Comparing the performances of three diffusion kurtosis tensor estimation algorithms via a ground truth diffusion template derived from HCP data
Daniel V. Olson1, Volkan Arpinar2, and L. Tugan Muftuler2

1Biophysics, Medical College of Wisconsin, Milwaukee, WI, United States, 2Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, United States

Synopsis

A diffusion weighted imaging template is proposed, which is derived from the in vivo data from the HCP database. Diffusion weighted signals are generated from this template using the DKI tensor model, and diffusion tensor and kurtosis tensor metric maps are produced. These maps established the ground truth, against which the outputs of different DKI tensor estimation algorithms were compared. Rician noise is added to simulate typical diffusion MRI acquisitions with different SNR levels. The performances of the algorithms are then compared via voxel-wise Mean Square Error and bias-plus-variance decomposition to determine the optimal algorithm for the desired application.

Purpose

Diffusion kurtosis imaging (DKI) is becoming increasingly popular in diffusion weighted imaging due to its higher sensitivity to tissue microstructure compared to conventional DTI1, while remaining within a clinically feasible scan time. However, the kurtosis tensor (KT) model is also more sensitive to noise and other data artifacts. To improve robustness and accuracy, several techniques have been proposed that address various components of the fitting process. However, there is no study to this date that compared the accuracy and robustness of these techniques against some gold standard. This work compares three KT algorithms: sparsity constraint2 (6KT), iterative outlier removal3 (REKINDLE), and directional weighting and regularization4 (DWAR). The data set to be used as the ground truth is derived from a subset of the HCP data using a novel approach.

For least-squares algorithms, the mean square error (MSE) can be decomposed into the squared bias (bias2), variance, and an irreducible error term. This is known as bias-plus-variance decomposition5. Bias2 error results from simplifying assumptions built into the method. Variance error is caused by the model fitting noise. By minimizing both of these terms, there is a balance between accuracy and robustness to noise. Here, our goal is to evaluate the effectiveness of the kurtosis tensor fitting algorithms with an in-vivo-based “ground truth” (GT) dataset and bias-plus-variance decomposition.

Methods

Preprocessed diffusion-weighted and T1-weighted images from 10 subjects (6 female, 22-30yrs) were retrieved from the WU-Minn HCP database6. Diffusion weighted images (DWI) were processed through software developed in-house to estimate DKI tensors based on the algorithm by Jensen et al7 to calculate 22 tensor elements. Then, the tensor images were registered to standard space in FSL8 and averaged across subjects to create GT tensor images. From these images, first the GT metric images (mean diffusivity (MD), fractional anisotropy (FA), and mean kurtosis (MK)) were generated (Figure 1). Then, a set of DWI images was generated using the same q-space sampling scheme as the HCP. A white matter (WM) mask was applied to reduce processing time.

Using the synthetically generated DWI dataset as a baseline, Rician noise was added to achieve SNR levels of 20 and 40 to simulate typical acquisitions. For both SNR levels, 100 noise iterations were performed. With singular value decomposition (SVD) as the control, the synthetic DWI signals were processed by the KT algorithms. Outputs were compared to the GT metric images at each voxel by MSE, bias2, and variance.

Results

MSE-Bias2-Variance analysis illuminates the differences between algorithm outputs and GT images. For MD (Figure 2), there is minimal contribution from algorithm bias. Noise has some effect on variance at low SNR, but the magnitude of the error is negligible (~1%). For FA (Figure 3), 6KT exhibits a bias in the white matter of about 5%; its variance is lower than the other algorithms but not enough to overcome the bias. The other three algorithms perform similarly. For MK (Figure 4), 6KT has bias and, to a lesser extent, variance error. In major WM tracts, its performance is comparable to SVD. REKINDLE has the lowest MSE of the algorithms, especially in the peripheral WM. The variance is the source of its error. DWAR performed similarly to SVD in nearly all regions.

Discussion

Based on the analysis, there is no clear winner for optimal KT fitting algorithm. The MSE was within about 5% of the actual metric value indicating good algorithm performance throughout the WM. 6KT has a slight advantage in robustness to noise, but the model simplicity underfits the data resulting in greater bias. With higher SNR, the robustness of 6KT no longer matters, yet the bias persists. REKINDLE is an improvement over SVD in terms of accuracy. Conversely, it is an iterative algorithm, and its processing time was approximately a 10-fold increase over SVD. DWAR performed similarly to SVD for all metrics in nearly all regions.

Although the GT metric images appear to be representative of typical in vivo values from human subjects, there could be errors from misregistration between subjects or partial volume effects. This might affect the GT tensor values at the edges of major fiber tracts and peripheral WM regions. These effects are mitigated to some extent by synthesizing the GT signal and metric images from the tensor images but should be considered when interpreting results. This GT template and analysis approach can be extended to arbitrary sampling schemes and other diffusion models.

Acknowledgements

Project funded by the MCW-AHW funds. Data were provided by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 NIH Institutes and Centers that support the NIH Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University.

References

1. Cheung MM, Hui ES, Chan KC, Helpern JA, Qi L, Wu EX. (2009). “Does diffusion kurtosis imaging lead to better neural tissue characterization? A rodent brain maturation study.” NeuroImage 45(2):386-392.

2. Tabesh A, Jensen JH, Fieremans E, Helpern JA. (2011). “Q-space undersampled diffusional kurtosis imaging.” Proc ISMRM 19 (Abstract 1990)

3. Tax CMW, Otte WM, Viergever MA, Dijkhuizen RM, Leemans A. (2015) “REKINDLE: robust extraction of kurtosis indices with linear estimation.” MRM 73(2):794-808

4. Kuder TA, Stieltjes B, Bachert P, Semmler W, Laun FB. (2012) “Advanced fit of the diffusion kurtosis tensor by directional weighting and regularization.” MRM 67(5):1401-11

5. Geman S, Bienenstock E, Doursat R. (1992) “Neural networks and the bias/variance dilemma.” Neural Computation 4:1-58

6. Van Essen DC, Smith SM, Barch DM, Behrens TEJ, Yacoub E, Ugurbil K. (2013). “The WU-Minn Human Connectome Project: An overview.” NeuroImage 80(2013):62-79

7. Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K. “Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging.” MRM. 53(6):1432-1440.

8. Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. (2012) “FSL.” NeuroImage 62: 782-90.

Figures

Figure 1: Ground truth metric images. Mean diffusivity (left), fractional anisotropy (center), mean kurtosis (right). MD has units of mm2/s. FA and MK are unitless.

Figure 2: Mean diffusivity MSE-Bias2-Variance analysis for kurtosis tensor fitting algorithms. All images have consistent window/levels according to the colorbar on the right. Overlaid on the ground truth MD image for reference.

Figure 3: Fractional anisotropy MSE-Bias2-Variance analysis for kurtosis tensor fitting algorithms. All images have consistent window/levels according to the colorbar on the right. Overlaid on the ground truth FA image for reference.

Figure 4: Mean kurtosis MSE-Bias2-Variance analysis for kurtosis tensor fitting algorithms. All images have consistent window/levels according to the colorbar on the right. Overlaid on the ground truth MK image for reference.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3072