我在用matlab画三维场强图时,找不到合适的colormap。最接近的colormap是matlab自带的jet,即随着函数值的增大颜色从蓝过渡到红。但由于场强图主要表征场的幅值特性,颜色区分只与函数值的绝对值相关,直接用jet函数就会导致场强最小位置的颜色为中间过渡色,而非蓝色。因此我利用colormap颜色RGB值渐变过渡的关系,编写了一个适用于场强表示的colormap。
程序效果见下图,本程序参考jet的色彩过渡关系,为了编写方便提升了色彩的纯度,采用blue[0,0,1]/cyan[0,1,1]/ yellow[1,1,0]/red[1,0,0]的纯色rgb值。
科研初学者,程序没啥含量,就是自己的一点小理解。
%% 电场colormap设定
n_step=9;
color_row=n_step*6+7;
color_map=zeros(color_row,3);
color_r = 0:1/(n_step+1):1;
color_g = 0:1/(n_step+1):1;
color_b = 0:1/(n_step+1):1;
% red to yellow
color_map(1:1+(n_step+1),1)=1;
color_map(1:1+(n_step+1),2)=color_g;
a=1+(n_step+1);
% yellow to cyan
color_map(a:a+(n_step+1),1)=fliplr(color_r);
color_map(a:a+(n_step+1),2)=1;
color_map(a:a+(n_step+1),3)=color_b;
a=a+(n_step+1);
% cyan to blue
color_map(a:a+(n_step+1),2)=fliplr(color_b);
color_map(a:a+(n_step+1),3)=1;
a=a+(n_step+1);
% blue to cyan
color_map(a:a+(n_step+1),2)=color_b;
color_map(a:a+(n_step+1),3)=1;
a=a+(n_step+1);
% cyan to yellow
color_map(a:a+(n_step+1),1)=color_r;
color_map(a:a+(n_step+1),2)=1;
color_map(a:a+(n_step+1),3)=fliplr(color_b);
a=a+(n_step+1);
% yellow to red
color_map(a:a+(n_step+1),1)=1;
color_map(a:a+(n_step+1),2)=fliplr(color_g);
colormap(color_map);