Estimating Kinetic Parameter Maps From Dynamic Contrast-Enhanced MRI Using Spatial Prior Knowledge

Abstract

Dynamic contrast-enhanced magnetic resonance (DCE-MR) imaging can be used to study microvascular structure in vivo by monitoring the abundance of an injected diffusible contrast agent over time. The resulting spatially resolved intensity-time curves are usually interpreted in terms of kinetic parameters obtained by fitting a pharmacokinetic model to the observed data. Least squares estimates of the highly nonlinear model parameters, however, can exhibit high variance and can be severely biased. As a remedy, we bring to bear spatial prior knowledge by means of a generalized Gaussian Markov random field (GGMRF). By using information from neighboring voxels and computing the maximum a posteriori solution for entire parameter maps at once, both bias and variance of the parameter estimates can be reduced thus leading to smaller root mean square error (RMSE). Since the number of variables gets very big for common image resolutions, sparse solvers have to be employed. To this end, we propose a generalized iterated conditional modes (ICM) algorithm operating on blocks instead of sites which is shown to converge considerably faster than the conventional ICM algorithm. Results on simulated DCE-MR images show a clear reduction of RMSE and variance as well as, in some cases, reduced estimation bias. The mean residual bias (MRB) is reduced on the simulated data as well as for all 37 patients of a prostate DCE-MRI dataset. Using the proposed algorithm, average computation times only increase by a factor of 1.18 (871 ms per voxel) for a Gaussian prior and 1.51 (1.12 s per voxel) for an edge-preserving prior compared to the single voxel approach (740 ms per voxel).

Dynamic Contrast Enhanced Magnetic Resonance Imaging in Rectal Cancer

“…This approach can potentially benefit from both ROI-based and pixel-by-pixel processing (Sourbron, 2010). Also, semi-automatic approaches for model-based segmentation of DCE-MRI images are currently being developed (Buonaccorsi et al, 2007;Kelm et al, 2009;Sansone et al, 2011;Schmid et al, 2006;Xiaohua et al, 2005) …”

Acceleration of dynamic fluorescence molecular tomography with principal component analysis

“…Meanwhile, compartmental modeling is a well-known method for pharmacokinetic analysis, and the pharmacokinetic parameters defined in the compartmental model can be utilized to describe the uptake and excretion of substances such as drugs, radiotracers and fluorophores in the body [6]. The boundary measurements of dynamic FMT can be transformed into three-dimensional (3-D) images of pharmacokinetic parameters through compartmental modeling [4][5][6][7], as has been done in positron emission tomography (PET) [8][9][10] and magnetic resonance imaging (MRI) [11][12][13]. Images of pharmacokinetic parameters are also known as parametric images, which can provide valuable physiological information for tumor detection [14][15][16], therapy assessment [17], and organ function evaluation [5][6][7]18].…”

An expectation-maximisation approach for simultaneous pixel classification and tracer kinetic modelling in dynamic contrast enhanced-magnetic resonance imaging

“…Kelm et al [12] proposed a block version of the iterated conditional modes (ICM) in order to reduce the long computational time of GMRF. Moreover in their approach, fewer assumptions on the underlying distribution are needed and the accuracy of the estimates is improved by using neighbourhood information.…”