【数字图像处理matlab】(乘积变换融合算法）

【数字图像处理matlab】(基于PCA的影像融合）

输入一张高分辨率的全色影像HighResultion，一张低分辨率的多光谱影像Multispectral，采用乘积变换影像融合算法，其中多光谱重采样采用imresize()函数的最近邻插值算法。

function F=multiplicative_fuse(HighResultion,Multispectral)
%调用代码----------------------------------------------
%Image1=imread('HR.jpg');Image2=imread('MS.jpg');
%multiplicative_fuse(Image1,Image2);
%-----------------------------------------------------

%乘积变换融合
B_m=im2double(Multispectral);
B_h=im2double(HighResultion);

[row,column,~]=size(B_h);
[~,~,bands]=size(B_m);
R= B_m(:,:,1);  r=imresize(R,[row,column]);
G= B_m(:,:,2);  g=imresize(G,[row,column]);
B= B_m(:,:,3);  b=imresize(B,[row,column]);

for i=1:row
    for j=1:column
        r1(i,j)=r(i,j)*B_h(i,j);
        b1(i,j)=b(i,j)*B_h(i,j);
        g1(i,j)=g(i,j)*B_h(i,j);
    end
end
B_new=cat(3,r1,g1,b1);
subplot(1,3,1),imshow(HighResultion);
title('原高分辨率图');
subplot(1,3,2),imshow(Multispectral);
title('原低分辨率图');
subplot(1,3,3);imshow(B_new);
title('乘积变换融合');

end

测试：