主题:求帮忙NC归一化函数
请高手帮忙看看是不是有错。
%计算nc(归一化相关系数)
function dNC = nc(ImageA,ImageB)
if (size(ImageA,1) ~= size(ImageB,1)) or (size(ImageA,2) ~= size(ImageB,2))
error('ImageA <> ImageB');
dNC = 0;
return ;
end
ImageA=double(ImageA);
ImageB=double(ImageB);
M = size(ImageA,1);
N = size(ImageA,2);
d1=0 ;
d2=0;
d3=0;
for i = 1:M
for j = 1:N
d1=d1+ImageA(i,j)*ImageB(i,j) ;
d2=d2+ImageA(i,j)*ImageA(i,j) ;
d3=d3+ImageB(i,j)*ImageB(i,j) ;
end
end
dNC=d1/(sqrt(d2)*sqrt(d3));
return
%计算nc(归一化相关系数)
function dNC = nc(ImageA,ImageB)
if (size(ImageA,1) ~= size(ImageB,1)) or (size(ImageA,2) ~= size(ImageB,2))
error('ImageA <> ImageB');
dNC = 0;
return ;
end
ImageA=double(ImageA);
ImageB=double(ImageB);
M = size(ImageA,1);
N = size(ImageA,2);
d1=0 ;
d2=0;
d3=0;
for i = 1:M
for j = 1:N
d1=d1+ImageA(i,j)*ImageB(i,j) ;
d2=d2+ImageA(i,j)*ImageA(i,j) ;
d3=d3+ImageB(i,j)*ImageB(i,j) ;
end
end
dNC=d1/(sqrt(d2)*sqrt(d3));
return