%????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
?????