makde lucasKanade Fcentre always
This requires principally changing computeDICoperators.cpp
around line 93
/* Variable to hold current coordinate (x_j(p)) in eq 12 */
int xjofp;
switch( jofp )
{
case 0: xjofp = z1; break;
case 1: xjofp = y1; break;
case 2: xjofp = x1; break;
case 3: xjofp = 1; break;
}
// later...
M[ p + (12 * m) ] += ( xjofm * gradim2iofm ) * ( xjofp * gradim2iofp );
So that z1, y1, and x1 are valued 0 in the middle of the image rather than the top corner.
This will make a number of things more consistent.
We can take the opportunity to apply Fpoint by default in the middle of images in applyTransformationOperator
Edited by Edward Andò