function?[]=Brovey_fusion（p1p2）
%%基于Brovey變換融合算法：p1代表彩色的多波段圖像；p2代表黑白的但波段圖像
%以d多光譜圖像作為圖像尺寸參考
%多光譜影像
p1=p1;p2=p2;
dp1=double（p1）;?
%高分辨率灰度影像
dp2=double（p2）;
[hwband]=size（p2）;
%將低分辨率的多光譜圖像進行重采樣，使其與高分影像像素數一致
B1=Resampling_Bilinearint（dp1（::1）hw）;
B2=Resampling_Bilinearint（dp1（::2）hw）;
B3=Resampling_Bilinearint（dp1（::3）hw）;
dp1=cat（3B1B2B3）;
dp1=double（dp1）;

xx=zeros（hwband）;
p=zeros（hwband）;

for?f=1:h
????for?e=1:w
????????xx（fe）=dp1（fe1）+dp1（fe2）+dp1（fe3）;???
????????p（fe1）=dp1（fe1）*dp2（fe）/xx（fe）;
????????p（fe2）=dp1（fe2）*dp2（fe）/xx（fe）;
????????p（fe3）=dp1（fe3）*dp2（fe）/xx（fe）;
????end
end

p=uint8（p）;
dp1=uint8（dp1）;
dp2=uint8（dp2）;
subplot（131）;imshow（dp1）;title（‘多光譜圖像1‘）;
subplot（132）;imshow（dp2）;title（‘高分辨率圖像2‘）;
subplot（133）;imshow（p）;title（‘1+2融合后圖像‘）;


end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????968??2019-11-09?15:04??fusion\Brovey_fusion.m

?????文件????????530??2019-11-15?08:14??fusion\HSI_fusion.m

?????文件???????1492??2019-11-09?15:58??fusion\PCA_fusion.m

?????目錄??????????0??2019-11-29?13:43??fusion

-----------?---------??----------?-----??----

?????????????????2990????????????????????4