We present in this paper a new registration and gain correction algorithm for 3D medical images. It is intensity based. The basic idea is to represent images by 4D points x j ; y j ; z j ; i j and to de ne a global energy function based on this representation. For minimisation, we propose a technique which does not require computing the derivatives of this criterion with respect to the parameters. It can be understood as an extension of the Iterative Closest Point algorithm 5, 56 or as an

more »

... ation of the formalism proposed in 13 . Two parameters enable us to develop a coarse-to-ne strategy both for resolution and for deformation. Our technique presents the advantage of minimising a well-de ned global criterion, to deal with various classes of transformations for example rigid, a ne, volume spline and radial basis functions, to be simple to implement, and to be e cient i n practice. Results on real brain and heart 3D images are presented to demonstrate the validity of our approach. We also explain how one can compute basic statistics on the deformation parameters to constrain the set of possible deformations by learning and to discriminate between di erent groups. 2