%首先讀入圖像文件
A=imread('myimage.tif');%A是一個unit8類型的S1×S2×3數組
A=double(A);
Xrgb=0.2990*A(:,:,1)+0.5870*A(:,:,2)+0.1140*A(:,:,3);
NbColors=255;
X=wcodemat(Xrgb,NbColors);
map=pink(NbColors);
image(X);
colormap(map);
colorbar;