1. A method for deformable registration of 2 digital images, comprising the steps of: providing a pair of digital images, each image comprising a plurality of intensities associated with a 2-dimensional grid of points, said pair including a fixed image and a moving image; extracting a set of edge images from each image of said pair of images, each edge set being extracted at a different resolution; selecting a pair of edge images with a lowest resolution; determining a mapping from edge points of said fixed image to edge points of moving image using a geodesic thin plate spline interpolation; applying said mapping to a next higher resolution edge point image of said moving image; selecting a pair of edge images at a next higher resolution, wherein a moving edge image is the moving edge image to which said mapping has been applied; repeating said steps of determining a mapping, applying said mapping, and selecting a pair of edge images at a next higher resolution for all edge images in said set of edge images; and applying said mapping to an entire moving image.