A Review of Organ Dose Calculation Tools for Patients Undergoing Computed Tomography Scans
Article information
Abstract
Background
Computed tomography (CT) is one of the crucial diagnostic tools in modern medicine. However, careful monitoring of radiation dose for CT patients is essential since the procedure involves ionizing radiation, a known carcinogen.
Materials and Methods
The most desirable CT dose descriptor for risk analysis is the organ absorbed dose. A variety of CT organ dose calculators currently available were reviewed in this article.
Results and Discussion
Key common elements included in CT dose calculators were discussed and compared, such as computational human phantoms, CT scanner models, organ dose database, effective dose calculation methods, tube current modulation modeling, and user interface platforms.
Conclusion
It is envisioned that more research needs to be conducted to more accurately map CT coverage on computational human phantoms, to automatically segment organs and tissues for patient-specific dose calculations, and to accurately estimate radiation dose in the cone beam computed tomography process during image-guided radiation therapy.
Introduction
Computed tomography (CT) is one of the crucial diagnostic tools in modern medicine. However, careful monitoring of radiation dose for CT patients is essential since the procedure involves ionizing radiation, a known carcinogen. There is direct evidence of adverse late effects reported from multiple epidemiological studies of pediatric CT patients [1–3]. Pediatric patients, in particular, have garnered increased attention since they are known to be more vulnerable to ionizing radiation due to their more active cell development and longer life expectancy compared to adults [4]. To make an informed decision on the use of CT, it is critical to define relevant dose metrics and to accurately estimate them.
Several dose quantities concerning CT, called CT dose descriptors, are available from CT scanners [5]: weighted CT dose index (CTDIw), volumetric CT dose index (CTDIvol), and dose-length product (DLP). However, these dose descriptors are basically derived from CTDI100, the radiation dose delivered to a 100 mm long ion chamber plugged in an acrylic cylinder phantom. However, since the cylindrical shape and homogeneity of the CTDI phantoms are different from those of the human anatomy, these CT dose descriptors cannot be considered as doses delivered to a patient. There is a “derived” dose descriptor, called effective dose, which is more relevant to patient doses. Effective dose is defined as the sum of the equivalent doses delivered to major radiosensitive organs weighted by tissue weighting factors [5]. To properly estimate the effective dose, absorbed doses to multiple radiosensitive organs are required as discussed later in this article. However, DLP-to-effective dose conversion factors, called k-factor [6], can be conveniently used to derive effective doses from DLP, which is readily available from CT scanners. Nonetheless, the effective dose is defined as ageand gender-averaged quantity, so cannot be used to estimate the risk of late effects from CT scans. The most desirable CT dose descriptor for risk analysis is organ absorbed dose [7].
There are two different approaches to estimating organ doses for CT patients. First, organ doses can be measured by using a surrogate patient anatomy model, called a physical anthropomorphic phantom. Different types of radiation dosimeters are embedded into the physical phantoms before CT scans and processed to read dose after scans. The procedure is often laborious and expensive. Second, organ doses can be calculated through the computer simulations of CT scans where CT scanner simulation models and the anatomical surrogate of a patient, called computational human phantoms, are involved. The simulation approach is more flexible and cost-effective than measurement so many researchers have developed different computational approaches, some of which then led to the development of user-friendly dose calculation tools for CT patients. To date, various CT dose calculation tools have been introduced to the radiology community, but little information is available concerning the comparison of their key elements.
In the current review article, a variety of CT dose calculation tools currently available will be reviewed. Common technical elements in CT dose calculators will be described and compared. Ongoing and future research on CT patients’ dosimetry will be discussed.
Elements of CT Dose Calculators
The summary of the elements in different CT dose calculators currently available is tabulated in Table 1 [8-15]. Several web applications simply providing effective doses or wholebody fetal doses are also available. However, the scope of the review was limited to CT dose calculation tools that calculate organ-level doses. CT dose calculators are commonly composed of multiple elements: computational human phantoms, CT scanner models, organ dose database, the algorithms of effective dose calculation, and tube current modulation modeling, which are housed in different types of platforms.
1. Computational Human Phantoms
Patient anatomy is one of the most critical elements of CT dose calculators. The most accurate anatomy model in CT dose calculations would be the anatomy of a patient, which possibly can be reconstructed from CT images. However, there are several challenges in this option including automatic segmentation of organs of interest, which will be discussed later. The next best anatomy models are computational human phantoms, one of the rapidly developing areas in medical physics over the past decades [16]. The computational phantoms are developed by contouring major organs of interest from CT or magnetic resonance (MR) images of patients and are adjusted to match the dimensions of the body and internal organs to those of the reference (or standard) human models established by a comprehensive literature search such as the International Commission on Radiological Protection (ICRP) Publication 89 [17]. As shown in Table 1, different libraries of computational human phantoms are implemented into CT dose calculators: age-dependent reference size pediatric and adult phantoms, body size-dependent phantoms, and pregnant women phantoms with fetal models. Depending on the availability of patient information, CT dosimetrists can choose different sets of phantoms for organ dose calculations.
When the only information available for a patient is age and gender, one can choose age and gender-dependent phantoms with reference body sizes. When this option is adopted, the impact of body size on organ dose would be ignored. Three categories of age-dependent phantoms have been introduced in the past decades. The first series of age-dependent phantoms were created using mathematical equations with Boolean operations, which are called mathematical or stylized phantoms [18]. ImPACT [8] and WinDose [9] employed adult male and female stylized phantoms developed at the Oak Ridge National Laboratory (ORNL). ImPACT program also provides a conversion table to scale the adult dose to that of pediatric patients. ImpactDose [19] and Radimetrics [10] adopted a series of pediatric and adult stylized phantoms. However, the mathematical equations and Boolean operations used to describe the anatomical structure of the stylized phantoms are limited to model the complicated structures of the human anatomy. To address the limitation, a new generation of phantoms based on tomographic images, called voxel phantoms [20, 21], were developed from CT images of pediatric and adult patients. CT-Expo [11] was the first CT dose calculator which included pediatric voxel phantoms. However, the program displays stylized-looking pediatric phantoms, whereas the software manual shows that the voxel phantoms, Baby and Child [22], were used for dose calculations. A more complete series of pediatric voxel phantoms were developed by the University of Florida research team [21] but were not implemented into CT dose calculators. Following voxel phantoms, a new phantom construction method using polygon mesh or non-uniform rational B-spline (NURBS) was used to develop a library of pediatric and adult phantoms, called hybrid phantoms (meaning combining the advantages of stylized and voxel phantoms) [23–25]. Different series of age-dependent hybrid phantoms were implemented in the DoseWatch [12], National Cancer Institute dosimetry system for CT (NCICT) [13], and VirtualDose CT [14]. Fig. 1 shows the frontal views of the stylized adult male phantom in CT-Expo 2.1, the ICRP adult male voxel phantom in ImpactDose 2.2, and the National Cancer Institute adult male hybrid phantom in NCICT 2.0.
The impact of the variation of body sizes even at the same age on organ dose in CT dosimetry has been largely investigated and highlighted the importance of implementing body size-dependent phantoms into CT dose calculators. The modeling flexibility of hybrid phantoms made it possible to modify reference size phantoms to create body size-dependent phantoms by adjusting the shape and dimension of phantoms composed of a polygon mesh and NURBS surfaces. DoseWatch [12] features a library of body-size-dependent adult XCAT phantoms [25]. VirtualDose CT also has a series of ten body size-dependent adult male and female phantoms [26]. NCICT 2.0 [27] implements a library of 351 pediatric and adult hybrid phantoms [28] with a range of height and weight.
In addition to the anatomical differences caused by age and body size, another dramatic change in human anatomy happens during pregnancy. Due to its high radiosensitivity and risk of adverse late effects, the fetus is of particular interest in CT dosimetry. Following the early stylized pregnant women phantoms [29] which are implemented in Radimetrics [10], the voxel [30] and hybrid [30, 31] pregnant women phantoms have evolved taking advantage of new phantom technologies. DoseWatch implements a series of XCAT hybrid pregnant phantoms for CT dose calculations. A library of eight pregnant women phantoms from the University of Florida research team is adopted in NCICT 3.0 [32]. VirtualDose CT calculates the dose to the fetus using three RPI hybrid pregnant women phantoms [33].
2. CT Scanner Models
Technical parameters in CT scanners such as tube potential, bowtie filter, collimation, fan beam angle, etc., make differences in the characteristics of X-ray beams, which result in differences in radiation dose delivered to patients. To accurately implement scanner-specific beam characteristics into organ dose calculations, detailed technical parameters should be used in the Monte Carlo radiation transport simulations of CT scanners. The modeling process of CT X-rays has been published by several researchers [34–37] but usually follows the steps below. The material composition and dimension of the head and body filters can be obtained from manufacturers. Different X-ray spectra corresponding to the filters are usually generated using an X-ray spectrum generation program such as SPECT78 [38]. Bowtie filters can be explicitly modeled using the detailed composition and dimensions obtained from manufacturers. However, to increase the efficiency of Monte Carlo simulations, weighting factors estimated from a lateral dose profile in a pencil ion chamber are applied to angular sampling within a fan beam angle [36]. The X-ray source is generated at a location defined by a focal-spot-to-axis distance, which is also obtained from manufacturers. Once CT X-ray source modeling is completed, the resulting CTDI values at the center and periphery of the 16 cm and 32 cm CTDI phantoms from the simulations are validated by comparing with measurements at different tube potentials and filtrations.
However, technical specifications of CT scanners are only available from scanner manufacturers and are often proprietary. ImPACT [8] contains a library of organ dose databases generated by modeling several CT scanners based on manufacturer-provided scanner-specific technical parameters. The software also has a “matching” feature to extend the simulated CT scanner models to new scanner models which were not included in the Monte Carlo simulations. However, Monte Carlo simulations of different scanners are time-consuming, and the approach is not flexible to cover a variety of CT scanners that are currently available at clinics and will be introduced.
To make CT dose calculators more flexible, it was suggested that the dependency of organ doses on CT scanners’ characteristics becomes minimal when organ doses calculated from a scanner model are normalized to CTDIvol measured from the given scanner. That is, although the organ doses from scanner A are different from those from scanner B, the organ doses per CTDIvol of scanner A become very close to those per CTDIvol of scanner B within coefficient of variation, 5% [39]. When organ doses for a new CT scanner are sought, the CTDIvol-normalized organ dose, also called organ dose conversion coefficients, can be multiplied by CTDIvol from the given scanner. The scanner-independent organ dose algorithm has been widely adopted in most CT dose calculators available to date [39].
3. Organ Dose Database
All CT dose calculators reviewed in the current article employ a precalculated organ dose database to perform rapid dose calculations. Organs and tissues to be included in the organ dose database are selected based on their radiosensitivity following the list of organs to which ICRP assigns tissue weighting factors. According to ICRP Publication 103 [40], the essential organs and tissues such as the red bone marrow, colon, lungs, stomach, and breast with the tissue weighting factor of 0.12 are always included in the list of organs in CT dose calculators. The primary organ group is followed by the gonads (testes for males and ovaries for females) with their tissue weighting factor, 0.08. The organs with a tissue weighting factor of 0.04 are then included such as the bladder wall, esophagus, liver, and thyroid. The remaining organs include the bone surface (also known as the endosteum), brain, salivary glands (parotid, sublingual, and submandibular), and skin, of which tissue weighting factor is 0.01. The ICRP groups several organs as “remainder tissues” and assigns the tissue weighting factor of 0.12: adrenals, extrathoracic region, gallbladder, heart, kidneys, lymphatic nodes, muscle, oral mucosa, pancreas, prostate for males, small intestine, spleen, thymus, uterus/cervix for females. The full list of the radiosensitive organs and tissues is not usually available in the stylized phantoms and not always available even in the modern voxel and hybrid phantoms so that some CT dose calculators use surrogate organs and tissues. Use of the surrogate organs in CT dose calculations may result in dosimetric errors [41].
CT dose calculators adopt an approach to approximate helical scans by the summation of doses from multiple axial slices included in the given scan range of interest. To create an organ dose database, a large amount of Monte Carlo calculations is required. For a given phantom, dose calculations for a series of axial slices must be conducted from the top of the head to the bottom of the feet of the phantom with a certain slice increment. For example, when an adult phantom 175 cm tall is adopted and the slice thickness 1 cm is desired, a total of 175 Monte Carlo calculations must be conducted. The number of calculations increases by the combination of the X-ray spectrum and the number of phantoms. As mentioned in the previous section, the resulting organ doses (mGy) are normalized by CTDIvol (mGy) to establish CTDIvolto-organ dose conversion coefficients (mGy/mGy or unitless), which are independent of CT scanner models.
4. Effective Dose Calculation Methods
As mentioned above, the effective dose is widely used as a derived CT dose descriptor in CT dosimetry. However, it is often unclear to the users what the quantity means and how it is derived [42]. There are two approaches to estimating effective doses for CT patients. First, the effective dose can be calculated by using DLP to effective dose conversion factors, called k-factor [6], which is adopted by many CT dose reporting systems. The effective dose for a chest CT scan for an adult patient, for example, can be estimated by multiplying the conversion factors defined for chest scans and adults by the DLP value displayed on the CT dose page. However, the conversion factors are often blindly used without acknowledging how the values are derived and what uncertainty factors are involved. It is critical to consider several factors affecting k-factors.
k-factors change depending on which computational human phantoms are used. The k-factors reported by the American Association of Physicists in Medicine (AAPM) [6] were calculated using the ORNL stylized phantoms, of which anatomical structures are not realistic compared to human anatomy as mentioned above. The organs included in the scan range of the chest scan in the ORNL stylized phantoms may not be the same as those covered in the chest scan. Since the torso of the stylized phantoms is reported to be narrower in anteroposterior direction compared to the human body, organ doses calculated from stylized phantoms tend to overestimate the doses estimated by using realistic voxel or hybrid phantoms. This implies that the effective dose derived from the stylized phantom-based k-factors may be greater than that from voxel or hybrid phantoms. k-factors also depend on which set of tissue weighting factors is used. The ICRP has published three sets of tissue weighting factors to date in ICRP Publication 26 [43], ICRP Publication 60 [44], and ICRP Publication 103 [40]. All CT dose calculators but ImPACT program estimate and display effective dose based on the latest tissue weighting factors from ICRP Publication 103 as well as effective dose based on the tissue weighting factors from ICRP Publication 60. Different sets of tissue weighting factors change the resulting effective dose by up to 40% between ICRP Publications 60 and 103 [16]. Lastly, it is suggested to account for body size in effective dose [45, 46]. It is found that reference size phantoms-based effective dose may underestimate effective dose received by underweight patients and overestimate effective dose received by overweight patients when not accounting for body size in the calculations of organ doses and corresponding effective doses [46].
The second method to derive effective dose for CT patients directly uses organ doses calculated for a given patient not relying on k-factors. All CT dose calculators reviewed in this article provide both k-factor-based and organ dose-based effective doses. To correctly follow the definition of effective dose by ICRP, gender-specific organ doses must be averaged before summing equivalent doses weighted by tissue weighting factors [40]. This gender neutralization was reflected in k-factors, which are not specified by gender. However, when effective dose is calculated from organ doses from a male or female phantom, that is a pseudo-effective dose different from the ICRP definition.
5. Tube Current Modulation Modeling
One of the dose reduction techniques adopted in modern CT scanners is tube current modulation (TCM), which is one branch under automatic exposure control techniques. The basic idea of TCM is to modulate the electric current of the x-ray tube depending on the signal-to-noise ratio which varies at different angles within a single slice due to the oval cross-sectional shape of the human body (e.g., the thickness difference between anteroposterior and lateral geometries) and at different longitudinal locations due to the difference of tissue density (e.g., the density around the lungs is smaller than that around the pelvic area). The purpose of this technique is to maintain a consistent image signal-to-noise ratio in the CT image across a patient’s anatomy with varying degrees of attenuation to reduce radiation dose to patients while maintaining image quality.
There are two approaches to implement modulated tube currents into organ dose calculations. First, modulated tube current can be extracted from the Digital Imaging and Communications in Medicine (DICOM) header or Radiation Dose Structured Report (RDSR) after CT scans [46], which is called a retrospective approach. Because the algorithm of TCM techniques varies by CT scanner manufacturers and is proprietary, the major advantage of this post-scan method is that the actual modulated current profiles can be used for organ dose calculations. Several CT dose calculators use this algorithm to account for TCM in organ dose calculations [10, 12–14] as shown in Table 1.
However, in some cases, dosimetrists may need to predict organ doses before CT scans are conducted or calculate generic organ doses not specific to a certain CT scanner model, which is called a prospective approach. Although the algorithm of TCM varies across different manufacturers, the fundamental goal is to reduce the dose while minimizing image noise by modulating the current of the X-ray tube depending on the signal-to-noise ratio. There are generic algorithms reported to derive TCM profiles from computational human phantoms [47–50]. The amount of attenuation is calculated following the X-ray beam path starting from the source through the phantom’s anatomy to the detector within the fan beam angle at a single source rotational angle, which is repeated from 0º to 360º following a rotating CT source. Angular modulation profiles can be averaged to create a modulation profile for a single slice, which is repeated from the top of the head to the bottom of the feet in a phantom. This algorithm was adopted in DoseWatch, ImpactDose, and NCICT 2.0. Some studies report the organ dose difference between TCM and fixed tube current [13, 47, 49, 51].
6. User Interface Platforms
CT dose calculators are available on different platforms and each platform has pros and cons. First, some programs are operating as stand-alone programs. CT-Expo and ImPACT are based on the Microsoft Excel program using Visual Basic for Applications (VBA) technology. The size of the program is small and the programs will not need to be installed. One drawback of the platform is a compatibility issue occurring whenever Microsoft Excel is upgraded. Other programs such as NCICT, ImpactDose, and WinDose are operating on Windows operating system as a stand-alone program without relying on additional programs such as Microsoft Excel. NCICT is also running on Mac and Linux. Second, VirtualDose CT is running on a web-based server. The users should upload their CT parameters or images to the web application to calculate organ doses. This platform is more efficient for update and maintenance compared to the stand-alone programs. However, the operation requires an internet connection and patient data transfer may be an issue for some clinical settings although they are anonymized. Third, other programs such as DoseWatch and Radimetrics are embedded in the Picture Archiving and Communication System (PACS). The patient data are readily transferred to the dose calculation modules within the PACS. The results are directly used for dose monitoring purposes and will help to make informed decisions.
Ongoing and Future Research
There are three major areas in CT organ dosimetry that are under active development and need more work in the future. First, although the organ dose calculation methods and tools are well developed as mentioned above, it is still challenging to map the patient’s scan range onto the computational human phantoms implemented in CT dose calculators. One method to do this is to use the anatomical locations defined in general scan protocols. For example, chest scan protocols usually start from the top point of the clavicle to the middle point of the liver. One can use the same protocol ignoring the actual scan ranges conducted on a patient whenever organ doses for chest CT patients are sought. This approach is generally used for a large-scale cohort of patients. However, it is reported that general scan protocols are not always strictly followed [46] even in a clinical trial where protocols are carefully designed and provided to participating hospitals. To avoid dose errors caused by a disagreement between a patient and phantom, automatic mapping methods have been investigated [52] that compare patient CT images with predefined anatomical landmarks in computational human phantoms to map the locations of scan start and end of a patient on those of phantoms. Machine learning could be used for mapping organs from CT images to computational human phantoms [53].
Second, although reference computational human phantoms and body size-matched phantoms are used to more accurately simulate a patient’s anatomy, there are still discrepancies between a phantom and a patient. As briefly mentioned above, the most accurate anatomical model would be the one constructed from the CT images of a given patient. To utilize CT images for organ dose calculations, however, the organs of interest need to be rapidly segmented. Recent studies used deep learning-based automatic segmentation methods to contour major radiosensitive organs from CT images [15, 54, 55]. Some organs such as the intestines are still challenging to segment due to their complicated and variable shapes so the pre-contoured organ models from computational phantoms are used as alternatives. The segmented patient models are coupled with fast radiation transport methods often based on graphical processing unit (GPU)-based Monte Carlo tools [56].
Lastly, as radiation dose delivered to patients from imaging procedures involved in the image-guided radiation therapy (IGRT) increases due to the advance of modern imaging techniques, more attention has been given to the needs of the accurate estimation of imaging dose and its implementation into treatment planning. The imaging dose from a single kV cone-beam computed tomography (CBCT) acquisition may be low compared to the therapy dose. However, the normal tissue dose can be substantially increased by repeated daily imaging procedures throughout many treatment fractions. Several researchers report algorithms and doses calculated for CBCT [57, 58], but no user-friendly rapid dose calculators are available yet.
Summary
A variety of CT organ dose calculators currently available were reviewed in this article. Key common elements included in CT dose calculators were discussed and compared such as computational human phantoms, CT scanner models, organ dose database, effective dose calculation methods, tube current modulation modeling, and user interface platforms. It is envisioned that more research needs to be conducted to more accurately map CT coverage on computational human phantoms, to automatically segment organs and tissues for patient-specific dose calculations, and to accurately estimate radiation dose in the CBCT process during IGRT.
Notes
Conflict of Interest
No potential conflict of interest relevant to this article was reported.
Acknowledgements
This research was funded by the intramural research program of the National Institutes of Health, National Cancer Institute, Division of Cancer Epidemiology and Genetics.