用matlab将灰度图转换为彩色图

function R=gray2rgb(img1,img2)
% img1 - Source Image  (gray image)   
% img2 - Selected color image for coloring the gray image. 

clc;
warning off;
imt=imread(img1);
ims=imread(img2);
[sx sy sz]=size(imt);
[tx ty tz]=size(ims);
if sz~=1
    imt=rgb2gray(imt);
end
if tz~=3
    disp ('img2 must be a color image (not indexed)');
else
    imt(:,:,2)=imt(:,:,1);
    imt(:,:,3)=imt(:,:,1);

% Converting to ycbcr color space
    nspace1=rgb2ycbcr(ims);
    nspace2= rgb2ycbcr(imt);
    ms=double(nspace1(:,:,1));
    mt=double(nspace2(:,:,1));
    m1=max(max(ms));
    m2=min(min(ms));
    m3=max(max(mt));
    m4=min(min(mt));
    d1=m1-m2;
    d2=m3-m4;
% Normalization
    dx1=ms;
    dx2=mt;
    dx1=(dx1*255)/(255-d1);
    dx2=(dx2*255)/(255-d2);
    [mx,my,mz]=size(dx2);
%Luminance Comparison
    disp('Please wait..................');
    for i=1:mx
        for j=1:my
             iy=dx2(i,j);
             tmp=abs(dx1-iy);
             ck=min(min(tmp));
             [r,c] = find(tmp==ck);
             ck=isempty(r);
             if (ck~=1)            
                 nimage(i,j,2)=nspace1(r(1),c(1),2);
                 nimage(i,j,3)=nspace1(r(1),c(1),3);
                 nimage(i,j,1)=nspace2(i,j,1);           
            end
         end
     end
    rslt=ycbcr2rgb(nimage)
    figure,imshow(uint8(imt));
    figure,imshow(uint8(rslt));
    R=uint8(rslt);
    toc
end  

subplot(221),imshow('pikaqiu.png');
subplot(222),imshow('maliao.png');
colorIm 
 = gray2rgb('pikaqiu.png','maliao.png');
subplot(223),imshow(colorIm); 

转载于： https://blog.csdn.net/u012150360/article/details/67656823