%?title：不采用任何算法，直接將來兩幅圖像進(jìn)行相加融合
%?Data：2019-1-6
function?DirectImageFusion
clc;
close?all;
%讀取兩幅圖像
Ima1?=?imread（‘pic_1.tif‘）;????????
Ima2?=?imread（‘pic_2.tif‘）;
%將來兩幅圖像直接進(jìn)行相加，得到融合后圖像
fusionImage?=?（Ima1?+?Ima2）?/?2;
set（0?‘defaultFigurePosition‘?[300?190?800?350]）;?%設(shè)置一個(gè)寬800，高350的窗口
set（0?‘defaultFigureColor‘?[1?1?1]）;?%設(shè)置窗口的顏色?為白色
%計(jì)算信息熵
H?=?entropy（fusionImage）;?%matlab自帶
%?計(jì)算標(biāo)準(zhǔn)差
SD_fusionImage?=??standardDeviation（fusionImage）;
%?計(jì)算清晰度
G_Ima1?=?sharpness（Ima1）;
G_Ima2?=?sharpness（Ima2）;
G_fusionImage?=?sharpness（fusionImage）;
%?計(jì)算空間頻率
SF_Ima1?=?kongjianpinlv（Ima1）;
SF_Ima2?=?kongjianpinlv（Ima2）;
SF_fusionImage?=?kongjianpinlv（fusionImage）;
subplot（131）imshow（Ima1）title（‘圖像1‘‘position‘[128305]）;
subplot（132）imshow（Ima2）title（‘圖像2‘‘position‘[128305]）;
subplot（133）imshow（fusionImage）title（‘圖像直接相加融合‘‘position‘[128305]）;
end

%?計(jì)算標(biāo)準(zhǔn)差
function?SD?=?standardDeviation（pic）
pic?=?double（pic）;
%?計(jì)算矩陣元素的平均值
average?=?mean2（pic（:））;
[mn]?=?size（pic）;
sum?=?0;
for?i=1:m
????for?j=1:n
????????sum?=?sum?+?（pic（ij）-average）.^2;
????end
end
SD?=?sqrt（sum?/?（m*n））;
end

%?計(jì)算清晰度
function?G?=?sharpness（pic）
pic?=?double（pic）;
[st]?=?size（pic）;
n?=?s?*?t;
z?=?0;
for?i?=?1?:?（s?-?1）
????for?j?=?1?:?（t?-?1）
????????x?=?pic（ij+1）?-?pic（ij）;
????????y?=?pic（i+1j）?-?pic（ij）;
????????z?=?z?+?sqrt（（x.^2+y.^2）/2）;
????end
end
G?=?z?/?n;
end

%?計(jì)算空間頻率
function?SF?=?kongjianpinlv（Ima）
I?=?double（Ima）;
[mn]?=?size（I）;???????????????????
RF?=?0;
CF?=?0;?
for?i?=?1?:?m
????for?j?=?2?:?n
????????RF?=?RF?+?（I（ij）?-?I（ij-1））.^2;?%求矩陣中每個(gè)元素的平方
????end
end
RF?=?RF?/?（m?*?n）;
for?i?=?2?:?m
????for?j?=?1?:?n
????????CF?=?CF?+?（I（ij）?-?I（i-1j））.^2;
????end
end
CF?=?CF?/?（m?*?n）;
SF?=?sqrt（RF?+?CF）;
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-04-17?10:11??Graduation?Project\
?????目錄???????????0??2019-04-17?10:12??Graduation?Project\Direct?image?fusion\
?????文件????????1953??2019-04-17?10:12??Graduation?Project\Direct?image?fusion\DirectImageFusion.m
?????文件???????65769??2019-04-17?10:12??Graduation?Project\Direct?image?fusion\pic_1.tif
?????文件???????65769??2019-04-17?10:11??Graduation?Project\Direct?image?fusion\pic_2.tif
?????目錄???????????0??2019-04-17?10:12??Graduation?Project\Simple?image?fusion\
?????文件????????1970??2019-04-17?10:11??Graduation?Project\Simple?image?fusion\WeightAvg.m
?????文件???????65769??2019-04-17?10:12??Graduation?Project\Simple?image?fusion\pic_1.tif
?????文件???????65769??2019-04-17?10:11??Graduation?Project\Simple?image?fusion\pic_2.tif
?????目錄???????????0??2019-04-17?10:11??Graduation?Project\Wavelet?image?fusion\
?????文件????????4295??2019-04-17?10:33??Graduation?Project\Wavelet?image?fusion\LWeightGEnergy.m
?????文件????????5374??2019-04-17?10:11??Graduation?Project\Wavelet?image?fusion\LWeightGEnergy1.m
?????文件????????4941??2019-04-17?10:33??Graduation?Project\Wavelet?image?fusion\LWeightGMax.m
?????文件???????65769??2019-04-17?10:11??Graduation?Project\Wavelet?image?fusion\pic_1.tif
?????文件???????65769??2019-04-17?10:11??Graduation?Project\Wavelet?image?fusion\pic_2.tif