brovey transform比值融合算法matlab,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

</