DT-REFinD: diffusion tensor registration with exact finite-strain differential

MMORF-FSL's MultiMOdal Registration Framework

We present MMORF-FSL's MultiMOdal Registration Framework-a newly released nonlinear image registration tool designed primarily for application to magnetic resonance imaging (MRI) images of the brain. MMORF is capable of simultaneously optimising both displacement and rotational transformations within a single registration framework by leveraging rich information from multiple scalar and tensor modalities. The regularisation employed in MMORF promotes local rigidity in the deformation, and we have previously demonstrated how this effectively controls both shape and size distortion, leading to more biologically plausible warps. The performance of MMORF is benchmarked against three established nonlinear registration methods-FNIRT, ANTs, and DR-TAMAS-across four domains: FreeSurfer label overlap, diffusion tensor imaging (DTI) similarity, task-fMRI cluster mass, and distortion. The evaluation is based on 100 unrelated subjects from the Human Connectome Project (HCP) dataset registered to the Oxford-MultiModal-1 (OMM-1) multimodal template via either the T1w contrast alone or in combination with a DTI/DTI-derived contrast. Results show that MMORF is the most consistently high-performing method across all domains-both in terms of accuracy and levels of distortion. MMORF is available as part of FSL, and its inputs and outputs are fully compatible with existing workflows. We believe that MMORF will be a valuable tool for the neuroimaging community, regardless of the domain of any downstream analysis, providing state-of-the-art registration performance that integrates into the rich and widely adopted suite of analysis tools in FSL.

Deep Diffusion MRI Registration (DDMReg): A Deep Learning Method for Diffusion MRI Registration

In this paper, we present a deep learning method, DDMReg, for accurate registration between diffusion MRI (dMRI) datasets. In dMRI registration, the goal is to spatially align brain anatomical structures while ensuring that local fiber orientations remain consistent with the underlying white matter fiber tract anatomy. DDMReg is a novel method that uses joint whole-brain and tract-specific information for dMRI registration. Based on the successful VoxelMorph framework for image registration, we propose a novel registration architecture that leverages not only whole brain information but also tract-specific fiber orientation information. DDMReg is an unsupervised method for deformable registration between pairs of dMRI datasets: it does not require nonlinearly pre-registered training data or the corresponding deformation fields as ground truth. We perform comparisons with four state-of-the-art registration methods on multiple independently acquired datasets from different populations (including teenagers, young and elderly adults) and different imaging protocols and scanners. We evaluate the registration performance by assessing the ability to align anatomically corresponding brain structures and ensure fiber spatial agreement between different subjects after registration. Experimental results show that DDMReg obtains significantly improved registration performance compared to the state-of-the-art methods. Importantly, we demonstrate successful generalization of DDMReg to dMRI data from different populations with varying ages and acquired using different acquisition protocols and different scanners.

Bayesian Fully Convolutional Networks for Brain Image Registration

The purpose of medical image registration is to find geometric transformations that align two medical images so that the corresponding voxels on two images are spatially consistent. Nonrigid medical image registration is a key step in medical image processing, such as image comparison, data fusion, target recognition, and pathological change analysis. Existing registration methods only consider registration accuracy but largely neglect the uncertainty of registration results. In this work, a method based on the Bayesian fully convolutional neural network is proposed for nonrigid medical image registration. The proposed method can generate a geometric uncertainty map to calculate the uncertainty of registration results. This uncertainty can be interpreted as a confidence interval, which is essential for judging whether the source data are abnormal. Moreover, the proposed method introduces group normalization, which is conducive to the network convergence of the Bayesian neural network. Some representative learning-based image registration methods are compared with the proposed method on different image datasets. Experimental results show that the registration accuracy of the proposed method is better than that of the methods, and its antifolding performance is comparable to that of fast image registration and VoxelMorph. Furthermore, the proposed method can evaluate the uncertainty of registration results.

A Patient-Specific 3Dt Coronary Artery Motion Modeling Method Using Hierarchical Deformation with Electrocardiogram. SENSORS (BASEL, SWITZERLAND) 2020;20:E5680. [PMID: 33027998 PMCID: PMC7582594 DOI: 10.3390/s20195680]

Cardiovascular-related diseases are one of the leading causes of death worldwide. An understanding of heart movement based on images plays a vital role in assisting postoperative procedures and processes. In particular, if shape information can be provided in real-time using electrocardiogram (ECG) signal information, the corresponding heart movement information can be used for cardiovascular analysis and imaging guides during surgery. In this paper, we propose a 3D+t cardiac coronary artery model which is rendered in real-time, according to the ECG signal, where hierarchical cage-based deformation modeling is used to generate the mesh deformation used during the procedure. We match the blood vessel's lumen obtained from the ECG-gated 3D+t CT angiography taken at multiple cardiac phases, in order to derive the optimal deformation. Splines for 3D deformation control points are used to continuously represent the obtained deformation in the multi-view, according to the ECG signal. To verify the proposed method, we compare the manually segmented lumen and the results of the proposed method for eight patients. The average distance and dice coefficient between the two models were 0.543 mm and 0.735, respectively. The required time for registration of the 3D coronary artery model was 23.53 s/model. The rendering speed to derive the model, after generating the 3D+t model, was faster than 120 FPS.

Multimodal MRI Template Creation in the Ring-Tailed Lemur and Rhesus Macaque

We present a multimodal registration algorithm for simultaneous alignment of datasets with both scalar and tensor MRI images. We employ a volumetric, cubic B-spline parametrised transformation model. Regularisation is based on the logarithm of the singular values of the local Jacobian and ensures diffeomorphic warps. Tensor registration takes reorientation into account during optimisation, through a finite-strain approximation of rotation due to the warp. The combination of scalar, tensor and regularisation cost functions allows us to optimise the deformations in terms of tissue matching, orientation matching and distortion minimisation simultaneously. We apply our method to creating multimodal T2 and DTI MRI brain templates of two small primates (the ring-tailed lemur and rhesus macaque) from high-quality, ex vivo, 0.5/0.6 mm isotropic data. The resulting templates are of very high quality across both modalities and species. Tissue contrast in the T2 channel is high indicating excellent tissue-boundary alignment. The DTI channel displays strong anisotropy in white matter, as well as consistent left/right orientation information even in relatively isotropic grey matter regions. Finally, we demonstrate where the multimodal templating approach overcomes anatomical inconsistencies introduced by unimodal only methods.

Image correction for diffusion tensor imaging of Rhesus monkey thoracic spinal cord

BACKGROUND

The relatively tiny spinal cord of non-human primate (NHP) causes increased challenge in diffusion tensor imaging (DTI) post-processing. This study aimed to establish a reliable correction strategy applied to clinical DTI images of NHP.

METHODS

Six normal and partial spinal cord injury (SCI) rhesus monkeys underwent 3T MR scanning. A correction strategy combining multiple iterations and non-rigid deformation was used for DTI image post-processing. Quantitative evaluations were then conducted to investigate effects of distortion correction.

RESULTS

After correction, longitudinal geometric distortion, global distortion, and residual distance errors were all significantly decreased (P < 0.05). Fractional anisotropy at the injured site was remarkably lower than that at the contralateral site (P = 0.0488) and was substantially lower than those at the adjacent superior (P = 0.0157) and inferior (P = 0.0128) areas at the same side.

CONCLUSIONS

Our image correction strategy can improve the quality of the DTI images of NHP thoracic cords, contributing to the development of SCI preclinical research.

Combining white matter diffusion and geometry for tract-specific alignment and variability analysis

We present a framework for along-tract analysis of white matter (WM) fiber bundles based on diffusion tensor imaging (DTI) and tractography. We introduce the novel concept of fiber-flux density for modeling fiber tracts' geometry, and combine it with diffusion-based measures to define vector descriptors called Fiber-Flux Diffusion Density (FFDD). The proposed model captures informative features of WM tracts at both the microscopic (diffusion-related) and macroscopic (geometry-related) scales, thus enabling improved sensitivity to subtle structural abnormalities that are not reflected by either diffusion or geometrical properties alone. A key step in this framework is the construction of an FFDD dissimilarity measure for sub-voxel alignment of fiber bundles, based on the fast marching method (FMM). The obtained aligned WM tracts enable meaningful inter-subject comparisons and group-wise statistical analysis. Moreover, we show that the FMM alignment can be generalized in a straight forward manner to a single-shot co-alignment of multiple fiber bundles. The proposed alignment technique is shown to outperform a well-established, commonly used DTI registration algorithm. We demonstrate the FFDD framework on the Human Connectome Project (HCP) diffusion MRI dataset, as well as on two different datasets of contact sports players. We test our method using longitudinal scans of a basketball player diagnosed with a traumatic brain injury, showing compatibility with structural MRI findings. We further perform a group study comparing mid- and post-season scans of 13 active football players exposed to repetitive head trauma, to 17 non-player control (NPC) subjects. Results reveal statistically significant FFDD differences (p-values<0.05) between the groups, as well as increased abnormalities over time at spatially-consistent locations within several major fiber tracts of football players.

Advances in computational and statistical diffusion MRI

Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole-brain connectivity information that describes the brain's wiring diagram and population-based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high-level overview of interest to diffusion MRI researchers, with a more in-depth treatment to illustrate selected computational advances.

Age-specific gray and white matter DTI atlas for human brain at 33, 36 and 39 postmenstrual weeks

During the 3rd trimester, dramatic structural changes take place in the human brain, underlying the neural circuit formation. The survival rate of premature infants has increased significantly in recent years. The large morphological differences of the preterm brain at 33 or 36 postmenstrual weeks (PMW) from the brain at 40PMW (full term) make it necessary to establish age-specific atlases for preterm brains. In this study, with high quality (1.5 × 1.5 × 1.6 mm3 imaging resolution) diffusion tensor imaging (DTI) data obtained from 84 healthy preterm and term-born neonates, we established age-specific preterm and term-born brain templates and atlases at 33, 36 and 39PMW. Age-specific DTI templates include a single-subject template, a population-averaged template with linear transformation and a population-averaged template with nonlinear transformation. Each of the age-specific DTI atlases includes comprehensive labeling of 126 major gray matter (GM) and white matter (WM) structures, specifically 52 cerebral cortical structures, 40 cerebral WM structures, 22 brainstem and cerebellar structures and 12 subcortical GM structures. From 33 to 39 PMW, dramatic morphological changes of delineated individual neural structures such as ganglionic eminence and uncinate fasciculus were revealed. The evaluation based on measurements of Dice ratio and L1 error suggested reliable and reproducible automated labels from the age-matched atlases compared to labels from manual delineation. Applying these atlases to automatically and effectively delineate microstructural changes of major WM tracts during the 3rd trimester was demonstrated. The established age-specific DTI templates and atlases of 33, 36 and 39 PMW brains may be used for not only understanding normal functional and structural maturational processes but also detecting biomarkers of neural disorders in the preterm brains.

Multimodal image registration for liver radioembolization planning and patient assessment

Purpose

Multimodal imaging plays a key role in patient assessment and treatment planning in liver radioembolization. It will reach its full potential for convenient use in combination with deformable image registration methods. A registration framework is proposed for multimodal liver image registration of multi-phase CT, contrast-enhanced late-phase T1, T2, and DWI MRI sequences.

Methods

A chain of four pair-wise image registrations based on a variational registration framework using normalized gradient fields as distance measure and curvature regularization is introduced. A total of 103 cases of 35 patients was evaluated based on anatomical landmarks and deformation characteristics.

Results

Good anatomical correspondence and physical plausibility of the deformation fields were attained. The global mean landmark errors vary from 3.20 to 5.36 mm, strongly influenced by low resolved images in z-direction. Moderate volume changes are indicated by mean minimum and maximum Jacobian determinants of 0.44 up to 1.88. No deformation foldings were detected. The mean average divergence of the deformation fields range from 0.08 to 0.16 and the mean harmonic energies vary from 0.08 to 0.58.

Conclusion

The proposed registration solutions enable the combined use of information from multimodal imaging and provide an excellent basis for patient assessment and primary planning for liver radioembolization.

Comparison of spatial normalization strategies of diffusion MRI data for studying motor outcome in subacute-chronic and acute stroke

A common means of studying motor recovery in stroke patients is to extract Diffusion Tensor Imaging (DTI) parameters from the corticospinal tract (CST) and correlate them with clinical outcome scores. To that purpose, conducting group-level analyses through spatial normalization has become a popular approach. However, the reliability of such analyses depends on the accuracy of the particular registration strategy employed. To date, most studies have employed scalar-based registration using either high-resolution T1 images or Fractional Anisotropy (FA) maps to warp diffusion data to a common space. However, more powerful registration algorithms exist for aligning major white matter structures, such as Fiber Orientation Distribution (FOD)-based registration. Regardless of the strategy chosen, automatic normalization algorithms are prone to distortions caused by stroke lesions. While lesion masking is a common means to lessen such distortions, the extent of its effect on tract-related DTI parameters and their correlation with motor outcome has yet to be determined. Here, we aimed to address these concerns by first investigating the effect of common T1 and FA-based registration as well as novel FOD-based registration algorithms with and without lesion masking on lesion load and DTI parameter extraction of the CST in datasets typically acquired for subacute-chronic and acute stroke patients. Second, we studied how differences in these procedures influenced correlation strength between CST damage (through DTI parameters) and motor outcome. Our results showed that, for high-quality subacute-chronic stroke data, FOD-based registration captured significantly higher lesion loads and significantly larger FA asymmetries in the CST. This was also associated with significantly stronger correlations in motor outcome with respect to T1 or FA-based registration methods. For acute data acquired in a clinical setting, there were few observed differences, suggesting that commonly employed FA-based registration is appropriate for group-level analyses.

Medical Image Enhancement Method Based on the Fractional Order Derivative and the Directional Derivative

In recent years, the fractional order derivative has been introduced for image enhancement. It was proved that the medical image enhancement method based on the fractional order derivative has better effect than the method based on the integral order calculus. However, a priori information such as texture surrounding a pixel is normally ignored by the traditional fractional differential operators with the same value in the eight directions. To address the above problem, this paper presents a new medical image enhancement method by taking the merits of fractional differential and directional derivative. The proposed method considers the surrounding information (such as the image edge, clarity and texture information) and structural features of different pixels, as well as the directional derivative of each pixel in constructing the masks. By proposing this method, it can not only improve the high frequency information, but also improve the low frequency information of the image. Ultimately, it enhances the texture information of the image. Extensive experiments on four kinds of medical image demonstrate that the proposed algorithm is in favor of preserving more texture details and superior to the existing fractional differential algorithms on medical image enhancement.

Tractography in the clinics: Implementing a pipeline to characterize early brain development

In imaging studies of neonates, particularly in the clinical setting, diffusion tensor imaging-based tractography is typically unreliable due to the use of fast acquisition protocols that yield low resolution and signal-to-noise ratio (SNR). These image acquisition protocols are implemented with the aim of reducing motion artifacts that may be produced by the movement of the neonate's head during the scanning session. Furthermore, axons are not yet fully myelinated in these subjects. As a result, the water molecules' movements are not as constrained as in older brains, making it even harder to define structure using diffusion profiles. Here, we introduce a post-processing method that overcomes the difficulties described above, allowing the determination of reliable tracts in newborns. We tested our method using neonatal data and successfully extracted some of the limbic, association and commissural fibers, all of which are typically difficult to obtain by direct tractography. Geometrical and diffusion based features of the tracts are then utilized to compare premature babies to term babies. Our results quantify the maturation of white matter fiber tracts in neonates.

•

The proposed method enables consistent tractography in clinical datasets.

•

The tractography is used to structural positioning purposes

•

Geometrical features and diffusion variables in the tracts' paths are analyzed.

•

The gestational age was predicted with regressions in term and preterm babies.

•

The extracted features can be used as indexes of early neurodevelopment.

DTI Image Registration under Probabilistic Fiber Bundles Tractography Learning

Diffusion Tensor Imaging (DTI) image registration is an essential step for diffusion tensor image analysis. Most of the fiber bundle based registration algorithms use deterministic fiber tracking technique to get the white matter fiber bundles, which will be affected by the noise and volume. In order to overcome the above problem, we proposed a Diffusion Tensor Imaging image registration method under probabilistic fiber bundles tractography learning. Probabilistic tractography technique can more reasonably trace to the structure of the nerve fibers. The residual error estimation step in active sample selection learning is improved by modifying the residual error model using finite sample set. The calculated deformation field is then registered on the DTI images. The results of our proposed registration method are compared with 6 state-of-the-art DTI image registration methods under visualization and 3 quantitative evaluation standards. The experimental results show that our proposed method has a good comprehensive performance.

A Method for Automated Classification of Parkinson's Disease Diagnosis Using an Ensemble Average Propagator Template Brain Map Estimated from Diffusion MRI

Parkinson's disease (PD) is a common and debilitating neurodegenerative disorder that affects patients in all countries and of all nationalities. Magnetic resonance imaging (MRI) is currently one of the most widely used diagnostic imaging techniques utilized for detection of neurologic diseases. Changes in structural biomarkers will likely play an important future role in assessing progression of many neurological diseases inclusive of PD. In this paper, we derived structural biomarkers from diffusion MRI (dMRI), a structural modality that allows for non-invasive inference of neuronal fiber connectivity patterns. The structural biomarker we use is the ensemble average propagator (EAP), a probability density function fully characterizing the diffusion locally at a voxel level. To assess changes with respect to a normal anatomy, we construct an unbiased template brain map from the EAP fields of a control population. Use of an EAP captures both orientation and shape information of the diffusion process at each voxel in the dMRI data, and this feature can be a powerful representation to achieve enhanced PD brain mapping. This template brain map construction method is applicable to small animal models as well as to human brains. The differences between the control template brain map and novel patient data can then be assessed via a nonrigid warping algorithm that transforms the novel data into correspondence with the template brain map, thereby capturing the amount of elastic deformation needed to achieve this correspondence. We present the use of a manifold-valued feature called the Cauchy deformation tensor (CDT), which facilitates morphometric analysis and automated classification of a PD versus a control population. Finally, we present preliminary results of automated discrimination between a group of 22 controls and 46 PD patients using CDT. This method may be possibly applied to larger population sizes and other parkinsonian syndromes in the near future.

DR-TAMAS: Diffeomorphic Registration for Tensor Accurate Alignment of Anatomical Structures

In this work, we propose DR-TAMAS (Diffeomorphic Registration for Tensor Accurate alignMent of Anatomical Structures), a novel framework for intersubject registration of Diffusion Tensor Imaging (DTI) data sets. This framework is optimized for brain data and its main goal is to achieve an accurate alignment of all brain structures, including white matter (WM), gray matter (GM), and spaces containing cerebrospinal fluid (CSF). Currently most DTI-based spatial normalization algorithms emphasize alignment of anisotropic structures. While some diffusion-derived metrics, such as diffusion anisotropy and tensor eigenvector orientation, are highly informative for proper alignment of WM, other tensor metrics such as the trace or mean diffusivity (MD) are fundamental for a proper alignment of GM and CSF boundaries. Moreover, it is desirable to include information from structural MRI data, e.g., T1-weighted or T2-weighted images, which are usually available together with the diffusion data. The fundamental property of DR-TAMAS is to achieve global anatomical accuracy by incorporating in its cost function the most informative metrics locally. Another important feature of DR-TAMAS is a symmetric time-varying velocity-based transformation model, which enables it to account for potentially large anatomical variability in healthy subjects and patients. The performance of DR-TAMAS is evaluated with several data sets and compared with other widely-used diffeomorphic image registration techniques employing both full tensor information and/or DTI-derived scalar maps. Our results show that the proposed method has excellent overall performance in the entire brain, while being equivalent to the best existing methods in WM.

Free-breathing diffusion tensor imaging and tractography of the human heart in healthy volunteers using wavelet-based image fusion

Free-breathing cardiac diffusion tensor imaging (DTI) is a promising but challenging technique for the study of fiber structures of the human heart in vivo. This work proposes a clinically compatible and robust technique to provide three-dimensional (3-D) fiber architecture properties of the human heart. To this end, 10 short-axis slices were acquired across the entire heart using a multiple shifted trigger delay (TD) strategy under free breathing conditions. Interscan motion was first corrected automatically using a nonrigid registration method. Then, two post-processing schemes were optimized and compared: an algorithm based on principal component analysis (PCA) filtering and temporal maximum intensity projection (TMIP), and an algorithm that uses the wavelet-based image fusion (WIF) method. The two methods were applied to the registered diffusion-weighted (DW) images to cope with intrascan motion-induced signal loss. The tensor fields were finally calculated, from which fractional anisotropy (FA), mean diffusivity (MD), and 3-D fiber tracts were derived and compared. The results show that the comparison of the FA values (FA(PCATMIP) = 0.45 ±0.10, FA(WIF) = 0.42 ±0.05, P=0.06) showed no significant difference, while the MD values ( MD(PCATMIP)=0.83 ±0.12×10(-3) mm (2)/s, MD(WIF)=0.74±0.05×10(-3) mm (2)/s, P=0.028) were significantly different. Improved helix angle variations through the myocardium wall reflecting the rotation characteristic of cardiac fibers were observed with WIF. This study demonstrates that the combination of multiple shifted TD acquisitions and dedicated post-processing makes it feasible to retrieve in vivo cardiac tractographies from free-breathing DTI acquisitions. The substantial improvements were observed using the WIF method instead of the previously published PCATMIP technique.

Diffeomorphic metric mapping and probabilistic atlas generation of hybrid diffusion imaging based on BFOR signal basis

We first propose a large deformation diffeomorphic metric mapping algorithm to align multiple b-value diffusion weighted imaging (mDWI) data, specifically acquired via hybrid diffusion imaging (HYDI). We denote this algorithm as LDDMM-HYDI. We then propose a Bayesian probabilistic model for estimating the white matter atlas from HYDIs. We adopt the work given in Hosseinbor et al. (2013) and represent the q-space diffusion signal with the Bessel Fourier orientation reconstruction (BFOR) signal basis. The BFOR framework provides the representation of mDWI in the q-space and the analytic form of the emsemble average propagator (EAP) reconstruction, as well as reduces memory requirement. In addition, since the BFOR signal basis is orthonormal, the L(2) norm that quantifies the differences in the q-space signals of any two mDWI datasets can be easily computed as the sum of the squared differences in the BFOR expansion coefficients. In this work, we show that the reorientation of the q-space signal due to spatial transformation can be easily defined on the BFOR signal basis. We incorporate the BFOR signal basis into the LDDMM framework and derive the gradient descent algorithm for LDDMM-HYDI with explicit orientation optimization. Additionally, we extend the previous Bayesian atlas estimation framework for scalar-valued images to HYDIs and derive the expectation-maximization algorithm for solving the HYDI atlas estimation problem. Using real HYDI datasets, we show that the Bayesian model generates the white matter atlas with anatomical details. Moreover, we show that it is important to consider the variation of mDWI reorientation due to a small change in diffeomorphic transformation in the LDDMM-HYDI optimization and to incorporate the full information of HYDI for aligning mDWI. Finally, we show that the LDDMM-HYDI outperforms the LDDMM algorithm with diffusion tensors generated from each shell of HYDI.

Large deformation diffeomorphic registration of diffusion-weighted imaging data

Registration plays an important role in group analysis of diffusion-weighted imaging (DWI) data. It can be used to build a reference anatomy for investigating structural variation or tracking changes in white matter. Unlike traditional scalar image registration where spatial alignment is the only focus, registration of DWI data requires both spatial alignment of structures and reorientation of local signal profiles. As such, DWI registration is much more complex and challenging than scalar image registration. Although a variety of algorithms has been proposed to tackle the problem, most of them are restricted by the diffusion model used for registration, making it difficult to fit to the registered data a different model. In this paper we describe a method that allows any diffusion model to be fitted after registration for subsequent multifaceted analysis. This is achieved by directly aligning DWI data using a large deformation diffeomorphic registration framework. Our algorithm seeks the optimal coordinate mapping by simultaneously considering structural alignment, local signal profile reorientation, and deformation regularization. Our algorithm also incorporates a multi-kernel strategy to concurrently register anatomical structures at different scales. We demonstrate the efficacy of our approach using in vivo data and report detailed qualitative and quantitative results in comparison with several different registration strategies.

Fast local trust region technique for diffusion tensor registration using exact reorientation and regularization

Diffusion tensor imaging is widely used in brain connectivity research. As more and more studies recruit large numbers of subjects, it is important to design registration methods which are not only theoretically rigorous, but also computationally efficient. However, the requirement of reorienting diffusion tensors complicates and considerably slows down registration procedures, due to the correlated impacts of registration forces at adjacent voxel locations. Based on the diffeomorphic Demons algorithm (Vercauteren , 2009), we propose a fast local trust region algorithm for handling inseparable registration forces for quadratic energy functions. The method guarantees that, at any time and at any voxel location, the velocity is always within its local trust region. This local regularization allows efficient calculation of the transformation update with numeric integration instead of completely solving a large linear system at every iteration. It is able to incorporate exact reorientation and regularization into the velocity optimization, and preserve the linear complexity of the diffeomorphic Demons algorithm. In an experiment with 84 diffusion tensor images involving both pair-wise and group-wise registrations, the proposed algorithm achieves better registration in comparison with other methods solving large linear systems (Yeo , 2009). At the same time, this algorithm reduces the computation time and memory demand tenfold.

Automatic clustering and population analysis of white matter tracts using maximum density paths

We introduce a framework for population analysis of white matter tracts based on diffusion-weighted images of the brain. The framework enables extraction of fibers from high angular resolution diffusion images (HARDI); clustering of the fibers based partly on prior knowledge from an atlas; representation of the fiber bundles compactly using a path following points of highest density (maximum density path; MDP); and registration of these paths together using geodesic curve matching to find local correspondences across a population. We demonstrate our method on 4-Tesla HARDI scans from 565 young adults to compute localized statistics across 50 white matter tracts based on fractional anisotropy (FA). Experimental results show increased sensitivity in the determination of genetic influences on principal fiber tracts compared to the tract-based spatial statistics (TBSS) method. Our results show that the MDP representation reveals important parts of the white matter structure and considerably reduces the dimensionality over comparable fiber matching approaches.

Deformable medical image registration: a survey

Deformable image registration is a fundamental task in medical image processing. Among its most important applications, one may cite: 1) multi-modality fusion, where information acquired by different imaging devices or protocols is fused to facilitate diagnosis and treatment planning; 2) longitudinal studies, where temporal structural or anatomical changes are investigated; and 3) population modeling and statistical atlases used to study normal anatomical variability. In this paper, we attempt to give an overview of deformable registration methods, putting emphasis on the most recent advances in the domain. Additional emphasis has been given to techniques applied to medical images. In order to study image registration methods in depth, their main components are identified and studied independently. The most recent techniques are presented in a systematic fashion. The contribution of this paper is to provide an extensive account of registration techniques in a systematic manner.

Patient-Specific Models of Cardiac Biomechanics

Patient-specific models of cardiac function have the potential to improve diagnosis and management of heart disease by integrating medical images with heterogeneous clinical measurements subject to constraints imposed by physical first principles and prior experimental knowledge. We describe new methods for creating three-dimensional patient-specific models of ventricular biomechanics in the failing heart. Three-dimensional bi-ventricular geometry is segmented from cardiac CT images at end-diastole from patients with heart failure. Human myofiber and sheet architecture is modeled using eigenvectors computed from diffusion tensor MR images from an isolated, fixed human organ-donor heart and transformed to the patient-specific geometric model using large deformation diffeomorphic mapping. Semi-automated methods were developed for optimizing the passive material properties while simultaneously computing the unloaded reference geometry of the ventricles for stress analysis. Material properties of active cardiac muscle contraction were optimized to match ventricular pressures measured by cardiac catheterization, and parameters of a lumped-parameter closed-loop model of the circulation were estimated with a circulatory adaptation algorithm making use of information derived from echocardiography. These components were then integrated to create a multi-scale model of the patient-specific heart. These methods were tested in five heart failure patients from the San Diego Veteran's Affairs Medical Center who gave informed consent. The simulation results showed good agreement with measured echocardiographic and global functional parameters such as ejection fraction and peak cavity pressures.

Improving alignment in Tract-based spatial statistics: evaluation and optimization of image registration

Anatomical alignment in neuroimaging studies is of such importance that considerable effort is put into improving the registration used to establish spatial correspondence. Tract-based spatial statistics (TBSS) is a popular method for comparing diffusion characteristics across subjects. TBSS establishes spatial correspondence using a combination of nonlinear registration and a “skeleton projection” that may break topological consistency of the transformed brain images. We therefore investigated feasibility of replacing the two-stage registration-projection procedure in TBSS with a single, regularized, high-dimensional registration.

To optimize registration parameters and to evaluate registration performance in diffusion MRI, we designed an evaluation framework that uses native space probabilistic tractography for 23 white matter tracts, and quantifies tract similarity across subjects in standard space. We optimized parameters for two registration algorithms on two diffusion datasets of different quality. We investigated reproducibility of the evaluation framework, and of the optimized registration algorithms. Next, we compared registration performance of the regularized registration methods and TBSS. Finally, feasibility and effect of incorporating the improved registration in TBSS were evaluated in an example study.

The evaluation framework was highly reproducible for both algorithms (R² 0.993; 0.931). The optimal registration parameters depended on the quality of the dataset in a graded and predictable manner. At optimal parameters, both algorithms outperformed the registration of TBSS, showing feasibility of adopting such approaches in TBSS. This was further confirmed in the example experiment.

Longitudinal change detection: inference on the diffusion tensor along white matter pathways

Diffusion weighted magnetic resonance imaging (DW-MRI) makes it possible to probe brain connections in vivo. This paper presents a change detection framework that relies on white matter pathways with application to neuromyelitis optica (NMO). The objective is to detect local or global fiber diffusion property modifications between two longitudinal DW-MRI acquisitions of a patient. To this end, we develop two frameworks based on statistical tests on tensor eigenvalues to detect local or global changes along the white matter pathways: a pointwise test that compares tensor populations extracted in bundles cross sections and a fiberwise test that compares paired tensors along all the fiber bundles. Experiments on both synthetic and real data highlight the benefit of considering fiber based statistical tests compared to standard voxelwise strategies.

A robust variational approach for simultaneous smoothing and estimation of DTI

Estimating diffusion tensors is an essential step in many applications - such as diffusion tensor image (DTI) registration, segmentation and fiber tractography. Most of the methods proposed in the literature for this task are not simultaneously statistically robust and feature preserving techniques. In this paper, we propose a novel and robust variational framework for simultaneous smoothing and estimation of diffusion tensors from diffusion MRI. Our variational principle makes use of a recently introduced total Kullback-Leibler (tKL) divergence for DTI regularization. tKL is a statistically robust dissimilarity measure for diffusion tensors, and regularization by using tKL ensures the symmetric positive definiteness of tensors automatically. Further, the regularization is weighted by a non-local factor adapted from the conventional non-local means filters. Finally, for the data fidelity, we use the nonlinear least-squares term derived from the Stejskal-Tanner model. We present experimental results depicting the positive performance of our method in comparison to competing methods on synthetic and real data examples.

Joint T1 and brain fiber log-demons registration using currents to model geometry

We present an extension of the diffeomorphic Geometric Demons algorithm which combines the iconic registration with geometric constraints. Our algorithm works in the log-domain space, so that one can efficiently compute the deformation field of the geometry. We represent the shape of objects of interest in the space of currents which is sensitive to both location and geometric structure of objects. Currents provides a distance between geometric structures that can be defined without specifying explicit point-to-point correspondences. We demonstrate this framework by registering simultaneously T1 images and 65 fiber bundles consistently extracted in 12 subjects and compare it against non-linear T1, tensor, and multi-modal T1 + Fractional Anisotropy (FA) registration algorithms. Results show the superiority of the Log-domain Geometric Demons over their purely iconic counterparts.

Large deformation diffeomorphic registration of diffusion-weighted images with explicit orientation optimization

We seek to compute a diffeomorphic map between a pair of diffusion-weighted images under large deformation. Unlike existing techniques, our method allows any diffusion model to be fitted after registration for subsequent multifaceted analysis. This is achieved by directly aligning the diffusion-weighted images using a large deformation diffeomorphic registration framework formulated from an optimal control perspective. Our algorithm seeks the optimal coordinate mapping by simultaneously considering structural alignment, local fiber reorientation, and deformation regularization. Our algorithm also incorporates a multi-kernel strategy to concurrently register anatomical structures of different scales. We demonstrate the efficacy of our approach using in vivo data and report on detailed qualitative and quantitative results in comparison with several different registration strategies.

Diffeomorphic metric mapping of hybrid diffusion imaging based on BFOR signal basis

In this paper, we propose a large deformation diffeomorphic metric mapping algorithm to align multiple b-value diffusion weighted imaging (mDWI) data, specifically acquired via hybrid diffusion imaging (HYDI), denoted as LDDMM-HYDI. We adopt the work given in Hosseinbor et al. (2012) and represent the q-space diffusion signal with the Bessel Fourier orientation reconstruction (BFOR) signal basis. The BFOR framework provides the representation of mDWI in the q-space and thus reduces memory requirement. In addition, since the BFOR signal basis is orthonormal, the L2 norm that quantifies the differences in q-space signals of any two mDWI datasets can be easily computed as the sum of the squared differences in the BFOR expansion coefficients. In this work, we show that the reorientation of the q-space signal due to spatial transformation can be easily defined on the BFOR signal basis. We incorporate the BFOR signal basis into the LDDMM framework and derive the gradient descent algorithm for LDDMM-HYDI with explicit orientation optimization. Using real HYDI datasets, we show that it is important to consider the variation of mDWI reorientation due to a small change in diffeomorphic transformation in the LDDMM-HYDI optimization.

Role of standardized and study-specific human brain diffusion tensor templates in inter-subject spatial normalization

PURPOSE

To investigate the effect of standardized and study-specific human brain diffusion tensor templates on the accuracy of spatial normalization, without ignoring the important roles of data quality and registration algorithm effectiveness.

MATERIALS AND METHODS

Two groups of diffusion tensor imaging (DTI) datasets, with and without visible artifacts, were normalized to two standardized diffusion tensor templates (IIT2, ICBM81) as well as study-specific templates, using three registration approaches. The accuracy of inter-subject spatial normalization was compared across templates, using the most effective registration technique for each template and group of data.

RESULTS

It was demonstrated that, for DTI data with visible artifacts, the study-specific template resulted in significantly higher spatial normalization accuracy than standardized templates. However, for data without visible artifacts, the study-specific template and the standardized template of higher quality (IIT2) resulted in similar normalization accuracy.

CONCLUSION

For DTI data with visible artifacts, a carefully constructed study-specific template may achieve higher normalization accuracy than that of standardized templates. However, as DTI data quality improves, a high-quality standardized template may be more advantageous than a study-specific template, because in addition to high normalization accuracy, it provides a standard reference across studies, as well as automated localization/segmentation when accompanied by anatomical labels.

Fully automatic stitching of diffusion tensor images in spinal cord

Diffusion tensor imaging (DTI) has become an important tool for studying the spinal cord pathologies. To enable high resolution imaging for modern studies, the DTI technique utilizes a small field of view (FOV) to capture partial human spinal cords. However, normal aging and many other diseases which affect the entire spinal cord increase the desire of acquiring the continuous full-view of the spinal cord. To overcome this problem, this paper presents a novel pipeline for automatic stitching of three-dimensional (3D) DTI of different portions of the spinal cord. The proposed technique consists of two operations, e.g. feature-based registration and adaptive composition to stitch every source image together to create a panoramic image. In the feature-based registration process, feature points are detected from the apparent diffusion coefficient map, and then a novel feature descriptor is designed to characterize feature points directly from a tensor neighborhood. 3D affine transforms are achieved by determining the correspondence matching. In the adaptive composition process, an effective feathering approach is presented to compute the tensors in the overlap region by the Log-Euclidean metrics. We evaluate the algorithm on real datasets from one healthy subject and one adolescent idiopathic scoliosis (AIS) patient. The colored FA maps and fiber tracking results show the effectiveness and accuracy of the proposed stitching framework.

Medical image registration: a review

This paper presents a review of automated image registration methodologies that have been used in the medical field. The aim of this paper is to be an introduction to the field, provide knowledge on the work that has been developed and to be a suitable reference for those who are looking for registration methods for a specific application. The registration methodologies under review are classified into intensity or feature based. The main steps of these methodologies, the common geometric transformations, the similarity measures and accuracy assessment techniques are introduced and described.

DTI analyses and clinical applications in Alzheimer's disease

DTI is one of the most effective MR tools for the investigation of the brain anatomy. In addition to the gray matter, histopathological studies indicate that white matter is also a good target for both the early diagnosis of AD and for monitoring disease progression, which motivates us to use DTI to study AD patients in vivo. There are already a large amount of studies reporting significant differences between AD patients and controls, as well as to predict progression of disease in symptomatic non-demented individuals. Application of these findings in clinical practice remains to be demonstrated.

Fast diffusion tensor registration with exact reorientation and regularization

Diffusion tensor imaging is widely used in brain connectivity study. As more and more group studies recruit a large number of subjects, it is important to design registration methods that are not only theoretically rigorous, but also computationally efficient, for processing large data sets. However, the requirement of reorienting diffusion tensors complicates and slows down the registration, especially for those methods whose scalar-image versions have linear complexity, for example, the Demons algorithm. In this paper, we propose an extension of the Demons algorithm that incorporates exact reorientation and regularization into the calculation of deforming velocity, yet preserving its linear complexity. This method restores the computational efficiency of the Demons algorithm to diffusion images, but does not sacrifice registration goodness. In our experiments, the new algorithm achieved state-of-art performance at a ten-fold decrease of computational time.

Registration and analysis of white matter group differences with a multi-fiber model

Diffusion magnetic resonance imaging has been used extensively to probe the white matter in vivo. Typically, the raw diffusion images are used to reconstruct a diffusion tensor image (DTI). The incapacity of DTI to represent crossing fibers leaded to the development of more sophisticated diffusion models. Among them, multi-fiber models represent each fiber bundle independently, allowing the direct extraction of diffusion features for population analysis. However, no method exists to properly register multi-fiber models, seriously limiting their use in group comparisons. This paper presents a registration and atlas construction method for multi-fiber models. The validity of the registration is demonstrated on a dataset of 45 subjects, including both healthy and unhealthy subjects. Morphometry analysis and tract-based statistics are then carried out, proving that multi-fiber models registration is better at detecting white matter local differences than single tensor registration.

Along-tract statistics allow for enhanced tractography analysis

Diffusion imaging tractography is a valuable tool for neuroscience researchers because it allows the generation of individualized virtual dissections of major white matter tracts in the human brain. It facilitates between-subject statistical analyses tailored to the specific anatomy of each participant. There is prominent variation in diffusion imaging metrics (e.g., fractional anisotropy, FA) within tracts, but most tractography studies use a "tract-averaged" approach to analysis by averaging the scalar values from the many streamline vertices in a tract dissection into a single point-spread estimate for each tract. Here we describe a complete workflow needed to conduct an along-tract analysis of white matter streamline tract groups. This consists of 1) A flexible MATLAB toolkit for generating along-tract data based on B-spline resampling and compilation of scalar data at different collections of vertices along the curving tract spines, and 2) Statistical analysis and rich data visualization by leveraging tools available through the R platform for statistical computing. We demonstrate the effectiveness of such an along-tract approach over the tract-averaged approach in an example analysis of 10 major white matter tracts in a single subject. We also show that these techniques easily extend to between-group analyses typically used in neuroscience applications, by conducting an along-tract analysis of differences in FA between 9 individuals with fetal alcohol spectrum disorders (FASDs) and 11 typically-developing controls. This analysis reveals localized differences between FASD and control groups that were not apparent using a tract-averaged method. Finally, to validate our approach and highlight the strength of this extensible software framework, we implement 2 other methods from the literature and leverage the existing workflow tools to conduct a comparison study.

Diffeomorphic image registration of diffusion MRI using spherical harmonics

Nonrigid registration of diffusion magnetic resonance imaging (MRI) is crucial for group analyses and building white matter and fiber tract atlases. Most current diffusion MRI registration techniques are limited to the alignment of diffusion tensor imaging (DTI) data. We propose a novel diffeomorphic registration method for high angular resolution diffusion images by mapping their orientation distribution functions (ODFs). ODFs can be reconstructed using q-ball imaging (QBI) techniques and represented by spherical harmonics (SHs) to resolve intra-voxel fiber crossings. The registration is based on optimizing a diffeomorphic demons cost function. Unlike scalar images, deforming ODF maps requires ODF reorientation to maintain its consistency with the local fiber orientations. Our method simultaneously reorients the ODFs by computing a Wigner rotation matrix at each voxel, and applies it to the SH coefficients during registration. Rotation of the coefficients avoids the estimation of principal directions, which has no analytical solution and is time consuming. The proposed method was validated on both simulated and real data sets with various metrics, which include the distance between the estimated and simulated transformation fields, the standard deviation of the general fractional anisotropy and the directional consistency of the deformed and reference images. The registration performance using SHs with different maximum orders were compared using these metrics. Results show that the diffeomorphic registration improved the affine alignment, and registration using SHs with higher order SHs further improved the registration accuracy by reducing the shape difference and improving the directional consistency of the registered and reference ODF maps.

ATLAS CONSTRUCTION FROM HIGH ANGULAR RESOLUTION DIFFUSION IMAGING DATA REPRESENTED BY GAUSSIAN MIXTURE FIELDS

Groupwise image registration is an essential part of atlas construction which is a very import and challenging task in medical image analysis. In this paper, we present a novel atlas construction technique using a groupwise registration of high angular resolution diffusion (MR) imaging datasets each of which is represented by a Gaussian Mixture field. To solve the registration problem, an L(2) distance is used to measure the similarity between two Gaussian Mixtures, which leads to an energy function whose gradient can be computed in closed form. A projection method is developed to construct a "sharp" (not blurred) atlas from the result of this groupwise registration. Synthetic and real data experiments are presented to demonstrate the efficacy of the proposed method.

Registration, atlas estimation and variability analysis of white matter fiber bundles modeled as currents

This paper proposes a generic framework for the registration, the template estimation and the variability analysis of white matter fiber bundles extracted from diffusion images. This framework is based on the metric on currents for the comparison of fiber bundles. This metric measures anatomical differences between fiber bundles, seen as global homologous structures across subjects. It avoids the need to establish correspondences between points or between individual fibers of different bundles. It can measure differences both in terms of the geometry of the bundles (like its boundaries) and in terms of the density of fibers within the bundle. It is robust to fiber interruptions and reconnections. In addition, a recently introduced sparse approximation algorithm allows us to give an interpretable representation of the fiber bundles and their variations in the framework of currents. First, we used this metric to drive the registration between two sets of homologous fiber bundles of two different subjects. A dense deformation of the underlying white matter is estimated, which is constrained by the bundles seen as global anatomical landmarks. By contrast, the alignment obtained from image registration is driven only by the local gradient of the image. Second, we propose a generative statistical model for the analysis of a collection of homologous bundles. This model consistently estimates prototype fiber bundles (called template), which capture the anatomical invariants in the population, a set of deformations, which align the geometry of the template to that of each subject and a set of residual perturbations. The statistical analysis of both the deformations and the residuals describe the anatomical variability in terms of geometry (stretching, torque, etc.) and "texture" (fiber density, etc.). Third, this statistical modeling allows us to simulate new synthetic bundles according to the estimated variability. This gives a way to interpret the anatomical features that the model detects consistently across the subjects. This may be used to better understand the bias introduced by the fiber extraction methods and eventually to give anatomical characterization of the normal or pathological variability of fiber bundles.

Group-wise diffeomorphic diffusion tensor image registration

We propose an unbiased group-wise diffeomorphic registration technique to normalize a group of diffusion tensor (DT) images. Our method uses an implicit reference group-wise registration framework to avoid bias caused by reference selection. Log-Euclidean metrics on diffusion tensors are used for the tensor interpolation and computation of the similarity cost functions. The overall energy function is constructed by a diffeomorphic demons approach. The tensor reorientation is performed and implicitly optimized during the registration procedure. The performance of the proposed method is compared with reference-based diffusion tensor imaging (DTI) registration methods. The registered DTI images have smaller shape differences in terms of reduced variance of the fractional anisotropy maps and more consistent tensor orientations. We demonstrate that fiber tract atlas construction can benefit from the group-wise registration by producing fiber bundles with higher overlaps.

Robust clustering of massive tractography datasets

This paper presents a clustering method that detects the fiber bundles embedded in any MR-diffusion based tractography dataset. Our method can be seen as a compressing operation, capturing the most meaningful information enclosed in the fiber dataset. For the sake of efficiency, part of the analysis is based on clustering the white matter (WM) voxels rather than the fibers. The resulting regions of interest are used to define subset of fibers that are subdivided further into consistent bundles using a clustering of the fiber extremities. The dataset is reduced from more than one million fiber tracts to about two thousand fiber bundles. Validations are provided using simulated data and a physical phantom. We see our approach as a crucial preprocessing step before further analysis of huge fiber datasets. An important application will be the inference of detailed models of the subdivisions of white matter pathways and the mapping of the main U-fiber bundles.