Intensity-based Non-rigid Registration Algorithm Research Articles

Automatic selection of landmarks in T1-weighted head MRI with regression forests for image registration initialization.

Medical image registration establishes a correspondence between images of biological structures and it is at the core of many applications. Commonly used deformable image registration methods are dependent on a good preregistration initialization. The initialization can be performed by localizing homologous landmarks and calculating a point-based transformation between the images. The selection of landmarks is however important. In this work, we present a learning-based method to automatically find a set of robust landmarks in 3D MR image volumes of the head to initialize non-rigid transformations. To validate our method, these selected landmarks are localized in unknown image volumes and they are used to compute a smoothing thin-plate splines transformation that registers the atlas to the volumes. The transformed atlas image is then used as the preregistration initialization of an intensity-based non-rigid registration algorithm. We show that the registration accuracy of this algorithm is statistically significantly improved when using the presented registration initialization over a standard intensity-based affine registration.

Ultrasound 2D Strain Estimator Based on Image Registration for Ultrasound Elastography.

In this paper, we present a new approach to calculate 2D strain through the registration of the pre- and post-compression (deformation) B-mode image sequences based on an intensity-based non-rigid registration algorithm (INRA). Compared with the most commonly used cross-correlation (CC) method, our approach is not constrained to any particular set of directions, and can overcome displacement estimation errors introduced by incoherent motion and variations in the signal under high compression. This INRA method was tested using phantom and in vivo data. The robustness of our approach was demonstrated in the axial direction as well as the lateral direction where the standard CC method frequently fails. In addition, our approach copes well under large compression (over 6%). In the phantom study, we computed the strain image under various compressions and calculated the signal-to-noise (SNR) and contrast-to-noise (CNS) ratios. The SNR and CNS values of the INRA method were much higher than those calculated from the CC-based method. Furthermore, the clinical feasibility of our approach was demonstrated with the in vivo data from patients with arm lymphedema.

MRI signal intensity based B‐Spline nonrigid registration for pre‐ and intraoperative imaging during prostate brachytherapy

To apply an intensity-based nonrigid registration algorithm to MRI-guided prostate brachytherapy clinical data and to assess its accuracy. A nonrigid registration of preoperative MRI to intraoperative MRI images was carried out in 16 cases using a Basis-Spline algorithm in a retrospective manner. The registration was assessed qualitatively by experts' visual inspection and quantitatively by measuring the Dice similarity coefficient (DSC) for total gland (TG), central gland (CG), and peripheral zone (PZ), the mutual information (MI) metric, and the fiducial registration error (FRE) between corresponding anatomical landmarks for both the nonrigid and a rigid registration method. All 16 cases were successfully registered in less than 5 min. After the nonrigid registration, DSC values for TG, CG, PZ were 0.91, 0.89, 0.79, respectively, the MI metric was -0.19 +/- 0.07 and FRE presented a value of 2.3 +/- 1.8 mm. All the metrics were significantly better than in the case of rigid registration, as determined by one-sided t-tests. The intensity-based nonrigid registration method using clinical data was demonstrated to be feasible and showed statistically improved metrics when compare to only rigid registration. The method is a valuable tool to integrate pre- and intraoperative images for brachytherapy.

Explicit Incorporation of Prior Anatomical Information Into a Nonrigid Registration of Thoracic and Abdominal CT and 18-FDG Whole-Body Emission PET Images

The aim of this paper is to develop a registration methodology in order to combine anatomical and functional information provided by thoracic/abdominal computed tomography (CT) and whole-body positron emission tomography (PET) images. The proposed procedure is based on the incorporation of prior anatomical information in an intensity-based nonrigid registration algorithm. This incorporation is achieved in an explicit way, initializing the intensity-based registration stage with the solution obtained by a nonrigid registration of corresponding anatomical structures. A segmentation algorithm based on a hierarchically ordered set of anatomy-specific rules is used to obtain anatomical structures in CT and emission PET scans. Nonrigid deformations are modeled in both registration stages by means of free-form deformations, the optimization of the control points being achieved by means of an original vector field-based approach instead of the classical gradient-based techniques, considerably reducing the computational time of the structure registration stage. We have applied the proposed methodology to 38 sets of images (33 provided by standalone machines and five by hybrid systems) and an assessment protocol has been developed to furnish a qualitative evaluation of the algorithm performance.

Three-dimensional average-shape atlas of the honeybee brain and its applications

The anatomical substrates of neural nets are usually composed from reconstructions of neurons that were stained in different preparations. Realistic models of the structural relationships between neurons require a common framework. Here we present 3-D reconstructions of single projection neurons (PN) connecting the antennal lobe (AL) with the mushroom body (MB) and lateral horn, groups of intrinsic mushroom body neurons (type 5 Kenyon cells), and a single mushroom body extrinsic neuron (PE1), aiming to compose components of the olfactory pathway in the honeybee. To do so, we constructed a digital standard atlas of the bee brain. The standard atlas was created as an average-shape atlas of 22 neuropils, calculated from 20 individual immunostained whole-mount bee brains. After correction for global size and positioning differences by repeatedly applying an intensity-based nonrigid registration algorithm, a sequence of average label images was created. The results were qualitatively evaluated by generating average gray-value images corresponding to the average label images and judging the level of detail within the labeled regions. We found that the first affine registration step in the sequence results in a blurred image because of considerable local shape differences. However, already the first nonrigid iteration in the sequence corrected for most of the shape differences among individuals, resulting in images rich in internal detail. A second iteration improved on that somewhat and was selected as the standard. Registering neurons from different preparations into the standard atlas reveals 1) that the m-ACT neuron occupies the entire glomerulus (cortex and core) and overlaps with a local interneuron in the cortical layer; 2) that, in the MB calyces and the lateral horn of the protocerebral lobe, the axon terminals of two identified m-ACT neurons arborize in separate but close areas of the neuropil; and 3) that MB-intrinsic clawed Kenyon cells (type 5), with somata outside the calycal cups, project to the peduncle and lobe output system of the MB and contact (proximate) the dendritic tree of the PE1 neuron at the base of the vertical lobe. Thus the standard atlas and the procedures applied for registration serve the function of creating realistic neuroanatomical models of parts of a neural net. The Honeybee Standard Brain is accessible at www.neurobiologie.fu-berlin.de/beebrain.

Reduction of motion artifacts in three-dimensional CT-DSA using constrained adaptive multilevel free-form registration

In this paper, we describe the application of an intensity-based nonrigid registration algorithm to (3-D) CT-digital subtraction angiography (CT-DSA). On a uniform control point grid (CPG), a B-spline deformation is computed, which maximizes the normalized mutual information (NMI) similarity measure between the 3-D pre- and post-contrast CT images. This transformation is then applied when computing subtraction images. We demonstrate that this algorithm can measure nonrigid patient motion between the image acquisitions and substantially reduces motion artifacts relative to results obtained with rigid-body registration.