Estimating radiation doses from multidetector CT using Monte Carlo simulations: effects of different size voxelized patient models on magnitudes of organ and effective dose

Automatic organ completion with image stitching for personalized radiation dosimetry in CT examinations

PURPOSE

Computed tomography (CT) image-based patient-specific voxel-based dosimetry has difficulties complementing missing tissues for organs located partially inside or completely outside the image volume. Previous studies constructed patient-specific whole-body models by rescaling reference phantoms or extending regional CT images with manually adjusted phantoms. This study proposes a methodology for automatic organ completion of regional CT images for CT dosimetry using a stitching approach.

METHODS

Virtual clinical trials were performed by truncating whole-body CT images to generate virtual clinical chest and abdominopelvic CT images. Corresponding anchor images for each patient were selected according to sex and similarity of the axial length and water equivalent diameter of the virtual regional CT images. Automatic image stitching was performed by transformation initialization and iteration, while the stitched CT images and organ atlas were used in GPU-based Geant4 Monte Carlo simulations to generate a radiation dose map and absorbed organ dose. To evaluate the performance of the stitching model in radiation dosimetry, organ mass differences and Jaccard's coefficient of stitched and rescaled anchor images were calculated, and the radiation doses were compared among the corresponding values from the VirtualDose®, original whole-body CT, stitching model, regional CT, registration-based rescaling method, and WED-based rescaling method.

RESULTS

The anatomical accuracy of stitched images was significantly improved. For organs partially inside the image volume, organ dose estimation from the stitching model could be more accurate than that reported in previous studies. The absolute differences in effective dose from the stitched images were 6.55% and 4.81% for chest and abdominopelvic CT scans, respectively.

CONCLUSION

The proposed automatic stitching model partially complements organs inside or outside the CT scan range and provides more accurate anatomical representations for radiation dosimetry than traditional phantom rescaling methods.

CT scanner-specific organ dose coefficients generated by Monte Carlo calculation for the ICRP adult male and female reference computational phantoms

Objective.Provide analyses of new organ dose coefficients (hereafter also referred to as normalized doses) for CT that have been developed to update the widely-utilized collection of data published 30 years ago in NRPB-SR250.Approach.In order to reflect changes in technology, and also ICRP recommendations concerning use of the computational phantoms adult male (AM) and adult female (AF), 102 series of new Monte Carlo simulations have been performed covering the range of operating conditions for 12 contemporary models of CT scanner from 4 manufacturers. Normalized doses (relative to free air on axis) have been determined for 39 organs, and for every 8 mm or 4.84 mm slab of AM and AF, respectively.Main results.Analyses of results confirm the significant influence (by up to a few tens of percent), on values of normalized organ (or contributions to effective dose (E_103,phan)), for whole body exposure arising from selection of tube voltage and beam shaping filter. Use of partial (when available) rather than a Full fan beam reduced both organ and effective dose by up to 7%. Normalized doses to AF were larger than corresponding figures for AM by up to 30% for organs and by 10% forE_103,phan. Additional simulations for whole body exposure have also demonstrated that: practical simplifications in the main modelling (point source, single slice thickness, neglect of patient couch and immobility of phantom arms) have sufficiently small (<5%) effect onE_103,phan; mis-centring of the phantom away from the axis of rotation by 5 mm (in any direction) leads to changes in normalized organ dose andE_103,phanby up to 20% and 6%, respectively; and angular tube current modulation can result in reductions by up to 35% and <15% in normalized organ dose andE_103,phan, respectively, for 100% cosine variation.Significance.These analyses help advance understanding of the influence of operational scanner settings on organ dose coefficients for contemporary CT, in support of improved patient protection. The results will allow the future development of a new dose estimation tool.

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.

MEASUREMENT OF HALF-VALUE LAYER IN COMPUTED TOMOGRAPHY SCANNERS USING LUMINESCENCE OF POLYETHERSULFONE RESIN BY X-RAY IRRADIATION

In this study, a method for estimating the half-value layer (HVL) and effective energy (Eeff) by imaging the luminescence from a polyethersulfone (PES) resin with rotating irradiation of X-rays in a computed tomography scanner was developed. The luminescence of the PES resin was imaged using a charge-coupled device camera. The PES-HVL was determined from the luminance attenuation profile corresponding to the X-ray attenuation within the resin. The PES-HVLs for tube potentials of 80-135 kVp were converted into Eeff values and were compared to those of a conventional lead-covered case method. The Eeff obtained using the proposed luminescence imaging method agreed within ~3.9% of that obtained using the conventional method. Moreover, dose simulations based on the X-ray spectrum calculated from the HVLs were performed using a poly(methyl methacrylate) phantom with a diameter of 16 cm. The simulated doses based on the luminescence imaging method agreed with the in-phantom dosimetry within ~9%.

A NOVEL METHOD FOR ESTIMATING PATIENT-SPECIFIC PRIMARY DOSE IN CONE-BEAM COMPUTED TOMOGRAPHY

For the purpose of real-time scan-protocol optimisation and patient-specific dose management in cone-beam computed tomography, we introduce a numerical algorithm that estimates the primary dose distributions in reconstructed images. The proposed algorithm is based on the ray-tracing technique and utilises reconstructed voxel data and scanning protocol. The algorithm is validated with the Monte Carlo (MC) and conventional model-based dose reconstruction methods for the simple cylindrical water and anthropomorphic head phantoms. The algorithm shows good agreement with both methods in terms of the zeroth-order x-ray interactions, which exclude the higher-order x-ray interactions at sites distant from the first interactions, and it consumes a significantly lower computational cost compared with the MC method. The differences between the proposed algorithm and the model-based dose reconstruction method as well as the improvement strategies of the algorithm are discussed in detail.

Investigation of Eye Lens Dose Estimate based on AAPM Report 293 in Head Computed Tomography

Background:

Estimation of eye lens dose is important in head computed tomography (CT) examination since the eye lens is a sensitive organ to ionizing radiation.

Objective:

The purpose of this study is to compare estimations of eye lens dose in head CT examinations using local size-specific dose estimate (SSDE) based on size-conversion factors of the American Association of Physicists in Medicine (AAPM) Report No. 293 with those based on size-conversion factors of the AAPM Report No. 220.

Material and Methods:

This experimental study is conducted on a group of patients who had undergone nasopharyngeal CT examination. Due to the longitudinal (z-axis) dose fluctuation, the average global SSDE and average local SSDE (i.e. particular slices where the eyes are located) were investigated. All estimates were compared to the measurement results using thermo-luminescent dosimeters (TLDs). The estimated and measured doses were implemented for 14 patients undergoing nasopharyngeal CT examination.

Results:

It was found that the percentage differences of the volume CT dose index (CTDI_vol), average global SSDE based on AAPM No. 220 (SSDE_o,g), average local SSDE based on AAPM No. 220 (SSDE_o,l), average global SSDE based on AAPM No. 293 (SSDE_n,g) and average local SSDE based on AAPM No. 293 (SSDE_n,l) against the measured TLD doses were 22.5, 21.7, 15.0, 9.3, and 2.1%, respectively. All comparisons between dose estimates and TLD measurements gave p-values less than 0.001, except for SSDE_n,l (p-value = 0.566).

Conclusion:

SSDE based on AAPM Report No. 293 can be used to accurately estimate eye lens radiation doses by performing the calculations on a number of specific slices containing the eyes.

Voxel model of a rabbit: assessment of absorbed doses in organs after CT examination performed by two different protocols

The objective of this work was to assess absorbed doses in organs and tissues of a rabbit, following computed tomography (CT) examinations, using a dedicated 3D voxel model. Absorbed doses in relevant organs were calculated using the MCNP5 Monte Carlo software. Calculations were perfomed for two standard CT protocols, using tube voltages of 110 kVp and 130 kVp. Absorbed doses were calculated in 11 organs and tissues, i.e., skin, bones, brain, muscles, heart, lungs, liver, spleen, kidney, testicles, and fat tissue. The doses ranged from 15.3 to 28.3 mGy, and from 40.2 to 74.3 mGy, in the two investigated protocols. The organs that received the highest dose were bones and kidneys. In contrast, brain and spleen were organs that received the smallest doses. Doses in organs which are stretched along the body did not change significantly with distance. On the other hand, doses in organs which are localized in the body showed maximums and minimums. Using the voxel model, it is possible to calculate the dose distribution in the rabbit's body after CT scans, and study the potential biological effects of CT doses in certain organs. The voxel model presented in this work can be used to calculated doses in all radiation experiments in which rabbits are used as experimental animals.

A comparison of breast and lung doses from chest CT scans using organ-based tube current modulation (OBTCM) vs. Automatic tube current modulation (ATCM)

PURPOSE

The purpose of this work was to estimate and compare breast and lung doses of chest CT scans using organ-based tube current modulation (OBTCM) to those from conventional, attenuation-based automatic tube current modulation (ATCM) across a range of patient sizes.

METHODS

Thirty-four patients (17 females, 17 males) who underwent clinically indicated CT chest/abdomen/pelvis (CAP) examinations employing OBTCM were collected from two multi-detector row CT scanners. Patient size metric was assessed as water equivalent diameter (Dw ) taken at the center of the scan volume. Breast and lung tissues were segmented from patient image data to create voxelized models for use in a Monte Carlo transport code. The OBTCM schemes for the chest portion were extracted from the raw projection data. ATCM schemes were estimated using a recently developed method. Breast and lung doses for each TCM scenario were estimated for each patient model. CTDIvol -normalized breast (nDbreast ) and lung (nDlung ) doses were subsequently calculated. The differences between OBTCM and ATCM normalized organ dose estimates were tested using linear regression models that included CT scanner and Dw as covariates.

RESULTS

Mean dose reduction from OBTCM in nDbreast was significant after adjusting for the scanner models and patient size (P = 0.047). When pooled with females and male patient, mean dose reduction from OBTCM in nDlung was observed to be trending after adjusting for the scanner model and patient size (P = 0.085).

CONCLUSIONS

One specific manufacturer's OBTCM was analyzed. OBTCM was observed to significantly decrease normalized breast relative to a modeled version of that same manufacturer's ATCM scheme. However, significant dose savings were not observed in lung dose over all. Results from this study support the use of OBTCM chest protocols for females only.

Reference dataset for benchmarking fetal doses derived from Monte Carlo simulations of CT exams

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.

Photon fluence and dose estimation in computed tomography using a discrete ordinates Boltzmann solver

In this study, cone-beam single projection and axial CT scans are modeled with a software package—DOCTORS, which solves the linear Boltzmann equation using the discrete ordinates method. Phantoms include a uniform 35 cm diameter water cylinder and a non-uniform abdomen phantom. Series simulations were performed with different simulation parameters, including the number of quadrature angles, the order of Legendre polynomial expansions, and coarse and fine mesh grid. Monte Carlo simulations were also performed to benchmark DOCTORS simulations. A quantitative comparison was made between the simulation results obtained using DOCTORS and Monte Carlo methods. The deterministic simulation was in good agreement with the Monte Carlo simulation on dose estimation, with a root-mean-square-deviation difference of around 2.87%. It was found that the contribution of uncollided photon fluence directly from the source dominates the local absorbed dose in the diagnostic X-ray energy range. The uncollided photon fluence can be calculated accurately using a ‘ray-tracing’ algorithm. The accuracy of collided photon fluence estimation is largely affected by the pre-calculated multigroup cross-sections. The primary benefit of DOCTORS lies in its rapid computation speed. Using DOCTORS, parallel computing with GPU enables the cone-beam CT dose estimation nearly in real-time.

Comparison of central, peripheral, and weighted size-specific dose in CT

The objective of this study is to determine X-ray dose distribution and the correlation between central, peripheral and weighted-centre peripheral doses for various phantom sizes and tube voltages in computed tomography (CT). We used phantoms developed in-house, with various water-equivalent diameters (Dw) from 8.5 up to 42.1 cm. The phantoms have one hole in the centre and four holes at the periphery. By using these five holes, it is possible to measure the size-specific central dose (Ds,c), peripheral dose (Ds,p), and weighted dose (Ds,w).The phantoms are scanned using a CT scanner (Siemens Somatom Definition AS), with the tube voltage varied from 80 up to 140 kVps. The doses are measured using a pencil ionization chamber (Ray safe X2 CT Sensor) in every hole for all phantoms. The relationships between Ds,c, Ds,p, and Ds,w, and the water-equivalent diameter are established. The size-conversion factors are calculated. Comparisons between Ds,c, Ds,p, and Ds,ware also established. We observe that the dose is relatively homogeneous over the phantom for water-equivalent diameters of 12-14 cm. For water-equivalent diameters less than 12 cm, the dose in the centre is higher than at the periphery, whereas for water-equivalent diameters greater than 14 cm, the dose at the centre is lower than that at the periphery. We also find that the distribution of the doses is influenced by the tube voltage. These dose distributions may be useful for calculating organ doses for specific patients using their CT images in future clinical practice.

SHAPE ESTIMATION OF BOWTIE FILTERS BASED ON THE LUMINESCENCE FROM POLYETHYLENE TEREPHTHALATE RESIN BY X-RAY IRRADIATION

In this study, we devised a novel method estimating the bowtie filter shapes by imaging luminescence from a polyethylene terephthalate (PET) resin with X-ray irradiation in a computed tomography (CT) scanner. The luminescence distribution of the PET resin corresponding to the thickness of bowtie filter was imaged using a charge-coupled device camera. On the assumption that the material of bowtie filter is aluminium (Al), the shape of bowtie filters was estimated from the correlation between Al attenuation curves and the angular-dependent luminance attenuation profiles according to the thickness of bowtie filters. Dose simulations based on the estimated bowtie filter shapes were performed using head and body PMMA phantoms with 16 and 32 cm in diameter. The simulated values of head and body weighted CT dose index (CTDIw) based on bowtie filter shape by the luminescence imaging method agreed within ~9% with the measured values by a dosemeter.

PATIENTS DOSE ESTIMATION IN CT EXAMINATIONS USING SIZE SPECIFIC DOSE ESTIMATES

The present work reports data of radiation exposure to the patients during head, chest, pelvis and abdomen CT examinations performed on a third-generation 16-slice CT machine. Radiation exposure was estimated using size specific dose estimates (SSDE) method, which takes into account patient's physical dimensions in phantom measured computed tomography dose index (CTDI) value. The reported median CT dose volume index CTDIvol values in head, chest, pelvis and abdomen examinations were 26.76, 16.27, 29.81 and 14.74 mGy, respectively. The median doses evaluated using SSDE methodology for the above mentioned procedure were 54.1, 23.1, 42.8 and 20.1 mGy, respectively. Our results showed variation in dose values estimated using CTDI and SSDE methods in all examinations. The evaluated SSDE values were also compared to the values derived from data reported by the American Association of Physicist in Medicine (AAPM). SSDE values in present measurements are 4-8% lower than AAPM values. The present results show that CTDI parameters recorded on CT console should not be used to specify patient dose during CT procedures.

ESTIMATION OF EFFECTIVE DOSE OF DENTAL X-RAY DEVICES

This study aims to estimate the effective doses of dental X-ray devices under common scanning protocols. After putting TLDs in the Alderson Radiation Therapy Phantom, we exposed the phantom under common scanning protocols of three dental X-ray devices, namely CBCT, dental panoramic machine and intraoral round cone device. Then effective doses were calculated using the measured absorbed doses of organs and tissues. Tissue weighting factors recommended by the ICRP were adopted in the calculation. Effective doses under common scanning protocols of three Dental X-ray devices were obtained. The effective dose of dental CT was 0.20 mSv, and that of dental panoramic machine and intraoral radiography were 0.013 and 0.0050 mSv, respectively. The tissue absorbed doses of dental CT scan were 0.63 mGy of brain, 7.7 mGy of salivary glands, 8.7 mGy of thyroid and 4.0 mGy of the lens of the eye. The tissue absorbed doses from dental panoramic machine are 0.62 mGy of salivary glands and 0.25 mGy of thyroid. And finally the tissue absorbed dose of intraoral radiography was 0.80 mGy of salivary gland. Among the three dental X-ray devices studied, dental CBCT scan can cause much higher effective dose than the other two. Brain, salivary glands, thyroid and the lens of the eye are tissues receiving relatively higher absorbed doses.

Estimating fetal dose from tube current-modulated (TCM) and fixed tube current (FTC) abdominal/pelvis CT examinations

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 ( n D fetus ¯ ) . n D fetus ¯ 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 n D fetus ¯ 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 n D fetus ¯ with respect to D_w and DE_f . A mean difference of 0.4% was observed between n D fetus ¯ 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.

Estimation of lifetime attributable risks (LARs) of cancer associated with abdominopelvic radiotherapy treatment planning computed tomography (CT) simulations

PURPOSE

The present study attempts to calculate organ-absorbed and effective doses for cancer patients to estimate the possible cancer induction and cancer mortality risks resulting from 64-slice abdominopelvic computed tomography (CT) simulations for radiotherapy treatment planning (RTTP).

MATERIAL AND METHODS

A group of 70 patients, who underwent 64-slice abdominopelvic CT scan for RTTP, voluntarily participated in the present study. To calculate organ and effective doses in a standard phantom of 70 kg, the collected dosimetric parameters were used with the ImPACT CT Patient Dosimetry Calculator. Patient-specific organ dose and effective dose were calculated by applying related correction factors. For the estimation of lifetime attributable risks (LARs) of cancer incidence and cancer-related mortality, doses in radiosensitive organs were converted to risks based on the data published in Biological Effects of Ionizing Radiation VII (BEIR VII).

RESULTS

The mean ± standard deviation (SD) of the effective dose for males and females were 13.87 ± 2.37 mSv (range: 9.25-18.82 mSv) and 13.04 ± 3.42 mSv (range: 6.99-18.37 mSv), respectively. The mean ± SD of LAR of cancer incidence was 35.34 ± 13.82 cases in males and 34.49 ± 9.63 cases in females per 100,000 persons. The LAR of cancer mortality had the mean ± SD value of 15.38 ± 4.25 and 16.72 ± 3.87 cases per 100,000 persons in males and females respectively.

CONCLUSION

Increase in the LAR of cancer occurrence and mortality due to abdominopelvic treatment planning CT simulation is noticeable and should be considered.

Determination and verification of the x-ray spectrum of a CT scanner

The accuracy of Monte Carlo (MC) simulations in estimating the computed tomography radiation dose is highly dependent on the proprietary x-ray source information. To address this, this study develops a method to precisely estimate the x-ray spectrum and bowtie (BT) filter thickness of the x-ray source based on physical measurements and calculations. The static x-ray source of the CT localizer radiograph was assessed to measure the total filtration at the isocenter for the x-ray spectrum characterization and the BT profile (air-kerma values as a function of fan angle). With these values, the utilized BT filter in the localizer radiograph was assessed by integrating the measured air kerma in a full 360-deg cycle. The consistency observed between the integrated BT filter profiles and the directly measured profiles pointed to the similarity in the utilized BT filter in terms of thickness and material between the static and rotating x-ray geometries. Subsequently, the measured air kerma was used to calculate the BT filter thickness and was verified using MC simulations by comparing the calculated and measured air-kerma values, where a very good agreement was observed. This would allow a more accurate computed tomography simulation and facilitate the estimation of the dose delivered to the patients.

New weighting factor of weighted CTDI equation for PMMA phantom diameter from 8 to 40 cm: A Monte Carlo study

PURPOSE

The weighted computed tomography dose index (CTDI_w ) uses measured CTDI values at the center and periphery of a cylindrical phantom. The CTDI_w value is calculated using conventional, Bakalyar's, and Choi's weighting factors. However, these weighting factors were produced from only 16- and 32-cm-diameter cylindrical phantoms. This study aims to devise new weighting factors to provide more accurate average dose in the central cross-sectional plane of cylindrical phantoms over a wide range of object diameters, by using Monte Carlo simulations.

METHODS

Simulations were performed by modeling a Toshiba Aquilion ONE CT scanner, in order to compute the cross-sectional dose profiles of polymethyl methacrylate (PMMA) cylindrical phantoms of each diameter (8-40 cm at 4-cm steps), for various tube voltages and longitudinal beam widths. Two phantom models were simulated, corresponding to the CTDI₁₀₀ method and the method recommended by American Association of Physicists in Medicine (AAPM) task group 111. The dose-computation PMMA cylinders of 1 mm diameter were located between the phantom surfaces and the centers at intervals of 1 mm, from which cross-sectional dose profiles were calculated. By using linear least-squares fits to the obtained cross-sectional dose profiles data, we determined new weighting factors to estimate more accurate average doses in the PMMA cylindrical phantoms by using the CTDI_w equation: CTDI_w = W_center ・ CTDI_center + W_periphery ・ CTDI_periphery . In order to demonstrate the validity of the devised new weighting factors, the percentage difference between average dose and CTDI_w value was evaluated for the weighting factors (conventional, Bakalyar's, Choi's, and devised new weighting factors) in each calculated cross-sectional dose profile.

RESULTS

With the use of linear least-squares techniques, new weighting factors (W_center = 3/8 and W_periphery = 5/8 where W_center and W_periphery are weighting factors for CTDI_center and CTDI_periphery ) were determined. The maximum percentage differences between average dose and CTDI_w value were 16, -12, -8, and -6% for the conventional, Bakalyar's, Choi's, and devised new weighting factors, respectively.

CONCLUSIONS

We devised new weighting factors (W_center = 3/8 and W_periphery = 5/8) to provide more accurate average dose estimation in PMMA cylindrical phantoms over a wide range of diameter. The CTDI_w equation with devised new weighting factors could estimate average dose in PMMA cylindrical phantoms with a maximum difference of -6%. The results of this study can estimate the average dose in PMMA cylindrical phantoms more accurately than the conventional weighting factors (W_center = 1/3 and W_periphery = 2/3).

Quantification of radiation dose reduction by reducing z-axis coverage in 320-detector coronary CT angiography

OBJECTIVE

To quantify the radiation dose reduction achievable by minimizing z-axis coverage in 320-detector coronary CT angiography (CCTA).

METHODS

We retrospectively reviewed 130 CCTAs performed on 320-detector CT that offers up to 16 cm z-axis coverage (adjustable in 2-cm increments), allowing complete coverage of the heart in a single gantry rotation. For each CT, we obtained the radiation dose [CT dose index and dose-length product (DLP)], measured the z-axis field of view and measured the craniocaudal cardiac size (distance from the left main coronary artery to the cardiac apex). We calculated the radiation dose savings achievable by reducing the z-axis coverage to the minimum necessary to cover the heart using 320 × 0.5-mm (maximum 16 cm) and 256 × 0.5-mm (maximum 12.8 cm) detector collimations.

RESULTS

Results are expressed as mean ± standard deviation. The mean craniocaudal cardiac size was 10.5 ± 1.0 cm, with 85% (n = 112) of CCTAs performed with 16 cm of z-axis coverage. The mean DLP was 417.6 ± 182.4 mGy cm, with the mean DLP saving achievable using the minimum z-axis coverage required to completely image the heart being 96.2 ± 47.4 mGy cm, an average dose reduction of 26.9 ± 7.0%. z-axis coverage of ≤12 cm was adequate for 92% and 12.8 cm for 98% of subjects.

CONCLUSION

Using the minimal z-axis coverage to adequately image the heart is a simple step that can reduce the DLP in 320-detector CCTA by approximately 27%. z-axis coverage of ≤12 cm is adequate for 92%, 12.8 cm for 98% and 14 cm for 100% of patients undergoing CCTA. Advances in knowledge: Reducing z-axis coverage in 320-detector CCTA can reduce DLP by approximately 27%.

Estimating organ doses from tube current modulated CT examinations using a generalized linear model

PURPOSE

Currently, available Computed Tomography dose metrics are mostly based on fixed tube current Monte Carlo (MC) simulations and/or physical measurements such as the size specific dose estimate (SSDE). In addition to not being able to account for Tube Current Modulation (TCM), these dose metrics do not represent actual patient dose. The purpose of this study was to generate and evaluate a dose estimation model based on the Generalized Linear Model (GLM), which extends the ability to estimate organ dose from tube current modulated examinations by incorporating regional descriptors of patient size, scanner output, and other scan-specific variables as needed.

METHODS

The collection of a total of 332 patient CT scans at four different institutions was approved by each institution's IRB and used to generate and test organ dose estimation models. The patient population consisted of pediatric and adult patients and included thoracic and abdomen/pelvis scans. The scans were performed on three different CT scanner systems. Manual segmentation of organs, depending on the examined anatomy, was performed on each patient's image series. In addition to the collected images, detailed TCM data were collected for all patients scanned on Siemens CT scanners, while for all GE and Toshiba patients, data representing z-axis-only TCM, extracted from the DICOM header of the images, were used for TCM simulations. A validated MC dosimetry package was used to perform detailed simulation of CT examinations on all 332 patient models to estimate dose to each segmented organ (lungs, breasts, liver, spleen, and kidneys), denoted as reference organ dose values. Approximately 60% of the data were used to train a dose estimation model, while the remaining 40% was used to evaluate performance. Two different methodologies were explored using GLM to generate a dose estimation model: (a) using the conventional exponential relationship between normalized organ dose and size with regional water equivalent diameter (WED) and regional CTDI_vol as variables and (b) using the same exponential relationship with the addition of categorical variables such as scanner model and organ to provide a more complete estimate of factors that may affect organ dose. Finally, estimates from generated models were compared to those obtained from SSDE and ImPACT.

RESULTS

The Generalized Linear Model yielded organ dose estimates that were significantly closer to the MC reference organ dose values than were organ doses estimated via SSDE or ImPACT. Moreover, the GLM estimates were better than those of SSDE or ImPACT irrespective of whether or not categorical variables were used in the model. While the improvement associated with a categorical variable was substantial in estimating breast dose, the improvement was minor for other organs.

CONCLUSIONS

The GLM approach extends the current CT dose estimation methods by allowing the use of additional variables to more accurately estimate organ dose from TCM scans. Thus, this approach may be able to overcome the limitations of current CT dose metrics to provide more accurate estimates of patient dose, in particular, dose to organs with considerable variability across the population.

Evaluation of Effective Dose from CT Scans for Overweight and Obese Adult Patients Using the VirtualDose Software

This paper evaluates effective dose (ED) of overweight and obese patients who undergo body computed tomography (CT) examinations. ED calculations were based on tissue weight factors in the International Commission on Radiological Protection Publication 103 (ICRP 103). ED per unit dose length product (DLP) are reported as a function of the tube voltage, body mass index (BMI) of patient. The VirtualDose software was used to calculate ED for male and female obese phantoms representing normal weight, overweight, obese 1, obese 2 and obese 3 patients. Five anatomic regions (chest, abdomen, pelvis, abdomen/pelvis and chest/abdomen/pelvis) were investigated for each phantom. The conversion factors were computed from the DLP, and then compared with data previously reported by other groups. It was observed that tube voltage and BMI are the major factors that influence conversion factors of obese patients, and that ED computed using ICRP 103 tissue weight factors were 24% higher for a CT chest examination and 21% lower for a CT pelvis examination than the ED using ICRP 60 factors. For body CT scans, increasing the tube voltage from 80 to 140 kVp would increase the conversion factors by as much as 19-54% depending on the patient's BMI. Conversion factor of female patients was ~7% higher than the factors of male patients. DLP and conversion factors were used to estimate ED, where conversion factors depended on tube voltage, sex, BMI and tissue weight factors. With increasing number of obese individuals, using size-dependence conversion factors will improve accuracy, in estimating patient radiation dose.

The importance of BMI in dosimetry of 153Sm-EDTMP bone pain palliation therapy: A Monte Carlo study

Using digital phantoms as an atlas compared to acquiring CT data for internal radionuclide dosimetry decreases patient overall radiation dose and reduces the required analysis effort and time for organ segmentation. The drawback is that the phantom may not match exactly with the patient. We assessed the effect of varying BMIs on dosimetry results for a bone pain palliation agent, ¹⁵³Sm-EDTMP. The simulation was done using the GATE Monte Carlo code. Female XCAT phantoms with the following different BMIs were employed: 18.6, 20.8, 22.1, 26.8, 30.3 and 34.7kg/m². S-factors (mGy/MBq.s) and SAFs (kg^-1) were calculated for the dosimetry of the radiation from major source organs including spine, ribs, kidney and bladder into different target organs as well as whole body dosimetry from spine. The differences in dose estimates from different phantoms compared to those from the phantom with BMI of 26.8kg/m² as the reference, were calculated for both gamma and beta radiations. The relative differences (RD) of the S-factors or SAFs from the values of reference phantom were calculated. RDs greater than 10% and 100% were frequent in radiations to organs for photon and beta particles, respectively. The relative differences in whole body SAFs from the reference phantom were 15.4%, 7%, 4.2%, -9.8% and -1.4% for BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m², respectively. The differences in whole body S-factors for the phantoms with BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m² were 39.5%, 19.4%, 8.8%, -7.9% and -4.3%, respectively. The dosimetry of the gamma photons and beta particles changes substantially with the use of phantoms with different BMIs. The change in S-factors is important for dose calculation and can change the prescribed therapeutic dose of ¹⁵³Sm-EDTMP. Thus a phantom with BMI better matched to the patient is suggested for therapeutic purposes where dose estimates closer to those in the actual patient are required.

A measurement-based generalized source model for Monte Carlo dose simulations of CT scans

The goal of this study is to develop a generalized source model for accurate Monte Carlo dose simulations of CT scans based solely on the measurement data without a priori knowledge of scanner specifications. The proposed generalized source model consists of an extended circular source located at x-ray target level with its energy spectrum, source distribution and fluence distribution derived from a set of measurement data conveniently available in the clinic. Specifically, the central axis percent depth dose (PDD) curves measured in water and the cone output factors measured in air were used to derive the energy spectrum and the source distribution respectively with a Levenberg-Marquardt algorithm. The in-air film measurement of fan-beam dose profiles at fixed gantry was back-projected to generate the fluence distribution of the source model. A benchmarked Monte Carlo user code was used to simulate the dose distributions in water with the developed source model as beam input. The feasibility and accuracy of the proposed source model was tested on a GE LightSpeed and a Philips Brilliance Big Bore multi-detector CT (MDCT) scanners available in our clinic. In general, the Monte Carlo simulations of the PDDs in water and dose profiles along lateral and longitudinal directions agreed with the measurements within 4%/1 mm for both CT scanners. The absolute dose comparison using two CTDI phantoms (16 cm and 32 cm in diameters) indicated a better than 5% agreement between the Monte Carlo-simulated and the ion chamber-measured doses at a variety of locations for the two scanners. Overall, this study demonstrated that a generalized source model can be constructed based only on a set of measurement data and used for accurate Monte Carlo dose simulations of patients' CT scans, which would facilitate patient-specific CT organ dose estimation and cancer risk management in the diagnostic and therapeutic radiology.

Mathematical modelling of scanner-specific bowtie filters for Monte Carlo CT dosimetry

The purpose of bowtie filters in CT scanners is to homogenize the x-ray intensity measured by the detectors in order to improve the image quality and at the same time to reduce the dose to the patient because of the preferential filtering near the periphery of the fan beam. For CT dosimetry, especially for Monte Carlo calculations of organ and tissue absorbed doses to patients, it is important to take the effect of bowtie filters into account. However, material composition and dimensions of these filters are proprietary. Consequently, a method for bowtie filter simulation independent of access to proprietary data and/or to a specific scanner would be of interest to many researchers involved in CT dosimetry. This study presents such a method based on the weighted computer tomography dose index, CTDI_w, defined in two cylindrical PMMA phantoms of 16 cm and 32 cm diameter. With an EGSnrc-based Monte Carlo (MC) code, ratios CTDI_w/CTDI_100,a were calculated for a specific CT scanner using PMMA bowtie filter models based on sigmoid Boltzmann functions combined with a scanner filter factor (SFF) which is modified during calculations until the calculated MC CTDI_w/CTDI_100,a matches ratios CTDI_w/CTDI_100,a, determined by measurements or found in publications for that specific scanner. Once the scanner-specific value for an SFF has been found, the bowtie filter algorithm can be used in any MC code to perform CT dosimetry for that specific scanner. The bowtie filter model proposed here was validated for CTDI_w/CTDI_100,a considering 11 different CT scanners and for CTDI_100,c, CTDI_100,p and their ratio considering 4 different CT scanners. Additionally, comparisons were made for lateral dose profiles free in air and using computational anthropomorphic phantoms. CTDI_w/CTDI_100,a determined with this new method agreed on average within 0.89% (max. 3.4%) and 1.64% (max. 4.5%) with corresponding data published by CTDosimetry (www.impactscan.org) for the CTDI HEAD and BODY phantoms, respectively. Comparison with results calculated using proprietary data for the PHILIPS Brilliance 64 scanner showed agreement on average within 2.5% (max. 5.8%) and with data measured for that scanner within 2.1% (max. 3.7%). Ratios of CTDI_100,c/CTDI₁₀₀, _p for this study and corresponding data published by CTDosimetry (www.impactscan.org) agree on average within about 11% (max. 28.6%). Lateral dose profiles calculated with the proposed bowtie filter and with proprietary data agreed within 2% (max. 5.9%), and both calculated data agreed within 5.4% (max. 11.2%) with measured results. Application of the proposed bowtie filter and of the exactly modelled filter to human phantom Monte Carlo calculations show agreement on the average within less than 5% (max. 7.9%) for organ and tissue absorbed doses.

The use of ultra-low-dose CT scans for the evaluation of limb fractures

Aims

To evaluate whether an ultra-low-dose CT protocol can diagnose selected limb fractures as well as conventional CT (C-CT).

Patients and Methods

We prospectively studied 40 consecutive patients with a limb fracture in whom a CT scan was indicated. These were scanned using an ultra-low-dose CT Reduced Effective Dose Using Computed Tomography In Orthopaedic Injury (REDUCTION) protocol. Studies from 16 selected cases were compared with 16 C-CT scans matched for age, gender and type of fracture. Studies were assessed for diagnosis and image quality. Descriptive and reliability statistics were calculated. The total effective radiation dose for each scanned site was compared.

Results

The mean estimated effective dose (ED) for the REDUCTION protocol was 0.03 milliSieverts (mSv) and 0.43 mSv (p < 0.005) for C-CT. The sensitivity (Sn), specificity (Sp), positive predictive value (PPV) and negative predictive value (NPV) of the REDUCTION protocol to detect fractures were 0.98, 0.89, 0.98 and 0.89 respectively when two occult fractures were excluded. Inter- and intra-observer reliability for diagnosis using the REDUCTION protocol (κ = 0.75, κ = 0.71) were similar to those of C-CT (κ = 0.85, κ = 0.82). Using the REDUCTION protocol, 3D CT reconstructions were equivalent in quality and diagnostic information to those generated by C-CT (κ = 0.87, κ = 0.94).

Conclusion

With a near 14-fold reduction in estimated ED compared with C-CT, the REDUCTION protocol reduces the amount of CT radiation substantially without significant diagnostic decay. It produces images that appear to be comparable with those of C-CT for evaluating fractures of the limbs. Cite this article: Bone Joint J 2016;98-B:1668-73.

Accuracy of patient-specific organ dose estimates obtained using an automated image segmentation algorithm

The overall goal of this work is to develop a rapid, accurate, and automated software tool to estimate patient-specific organ doses from computed tomography (CT) scans using simulations to generate dose maps combined with automated segmentation algorithms. This work quantified the accuracy of organ dose estimates obtained by an automated segmentation algorithm. We hypothesized that the autosegmentation algorithm is sufficiently accurate to provide organ dose estimates, since small errors delineating organ boundaries will have minimal effect when computing mean organ dose. A leave-one-out validation study of the automated algorithm was performed with 20 head-neck CT scans expertly segmented into nine regions. Mean organ doses of the automatically and expertly segmented regions were computed from Monte Carlo-generated dose maps and compared. The automated segmentation algorithm estimated the mean organ dose to be within 10% of the expert segmentation for regions other than the spinal canal, with the median error for each organ region below 2%. In the spinal canal region, the median error was [Formula: see text], with a maximum absolute error of 28% for the single-atlas approach and 11% for the multiatlas approach. The results demonstrate that the automated segmentation algorithm can provide accurate organ dose estimates despite some segmentation errors.

Source modeling for Monte Carlo dose calculation of CT examinations with a radiotherapy treatment planning system

PURPOSE

Radiation dose to patients undergoing examinations with Multislice Computed Tomography (MSCT) as well as Cone Beam Computed Tomography (CBCT) is a matter of concern. Risk management could benefit from efficient replace rational dose calculation tools. The paper aims to verify MSCT dose calculations using a Treatment Planning System (TPS) for radiotherapy and to evaluate four different variations of bow-tie filter characterizations for the beam model used in the dose calculations.

METHODS

A TPS (RayStation™, RaySearch Laboratories, Stockholm, Sweden) was configured to calculate dose from a MSCT (GE Healthcare, Wauwatosa, WI, USA). The x-ray beam was characterized in a stationary position the by measurements of the Half-Value Layer (HVL) in aluminum and kerma along the principal axes of the isocenter plane perpendicular to the beam. A Monte Carlo source model for the dose calculation was applied with four different variations on the beam-shaping bow-tie filter, taking into account the different degrees of HVL information but reconstructing the measured kerma distribution after the bow-tie filter by adjusting the photon sampling function. The resulting dose calculations were verified by comparison with measurements in solid water as well as in an anthropomorphic phantom.

RESULTS

The calculated depth dose in solid water as well as the relative dose profiles was in agreement with the corresponding measured values. Doses calculated in the anthropomorphic phantom in the range 26-55 mGy agreed with the corresponding thermo luminescence dosimeter (TLD) measurements. Deviations between measurements and calculations were of the order of the measurement uncertainties. There was no significant difference between the different variations on the bow-tie filter modeling.

CONCLUSIONS

Under the assumption that the calculated kerma after the bow-tie filter replicates the measured kerma, the central specification of the HVL of the x-ray beam together with the kerma distribution can be used to characterize the beam. Thus, within the limits of the study, a flat bow-tie filter with an HVL specified by the vendor suffices to calculate the dose distribution. The TPS could be successfully configured to replicate the beam movement and intensity modulation of a spiral scan with dose modulation, on the basis of the specifications available in the metadata of the digital images and the log file of the CT.

Radiation Dose to the Breast by 64-slice CT: Effects of Scanner Model and Study Protocol

RATIONALE AND OBJECTIVES

This work aimed to study the effects of scanner model and study protocol on radiation dose received by breast tissues from 64-slice computed tomography (CT) studies.

MATERIALS AND METHODS

Four scanner models and three study protocols were used in scanning an anthropomorphic phantom with breast modules. Each protocol follows recommendations or guidelines from the American Association of Physicists in Medicine and the American College of Radiology. Twenty thermoluminescent dosimeters were placed inside the breast modules to measure breast tissue doses. Both the absolute and the normalized breast tissue doses were analyzed.

RESULTS

The mean glandular doses of a lung cancer screening CT, a chest/abdomen/pelvis CT, and a virtual colonoscopy CT are equivalent to less than 1, 5-7, and 1-3 two-view digital mammograms, respectively, for a standard-sized patient. The normalized breast dose differs significantly (P < 0.01) between lung cancer screening CT and chest/abdomen/pelvis CT; however, it shows less than ±10% variation among scanner models for the same protocol. In virtual colonoscopy CT, breast tissue dose decreases with the distance between local tissues to the edge of the x-ray field, although the decreasing trend varies for different scanner models and protocol settings.

CONCLUSIONS

When breasts are entirely included in the primary x-ray field, breast dose by 64-slice CT is mainly protocol dependent, with the normalized breast dose about 15% lower for protocols with modulated mA than for those with constant mA; when breasts are only partially included in the primary beam field, breast dose by 64-slice CT is dependent on both the scanner model and the protocol settings.

Influence of 320-detector-row volume scanning and AAPM report 111 CT dosimetry metrics on size-specific dose estimate: a Monte Carlo study

The American Association of Physicists in Medicine (AAPM) task group 204 has recommended the use of size-dependent conversion factors to calculate size-specific dose estimate (SSDE) values from volume computed tomography dose index (CTDIvol) values. However, these conversion factors do not consider the effects of 320-detector-row volume computed tomography (CT) examinations or the new CT dosimetry metrics proposed by AAPM task group 111. This study aims to investigate the influence of these examinations and metrics on the conversion factors reported by AAPM task group 204, using Monte Carlo simulations. Simulations were performed modelling a Toshiba Aquilion ONE CT scanner, in order to compute dose values in water for cylindrical phantoms with 8-40-cm diameters at 2-cm intervals for each scanning parameter (tube voltage, bow-tie filter, longitudinal beam width). Then, the conversion factors were obtained by applying exponential regression analysis between the dose values for a given phantom diameter and the phantom diameter combined with various scanning parameters. The conversion factors for each scanning method (helical, axial, or volume scanning) and CT dosimetry method (i.e., the CTDI100 method or the AAPM task group 111 method) were in agreement with those reported by AAPM task group 204, within a percentage error of 14.2 % for phantom diameters ≥11.2 cm. The results obtained in this study indicate that the conversion factors previously presented by AAPM task group 204 can be used to provide appropriate SSDE values for 320-detector-row volume CT examinations and the CT dosimetry metrics proposed by the AAPM task group 111.

Organ doses can be estimated from the computed tomography (CT) dose index for cone-beam CT on radiotherapy equipment

Cone beam computed tomography (CBCT) systems are fitted to radiotherapy linear accelerators and used for patient positioning prior to treatment by image guided radiotherapy (IGRT). Radiotherapists' and radiographers' knowledge of doses to organs from CBCT imaging is limited. The weighted CT dose index for a reference beam of width 20 mm (CTDIw,ref) is displayed on Varian CBCT imaging equipment known as an On-Board Imager (OBI) linked to the Truebeam linear accelerator. This has the potential to provide an indication of organ doses. This knowledge would be helpful for guidance of radiotherapy clinicians preparing treatments. Monte Carlo simulations of imaging protocols for head, thorax and pelvic scans have been performed using EGSnrc/BEAMnrc, EGSnrc/DOSXYZnrc, and ICRP reference computational male and female phantoms to derive the mean absorbed doses to organs and tissues, which have been compared with values for the CTDIw,ref displayed on the CBCT scanner console. Substantial variations in dose were observed between male and female phantoms. Nevertheless, the CTDIw,ref gave doses within  ±21% for the stomach and liver in thorax scans and 2  ×  CTDIw,ref can be used as a measure of doses to breast, lung and oesophagus. The CTDIw,ref could provide indications of doses to the brain for head scans, and the colon for pelvic scans. It is proposed that knowledge of the link between CTDIw for CBCT should be promoted and included in the training of radiotherapy staff.

Comparison of Measured and Estimated CT Organ Doses for Modulated and Fixed Tube Current:: A Human Cadaver Study

RATIONALE AND OBJECTIVES

The aim of this study was to compare the directly measured and the estimated computed tomography (CT) organ doses obtained from commercial radiation dose-tracking (RDT) software for CT performed with modulated tube current or automatic exposure control (AEC) technique and fixed tube current (mAs).

MATERIALS AND METHODS

With the institutional review board (IRB) approval, the ionization chambers were surgically implanted in a human cadaver (88 years old, male, 68 kg) in six locations such as liver, stomach, colon, left kidney, small intestine, and urinary bladder. The cadaver was scanned with routine abdomen pelvis protocol on a 128-slice, dual-source multidetector computed tomography (MDCT) scanner using both AEC and fixed mAs. The effective and quality reference mAs of 100, 200, and 300 were used for AEC and fixed mAs, respectively. Scanning was repeated three times for each setting, and measured and estimated organ doses (from RDT software) were recorded (N = 3*3*2 = 18).

RESULTS

Mean CTDIvol for AEC and fixed mAs were 4, 8, 13 mGy and 7, 14, 21 mGy, respectively. The most estimated organ doses were significantly greater (P < 0.01) than the measured organ doses for both AEC and fixed mAs. At AEC, the mean estimated organ doses (for six organs) were 14.7 mGy compared to mean measured organ doses of 12.3 mGy. Similarly, at fixed mAs, the mean estimated organ doses (for six organs) were 24 mGy compared to measured organ doses of 22.3 mGy. The differences among the measured and estimated organ doses were higher for AEC technique compared to the fixed mAs for most organs (P < 0.01).

CONCLUSIONS

The most CT organ doses estimated from RDT software are greater compared to directly measured organ doses, particularly when AEC technique is used for CT scanning.

Diagnostic Reference Levels of CT Radiation Dose in Whole-Body PET/CT

The role of CT in PET/CT imaging includes acquisition techniques for diagnostic, anatomic localization, and attenuation correction purposes. Diagnostic reference levels of the volumetric CT dose index (CTDIvol) are available for dedicated CT procedures on selected body regions, but similar reference levels for whole-body CT used in PET/CT examinations are limited. This work reports CTDIvol values from sites that conduct whole-body oncologic PET/CT examinations and participated in the scanner validation program of the Society of Nuclear Medicine and Molecular Imaging Clinical Trials Network.

METHODS

From 2010 to 2014, a total of 154 sites submitted CT acquisition parameters used in their clinical (18)F-FDG PET/CT oncology protocols. From these parameters, the CTDIvol was estimated using the ImPACT CTDI dosimetry tables. Histograms of CTDIvol values were created for each year, and descriptive statistics, including mean, median, and 75th percentile, were reported. Repeated-measures ANOVA was performed to determine whether significant differences occurred between reporting years.

RESULTS

A wide range of technical parameters was reported, most notably in tube current. Between 2010 and 2014, the median CTDIvol ranged from 4.9 to 6.2 mGy and the 75th percentile from 9.7 to 10.2 mGy. There was no significant change in CTDIvol between reporting years (repeated-measures ANOVA, P = 0.985).

CONCLUSION

The 75th percentile CTDIvol reported in this work was 9.8 mGy averaged over all reporting years. These data provide a resource for establishing CTDIvol reference values specific to performing CT in PET/CT whole-body examinations. The wide ranges of CT acquisition parameters reported by sites suggest that CTDIvol reference levels may be beneficial for optimization of CT protocols.

Patient-specific dose calculations for pediatric CT of the chest, abdomen and pelvis

BACKGROUND

Organ dose is essential for accurate estimates of patient dose from CT.

OBJECTIVE

To determine organ doses from a broad range of pediatric patients undergoing diagnostic chest-abdomen-pelvis CT and investigate how these relate to patient size.

MATERIALS AND METHODS

We used a previously validated Monte Carlo simulation model of a Philips Brilliance 64 multi-detector CT scanner (Philips Healthcare, Best, The Netherlands) to calculate organ doses for 40 pediatric patients (M:F = 21:19; range 0.6-17 years). Organ volumes and positions were determined from the images using standard segmentation techniques. Non-linear regression was performed to determine the relationship between volume CT dose index (CTDIvol)-normalized organ doses and abdominopelvic diameter. We then compared results with values obtained from independent studies.

RESULTS

We found that CTDIvol-normalized organ dose correlated strongly with exponentially decreasing abdominopelvic diameter (R(2) > 0.8 for most organs). A similar relationship was determined for effective dose when normalized by dose-length product (R(2) = 0.95). Our results agreed with previous studies within 12% using similar scan parameters (e.g., bowtie filter size, beam collimation); however results varied up to 25% when compared to studies using different bowtie filters.

CONCLUSION

Our study determined that organ doses can be estimated from measurements of patient size, namely body diameter, and CTDIvol prior to CT examination. This information provides an improved method for patient dose estimation.

Attenuation-based size metric for estimating organ dose to patients undergoing tube current modulated CT exams

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.

Organ and effective dose evaluation in coronary angiography by using a 320 MDCT based on in-phantom dose measurements with TLDs

The purpose of this study is to estimate the organ and effective dose (E) from computed tomography coronary angiography (CTCA) on a 320-MDCT scanner. Radiation dose was estimated for the prospectively ECG-gated CTCA in a male phantom. A total of 451 chips of thermoluminescent dosimeter were implanted in the phantom for measuring the organ doses. The effective doses were calculated using measured organ doses and the tissue-weighting factors. The dose length product (DLP) values were recorded and used to develop the conversion coefficient k = 0.017 mSv•(mGy•cm)(-1) (E/DLP). In a 3-beat acquisition, the organ doses ranged from 0.24 to 71.55 mGy, and the doses in breast, bone surface, oesophagus, and lung were higher than 20 mGy. The effective doses in 2-beat and 3-beat acquisition were estimated to be 14.3 and 24.3 mSv. More beats of acquisition led to higher radiation dose. The reported k values for chest CT scan can be used to roughly estimate the E value from CTCA for 320 MDCT.

Deterministic absorbed dose estimation in computed tomography using a discrete ordinates method

PURPOSE

Organ dose estimation for a patient undergoing computed tomography (CT) scanning is very important. Although Monte Carlo methods are considered gold-standard in patient dose estimation, the computation time required is formidable for routine clinical calculations. Here, the authors instigate a deterministic method for estimating an absorbed dose more efficiently.

METHODS

Compared with current Monte Carlo methods, a more efficient approach to estimating the absorbed dose is to solve the linear Boltzmann equation numerically. In this study, an axial CT scan was modeled with a software package, Denovo, which solved the linear Boltzmann equation using the discrete ordinates method. The CT scanning configuration included 16 x-ray source positions, beam collimators, flat filters, and bowtie filters. The phantom was the standard 32 cm CT dose index (CTDI) phantom. Four different Denovo simulations were performed with different simulation parameters, including the number of quadrature sets and the order of Legendre polynomial expansions. A Monte Carlo simulation was also performed for benchmarking the Denovo simulations. A quantitative comparison was made of the simulation results obtained by the Denovo and the Monte Carlo methods.

RESULTS

The difference in the simulation results of the discrete ordinates method and those of the Monte Carlo methods was found to be small, with a root-mean-square difference of around 2.4%. It was found that the discrete ordinates method, with a higher order of Legendre polynomial expansions, underestimated the absorbed dose near the center of the phantom (i.e., low dose region). Simulations of the quadrature set 8 and the first order of the Legendre polynomial expansions proved to be the most efficient computation method in the authors' study. The single-thread computation time of the deterministic simulation of the quadrature set 8 and the first order of the Legendre polynomial expansions was 21 min on a personal computer.

CONCLUSIONS

The simulation results showed that the deterministic method can be effectively used to estimate the absorbed dose in a CTDI phantom. The accuracy of the discrete ordinates method was close to that of a Monte Carlo simulation, and the primary benefit of the discrete ordinates method lies in its rapid computation speed. It is expected that further optimization of this method in routine clinical CT dose estimation will improve its accuracy and speed.

Assessment of patient dose reduction by bismuth shielding in CT using measurements, GEANT4 and MCNPX simulations

This work reports on the use of two different Monte Carlo codes (GEANT4 and MCNPX) for assessing the dose reduction using bismuth shields in computer tomography (CT) procedures in order to protect radiosensitive organs such as eye lens, thyroid and breast. Measurements were performed using head and body PMMA phantoms and an ionisation chamber placed in five different positions of the phantom. Simulations were performed to estimate Computed Tomography Dose Index values using GEANT4 and MCNPX. The relative differences between measurements and simulations were <10 %. The dose reduction arising from the use of bismuth shielding ranges from 2 to 45 %, depending on the position of the bismuth shield. The percentage of dose reduction was more significant for the area covered by the bismuth shielding (36 % for eye lens, 39 % for thyroid and 45 % for breast shields).

Validation of a Monte Carlo model used for simulating tube current modulation in computed tomography over a wide range of phantom conditions/challenges

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.

New absorbed dose measurement with cylindrical water phantoms for multidetector CT

The aim of this study was to develop new dosimetry with cylindrical water phantoms for multidetector computed tomography (MDCT). The ionization measurement was performed with a Farmer ionization chamber at the center and four peripheral points in the body-type and head-type cylindrical water phantoms. The ionization was converted to the absorbed dose using a (60)Co absorbed-dose-to-water calibration factor and Monte Carlo (MC) -calculated correction factors. The correction factors were calculated from MDCT (Brilliance iCT, 64-slice, Philips Electronics) modeled with GMctdospp (IMPS, Germany) software based on the EGSnrc MC code. The spectrum of incident x-ray beams and the configuration of a bowtie filter for MDCT were determined so that calculated photon intensity attenuation curves for aluminum (Al) and calculated off-center ratio (OCR) profiles in air coincided with those measured. The MC-calculated doses were calibrated by the absorbed dose measured at the center in both cylindrical water phantoms. Calculated doses were compared with measured doses at four peripheral points and the center in the phantom for various beam pitches and beam collimations. The calibration factors and the uncertainty of the absorbed dose determined using this method were also compared with those obtained by CTDIair (CT dose index in air). Calculated Al half-value layers and OCRs in air were within 0.3% and 3% agreement with the measured values, respectively. Calculated doses at four peripheral points and the centers for various beam pitches and beam collimations were within 5% and 2% agreement with measured values, respectively. The MC-calibration factors by our method were 44-50% lower than values by CTDIair due to the overbeaming effect. However, the calibration factors for CTDIair agreed within 5% with those of our method after correction for the overbeaming effect. Our method makes it possible to directly measure the absorbed dose for MDCT and is more robust and accurate than the CTDIair measurement.

Development and validation of a GEANT4 radiation transport code for CT dosimetry

The authors have created a radiation transport code using the GEANT4 Monte Carlo toolkit to simulate pediatric patients undergoing CT examinations. The focus of this paper is to validate their simulation with real-world physical dosimetry measurements using two independent techniques. Exposure measurements were made with a standard 100-mm CT pencil ionization chamber, and absorbed doses were also measured using optically stimulated luminescent (OSL) dosimeters. Measurements were made in air with a standard 16-cm acrylic head phantom and with a standard 32-cm acrylic body phantom. Physical dose measurements determined from the ionization chamber in air for 100 and 120 kVp beam energies were used to derive photon-fluence calibration factors. Both ion chamber and OSL measurement results provide useful comparisons in the validation of the Monte Carlo simulations. It was found that simulated and measured CTDI values were within an overall average of 6% of each other.

Compliance with AAPM Practice Guideline 1.a: CT Protocol Management and Review - from the perspective of a university hospital

The purpose of this paper is to describe our experience with the AAPM Medical Physics Practice Guideline 1.a: “CT Protocol Management and Review Practice Guideline”. Specifically, we will share how our institution's quality management system addresses the suggestions within the AAPM practice report. We feel this paper is needed as it was beyond the scope of the AAPM practice guideline to provide specific details on fulfilling individual guidelines. Our hope is that other institutions will be able to emulate some of our practices and that this article would encourage other types of centers (e.g., community hospitals) to share their methodology for approaching CT protocol optimization and quality control. Our institution had a functioning CT protocol optimization process, albeit informal, since we began using CT. Recently, we made our protocol development and validation process compliant with a number of the ISO 9001:2008 clauses and this required us to formalize the roles of the members of our CT protocol optimization team. We rely heavily on PACS‐based IT solutions for acquiring radiologist feedback on the performance of our CT protocols and the performance of our CT scanners in terms of dose (scanner output) and the function of the automatic tube current modulation. Specific details on our quality management system covering both quality control and ongoing optimization have been provided. The roles of each CT protocol team member have been defined, and the critical role that IT solutions provides for the management of files and the monitoring of CT protocols has been reviewed. In addition, the invaluable role management provides by being a champion for the project has been explained; lack of a project champion will mitigate the efforts of a CT protocol optimization team. Meeting the guidelines set forth in the AAPM practice guideline was not inherently difficult, but did, in our case, require the cooperation of radiologists, technologists, physicists, IT, administrative staff, and hospital management. Some of the IT solutions presented in this paper are novel and currently unique to our institution.

PACS number: 87.57.Q

Comprehensive assessment of radiation dose estimates for the CORE320 study

OBJECTIVE. The purpose of this study was to comprehensively study estimated radiation doses for subjects included in the main analysis of the Combined Non-invasive Coronary Angiography and Myocardial Perfusion Imaging Using 320 Detector Computed Tomography (CORE320) study ( ClinicalTrials.gov identifier NCT00934037), a clinical trial comparing combined CT angiography (CTA) and perfusion CT with the reference standard catheter angiography plus myocardial perfusion SPECT. SUBJECTS AND METHODS. Prospectively acquired data on 381 CORE320 subjects were analyzed in four groups of testing related to radiation exposure. Radiation dose estimates were compared between modalities for combined CTA and perfusion CT with respect to covariates known to influence radiation exposure and for the main clinical outcomes defined by the trial. The final analysis assessed variations in radiation dose with respect to several factors inherent to the trial. RESULTS. The mean radiation dose estimate for the combined CTA and perfusion CT protocol (8.63 mSv) was significantly (p < 0.0001 for both) less than the average dose delivered from SPECT (10.48 mSv) and the average dose from diagnostic catheter angiography (11.63 mSv). There was no significant difference in estimated CTA-perfusion CT radiation dose for subjects who had false-positive or false-negative results in the CORE320 main analyses in a comparison with subjects for whom the CTA-perfusion CT findings were in accordance with the reference standard SPECT plus catheter angiographic findings. CONCLUSION. Radiation dose estimates from CORE320 support clinical implementation of a combined CT protocol for assessing coronary anatomy and myocardial perfusion.

Accuracy of Monte Carlo simulations compared to in-vivo MDCT dosimetry

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.

Calculating patient-specific organ doses from adult body CT scans by Monte Carlo analysis using male-individual voxel phantoms

Computed tomography (CT) dose calculators such as the WAZA-ARI are useful for estimating the radiation dose from CT examination. This study determined correction coefficients for estimating organ doses to patients of any size attended to in daily clinical practice. To this end, the authors constructed voxel phantoms based on the CT images of patients of different size and simulated radiation transport in CT examinations to obtain organ doses using Monte Carlo simulation. The results show that the linear relationship between effective diameter and organ dose can predict patient-specific organ doses. The effective diameter-based prediction can provide accuracy within an error of ±10%, whereas an error of >20% was observed only in the liver and bladder. These results may contribute to practical irradiation dose calculation from a CT examination depending on individual patient size within a certain degree of accuracy.