Starting from the sub-pixel de-registration information at each GCP, a deformation model is sought, describing de-registration at every pixel. The image is tessellated by means of Delaunay triangulation into polygons containing N GCP. For each of these polygons the coefficients of a 2D polynomial are estimated. These polynomials describe de-registration within the respective polygon.