Abstract—An overview is presented of the medical image processing literature on mutual-information-based registration. The aim of the survey is threefold: an introduction for those new to the field, an overview for those working in the field, and a reference for those searching for literature on a specific application. Methods are classified according to the different aspects of mutual-information-based registration. The main division is in aspects of the methodology and of the application. The part on methodology describes choices made on facets such as preprocessing of images, gray value interpolation, optimization, adaptations to the mutual information measure, and different types of geometrical transformations. The part on applications is a reference of the literature available on different modalities, on interpatient registration and on different anatomical objects. Comparison studies including mutual information are also considered. The paper starts with a description of entropy and mutual information and it closes with a discussion on past achievements and some future challenges. Index Terms—Image registration, literature survey, matching, mutual information. I.

One major problem with nonrigid image registration techniques is their high computational cost. Because of this, these methods have found limited application to clinical situations where fast execution is required, e.g., intraoperative imaging. This paper presents a parallel implementation of a nonrigid image registration algorithm. It takes advantage of shared-memory multiprocessor computer architectures using multithreaded programming by partitioning of data and partitioning of tasks, depending on the computational subproblem. For three different biomedical applications (intraoperative brain deformation, contrast-enhanced MR mammography, intersubject brain registration), the scaling behavior of the algorithm is quantitatively analyzed. The method is demonstrated to perform the computation of intra-operative brain deformation in less than a minute using 64 CPUs on a 128-CPU shared-memory supercomputer (SGI Origin 3800). It is shown that its serial component is no more than 2% of the total computation time, allowing a speedup of at least a factor of 50. In most cases, the theoretical limit of the speedup is substantially higher (up to 132-fold in the application examples presented in this paper). The parallel implementation of our algorithm is, therefore, capable of solving nonrigid registration problems with short execution time requirements and may be considered an important step in the application of such techniques to clinically important problems such as the computation of brain deformation during cranial imageguided surgery.

In this paper, we extend a previously reported intensity-based nonrigid registration algorithm by using a novel regularization term to constrain the deformation. Global motion is modeled by a rigid transformation while local motion is described by a free-form deformation based on B-splines. An information theoretic measure, normalized mutual information, is used as an intensity-based image similarity measure. Registration is performed by searching for the deformation that minimizes a cost function consisting of a weighted combination of the image similarity measure and a regularization term. The novel regularization term is a local volume-preservation (incompressibility) constraint, which is motivated by the assumption that soft tissue is incompressible for small deformations and short time periods. The incompressibility constraint is implemented by penalizing deviations of the Jacobian determinant of the deformation from unity. We apply the nonrigid registration algorithm with and without the incompressibility constraint to precontrast and postcontrast magnetic resonance (MR) breast images from 17 patients. Without using a constraint, the volume of contrast-enhancing lesions decreases by 1%--78% (mean 26%). Image improvement (motion artifact reduction) obtained using the new constraint is compared with that obtained using a smoothness constraint based on the bending energy of the coordinate grid by blinded visual assessment of maximum intensity projections of subtraction images. For both constraints, volume preservation improves, and motion artifact correction worsens, as the weight of the constraint penalty term increases. For a given volume change of the contrast-enhancing lesions (2% of the original volume), the incompressibility constraint reduces motion artifacts ...

Image registration refers to the process of finding the spatial correspondence between two or more images. This is usually done by applying a spatial transformation, computed automatic or manually, to a given image using a continuous image model computed either with interpolation or approximation methods. We show that noise induced signal variance in interpolated images differs significantly from the signal variance of the original images in native space. We describe a simple approach to compute the signal variance in registered images based on the signal variance and covariance of the original images, the spatial transformations computed by the registration procedure, and the interpolation or approximation kernel chosen. Our approach is applied to diffusion tensor (DT) MRI data. We show that incorrect noise variance estimates in registered diffusion weighted images can affect the estimated DT parameters, their estimated uncertainty, as well as indices of goodness of fit such as chi-square maps. In addition to DT-MRI, we believe that this methodology would be useful any time parameter extraction methods are applied to registered or interpolated data.

Abstract. Diffusion weighted images (DWIs) are commonly acquired with Echo-planar imaging (EPI). B0 inhomogeneities affect EPI by producing spatially nonlinear image distortions. Several strategies have been proposed to correct EPI distortions including B 0 field mapping (B 0M) and image registration. In this study, an experimental framework is proposed to evaluation the performance of different EPI distortion correction methods in improving DTderived quantities. A deformable registration based method with mutual information metric and cubic B-spline modeled constrained deformation field (BSP) is proposed as an alternative when B 0 mapping data are not available. BSP method is qualitatively and quantitatively compared to B 0M method using the framework. Both methods can successful reduce EPI distortions and significantly improve the quality of DT-derived quantities. Overall, B 0M was clearly superior in infratentorial regions including brainstem and cerebellum, as well as in the ventral areas of the temporal lobes while BSP was better in all rostral brain regions. Keywords: DTI, EPI distortion correction, B 0 field mapping, B-spline image registration. 1

Passband balanced-steady-state free precession (b-SSFP) fMRI is a recently developed method that utilizes the passband (flat portion) of the b-SSFP off-resonance response to measure MR signal changes elicited by changes in tissue oxygenation following increases in neuronal activity. Rapid refocusing and short readout durations of b-SSFP, combined with the relatively large flat portion of the b-SSFP off-resonance spectrum allows distortion-free full-brain coverage with only two acquisitions. This allows for high-resolution functional imaging, without the spatial distortion frequently encountered in conventional highresolution functional images. Finally, the 3D imaging compatibility of the b-SSFP acquisitions permits isotropic-voxel-size high-resolution acquisitions. In this study we address some of the major technical issues involved in obtaining passband b-SSFP-based functional brain images with practical imaging parameters and demonstrate the advantages through breathholding and visual field mapping experiments. Magn Reson

Estimating intensity variance due to noise in registered images:

Abstract. This paper presents a method for correcting the geometric and greyscale distortions in diffusion-weighted MRI that result from inhomogeneities in the static magnetic field. These inhomogeneities may due to imperfections in the magnet or to spatial variations in the magnetic susceptibility of the object being imaged—so called susceptibility artifacts. Echo-planar imaging (EPI), used in virtually all diffusion weighted acquisition protocols, assumes a homogeneous static field, which generally does not hold for head MRI. The resulting distortions are significant, sometimes more than ten millimeters. These artifacts impede accurate alignment of diffusion images with structural MRI, and are generally considered an obstacle to the joint analysis of connectivity and structure in head MRI. In principle, susceptibility artifacts can be corrected by acquiring (and applying) a field map. However, as shown in the literature and demonstrated in this paper, field map corrections of susceptibility artifacts are not entirely accurate and reliable, and thus field maps do not produce reliable alignment of EPIs with corresponding structural images. This paper presents a new, image-based method for correcting susceptibility artifacts. The method relies on a variational formulation of the match between an EPI baseline image and a corresponding T2-weighted structural image but also specifically accounts for the physics of susceptibility artifacts. We derive a set of partial differential equations associated with the optimization, describe the numerical methods for solving these equations, and present results that demonstrate the effectiveness of the proposed method compared with field-map correction. 1