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

piday123

我们有一张图片，读进来后，转换为gray，然后打算画两张图，一张用gray的色图，一张用jet的色图。
首先我们读取图片：
clear all;
close all;
clc;
%%
) C/ m. i6 P% C2 j3 ], Pimgname='.\frog.bmp';
[imagRGB]     = imread( imgname );

  P' M4 `6 O# m9 I7 q& }) sif (length(size(imagRGB)) > 2)
    imagGray = rgb2gray(imagRGB);
else
    imagGray = imagRGB;
end
clear('imagRGB');
data =  double(imagGray);
data =flipdim(data,1);  
* R0 A  a6 G; h% k, W   % Flip linear dimensions (image coordinate system starts at upper left
   % corner )
7 B7 y5 }. I9 U: w- Y7 q) F7 ?

& o2 }2 J# r+ N5 J& g, z7 c方法1： 我们生成两个figure，分开画两个图
1 I1 q5 I0 `6 i" Efig1=figure;
# r+ X+ O8 t  S: E6 \2 W, u7 splotHandles(1) = pcolor( data );
set( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
* u) b/ \. g) g- S  xcolormap gray;
4 _5 L' x5 t) I- W' l8 K
fig2=figure;
% _& x( g. ^8 mplotHandles(1) = pcolor( data );
4 t7 j% x& U2 D7 s/ J9 l7 C1 O+ |set( plotHandles(1) , 'EdgeColor' , 'none');
/ h0 o' ]7 r; V2 @1 v$ {caxis([0,255])
9 `: p! S, |5 O6 T- T& Mcolormap jet;


2 @$ l& _6 X3 z" I: t! f/ @
" |+ z- I, G7 J/ {; T结果如下：

注意：这是两个figure，每个figure可以使用一个唯一的colormap
# u, e& z& a  r

方法二： 一个figure，subplot
fig3=figure;

subplot(1,2,1);
; p% Q/ b2 |. d. A# V1 @" ?6 O% rplotHandles(1) = pcolor( data );
, k0 Y" z" ^! M7 s2 w8 eset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255]);
colormap gray;

- n4 ]" W) L4 I5 g2 \5 B. Tsubplot(1,2,2);
plotHandles(1) = pcolor( data );
! a- v' v8 ?5 s7 w. V4 Xset( plotHandles(1) , 'EdgeColor' , 'none');
caxis([0,255])
3 K  t! }$ N8 q+ a5 H7 hcolormap jet;
/ D& ?7 t! t, m$ n$ H
8 a1 n/ t1 ?% j! a% M3 M
( S: F5 [( [( g5 I6 @结果如何呢？

) d' G2 o; n+ F' l4 B) j
, T" c( w7 Z" N1 Y, d) ?% |4 o
& Z+ _  r7 z! g1 r后面定义的colormap会把前面的覆盖掉。因此最终两张图片都使用了jet的色图。

方法三： 既然只能使用一个colormap，那么我们就自己动手组装我们需要的色图
options.colormaps{1}.offset     = 0;
. h1 D1 R( ?1 {) }+ qoptions.colormaps{1}.num        = 256;  % Number of colors
options.colormaps{1}.cmin       = 0    + options.colormaps{1}.offset;
0 o4 }/ u# g0 O/ N( ~& @" z0 ooptions.colormaps{1}.cmax       = 255  + options.colormaps{1}.offset;
options.colormaps{1}.map        = gray(256);
' ]2 P) u1 P# [' z; z* ^* L
options.colormaps{2}.offset     = 256;
0 V5 }6 N; K, D3 o, Noptions.colormaps{2}.num        = 256;  % Number of colors
options.colormaps{2}.cmin       = 0    + options.colormaps{1}.offset;
options.colormaps{2}.cmax       = 255  + options.colormaps{1}.offset;
" m; T# ?( M7 n" B: eoptions.colormaps{2}.map        = jet(256);
5 P& b& m5 U! k4 B) K
# g: E2 ~8 {& w$ K: U) dfor iCMap = 1 : length(options.colormaps)
    if iCMap == 1
  L- {$ F) z( Q( p8 m        cmap = options.colormaps{iCMap}.map;
    else
        cmap = [cmap ; options.colormaps{iCMap}.map ];
    end
. s! C* l. Y* uend

: n9 }3 Q& u: b! W, f! O
2 W* u$ o3 a  k3 y我们定义了两个色图，然后把两个色图拼装起来。
2 n" I, p5 k' @' p3 g& B
fig4 = figure;
' W3 w$ L' ?# Hcolormap(cmap);
2 d" f) _3 M' D+ j  |
; f) E* k* K8 d6 z, Q% @# K% {7 |1 Esubplot(1,2,1);
data1 = data +options.colormaps{1}.offset;
) ^/ a5 r; i2 R6 b$ J, ZplotHandles(1) = pcolor( data1 );
1 e$ B$ {) T& Q8 A' Eset( plotHandles(1) , 'EdgeColor' , 'none');
. _% S* ^, l4 V; f. ?1 kcaxis([0 , size(cmap,1) - 1]);
& ~: B  t; v, _& H

$ l' ~- \- q; J7 n' Asubplot(1,2,2);
data2 = data + options.colormaps{2}.offset;
: Q" G0 p9 w( S+ KplotHandles(2) = pcolor( data2 );
% h6 u, _/ j8 i  \
' s. [/ d8 [+ c3 v* h) O3 u3 x) wset( plotHandles(2) , 'EdgeColor' , 'none');
& K0 y5 l7 @$ A: p" q1 {( kcaxis([0 , size(cmap,1) - 1]);

) W2 ~1 a1 z7 W" T* N: Y
& k" ~% t  J& B6 N* m结果如下：

6 P  j1 L  D# R6 T: L$ Q0 k
6 f/ y' \: s" g- G( j实际上我们只有一个colormap，但达到了多个colormap的效果。
# a$ M3 U; ]0 F, ~+ [* w

SsaaM7

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