A model-constrained Monte Carlo method for blind arterial input function estimation in dynamic contrast-enhanced MRI: II. In vivo results

Accurate quantification of pharmacokinetic model parameters in tracer kinetic imaging experiments requires correspondingly accurate determination of the arterial input function (AIF). Despite significant effort expended on methods of directly measuring patient-specific AIFs in modalities as diverse as dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), dynamic positron emission tomography (PET), and perfusion computed tomography (CT), fundamental and technical difficulties have made consistent and reliable achievement of that goal elusive. Here, we validate a new algorithm for AIF determination, the Monte Carlo blind estimation (MCBE) method (which is described in detail and characterized by extensive simulations in a companion paper), by comparing AIFs measured in DCE-MRI studies of eight brain tumor patients with results of blind estimation. Blind AIFs calculated with the MCBE method using a pool of concentration-time curves from a region of normal brain tissue were found to be quite similar to the measured AIFs, with statistically significant decreases in fit residuals observed in six of eight patients. Biases between the blind and measured pharmacokinetic parameters were the dominant source of error. Averaged over all eight patients, the mean biases were +7% in K^trans, 0% in k_ep,-11% in v_p and +10% in v_e. Corresponding uncertainties (median absolute deviation from the best fit line) were 0.0043 min^-1 in K^trans, 0.0491 min^-1 in k _ep, 0.29% in v_p and 0.45% in v_e. The use of a published population-averaged AIF resulted in larger mean biases in three of the four parameters (-23% in K^trans,-22% in k_ep,-63% in v_p), with the bias in v_e unchanged, and led to larger uncertainties in all four parameters (0.0083 min^-1 in K ^trans, 0.1038 min^-1 in k_ep, 0.31% in v _p and 0.95% in v_e). When blind AIFs were calculated from a region of tumor tissue, statistically significant decreases in fit residuals were observed in all eight patients despite larger deviations of these blind AIFs from the measured AIFs. The observed