Symmetric diffeomorphic registration of fibre orientation distributions

• Published in: NeuroImage (Orlando, Fla.) Vol. 56; no. 3; pp. 1171 - 1180
• Main Authors: Raffelt, David, Tournier, J-Donald, Fripp, Jurgen, Crozier, Stuart, Connelly, Alan, Salvado, Olivier
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01.06.2011; Elsevier Limited
• Subjects: Algorithms; Anisotropy; Attention deficit hyperactivity disorder; Brain - anatomy & histology; Brain research; Data Interpretation, Statistical; Diffusion Magnetic Resonance Imaging - methods; Diffusion MRI; Fibre orientation distribution; HARDI; Humans; Image Processing, Computer-Assisted - methods; Methods; Nerve Fibers - physiology; Nonlinear Dynamics; Pattern Recognition, Automated; Registration; Reproducibility of Results; Software; Spatial normalisation; Studies; Water - chemistry; CC; PVF; PSF; Fibre orientation distribution; Spatial normalisation; SyN; HARDI; DT; MAE; ODF; ANTS; CSD; F3D; Registration; DWI; MD; Diffusion MRI; FOD; FA; T1W; MSD; AFD
• ISSN: 1053-8119; 1095-9572; 1095-9572
• DOI: 10.1016/j.neuroimage.2011.02.014

Registration of diffusion-weighted images is an important step in comparing white matter fibre bundles across subjects, or in the same subject at different time points. Using diffusion-weighted imaging, Spherical Deconvolution enables multiple fibre populations within a voxel to be resolved by computing the fibre orientation distribution (FOD). In this paper, we present a novel method that employs FODs for the registration of diffusion-weighted images. Registration was performed by optimising a symmetric diffeomorphic non-linear transformation model, using image metrics based on the mean squared difference, and cross-correlation of the FOD spherical harmonic coefficients. The proposed method was validated by recovering known displacement fields using FODs represented with maximum harmonic degrees (lmax) of 2, 4 and 6. Results demonstrate a benefit in using FODs at lmax=4 compared to lmax=2. However, a decrease in registration accuracy was observed when lmax=6 was used; this was likely caused by noise in higher harmonic degrees. We compared our proposed method to fractional anisotropy driven registration using an identical code base and parameters. FOD registration was observed to perform significantly better than FA in all experiments. The cross-correlation metric performed significantly better than the mean squared difference. Finally, we demonstrated the utility of this method by computing an unbiased group average FOD template that was used for probabilistic fibre tractography. This work suggests that using crossing fibre information aids in the alignment of white matter and could therefore benefit several methods for investigating population differences in white matter, including voxel based analysis, tensor based morphometry, atlas based segmentation and labelling, and group average fibre tractography. ► A new method for registration of fibre orientation distributions (FOD) is presented. ► Quantitative validation is performed using two FOD image metrics on human brain data. ► The cross-correlation metric outperforms the mean squared difference metric. ► Higher-order information in FODs improves white matter alignment. ► The method is applied to a potential application: group average FOD tractography.