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)