非常好用的融合图像评价方法

%    objective Image Fusion Performance Measure
%    C.S.Xydeas and V.Petrovic
%    Electronics Letters

%    edit by Richang Hong; May 27 07

function output = Qabf（strA strB strF）
% strA and strB is the source images and strF is the fusion result

% model parameters
    L=1; Tg=0.9994;kg=-15;Dg=0.5;Ta=0.9879;ka=-22;Da=0.8;    
    
% Sobel Operator
h1=[1 2 1;0 0 0;-1 -2 -1]; h2=[0 1 2;-1 0 1;-2 -1 0]; h3=[-1 0 1;-2 0 2;-1 0 1];
% if y is the response to h1 and x is the response to h3;
% then the intensity is sqrt（x^2+y^2） and orientation is arctan（y/x）;
pA = imread（strA）; pA = double（pA）;
pB = imread（strB）; pB = double（pB）;
pF = imread（strF）; pF = double（pF）;

SAx = conv2（pAh3‘same‘）; SAy = conv2（pAh1‘same‘）;
gA = sqrt（SAx.^2 + SAy.^2）; 
[MN] = size（SAx）; aA = zeros（MN）;
for i=1:M
    for j=1:N
        if （ SAx（ij） == 0 ） aA（ij） = pi/2;
        else
            aA（ij） = atan（SAy（ij）/SAx（ij））;
        end
    end
end

SBx = conv2（pBh3‘same‘）; SBy = conv2（pBh1‘same‘）;
gB = sqrt（SBx.^2 + SBy.^2）; 
[MN] = size（SBx）; aB = zeros（MN）;
for i=1:M
    for j=1:N
        if （ SBx（ij） == 0 ） aB（ij） = pi/2;
        else
            aB（ij） = atan（SBy（ij）/SBx（ij））;
        end
    end
end

SFx = conv2（pFh3‘same‘）; SFy = conv2（pFh1‘same‘）;
gF = sqrt（SFx.^2 + SFy.^2）; 
[MN] = size（SAx）; aF = zeros（MN）;
for i=1:M
    for j=1:N
        if （ SFx（ij） == 0 ） aF（ij） = pi/2;
        else