Patient-specific Monte Carlo-based organ dose estimates in spiral CT via optical 3D body scanning and adaptation of a voxelized phantom dataset: proof-of-principle

Objective. We present a method for personalized organ dose estimates obtained before the CT exam, via 3D optical body scanning and Monte Carlo (MC) simulations.Approach. A voxelized phantom is derived by adapting a reference phantom to the body size and shape measured with a portable 3D optical scanner, which returns the 3D silhouette of the patient. This was used as an external rigid envelope for incorporating a tailored version of the internal body anatomy derived from a phantom dataset (National Cancer Institute, NIH, USA) matched for gender, age, weight, and height. The proof-of-principle was conducted on adult head phantoms. The Geant4 MC code provided estimates of the organ doses from 3D absorbed dose maps in voxelized body phantom.Main results. We applied this approach for head CT scanning using an anthropomorphic voxelized head phantom derived from 3D optical scans of mannequins. We compared the estimates of head organ doses with those provided by the NCICT3.0 software (NCI, NIH, USA). Head organ doses differed up to 38% using the proposed personalized estimate and MC code, with respect to corresponding estimates calculated for the standard (non-personalized) reference head phantom. Preliminary application of the MC code to chest CT scans is shown. Real-time pre-exam personalized CT dosimetry is envisaged with adoption of a GPU-based fast MC code.Significance. The developed procedure for personalized organ dose estimates before the CT exam, introduces a new approach for realistic description of size and shape of patients via voxelized phantoms specific for each patient.

A survey of local diagnostic reference levels for the head, thorax, abdomen and pelvis computed tomography in Norway and Canada

Background

Computed tomography (CT) contributes to 60% of the collective dose in medical imaging. Literature has demonstrated that patient dose varies across regions and countries. Establishing diagnostic reference levels (DRLs) contributes to the optimization of clinical practices and radiation protection.

Purpose

To survey the dose indices (CTDIvol and dose-length product) for frequently performed CT examinations from the chosen hospitals in Norway and Canada and to determine local DRLs (LDRLs) based on the collected data.

Material and Methods

The survey included eight scanners from two Norwegian hospitals and four scanners from four Canadian hospitals. Dosimetry data were collected for the following routine CT examinations: head, contrast-enhanced thorax, and abdomen and pelvis. Overall 480 adult average-sized patients from Norway and 360 from Canada were included in the survey. The LDRLs were determined as the 75th percentile of distributions of median values of dose indicators from different CT scanners. The differences in dose between scanners were determined using single-factor ANOVA.

Results

The LDRLs determined in Norway were higher overall than in Canada. The obtained values were compared to the national DRLs. The dose from several scanners in Norway exceeded national Norwegian DRLs, while Canadian LDRLs were below the Canadian reference levels. The differences between the means of the dose distributions from each scanner were statistically significant (p < 0.05) for all examinations with exception of identical scanners located in the same hospital and using the same protocols.

Conclusion

Observed dose variations even in the same hospital, or from the same scanner model confirmed the need for CT protocol optimization.

A comprehensive Monte Carlo study of CT dose metrics proposed by the AAPM Reports 111 and 200

PURPOSE

A Monte Carlo (MC) modeling of single axial and helical CT scan modes has been developed to compute single and accumulated dose distributions. The radiation emission characteristics of an MDCT scanner has been modeled and used to evaluate the dose deposition in infinitely long head and body PMMA phantoms. The simulated accumulated dose distributions determined the approach to equilibrium function, H(L). From these H ( L ) curves, dose-related information was calculated for different head and body clinical protocols.

METHODS

The PENELOPE/penEasy package has been used to model the single axial and helical procedures and the radiation transport of photons and electrons in the phantoms. The bowtie filters, heel effect, focal-spot angle, and fan-beam geometry were incorporated. Head and body protocols with different pitch values were modeled for x-ray spectra corresponding to 80, 100, 120, and 140 kV. The analytical formulation for the single dose distributions and experimental measurements of single and accumulated dose distributions were employed to validate the MC results. The experimental dose distributions were measured with OSLDs and a thimble ion chamber inserted into PMMA phantoms. Also, the experimental values of the C T D I 100 along the center and peripheral axes of the CTDI phantom served to calibrate the simulated single and accumulated dose distributions.

RESULTS

The match of the simulated dose distributions with the reference data supports the correct modeling of the heel effect and the radiation transport in the phantom material reflected in the tails of the dose distributions. The validation of the x-ray source model was done comparing the CTDI ratios between simulated, measured and CTDosimetry data. The average difference of these ratios for head and body protocols between the simulated and measured data was in the range of 13-17% and between simulated and CTDosimetry data varied 10-13%. The distributions of simulated doses and those measured with the thimble ion chamber are compatible within 3%. In this study, it was demonstrated that the efficiencies of the C T D I 100 measurements in head phantoms with nT = 20 mm and 120 kV are 80.6% and 87.8% at central and peripheral axes, respectively. In the body phantoms with n T = 40 mm and 120 kV, the efficiencies are 56.5% and 86.2% at central and peripheral axes, respectively. In general terms, the clinical parameters such as pitch, beam intensity, and voltage affect the D_eq values with the increase of the pitch decreasing the D_eq and the beam intensity and the voltage increasing its value. The H(L) function does not change with the pitch values, but depends on the phantom axis (central or peripheral).

CONCLUSIONS

The computation of the pitch-equilibrium dose product, D ̂ eq , evidenced the limitations of the C T D I 100 method to determine the dose delivered by a CT scanner. Therefore, quantities derived from the C T D I 100 propagate this limitation. The developed MC model shows excellent compatibility with both measurements and literature quantities defined by AAPM Reports 111 and 200. These results demonstrate the robustness and versatility of the proposed modeling method.

Dose quantities for measurement and comparison of doses to individual patients in computed tomography (CT)

The dose quantities displayed routinely on CT scanners, the volume averaged CT dose index (CTDI_vol) and dose length product, provide measures of doses calculated for standard phantoms. The American Association of Medical Physics has published conversion factors for the adjustment of CTDI_volto take account of variations in patient size, the results being termed size-specific dose estimate (SSDE). However, CTDI_voland SSDE, while useful in comparing and optimising doses from a set procedure, do not provide risk-related information that takes account of the organs and tissues irradiated and associated cancer risks. A derivative of effective dose that takes account of differences in body and organ sizes and masses, referred to here as size-specific effective dose (SED), can provide such information. Data on organ doses from NCICT software that is based on Monte Carlo simulations of CT scans for 193 adult phantoms have been used to compute values of SED for CT examinations of the trunk and results compared with corresponding values of SSDE. Relationships within ±8% were observed between SED and SSDE for scans extending over similar regions for phantoms with a wide range of sizes. Coefficients have been derived from fits of the data to estimate SED values from SSDEs for different regions of the body for scans of standard lengths based on patient height. A method developed to take account of differences in scan length gave SED results within ±5% of values calculated using the NCI phantom library. This approach could potentially be used to estimate SED from SSDE values, allowing their display at the time a CT scan is performed.

Validation of a method for estimating peak skin dose from CT-guided procedures

A method for estimating peak skin dose (PSD) from CTDI_vol has been published but not validated. The objective of this study was to validate this method during CT‐guided ablation procedures. Radiochromic film was calibrated and used to measure PSD. Sixty‐eight patients were enrolled in this study, and measured PSD were collected for 46 procedures. CTDI_vol stratified by axial and helical scanning was used to calculate an estimate of PSD using the method [1.2 × CTDI_vol(helical) + 0.6 × CTDI_vol(axial)], and both calculated PSD and total CTDI_vol were compared to measured PSD using paired t‐tests on the log‐transformed data and Bland‐Altman analysis. Calculated PSD were significantly different from measured PSD (P < 0.0001, bias, 18.3%, 95% limits of agreement, −63.0% to 26.4%). Measured PSD were not significantly different from total CTDI_vol (P = 0.27, bias, 3.97%, 95% limits of agreement, −51.6% to 43.7%). Considering that CTDI_vol is reported on the console of all CT scanners, is not stratified by axial and helical scanning modes, and is immediately available to the operator during CT‐guided interventional procedures, it may be reasonable to use the scanner‐reported CTDI_vol as an indicator of PSD during CT‐guided procedures. However, further validation is required for other models of CT scanner.

Accuracy of weighted CTDI in estimating average dose delivered to CTDI phantoms: An experimental study

PURPOSE

The concept of the weighted computed tomography dose index ( CTDI w ) was proposed in 1995 to represent the average CTDI across an axial section of a cylindrical phantom. The purpose of this work was to experimentally re-examine the validity of the underlying assumptions behind CTDI w for modern MDCT systems.

METHODS

To enable experimental mapping of CTDI 100 in the axial plane, in-house 16 and 32 cm cylindrical phantoms were fabricated to allow the pencil chamber to reach any arbitrary axial location within the phantoms. The phantoms were scanned on a clinical MDCT with five beam collimation widths, three bowtie filters, and four kV levels. To evaluate the linearity and rotational invariance assumptions implicitly made when the weighting factors of 1/3 and 2/3 in the CTDI w formula were originally derived, CTDI 100 was measured at different radial and angular locations within the phantom for different collimation, bowtie, and kV combinations. The average CTDI ( CTDI avg ) across the axial plane was calculated from the experimental two-dimensional (2D) dose distribution and was compared with the traditional CTDI w .

RESULTS

For both phantoms under all scan conditions, the axial dose distributions were found to have significant angular dependence, potentially due to the x-ray attenuation by the patient couch or the head holder. The radial dose profiles were also found to significantly deviate from linearity in many cases due to the presence of the bowtie filter. When only the 12 o'clock peripheral CTDI 100 and the traditional weighting factors were used to calculate CTDI w , the average dose was overestimated in the 16 cm phantom by up to 8.4% at isocenter and up to 35.3% when the phantom was off-centered by 6 cm; in the 32 cm phantom at isocenter, the average dose was overestimated by up to 12.8%. Using an average of the four peripheral CTDI 100 measurements at the 12, 3, 6, and 9 o'clock locations reduced the error of CTDI w to within 1.2% in the 16 cm phantom. For the 32 cm phantom, even by using the average of the peripheral measurements, the traditional CTDI w underestimated the average dose by up to 4.3% due to aggressive drop-off of the CTDI 100 at the phantom periphery.

CONCLUSIONS

The linearity and rotational-invariance assumptions behind the traditional CTDI w formalism may not be valid for modern CT systems and thus CTDI w may not accurately represent the average dose or radiation output within a CTDI phantom. Utilizing data from all four peripheral locations always improves accuracy of CTDI w in representing the true average dose. For the large (32 cm) phantom, nonlinear models and more measurement points are needed if a more precise estimation of the average axial dose is required.

Dosimetric and radiation cancer risk evaluation of high resolution thorax CT during COVID-19 outbreak

PURPOSE

The aim of this work was to evaluate the dosimetric impact of high-resolution thorax CT during COVID-19 outbreak in the University Hospital of Parma. In two months we have performed a huge number of thorax CT scans collecting effective and equivalent organ doses and evaluating also the lifetime attributable risk (LAR) of lung and other major cancers.

MATERIALS AND METHOD

From February 24th to April 28th, 3224 high-resolution thorax CT were acquired. For all patients we have examined the volumetric computed tomography dose index (CTDIvol), the dose length product (DLP), the size-specific dose estimate (SSDE) and effective dose (E103) using a dose tracking software (Radimetrics Bayer HealthCare). From the equivalent dose to organs for each patient, LAR for lung and major cancers were estimated following the method proposed in BEIR VII which considers age and sex differences.

RESULTS

Study population included 3224 patients, 1843 male and 1381 female, with an average age of 67 years. The average CTDIvol, SSDE and DLP, and E103 were 6.8 mGy, 8.7 mGy, 239 mGy·cm and 4.4 mSv respectively. The average LAR of all solid cancers was 2.1 cases per 10,000 patients, while the average LAR of leukemia was 0.2 cases per 10,000 patients. For both male and female the organ with a major cancer risk was lung.

CONCLUSIONS

Despite the impressive increment in thoracic CT examinations due to COVID-19 outbreak, the high resolution low dose protocol used in our hospital guaranteed low doses and very low risk estimation in terms of LAR.

Patient-adapted organ absorbed dose and effective dose estimates in pediatric 18F-FDG positron emission tomography/computed tomography studies

BACKGROUND

Organ absorbed doses and effective doses can be used to compare radiation exposure among medical imaging procedures, compare alternative imaging options, and guide dose optimization efforts. Individual dose estimates are important for relatively radiosensitive patient populations such as children and for radiosensitive organs such as the eye lens. Software-based dose calculation methods conveniently calculate organ dose using patient-adjusted and examination-specific inputs.

METHODS

Organ absorbed doses and effective doses were calculated for 429 pediatric 18F-FDG PET-CT patients. Patient-adjusted and scan-specific information was extracted from the electronic medical record and scanner dose-monitoring software. The VirtualDose and OLINDA/EXM (version 2.0) programs, respectively, were used to calculate the CT and the radiopharmaceutical organ absorbed doses and effective doses. Patients were grouped according to age at the time of the scan as follows: less than 1 year old, 1 to 5 years old, 6 to 10 years old, 11 to 15 years old, and 16 to 17 years old.

RESULTS

The mean (+/- standard deviation, range) total PET plus CT effective dose was 14.5 (1.9, 11.2-22.3) mSv. The mean (+/- standard deviation, range) PET effective dose was 8.1 (1.2, 5.7-16.5) mSv. The mean (+/- standard deviation, range) CT effective dose was 6.4 (1.8, 2.9-14.7) mSv. The five organs with highest PET dose were: Urinary bladder, heart, liver, lungs, and brain. The five organs with highest CT dose were: Thymus, thyroid, kidneys, eye lens, and gonads.

CONCLUSIONS

Organ and effective dose for both the CT and PET components can be estimated with actual patient and scan data using commercial software. Doses calculated using software generally agree with those calculated using dose conversion factors, although some organ doses were found to be appreciably different. Software-based dose calculation methods allow patient-adjusted dose factors. The effort to gather the needed patient data is justified by the resulting value of the characterization of patient-adjusted dosimetry.

Impact of patient centering in CT on organ dose and the effect of using a positioning compensation system: Evidence from OSLD measurements in postmortem subjects

The purpose of this study was to investigate the frequency and impact of vertical mis-centering on organ doses in computed tomography (CT) exams and evaluate the effect of a commercially available positioning compensation system (PCS). Mis-centering frequency and magnitude was retrospectively measured in 300 patients examined with chest-abdomen-pelvis CT. Organ doses were measured in three postmortem subjects scanned on a CT scanner at nine different vertical table positions (maximum shift ± 4 cm). Organ doses were measured with optically stimulated luminescent dosimeters inserted within organs. Regression analysis was performed to determine the correlation between organ doses and mis-centering. Methods were repeated using a PCS that automatically detects the table offset to adjust tube current output accordingly. Clinical mis-centering was >1 cm in 53% and 21% of patients in the vertical and lateral directions, respectively. The 1-cm table shifts resulted in organ dose differences up to 8%, while 4-cm shifts resulted in organ dose differences up to 35%. Organ doses increased linearly with superior table shifts for the lung, colon, uterus, ovaries, and skin (R²  = 0.73-0.99, P < 0.005). When the PCS was utilized, organ doses decreased with superior table shifts and dose differences were lower (average 5%, maximum 18%) than scans performed without PCS (average 9%, maximum 35%) at all table shifts. Mis-centering occurs frequently in the clinic and has a significant effect on patient dose. While accurate patient positioning remains important for maintaining optimal imaging conditions, a PCS has been shown to reduce the effects of patient mis-centering.

CT dosimetry at the Australian Synchrotron for 25-100 keV photons and 35-160 mm-diameter biological specimens

The dose length product (DLP) method for medical computed tomography (CT) dosimetry is applied on the Australian Synchrotron Imaging and Medical Beamline (IMBL). Beam quality is assessed from copper transmission measurements using image receptors, finding near 100% (20 keV), 3.3% (25 keV) and 0.5% (30-40 keV) relative contributions from third-harmonic radiation. The flat-panel-array medical image receptor is found to have a non-linear dose response curve. The amount of radiation delivered during an axial CT scan is measured as the dose in air alone, and inside cylindrical PMMA phantoms with diameters 35-160 mm for mono-energetic radiation 25-100 keV. The radiation output rate for the IMBL is comparable with that used for medical CT. Results are presented as the ratios of CT dose indices (CTDI) inside phantoms to in air with no phantom. Ratios are compared for the IMBL against medical CT where bow-tie filters shape the beam profile to reduce the absorbed dose to surface organs. CTDI ratios scale measurements in air to estimate the volumetric CTDI representing the average dose per unit length, and the dose length product representing the absorbed dose to the scanned volume. Medical CT dose calculators use the DLP, beam quality, axial collimation and helical pitch to estimate organ doses and the effective dose. The effective dose per unit DLP for medical CT is presented as a function of body region, beam energy and sample sizes from neonate to adult.

Measurements of air kerma index in computed tomography: A comparison among methodologies

As CT exams impart high doses to patients in comparison to other radiologist techniques, reliable dosimetry is required. In this work, dosimetry in CT beams was carried out in terms of air kerma index in air or in a phantom measured by a pencil ionization chamber, thermoluminescent (TL) detectors and radiochromic film. Calibration results showed the low energy dependence of all three dosimetric systems for the 100-120kV range, the very high uncertainty of the TL dosimeters in comparison to the other systems and high nonlinearity response in terms of air kerma of the radiochromic film. Measurements with the three systems in a 120kV CT protocol showed an acceptable agreement among the weighted air kerma index values, but TL dosimeters presented the highest uncertainties in the values.

A comparison of pediatric and adult CT organ dose estimation methods

BACKGROUND

Computed Tomography (CT) contributes up to 50% of the medical exposure to the United States population. Children are considered to be at higher risk of developing radiation-induced tumors due to the young age of exposure and increased tissue radiosensitivity. Organ dose estimation is essential for pediatric and adult patient cancer risk assessment. The objective of this study is to validate the VirtualDose software in comparison to currently available software and methods for pediatric and adult CT organ dose estimation.

METHODS

Five age groups of pediatric patients and adult patients were simulated by three organ dose estimators. Head, chest, abdomen-pelvis, and chest-abdomen-pelvis CT scans were simulated, and doses to organs both inside and outside the scan range were compared. For adults, VirtualDose was compared against ImPACT and CT-Expo. For pediatric patients, VirtualDose was compared to CT-Expo and compared to size-based methods from literature. Pediatric to adult effective dose ratios were also calculated with VirtualDose, and were compared with the ranges of effective dose ratios provided in ImPACT.

RESULTS

In-field organs see less than 60% difference in dose between dose estimators. For organs outside scan range or distributed organs, a five times' difference can occur. VirtualDose agrees with the size-based methods within 20% difference for the organs investigated. Between VirtualDose and ImPACT, the pediatric to adult ratios for effective dose are compared, and less than 21% difference is observed for chest scan while more than 40% difference is observed for head-neck scan and abdomen-pelvis scan. For pediatric patients, 2 cm scan range change can lead to a five times dose difference in partially scanned organs.

CONCLUSIONS

VirtualDose is validated against CT-Expo and ImPACT with relatively small discrepancies in dose for organs inside scan range, while large discrepancies in dose are observed for organs outside scan range. Patient-specific organ dose estimation is possible using the size-based methods, and VirtualDose agrees with size-based method for the organs investigated. Careful range selection for CT protocols is necessary for organ dose optimization for pediatric and adult patients.

Imaging wisely: patient safety in CT

The past decade has seen a significant growth in diagnostic CT imaging as a direct result of the clinical value provided by CT imaging. At the same time, many new techniques and resources are now available to make CT imaging safe. This article presents the basics of CT dosimetry and their usage in clinical practices, methods to implement CT dose reduction, followed by a summary of legislation, and guidelines related to patient safety in diagnostic CT imaging. Also, CT radiation dose diagnostic reference levels from published regional and national surveys are reviewed and applied in a CT dose tracking and monitoring program.

Comparison of pencil-type ionization chamber calibration results and methods between dosimetry laboratories

A comparison of calibration results and procedures in terms of air kerma length product, PKL, and air kerma, K, was conducted between eight dosimetry laboratories. A pencil-type ionization chamber (IC), generally used for computed tomography dose measurements, was calibrated according to three calibration methods, while its residual signal and other characteristics (sensitivity profile, active length) were assessed. The results showed that the "partial irradiation method" is the preferred method for the pencil-type IC calibration in terms of PKL and it could be applied by the calibration laboratories successfully. Most of the participating laboratories achieved high level of agreement (>99%) for both dosimetry quantities (PKL and K). Estimated relative standard uncertainties of comparison results vary among laboratories from 0.34% to 2.32% depending on the quantity, beam quality and calibration method applied. Detailed analysis of the assigned uncertainties have been presented and discussed.

Dosimetric characterization and image quality evaluation of the AIRO mobile CT scanner

Radiation dose and image quality from a recently introduced mobile CT imaging system are presented. Radiation dose was measured using a conventional 100 mm pencil ionization chamber and CT polymethylmetacrylate (PMMA) body and head phantoms. Image quality was evaluated with a CATPHAN 500 phantom. Spatial resolution, low contrast resolution, Modulation Transfer Function (MTF), and Normalized Noise Power Spectrum (NNPS) were analyzed. Radiation dose and image quality were compared to those from a multi-detector CT scanner (Siemens Sensation 64). Under identical technique factors radiation dose (mGy/mAs) from the AIRO mobile CT system (AIRO) is higher than that from a 64 slice CT scanner. Based on MTF analysis, both Soft and Standard filters of the AIRO system lost resolution quickly compared to the Sensation 64 slice CT. The Siemens scanner had up to 7 lp/cm for the head FOV and H40 kernel and up to 5 lp/cm at body FOV for the B40f kernel. The Standard kernel in the AIRO system was evaluated to have 3 lp/cm and 4 lp/cm for the body and head FOVs respectively. NNPS of the AIRO shows low frequency noise due to ring-like artifacts which may be caused by detector calibration or lack of artifact reducing image post-processing. Due to a higher dose in terms of mGy/mAs at both head and body FOV, the contrast to noise ratio is higher in the AIRO system than in the Siemens scanner. However detectability of the low contrast objects is poorer in the AIRO due to the presence of ring artifacts in the location of the targets.