Tissue Partial Volume Correction of Perfusion Maps in Dynamic Susceptibility Contrast MRI
André Ahlgren1, Ronnie Wirestam1, Freddy Ståhlberg1,2,3, and Linda Knutsson1

1Department of Medical Radiation Physics, Lund University, Lund, Sweden, 2Department of Diagnostic Radiology, Lund University, Lund, Sweden, 3Lund University Bioimaging Center, Lund University, Lund, Sweden

Synopsis

Partial volume effects (PVE) can significantly affect parameter estimates in perfusion MRI. In contrast to arterial spin labeling (ASL), the impact of PVEs in dynamic susceptibility contrast MRI (DSC-MRI) has not yet been well established. In this work, we assess and compare partial volume correction (PVC) of DSC-MRI and ASL data in 20 healthy subjects. PVC reduced the tissue volume dependence of perfusion estimates in DSC-MRI and ASL. White matter perfusion maps were of higher quality for DSC-MRI. However, for PVC of DSC-MRI we used several assumptions which need further evaluation.

Purpose

To assess the possibility and the potential impact of partial volume correction (PVC) in dynamic susceptibility contrast MRI (DSC-MRI).

Introduction

Tissue partial volume effects (PVEs) constitute an important source of bias in measurements of absolute brain perfusion. In particular, low spatial resolution tends to result in under- and overestimation of cerebral blood flow (CBF) estimates in gray matter (GM) and white matter (WM), respectively. Thus, disentangling the effects of tissue volume and perfusion is important, especially in diseases in which atrophy and perfusion alterations are common.

The impact of tissue PVEs in perfusion measurements with DSC-MRI has not previously been thoroughly investigated. In this work, two PVC methods were applied to DSC-MRI data, and, for comparison, to the corresponding arterial spin labeling (ASL) maps.

Methods

The DSC-MRI signal was assumed to be given by $$$S(t)=S_0e^{-t_E\cdot r_2^*\cdot c(t)}$$$, where $$$S_0$$$ is the baseline signal, $$$t_E$$$ is the echo time, $$$r_2^*$$$ is the transverse relaxivity of the contrast agent, and $$$c(t)$$$ is the tracer concentration. The signal can also be described as the sum of signal contributions from different components according to $$$S(t)=S_0\sum_ip_{v,i}e^{-t_E\cdot r_2^*\cdot c_i(t)}$$$, where $$$p_{v,i}$$$ is the partial volume of component $$$i$$$, $$$c_i(t)$$$ is the tracer concentration in component $$$i$$$, and $$$i=$$$ {GM,WM,CSF}. By combining the two equations above, acknowledging that cerebrospinal fluid (CSF) is non-perfused, and that partial volumes sum up to 1, Taylor expansion to the first order yields: $$$c(t)\approx p_{v,GM}\cdot c_{GM}(t)+p_{v,WM}\cdot c_{WM}(t)$$$. Preliminary simulations with plausible assumptions about $$$t_E$$$, $$$r_2^*$$$ and $$$c(t)$$$ indicate that a first-order Taylor expansion is not unreasonable (data not shown). Still, when second- and higher-order terms are non-negligible, this approximation fails. By further assuming similar temporal tracer-kinetic patterns in GM and WM we can assume that the measured perfusion $$$f\approx p_{v,GM}\cdot f_{GM}+p_{v,WM}\cdot f_{WM}$$$, and this can be used for post-quantification PVC (this relationship was also used for ASL). Here, $$$f_{GM}$$$ denotes GM perfusion and $$$f_{WM}$$$ denotes WM perfusion.

Twenty healthy volunteers underwent prebolus DSC-MRI1, pseudo-continuous ASL, and quantitative T1-based (IR Look-Locker) partial volume (PV) estimation using the fractional signal approach2. DSC-MRI and PV data had a voxel size of 1.72×1.72×5 mm3 and ASL data 2.3×2.3×5 mm3. (See Ref. 3 for details on protocols.)

PVC was performed using (i) linear regression4 (LR) and (ii) modified least trimmed squares5 (mLTS), with circular kernels (diameter 5 pixels), yielding partial-volume-free GM and WM perfusion maps. Uncorrected and PV-corrected CBF maps were analyzed, where GM and WM ROIs were extracted from the PV estimates ($$$p_{v,i}=$$$ 10–20%,...,90–100%). Visual inspection and calculation of the relative range (rR) of CBF values over the range of partial volumes (for no PVC, mLTS and LR) allowed for assessment of the degree of volume dependence of estimated CBF. Spuriously high values from large vessels were removed by masking out voxels with CBV>10 ml/100g.

Results

Example PV and PVC maps are shown for DSC-MRI (Fig. 1) and ASL (Fig. 2). (Note that PVC CBF maps are masked at PV>10%.) LR yielded more apparent smoothing in the PVC maps than mLTS, and WM perfusion maps were less noisy in DSC-MRI compared to ASL. Figure 3 shows results of the ROI analysis. PVC decreased the volume dependence of CBF estimates for DSC-MRI as well as ASL, and for both tissue types. Notably, for DSC-MRI, the LR method decreased the rR from 44.4% (for uncorrected data) to 4.2% in GM and from 55.3% to 14.2% in WM. As expected, the PVEs were larger in the ASL data, with a rR of 60-80% for uncorrected data.

Discussion and Conclusion

In this work, a linear approximation of the DSC-MRI signal was applied to allow for post-quantification PVC of perfusion estimates. PVEs were clearly identified in DSC-MRI data, and PVC decreased the volume dependence of CBF estimates. However, the first-order Taylor approximation needs to be further evaluated, and it might be preferable to expand the concept by, for example, nonlinear regression or, alternatively, by modelling the tissue signal as a combination of nonlinear tracer responses from GM and WM (i.e., modelling the PVE within the quantification process).

The mLTS algorithm performed less convincing than LR with regard to removing volume dependence, although it preserved more detail in the corrected maps due to the use of adaptive matrix trimming. However, the low volume dependence in LR could, at least partially, be an effect of the blurring. PVC can also be applied to CBV estimates, and therefore, in principle, allow for separate GM and WM MTT estimates. In conclusion, PVEs significantly affected CBF estimates in DSC-MRI, and application of PVC should be considered when absolute values are of interest.

Acknowledgements

This work was supported by the Swedish Research Council.

References

[1] Knutsson et al. MRM 2014;72:996-1006. [2] Shin et al. Neuroimage 2010;52:1347-1354. [3] Lindgren et al. MAGMA 2014;27:487-499. [4] Asllani et al. MRM 2008;60:1362-1371. [5] Liang et al. MRM 2013;69:531-537.

Figures

Figure 1. DSC-MRI parameter maps in one slice for one subject. (a) CBF map; (b) GM PV; (c) with mLTS; (d) with LR; (e) WM PV; (f) with mLTS; (g) with LR. Displayed ranges are 0-100 ml/100g/min (a,c,d), 0-30 ml/100g/min (f,g) and 0-1 (b,e).

Figure 2. ASL parameter maps in one slice for one subject (same as Fig. 1). (a) CBF map; (b) GM PV; (c) with mLTS; (d) with LR; (e) WM PV; (f) with mLTS; (g) with LR. Displayed ranges are 0-100 ml/100g/min (a,c,d), 0-30 ml/100g/min (f,g) and 0-1 (b,e).

Figure 3. Uncorrected and PV-corrected CBF estimates in GM and WM, for DSC-MRI and ASL, using LR and mLTS, for all 20 subjects. (a) DSC-MRI GM, (b) ASL GM, (c) DSC-MRI WM, (d) ASL WM.



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