Purpose
To present a method for estimating the local arterial input function (AIF) within a dynamic contrast‐enhanced MRI scan, based on the alternating minimization with model (AMM) method.
Materials and Methods
This method clusters a subset of data into representative curves, which are then input to the AMM algorithm to return a parameterized AIF and pharmacokinetic parameters. Computer simulations are used to investigate the accuracy with which the AMM is able to estimate the true AIF as a function of the input tissue curves.
Results
Simulations show that a power law relates uncertainty in kinetic parameters and SNR and heterogeneity of the input. Kinetic parameters calculated with the measured AIF are significantly different from those calculated with either a global (P < 0.005) or a local input function (P = 0.0). The use of local AIFs instead of measured AIFs yield mean lesion‐averaged parameter changes: Ktrans: +24% [+15%, +70%], kep: +13% [−36%, +300%]. Globally estimated input functions yield mean lesion‐averaged changes: Ktrans: +9% [−38%, +65%], kep: +13% [−100%, +400%]. The observed improvement in fit quality with local AIFs was found to be significant when additional free parameters were accounted for using the Akaike information criterion.
Conclusion
Local AIFs result in significantly different kinetic parameter values. The statistically significant improvement in fit quality suggests that changes in parameter estimates using local AIFs reflect differences in underlying tissue physiology. J. Magn. Reson. Imaging 2010;32:924–934. © 2010 Wiley‐Liss, Inc.