I = imread('saturn.tif');
%把灰度圖像I分成8層
I=double(I) ;
J=floor(I/32);
%對圖像數組進行等分層處理
J1=floor(J/4);
J2=rem(floor(J/2),2);
J3=rem(J,2);
%灰度與彩色變換關系
PC(:,:,1)=J1;
PC(:,:,2)=J2;
PC(:,:,3)=J3;
%組成m×n×3圖像數組
imshow(PC)