-
Notifications
You must be signed in to change notification settings - Fork 15
/
Copy pathanalysis_Reference.m
40 lines (29 loc) · 917 Bytes
/
analysis_Reference.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
function [EN,Qabf,SCD,FMI_w,FMI_dct,SSIM, MS_SSIM,FMI_pixel,Nabf,MI] = analysis_Reference(image_f,image_ir,image_vis)
[s1,s2] = size(image_ir);
imgSeq = zeros(s1, s2, 2);
imgSeq(:, :, 1) = image_ir;
imgSeq(:, :, 2) = image_vis;
image1 = im2double(image_ir);
image2 = im2double(image_vis);
image_fused = im2double(image_f);
%EN
EN = entropy(image_fused);
%MI
MI = analysis_MI(image_ir,image_vis,image_f);
%Qabf
Qabf = analysis_Qabf(image1,image2,image_fused);
%FMI
FMI_pixel = analysis_fmi(image1,image2,image_fused);
FMI_dct = analysis_fmi(image1,image2,image_fused,'dct');
FMI_w = analysis_fmi(image1,image2,image_fused,'wavelet');
%Nabf
Nabf = analysis_nabf(image_fused,image1,image2);
%SCD
SCD = analysis_SCD(image1,image2,image_fused);
% SSIM_a
SSIM1 = ssim(image_fused,image1);
SSIM2 = ssim(image_fused,image2);
SSIM = (SSIM1+SSIM2)/2;
%MS_SSIM
[MS_SSIM,t1,t2]= analysis_ms_ssim(imgSeq, image_f);
end