Computational analysis of LDDMM for brain mapping

One goal of computational anatomy (CA) is to develop tools to accurately segment brain structures in healthy and diseased subjects. In this paper, we examine the performance and complexity of such segmentation in the framework of the large deformation diffeomorphic metric mapping (LDDMM) registration method with reference to atlases and parameters. First we report the application of a multi-atlas segmentation approach to define basal ganglia structures in healthy and diseased kids' brains. The segmentation accuracy of the multi-atlas approach is compared with the single atlas LDDMM implementation and two state-of-the-art segmentation algorithms—Freesurfer and FSL—by computing the overlap errors between automatic and manual segmentations of the six basal ganglia nuclei in healthy subjects as well as subjects with diseases including ADHD and Autism. The high accuracy of multi-atlas segmentation is obtained at the cost of increasing the computational complexity because of the calculations necessary between the atlases and a subject. Second, we examine the effect of parameters on total LDDMM computation time and segmentation accuracy for basal ganglia structures. Single atlas LDDMM method is used to automatically segment the structures in a population of 16 subjects using different sets of parameters. The results show that a cascade approach and using fewer time steps can reduce computational complexity as much as five times while maintaining reliable segmentations.

Introduction

Computational anatomy (CA) methods have been at the forefront of neuroimaging studies of neurodevelopment and neurodegeneration. Software such as FSL (Jenkinson et al., 2012), AFNI (Cox, 2012), BrainVoyager (Goebel, 2012), Caret (Van Essen, 2012), FIASCO (Aguirre, 2012), SUMA (Saad and Reynolds, 2012), FreeSurfer (Fischl, 2012), and SPM (Ashburner, 2012) have been used to analyze structural and functional properties of the human brain via magnetic resonance images (MRI), diffusion tensor images (DTI), and functional MRI (fMRI) at a resolution of 1 mm³. In addition, pipelines for registering and visualizing hundreds of images have emerged. These include DTIStudio (Jiang et al., 2006), LONI (Dinov et al., 2006), ANTs (Avants et al., 2008), and 3D Slicer (Pieper et al., 2004). At the heart of these methods is the modeling of anatomical structures and their huge variation, an active field originally formalized as CA (Grenander and Miller, 1998). The main difficulty is the complexity of anatomical structures and the large variation between individuals. Here anatomical structures are represented as a collection of coordinate systems: landmark points, curves, surfaces, and sub-volumes. These structures are represented as deformable templates, with the space of anatomical images being the set generated by the group of diffeomorphic transformations acting on the template with associated probability laws, which describe their variation. The transformations are detailed so that a large family of shapes may be generated with the precise topology of the template maintained. In particular, our methods have demonstrated localized shape differences in multiple sub-cortical structures in neuroimaging studies of Alzheimer's Disease (Qiu et al., 2008a, 2009c), ADHD (Qiu et al., 2009b), Autism (Qiu et al., 2010), schizophrenia (Wang et al., 2008), and Tourette Syndrome (Wang et al., 2007).

The theoretical framework we have adopted is based on the large deformation diffeomorphic metric mapping (LDDMM) algorithm (Beg et al., 2005) and advances have been developed by others as well (Risser et al., 2010, 2011; Auzias et al., 2011). It also has been demonstrated that the best registration methods have high dimensional and diffeomorphic properties that have in common many of the properties incorporated into LDDMM (Klein et al., 2009, 2010). The LDDMM approach allows shape to be uniquely encoded by the vectors normal to the outline of the template (Miller et al., 2006). This property provides the basis of our shape analysis projects and is thus crucial in locating changes in multiple structures caused by disease leading to a better understanding of the effect of disease on neighboring structures.

To date, many MRI studies of subcortical gray matter nuclei have defined a single measure of structural volume. While this has the advantage of being quantitative, it does not give specific information about subregions of atrophying nuclei. This information would be useful in order to determine whether MRI morphometric results correlate with neuropathologic studies, to define better the subregional distribution of atrophy, and for correlation of pathologic changes with clinical features of the disease. Methods of statistical shape analysis have proved useful for studying normal age related changes in subcortical nuclei and for studying a number of other diseases (Bansal et al., 2007; Brun et al., 2009; Chung et al., 2010; Stein et al., 2011; van den Bogaard et al., 2011). Such analysis based on our LDDMM framework (Miller and Qiu, 2009) has proved useful for studying changes in subcortical nuclei in normal aging and several diseases (Csernansky et al., 1998, 2000, 2002; Wang et al., 2007; Qiu and Miller, 2008; Qiu et al., 2008b, 2009a,b,c,d, 2010). For example, analysis with the LDDMM based surface based morphometry (SBM) pipeline revealed heterogeneity of atrophy in subcortical structures in studies of Huntington's disease (Younes et al., 2012) and dementia (Miller et al., 2012).

So as the usage of neuroimaging software expands, it becomes critically important to provide novel analyses that go beyond volumes and/or other scalar quantities at the voxel or regional level. Specifically, the ability to segment subcortical structures in hundreds of scans which can then be used for shape analysis. However, segmentation accuracy and performance is limited by two complementary factors.

The first is the restriction to a single atlas as in MRIStudio which has become widely used in neuroimaging studies of MRI and DTI data as evidenced by meta-data search in Google Scholar. Currently, MRIStudio consists of three programs: DTIStudio (Jiang et al., 2006), DiffeoMap and ROIEditor. In particular, DiffeoMap is a program for image transformation based on LDDMM and ROIEditor uses the results of DiffeoMap to perform image analysis with respect to a single atlas both at the voxel and regional level. The restriction can be compounded by the wide anatomical variability in structures. For example, ventricles differ in shape and geometry in neuroimaging studies of aging and dementia (Reig et al., 2007). Further, proximity of sub-cortical structures can affect the accuracy and performance of automated segmentation algorithms resulting in conflicting results at the volumetric level (Heckemann et al., 2010). Here we explore the incorporation of multi-atlases to drive the segmentation of subcortical structures which should lead to improved statistical significance in shape analysis even from a small sample.

The segmentation problem, which relies on learning anatomical information from pre-segmented training datasets, is usually handled in the Bayesian framework by solving a maximum a posteriori (MAP) estimation problem. The information of multiple pre-segmented training datasets can be combined naturally in the Bayesian framework. Compared with single atlas based segmentation, multi-atlas based segmentation has been shown to be more powerful and more accurate (Rohlfing et al., 2004; Heckemann et al., 2010; Langerak et al., 2010; Lotjonen et al., 2010). More recently we have demonstrated the use of LDDMM with multiple atlases in which both the diffeomorphic coordinate change as well as the atlas being selected is unknown. Such an approach has been shown to be efficient in mediating large deformations in the context of building reliable segmentations of both subcortical and cortical structures (Tang et al., 2012).

This leads to the second problem which is the computational complexity of LDDMM which can be compounded by multiple atlases. By default, LDDMM computes the geodesic along the Riemannian manifold from the template to the target in ten uniform time steps (Beg et al., 2005; Ceritoglu et al., 2009; Tward et al., 2011). As such it takes about 30 min to map one whole brain at 1 mm³ resolution using a parallelized implementation on Intel Xeon CPU E5530 at 2.40GHz with 16 cores. By evolving the flow field over time, it allows for large deformations which come at a cost in terms of computational resources. In the multiple atlas setting, the computational requirements increase as order number of atlases. In an attempt to overcome the computational complexity, LDDMM has inspired alternative diffeomorphic implementations image registration algorithms such as DARTEL (Ashburner, 2007), Diffeo DEMONS (Vercauteren et al., 2009) and Spherical Demons (Yeo et al., 2010) to name but a few. These methods essentially compute the mapping in one time step to generate a “stationary velocity field” in the spirit of exponential maps for finite dimensional groups. While these mappings no longer satisfy the conservation law and therefore metric property giving rise to complete encoding of the diffeomorphic flow of the shape in the initial tangent vector (Miller et al., 2006), their diffeomorphic features make them powerful alternatives for generating segmentations which are accurate. Then metric structure can be derived from the segmentations as demonstrated (Miller and Qiu, 2009). By reformulating the optimization problem, LDDMM can be recast as EPDiff, i.e., an initial value problem (IVP) in which geodesic shooting is used to evolve the initial momenta to match with the target shape (Miller et al., 2006; Marsland and McLachlan, 2007; Ashburner and Friston, 2011; Vialard et al., 2012).

This paper examines performance and complexity of LDDMM with reference to atlases and parameters. A natural way is to examine tradeoff of computational complexity vs. segmentation reliability. Hence we examine segmentations of multiple sub-cortical structures from ongoing large scale neuroimaging studies which are different and more complex than those previously studied (Tang et al., 2012, 2013).

Methods and Data

Data

The first dataset included thirty whole brain, high resolution T1-weighted 3D-volume MPRAGE images (matrix size = 256 × 256, echo time = 3.76 ms, repetition time = 7.99 ms, field of view = 256 mm, slice thickness = 1.0 mm) acquired from a 3T Philips Gyroscan NT scanner (Royal Philips Electronics, Amsterdam, The Netherlands). The dataset included 13 healthy subjects (mean age 10.42 years old; 5 males, and 8 females), 6 male subjects with autism spectrum disorder (ASD) (mean age 9.74 years old) and 11 subjects diagnosed with Attention Deficit/Hyperactivity Disorder (ADHD) (mean age 10.2 years old; 4 males, and 7 females). The second dataset included sixteen whole brain, high resolution T1-weighted 3D-volume MPRAGE images. Eight were healthy subjects (mean age 20.75 years old; all female) acquired from a Philips Intera 3T scanner (Royal Philips Electronics, Amsterdam, The Netherlands; matrix size = 256 × 256, echo time = 3.8 ms, repetition time = 8.1 ms, field of view = 256 mm, slice thickness = 1.0 mm). Four were diagnosed with primary progressive aphasia (mean age 68 years old) acquired from a 1.5T Philips Gyroscan NT (Royal Philips Electronics, Amsterdam, The Netherlands; matrix size = 256 × 256, echo time = 32 ms, repetition time = 6.85 ms, field of view = 230 mm, slice thickness = 1.0 mm). The other 4 were subjects with Alzheimer's disease (AD) (mean 75.6 years old) acquired from a 1.5T Philips Gyroscan NT (Royal Philips Electronics, Amsterdam, The Netherlands; matrix size = 256 × 256, echo time = 3.2 ms, repetition time = 6.9 ms, field of view = 240 mm, slice thickness = 1.2 mm). All data were resampled to isotropic 1 mm³ voxel resolution and image size of 181 × 217 × 181.

LDDMM

Given an atlas image I₀ and a target image I₁, which can be represented as functions I₀,I₁:Ω ⊆ R³ → R on the spatial domain Ω ⊆ R³, the LDDMM algorithm (Beg et al., 2005; Ceritoglu et al., 2009) computes a diffeomorphic transformation φ:Ω → Ω between these images such that I₁ = I₀·φ⁻¹. The diffeomorphism φ = ϕ₁ is defined as the end point of the curve φ = ϕ_t, t ∈ [0,1] satisfying the ordinary differential equation ${\dot{\varphi }}_t=v_t\left({\varphi }_t\right)$, where ϕ₀ = Id is the identity transformation and v_t:Ω → R³, t ∈ [0, 1] is the time dependant velocity vector field of the flow of deformation.

The diffeomorphism φ is calculated as: φ = ϕ₁ = ∫¹₀v_t (ϕ_t)dt with ϕ₀ = Id, where the optimal v_t is estimated by solving the variational problem:

$\begin{array}{cc}\widehat{v}=\underset{v:{\dot{\varphi }}_t=v_t\left({\varphi }_t\right)}{\arg \min }\left({\displaystyle {\int }_0^1{\Vert L{v}_t\Vert }_{L^2}^{2}dt+\frac{1}{{\sigma }^2}{\Vert I_0·{\varphi }_1^{-1}-I_1\Vert }_{L^2}^2}\right)& (1)\end{array}$

Theoretically, the minimizer of equation (1) results in a geodesic path for the curve ϕ_t in the space of diffeomorphisms. To ensure that the solution lies in the space of diffeomorphisms, smoothness is achieved by defining the operator L as: L = −α∇² + γI_{3 × 3}, where ∇² is the Laplacian operator.

In LDDMM, steepest gradient descent approach is used to perform the minimization in equation (1) and the velocity field at each gradient descent iteration k is updated with

$\begin{array}{cc}{v}_t^{k+1}=v_t^k-\epsilon ({\nabla }_{v_t^k}{E}_t),& (2)\end{array}$

where ∇_v E_t is the gradient of the cost in equation (1).

$\begin{array}{cc}{\nabla }_v{E}_t=2v_t-K\text{*}\left(\frac{2}{{\sigma }^2}\left|D{\varphi }_{t,1}^v\right|\nabla J_t^0\left(J_t^0-J_t^1\right)\right),& (3)\end{array}$

where ϕ_{s, t} = ϕ_t · ϕ⁻¹_s, J⁰_t = I₀ · ϕ_{t, 0}, J¹_t = I₁ · ϕ_{t, 1}, |Dϕ^v_{t, 1}| is the determinant of the Jacobian matrix, K = (L^†L)⁻¹ and * is the convolution operation.

In the numerical implementation of LDDMM, the time parameter t of the flow is discretized with a fixed total number of timesteps T, where T = 10 is selected as the default value in general. Selection of a smaller T causes gradient descent to terminate with a higher final mismatch error ‖I₀ · φ⁻¹ − I₁‖²_L² between the registered atlas image and the target image. In this article T = 1 is chosen to approximate the small deformation setup for comparison with T = 10.

The background space of the images is represented with Ω = [0, 1]³ and the convolution operation in equation (3) is calculated in Fourier domain. The operator K acts as a low pass filter at each iteration of gradient descent and the parameters α and γ controls the amount of smoothing and the elasticity of the deformation. Selection of these parameters depends on the size of the deformation (in pixels) necessary to register the features of the atlas image to the features of the target image. Resampling the atlas and target images to same resolution and doing an initial rigid or affine alignment as a preprocessing step for LDDMM allows the selection of α/γ in the range of 0.01–0.001 to accurately register the features of MR brain data with a typical resolution in the range of 0.5–2 mm/voxel. As α/γ decreases in the 0.01–0.001 range, the image matching quality increases and the deformations become more local and elastic (Beg et al., 2005; Ceritoglu et al., 2009; Tward et al., 2011).

In the cascading implementation of LDDMM (Ceritoglu et al., 2009), The final diffeomorphic deformation φ between images I₀ andI₁ is calculated by combining deformations φ₁, φ₂, … , φ_n and φ = φ₁⁻¹ · φ⁻¹₂ · … · φ⁻¹_n where (i) φ₁ is calculated between I₀ and I₁, (ii) φ₂ is calculated between I₀ · φ⁻¹₁ and I₁, (iii) φ₃ is calculated between I₀ · φ⁻¹₁ · φ⁻¹₂ and I₁. In each step, the regularization parameter γ is fixed to be 1 and α is decreased to ensure that final registration is more robust and does not converge to an apparently incorrect solution compared to computing φ directly between I₀ and I₁ using very small α value (Ceritoglu et al., 2009) directly without cascading. The cascading implementation does not calculate a geodesic trajectory between the images and the mathematical properties of LDDMM such as metric distance and the encoding of the diffeomorphic flow of the shape with the initial tangent vector are not retained.

It is worth noting that geodesic trajectories can be obtained via multi-scale approaches (Risser et al., 2011; Vialard et al., 2012) where a weighted sum of Gaussian kernels is used to define the operator K in equation (3) control the elasticity of the deformation at each scale.

Segmentation Via Single-Atlas

The single-subject atlas is the JHU-DTI-MNI atlas (Oishi et al., 2009), which is a single-subject template in the ICBM-DTI-81 space. The template has an isotropic 1 mm³ voxel resolution and image size of 181 × 217 × 181. Histogram matching was applied to the atlas and subject image, followed by initial affine alignment of the atlas onto the subject via AIR (Woods et al., 1998a,b). MRIStudio was then used to register the subject to the atlas. The resulting diffeomorphism was used to transfer the ROIs defined in the atlas to the subject space resulting in ROIs segmented in the subject space. This procedure is illustrated in Figure 1. All experiments were performed using Intel Xeon CPU E5530 at 2.40GHz with 16 cores.

FIGURE 1

Figure 1. Example axial slices of atlas and subject T1 weighted MR image and the boundaries of their ROI labels at each registration step. (A) Original atlas image and atlas ROIs; (B) atlas image and ROIs after linear AIR transformation; (C) atlas image and ROIs after AIR and LDDMM transformations; (D) subject image and subject ROIs.

Segmentation Via Multi-Atlas Fusion

The algorithm for multi-atlas segmentation using LDDMM (Tang et al., 2012, 2013) is briefly summarized as follows. Assume there are multiple MRI atlases, each of which contains a collection of locally-defined charts obtained from manually segmented structures. A MAP approach is used to estimate the high dimensional segmentations from the class of generative models representing the observed MRI which is assumed to be a Gaussian random fields conditioned on the atlas charts and diffeomorphic change of coordinates of each chart. The charts and their diffeomorphic correspondences are unknown and viewed as latent or hidden variables. The expectation-maximization (EM) algorithm yields the likelihood-fusion equation which is maximized by the a posteriori estimator of the segmentation labels. The fused likelihoods are modeled as conditional Gaussian random fields with mean fields a function of each atlas chart under its diffeomorphic change of coordinates onto the target. The conditional-mean in the EM algorithm specifies the convex weights with which the chart-specific likelihoods are fused. The multiple atlases with the associated convex weights imply that the posterior distribution is a multi-modal representation of the measured MRI. As with the single-atlas segmentation, the resulting diffeomorphism was used to transfer the ROIs to the subject space resulting in ROIs segmented in the subject space.

Comparison Metrics

Automated segmentation of subcortical structures was compared with manual segmentations of the same subjects. The manual segmentations are based on anatomical definitions (Wang et al., 2007; Qiu et al., 2009b, 2010) with the assistance of an atlas (Mai et al., 1997); for details see http://caportal.cis.jhu.edu/protocols. Manual segmentation was performed with the open source software Seg3D which follows the radiological convention for displaying images.

Accuracy for each ROI was quantified by three metrics. The first is the kappa statistic (Cohen, 1960) is defined by κ = (p_agree − p_random)/(1 − p_random)where p_agree is the fraction of voxels in which the given segmentation agrees with the manual segmentation, and p_random is the fraction one would expect by random chance (based only on the volumes of foreground and background). Note that κ is biased by the volume size of the structure; generally, the bigger the structure, the higher the kappa statistic is. For brain structures, a value of κ = 0.8is considered quite good. The second is the volume error which quantifies volume difference between two label defined by VD(L_A, L_M) = 100|V(L_A) − V(L_M)|/V(L_M) where V(L_A)and V(L_M) are respectively the volume size of the automated and manual segmentation. The third is the L₁ misclassification error (Miller et al., 2000) defined by L₁ = |A ∪ M| − |A ∩ M|/|A ∪ M|/2 where A and Mrespectively are the set of automated and manual labeled voxels.

Statistics

In the first experiment, One-Way ANOVA was performed on κ and L₁ from the automated and manual segmentations. Multiple comparisons were conducted using Tukey's HSD (honestly significant difference) test between pair wise groups. In the second experiment, two statistical analyses were performed to examine the effects of LDDMM parameters—T and α —on the segmentation accuracy of single-atlas based registration. In the first analysis, meanL₁ and volume errors were calculated for each subjects by averaging the errors of the segmented ROIs for that subject. Then Two-Way ANOVA was conducted to examine the effects of changing T and α on mean L₁ and volume error. Multiple comparison tests using Tukey's HSD was carried out after ANOVA. In the second analysis, the effect of T on the errors for each ROI was analyzed separately via one sided paired samples t-test between T = 1 and T = 10 results for given α.

Results

For the first experiment, multi-atlas segmentation with cascading values of α = 0.01, 0.005, 0.002 and T = 10 was compared with single-atlas segmentation with the same values of α and T, FSL (version 5.0) and FreeSurfer (version 5.2.0). Six structures were analyzed (left and right pairs of caudate, putamen and globus pallidus). Mean and standard deviation of kappa overlap and L₁ errors are shown in Figure 2. One-Way ANOVA results (F_L.caudate = 155.57, F_R.caudate = 196.8, F_{L.globuspallidus} = 95.45, F_R.pallidus = 96.85, F_L.putamen = 182.73, F_R.putamen = 228.23, df = 3, P < 0.001) for kappa overlap and (F_L.caudate = 163.08, F_R.caudate = 187.79, F_{L.globuspallidus} = 73.89, F_{R.globuspallidus} = 75.90, F_L.putamen = 186.88, F_R.putamen = 244.95, df = 3, P < 0.001) and the multiple comparisons test with Tukey's HSD showed that for each of six structures considered, the automated segmentations obtained from multi-atlas LDDMM are statistically more accurate than those from the other three methods. For the second experiment, the following combinations were used: α = {0.01, 0.005, 0.002, 0.002with cascading} and T = {1, 10}. So, a total of (4 α-values) × (2 T-values) = 8 different automated segmentations were calculated for each subject using LDDMM single-atlas segmentation. Fourteen structures (left and right pairs of caudate, putamen, globus pallidus, hippocampus, amygdala, thalamus, ventricle) were analyzed. Figure 3 (top and middle rows) shows the mean and standard deviation of these errors computed over 16 subjects. There was no significant interaction between the effects of T and α, (F = 1.6785, df = 3, P = 0.17527). Simple main effects analysis showed that the effects of both α (F = 7.5801, df = 3, P = 0.00011) and T (F = 9.3323, df = 1, P = 0.0028) are statistically significant. A similar analysis on mean volume error showed that there was no significant interaction between the effects of T and α, (F = 0.16, df = 3, P = 0.9230). Simple main effects analysis of the mean volume error showed that the effect of α (F = 4.1938, df = 3, P = 0.0073) was statistically significant but the effect of T (F = 2.7054, df = 1, P = 0.1026) was not statistically significant. Multiple comparisons test with Tukey's HSD, showed that there is a significant difference (P < 0.01) between using a high α (0.01) and low α (0.002 with and without cascading). Tables 1, 2 shows the statistically significant differences (P < 0.01) for L₁ error and volume error respectively for the results at T = 1 and T = 10 at a significance alpha level of 0.01.

FIGURE 2

Figure 2. Mean and standard deviation of mean Kappa statistics (top) and mean L₁ errors (bottom) for left and right pairs of caudate, putamen and pallidus in 30 subjects via four different segmentation methods: multi-atlas (blue), single-atlas (red), FSL (green), and FreeSurfer (purple).

FIGURE 3

Figure 3. Mean and standard deviation of mean L₁ error (top), mean volume error (middle) and LDDMM running times (bottom) calculated from a population of 16 subjects for different values of α, T = 1 (blue) and T = 10 (red). Mean errors were calculated using 7 pairs of caudate, putamen, pallidus, amygdala, hippocampus, thalamus, and ventricle for each subject.

TABLE 1

Table 1. For each ROI and for each α in LDDMM, the L₁ error results for T = 1 and T = 10 were compared using one sided paired sample t-test.

TABLE 2

Table 2. For each ROI and for each α in LDDMM, the volume error results for T = 1 and T = 10 were compared using one sided paired samples t-test.

Figure 3 (bottom row) shows the average running times of LDDMM calculations for any selection of α and T for the populaton of 16 subjects. Computational complexity of each gradient descent iteration in LDDMM depends on T linearly. Smaller α also increases the computation time.

Discussion

In this paper, we examined the computational complexity in mapping multiple atlases for segmentation of subcortical structures—caudate, putamen, and globus pallidus. Our results show that a cascade of α parameters and selecting fewer time steps can reduce computational complexity as much as five times while obtaining reliable segmentations.

We also compared the segmentations with two other methods—FSL and FreeSurfer. LDDMM based segmentation is capable of achieving comparable or superior accuracy as measured by either the kappa statistic or the L₁ error. Extending the single-atlas LDDMM based segmentation to multiple atlases, we found that the segmentation accuracy was significantly increased. The even split between male and female subjects in the dataset of 30 subjects suggests that our algorithm is robust with respect to gender. Then in terms of pathological variability, the datasets consisted of control subjects, subjects with autism, as well as subjects with ADHD. So our multi-atlas LDDMM segmentation algorithm can be applied to a wide variety of subjects. This is crucial for clinical neuroimaging studies since the ultimate goal of any segmentation algorithm is to be able to accurately delineate brain structures with different diseases. The performance of multi-atlas LDDMM, in terms of segmenting subcortical structures of a different dataset, has been compared with STAPLE which is a classic label-fusion based segmentation method (Warfield et al., 2004). The segmentations were demonstrated to be superior to those from STAPLE (Tang et al., 2012). The high segmentation accuracy obtained from multi-atlas LDDMM is at the cost of more running time, given that we need to first segment multiple atlases and we need to do multiple registrations. This issue led to the second experiment—analyzing the two parameters α and T that determines the speed, memory requirement and accuracy of LDDMM, trying to find the optimal α and T that will speed up the registration without affecting segmentation accuracy.

In LDDMM, α and T are the two most important parameters that affect the smoothness of the computed diffeomorphism and the registration accuracy between the images that are registered. The results of the second experiment showed that decreasing α, decreased the mean and standard deviation of the errors and the results are statistically significant. Increasing T and using a large deformation setup instead of a small deformation setup also decreased the errors (Table 1) but results are less significant. Therefore, in general a small α should be selected for better registration. When each ROI is considered separately (first 3 columns of Table 1 and Table 2), using a large deformation map (T = 10) instead of a small deformation map (T = 1) decreased the errors further with smaller α for given ROIs. Instead of using a single α to compute the final transformation with LDDMM, using a cascading scheme and decreasing α gradually to compute the final transformation allows the selection of a smaller T (last column of Table 2 and Table 2). This result is very important especially for multi-atlas LDDMM where a large number of LDDMM computations are necessary for superior segmentation accuracy. Thus in the extreme case of selecting T = 1 leads to the mappings that capture only small deformations while increasing T with finer discretization of time flow leads to the mappings that capture larger deformations and are more smooth.

Considered by some to be the “state of art” registration method (Ashburner and Friston, 2011), LDDMM has been shown to yield accurate segmentation of subcortical structures in diverse datasets. Accuracy is obtained at a cost of computational complexity. As in the original pipeline (Miller and Qiu, 2009) in which segmentations were used to drive geodesic representations from which momentum and statistics were performed, the cascade approach described here allows us to generate intermediate analyses such as segmentations in a computationally efficient way. This is extremely important in our multi-atlas formulation with complexity O(N), where N is the number of atlases. Thus via MRIStudio, LDDMM offers the potential for more precise and sensitive analysis of anatomical structures in neuroimaging studies.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The research here was supported by NIH grants R01 EB000975, P41 EB015909, R01 EB008171, R01 MH084803, S10 RR025053, R01 MH085328, R01 NS048527, P30 HD024061, R01 MH078160, UL1 TR 000424, P50-AG005146 and P50-AG 021334; The images acquired from Alzheimer's patients and age-matched controls were obtained with the help of a methods development grant from Glaxo-Smith-Kline, awarded to Drs. Marilyn Albert and Constantine G. Lyketsos.

References

Aguirre, G. K. (2012). FIASCO, VoxBo, and MEDx: behind the code. Neuroimage 62, 765–767. doi: 10.1016/j.neuroimage.2012.02.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113. doi: 10.1016/j.neuroimage.2007.07.007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ashburner, J. (2012). SPM: a history. Neuroimage 62, 791–800. doi: 10.1016/j.neuroimage.2011.10.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ashburner, J., and Friston, K. J. (2011). Diffeomorphic registration using geodesic shooting and Gauss-Newton optimisation. Neuroimage 55, 954–967. doi: 10.1016/j.neuroimage.2010.12.049

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Auzias, G., Colliot, O., Glaunes, J. A., Perrot, M., Mangin, J. F., Trouve, A., et al. (2011). Diffeomorphic brain registration under exhaustive sulcal constraints. IEEE Trans. Med. Imaging 30, 1214–1227. doi: 10.1109/TMI.2011.2108665

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Avants, B. B., Epstein, C. L., Grossman, M., and Gee, J. C. (2008). Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med. Image Anal. 12, 26–41. doi: 10.1016/j.media.2007.06.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bansal, R., Staib, L. H., Xu, D., Zhu, H., and Peterson, B. S. (2007). Statistical analyses of brain surfaces using Gaussian random fields on 2-D manifolds. IEEE Trans. Med. Imaging 26, 46–57. doi: 10.1109/TMI.2006.884187

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beg, M. F., Miller, M. I., Trouve, A., and Younes, L. (2005). Computing large deformation metric mapping via geodesic flows of diffeomorphisms. Int. J. Comput. Vis. 61, 139–157. doi: 10.1023/B:VISI.0000043755.93987.aa

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Brun, C. C., Nicolson, R., Lepore, N., Chou, Y. Y., Vidal, C. N., Devito, T. J., et al. (2009). Mapping brain abnormalities in boys with autism. Hum. Brain Mapp. 30, 3887–3900. doi: 10.1002/hbm.20814

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ceritoglu, C., Oishi, K., Li, X., Chou, M. C., Younes, L., Albert, M., et al. (2009). Multi-contrast large deformation diffeomorphic metric mapping for diffusion tensor imaging. Neuroimage 47, 618–627. doi: 10.1016/j.neuroimage.2009.04.057

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chung, M. K., Worsley, K. J., Nacewicz, B. M., Dalton, K. M., and Davidson, R. J. (2010). General multivariate linear modeling of surface shapes using SurfStat. Neuroimage 53, 491–505. doi: 10.1016/j.neuroimage.2010.06.032

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cohen, J. (1960). A coefficient of agreement for nominal scales. Educ. Psychol. Meas. 20, 37–46. doi: 10.1177/001316446002000104

CrossRef Full Text

Cox, R. W. (2012). AFNI: what a long strange trip it's been. Neuroimage 62, 743–747. doi: 10.1016/j.neuroimage.2011.08.056

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Csernansky, J. G., Joshi, S., Wang, L., Haller, J. W., Gado, M., Miller, J. P., et al. (1998). Hippocampal morphometry in schizophrenia by high dimensional brain mapping. Proc. Natl. Acad. Sci. U.S.A. 95, 11406–11411. doi: 10.1073/pnas.95.19.11406

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Csernansky, J. G., Wang, L., Jones, D., Rastogi-Cruz, D., Posener, J. A., Heydebrand, G., et al. (2002). Hippocampal deformities in schizophrenia characterized by high dimensional brain mapping. Am. J. Psychiatry 159, 2000–2006. doi: 10.1176/appi.ajp.159.12.2000

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Csernansky, J. G., Wang, L., Joshi, S., Miller, J. P., Gado, M., Kido, D., et al. (2000). Early DAT is distinguished from aging by high-dimensional mapping of the hippocampus. Neurology 55, 1636–1643. doi: 10.1212/WNL.55.11.1636

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Dinov, I. D., Valentino, D., Shin, B. C., Konstantinidis, F., Hu, G., Mackenzie-Graham, A., et al. (2006). LONI visualization environment. J. Digit. Imaging 19, 148–158. doi: 10.1007/s10278-006-0266-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fischl, B. (2012). FreeSurfer. Neuroimage 62, 774–781. doi: 10.1016/j.neuroimage.2012.01.021

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Goebel, R. (2012). BrainVoyager–past, present, future. Neuroimage 62, 748–756. doi: 10.1016/j.neuroimage.2012.01.083

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Grenander, U., and Miller, M. I. (1998). Computational anatomy: an emerging discipline. Quart. Appl. Math. 56, 617–694.

Heckemann, R. A., Keihaninejad, S., Aljabar, P., Rueckert, D., Hajnal, J. V., Hammers, A., et al. (2010). Improving intersubject image registration using tissue-class information benefits robustness and accuracy of multi-atlas based anatomical segmentation. Neuroimage 51, 221–227. doi: 10.1016/j.neuroimage.2010.01.072

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Jenkinson, M., Beckmann, C. F., Behrens, T. E., Woolrich, M. W., and Smith, S. M. (2012). Fsl. Neuroimage 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Jiang, H., Van Zijl, P. C., Kim, J., Pearlson, G. D., and Mori, S. (2006). DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking. Comput. Methods Programs Biomed. 81, 106–116. doi: 10.1016/j.cmpb.2005.08.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Klein, A., Andersson, J., Ardekani, B. A., Ashburner, J., Avants, B., Chiang, M. C., et al. (2009). Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration. Neuroimage 46, 786–802. doi: 10.1016/j.neuroimage.2008.12.037

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Klein, A., Ghosh, S. S., Avants, B., Yeo, B. T., Fischl, B., Ardekani, B., et al. (2010). Evaluation of volume-based and surface-based brain image registration methods. Neuroimage 51, 214–220. doi: 10.1016/j.neuroimage.2010.01.091

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Langerak, T. R., van der Heide, U. A., Kotte, A. N. T. J., Viergever, M. A., Van Vulpen, M., and Pluim, J. P. W. (2010). Label fusion in atlas-based segmentation using a selective and iterative method for performance level estimation (SIMPLE). IEEE Trans. Med. Imaging 29, 2000–2008. doi: 10.1109/TMI.2010.2057442

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lotjonen, J. M. P., Wolz, R., Koikkalainen, J. R., Thurfjell, L., Waldemar, G., Soininen, H., et al. (2010). Fast and robust multi-atlas segmentation of brain magnetic resonance images. Neuroimage 49, 2352–2365. doi: 10.1016/j.neuroimage.2009.10.026

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mai, J. K., Assheuer, J., and Paxinos, G. (1997). Atlas of the Human Brain. San Diego, CA: Academic Press.

Marsland, S., and McLachlan, R. (2007). “A hamiltonian particle method for diffeomorphic image registration,” in Proceedings of the 20th International Conference on Information Processing in Medical Imaging (Kerkrade). eds N. Karssemeijer, and B. Lelieveldt (Berlin; New York, NY: Springer-Verlag), 396–407.

Miller, M. I., Massie, A. B., Ratnanather, J. T., Botteron, K. N., and Csernansky, J. G. (2000). Bayesian construction of geometrically based cortical thickness metrics. Neuroimage 12, 676–687. doi: 10.1006/nimg.2000.0666

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miller, M. I., and Qiu, A. (2009). The emerging discipline of computational functional anatomy. Neuroimage 45, S16–S39. doi: 10.1016/j.neuroimage.2008.10.044

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miller, M. I., Trouve, A., and Younes, L. (2006). Geodesic Shooting for Computational Anatomy. J. Math. Imaging Vis. 24, 209–228. doi: 10.1007/s10851-005-3624-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miller, M. I., Younes, L., Ratnanather, J. T., Brown, T., Reigel, T., Trinh, H., et al. (2012). “Amygdala atrophy in MCI/Alzheimer's disease in the BIOCARD cohort based on diffeomorphic morphometry,” in MICCAI 2012 Workshop on Novel Biomarkers for Alzheimer's Disease and Related Disorders, eds L. Wang, P. A. Yushkevich, and S. Ourselin. Available on http://picsl.upenn.edu/wiki/pmwiki.php?n = Miccai2012/program, 155–166.

Oishi, K., Faria, A., Jiang, H., Li, X., Akhter, K., Zhang, J., et al. (2009). Atlas-based whole brain white matter analysis using large deformation diffeomorphic metric mapping: application to normal elderly and Alzheimer's disease participants. Neuroimage 46, 486–499. doi: 10.1016/j.neuroimage.2009.01.002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pieper, S., Halle, M., and Kikinis, R. (2004). “3D Slicer,” in IEEE International Symposium on Biomedical Imaging: Nano to Macro. Arlington, VA: IEEE Society.

Qiu, A., Adler, M., Crocetti, D., Miller, M., and Mostofsky, S. (2010). Basal ganglia shapes predict social, communication, and motor dysfunctions in boys with autism spectrum disorder. J. Am. Acad. Child Adolesc. Psychiatry 49, 539–551.

Pubmed Abstract | Pubmed Full Text

Qiu, A., Albert, M., Younes, L., and Miller, M. I. (2009a). Time sequence diffeomorphic metric mapping and parallel transport track time-dependent shape changes. Neuroimage 45, S51–S60. doi: 10.1016/j.neuroimage.2008.10.039

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., Crocetti, D., Adler, M., Mahone, E., Mb, D., Miller, M., et al. (2009b). Basal ganglia volume and shape in children with attention deficit hyperactivity disorder. Am. J. Psychiatry 166, 74–82. doi: 10.1176/appi.ajp.2008.08030426

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., Fennema-Notestine, C., Dale, A., and Mi, M. (2009c). Regional shape abnormalities in mild cognitive impairment and Alzheimer's disease. Neuroimage 45, 656–661. doi: 10.1016/j.neuroimage.2009.01.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., Wang, L., Younes, L., Harms, M., Ratnanather, J., Miller, M., et al. (2009d). Neuroanatomical asymmetry patterns in individuals with schizophrenia and their non-psychotic siblings. Neuroimage 47, 1221–1229. doi: 10.1016/j.neuroimage.2009.05.054

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., and Miller, M. I. (2008). Multi-structure network shape analysis via normal surface momentum maps. Neuroimage 42, 1430–1438. doi: 10.1016/j.neuroimage.2008.04.257

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., Younes, L., Miller, M., and Csernansky, J. (2008a). Parallel transport in diffeomorphisms distinguishes the time-dependent pattern of hippocampal surface deformation due to healthy aging and the dementia of the Alzheimer's type. Neuroimage 40, 68–76. doi: 10.1016/j.neuroimage.2007.11.041

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Qiu, A., Younes, L., and Miller, M. I. (2008b). Intrinsic and extrinsic analysis in computational anatomy. Neuroimage 39, 1803–1814. doi: 10.1016/j.neuroimage.2007.08.043

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Reig, S., Penedo, M., Gispert, J. D., Pascau, J., Sanchez-Gonzalez, J., Garcia-Barreno, P., et al. (2007). Impact of ventricular enlargement on the measurement of metabolic activity in spatially normalized PET. Neuroimage 35, 748–758. doi: 10.1016/j.neuroimage.2006.12.015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Risser, L., Vialard, F. X., Wolz, R., Holm, D. D., and Rueckert, D. (2010). Simultaneous fine and coarse diffeomorphic registration: application to atrophy measurement in Alzheimer's disease. Med. Image Comput. Comput. Assist. Interv. 13, 610–617.

Pubmed Abstract | Pubmed Full Text

Risser, L., Vialard, F. X., Wolz, R., Murgasova, M., Holm, D. D., Rueckert, D., et al. (2011). Simultaneous multi-scale registration using large deformation diffeomorphic metric mapping. IEEE Trans. Med. Imaging 30, 1746–1759. doi: 10.1109/TMI.2011.2146787

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rohlfing, T., Brandt, R., Menzel, R., and Maurer, C. R. (2004). Evaluation of atlas selection strategies for atlas-based image segmentation with application to confocol microscopy images of bee brains. Neuroimage 21, 1428–1442. doi: 10.1016/j.neuroimage.2003.11.010

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Saad, Z. S., and Reynolds, R. C. (2012). Suma. Neuroimage 62, 768–773. doi: 10.1016/j.neuroimage.2011.09.016

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Stein, J. L., Hibar, D. P., Madsen, S. K., Khamis, M., McMahon, K. L., De Zubicaray, G. I., et al. (2011). Discovery and replication of dopamine-related gene effects on caudate volume in young and elderly populations (N = 1198) using genome-wide search. Mol. Psychiatry 16, 927–37, 881.

Pubmed Abstract | Pubmed Full Text

Tang, X., Mori, S., and Miller, M. I. (2012). “Segmentation via multi-atlas LDDMM,” in MICCAI 2012 Workshop on Multi-Atlas Labeling, 2012 Nice, eds B. A. Landman and S. K. Warfield (Nice: CreateSpace Independent Publishing Platform), 123–133.

Tang, X., Oishi, K., Faria, A. V., Hillis, A. E., Albert, M. S., Mori, S., et al. (2013). Bayesian parameter estimation and segmentation in the multi-atlas random orbit model. PLoS ONE 8:e65591. doi: 10.1371/journal.pone.0065591

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tward, D. J., Ceritoglu, C., Kolasny, A., Sturgeon, G. M., Segars, W. P., Miller, M. I., et al. (2011). Patient specific dosimetry phantoms using multichannel LDDMM of the whole body. Int. J. Biomed. Imaging 2011:481064. doi: 10.1155/2011/481064

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

van den Bogaard, S. J. A., Dumas, E. M., Acharya, T. P., Johnson, H., Langbehn, D. R., Scahill, R. I., et al. (2011). Early atrophy of pallidum and accumbens nucleus in Huntington's disease. J. Neurol. 258, 412–420. doi: 10.1007/s00415-010-5768-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Van Essen, D. C. (2012). Cortical cartography and Caret software. Neuroimage 62, 757–764. doi: 10.1016/j.neuroimage.2011.10.077

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Vercauteren, T., Pennec, X., Perchant, A., and Ayache, N. (2009). Diffeomorphic demons: efficient non-parametric image registration. Neuroimage 45, S61–S72. doi: 10.1016/j.neuroimage.2008.10.040

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Vialard, F.-X., Risser, L., Rueckert, D., and Cotter, C. J. (2012). Diffeomorphic 3D image registration via geodesic shooting using an efficient adjoint calculation. Int. J. Comput. Vis. 97, 229–241. doi: 10.1007/s11263-011-0481-8

CrossRef Full Text

Wang, L., Lee, D. Y., Bailey, E., Hartlein, J. M., Gado, M. H., Miller, M. I., et al. (2007). Validity of large-deformation high dimensional brain mapping of the basal ganglia in adults with Tourette syndrome. Psychiatry Res. 154, 181–190. doi: 10.1016/j.pscychresns.2006.08.006

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wang, L., Mamah, D., Harms, M. P., Karnik, M., Price, J. L., Gado, M. H., et al. (2008). Progressive deformation of deep brain nuclei and hippocampal-amygdala formation in schizophrenia. Biol. Psychiatry 64, 1060–1068. doi: 10.1016/j.biopsych.2008.08.007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Warfield, S. K., Zou, K. H., and Wells, W. M. (2004). Simultaneous truth and performance level estimation (STAPLE): an algorithm for the validation of image segmentation. IEEE Trans. Med. Imaging 23, 903–921.

Pubmed Abstract | Pubmed Full Text

Woods, R. P., Grafton, S. T., Holmes, C. J., Cherry, S. R., and Mazziotta, J. C. (1998a). Automated image registration: I. General methods and intrasubject, intramodality validation. J. Comput. Assist. Tomogr. 22, 139–152. doi: 10.1097/00004728-199801000-00027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Woods, R. P., Grafton, S. T., Watson, J. D., Sicotte, N. L., and Mazziotta, J. C. (1998b). Automated image registration: II. Intersubject validation of linear and nonlinear models. J. Comput. Assist. Tomogr. 22, 153–165. doi: 10.1097/00004728-199801000-00028

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yeo, B. T., Sabuncu, M. R., Vercauteren, T., Ayache, N., Fischl, B., and Golland, P. (2010). Spherical demons: fast diffeomorphic landmark-free surface registration. IEEE Trans. Med. Imaging 29, 650–668. doi: 10.1109/TMI.2009.2030797

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Younes, L., Ratnanather, J. T., Brown, T., Aylward, E., Nopoulos, P., Pierson, R., et al. (2012). Regionally selective atrophy of subcortical structures in prodromal HD as revealed by statistial shape analysis. Hum. Brain Mapp. doi: 10.1002/hbm.22214. [Epub ahead of print].

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text