求帮!立体匹配中NCC算法程序报错
程序如下:
function res=youhuaNCC(a1,r1,c1,a2,r2,c2) %a1为左灰度图像 r1,c1为左图像角点横纵坐标;
%a2为右灰度图像,r2,c2为右图像角点坐标
N=2; %窗口大小
[ligne1,col1]= size(c1); %源图像特征点数
%左图窗口内平均灰度值
realpt1=1;
tic
for i=1:ligne1
somme = 0;
ZNCC1 = 0;
num=0;
denom=0;
denom1=0;
denom2=0;
moyenne1=0;
for l=-N:1:N
for k=-N:1:N
somme=somme+double(a1(r1(i)+l,c1(i)+k));%求a1\2N+1 2N+1串口的灰度值的和
end
end
moyenne1=somme/(2N+1)^2; %a1\单个特征点窗口灰度值的均值T~X