depnormcorr.m

Lukas Kanade Multi resolution program in matlab language.

function [dx,dy]=depnormcorr(I1,I2);

[L,C]=size(I1);

indMin=(L/4)+1;
indMax=(3*L/4);
M=I1(indMin:indMax,indMin:indMax);
F = normxcorr2(M,I2);

% max de correlation au pixel pres
[valMax,xmax]=max(max(F));
[valMax,ymax]=max(max(F'));

% max de correlation subpixel
 vois=3;
 Profx=F(ymax,(xmax-vois):(xmax+vois));
 Profy=F((ymax-vois):(ymax+vois),xmax)';
  
 [Px,Sx] = polyfit((xmax-vois):(xmax+vois),Profx,2);
 [Py,Sy] = polyfit((ymax-vois):(ymax+vois),Profy,2);
 
 xmaxs=-Px(2)/2/Px(1);
 ymaxs=-Py(2)/2/Py(1);
 
 [LF,CF]=size(F);
 dx=xmaxs-(LF/2)-0.5;
 dy=ymaxs-(CF/2)-0.5;