A fast, linear Boltzmann transport equation solver for computed tomography dose calculation (Acuros CTD)

Monte Carlo methods for medical imaging research

In radiation-based medical imaging research, computational modeling methods are used to design and validate imaging systems and post-processing algorithms. Monte Carlo methods are widely used for the computational modeling as they can model the systems accurately and intuitively by sampling interactions between particles and imaging subject with known probability distributions. This article reviews the physics behind Monte Carlo methods, their applications in medical imaging, and available MC codes for medical imaging research. Additionally, potential research areas related to Monte Carlo for medical imaging are discussed.

Is deep learning-enabled real-time personalized CT dosimetry feasible using only patient images as input?

PURPOSE

To propose a novel deep-learning based dosimetry method that allows quick and accurate estimation of organ doses for individual patients, using only their computed tomography (CT) images as input.

METHODS

Despite recent advances in medical dosimetry, personalized CT dosimetry remains a labour-intensive process. Current state-of-the-art methods utilize time-consuming Monte Carlo (MC) based simulations for individual organ dose estimation in CT. The proposed method uses conditional generative adversarial networks (cGANs) to substitute MC simulations with fast dose image generation, based on image-to-image translation. The pix2pix architecture in conjunction with a regression model was utilized for the generation of the synthetic dose images. The lungs, heart, breast, bone and skin were manually segmented to estimate and compare organ doses calculated using both the original and synthetic dose images, respectively.

RESULTS

The average organ dose estimation error for the proposed method was 8.3% and did not exceed 20% for any of the organs considered. The performance of the method in the clinical environment was also assessed. Using segmentation tools developed in-house, an automatic organ dose calculation pipeline was set up. Calculation of organ doses for heart and lung for each CT slice took about 2 s.

CONCLUSIONS

This work shows that deep learning-enabled personalized CT dosimetry is feasible in real-time, using only patient CT images as input.

A fully automated machine learning-based methodology for personalized radiation dose assessment in thoracic and abdomen CT

PURPOSE

To develop a machine learning-based methodology for patient-specific radiation dosimetry in thoracic and abdomen CT.

METHODS

Three hundred and thirty-one thoracoabdominal radiotherapy-planning CT examinations with the respective organ/patient contours were collected retrospectively for the development and validation of segmentation 3D-UNets. Moreover, 97 diagnostic thoracic and 89 diagnostic abdomen CT examinations were collected retrospectively. For each of the diagnostic CT examinations, personalized MC dosimetry was performed. The data derived from MC simulations along with the respective CT data were used for the training and validation of a dose prediction deep neural network (DNN). An algorithm was developed to utilize the trained models and perform patient-specific organ dose estimates for thoracic and abdomen CT examinations. The doses estimated with the DNN were compared with the respective doses derived from MC simulations. A paired t-test was conducted between the DNN and MC results. Furthermore, the time efficiency of the proposed methodology was assessed.

RESULTS

The mean percentage differences (range) between DNN and MC dose estimates for the lungs, liver, spleen, stomach, and kidneys were 7.2 % (0.2-24.1 %), 5.5 % (0.4-23.0 %), 7.9 % (0.6-22.3 %), 6.9 % (0.0-23.0 %) and 6.7 % (0.3-22.6 %) respectively. The differences between DNN and MC dose estimates were not significant (p-value = 0.12). Moreover, the mean processing time of the proposed workflow was 99 % lower than the respective time needed for MC-based dosimetry.

CONCLUSIONS

The proposed methodology can be used for rapid and accurate patient-specific dosimetry in chest and abdomen CT.

Real-time, acquisition parameter-free voxel-wise patient-specific Monte Carlo dose reconstruction in whole-body CT scanning using deep neural networks

OBJECTIVE

We propose a deep learning-guided approach to generate voxel-based absorbed dose maps from whole-body CT acquisitions.

METHODS

The voxel-wise dose maps corresponding to each source position/angle were calculated using Monte Carlo (MC) simulations considering patient- and scanner-specific characteristics (SP_MC). The dose distribution in a uniform cylinder was computed through MC calculations (SP_uniform). The density map and SP_uniform dose maps were fed into a residual deep neural network (DNN) to predict SP_MC through an image regression task. The whole-body dose maps reconstructed by the DNN and MC were compared in the 11 test cases scanned with two tube voltages through transfer learning with/without tube current modulation (TCM). The voxel-wise and organ-wise dose evaluations, such as mean error (ME, mGy), mean absolute error (MAE, mGy), relative error (RE, %), and relative absolute error (RAE, %), were performed.

RESULTS

The model performance for the 120 kVp and TCM test set in terms of ME, MAE, RE, and RAE voxel-wise parameters was  - 0.0302 ± 0.0244 mGy, 0.0854 ± 0.0279 mGy,  - 1.13 ± 1.41%, and 7.17 ± 0.44%, respectively. The organ-wise errors for 120 kVp and TCM scenario averaged over all segmented organs in terms of ME, MAE, RE, and RAE were  - 0.144 ± 0.342 mGy, and 0.23 ± 0.28 mGy,  - 1.11 ± 2.90%, 2.34 ± 2.03%, respectively.

CONCLUSION

Our proposed deep learning model is able to generate voxel-level dose maps from a whole-body CT scan with reasonable accuracy suitable for organ-level absorbed dose estimation.

CLINICAL RELEVANCE STATEMENT

We proposed a novel method for voxel dose map calculation using deep neural networks. This work is clinically relevant since accurate dose calculation for patients can be carried out within acceptable computational time compared to lengthy Monte Carlo calculations.

KEY POINTS

• We proposed a deep neural network approach as an alternative to Monte Carlo dose calculation. • Our proposed deep learning model is able to generate voxel-level dose maps from a whole-body CT scan with reasonable accuracy, suitable for organ-level dose estimation. • By generating a dose distribution from a single source position, our model can generate accurate and personalized dose maps for a wide range of acquisition parameters.

Reduced Chest Computed Tomography Scan Length for Patients Positive for Coronavirus Disease 2019: Dose Reduction and Impact on Diagnostic Utility

METHODS

This study used the Personalized Rapid Estimation of Dose in CT (PREDICT) tool to estimate patient-specific organ doses from CT image data. The PREDICT is a research tool that combines a linear Boltzmann transport equation solver for radiation dose map generation with deep learning algorithms for organ contouring. Computed tomography images from 74 subjects in the Medical Imaging Data Resource Center-RSNA International COVID-19 Open Radiology Database data set (chest CT of adult patients positive for COVID-19), which included expert annotations including "infectious opacities," were analyzed. First, the full z-scan length of the CT image data set was evaluated. Next, the z-scan length was reduced from the left hemidiaphragm to the top of the aortic arch. Generic dose reduction based on dose length product (DLP) and patient-specific organ dose reductions were calculated. The percentage of infectious opacities excluded from the reduced z-scan length was used to quantify the effect on diagnostic utility.

RESULTS

Generic dose reduction, based on DLP, was 69%. The organ dose reduction ranged from approximately equal to 18% (breasts) to approximately equal to 64% (bone surface and bone marrow). On average, 12.4% of the infectious opacities were not included in the reduced z-coverage, per patient, of which 5.1% were above the top of the arch and 7.5% below the left hemidiaphragm.

CONCLUSIONS

Limiting z-scan length of chest CTs reduced radiation dose without significantly compromising diagnostic utility in COVID-19 patients. The PREDICT demonstrated that patient-specific organ dose reductions varied from generic dose reduction based on DLP.

Real-time estimation of patient-specific dose distributions for medical CT using the deep dose estimation

PURPOSE

With the rising number of computed tomography (CT) examinations and the trend toward personalized medicine, patient-specific dose estimates are becoming more and more important in CT imaging. However, current approaches are often too slow or too inaccurate to be applied routinely. Therefore, we propose the so-called deep dose estimation (DDE) to provide highly accurate patient dose distributions in real time METHODS: To combine accuracy and computational performance, the DDE algorithm uses a deep convolutional neural network to predict patient dose distributions. To do so, a U-net like architecture is trained to reproduce Monte Carlo simulations from a two-channel input consisting of a CT reconstruction and a first-order dose estimate. Here, the corresponding training data were generated using CT simulations based on 45 whole-body patient scans. For each patient, simulations were performed for different anatomies (pelvis, abdomen, thorax, head), different tube voltages (80 kV, 100 kV, 120 kV), different scan trajectories (circle, spiral), and with and without bowtie filtration and tube current modulation. Similar simulations were performed using a second set of eight whole-body CT scans from the Visual Concept Extraction Challenge in Radiology (Visceral) project to generate testing data. Finally, the DDE algorithm was evaluated with respect to the generalization to different scan parameters and the accuracy of organ dose and effective dose estimates based on an external organ segmentation.

RESULTS

DDE dose distributions were quantified in terms of the mean absolute percentage error (MAPE) and a gamma analysis with respect to the ground truth Monte Carlo simulation. Both measures indicate that DDE generalizes well to different scan parameters and different anatomical regions with a maximum MAPE of 6.3% and a minimum gamma passing rate of 91%. Evaluating the organ dose values for all organs listed in the International Commission on Radiological Protection (ICRP) recommendation, shows an average error of 3.1% and maximum error of 7.2% (bone surface).

CONCLUSIONS

The DDE algorithm provides an efficient approach to determine highly accurate dose distributions. Being able to process a whole-body CT scan in about 1.5 s, it provides a valuable alternative to Monte Carlo simulations on a graphics processing unit (GPU). Here, the main advantage of DDE is that it can be used on top of any existing Monte Carlo code such that real-time performance can be achieved without major adjustments. Thus, DDE opens up new options not only for dosimetry but also for scan and protocol optimization.

Technical note: Evaluation of a V‐Net autosegmentation algorithm for pediatric CT scans: Performance, generalizability and application to patient‐specific CT dosimetry

PURPOSE

This study developed and evaluated a fully convolutional network (FCN) for pediatric CT organ segmentation and investigated the generalizability of the FCN across image heterogeneities such as CT scanner model protocols and patient age. We also evaluated the autosegmentation models as part of a software tool for patient-specific CT dose estimation.

METHODS

A collection of 359 pediatric CT datasets with expert organ contours were used for model development and evaluation. Autosegmentation models were trained for each organ using a modified FCN 3D V-Net. An independent test set of 60 patients was withheld for testing. To evaluate the impact of CT scanner model protocol and patient age heterogeneities, separate models were trained using a subset of scanner model protocols and pediatric age groups. Train and test sets were split to answer questions about the generalizability of pediatric FCN autosegmentation models to unseen age groups and scanner model protocols, as well as the merit of scanner model protocol or age-group-specific models. Finally, the organ contours resulting from the autosegmentation models were applied to patient-specific dose maps to evaluate the impact of segmentation errors on organ dose estimation.

RESULTS

Results demonstrate that the autosegmentation models generalize to CT scanner acquisition and reconstruction methods which were not present in the training dataset. While models are not equally generalizable across age groups, age-group-specific models do not hold any advantage over combining heterogeneous age groups into a single training set. Dice similarity coefficient (DSC) and mean surface distance results are presented for 19 organ structures, for example, median DSC of 0.52 (duodenum), 0.74 (pancreas), 0.92 (stomach), and 0.96 (heart). The FCN models achieve a mean dose error within 5% of expert segmentations for all 19 organs except for the spinal canal, where the mean error was 6.31%.

CONCLUSIONS

Overall, these results are promising for the adoption of FCN autosegmentation models for pediatric CT, including applications for patient-specific CT dose estimation.

Validation of a deterministic linear Boltzmann transport equation solver for rapid CT dose computation using physical dose measurements in pediatric phantoms

PURPOSE

The risk of inducing cancer to patients undergoing CT examinations has motivated efforts for CT dose estimation, monitoring, and reduction, especially among pediatric population. The method investigated in this study is Acuros CTD (Varian Medical Systems, Palo Alto, CA), a deterministic linear Boltzmann transport equation (LBTE) solver aimed at generating rapid and reliable dose maps of CT exams. By applying organ contours, organ doses can also be obtained, thus patient-specific organ dose estimates can be provided. This study experimentally validated Acuros against measurements performed on a clinical CT system using a range of physical pediatric anthropomorphic phantoms and acquisition protocols.

METHODS

The study consisted of (1) the acquisition of dose measurements on a clinical CT scanner through thermoluminescent dosimeters (TLDs), and (2) the modeling in the Acuros platform of the measurement set up, which includes the modeling of the CT scanner and of the anthropomorphic phantoms. For the measurements, 1-year-old, 5-year-old, and 10-year-old anthropomorphic phantoms of the CIRS ATOM family were used. TLDs were placed in selected organ locations such as stomach, liver, lungs, and heart. The pediatric phantoms were scanned helically with the GE Discovery 750 HD clinical scanner for several examination protocols. For the simulations in Acuros, scanner-specific input, such as bowtie filters, overrange collimation, and tube current modulation schemes, were modeled. These scanner complexities were implemented by defining discretized X-ray beams whose spectral distribution, defined in Acuros by only six energy bins, varied across fan angle, cone angle, and slice position. The images generated during the CT acquisitions were used to create the geometrical models, by applying thresholding algorithms and assigning materials to the HU values. The TLDs were contoured in the phantom models as sensitive cylindrical volumes at the locations selected for dosimeters placement, to provide dose estimates, in terms of dose per unit photon. To compare measured doses with dose estimates, a calibration factor was derived from the CTDI_vol displayed by the scanner, to account for the number of photons emitted by the X-ray tube during the procedure.

RESULTS

The differences of the measured and estimated doses, in terms of absolute % errors, were within 13% for 153 TLD locations, with an error of 17% at the stomach for one study with the 10-year-old phantom. Root-mean-squared-errors (RMSE) across all TLD locations for all configurations were in the range of 3%-8%, with Acuros providing dose estimates in a time range of a few seconds up to 2 min.

CONCLUSIONS

An overall good agreement between measurements and simulations was achieved, with average RMSE of 6% across all cases. The results demonstrate that Acuros can model a specific clinical scanner despite the required discretization in spatial and energy domains. The proposed deterministic tool has the potential to be part of a near real-time individualized dosimetry monitoring system for CT applications, providing patient-specific organ dose estimates.

A GPU-accelerated framework for rapid estimation of scanner-specific scatter in CT for virtual imaging trials

Virtual imaging trials (VITs), defined as the process of conducting clinical imaging trials using computer simulations, offer a time- and cost-effective alternative to traditional imaging trials for CT. The clinical potential of VITs hinges on the realism of simulations modeling the image acquisition process, where the accurate scanner-specific simulation of scatter in a time-feasible manner poses a particular challenge. To meet this need, this study proposes, develops, and validates a rapid scatter estimation framework, based on GPU-accelerated Monte Carlo (MC) simulations and denoising methods, for estimating scatter in single source, dual-source, and photon-counting CT. A CT simulator was developed to incorporate parametric models for an anti-scatter grid and a curved energy integrating detector with an energy-dependent response. The scatter estimates from the simulator were validated using physical measurements acquired on a clinical CT system using the standard single-blocker method. The MC simulator was further extended to incorporate a pre-validated model for a PCD and an additional source-detector pair to model cross scatter in dual-source configurations. To estimate scatter with desirable levels of statistical noise using a manageable computational load, two denoising methods using a (1) convolutional neural network and an (2) optimized Gaussian filter were further deployed. The viability of this framework for clinical VITs was assessed by integrating it with a scanner-specific ray-tracer program to simulate images for an image quality (Mercury) and an anthropomorphic phantom (XCAT). The simulated scatter-to-primary ratios agreed with physical measurements within 4.4% ± 10.8% across all projection angles and kVs. The differences of ∼121 HU between images with and without scatter, signifying the importance of scatter for simulating clinical images. The denoising methods preserved the magnitudes and trends observed in the reference scatter distributions, with an averaged rRMSE value of 0.91 and 0.97 for the two methods, respectively. The execution time of ∼30 s for simulating scatter in a single projection with a desirable level of statistical noise indicates a major improvement in performance, making our tool an eligible candidate for conducting extensive VITs spanning multiple patients and scan protocols.

XDose: toward online cross-validation of experimental and computational X-ray dose estimation

PURPOSE

As the spectrum of X-ray procedures has increased both for diagnostic and for interventional cases, more attention is paid to X-ray dose management. While the medical benefit to the patient outweighs the risk of radiation injuries in almost all cases, reproducible studies on organ dose values help to plan preventive measures helping both patient as well as staff. Dose studies are either carried out retrospectively, experimentally using anthropomorphic phantoms, or computationally. When performed experimentally, it is helpful to combine them with simulations validating the measurements. In this paper, we show how such a dose simulation method, carried out together with actual X-ray experiments, can be realized to obtain reliable organ dose values efficiently.

METHODS

A Monte Carlo simulation technique was developed combining down-sampling and super-resolution techniques for accelerated processing accompanying X-ray dose measurements. The target volume is down-sampled using the statistical mode first. The estimated dose distribution is then up-sampled using guided filtering and the high-resolution target volume as guidance image. Second, we present a comparison of dose estimates calculated with our Monte Carlo code experimentally obtained values for an anthropomorphic phantom using metal oxide semiconductor field effect transistor dosimeters.

RESULTS

We reconstructed high-resolution dose distributions from coarse ones (down-sampling factor 2 to 16) with error rates ranging from 1.62 % to 4.91 %. Using down-sampled target volumes further reduced the computation time by 30 % to 60 %. Comparison of measured results to simulated dose values demonstrated high agreement with an average percentage error of under [Formula: see text] for all measurement points.

CONCLUSIONS

Our results indicate that Monte Carlo methods can be accelerated hardware-independently and still yield reliable results. This facilitates empirical dose studies that make use of online Monte Carlo simulations to easily cross-validate dose estimates on-site.

Deterministic linear 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 computed tomography (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 s. 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 dose 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.