A method to generate equivalent energy spectra and filtration models based on measurement for multidetector CT Monte Carlo dosimetry simulations

Accuracy of patient-specific CT organ doses from Monte Carlo simulations: influence of CT-based voxel models

Monte Carlo simulations using patient CT images as input are the gold standard to perform patient-specific dosimetry. However, in standard clinical practice patient's CT images are limited to the reconstructed CT scan range. In this study, organ dose calculations were performed with ImpactMC for chest and cardiac CT using whole-body and anatomy-specific voxel models to estimate the accuracy of CT organ doses based on the latter model. When the 3D patient model is limited to the CT scan range, CT organ doses from Monte Carlo simulations are the most accurate for organs entirely in the field of view. For these organs only the radiation dose related to scatter from the rest of the body is not incorporated. For organs lying partially outside the field of view organ doses are overestimated by not accounting for the non-irradiated tissue mass. This overestimation depends strongly on the amount of the organ volume located outside the field of view. To get a more accurate estimation of the radiation dose to these organs, the ICRP reference organ masses and densities could form a solution. Except for the breast, good agreement in dose was found for most organs. Voxel models generated from clinical CT examinations do not include the overscan in the z-direction. The availability of whole-body voxel models allowed to study this influence as well. As expected, overscan induces slightly higher organ doses.

PATIENT-SPECIFIC ORGAN DOSE EVALUATION BASED ON MONTE CARLO SIMULATION AND DOSE METRICS IN PAEDIATRIC CHEST-ABDOMEN-PELVIS CT EXAMINATIONS

This study aimed to determine organ doses based on Monte Carlo (MC) simulations for individual paediatric patients undergoing chest-abdomen-pelvis computed tomography (CT) examinations and to evaluate correlations of organ doses with dose metrics. MC simulations were performed by inputting detailed descriptions of a CT scanner, scanning parameters and CT images of 51 paediatric patients aged from 0 to 10 years into the simulation software. Organ doses for six radiosensitive organs were determined from dose distribution images obtained as the simulation results. The correlations of organ doses with dose metrics such as volume CT dose index, size-specific dose estimates (SSDEs), and organ-specific SSDEs were evaluated from the corresponding coefficients of determination. Organ doses for ages of 0-1 and 1-5 years were 22%-32% lower than those for ages of 5-10 years. Organ doses exhibited higher linear correlations with SSDEs and organ-specific SSDEs, and can be easily estimated using the linear regression.

Evaluating Size-Specific Dose Estimate (SSDE) as an estimate of organ doses from routine CT exams derived from Monte Carlo simulations

PURPOSE

Size-specific dose estimate (SSDE) is a metric that adjusts CTDI_vol to account for patient size. While not intended to be an estimate of organ dose, AAPM Report 204 notes the difference between the patient organ dose and SSDE is expected to be 10-20%. The purpose of this work was therefore to evaluate SSDE against estimates of organ dose obtained using Monte Carlo (MC) simulation techniques applied to routine exams across a wide range of patient sizes.

MATERIALS AND METHODS

Size-specific dose estimate was evaluated with respect to organ dose based on three routine protocols taken from Siemens scanners: (a) brain parenchyma dose in routine head exams, (b) lung and breast dose in routine chest exams, and (c) liver, kidney, and spleen dose in routine abdomen/pelvis exams. For each exam, voxelized phantom models were created from existing models or derived from clinical patient scans. For routine head exams, 15 patient models were used which consisted of 10 GSF/ICRP voxelized phantom models and five pediatric voxelized patient models created from CT image data. For all exams, the size metric used was water equivalent diameter (D_w ). For the routine chest exams, data from 161 patients were collected with a D_w range of ~16-44 cm. For the routine abdomen/pelvis exams, data from 107 patients were collected with a range of D_w from ~16 to 44 cm. Image data from these patients were segmented to generate voxelized patient models. For routine head exams, fixed tube current (FTC) was used while tube current modulation (TCM) data for body exams were extracted from raw projection data. The voxelized patient models and tube current information were used in detailed MC simulations for organ dose estimation. Organ doses from MC simulation were normalized by CTDI_vol and parameterized as a function of D_w . For each patient scan, the SSDE was obtained using D_w and CTDI_vol values of each scan, according to AAPM Report 220 for body scans and Report 293 for head scans. For each protocol and each patient, normalized organ doses were compared with SSDE. A one-sided tolerance limit covering 95% (P = 0.95) of the population with 95% confidence (α = 0.05) was used to assess the upper tolerance limit (T_U ) between SSDE and normalized organ dose.

RESULTS

For head exams, the T_U between SSDE and brain parenchyma dose was observed to be 12.5%. For routine chest exams, the T_U between SSDE and lung and breast dose was observed to be 35.6% and 68.3%, respectively. For routine abdomen/pelvis exams, the T_U between SSDE and liver, spleen, and kidney dose was observed to be 30.7%, 33.2%, and 33.0%, respectively.

CONCLUSIONS

The T_U of 20% between SSDE and organ dose was found to be insufficient to cover 95% of the sampled population with 95% confidence for all of the organs and protocols investigated, except for brain parenchyma dose. For the routine body exams, excluding the breasts, a wider threshold difference of ~30-36% would be needed. These results are, however, specific to Siemens scanners.

Estimation of patient skin dose in fluoroscopy: summary of a joint report by AAPM TG357 and EFOMP

BACKGROUND

Physicians use fixed C-arm fluoroscopy equipment with many interventional radiological and cardiological procedures. The associated effective dose to a patient is generally considered low risk, as the benefit-risk ratio is almost certainly highly favorable. However, X-ray-induced skin injuries may occur due to high absorbed patient skin doses from complex fluoroscopically guided interventions (FGI). Suitable action levels for patient-specific follow-up could improve the clinical practice. There is a need for a refined metric regarding follow-up of X-ray-induced patient injuries and the knowledge gap regarding skin dose-related patient information from fluoroscopy devices must be filled. The most useful metric to indicate a risk of erythema, epilation or greater skin injury that also includes actionable information is the peak skin dose, that is, the largest dose to a region of skin.

MATERIALS AND METHODS

The report is based on a comprehensive review of best practices and methods to estimate peak skin dose found in the scientific literature and situates the importance of the Digital Imaging and Communication in Medicine (DICOM) standard detailing pertinent information contained in the Radiation Dose Structured Report (RDSR) and DICOM image headers for FGI devices. Furthermore, the expertise of the task group members and consultants have been used to bridge and discuss different methods and associated available DICOM information for peak skin dose estimation.

RESULTS

The report contributes an extensive summary and discussion of the current state of the art in estimating peak skin dose with FGI procedures with regard to methodology and DICOM information. Improvements in skin dose estimation efforts with more refined DICOM information are suggested and discussed.

CONCLUSIONS

The endeavor of skin dose estimation is greatly aided by the continuing efforts of the scientific medical physics community, the numerous technology enhancements, the dose-controlling features provided by the FGI device manufacturers, and the emergence and greater availability of the DICOM RDSR. Refined and new dosimetry systems continue to evolve and form the infrastructure for further improvements in accuracy. Dose-related content and information systems capable of handling big data are emerging for patient dose monitoring and quality assurance tools for large-scale multihospital enterprises.

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.

A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PURPOSE

The conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.

METHODS

First, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8-40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba's CTDI_w methods, respectively.

RESULTS

The relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba's CTDI_w methods, respectively.

CONCLUSIONS

The devised f(0)^CB_w value was calculated by averaging four "point doses" at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.

EVALUATION OF ORGAN DOSE BASED ON PATIENT DATA IN EAST ASIAN DESCENT PEDIATRIC HEAD CT EXAMINATIONS

The purpose of this study was to investigate other indices estimating absorbed dose for eye lens and brain, using clinical images of East Asian pediatric patients. We simulated head computed tomography (CT ) examinations in 104 pediatric patients. Effective diameter (deff) and water equivalent diameter (dw) were measured on clinical images. Various size metrics and age were compared with absorbed dose normalised by CTDIvol (nD). The nD was estimated for eye and brain. The nD tended to decrease with advancing age. R2 between age and nD were 0.38 and 0.31 for eye and brain, respectively. Increasing head diameters decreased each nD. R2 between deff and dw, and nD were 0.20-0.24 and 0.51-0.53 for eye and brain, respectively. Head sizes allowed us to estimate absorbed dose in brain CT on East Asian pediatric patients. Scanning parameters for pediatric head CT may need to be based on individual patient information.

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.

COMPARISON OF EXPERIMENTAL AND NUMERICAL METHODS OF PATIENT DOSE ESTIMATIONS IN CT USING ANTHROPOMORPHIC MODELS

The common methods for patient dose estimations in computed tomography (CT) are thermoluminescence dosemeter (TLD) measurements or the usage of software packages based on Monte Carlo simulations like CT-Expo or the newer CTVoxDos, which uses the ICRP Reference Adult Male (ICRP 110). Organ (OD) and effective doses of a CT protocol of the upper abdomen are compared. Compared to CTVoxDos, ODs inferred by TLD measurement using an anthropomorphic phantom differ by $\mathbf{(19\pm 16)\,\%}$ inside the primary radiation field, $\mathbf{(14\pm 2)\,\%}$ for partially primary irradiated organs and $\mathbf{(34\pm 38)\,\%}$ in the scattered radiation field. ODs estimated by CT-Expo show a mean deviation of $\mathbf{(16\pm 9)\,\%}$ (primary irradiated) and $\mathbf{(28\pm 31)\,\%}$ (scatter irradiated) from ODs estimated by CTVoxDos.

Peak skin and eye lens radiation dose from brain perfusion CT: CTDIvol and Monte Carlo based estimations

PURPOSE

To quantify the eye lens, peak skin and brain doses associated with head CT perfusion exam by means of thermoluminescent dosimeters (TLDs) measurements in a cadaver and compare them to Monte Carlo (MC) dose estimations as well as to the CTDI_vol.

METHOD

18 TLDs were inserted in the brain, skin, and eye lenses of a female cadaver head, who underwent a CT brain perfusion scan using a Siemens Definition Flash. The table-toggling protocol used 80 kVp, 200 mAs, 32 × 1.2 mm collimation and 30 sequences. From the CT images, a voxel model was created. Doses were calculated with a MC framework (EGSnrc) and compared to TLD measurements. TLD measurements were also compared to the displayed CTDI_vol.

RESULTS

The average measured doses were: 185 mGy for the eyes lenses, 107 mGy for the skin, 172 mGy for the brain and 273 mGy for the peak skin. The reported CTDI_vol of 259 mGy overestimated the averaged organ doses but not the peak skin dose. MC estimated organ doses were 147 mGy for the eyes (average), 104 mGy for the skin and 178 mGy for the brain (-20 %, -3% and 4% difference respect to the TLDs measurements, respectively).

CONCLUSIONS

CTDI_vol remains a conservative metric for average brain, skin and eyes lenses doses. For accurate eye lens and skin dose estimates MC simulations can be used. CTDI_vol should be used with caution as it was of the same order of magnitude as the peak skin dose for this protocol and this particular CT scanner.

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.

Size-specific dose estimates for various weighting factors of CTDI equation

Size-specific dose estimate (SSDE) was proposed by the American Association of Physicists in Medicine (AAPM) Task Group 204 to consider the effect of patient size in the x-ray CT dose estimation. Size correction factors to calculate SSDE were derived based on the conventional weighted CT dose index (CTDIw) equation. This study aims to investigate the influence of Bakalyar's and the authors' own CTDIw equations on the size correction factors described by the AAPM Task Group 204, using Monte Carlo simulations. The simulations were performed by modeling four types of x-ray CT scanner designs, to compute the dose values in water for cylindrical phantoms with 8-40 cm diameters. CTDI100 method and the AAPM Task Group 111's proposed method were employed as the CT dosimetry models. Size correction factors were obtained for the computed dose values of various phantom diameters for the conventional, Bakalyar's, and the authors' weighting factors. Maximum difference between the size correction factors for the Bakalyar's weighting factor and those of the AAPM Task Group 204 was 27% for a phantom diameter of 11.2 cm. On the other hand, the size correction factors calculated for the authors' weighting factor were in good agreement with those from the AAPM Task Group 204 report with a maximum difference of 17%. The results indicate that the SSDE values obtained with the authors' weighting factor can be evaluated by using the size correction factors reported by the AAPM Task Group 204, which is currently accepted as a standard.

A Scalable Database of Organ Doses for Common Diagnostic Fluoroscopy Procedures of Children: Procedures of Historical Practice for Use in Radiation Epidemiology Studies

Assessment of health effects from low-dose radiation exposures in patients undergoing diagnostic imaging is an active area of research. High-quality dosimetry information pertaining to these medical exposures is generally not readily available to clinicians or epidemiologists studying radiation-related health risks. The purpose of this study was to provide methods for organ dose estimation in pediatric patients undergoing four common diagnostic fluoroscopy procedures: the upper gastrointestinal (UGI) series, the lower gastrointestinal (LGI) series, the voiding cystourethrogram (VCUG) and the modified barium swallow (MBS). Abstracted X-ray film data and physician interviews were combined to generate procedure outlines detailing X-ray beam projections, imaged anatomy, length of X-ray exposure, and presence and amount of contrast within imaged anatomy. Monte Carlo radiation transport simulations were completed for each of the four diagnostic fluoroscopy procedures across the 162-member (87 males and 75 females) University of Florida/National Cancer Institute pediatric phantom library, which covers variations in both subject height and weight. Absorbed doses to 28 organs, including the active marrow and bone endosteum, were assigned for all 162 phantoms by procedure. Additionally, we provide dose coefficients (DCs) in a series of supplementary tables. The DCs give organ doses normalized to procedure-specific dose metrics, including: air kerma-area product (µGy/mGy · cm2), air kerma at the reference point (µGy/µGy), number of spot films (SF) (µGy/number of SFs) and total fluoroscopy time (µGy/s). Organs accumulating the highest absorbed doses per procedure were as follows: kidneys between 0.9-25.4 mGy, 1.1-16.6 mGy and 1.1-9.7 mGy for the UGI, LGI and VCUG procedures, respectively, and salivary glands between 0.2-3.7 mGy for the MBS procedure. Average values of detriment-weighted dose, a phantom-specific surrogate for the effective dose based on ICRP Publication 103 tissue-weighting factors, were 0.98 mSv, 1.16 mSv, 0.83 mSv and 0.15 mSv for the UGI, LGI, VCUG and MBS procedures, respectively. Scalable database of organ dose coefficients by patient sex, height and weight, and by procedure exposure time, reference point air kerma, kerma-area product or number of spot films, allows clinicians and researchers to compute organ absorbed doses based on their institution-specific and patient-specific dose metrics. In addition to informing on patient dosimetry, this work has the potential to facilitate exposure assessments in epidemiological studies designed to investigate radiation-related risks.

Size-specific dose estimations for pediatric chest, abdomen/pelvis and head CT scans with the use of GATE

PURPOSE

The purpose of this study is to create an organ dose database for pediatric individuals undergoing chest, abdomen/pelvis, and head computed tomography (CT) examinations, and to report the differences in absorbed organ doses, when anatomical differences exist for pediatric patients.

METHODS

The GATE Monte Carlo (MC) toolkit was used to model the GE BrightSpeed Elite CT model. The simulated scanner model was validated with the standard Computed Tomography Dose Index (CTDI) head phantom. Twelve computational models (2.1-14 years old) were used. First, contributions to effective dose and absorbed doses per CTDI_vol and per 100 mAs were estimated for all organs. Then, doses per CTDI_vol were correlated with patient model weight for the organs inside the scan range for chest and abdomen/pelvis protocols. Finally, effective doses per dose-length product (DLP) were estimated and compared with the conventional conversion k-factors.

RESULTS

The system was validated against experimental CTDI_w measurements. The doses per CTDI_vol and per 100 mAs for selected organs were estimated. The magnitude of the dependency between the dose and the anatomical characteristics was calculated with the coefficient of determination at 0.5-0.7 for the internal scan organs for chest and abdomen/pelvis protocols. Finally, effective doses per DLP were compared with already published data, showing discrepancies between 13 and 29% and were correlated strongly with the total weight (R² > 0.8) for the chest and abdomen protocols.

CONCLUSIONS

Big differences in absorbed doses are reported even for patients of similar age or same gender, when anatomical differences exist on internal organs of the body.

A scalable database of organ doses for common diagnostic fluoroscopy examinations of children: procedures of current practice at the University of Florida

Of all the medical imaging modalities that utilize ionizing radiation, fluoroscopy proves to be the most difficult to assess values of patient organ dose owing to the dynamic and patient-specific nature of the irradiation geometry and its associated x-ray beam characteristics. With the introduction of the radiation dose structured report (RDSR) in the mid-2000s, however, computational tools have been developed to extract patient and procedure-specific data for each irradiation event of the study, and when coupled to a computational phantom of the patient, values of skin and internal organ dose may be assessed. Unfortunately, many legacy and even current diagnostic fluoroscopy units do not have RDSR reporting capabilities, thus limiting these dosimetry reporting advances. Nevertheless, knowledge of patient organ doses for patient care, as well as for radiation epidemiology studies, remains a research and regulatory priority. In this study, we created procedural outlines which document all radiation exposure information required for organ dose assessment, akin to a reference RDSR, for six common diagnostic fluoroscopy procedures performed at the University of Florida (UF) Shands Pediatric Hospital. These procedures include the voiding cystourethrogram, the gastrostomy-tube placement, the lower gastrointestinal study, the rehabilitation swallow, the upper gastrointestinal study, and the upper gastrointestinal study with follow through. These procedural outlines were used to develop an extensive database of organ doses for the 162-member UF/NCI (National Cancer Institute) library of pediatric hybrid phantoms, with each member varying combinations of sex, height, and weight. The organ dose assessment accounts for the varying x-ray fields, fluoroscopy time, relative concentration of x-ray contrast in the organs, and changes in the fluoroscope output due to patient size. Furthermore, we are also reporting organ doses normalized to total fluoroscopy time, reference point air kerma, and kerma-area product, effectively providing procedure dose coefficients. The extensive organ dose library produced in this study may be used prospectively for patient organ dose reporting or retrospectively in epidemiological studies of radiation-associated health risks.

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.

CT beam dosimetric characterization procedure for personalized dosimetry

Personalized dosimetry in computed tomography (CT) can be realized by a full Monte Carlo (MC) simulation of the scan procedure. Essential input data needed for the simulation are appropriate CT x-ray source models and a model of the patient's body which is based on the CT image. The purpose of this work is to develop comprehensive procedures for the determination of CT x-ray source models and their verification by comparison of calculated and measured dose distributions in physical phantoms. Mobile equipment together with customized software was developed and used for non-invasive determination of equivalent source models of CT scanners under clinical conditions. Standard and physical anthropomorphic CT dose phantoms equipped with real-time CT dose probes at five representative positions were scanned. The accumulated dose was measured during the scan at the five positions. ImpactMC, an MC-based CT dose software program, was used to simulate the scan. The necessary inputs were obtained from the scan parameters, from the equivalent source models and from the material-segmented CT images of the phantoms. 3D dose distributions in the phantoms were simulated and the dose values calculated at the five positions inside the phantom were compared to measured dose values. Initial results were obtained by means of a General Electric Optima CT 660 and a Toshiba (Canon) Aquilion ONE. In general, the measured and calculated dose values were within relative uncertainties that had been estimated to be less than 10%. The procedures developed were found to be viable and rapid. The procedures are applicable to any scanner type under clinical conditions without making use of the service mode with stationary x-ray tube position. Results show that the procedures are well suited for determining and verifying the equivalent source models needed for personalized CT dosimetry based on post-scan MC calculations.

Development and validation of an open source Monte Carlo dosimetry model for wide-beam CT scanners using Fluka

PURPOSE

Development and validation of an open source Fluka-based Monte Carlo source model for diagnostic patient dose calculations.

METHODS

A framework to simulate a computed tomography (CT) scanner using Fluka Monte Carlo particle transport code was developed. The General Electric (GE) Revolution scanner with the large body filter and 120 kV tube potential was characterized using measurements. The model was validated on benchmark CT test problems and on dose measurements in computed tomography dose index (CTDI) and anthropomorphic phantoms. Axial and helical operation modes with provision for tube current modulation (TCM) were implemented. The particle simulations in Fluka were accelerated by executing them on a high-performance computing cluster.

RESULTS

The simulation results agreed to better than an average of 4% of the reference simulation results from the AAPM Report 195 test scenarios, namely: better than 2% for both test problems in case 4 using the PMMA phantom, and better than 5% of the reference result for 14 of 17 organs in case 5, and within 10% for the three remaining organs. The Fluka simulation results agreed to better than 2% of the air kerma measured in-air at isocenter of the GE Revolution scanner. The simulated air kerma in the center of the CTDI phantom overestimated the measurement by 7.5% and a correction factor was introduced to account for this. The simulated mean absorbed doses for a chest scan of the pediatric anthropomorphic phantom was completed in ~9 min and agreed to within the 95% CI for bone, soft tissue, and lung measurements made using MOSFET detectors for fixed current axial and helical scans as well as helical scan with TCM.

CONCLUSION

A Fluka-based Monte Carlo simulation model of axial and helical acquisition techniques using a wide-beam collimation CT scanner demonstrated good agreement between measured and simulated results for both fixed current and TCM in complex and simple geometries. Code and dataset will be made available at https://github.com/chezhia/FLUKA_CT.

Calibration strategies for use of the nanoDot OSLD in CT applications

Aluminum oxide based optically stimulated luminescent dosimeters (OSLD) have been recognized as a useful dosimeter for measuring CT dose, particularly for patient dose measurements. Despite the increasing use of this dosimeter, appropriate dosimeter calibration techniques have not been established in the literature; while the manufacturer offers a calibration procedure, it is known to have relatively large uncertainties. The purpose of this work was to evaluate two clinical approaches for calibrating these dosimeters for CT applications, and to determine the uncertainty associated with measurements using these techniques. Three unique calibration procedures were used to calculate dose for a range of CT conditions using a commercially available OSLD and reader. The three calibration procedures included calibration (a) using the vendor-provided method, (b) relative to a 120 kVp CT spectrum in air, and (c) relative to a megavoltage beam (implemented with 60 Co). The dose measured using each of these approaches was compared to dose measured using a calibrated farmer-type ion chamber. Finally, the uncertainty in the dose measured using each approach was determined. For the CT and megavoltage calibration methods, the dose measured using the OSLD nanoDot was within 5% of the dose measured using an ion chamber for a wide range of different CT scan parameters (80-140 kVp, and with measurements at a range of positions). When calibrated using the vendor-recommended protocol, the OSLD measured doses were on average 15.5% lower than ion chamber doses. Two clinical calibration techniques have been evaluated and are presented in this work as alternatives to the vendor-provided calibration approach. These techniques provide high precision for OSLD-based measurements in a CT environment.

Estimating a size-specific dose for helical head CT examinations using Monte Carlo simulation methods

PURPOSE

Size-specific dose estimates (SSDE) conversion factors have been determined by AAPM Report 204 to adjust CTDI_vol to account for patient size but were limited to body CT examinations. The purpose of this work was to determine conversion factors that could be used for an SSDE for helical, head CT examinations for patients of different sizes.

METHODS

Validated Monte Carlo (MC) simulation methods were used to estimate dose to the center of the scan volume from a routine, helical head examination for a group of patient models representing a range of ages and sizes. Ten GSF/ICRP voxelized phantom models and five pediatric voxelized patient models created from CT image data were used in this study. CT scans were simulated using a Siemens multidetector row CT equivalent source model. Scan parameters were taken from the AAPM Routine Head protocols for a fixed tube current (FTC), helical protocol, and scan lengths were adapted to the anatomy of each patient model. MC simulations were performed using mesh tallies to produce voxelized dose distributions for the entire scan volume of each model. Three tally regions were investigated: (1) a small 0.6 cc volume at the center of the scan volume, (2) 0.8-1.0 cm axial slab at the center of the scan volume, and (3) the entire scan volume. Mean dose to brain parenchyma for all three regions was calculated. Mean bone dose and a mass-weighted average dose, consisting of brain parenchyma and bone, were also calculated for the slab in the central plane and the entire scan volume. All dose measures were then normalized by CTDI_vol for the 16 cm phantom (CTDI_vol,16 ). Conversion factors were determined by calculating the relationship between normalized doses and water equivalent diameter (D_w ).

RESULTS

CTDI_vol,16 -normalized mean brain parenchyma dose values within the 0.6 cc volume, 0.8-1.0 cm central axial slab, and the entire scan volume, when parameterized by D_w , had an exponential relationship with a coefficient of determination (R² ) of 0.86, 0.84, and 0.88, respectively. There was no statistically significant difference between the conversion factors resulting from these three different tally regions. Exponential relationships between CTDI_vol,16 -normalized mean bone doses had R² values of 0.83 and 0.87 for the central slab and for the entire scan volume, respectively. CTDI_vol,16 -normalized mass-weighted average doses had R² values of 0.39 and 0.51 for the central slab and for the entire scan volume, respectively.

CONCLUSIONS

Conversion factors that describe the exponential relationship between CTDI_vol,16 -normalized mean brain dose and a size metric (D_w ) for helical head CT examinations have been reported for two different interpretations of the center of the scan volume. These dose descriptors have been extended to describe the dose to bone in the center of the scan volume as well as a mass-weighted average dose to brain and bone. These may be used, when combined with other efforts, to develop an SSDE dose coefficients for routine, helical head CT examinations.

VALIDATION OF MONTE CARLO DOSE CALCULATION FOR PAEDIATRIC CT EXAMINATIONS USING TUBE CURRENT MODULATION BASED ON IN-PHANTOM DOSIMETRY

The aim of this study is to estimate tube current modulation (TCM) profiles in paediatric computed tomography (CT) examinations with a TCM scheme (Volume-EC) and evaluate the estimation accuracy of TCM profiles. Another aim is to validate organ doses calculated using Monte Carlo-based CT dosimetry software and estimated TCM profiles by comparing them with those measured using 5-year-old and 10-year-old anthropomorphic phantoms and radio-photoluminescence glass dosemeters. Dose calculations were performed by inputting detailed descriptions of a CT scanner, scan parameters and CT images of the phantoms into the software. Organ doses were evaluated from the calculated dose distribution images. Average relative differences (RDs) between the estimated and actual TCM profiles ranged from -3.6 to 5.6%. RDs between the calculated and measured organ doses ranged from -4.2 to 13.0% and -18.1 to 4.9% for 5-year-old and 10-year-old phantoms, respectively. These results validate dose calculations for paediatric CT scans using TCM.

Estimating lung, breast, and effective dose from low-dose lung cancer screening CT exams with tube current modulation across a range of patient sizes

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Two-step validation of a Monte Carlo dosimetry framework for general radiology

The Monte Carlo technique is considered gold standard when it comes to patient-specific dosimetry. Any newly developed Monte Carlo simulation framework, however, has to be carefully calibrated and validated prior to its use. For many researchers this is a tedious work. We propose a two-step validation procedure for our newly built Monte Carlo framework and provide all input data to make it feasible for future related application by the wider community. The validation was at first performed by benchmarking against simulation data available in literature. The American Association of Physicists in Medicine (AAPM) report of task group 195 (case 2) was considered most appropriate for our application. Secondly, the framework was calibrated and validated against experimental measurements for trunk X-ray imaging protocols using a water phantom. The dose results obtained from all simulations and measurements were compared. Our Monte Carlo framework proved to agree with literature data, by showing a maximal difference below 4% to the AAPM report. The mean difference with the water phantom measurements was around 7%. The statistical uncertainty for clinical applications of the dosimetry model is expected to be within 10%. This makes it reliable for clinical dose calculations in general radiology. Input data and the described procedure allow for the validation of other Monte Carlo frameworks.

Consideration of the usefulness of a size-specific dose estimate in pediatric CT examination

Computed tomography (CT) has recently been utilized in various medical settings, and technological advances have resulted in its widespread use. However, medical radiation exposure associated with CT scans accounts for the largest share of examinations using radiation; thus, it is important to understand the organ dose and effective dose in detail. The CT dose index and dose-length product are used to evaluate the organ dose. However, evaluations using these indicators fail to consider the age and body type of patients. In this study, we evaluated the effective dose based on the CT examination data of 753 patients examined at our hospital using the size-specific dose estimate (SSDE) method, which can calculate the exposure dose with consideration of the physique of a patient. The results showed a large correlation between the SSDE conversion factor and physique, with a larger exposure dose in patients with a small physique when a single scan is considered. Especially for children, the SSDE conversion factor was found to be 2 or more. In addition, the patient exposed to the largest dose in this study was a 10-year-old, who received 40.4 mSv (five series/examination). In the future, for estimating exposure using the SSDE method and in cohort studies, the diagnostic reference level of SSDE should be determined and a low-exposure imaging protocol should be developed to predict the risk of CT exposure and to maintain the quality of diagnosis with better radiation protection of patients.

Organ doses in pediatric patients undergoing cardiac-centered fluoroscopically guided interventions: Comparison of three methods for computational phantom alignment

PURPOSE

To assess various computational phantom alignment techniques within Monte Carlo radiation transport models of pediatric fluoroscopically guided cardiac interventional studies.

METHODS

Logfiles, including all procedure radiation and machine data, were extracted from a Toshiba infinix-I unit in the University of Florida Pediatric Catheterization Laboratory for a cohort of 10 patients. Two different alignment methods were then tested against a ground truth standard based upon identification of a unique anatomic reference point within images co-registered to specific irradiation events within each procedure. The first alignment method required measurement of the distance from the edge of the exam table to the top of the patient's head (table alignment method). The second alignment method fixed the anatomic reference point to be the geometric center of the heart muscle, as all 10 studies were cardiac in nature. Monte Carlo radiation transport simulations were performed for each patient and intervention using morphometry-matched hybrid computational phantoms for the reference and two tested alignment methods. For each combination, absorbed doses were computed for 28 organs and root mean square organ doses were assessed and compared across the alignment methods.

RESULTS

The percent error in root mean square organ dose ranged from -57% to +41% for the table alignment method, and from -27% to +22% for the heart geometric centroid alignment method. Absorbed doses to specific organs, such as the heart and lungs, demonstrated higher accuracy in the heart geometric centroid alignment method, with average percent errors of 10% and 1.4%, respectively, compared to average percent errors of -32% and 24%, respectively, using the table alignment method.

CONCLUSIONS

Of the two phantom alignment methods investigated in this study, the use of an anatomical reference point - in this case the geometric centroid of the heart - provided a reliable method for radiation transport simulations of organ dose in pediatric interventional cardiac studies. This alignment method provides the added benefit of requiring no physician input, making retrospective calculations possible. Moving forward, additional anatomical reference methods can be tested to assess the reliability of anatomical reference points beyond cardiac centered procedures.

Physical validation of UF-RIPSA: A rapid in-clinic peak skin dose mapping algorithm for fluoroscopically guided interventions

Purpose

The purpose of this study was to experimentally validate UF‐RIPSA, a rapid in‐clinic peak skin dose mapping algorithm developed at the University of Florida using optically stimulated luminescent dosimeters (OSLDs) and tissue‐equivalent phantoms.

Methods

The OSLDs used in this study were InLight^TM Nanodot dosimeters by Landauer, Inc. The OSLDs were exposed to nine different beam qualities while either free‐in‐air or on the surface of a tissue equivalent phantom. The irradiation of the OSLDs was then modeled using Monte Carlo techniques to derive correction factors between free‐in‐air exposures and more complex irradiation geometries. A grid of OSLDs on the surface of a tissue equivalent phantom was irradiated with two fluoroscopic x ray fields generated by the Siemens Artis zee bi‐plane fluoroscopic unit. The location of each OSLD within the grid was noted and its dose reading compared with UF‐RIPSA results.

Results

With the use of Monte Carlo correction factors, the OSLD's response under complex irradiation geometries can be predicted from its free‐in‐air response. The predicted values had a percent error of −8.7% to +3.2% with a predicted value that was on average 5% below the measured value. Agreement within 9% was observed between the values of the OSLDs and RIPSA when irradiated directly on the phantom and within 14% when the beam first traverses the tabletop and pad.

Conclusions

The UF‐RIPSA only computes dose values to areas of irradiated skin determined to be directly within the x ray field since the algorithm is based upon ray tracing of the reported reference air kerma value, with subsequent corrections for air‐to‐tissue dose conversion, x ray backscatter, and table/pad attenuation. The UF‐RIPSA algorithm thus does not include the dose contribution of scatter radiation from adjacent fields. Despite this limitation, UF‐RIPSA is shown to be fairly robust when computing skin dose to patients undergoing fluoroscopically guided interventions.

Assessment of PCXMC for patients with different body size in chest and abdominal x ray examinations: a Monte Carlo simulation study

A PC Program for x ray Monte Carlo (PCXMC) has been used to calculate organ doses in patient dosimetry and for the exposure assessment in epidemiological studies of radiogenic health related risks. This study compared the dosimetry from using the built-in stylized phantoms in the PCXMC to that of a newer hybrid phantom library with improved anatomical realism. We simulated chest and abdominal x ray projections for 146 unique body size computational phantoms, 77 males and 69 females, with different combinations of height (125-180 cm) and weight (20-140 kg) using the built-in stylized phantoms in the PCXMC version 2.0.1.4 and the hybrid phantom library using the Monte Carlo N-particle eXtended transport code 2.7 (MCNPX). Unfortunately, it was not possible to incorporate the hybrid phantom library into the PCXMC. We compared 14 organ doses, including dose to the active bone marrow, to evaluate differences between the built-in stylized phantoms in the PCXMC and the hybrid phantoms (Cristy and Eckerman 1987 Technical Report ORNL/TM-8381/V1, Oak Ridge National Laboratory, Eckerman and Ryman 1993 Technical Report 12 Oak Ridge, TN, Geyer et al 2014 Phys. Med. Biol. 59 5225-42). On average, organ doses calculated using the built-in stylized phantoms in the PCXMC were greater when compared to the hybrid phantoms. This is most prominent in AP abdominal exams by an average factor of 2.4-, 2.8-, and 2.8-fold for the 10-year-old, 15-year-old, and adult phantoms, respectively. For chest exams, organ doses are greater by an average factor of 1.1-, 1.4-, and 1.2-fold for the 10-year-old, 15-year-old, and adult phantoms, respectively. The PCXMX, due to its ease of use, is often selected to support dosimetry in epidemiological studies; however, it uses simplified models of the human anatomy that fail to account for variations in body morphometry for increasing weight. For epidemiological studies that use PCXMC dosimetry, associations between radiation-related disease risks and organ doses may be underestimated, and to a greater degree in pediatric, especially obese pediatric, compared to adult patients.

Evaluation of the UF/NCI hybrid computational phantoms for use in organ dosimetry of pediatric patients undergoing fluoroscopically guided cardiac procedures

Epidemiologic data demonstrate that pediatric patients face a higher relative risk of radiation induced cancers than their adult counterparts at equivalent exposures. Infants and children with congenital heart defects are a critical patient population exposed to ionizing radiation during life-saving procedures. These patients will likely incur numerous procedures throughout their lifespan, each time increasing their cumulative radiation absorbed dose. As continued improvements in long-term prognosis of congenital heart defect patients is achieved, a better understanding of organ radiation dose following treatment becomes increasingly vital. Dosimetry of these patients can be accomplished using Monte Carlo radiation transport simulations, coupled with modern anatomical patient models. The aim of this study was to evaluate the performance of the University of Florida/National Cancer Institute (UF/NCI) pediatric hybrid computational phantom library for organ dose assessment of patients that have undergone fluoroscopically guided cardiac catheterizations. In this study, two types of simulations were modeled. A dose assessment was performed on 29 patient-specific voxel phantoms (taken as representing the patient's true anatomy), height/weight-matched hybrid library phantoms, and age-matched reference phantoms. Two exposure studies were conducted for each phantom type. First, a parametric study was constructed by the attending pediatric interventional cardiologist at the University of Florida to model the range of parameters seen clinically. Second, four clinical cardiac procedures were simulated based upon internal logfiles captured by a Toshiba Infinix-i Cardiac Bi-Plane fluoroscopic unit. Performance of the phantom library was quantified by computing both the percent difference in individual organ doses, as well as the organ dose root mean square values for overall phantom assessment between the matched phantoms (UF/NCI library or reference) and the patient-specific phantoms. The UF/NCI hybrid phantoms performed at percent differences of between 15% and 30% for the parametric set of irradiation events. Among internal logfile reconstructed procedures, the UF/NCI hybrid phantoms performed with RMS organ dose values between 7% and 29%. Percent improvement in organ dosimetry via the use of hybrid library phantoms over the reference phantoms ranged from 6.6% to 93%. The use of a hybrid phantom library, Monte Carlo radiation transport methods, and clinical information on irradiation events provide a means for tracking organ dose in these radiosensitive patients undergoing fluoroscopically guided cardiac procedures.

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).

Physical validation of a Monte Carlo-based, phantom-derived approach to computed tomography organ dosimetry under tube current modulation

PURPOSE

To physically validate the accuracy of a Monte Carlo-based, phantom-derived methodology for computed tomography (CT) dosimetry that utilizes organ doses from precomputed axial scans and that accounts for tube current modulation (TCM).

METHODS

The output of a Toshiba Aquilion ONE CT scanner was modeled, based on physical measurement, in the Monte Carlo radiation transport code MCNPX (v2.70). CT examinations were taken of two anthropomorphic phantoms representing pediatric and adult patients (15-yr-old female and adult male) at various energies, in which physical organ dose measurements were made using optically stimulated luminescence dosimeters (OSLDs). These exams (chest-abdomen-pelvis) were modeled using organ dose data obtained from the computationally equivalent phantom of each anthropomorphic phantom. TCM was accounted for by weighting all organ dose contributions by both the relative attenuation of the phantom and the image-derived mA value (from the DICOM header) at the same z-extent (cranial-caudal direction) of the axial dose data.

RESULTS

The root mean squares of percent difference in organ dose when comparing the physical organ dose measurements to the computational estimates were 21.2, 12.1, and 15.1% for the uniform (no attenuation weighting), weighted (computationally derived), and image-based methodologies, respectively.

CONCLUSIONS

Overall, these data suggest that the Monte Carlo-based dosimetry presented in this work is viable for CT dosimetry. Additionally, for CT exams with TCM, local attenuation weighting of organ dose contributions from precomputed axial dosimetry libraries increases organ dose accuracy.

Assessment of different patient-to-phantom matching criteria applied in Monte Carlo-based computed tomography dosimetry

PURPOSE

To quantify differences in computationally estimated computed tomography (CT) organ doses for patient-specific voxel phantoms to estimated organ doses in matched computational phantoms using different matching criteria.

MATERIALS AND METHODS

Fifty-two patient-specific computational voxel phantoms were created through CT image segmentation. In addition, each patient-specific phantom was matched to six computational phantoms of the same gender based, respectively, on age and gender (reference phantoms), height and weight, effective diameter (both central slice and exam range average), and water equivalent diameter (both central slice and exam range average). Each patient-specific phantom and matched computational phantom were then used to simulate six different torso examinations using a previously validated Monte Carlo CT dosimetry methodology that accounts for tube current modulation. Organ doses for each patient-specific phantom were then compared with the organ dose estimates of each of the matched phantoms.

RESULTS

Relative to the corresponding patient-specific phantoms, the root mean square of the difference in organ dose was 39.1%, 20.3%, 22.7%, 21.6%, 20.5%, and 17.6%, for reference, height and weight, effective diameter (central slice and scan average), and water equivalent diameter (central slice and scan average), respectively. The average magnitude of difference in organ dose was 24%, 14%, 16.9%, 16.2%, 14%, and 11.9%, respectively.

CONCLUSION

Overall, these data suggest that matching a patient to a computational phantom in a library is superior to matching to a reference phantom. Water equivalent diameter is the superior matching metric, but it is less feasible to implement in a clinical and retrospective setting. For these reasons, height-and-weight matching is an acceptable and reliable method for matching a patient to a member of a computational phantom library with regard to CT dosimetry.

A hybrid phantom system for patient skin and organ dosimetry in fluoroscopically guided interventions

PURPOSE

The purpose of this study was to investigate calibrations for improved estimates of skin dose and to develop software for computing absorbed organ doses for fluoroscopically guided interventions (FGIs) with the use of radiation dose structured reports (RDSR) and the UF/NCI family of hybrid computational phantoms.

METHODS AND MATERIALS

Institutional review board approval was obtained for this retrospective study in which ten RDSRs were selected for their high cumulative reference air kerma values. Skin doses were computed using the University of Florida's rapid in-clinic peak skin dose algorithm (or UF-RIPSA). Kerma-area product (KAP) meter calibrations and attenuation of the tabletop with pad were incorporated into the UF-RIPSA. To compute absorbed organ doses the RDSRs were coupled with software to develop Monte Carlo input decks for each irradiation event. The effects of spectrum matching were explored by modeling (a) a polychromatic x-ray energy beam made to match measured first half-value layers of aluminum, (b) an unmatched spectrum, (c) and a mono-energetic beam equivalent to the effective x-ray energy. The authors also considered the practicality of computing organ doses for each irradiation event within a RDSR.

RESULTS

The KAP meter is highly dependent on the quality of the x-ray spectra. Monte Carlo based attenuation coefficients for configurations in which the beam is transmitted through the tabletop with pad reduced the amount by which the software overestimated skin doses. For absorbed organ dose computations, the average ratios of computed organ doses for a non-fitted to fitted spectrum and effective energy to fitted spectrum were 0.45 and 0.03, respectively. Monte Carlo simulations on average took 38 min per patient. All in-field organ tallies converged with a relative error of less than 1% and out-of-field organs tallies within 10% relative error.

CONCLUSIONS

This work details changes to the UF-RIPSA software that include an expanded library of computational phantoms, attenuation coefficients for tabletop with pad, and calibration curves for the KAP meter. For the computation of absorbed organ dose, it is possible to model each irradiation event separately on a patient-dependent model that best morphometrically matches the patient, thus providing a full report of internal organ doses for FGI patients.

Estimating patient dose from CT exams that use automatic exposure control: Development and validation of methods to accurately estimate tube current values

PURPOSE

The vast majority of body CT exams are performed with automatic exposure control (AEC), which adapts the mean tube current to the patient size and modulates the tube current either angularly, longitudinally or both. However, most radiation dose estimation tools are based on fixed tube current scans. Accurate estimates of patient dose from AEC scans require knowledge of the tube current values, which is usually unavailable. The purpose of this work was to develop and validate methods to accurately estimate the tube current values prescribed by one manufacturer's AEC system to enable accurate estimates of patient dose.

METHODS

Methods were developed that took into account available patient attenuation information, user selected image quality reference parameters and x-ray system limits to estimate tube current values for patient scans. Methods consistent with AAPM Report 220 were developed that used patient attenuation data that were: (a) supplied by the manufacturer in the CT localizer radiograph and (b) based on a simulated CT localizer radiograph derived from image data. For comparison, actual tube current values were extracted from the projection data of each patient. Validation of each approach was based on data collected from 40 pediatric and adult patients who received clinically indicated chest (n = 20) and abdomen/pelvis (n = 20) scans on a 64 slice multidetector row CT (Sensation 64, Siemens Healthcare, Forchheim, Germany). For each patient dataset, the following were collected with Institutional Review Board (IRB) approval: (a) projection data containing actual tube current values at each projection view, (b) CT localizer radiograph (topogram) and (c) reconstructed image data. Tube current values were estimated based on the actual topogram (actual-topo) as well as the simulated topogram based on image data (sim-topo). Each of these was compared to the actual tube current values from the patient scan. In addition, to assess the accuracy of each method in estimating patient organ doses, Monte Carlo simulations were performed by creating voxelized models of each patient, identifying key organs and incorporating tube current values into the simulations to estimate dose to the lungs and breasts (females only) for chest scans and the liver, kidney, and spleen for abdomen/pelvis scans. Organ doses from simulations using the actual tube current values were compared to those using each of the estimated tube current values (actual-topo and sim-topo).

RESULTS

When compared to the actual tube current values, the average error for tube current values estimated from the actual topogram (actual-topo) and simulated topogram (sim-topo) was 3.9% and 5.8% respectively. For Monte Carlo simulations of chest CT exams using the actual tube current values and estimated tube current values (based on the actual-topo and sim-topo methods), the average differences for lung and breast doses ranged from 3.4% to 6.6%. For abdomen/pelvis exams, the average differences for liver, kidney, and spleen doses ranged from 4.2% to 5.3%.

CONCLUSIONS

Strong agreement between organ doses estimated using actual and estimated tube current values provides validation of both methods for estimating tube current values based on data provided in the topogram or simulated from image data.

Realistic phantoms to characterize dosimetry in pediatric CT

BACKGROUND

The estimation of organ doses and effective doses for children receiving CT examinations is of high interest. Newer, more realistic anthropomorphic body models can provide information on individual organ doses and improved estimates of effective dose.

MATERIALS AND METHODS

Previously developed body models representing 50th-percentile individuals at reference ages (newborn, 1, 5, 10 and 15 years) were modified to represent 10th, 25th, 75th and 90th height percentiles for both genders and an expanded range of ages (3, 8 and 13 years). We calculated doses for 80 pediatric reference phantoms from simulated chest-abdomen-pelvis exams on a model of a Philips Brilliance 64 CT scanner. Individual organ and effective doses were normalized to dose-length product (DLP) and fit as a function of body diameter.

RESULTS

We calculated organ and effective doses for 80 reference phantoms and plotted them against body diameter. The data were well fit with an exponential function. We found DLP-normalized organ dose to correlate strongly with body diameter (R2>0.95 for most organs). Similarly, we found a very strong correlation with body diameter for DLP-normalized effective dose (R2>0.99). Our results were compared to other studies and we found average agreement of approximately 10%.

CONCLUSION

We provide organ and effective doses for a total of 80 reference phantoms representing normal-stature children ranging in age and body size. This information will be valuable in replacing the types of vendor-reported doses available. These data will also permit the recording and tracking of individual patient doses. Moreover, this comprehensive dose database will facilitate patient matching and the ability to predict patient-individualized dose prior to examination.

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.

Organ Dose Evaluations Based on Monte Carlo Simulation for CT Examinations Using Tube Current Modulation

The aims of this study were to estimate tube current values for each X-ray projection angle used in adult chest computed tomography (CT) and abdomen-pelvis CT examinations with tube current modulation (TCM) and to validate organ doses determined using Monte Carlo (MC) simulations through comparisons with the doses measured using in-phantom dosimetry. For dose simulations, dose distribution images were obtained by inputting the geometry of a CT scanner, scan parameters including estimated TCM curves and CT images of an adult anthropomorphic phantom into MC simulation software. Organ doses were then determined from the dose distribution images. For dose measurements, organ doses were evaluated using radio-photoluminescence glass dosemeters located at various organ positions within the phantom. Relative differences between the simulated and measured organ doses were -2.5 to 11.0% and -1.5 to 10.5% for organs in chest and abdomen-pelvis CT scan ranges, respectively. Thus, the simulated and measured doses agreed well.

Germain J, Pandit-Taskar N, Xu XG, Bolch WE, Dauer LT. 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.

Direct and fast measurement of CT beam filter profiles with simultaneous geometrical calibration

PURPOSE

To accurately measure the beam filter profiles from a variety of CT scanner models and to provide reference data for Monte Carlo simulations of CT scanners.

METHODS

This study proposed a new method to measure CT beam filter profiles using a linear-array x-ray detector (X-Scan 0.8f3-512; Detection Technology Inc., Espoo, Finland) under gantry rotation mode. A robust geometrical calibration approach was developed to determine key geometrical parameters by considering the x-ray focal spot location relative to the linear-array detector and the gantry's angular increment at each acquisition point. CT beam intensity profiles were synthesized from continuously measured data during a 10° gantry rotation range with calibrated detector response and system geometry information. Relative transmission profiles of nineteen sets of beam filters were then derived for nine different CT scanner models from three different manufacturers. Equivalent aluminum thickness profiles of these beam filters were determined by analytical calculation using the Spektr Matlab software package to match the measured transmission profiles. Three experiments were performed to validate the accuracy of the geometrical calibration, detector response modeling, and the derived equivalent aluminum thickness profiles.

RESULTS

The beam intensity profiles measured from gantry rotation mode showed very good agreement with those measured with gantry stationary mode, with a maximal difference of 3%. The equivalent aluminum thickness determined by this proposed method agreed well with what was measured by an ion chamber, with a mean difference of 0.4%. The determined HVL profiles matched well with data from a previous study (max difference of 4.7%).

CONCLUSIONS

An accurate and robust method to directly measure profiles from a broad list of beam filters and CT scanner models was developed, implemented, and validated. Useful reference data was provided for future research on CT system modeling.

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.

Organ Doses to Airline Passengers Screened by X-Ray Backscatter Imaging Systems

Advanced imaging technologies (AIT) are being developed for passenger airline transportation. They are designed to provide enhanced security benefits by identifying objects on passengers that would not be detected by methodologies now used for routine surveillance. X-ray backscatter imaging is one AIT system being considered. Since this technology is based on scanning passengers with ionizing radiation, concern has been raised relating to the health risks associated with these exposures. Recommendations for standards of radiation safety have been proposed by the American National Standards Institute published in ANSI/HPS N43.17-2009. A Monte Carlo based methodology for estimating organ doses received from an X-ray backscatter AIT system is presented. Radiological properties of a reference scanner including beam intensity, geometry and energy spectra were modeled based on previous studies and physical measurements. These parameters were incorporated into a Monte Carlo source subroutine and validated with comparison of simulated versus measured data. One extension of this study was to calculate organ and effective dose on a wide range of potential passengers. Computational phantoms with realistic morphologies were used including adults of 5th, 25th, 50th, 75th and 95th percentile weight, children of 5th, 50th and 95th percentile weight, and the developing fetus of 15, 25, and 38 weeks after conception. Additional sensitivity studies were performed to evaluate effects of passenger positioning within the scanner, energy spectrum and beam geometry, as well as failure mode analyses. Results for routine operations yielded a maximum effective dose to the adult and pediatric passengers of 15 and 25 nSv per screen, respectively. The developing fetus received a maximum organ dose and whole body dose of 16 nGy and 8.5 nGy per screen, respectively. The sensitivity analyses indicated that variations in positioning, energy spectra, and beam geometry yielded a range of effective doses per screen that were an order of magnitude below the ANSI recommendation.