function out = aaltoGradient(numLevels) if nargin<1 numLevels=128; end col1=[0 hex2dec('65') hex2dec('bd')]; col2=255*[1 1 1]; col3=[hex2dec('ed') hex2dec('29') hex2dec('39')]; levs=numLevels/2; colGrad=zeros(levs,3); colGrad2=zeros(levs,3); for k=1:3 colGrad(:,k)=col1(k):(col2(k)-col1(k))/(levs-1):col2(k); end for k=1:3 colGrad2(:,k)=col2(k):(col3(k)-col2(k))/(levs-1):col3(k); end out=([colGrad; colGrad2])/255; end