bone(varargin)

%BONE  Black-blue-white colormap
%
% Examples:
%   map = bone;
%   map = bone(len);
%   B = bone(A);
%   B = bone(A, lims);
%
% Similar to MATLAB's bone function, but also able to return a concise
% colormap table.
%
% The function can additionally be used to convert a real-valued array into
% a truecolor array using the colormap.
%
% IN:
%   len - Scalar length of the output colormap. If len == Inf the concise
%         table is returned. Default: len = size(get(gcf, 'Colormap'), 1);
%   A - Non-scalar numeric array of real values to be converted into
%       truecolor.
%   lims - 1x2 array of saturation limits to be used on A. Default:
%          [min(A(:)) max(A(:))].
%
% OUT:
%   map - (len)xJ colormap table. J = 3, except in the concise case, when
%         J = 4, map(1:end-1,4) giving the relative sizes of the 
%         inter-color bins.
%   B - size(A)x3 truecolor array.

% Copyright: Oliver Woodford, 2009

function map = bone(varargin)
map = [0 0 0 3; 21 21 29 3; 42 50 50 2; 64 64 64 1]/64;
map = colormap_helper(map, varargin{:});