1. A method for registering a source image with a target image using multiple scales, wherein the images are deformable, the method comprising: measuring dissimilarity between the source image and the target image; using the measured dissimilarity for each scale of the multiple scales to minimize a discrete energy function, wherein the discrete energy function for each scale contains a data term and a smoothness term, the data term represents deformation vectors as discrete labels, and the smoothness term imposes prior probabilities on the deformation to preserve discontinuities, and wherein prior to warping the source image, each scale includes: applying multi-scale Markov random field registration to determine a labeling which corresponds to a deformation vector field; and accumulating the deformation vector field via deformation vector addition to an accumulated deformation vector field determined at a previous scale, wherein the deformation vectors are added on the basis that the deformation vector field at a final scale corresponds to the deformation vector field at a desired resolution, such that the final deformation vectors are aggregated at each scale prior to warping; warping the source image, according to the accumulated deformation vector field at the desired resolution, to obtain the target image, wherein the steps are performed in a processor.