MATLAB 图像融合评价算法

function laplacian(ori_A,ori_B)
path(path,'FusionEvaluation/')
%Ori_A = imread('clockA_tif2bmp (256x256).BMP');
%Ori_B = imread('clockB_tif2bmp (256x256).BMP');
ori_A = imread('clockA.BMP');
ori_B = imread('clockB.BMP');
ori_A=rgb2gray(ori_A);
ori_B=rgb2gray(ori_B);
A_a = double(ori_A)/255;
B_b = double(ori_B)/255;
[p1,q1] = size(ori_A);
[p2,q2] = size(ori_B);
if(p1~=p2||q1~=q2)
    error('The size of image A and image B must be the same......');
end
F = fuse_lap(A_a,B_b,4,1,3);
%%
F=F*255;
F(F<0)=0;
disp('拉普拉斯指标:')
F(F>255)=255;
F=round(F);
%%
grey_level=256;
Criteria=Evaluation(double(ori_A),double(ori_B),F,grey_level);
disp(Criteria)
figure,imshow(F-double(ori_A),[])
F_DWT = fuse_dwb(A_a,B_b,4,1,3);
F_DWT=F_DWT*255;
F_DWT(F_DWT<0)=0;
disp('离散小波变换指标:')
F_DWT(F_DWT>255)=255;
F_DWT=round(F_DWT);
%%
grey_level=256;
Criteria=Evaluation(double(ori_A),double(ori_B),F_DWT,grey_level);
disp(Criteria)
figure,imshow(F_DWT-double(ori_A),[])
F_CON = fuse_con(A_a,B_b,4,1,3);


F_CON=F_CON*255;
F_CON(F_CON<0)=0;
disp('Contrast Pyramid变换指标:')
F_CON(F_CON>255)=255;
F_CON=round(F_CON);
%%
grey_level=256;
figure,imshow(F_CON-double(ori_A),[])
Criteria=Evaluation(double(ori_A),double(ori_B),F_CON,grey_level);
disp(Criteria)

F_SIH = fuse_sih(A_a,B_b,4,1,3);

F_SIH=F_SIH*255;
F_SIH(F_SIH<0)=0;
disp('FSD Pyramid变换指标:')
F_SIH(F_SIH>255)=255;
F_SIH=round(F_SIH);
%%
grey_level=256;
figure,imshow(F_SIH-double(ori_A),[])
Criteria=Evaluation(double(ori_A),double(ori_B),F_SIH,grey_level);
disp(Criteria)