Matlab的colormap函数详解(多个colormap)

piday123

我们有一张图片，读进来后，转换为gray，然后打算画两张图，一张用gray的色图，一张用jet的色图。
8 m: @6 e+ ~0 i首先我们读取图片：
0 f, ]4 ?# z8 L- D" Tclear all;
( c; t- I' H% N! I& X+ Iclose all;
clc;
%%
imgname='.\frog.bmp';
[imagRGB]     = imread( imgname );
( L4 Y1 X& r0 x9 Z, V4 P
5 L. @+ s4 [6 y* m9 Vif (length(size(imagRGB)) > 2)
    imagGray = rgb2gray(imagRGB);
else
9 {" q3 q: s* n" \7 G  x; n6 ^6 [6 e    imagGray = imagRGB;
5 t. s- B9 Y6 k( o0 l2 X8 nend
clear('imagRGB');
# r. i4 Q* E: `& l9 ydata =  double(imagGray);
data =flipdim(data,1);  
   % Flip linear dimensions (image coordinate system starts at upper left
   % corner )

7 r5 `7 \4 D7 Y2 T
5 \/ P, l5 t4 i, J0 x2 h6 V8 m4 k方法1： 我们生成两个figure，分开画两个图
4 f; ]  \, y$ \" t* S6 k, Rfig1=figure;
plotHandles(1) = pcolor( data );
1 ?5 _  U6 z7 N  H+ Jset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
% I' j: H! R) Ecolormap gray;
  K- a# r1 I( s) ~4 f
! |- C3 B$ b2 ]fig2=figure;
plotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255])
colormap jet;

! C# M& i( j) H7 ]/ M& Y. W0 d( _
2 h/ }: `- C+ j1 {8 `+ E
结果如下：
1 r1 B& I/ a" |  x
2 ?  r% k/ m: p1 X+ Q注意：这是两个figure，每个figure可以使用一个唯一的colormap
5 k( A+ z5 k# X: E- l1 x

, o. Z0 H) y9 L! j方法二： 一个figure，subplot
fig3=figure;
! G& `0 g! W6 K2 q! h
subplot(1,2,1);
plotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
4 \  B3 z: D% J) _colormap gray;
' M7 o: d/ U1 h# w) J' E
subplot(1,2,2);
4 u/ U/ g. E( p, j! t: wplotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
; \) w: s8 P, R  Ycaxis([0,255])
colormap jet;


# `. X1 k! A) v+ b5 {) _结果如何呢？

) Y% _/ L) Z$ n/ h# q$ S; @- W

后面定义的colormap会把前面的覆盖掉。因此最终两张图片都使用了jet的色图。
/ Z1 Q% K/ h, l
4 m* L* @0 I- A  H" O1 ~7 n方法三： 既然只能使用一个colormap，那么我们就自己动手组装我们需要的色图
options.colormaps{1}.offset     = 0;
0 p. |. [" D/ T7 U5 K' Q9 moptions.colormaps{1}.num        = 256;  % Number of colors
% }6 N2 j& n2 E+ M$ N: |, Qoptions.colormaps{1}.cmin       = 0    + options.colormaps{1}.offset;
options.colormaps{1}.cmax       = 255  + options.colormaps{1}.offset;
2 [# [* [6 t, Z; L( noptions.colormaps{1}.map        = gray(256);

options.colormaps{2}.offset     = 256;
options.colormaps{2}.num        = 256;  % Number of colors
options.colormaps{2}.cmin       = 0    + options.colormaps{1}.offset;
2 M0 ~7 ?/ T( s! Noptions.colormaps{2}.cmax       = 255  + options.colormaps{1}.offset;
4 p3 }5 O* o+ l1 [4 ^% V- soptions.colormaps{2}.map        = jet(256);
+ G! k, g, m/ {
for iCMap = 1 : length(options.colormaps)
/ K5 e. l1 v, `' p" H    if iCMap == 1
        cmap = options.colormaps{iCMap}.map;
    else
! Q2 U7 |, t" a; C3 h6 x        cmap = [cmap ; options.colormaps{iCMap}.map ];
    end
! P9 A# F9 _: c2 K5 t: ~  Send
1 ?$ K* f+ f; E4 N% |- N
8 ~( N; O7 Y7 z0 u) S: p
我们定义了两个色图，然后把两个色图拼装起来。

fig4 = figure;
1 ?9 S* ?' o( ]colormap(cmap);

subplot(1,2,1);
data1 = data +options.colormaps{1}.offset;
plotHandles(1) = pcolor( data1 );
set( plotHandles(1) , 'EdgeColor' , 'none');
$ @$ J: u# N& n( s. s0 c' [caxis([0 , size(cmap,1) - 1]);


subplot(1,2,2);
data2 = data + options.colormaps{2}.offset;
) o$ T4 l0 C# z- l6 B8 qplotHandles(2) = pcolor( data2 );

set( plotHandles(2) , 'EdgeColor' , 'none');
caxis([0 , size(cmap,1) - 1]);


4 Z! K- I( C+ W. ~# a) I结果如下：

9 u3 D0 n/ g* \: W4 E* ]: d
( f  {1 |6 O" W实际上我们只有一个colormap，但达到了多个colormap的效果。

7 R6 }) m( q$ b- e4 K  K

SsaaM7

Matlab的colormap函数详解(多个colormap)