Deterministic Boltzmann Transport Equation Solver for Patient-Specific CT Dose Estimation: Comparison Against a Monte Carlo Benchmark for Realistic Scanner Configurations and Patient Models.

PURPOSE Epidemiological evidence suggests an increased risk of cancer related to CT scans, with children exposed to greater risk. The purpose of this work is to test the reliability of a linear Boltzmann Transport Equation (LBTE) solver for rapid and patient-specific CT dose estimation. This includes building a flexible LBTE framework for modeling modern clinical CT scanners and to validate the resulting dose maps across a range of realistic scanner configurations and patient models. METHODS In this study, computational tools were developed for modeling CT scanners, including a bowtie filter, overrange collimation, and tube current modulation. The LBTE solver requires discretization in the spatial, angular, and spectral dimensions, which may affect the accuracy of scanner modeling. To investigate these effects, this study evaluated the LBTE dose accuracy for different discretization parameters, scanner configurations, and patient models (male, female, adults, pediatric). The method used to validate the LBTE dose maps was the Monte Carlo code Geant4, which provided ground truth dose maps. LBTE simulations were implemented on a GeForce GTX 1080 graphic unit, while Geant4 was implemented on a distributed cluster of CPUs. RESULTS The agreement between Geant4 and the LBTE solver quantifies the accuracy of the LBTE, which was similar across the different protocols and phantoms. The results suggest that 18 views per rotation provides sufficient accuracy, as no significant improvement in the accuracy was observed by increasing the number of projection views. Considering this discretization, the LBTE solver average simulation time was approximately 30 seconds. However, in the LBTE solver the phantom model was implemented with a lower voxel resolution with respect to Geant4, as it is limited by the memory of the GPU. Despite this discretization the results showed a good agreement between the LBTE and Geant4, with root mean square error of the deposited energy in organs of approximately 3.5% for most of the studied configurations. CONCLUSIONS The LBTE solver is proposed as an alternative to Monte Carlo for patient-specific organ dose estimation. This study demonstrated accurate organ dose estimates for the rapid LBTE solver when considering realistic aspects of CT scanners and a range of phantom models. Future plans will combine the LBTE framework with deep-learning autosegmentation algorithms to provide near real-time patient-specific organ dose estimation.

[1]  J. Hsieh,et al.  Deterministic absorbed dose estimation in computed tomography using a discrete ordinates method. , 2015, Medical physics.

[2]  K. D. Lathrop,et al.  Numerical Solution of the Boltzmann Transport Equation , 1966 .

[3]  Firas Mourtada,et al.  Validation of a new grid-based Boltzmann equation solver for dose calculation in radiotherapy with photon beams , 2010, Physics in medicine and biology.

[4]  A. Dell'Acqua,et al.  Geant4 - A simulation toolkit , 2003 .

[5]  Kai Yang,et al.  Characterization of dynamic collimation mechanisms for helical CT scans with direct measurements , 2019, Physics in medicine and biology.

[6]  D. Brenner,et al.  Estimated risks of radiation-induced fatal cancer from pediatric CT. , 2001, AJR. American journal of roentgenology.

[7]  Adam S. Wang,et al.  A fast, linear Boltzmann transport equation solver for computed tomography dose calculation (Acuros CTD) , 2018, Medical physics.

[8]  J. Boone,et al.  CT dose index and patient dose: they are not the same thing. , 2011, Radiology.

[9]  Dose is not always what it seems: where very misleading values can result from volume CT dose index and dose length product. , 2014, Journal of the American College of Radiology : JACR.

[10]  W. Segars,et al.  4D XCAT phantom for multimodality imaging research. , 2010, Medical physics.

[11]  Daniel Kolditz,et al.  Fast on-site Monte Carlo tool for dose calculations in CT applications. , 2012, Medical physics.

[12]  A. Einstein,et al.  Estimating risk of cancer associated with radiation exposure from 64-slice computed tomography coronary angiography. , 2007, JAMA.

[13]  J. Boone,et al.  Dose equations for tube current modulation in CT scanning and the interpretation of the associated CTDIvol. , 2013, Medical physics.

[14]  Aldo Badano,et al.  Accelerating Monte Carlo simulations of photon transport in a voxelized geometry using a massively parallel graphics processing unit. , 2009, Medical physics.

[15]  John M Boone,et al.  Radiation exposure from CT scans: how to close our knowledge gaps, monitor and safeguard exposure--proceedings and recommendations of the Radiation Dose Summit, sponsored by NIBIB, February 24-25, 2011. , 2012, Radiology.

[16]  Taly Gilat Schmidt,et al.  Accuracy of patient-specific organ dose estimates obtained using an automated image segmentation algorithm , 2016, Journal of medical imaging.

[17]  Taly Gilat Schmidt,et al.  Evaluation of deep learning segmentation for rapid, patient-specific CT organ dose estimation using an LBTE solver , 2020 .

[18]  K. Bush,et al.  Dosimetric validation of Acuros XB with Monte Carlo methods for photon dose calculations. , 2011, Medical physics.

[19]  Ehsan Samei,et al.  A real-time Monte Carlo tool for individualized dose estimations in clinical CT , 2019, Physics in medicine and biology.

[20]  Xun Jia,et al.  Point/counterpoint. GPU technology is the hope for near real-time Monte Carlo dose calculations. , 2015, Medical physics.

[21]  Xun Jia,et al.  Comprehensive evaluations of cone-beam CT dose in image-guided radiation therapy via GPU-based Monte Carlo simulations. , 2014, Physics in medicine and biology.

[22]  Ehsan Samei,et al.  Monte Carlo reference data sets for imaging research: Executive summary of the report of AAPM Research Committee Task Group 195. , 2015, Medical physics.

[23]  Taly Gilat Schmidt,et al.  Technical Note: Phantom study to evaluate the dose and image quality effects of a computed tomography organ-based tube current modulation technique. , 2015, Medical physics.

[24]  Alexander Maslowski,et al.  Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part II: System modeling, scatter correction, and optimization. , 2018, Medical physics.

[25]  D. Cullen,et al.  Tables and graphs of photon interaction cross sections from 1.0 keV to 100 MeV derived from the LLL Evaluated Nuclear Data Library , 1975 .

[26]  Josh Star-Lack,et al.  Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part I: Core algorithms and validation. , 2018, Medical physics.

[27]  S. T. Perkins,et al.  Tables and graphs of electron-interaction cross sections from 10 eV to 100 GeV derived from the LLNL Evaluated Electron Data Library (EEDL), Z = 1--100 , 1991 .

[28]  M. Kalra,et al.  A Method of Rapid Quantification of Patient-Specific Organ Doses for CT Using Deep-Learning based Multi-Organ Segmentation and GPU-accelerated Monte Carlo Dose Computing. , 2019, Medical physics.