灰度图像转伪彩色

先只贴出代码

clear all;close all;clc;
for k=1:6
    I=imread(['H:\mwtmodel\ths_featuremap\Resize_2\Original_' num2str(k) '.tiff']);
    % figure;
    % imshow(I);
    I=double(I);
    [m,n]=size(I);
    L=256;
    for i=1:m
        for j=1:n
            if I(i,j)<L/4
                R(i,j)=0;
                G(i,j)=4*I(i,j);
                B(i,j)=L;
            else if I(i,j)<=L/2
                    R(i,j)=0;
                    G(i,j)=L;
                    B(i,j)=-4*I(i,j)+2*L;
                else if I(i,j)<=3*L/4
                        R(i,j)=4*I(i,j)-2*L;
                        G(i,j)=L;
                        B(i,j)=0;
                    else
                        R(i,j)=L;
                        G(i,j)=-4*I(i,j)+4*L;
                        B(i,j)=0;
                    end
                end
            end
        end
    end
    for i=1:m
        for j=1:n
            G2C(i,j,1)=R(i,j);
            G2C(i,j,2)=G(i,j);
            G2C(i,j,3)=B(i,j);
        end
    end
    G2C=G2C/256;
    imwrite(G2C,['C:\Users\Administrator\Desktop\ths_featuremap\Resize_2\Original_' num2str(k) '.tiff']);
end