Deformation and resolution issues in partial volume correction of 2D arterial spin labeling data
Jan Petr1, Henri JMM Mutsaerts2, Enrico De Vita3,4, Jens Maus1, Jörg van den Hoff1, and Iris Asllani5

1Institute of Radiopharmaceutical Cancer Research, Helmholtz-Zentrum Dresden-Rossendorf, Dresden, Germany, 2Sunnybrook Research Institute, Toronto, ON, Canada, 3Lysholm Department of Neuroradiology, National Hospital for Neurology and Neurosurgery, London, United Kingdom, 4Academic Neuroradiological Unit, Department of Brain Repair and Rehabilitation, UCL Institute of Neurology, London, United Kingdom, 5Department of Biomedical Engineering, Rochester Institute of Technology, Rochester, NY, United States

Synopsis

Partial volume (PV) effects are a well-recognized confounder in arterial spin labeling due to its limited spatial resolution. Several algorithms exist to correct for these errors. Nevertheless, PV-correction is rarely used, mainly because the PV maps obtained from segmented T1-weighted images are regarded as not being sufficiently reliable when transformed into ASL space. Here, we show the impact of spatial deformation and resolution in the PV-maps used for PV-correction in the calculation of mean total gray matter (GM) cerebral blood flow (CBF). We also show how the deformations affect the calculation of PV-uncorrected mean GM CBF.

Purpose

Partial volume (PV) effects are a well-recognized challenge in arterial spin labeling (ASL), and can have profound effects on the accuracy of the ASL-derived cerebral blood flow (CBF) maps1,2,3. Several algorithms exist to correct for these errors4,5,6,7. A limitation of the current strategies is the use of PV-maps from segmented T1-weighted images (pGM), which differ geometrically from the ASL acquisition8. This study aims to quantify the impact of common errors in PV maps caused by mismatches in geometric distortion and spatial resolution between the low-resolution ASL and high-resolution T1-weighted acquisitions. Furthermore, we show how PV-correction relate to the commonly used calculation of the mean gray matter (GM) CBF within a GM mask, which is essentially a form of PV-correction. Lastly, the CBF maps were smoothed to simulate a larger acquisition point-spread-function (PSF) and sensitivity to motion.

Methods

Three methods to estimate the mean GM CBF (within a pGM-mask obtained by thresholding the pGM map) were compared: 1) GM-Mask: mean CBF within a pGM mask; 2) GM-Weighted: within the pGM mask, the mean CBF is divided by the mean pGM; 3) PVEC: within the pGM mask, the mean of the PV-corrected GM-CBF values (the Asllani’s method with a 3x3x1 kernel3). Both 2D-EPI pseudo-continuous ASL9 and T1-weighted (T1w) images were acquired for a young healthy volunteer (Table 1). The ASL and M0 images were motion corrected and the perfusion-weighted difference (PWI) was calculated. The PWI and M0 were upsampled to the T1w space and the T1w volume was segmented into GM/WM using SPM12. To estimate the deformations between ASL and T1w volumes they were co-registered using four different methods: A) T1w-Rigid: rigid transformation (between T1w and M0); B) GM-Rigid: rigid transformation (between PWI-GM); C) Nonlin: SPM’s affine transformation followed by nonlinear transformation10 (between PWI-GM); D) DARTEL: creating a DARTEL11 template from the PWI and GM image (Figure 1). Each of these transformations was separately applied to the original GM/WM maps to create four deformed volumes, which were subsequently used to simulate the deformation between the ASL-image and T1-based PV-maps. The results were downsampled using an anisotropic Gaussian kernel identical to the acquisition resolution. The downsampled GM/WM maps from method D (DARTEL) were used to generate a simulated CBF (sCBF) image by assuming a uniform GM-CBF of 80 ml/min/100g, GM/WM ratio of 3, and SD-4 Gaussian noise. The mean CBF in sCBF was evaluated using the four different downsampled PV-maps for the three PV-correction methods. A second set of sCBF maps was created from method D by increasing FWHM from 3x3x7mm3 to 7x7x7mm3. The mean CBF was evaluated using PV-maps with FWHM 3x3x7mm3 to simulate the effect of having sharp PV-maps and blurred CBF (Figure 2). The globus pallidus and thalamus were excluded from the analysis due to T1-segmentation issues.

Results

To assess the differences across the 4 methods, for each PV-correction algorithm, we computed: the mean relative voxel-wise difference between all pGMs and the ground-truth pGM (D), the relative error in the mean GM-CBF, and the mean relative error in the GM-CBF calculated on a local neighborhood of 15x15x7 mm3. Deformation methods A-D (Figure 3), and FWHMs from 3x3x7mm3 to 7x7x7mm3 for the method D (Figure 4) were compared. For a typical threshold (pGM>0.7), the mean pixel-wise error in pGM-maps was between 10-15%. Without PV-correction, global and local CBF errors were 16-21% for all methods, including D. With both PV-corrections, the local error was under 8%. A local error of 4% with perfect PV-maps was achieved for GM-Weighted PV-correction. The global CBF error was <5%, however, this is interpreted carefully as a uniform GM-CBF value was used. The CBF-map blurring caused 3-13% error in pGM. This added 2-10% to the error in the mean global/local CBF calculation without PV-correction. This effect persisted with PV-correction regardless of deformation for the global error. However, the blurring did not affect the local CBF error for the deformations methods A-C with PV-correction.

Conclusions

Without any PV-correction, the ASL-T1w deformations had only a small effect on the calculated mean CBF values, however, blurring further underestimated the CBF. Blurring of data caused by acquisition, motion or post-processing thus presents an important factor for the mean CBF evaluation. This needs to be taken into account in multi-center ASL studies, especially for 3D sequences with large PSF and in patients that are prone to move. Even on a local scale, both PV-correction methods decreased the errors in mean CBF calculation. The use of PV-correction did always improve the results even in the presence of deformations and data blurring, despite a decrease of the accuracy of the PV-maps.

Acknowledgements

The authors would like to acknowledge networking support by the COST Action BM1103. Part of this work was undertaken at UCLH/UCL who received a proportion of funding from the Department of Health’s NIHR Biomedical Research Centres funding scheme.

References

1. Asllani, I., Habeck, C., Borogovac, A., Brown, T. R., Brickman, A. M., and Stern, Y. Separating function from structure in perfusion imaging of the aging brain. Human Brain Mapping, 30(9), 2927–35, 2009.

2. Chen, J. J., Rosas, H. D., and Salat, D. H. Age-associated reductions in cerebral blood flow are independent from regional atrophy. NeuroImage, 55(2), 468–78, 2011.

3. Mutsaerts, H. J., Richard, E., Heijtel, D. F., van Osch, M. J., Majoie, C. B., and Nederveen, A. J. Gray matter contamination in arterial spin labeling white matter perfusion measurements in patients with dementia. NeuroImage: Clinical,4:139-44, 2014.

4. Asllani, I., Borogovac, A., & Brown, T. R. Regression algorithm correcting for partial volume effects in arterial spin labeling MRI. Magnetic Resonance in Medicine, 60(6), 1362–71, 2008.

5. Liang, X., Connelly, A., and Calamante, F. Improved partial volume correction for single inversion time arterial spin labeling data. Magnetic Resonance in Medicine, 69(2), 531–7, 2013.

6. Chappell, M. A., Groves, A. R., MacIntosh, B. J., Donahue, M. J., Jezzard, P., and Woolrich, M. W. Partial volume correction of multiple inversion time arterial spin labeling MRI data. Magnetic Resonance in Medicine, 65(4), 1173–83, 2011.

7. Petr, J., Schramm, G., Hofheinz, F., Langner, J., and van den Hoff, J. Partial volume correction in arterial spin labeling using a Look-Locker sequence. Magnetic Resonance in Medicine, 70(6), 1535–43, 2013.

8. Chen, N.-K., and Wyrwicz, A. M. Optimized distortion correction technique for echo planar imaging. Magnetic Resonance in Medicine, 45:525–8, 2001.

9. Alsop, D. C., Detre, J.A., Golay, X., Günther, M., Hendrikse, J., Hernandez-Garcia, L., Hanzhang, L., MacIntosh, B.J., Parkes, L. M., Smits, M., van Osch, M. J. P., Wang, D. J. J., Wong, E. C., Zaharchuk, G. Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: A consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia. Magnetic Resonance in Medicine, 73(1):102-16, 2015.

10. Ashburner, J., and Friston, K.J. Nonlinear spatial normalization using basis functions. Human Brain Mapping, 7(4):254-266, 1999.

11. Ashburner, J. A fast diffeomorphic image registration algorithm. Neuroimage, 38(1):95-113, 2007.

Figures

Table 1. Acquisition parameters of the 3D T1-weighted and arterial spin labeling (ASL) sequences. A 3T MRI scanner with an 8-channel head coil was used for all measurements.

Figure 1. ASL images were motion-corrected, upsampled, T1-weighted image was segmented, and T1-weighted and ASL were co-registered (above). The inverse transformations were applied to the segmented GM and WM. Results were downsampled to create PV-maps. CBF image was generated from the reference (DARTEL) GM and WM PV-maps.

Figure 2. The original PWI image (shown for comparison only) and the downsampled GM maps. Deformations are visible when different co-registration methods were used (upper row, PSF 3x3x7mm3). The DARTEL GM map is displayed for downsampling PSFs between 4x4x7 mm3 and 7x7x7 mm3 (lower row).

Figure 3. Comparison of PV-correction methods with different PV-map deformations. The mean relative GM error, the local and global CBF errors are shown for different types of PV-correction (GM-Mask, GM-Weighted, and PVEC).

Figure 4. Comparison of PV-corrections methods with different FWHMs used for the PV-map downsampling. The mean relative GM error, the local and global relative CBF errors are shown. All images for this comparison were created based on the registration D (DARTEL).



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