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

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.

Fast prediction of personalized abdominal organ doses from CT examinations by radiomics feature-based machine learning models

The X-rays emitted during CT scans can increase solid cancer risks by damaging DNA, with the risk tied to patient-specific organ doses. This study aims to establish a new method to predict patient specific abdominal organ doses from CT examinations using minimized computational resources at a fast speed. The CT data of 247 abdominal patients were selected and exported to the auto-segmentation software named DeepViewer to generate abdominal regions of interest (ROIs). Radiomics feature were extracted based on the selected CT data and ROIs. Reference organ doses were obtained by GPU-based Monte Carlo simulations. The support vector regression (SVR) model was trained based on the radiomics features and reference organ doses to predict abdominal organ doses from CT examinations. The prediction performance of the SVR model was tested and verified by changing the abdominal patients of the train and test sets randomly. For the abdominal organs, the maximal difference between the reference and the predicted dose was less than 1 mGy. For the body and bowel, the organ doses were predicted with a percentage error of less than 5.2%, and the coefficient of determination (R2) reached up to 0.9. For the left kidney, right kidney, liver, and spinal cord, the mean absolute percentage error ranged from 5.1 to 8.9%, and the R2 values were more than 0.74. The SVR model could be trained to achieve accurate prediction of personalized abdominal organ doses in less than one second using a single CPU core.

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.

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.

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.