We present a method to automatically register images presenting both global and local deformations. The image registration process is performed by exploiting a multi-level/multi-image approach whereby after having wavelet-transformed the images, the subband images at different levels are used in a non-feature-based way to determine the motion vectors between the reference and the target images. The crude motion field determined by block matching at the coarsest level of the pyramid is successively refined by taking advantage of both the orientation sensitivity of the different subbands and the contribution of the adjacent levels. The final registered image is obtained by applying the motion field to the lowest level of the pyramid and by inversely transforming it.