MATLAB实现Brovey图像融合

自定义函数：

function BF=Brovey_fuse(Hyperspectral_image,High_resolution_image)
x0=imread(Hyperspectral_image); %高光谱影像
[a,b,c]=size(x0);
x=double(x0)/;
y=imread(High_resolution_image); %高分辨率灰度影像
y1=double(y)/;
xx=zeros(a,b,c);
p=zeros(a,b,c);
for f=:a
    for e=:b
        xx(f,e)=x(f,e,)+x(f,e,)+x(f,e,);
        p(f,e,)=x(f,e,)*y1(f,e)/xx(f,e);
        p(f,e,)=x(f,e,)*y1(f,e)/xx(f,e);
        p(f,e,)=x(f,e,)*y1(f,e)/xx(f,e);
    end
end
BF=uint8(round(p*));
%imshow(BF);

在命令行中输入：

>> a=Brovey_fuse('f1.jpg','f2.bmp');
>> imshow(a)

f1.jpg

f2.bmp

融合结果：