%-----------------------------------------------------------------------------
%灰度圖像的偽彩色處理：給定一個灰度－彩色轉(zhuǎn)換函數(shù)，對一灰度圖像進行偽彩色處理。
%轉(zhuǎn)換函數(shù)如下表：
%輸入灰度級	輸出彩色
%0～31	藍色
%32～63	綠色
%64～95	淡藍色
%96～127	紫色
%128～159	紅色
%160～191	橙色
%192～223	黃色
%224～255	淡黃色
%--------------------------------------------------------------------------

clear;

%讀入源圖像
A=imread('DIP_070729 趙明_灰度圖像的偽彩色處理_02-lena.tif');
[M,N]=size(A);%讀取圖像的大小

for i=1:M;
    for j=1:N;
        if A(i,j)<32                 %轉(zhuǎn)化為藍色
            B(i,j,:,:,:)=[0,0,1];
        elseif A(i,j)<64             %轉(zhuǎn)化為綠色
            B(i,j,:,:,:)=[0,1,0];
        elseif A(i,j)<96             %轉(zhuǎn)化為淡藍色
            B(i,j,:,:,:)=[0.4,0.9,1];
        elseif A(i,j)<128            %轉(zhuǎn)化為紫色
            B(i,j,:,:,:)=[1,0,1];
        elseif A(i,j)<160            %轉(zhuǎn)化為紅色
            B(i,j,:,:,:)=[1,0,0];
        elseif A(i,j)<192            %轉(zhuǎn)化為橙色
            B(i,j,:,:,:)=[1,0.8,0];
        elseif A(i,j)<224            %轉(zhuǎn)化為黃色
            B(i,j,:,:,:)=[1,1,0];
        elseif A(i,j)<255            %轉(zhuǎn)化為淡黃色
            B(i,j,:,:,:)=[1,0.8,0];
        end
    end
end

%--------------------------------------------------------------------------
%將處理結(jié)果存入目標(biāo)文件夾并畫圖
imwrite(B,'elaine_1.tif','tiff');
clf;
subplot(1,2,1),imshow(A);title('原始的灰度圖像');
subplot(1,2,2),imshow(B);title('偽彩色圖像');