Purpose

Monte Carlo (MC) simulation methods have been widely used in patient dosimetry in computed tomography (CT), including estimating patient organ doses. However, most simulation methods have undergone a limited set of validations, often using homogeneous phantoms with simple geometries. As clinical scanning has become more complex and the use of tube current modulation (TCM) has become pervasive in the clinic, MC simulations should include these techniques in their methodologies and therefore should also be validated using a variety of phantoms with different shapes and material compositions to result in a variety of differently modulated tube current profiles. The purpose of this work is to perform the measurements and simulations to validate a Monte Carlo model under a variety of test conditions where fixed tube current (FTC) and TCM were used.

Methods

A previously developed MC model for estimating dose from CT scans that models TCM, built using the platform of mcnpx, was used for CT dose quantification. In order to validate the suitability of this model to accurately simulate patient dose from FTC and TCM CT scan, measurements and simulations were compared over a wide range of conditions. Phantoms used for testing range from simple geometries with homogeneous composition (16 and 32 cm computed tomography dose index phantoms) to more complex phantoms including a rectangular homogeneous water equivalent phantom, an elliptical shaped phantom with three sections (where each section was a homogeneous, but different material), and a heterogeneous, complex geometry anthropomorphic phantom. Each phantom requires varying levels of x-, y- and z-modulation. Each phantom was scanned on a multidetector row CT (Sensation 64) scanner under the conditions of both FTC and TCM. Dose measurements were made at various surface and depth positions within each phantom. Simulations using each phantom were performed for FTC, detailed x-y-z TCM, and z-axis-only TCM to obtain dose estimates. This allowed direct comparisons between measured and simulated dose values under each condition of phantom, location, and scan to be made.

Results

For FTC scans, the percent root mean square (RMS) difference between measurements and simulations was within 5% across all phantoms. For TCM scans, the percent RMS of the difference between measured and simulated values when using detailed TCM and z-axis-only TCM simulations was 4.5% and 13.2%, respectively. For the anthropomorphic phantom, the difference between TCM measurements and detailed TCM and z-axis-only TCM simulations was 1.2% and 8.9%, respectively. For FTC measurements and simulations, the percent RMS of the difference was 5.0%.

Conclusions

This work demonstrated that the Monte Carlo model developed provided good agreement between measured and simulated values under both simple and complex geometries including an anthropomorphic phantom. This work also showed the increased dose differences for z-axis-only TCM simulations, where considerable modulation in the x-y plane was present due to the shape of the rectangular water phantom. Results from this investigation highlight details that need to be included in Monte Carlo simulations of TCM CT scans in order to yield accurate, clinically viable assessments of patient dosimetry.

Objective

The U.S. Centers for Medicare & Medicaid Services (CMS) recently approved the use of low-dose CT for lung cancer screening and described volumetric CT dose index (CTDI_vol) requirements. These were based on the National Lung Screening Trial, which used only fixed-tube-current techniques. The aim of this study was to evaluate dose index data from a lung cancer screening program using automatic exposure control (AEC) techniques to ensure compliance with requirements and to correlate dose index values with patient size.

Materials and methods

CTDI_vol, dose-length product (DLP), and body mass index (BMI) data were collected for 563 lung cancer screening examinations performed with AEC between January 1, 2014, through August 31, 2015. CTDI_vol and DLP were analyzed according to the patient's BMI classification. Results were compared with the CMS requirement that the CTDI_vol for a standard-sized patient (height, 170 cm; weight, 70 kg) be 3.0 mGy or less, with adjustments for patients of different sizes. For a subset of patients, the average water-equivalent diameter and size-specific dose estimate were estimated.

Results

The average CTDI_vol for a standard-sized patient was 1.8 mGy, which meets CMS requirements. CTDI_vol values were lower for smaller patients and higher for larger patients. Overall, the mean CTDI_vol and DLP were 2.1 mGy and 74 mGy⋅cm, respectively. The size-specific dose estimate for the average water-equivalent diameter (27.5 cm) of the patient subset was 2.6 mGy.

Conclusion

The screening protocols using AEC resulted in CTDI_vol values that were compliant with CMS requirements. CTDI_vol values greater than 3.0 mGy were only observed for overweight or obese patients.

Purpose

AAPM Task Group 204 introduced size-specific dose estimates for pediatric and adult patients undergoing body CT examinations. This investigation extends that work to head CT exams by using Monte Carlo simulations to develop size-specific, scanner-independent CTDIvol-to-organ-dose conversion coefficients.

Methods

Using eight patient models from the GSF family of voxelized phantoms, dose to the brain and lens of the eye was estimated using Monte Carlo simulations of contiguous axial and helical scans for 64-slice multidetector CT scanners from four major manufacturers. For each patient model and scan mode, scanner-independent CTDIvol-to-organ-dose conversion coefficients were calculated by normalizing organ dose by scanner-specific 16 cm CTDIvol values and averaging across all scanners. Head size was measured using both geometric and attenuation-based size metrics. Head perimeter and effective diameter (ED), both geometric size metrics, were measured directly from the GSF data at the first slice superior to the eyes. Because the GSF models' pixel data are provided in terms of organ identification numbers instead of CT numbers, an indirect estimate of water equivalent diameter (WED), an attenuation-based size metric, was determined based on the relationships between WED and both ED and perimeter for a sample of patient data. Correlations between CTDIvol-to-organ-dose conversion coefficients and the various patient size metrics were then explored.

Results

The analysis of the patient data revealed a best-fit linear relationship (R(2) of 0.87) between ED and WED across a wide variety of patient sizes. Using this relationship along with ED determined from the GSF data, WED was estimated for each GSF model. An exponential relationship between CTDIvol normalized organ dose and WED was observed for both contiguous axial and helical scanning. For head perimeter and ED measured directly from the GSF data, an exponential relationship between CTDIvol normalized organ dose and patient size was also observed for each scan mode. For all patient size metrics and scan modes, R(2) of the exponential fits ranged from 0.92 to 0.93 and 0.73 to 0.85 for the brain and lens of the eye, respectively.

Conclusions

For all scan modes, strong correlation exists between CTDIvol normalized brain dose and both geometric and attenuation-based patient size metrics. A slightly lower correlation between CTDIvol normalized organ dose and patient size was observed for the lens of the eye. This may be due to the combination of the eye lens being a small peripheral organ and the presence of surface dose variation in both contiguous axial and helical scanning. Results indicate that robust estimates of patient-specific head CT dose may be provided using the size-specific, scanner-independent CTDIvol-to-organ-dose conversion coefficients described in this work.

Purpose

The purpose of this study was to estimate the radiation dose to the lung and breast as well as the effective dose from tube current modulated (TCM) lung cancer screening (LCS) scans across a range of patient sizes.

Methods

Monte Carlo (MC) methods were used to calculate lung, breast, and effective doses from a low-dose LCS protocol for a 64-slice CT that used TCM. Scanning parameters were from the protocols published by AAPM's Alliance for Quality CT. To determine lung, breast, and effective doses from lung cancer screening, eight GSF/ICRP voxelized phantom models with all radiosensitive organs identified were used to estimate lung, breast, and effective doses. Additionally, to extend the limited size range provided by the GSF/ICRP phantom models, 30 voxelized patient models of thoracic anatomy were generated from LCS patient data. For these patient models, lung and breast were semi-automatically segmented. TCM schemes for each of the GSF/ICRP phantom models were generated using a validated method wherein tissue attenuation and scanner limitations were used to determine the TCM output as a function of table position and source angle. TCM schemes for voxelized patient models were extracted from the raw projection data. The water equivalent diameter, Dw, was used as the patient size descriptor. Dw was estimated for the GSF/ICRP models. For the thoracic patient models, Dw was extracted from the DICOM header of the CT localizer radiograph. MC simulations were performed using the TCM scheme for each model. Absolute organ doses were tallied and effective doses were calculated using ICRP 103 tissue weighting factors for the GSF/ICRP models. Metrics of scanner radiation output were determined based on each model's TCM scheme, including CTDI_vol , dose length product (DLP), and CTDI_{vol, Low Att} , a previously described regional metric of scanner output covering most of the lungs and breast. All lung and breast doses values were normalized by scan-specific CTDI_vol and CTDI_{vol, Low Att} . Effective doses were normalized by scan-specific CTDI_vol and DLP. Absolute and normalized doses were reported as a function of Dw.

Results

Lung doses normalized by CTDI_{vol, Low Att} were modeled as an exponential relationship with respect to Dw with coefficients of determination (R² ) of 0.80. Breast dose normalized by CTDI_{vol, Low Att} was modeled with an exponential relationship to Dw with an R² of 0.23. For all eight GSF/ICRP phantom models, the effective dose using TCM protocols was below 1.6 mSv. Effective doses showed some size dependence but when normalized by DLP demonstrated a constant behavior.

Conclusion

Lung, breast, and effective doses from LCS CT exams with TCM were estimated with respect to patient size. Normalized lung dose can be reasonably estimated with a measure of a patient size such as Dw and regional metric of CTDI_vol covering the thorax such as CTDI_{vol, Low Att} , while normalized breast dose can also be estimated with a regional metric of CTDI_vol but with a larger degree of variability than observed for lung. Effective dose normalized by DLP can be estimated with a constant multiplier.

Purpose

Task Group Report 195 of the American Association of Physicists in Medicine contains reference datasets for the direct comparison of results among different Monte Carlo (MC) simulation tools for various aspects of imaging research that employs ionizing radiation. While useful for comparing and validating MC codes, that effort did not provide the information needed to compare absolute dose estimates from CT exams. Therefore, the purpose of this work is to extend those efforts by providing a reference dataset for benchmarking fetal dose derived from MC simulations of clinical CT exams.

Acquisition and validation methods

The reference dataset contains the four necessary elements for validating MC engines for CT dosimetry: (a) physical characteristics of the CT scanner, (b) patient information, (c) exam specifications, and (d) fetal dose results from previously validated and published MC simulations methods in tabular form. Scanner characteristics include non-proprietary descriptions of equivalent source cumulative distribution function (CDF) spectra and bowtie filtration profiles, as well as scanner geometry information. Additionally, for the MCNPX MC engine, normalization factors are provided to convert raw simulation results to absolute dose in mGy. The patient information is based on a set of publicly available fetal dose models and includes de-identified image data; voxelized MC input files with fetus, uterus, and gestational sac identified; and patient size metrics in the form of water equivalent diameter (D_w ) z-axis distributions from a simulated topogram (D_w,topo ) and from the image data (D_w,image ). Exam characteristics include CT scan start and stop angles and table and patient locations, helical pitch, nominal collimation and measured beam width, and gantry rotation time for each simulation. For simulations involving estimating doses from exams using tube current modulation (TCM), a realistic TCM scheme is presented that is estimated based upon a validated method. (d) Absolute and CTDI_vol -normalized fetal dose results for both TCM and FTC simulations are given for each patient model under each scan scenario.

Data format and usage notes

Equivalent source CDFs and bowtie filtration profiles are available in text files. Image data are available in DICOM format. Voxelized models are represented by a header followed by a list of integers in a text file representing a three-dimensional model of the patient. Size distribution metrics are also given in text files. Results of absolute and normalized fetal dose with associated MC error estimates are presented in tabular form in an Excel spreadsheet. All data are stored on Zenodo and are publicly accessible using the following link: https://zenodo.org/record/3959512.

Potential applications

Similar to the work of AAPM Report 195, this work provides a set of reference data for benchmarking fetal dose estimates from clinical CT exams. This provides researchers with an opportunity to compare MC simulation results to a set of published reference data as part of their efforts to validate absolute and normalized fetal dose estimates. This could also be used as a basis for comparison to other non-MC approaches, such as deterministic approaches, or to commercial packages that provide estimates of fetal doses from clinical CT exams.

PURPOSE: The vast majority of body CT exams are performed with automatic exposure control (AEC), which adapts the mean tube current to the patient size and modulates the tube current either angularly, longitudinally or both. However, most radiation dose estimation tools are based on fixed tube current scans. Accurate estimates of patient dose from AEC scans require knowledge of the tube current values, which is usually unavailable. The purpose of this work was to develop and validate methods to accurately estimate the tube current values prescribed by one manufacturers AEC system to enable accurate estimates of patient dose. METHODS: Methods were developed that took into account available patient attenuation information, user selected image quality reference parameters and x-ray system limits to estimate tube current values for patient scans. Methods consistent with AAPM Report 220 were developed that used patient attenuation data that were: (a) supplied by the manufacturer in the CT localizer radiograph and (b) based on a simulated CT localizer radiograph derived from image data. For comparison, actual tube current values were extracted from the projection data of each patient. Validation of each approach was based on data collected from 40 pediatric and adult patients who received clinically indicated chest (n = 20) and abdomen/pelvis (n = 20) scans on a 64 slice multidetector row CT (Sensation 64, Siemens Healthcare, Forchheim, Germany). For each patient dataset, the following were collected with Institutional Review Board (IRB) approval: (a) projection data containing actual tube current values at each projection view, (b) CT localizer radiograph (topogram) and (c) reconstructed image data. Tube current values were estimated based on the actual topogram (actual-topo) as well as the simulated topogram based on image data (sim-topo). Each of these was compared to the actual tube current values from the patient scan. In addition, to assess the accuracy of each method in estimating patient organ doses, Monte Carlo simulations were performed by creating voxelized models of each patient, identifying key organs and incorporating tube current values into the simulations to estimate dose to the lungs and breasts (females only) for chest scans and the liver, kidney, and spleen for abdomen/pelvis scans. Organ doses from simulations using the actual tube current values were compared to those using each of the estimated tube current values (actual-topo and sim-topo). RESULTS: When compared to the actual tube current values, the average error for tube current values estimated from the actual topogram (actual-topo) and simulated topogram (sim-topo) was 3.9% and 5.8% respectively. For Monte Carlo simulations of chest CT exams using the actual tube current values and estimated tube current values (based on the actual-topo and sim-topo methods), the average differences for lung and breast doses ranged from 3.4% to 6.6%. For abdomen/pelvis exams, the average differences for liver, kidney, and spleen doses ranged from 4.2% to 5.3%. CONCLUSIONS: Strong agreement between organ doses estimated using actual and estimated tube current values provides validation of both methods for estimating tube current values based on data provided in the topogram or simulated from image data.

Purpose

The purpose of this work was to estimate scanner-independent CTDI_vol -to-fetal-dose coefficients for tube current-modulated (TCM) and fixed tube current (FTC) computed tomography (CT) examinations of pregnant patients of various gestational ages undergoing abdominal/pelvic CT examinations.

Methods

For 24 pregnant patients of gestational age from <5 to 36 weeks who underwent clinically indicated CT examinations, voxelized models of maternal and fetal (or embryo) anatomy were created from abdominal/pelvic image data. Absolute fetal dose (D_fetus ) was estimated using Monte Carlo (MC) simulations of helical scans covering the abdomen and pelvis for TCM and FTC scans. Estimated TCM schemes were generated for each patient model using a validated method that accounts for patient attenuation and scanner output limits for one scanner model and were incorporated into MC simulations. FTC scans were also simulated for each patient model with multidetector row CT scanners from four manufacturers. Normalized fetal dose estimates, nD_fetus , was obtained by dividing D_fetus from the MC simulations by CTDI_vol . Patient size was described using water equivalent diameter (D_w ) measured at the three-dimensional geometric centroid of the fetus. Fetal depth (DE_f ) was measured from the anterior skin surface to the anterior part of the fetus. nD_fetus and D_w were correlated using an exponential model to develop equations for fetal dose conversion coefficients for TCM and FTC abdominal/pelvic CT examinations. Additionally, bivariate linear regression was performed to analyze the correlation of nD_fetus with D_w and fetal depth (DE_f ). For one scanner model, nD_fetus from TCM was compared to FTC and the size-specific dose estimate (SSDE) conversion coefficients (f-factors) from American Association of Physicists in Medicine (AAPM) Report 204. nD_fetus from FTC simulations was averaged across all scanners for each patient (nDfetus¯) . nDfetus¯ was then compared with SSDE f-factors and correlated with D_w using an exponential model and with D_w and DE_f using a bivariate linear model.

Results

For TCM, the coefficient of determination (R² ) of nD_fetus and D_w was observed to be 0.73 using an exponential model. Using the bivariate linear model with D_w and DE_f , an R² of 0.78 was observed. For the TCM technology modeled, TCM yielded nD_fetus values that were on average 6% and 17% higher relative to FTC and SSDE f-factors, respectively. For FTC, the R² of nDfetus¯ with respect to D_w was observed to be 0.64 using an exponential model. Using the bivariate linear model, an R² of 0.75 was observed for nDfetus¯ with respect to D_w and DE_f . A mean difference of 0.4% was observed between nDfetus¯ and SSDE f-factors.

Conclusion

Good correlations were observed for nD_fetus from TCM and FTC scans using either an exponential model with D_w or a bivariate linear model with both D_w and DE_f . These results indicate that fetal dose from abdomen/pelvis CT examinations of pregnant patients of various gestational ages may be reasonably estimated with models that include (a) scanner-reported CTDI_vol and (b) D_w as a patient size metric, in addition to (c) DE_f if available. These results also suggest that SSDE f-factors may provide a reasonable (within ±25%) estimate of nD_fetus for TCM and FTC abdomen/pelvis CT exams.

Although there has been extensive research done on the biological response to doses of ionizing radiation relevant to radiodiagnostic procedures, very few studies have examined radiation schemes similar to those frequently utilized in CT exams. Instead of a single exposure, CT exams are often made up of a series of scans separated on the order of minutes. DNA damage dose-response kinetics after radiation doses and schemes similar to CT protocols were established in both cultured (ESW-WT3) and whole blood lymphocytes and compared to higher dose exposures. Both the kinetics and extent of H2AX phosphorylation were found to be dose dependent. Damage induction and detection showed a clear dose response, albeit different, at all time points and differences in the DNA repair kinetics of ESW-WT3 and whole blood lymphocytes were characterized. Moreover, using a modified split-dose in vitro experiment, we show that phosphorylation of H2AX is significantly reduced after exposure to CT doses fractionated over a few minutes compared to the same total dose delivered as a single exposure. Because the split-dose exposures investigated here are more similar to those experienced during a CT examination, it is essential to understand why and how these differences occur. This work provides compelling evidence supporting differential biological responses not only between high and low doses, but also between single and multiple exposures to low doses of ionizing radiation.

Purpose

Task Group 204 introduced effective diameter (ED) as the patient size metric used to correlate size-specific-dose-estimates. However, this size metric fails to account for patient attenuation properties and has been suggested to be replaced by an attenuation-based size metric, water equivalent diameter (DW). The purpose of this study is to investigate different size metrics, effective diameter, and water equivalent diameter, in combination with regional descriptions of scanner output to establish the most appropriate size metric to be used as a predictor for organ dose in tube current modulated CT exams.

Methods

101 thoracic and 82 abdomen/pelvis scans from clinically indicated CT exams were collected retrospectively from a multidetector row CT (Sensation 64, Siemens Healthcare) with Institutional Review Board approval to generate voxelized patient models. Fully irradiated organs (lung and breasts in thoracic scans and liver, kidneys, and spleen in abdominal scans) were segmented and used as tally regions in Monte Carlo simulations for reporting organ dose. Along with image data, raw projection data were collected to obtain tube current information for simulating tube current modulation scans using Monte Carlo methods. Additionally, previously described patient size metrics [ED, DW, and approximated water equivalent diameter (DWa)] were calculated for each patient and reported in three different ways: a single value averaged over the entire scan, a single value averaged over the region of interest, and a single value from a location in the middle of the scan volume. Organ doses were normalized by an appropriate mAs weighted CTDIvol to reflect regional variation of tube current. Linear regression analysis was used to evaluate the correlations between normalized organ doses and each size metric.

Results

For the abdominal organs, the correlations between normalized organ dose and size metric were overall slightly higher for all three differently (global, regional, and middle slice) reported DW and DWa than they were for ED, but the differences were not statistically significant. However, for lung dose, computed correlations using water equivalent diameter calculated in the middle of the image data (DW,middle) and averaged over the low attenuating region of lung (DW,regional) were statistically significantly higher than correlations of normalized lung dose with ED.

Conclusions

To conclude, effective diameter and water equivalent diameter are very similar in abdominal regions; however, their difference becomes noticeable in lungs. Water equivalent diameter, specifically reported as a regional average and middle of scan volume, was shown to be better predictors of lung dose. Therefore, an attenuation-based size metric (water equivalent diameter) is recommended because it is more robust across different anatomic regions. Additionally, it was observed that the regional size metric reported as a single value averaged over a region of interest and the size metric calculated from a single slice/image chosen from the middle of the scan volume are highly correlated for these specific patient models and scan types.

Purpose

The purpose of this study was to assess the accuracy of a Monte Carlo simulation-based method for estimating radiation dose from multidetector computed tomography (MDCT) by comparing simulated doses in ten patients to in-vivo dose measurements.

Methods

MD Anderson Cancer Center Institutional Review Board approved the acquisition of in-vivo rectal dose measurements in a pilot study of ten patients undergoing virtual colonoscopy. The dose measurements were obtained by affixing TLD capsules to the inner lumen of rectal catheters. Voxelized patient models were generated from the MDCT images of the ten patients, and the dose to the TLD for all exposures was estimated using Monte Carlo based simulations. The Monte Carlo simulation results were compared to the in-vivo dose measurements to determine accuracy.

Results

The calculated mean percent difference between TLD measurements and Monte Carlo simulations was -4.9% with standard deviation of 8.7% and a range of -22.7% to 5.7%.

Conclusions

The results of this study demonstrate very good agreement between simulated and measured doses in-vivo. Taken together with previous validation efforts, this work demonstrates that the Monte Carlo simulation methods can provide accurate estimates of radiation dose in patients undergoing CT examinations.