motion estimation and compensation for very low bitrate video coding

D.2 Pseudo-code of the implementation 201
// normal case
else
{
(*current).feature = O;
// next point
(*current).next = new point;
(*(*current).next).contour = num;
(*(*current).next).vpos = posv;
(*(*current).next).hpos = posh;
(*(*current).next).next = NULL;
(*(*current).next).previous = current;
return track_forward(edged,type,typer,num,pDouble,pTraced,
(*current).next);
}
}
else
// the point is a gap
{
(*current).feature = G;
// no next point
(*current).next = NULL;
return G;
}
return 0;
}
/* Recursive tracking of half-boundaries in backwardd direction */
short track_backward(contour *edged,short type,short typer,short num,
picture pDouble, picture pTraced, Point current)
{
// similar to track_forward routine but for the
// (*current).previous point
}
D.2.3 Inverse Kriging System
/* Determine the system matrix A and the solution vector B
for the Inverse Kriging Interpolation, as exposed in Chapter 5 */
determine_system(wireframe wF, motion mField, double **A, double *B);
{
// Initialization
A = 0;
for (every block #i of the motion vector field)
{
/* Backwards motion field to be reversed
Determine origin of the motion vector according to the
center position (posv[i],posh[i]) of the block */
VertOrigin = MaxVector+posv[i]+vertical_component(mField,i);
HoriOrigin = MaxVector+posh[i]+horizontal_component(mField,i);
/* Determine into wich triangle, made out of vertices
(Top1V,Top1H), (Top2V,Top2H), and (Top3V,Top3H), is located
the origin of the motion vector */