Brace indexing is not supported for variables of this type.

1 次查看(过去 30 天)
Error in objective_fusion_perform_fn (line 38)
[gvA,ghA]=sobel_fn(x{1});
function [gv gh]=sobel_fn(x)
%%% sobel_fn: Computes the vertical & horizontal edges of x using sobel operator.
%%%
%%% [gv gh]=sobel_fn(x)
%%%
%%% Author : B. K. SHREYAMSHA KUMAR
%%% Created on 28-10-2011.
%%% Updated on 28-10-2011.
vtemp=[-1 0 1;-2 0 2;-1 0 1]/8;
htemp=[-1 -2 -1;-0 0 0;1 2 1]/8;
[a b]=size(htemp);
x_ext=per_extn_im_fn(x,a);
[p,q]=size(x_ext);
for ii=2:p-1
for jj=2:q-1
gv(ii-1,jj-1)=sum(sum(x_ext(ii-1:ii+1,jj-1:jj+1).*vtemp));
gh(ii-1,jj-1)=sum(sum(x_ext(ii-1:ii+1,jj-1:jj+1).*htemp));
end
end
function [QABF,LABF,NABF,NABF1]=objective_fusion_perform_fn(xrcw,x)
%%% objective_fusion_perform_fn: Computes the Objective Fusion Performance Parameters proposed by Petrovic
%%% and modified Fusion Artifacts (NABF) measure proposed by B. K. Shreyamsha Kumar
%%%
%%% Inputs:
%%% xrcw -> fused image
%%% x -> source images, x{1}, x{2}
%%%
%%% Outputs:
%%% QABF -> Total information transferred from source images to fused image measure proposed by Petrovic
%%% LABF -> Total loss of information measure proposed by Petrovic
%%% NABF1 -> Fusion Artifacts measure proposed by Petrovic
%%% NABF -> Modified Fusion Artifacts measure proposed by B. K. Shreyamsha Kumar
%%%
%%% Author : B. K. SHREYAMSHA KUMAR
%%% Created on 28-10-2011.
%%% Updated on 08-11-2011.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Petrovic Metrics %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Parameters for Petrovic Metrics Computation.
Td=2;
wt_min=0.001;
P=1;
Lg=1.5;
Nrg=0.9999;
kg=19;
sigmag=0.5;
Nra=0.9995;
ka=22;
sigmaa=0.5;
%%% Edge Strength & Orientation.
[gvA,ghA]=sobel_fn(x{1});
gA=sqrt(ghA.^2+gvA.^2);
[gvB,ghB]=sobel_fn(x{2});
gB=sqrt(ghB.^2+gvB.^2);
[gvF,ghF]=sobel_fn(xrcw);
gF=sqrt(ghF.^2+gvF.^2);
%%% Relative Edge Strength & Orientation.
[p,q]=size(xrcw);
for ii=1:p
for jj=1:q
if(gA(ii,jj)==0 | gF(ii,jj)==0)
gAF(ii,jj)=0;
elseif(gA(ii,jj)>gF(ii,jj))
gAF(ii,jj)=gF(ii,jj)/gA(ii,jj);
else
gAF(ii,jj)=gA(ii,jj)/gF(ii,jj);
end
if(gB(ii,jj)==0 | gF(ii,jj)==0)
gBF(ii,jj)=0;
elseif(gB(ii,jj)>gF(ii,jj))
gBF(ii,jj)=gF(ii,jj)/gB(ii,jj);
else
gBF(ii,jj)=gB(ii,jj)/gF(ii,jj);
end
if(gvA(ii,jj)==0 & ghA(ii,jj)==0)
aA(ii,jj)=0;
else
aA(ii,jj)=atan(gvA(ii,jj)/ghA(ii,jj));
end
if(gvB(ii,jj)==0 & ghB(ii,jj)==0)
aB(ii,jj)=0;
else
aB(ii,jj)=atan(gvB(ii,jj)/ghB(ii,jj));
end
if(gvF(ii,jj)==0 & ghF(ii,jj)==0)
aF(ii,jj)=0;
else
aF(ii,jj)=atan(gvF(ii,jj)/ghF(ii,jj));
end
end
end
aAF=abs(abs(aA-aF)-pi/2)*2/pi;
aBF=abs(abs(aB-aF)-pi/2)*2/pi;
%%% Edge Preservation Coefficient.
QgAF=Nrg./(1+exp(-kg*(gAF-sigmag)));
QaAF=Nra./(1+exp(-ka*(aAF-sigmaa)));
QAF=sqrt(QgAF.*QaAF);
QgBF=Nrg./(1+exp(-kg*(gBF-sigmag)));
QaBF=Nra./(1+exp(-ka*(aBF-sigmaa)));
QBF=sqrt(QgBF.*QaBF);
%%% Total Fusion Performance (QABF).
wtA=wt_min*ones(p,q);
wtB=wt_min*ones(p,q);
cA=ones(p,q); cB=ones(p,q);
for ii=1:p
for jj=1:q
if(gA(ii,jj)>=Td)
wtA(ii,jj)=cA(ii,jj)*gA(ii,jj)^Lg;
end
if(gB(ii,jj)>=Td)
wtB(ii,jj)=cB(ii,jj)*gB(ii,jj)^Lg;
end
end
end
wt_sum=sum(sum(wtA+wtB));
QAF_wtsum=sum(sum(QAF.*wtA))/wt_sum; %% Information Contributions of A.
QBF_wtsum=sum(sum(QBF.*wtB))/wt_sum; %% Information Contributions of B.
QABF=QAF_wtsum+QBF_wtsum %% QABF=sum(sum(QAF.*wtA+QBF.*wtB))/wt_sum -> Total Fusion Performance.
%%% Fusion Gain (QdeltaABF).
Qdelta=abs(QAF-QBF);
QCinfo=(QAF+QBF-Qdelta)/2;
QdeltaAF=QAF-QCinfo;
QdeltaBF=QBF-QCinfo;
QdeltaAF_wtsum=sum(sum(QdeltaAF.*wtA))/wt_sum;
QdeltaBF_wtsum=sum(sum(QdeltaBF.*wtB))/wt_sum;
QdeltaABF=QdeltaAF_wtsum+QdeltaBF_wtsum; %% Total Fusion Gain.
QCinfo_wtsum=sum(sum(QCinfo.*(wtA+wtB)))/wt_sum;
QABF11=QdeltaABF+QCinfo_wtsum; %% Total Fusion Performance.
%%% Fusion Loss (LABF).
rr=zeros(p,q);
for ii=1:p
for jj=1:q
if(gF(ii,jj)<=gA(ii,jj) | gF(ii,jj)<=gB(ii,jj))
rr(ii,jj)=1;
else
rr(ii,jj)=0;
end
end
end
LABF=sum(sum(rr.*((1-QAF).*wtA+(1-QBF).*wtB)))/wt_sum
%%% Fusion Artifacts (NABF) by Petrovic.
for ii=1:p
for jj=1:q
if(gF(ii,jj)>gA(ii,jj) & gF(ii,jj)>gB(ii,jj))
na1(ii,jj)=2-QAF(ii,jj)-QBF(ii,jj);
else
na1(ii,jj)=0;
end
end
end
NABF1=sum(sum(na1.*(wtA+wtB)))/wt_sum
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Fusion Artifacts (NABF) changed by B. K. Shreyamsha Kumar .
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for ii=1:p
for jj=1:q
if(gF(ii,jj)>gA(ii,jj) & gF(ii,jj)>gB(ii,jj))
na(ii,jj)=1;
else
na(ii,jj)=0;
end
end
end
NABF=sum(sum(na.*((1-QAF).*wtA+(1-QBF).*wtB)))/wt_sum
QABF+LABF+NABF1
QABF+LABF+NABF

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Image Analysis 的更多信息

标签

产品


版本

R2018a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by