clear
I1=imread（‘C:\Users\攀攀\Desktop\融合\spot.bmp‘）;
I2=imread（‘C:\Users\攀攀\Desktop\融合\tm_743.bmp‘）;

I1=double（I1）;
I2=double（I2）;
mI1=mean2（I1）;
I2r=I2（::1）;
I2g=I2（::2）;
I2b=I2（::3）;
mI2r=mean2（I2r）;
mI2g=mean2（I2g）;
mI2b=mean2（I2b）;

%加權(quán)融合
r_r=sum（sum（（I1-mI1）.*（I2r-mI2r）））/（sqrt（sum（sum（（I1-mI1）.^2））*sum（sum（（I2r-mI2r）.^2））））;
r_g=sum（sum（（I1-mI1）.*（I2g-mI2g）））/（sqrt（sum（sum（（I1-mI1）.^2））*sum（sum（（I2g-mI2g）.^2））））;
r_b=sum（sum（（I1-mI1）.*（I2b-mI2b）））/（sqrt（sum（sum（（I1-mI1）.^2））*sum（sum（（I2b-mI2b）.^2））））;
RHr=（（1+abs（r_r））.*I1+（1-abs（r_r））.*I2r）/2;
RHg=（（1+abs（r_g））.*I1+（1-abs（r_g））.*I2g）/2;
RHb=（（1+abs（r_b））.*I1+（1-abs（r_b））.*I2b）/2;
RH1=cat（3RHrRHgRHb）;%融合成三維
RH1=uint8（RH1）;

%比值變換
B1=I2r;
B2=I2g;
B3=I2b;
D=I1;
DB1=（B1./（B1+B2+B3））.*D;
DB2=（B2./（B1+B2+B3））.*D;
DB3=（B3./（B1+B2+B3））.*D;
RH2=cat（3DB1DB2DB3）;
RH2=uint8（RH2）;

%乘積變換
db1=D.*B1/mI1;
db2=D.*B2/mI1;
db3=D.*B3/mI1;
RH3=cat（3db1db2db3）;
RH3=uint8（RH3）;

%融合后圖像顯示
subplot（131）;imshow（RH1）;title（‘加權(quán)融合后圖像‘）;
subplot（132）;imshow（RH2）;title（‘比值變換后圖像‘）;
subplot（133）;imshow（RH3）;title（‘乘積變換后圖像‘）;

%融合效果評(píng)價(jià)
RH1=double（RH1）;
RH2=double（RH2）;
RH3=double（RH3）;
hdr?=?mI2r;
hdg?=?mI2g;
hdb?=?mI2b;??%原始多光譜影像
hd1r?=?mean2（RH1（::1））;
hd1g?=?mean2（RH1（::2））;
hd1b?=?mean2（RH1（::3））;???%加權(quán)融合后
hd2r?=?mean2（RH2（::1））;
hd2g?=?mean2（RH2（::2））;
hd2b?=?mean2（RH2（::3））;????%比值變換后
hd3r?=?mean2（RH3（::1））;
hd3g?=?mean2（RH3（::2））;
hd3b?=?mean2（RH3（::3））;????%乘積變換后
%相關(guān)系數(shù)
r1R?=?sum（sum（（I2r?-?hdr）.*（RH1（::1）?-?hd1r）））?/?sqrt（（sum（sum（（I2r?-?hdr）.*（I2r?-?hdr））））?*?（sum（sum（（RH1（::1）?-?hd1r）.*（RH1（::1）?-?hd1r）））））;
r1G?=?sum（sum（（I2g?-?hdg）.*（RH1（::2）?-?hd1g）））?/?sqrt（（sum（sum（（I2g?-?hdg）.*（I2g?-?hdg））））?*?（sum（sum（（RH1（::2）?-?hd1g）.*（RH1（::2）?-?hd1g）））））;
r1B?=?sum（sum（（I2b?-?hdb）.*（RH1（::3）?-?hd1b）））?/?sqrt（（sum（sum（（I2b?-?hdb）.*（I2b?-?hdb））））?*?（sum（sum（（RH1（::3）?-?hd1b）.*（RH1（::3）?-?hd1b）））））;
r2R?=?sum（sum（（I2r?-?hdr）.*（RH2（::1）?-?hd2r）））?/?sqrt（（sum（sum（（I2r?-?hdr）.*（I2r?-?hdr））））?*?（sum（sum（（RH2（::1）?-?hd2r）.*（RH2（::1）?-?hd2r）））））;
r2G?=?sum（sum（（I2g?-?hdg）.*（RH2（::2）?-?hd2g）））?/?sqrt（（sum（sum（（I2g?-?hdg）.*（I2g?-?hdg））））?*?（sum（sum（（RH2（::2）?-?hd2g）.*（RH2（::2）?-?hd2g）））））;
r2B?=?sum（sum（（I2b?-?hdb）.*（RH2（::3）?-?hd2b）））?/?sqrt（（sum（sum（（I2b?-?hdb）.*（I2b?-?hdb））））?*?（sum（sum（（RH2（::3）?-?hd2b）.*（RH2（::3）?-?hd2b）））））;
r3R?=?sum（sum（（I2r?-?hdr）.*（RH3（::1）?-?hd3r）））?/?sqrt（（sum（sum（（I2r?-?hdr）.*（I2r?-?hdr））））?*?（sum（sum（（RH3（::1）?-?hd3r）.*（RH3（::1）?-?hd3r）））））;
r3G?=?sum（sum（（I2g?-?hdg）.*（RH3（::2）?-?hd3g）））?/?sqrt（（sum（sum（（I2g?-?hdg）.*（I2g?-?hdg））））?*?（sum（sum（（RH3（::2）?-?hd3g）.*（RH3（::2）?-?hd3g）））））;
r3B?=?sum（sum（（I2b?-?hdb）.*（RH3（::3）?-?hd3b）））?/?sqrt（（sum（sum（（I2b?-?hdb）.*（I2b?-?hdb））））?*?（sum（sum（（RH3（::3）?-?hd3b）.*（RH3（::3）?-?hd3b）））））;

%均值偏差：
b1R?=?abs（hdr?-?hd1r）?/?hd1r;
b1G?=?abs（hdg?-?hd1g）?/?hd1g;
b1B?=?abs（hdb?-?hd1b）?/?hd1b;
b2R?=?abs（hdr

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????3966??2015-10-29?13:57??融合\Fusion.m
?????文件??????187904??2015-10-29?14:54??融合\Thumbs.db
?????文件????17241982??2015-10-21?20:01??融合\spot.bmp
?????文件????51689158??2015-10-22?15:03??融合\tm_743.bmp
?????文件?????????305??2015-10-29?16:49??融合\融合精度.txt
?????目錄???????????0??2015-10-29?13:58??融合\