clear;

A=imread('Einstein.tif');
[M,N]=size(A);

for i=1:M;
    for j=1:N;
        if A(i,j)<32                 
            B(i,j,:,:,:)=[0,0,1];
        elseif A(i,j)<64            
            B(i,j,:,:,:)=[0,1,0];
        elseif A(i,j)<96            
            B(i,j,:,:,:)=[0.4,0.9,1];
        elseif A(i,j)<128            
            B(i,j,:,:,:)=[1,0,1];
        elseif A(i,j)<160           
            B(i,j,:,:,:)=[1,0,0];
        elseif A(i,j)<192            
            B(i,j,:,:,:)=[1,0.8,0];
        elseif A(i,j)<224            
            B(i,j,:,:,:)=[1,1,0];
        elseif A(i,j)<255            
            B(i,j,:,:,:)=[1,0.8,0];
        end
    end
end

imwrite(B,'Einstein2.tif','tiff');
clf;
subplot(1,2,1),imshow(A);title('原始的灰度圖像');
subplot(1,2,2),imshow(B);title('偽彩色圖像');