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

piday123

我们有一张图片，读进来后，转换为gray，然后打算画两张图，一张用gray的色图，一张用jet的色图。
首先我们读取图片：
* i. |7 J/ D0 U, |# }7 bclear all;
7 {4 \) W2 P" E8 \/ {' O; oclose all;
clc;
%%
imgname='.\frog.bmp';
, k/ H- n& }8 u' Z( S+ b[imagRGB]     = imread( imgname );

( R, V- k( b7 ^0 ]7 g$ j" G, `if (length(size(imagRGB)) > 2)
    imagGray = rgb2gray(imagRGB);
+ p5 |  }$ G  X1 V* ?: l8 welse
9 f: r4 e6 o+ W) Q    imagGray = imagRGB;
% M/ x7 M5 Y  v* w6 r4 j/ w! Wend
clear('imagRGB');
data =  double(imagGray);
data =flipdim(data,1);  
   % Flip linear dimensions (image coordinate system starts at upper left
3 U, f9 z/ h1 c! z/ O   % corner )
- M+ l- z  L2 s1 a* C% y

方法1： 我们生成两个figure，分开画两个图
0 L' T7 S, c( W4 O, a2 Z( pfig1=figure;
plotHandles(1) = pcolor( data );
* T: z: A) p3 yset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
" y# @2 |" g3 y1 h2 ^colormap gray;
! ?( W, w2 o4 i, [+ M  Y- `
' g# U. t, C# V$ s7 A" b# Kfig2=figure;
8 P. k& _. U7 `+ l" B- hplotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
/ H, E, q/ S+ X5 q7 W* O2 ]caxis([0,255])
colormap jet;
/ o; y/ b& M% N/ [' J9 n1 `) i


结果如下：
+ g$ L  k3 S, c: a& P" z
  |, N; R" [2 G% a( I注意：这是两个figure，每个figure可以使用一个唯一的colormap

! M( ~6 D0 e3 H, R2 {9 ]
方法二： 一个figure，subplot
& r* s7 N0 f9 b9 q5 h+ Z2 efig3=figure;
7 w% a+ C9 L5 \9 v2 _  ^
subplot(1,2,1);
plotHandles(1) = pcolor( data );
2 J3 ^6 L3 y9 Lset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
colormap gray;

subplot(1,2,2);
plotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
& d, [5 k9 L, C( X& N. n, X; Ecaxis([0,255])
colormap jet;
/ }. E% [& }3 ~' @; u

结果如何呢？

" r' }9 P- V, G$ h# ~/ @+ i' p1 v7 K
- o6 d; ]) P' a, x
后面定义的colormap会把前面的覆盖掉。因此最终两张图片都使用了jet的色图。
  V; z3 ^7 b) l! N
0 c' O+ G) r3 t5 _0 `方法三： 既然只能使用一个colormap，那么我们就自己动手组装我们需要的色图
options.colormaps{1}.offset     = 0;
6 ~3 l- ^& V- V& J, Uoptions.colormaps{1}.num        = 256;  % Number of colors
2 O9 S4 y9 n7 Z3 b: ]options.colormaps{1}.cmin       = 0    + options.colormaps{1}.offset;
. |, `  Y$ L) }options.colormaps{1}.cmax       = 255  + options.colormaps{1}.offset;
options.colormaps{1}.map        = gray(256);
1 A" \+ }% V. M' D- N
" W( `2 ~7 v  t) P3 w" Toptions.colormaps{2}.offset     = 256;
options.colormaps{2}.num        = 256;  % Number of colors
1 P- I3 U7 [' f" c8 J9 @4 W7 L, ]  Loptions.colormaps{2}.cmin       = 0    + options.colormaps{1}.offset;
2 x. G8 r. K7 A3 w) ?7 V& k; Woptions.colormaps{2}.cmax       = 255  + options.colormaps{1}.offset;
options.colormaps{2}.map        = jet(256);
- [' v5 x* G. ^+ _% v/ b
5 [: V8 ^# d6 U. I" X8 Ofor iCMap = 1 : length(options.colormaps)
    if iCMap == 1
: _/ }7 b: g: U: c3 @1 N+ X" d        cmap = options.colormaps{iCMap}.map;
    else
+ E/ y7 {& l  M" u' K" n        cmap = [cmap ; options.colormaps{iCMap}.map ];
3 U. T: m! ]/ @1 _: n1 o! ^% e3 H    end
2 w) h: `! U$ x+ p  Oend
6 V' c0 z2 G! R- n6 q7 s& z
  N% ~; V/ \( M9 {8 P: w5 T
我们定义了两个色图，然后把两个色图拼装起来。

+ Q, u' `) z) s( ~  b0 rfig4 = figure;
6 U! J4 g3 g1 {% H# Lcolormap(cmap);

. }' p& u' ^5 R. hsubplot(1,2,1);
data1 = data +options.colormaps{1}.offset;
plotHandles(1) = pcolor( data1 );
- v; [! o  v& H5 w6 jset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0 , size(cmap,1) - 1]);


subplot(1,2,2);
0 N% J0 y& d, L. d9 y. x& cdata2 = data + options.colormaps{2}.offset;
plotHandles(2) = pcolor( data2 );
8 w( x" C2 H' P  F& a. H% u
set( plotHandles(2) , 'EdgeColor' , 'none');
caxis([0 , size(cmap,1) - 1]);
, k2 T' \8 a2 _7 I8 S" P6 I; V
0 Y5 N+ c9 I' {6 l) j
, J* z5 C5 r) W7 K结果如下：

/ I, z! o% u0 B1 R0 |9 \
5 O! {9 D1 \+ s, ]1 `$ E. {' T: k, e实际上我们只有一个colormap，但达到了多个colormap的效果。
9 h3 }( _5 R# `( s+ ^. F/ j! N

SsaaM7

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