自定义函数:
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
</