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

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

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

function F=Broveymelt(HighResultion,Multispectral)
%调用代码-------------------------------------------------
%Image1=imread('HR.jpg');Image2=imread('MS.jpg');
%%Image1为高分辨率，Image2为多光谱
%Broveymelt(Image1,Image2);
%------------------------------------------------------

subplot(1,3,1),imshow(HighResultion);
title('原高分辨率图');
subplot(1,3,2),imshow(Multispectral);
title('原低分辨率图');
HighResultion=im2double(HighResultion);
Multispectral=im2double(Multispectral);
[row,column,bands]=size(HighResultion);   %求高分辨率的size
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  % B = imresize(A,[nrows ncols],method)
  %method的几种可选值：
          %'nearest'最近邻插值（默认值）
          % 'bilinear'双线性插值
          % 'bicubic'三次卷积插值
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
r=Multispectral(:,:,1);   R=imresize(r,[row,column]);
g=Multispectral(:,:,2);   G=imresize(g,[row,column]);
b=Multispectral(:,:,3);   B=imresize(b,[row,column]);

Xr=zeros(row,column);Xg=zeros(row,column);Xb=zeros(row,column);
for i=1:row
    for j=1:column
        sum=R(i,j)+G(i,j)+B(i,j);
        Xr(i,j)=HighResultion(i,j)*R(i,j)/sum;
        Xg(i,j)=HighResultion(i,j)*G(i,j)/sum;
        Xb(i,j)=HighResultion(i,j)*B(i,j)/sum;
    end
end

X=cat(3,Xr,Xg,Xb);
subplot(1,3,3),imshow(X);
title('Brovey变换融合后的图像');

end

测试：