Accuracy and Performance of Functional Parameter Estimation Using a Novel Numerical Optimization Approach for GPU-Based Kinetic Compartmental Modeling

Quantitative kinetic parameters derived from dynamic contrast-enhanced (DCE) data are dependent on signal measurement quality and choice of pharmacokinetic model. However, the fundamental optimization analysis method is equally important and its impact on pharmacokinetic parameters has been mostly overlooked. We examine the effects of those choices on accuracy and performance of parameter estimation using both computer processing unit and graphical processing unit (GPU) numerical optimization implementations and evaluate the improvements offered by a novel optimization approach. A test framework was developed where experimentally derived population-average arterial input function and randomly sampled parameter sets {Ktrans, Kep, Vb, τ} were used to generate known tissue curves. Five numerical optimization algorithms were evaluated: sequential quadratic programming, downhill simplex (Nelder–Mead), pattern search, simulated annealing, and differential evolution. This was combined with various objective function implementation details: delay approximation, discretization and varying sampling rates. Then, impact of noise and CPU/GPU implementation was tested for speed and accuracy. Finally, the optimal method was compared to conventional implementation as applied to clinical DCE computed tomography. Nelder–Mead, differential evolution and sequential quadratic programming produced good results on clean and noisy input data outperforming simulated annealing and pattern search in terms of speed and accuracy in the respective order of 10−8%, 10−7%, and ×10−6%). A novel approach for DCE numerical optimization (infinite impulse response with fractional delay approximation) was implemented on GPU for speed increase of at least 2 orders of magnitude. Applied to clinical data, the magnitude of overall parameter error was <10%.

[1]  David Jaffray,et al.  Automated voxel-based analysis of volumetric dynamic contrast-enhanced CT data improves measurement of serial changes in tumor vascular biomarkers. , 2015, International journal of radiation oncology, biology, physics.

[2]  H Wang,et al.  GPU-accelerated voxelwise hepatic perfusion quantification. , 2012, Physics in medicine and biology.

[3]  D. Jaffray,et al.  Implementation and characterization of a 320-slice volumetric CT scanner for simulation in radiation oncology. , 2009, Medical physics.

[4]  Thomas E. Yankeelov,et al.  The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge , 2016, Tomography.

[5]  M. Knopp,et al.  Estimating kinetic parameters from dynamic contrast‐enhanced t1‐weighted MRI of a diffusable tracer: Standardized quantities and symbols , 1999, Journal of magnetic resonance imaging : JMRI.

[6]  Charles Audet,et al.  Analysis of Generalized Pattern Searches , 2000, SIAM J. Optim..

[7]  Catherine Coolens,et al.  Unified platform for multimodal voxel-based analysis to evaluate tumour perfusion and diffusion characteristics before and after radiation treatment evaluated in metastatic brain cancer. , 2019, The British journal of radiology.

[8]  Unto K. Laine,et al.  Splitting the unit delay [FIR/all pass filters design] , 1996, IEEE Signal Process. Mag..

[9]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[10]  A. Jackson,et al.  Experimentally‐derived functional form for a population‐averaged high‐temporal‐resolution arterial input function for dynamic contrast‐enhanced MRI , 2006, Magnetic resonance in medicine.

[11]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[12]  Catherine Coolens,et al.  Quantitative Imaging in Radiation Oncology: An Emerging Science and Clinical Service. , 2015, Seminars in radiation oncology.

[13]  C Coolens,et al.  Development of a dynamic quality assurance testing protocol for multisite clinical trial DCE-CT accreditation. , 2013, Medical physics.

[14]  William H. Press,et al.  Numerical recipes: the art of scientific computing, 3rd Edition , 2007 .

[15]  Catherine Coolens,et al.  Comparison of Voxel-Wise Tumor Perfusion Changes Measured With Dynamic Contrast-Enhanced (DCE) MRI and Volumetric DCE CT in Patients With Metastatic Brain Cancer Treated with Radiosurgery , 2016, Tomography.