图像融合小波变换

资源包含加权平均、简单图像融合、基于小波变换.

% title：不采用任何算法，直接将来两幅图像进行相加融合
% Data：2019-1-6
function DirectImageFusion
clc;
close all;
%读取两幅图像
Ima1 = imread（‘pic_1.tif‘）;        
Ima2 = imread（‘pic_2.tif‘）;
%将来两幅图像直接进行相加，得到融合后图像
fusionImage = （Ima1 + Ima2） / 2;
set（0 ‘defaultFigurePosition‘ [300 190 800 350]）; %设置一个宽800，高350的窗口
set（0 ‘defaultFigureColor‘ [1 1 1]）; %设置窗口的颜色 为白色
%计算信息熵
H = entropy（fusionImage）; %matlab自带
% 计算标准差
SD_fusionImage =  standardDeviation（fusionImage）;
% 计算清晰度
G_Ima1 = sharpness（Ima1）;
G_Ima2 = sharpness（Ima2）;
G_fusionImage = sharpness（fusionImage）;
% 计算空间频率
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

% 计算标准差
function SD = standardDeviation（pic）
pic = double（pic）;
% 计算矩阵元素的平均值
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

% 计算清晰度
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

% 计算空间频率
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; %求矩阵中每个元素的平方
    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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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