utils/matlab/cm_landwater.m

function [c] = cm_landwater(n,varargin)
% [c] = cm_landwater(n)
% [c] = cm_landwater(n,f)
%
%  n - is size of color map
%  f is scalar fraction that should be water (default=0.5)
%    or vector [hmin hmax] from which f will be calculated
%
% e.g.
% >> cm_landwater(64,caxis);
%
% Green-yellor land/blue water colormap

if nargin==1
 f=0.5;
else
 f=varargin{1};
end

if prod(size(f))==2
 hmin=f(1);hmax=f(2);
 f=-hmin/(hmax-hmin);
 f=max(0,f);f=min(1,f);
end

if mod(n,2)==0
 nb=round(n*f(1));
 ng=n-nb;
 nw=0;
else
 nb=floor(n*f(1));
 nw=1;
 ng=n-nb-nw;
end
b=bluewater(nb);
g=greenland(ng);
w=ones(nw,1)*[1 1 0.55];

c=[b' w' g']';
if nargout==0
 colormap(c);
end

%subplot(211)
%plot((0:(n-1))/(n-1),c(:,3:-1:1))
%subplot(212)
%pcol([-1:.1/n:1])
%axis off

% ----------------------------------------------------
function [c] = bluewater(n)
% Blue colormap

x=(0:n-1)'/(n-1);
c=((.7-.7*x)*[1 1 0]+max(0,(1-.65*x))*[0 0 1]);
c=c(end:-1:1,:);
%colormap(c);
% ----------------------------------------------------
function [c] = greenland(n)
% Green land colormap

x=(0:n-1)'/(n-1);
r=max(0,(-.25+1.8*x));
g=.4*(1+2.5*x);
b=0.5*max(0,(-.25+2.0*x));
i=find(r>1);
r(i)=max( 1.7-1*x(i), b(i));
i=find(g>1);
g(i)=max( 1.5-1*x(i), b(i));
c=r*[1 0 0]+g*[0 1 0]+b*[0 0 1];
c=min(max(0,c),1);
%colormap(c);
% ----------------------------------------------------
1	adcroft	1.1	function [c] = cm_landwater(n,varargin)
2			% [c] = cm_landwater(n)
3			% [c] = cm_landwater(n,f)
4			%
5			% n - is size of color map
6			% f is scalar fraction that should be water (default=0.5)
7			% or vector [hmin hmax] from which f will be calculated
8			%
9			% e.g.
10			% >> cm_landwater(64,caxis);
11			%
12			% Green-yellor land/blue water colormap
13
14			if nargin==1
15			f=0.5;
16			else
17			f=varargin{1};
18			end
19
20			if prod(size(f))==2
21			hmin=f(1);hmax=f(2);
22			f=-hmin/(hmax-hmin);
23			f=max(0,f);f=min(1,f);
24			end
25
26			if mod(n,2)==0
27			nb=round(n*f(1));
28			ng=n-nb;
29			nw=0;
30			else
31			nb=floor(n*f(1));
32			nw=1;
33			ng=n-nb-nw;
34			end
35			b=bluewater(nb);
36			g=greenland(ng);
37			w=ones(nw,1)*[1 1 0.55];
38
39			c=[b' w' g']';
40			if nargout==0
41			colormap(c);
42			end
43
44			%subplot(211)
45			%plot((0:(n-1))/(n-1),c(:,3:-1:1))
46			%subplot(212)
47			%pcol([-1:.1/n:1])
48			%axis off
49
50			% ----------------------------------------------------
51			function [c] = bluewater(n)
52			% Blue colormap
53
54			x=(0:n-1)'/(n-1);
55			c=((.7-.7x)[1 1 0]+max(0,(1-.65x))[0 0 1]);
56			c=c(end:-1:1,:);
57			%colormap(c);
58			% ----------------------------------------------------
59			function [c] = greenland(n)
60			% Green land colormap
61
62			x=(0:n-1)'/(n-1);
63			r=max(0,(-.25+1.8*x));
64			g=.4(1+2.5x);
65			b=0.5max(0,(-.25+2.0x));
66			i=find(r>1);
67			r(i)=max( 1.7-1*x(i), b(i));
68			i=find(g>1);
69			g(i)=max( 1.5-1*x(i), b(i));
70			c=r[1 0 0]+g[0 1 0]+b*[0 0 1];
71			c=min(max(0,c),1);
72			%colormap(c);
73			% ----------------------------------------------------