A method for elastic registration of a two-dimensional source digital image of an object of interest with a slice of a three dimensional target volume of the object of interest, the method being characterized in that it comprises defining a Markov Random Field framework comprising at least one undirected pairwise graph superimposed on the two-dimensional image domain comprising one or two sets of regular vertices and one, two or three sets of edges, defining at least a grid of control points ( p i), each control point corresponding to either a vertex of a set of vertices, or to two vertices each of a set of vertices, and a neighborhood system of edges associated with vertices, further comprising: - defining a set of multi-dimensional labels of a discrete space to associate to each control point ( p i), wherein a spatially displaced control point ( p' i) is obtained by assigning a given label to said control point, and wherein a plane transformation adapted to obtain a plane slice of the target volume and an in-plane deformation transformation of the source digital image are determined simultaneously by said labels - obtaining a set of labels defining a set of displaced control points ( p' i), by minimizing an energy target function, the set of displaced control points defining a plane slice (Π i) of the target volume that elastically matches the two-dimensional source digital image according to the energy target function, - the energy target function being expressed as a sum of terms comprising: at least one first order term computed as a sum, over the control points, of unary potentials depending on the vertex multi-dimensional label, and/or at least one second order term computed as a sum, over the pairs of control points according to a second order neighborhood system, of pairwise potentials depending on pairs of multi-dimensional labels of two connected vertices, and/or at least one higher order term, over higher order cliques of control points according to a higher order neighborhood system, of higher order potentials depending on a number of labels depending on the order of said cliques. - defining a set of multi-dimensional labels of a discrete space to associate to each control point ( p i), wherein a spatially displaced control point ( p' i) is obtained by assigning a given label to said control point, and wherein a plane transformation adapted to obtain a plane slice of the target volume and an in-plane deformation transformation of the source digital image are determined simultaneously by said labels - obtaining a set of labels defining a set of displaced control points ( p' i), by minimizing an energy target function, the set of displaced control points defining a plane slice (Π i) of the target volume that elastically matches the two-dimensional source digital image according to the energy target function, - the energy target function being expressed as a sum of terms comprising: at least one first order term computed as a sum, over the control points, of unary potentials depending on the vertex multi-dimensional label, and/or at least one second order term computed as a sum, over the pairs of control points according to a second order neighborhood system, of pairwise potentials depending on pairs of multi-dimensional labels of two connected vertices, and/or at least one higher order term, over higher order cliques of control points according to a higher order neighborhood system, of higher order potentials depending on a number of labels depending on the order of said cliques. at least one first order term computed as a sum, over the control points, of unary potentials depending on the vertex multi-dimensional label, and/or at least one second order term computed as a sum, over the pairs of control points according to a second order neighborhood system, of pairwise potentials depending on pairs of multi-dimensional labels of two connected vertices, and/or at least one higher order term, over higher order cliques of control points according to a higher order neighborhood system, of higher order potentials depending on a number of labels depending on the order of said cliques.