Abstract
A new generation of reference computational phantoms, based on image-based models tied to the reference masses defined by the International Commission on Radiological Protection (ICRP) for dose calculations, is presented. Methods: Anatomic models based on nonuniform rational b-spline modeling techniques were used to define reference male and female adults, 15-y-olds, 10-y-olds, 5-y-olds, 1-y-olds, newborns, and pregnant women at 3 stages of gestation, using the defined reference organ masses in ICRP publication 89. Absorbed fractions and specific absorbed fractions for internal emitters were derived using standard Monte Carlo radiation transport simulation codes. Results: Differences were notable between many pairs of organs in specific absorbed fractions because of the improved realism of the models, with adjacent organs usually closer and sometimes touching. Final estimates of absorbed dose for radiopharmaceuticals, for example, were only slightly different overall, as many of the differences were small and most pronounced at low radiation energies. Some new important organs were defined (salivary glands, prostate, eyes, and esophagus), and the identity of a few gastrointestinal tract organs changed. Conclusion: A new generation of reference models for standardized internal and external dose calculations has been defined. The models will be implemented in standardized software for internal dose calculations and be used to produce new standardized dose estimates for radiopharmaceuticals and other applications.
Anthropomorphic models—also called computational phantoms—applied for dose calculations in nuclear medicine for the past 30 y have used the stylized geometric definitions that were developed for reference adults and children (1) and pregnant women (2). Figure 1 shows an exterior view of one of the Cristy/Eckerman stylized adult and pediatric models. The models’ exteriors were defined in 3 sections: an elliptic cylinder representing the arm, torso, and hips; a truncated elliptic cone representing the legs and feet; and an elliptic cylinder representing the head and neck. Several organs and tissues were mathematically defined as occupying finite spaces within the body space and comprised 3 types of tissue: soft tissue, bone, and lung. The mathematic descriptions of the organs were formulated on the basis of descriptive and schematic materials from general anatomy references. This series of stylized models was implemented in the MIRDOSE (3) and OLINDA/EXM 1.0 (4) personal computer codes to facilitate calculation of standardized internal dose calculations.
Exterior view of one of the Cristy/Eckerman 1987 phantoms (1).
The stylized models that originated at Oak Ridge National Laboratory some 40 y ago have now been replaced with realistic body models based mostly on human image data (5). The image-based models are now considered to be the second-generation anatomic models that use the so-called voxel geometry. Petoussi-Henss et al. created tomographic voxelized models from CT images (6) and found that the stylized computational models could over- and underestimate organ absorbed fractions by several tens of percentage points when compared against the more realistic, image-based models. Similar discrepancies between the stylized and voxelized models were reported using the VIP-Man model, which is based on the Visible Human Project cadaver cross-sectional anatomic images (7,8). Hurtado et al. (9) reported on hybrid realistic computational phantoms for adults and used them in retrospective dosimetry analyses for contaminated radiation workers. Many others have developed individual realistic phantoms of adults (10,11), children (12), and pregnant women (13,14), using various technologies, as summarized by Xu and Eckerman (5). In particular, representing the third generation of anatomic model is the work by Segars, now with Duke University, who used the nonuniform rational b-spline (NURBS) modeling technique (15) that has allowed not only a more realistic rendering of the human body but also the attractive capability of rapid and easy scaling of organs and the whole body. Figure 2 shows the adult male and female NURBS models originally developed by Segars.
Anterior views of Segars (15) NURBS models of adult male (A) and adult female (B).
This work describes the development of a complete male and female adult and pediatric (15-y-old, 10-y-old, 5-y-old, 1-y-old, and newborn) and pregnant woman phantom series for the purposes of internal radiation dosimetry in nuclear medicine using a realistic, NURBS-based body modeling technique, based on the recommended organ, body, and fetal masses in publication 89 of the International Commission on Radiological Protection (ICRP) (16). We have calculated photon- and electron-specific absorbed fractions (SAFs) and dose factors for more than 1,000 radionuclides for this phantom series, as defined in the unified internal and external dose assessment methodology of the RAdiation Dose Assessment Resource (RADAR) for nuclear medicine patient and occupational radiation worker dosimetry (17):
In certain cases (e.g., hollow organs and bone or marrow), results from other models are used to provide a complete dosimetric modeling environment that has been implemented for internal dose assessment in an update of the OLINDA/EXM personal computer code (4).
MATERIALS AND METHODS
Adult and Pediatric Models
NURBS models may be easily deformed to represent different sizes and shapes. We developed a series of models representing adults and children of different ages scaled from these original adult models to match the recommended age-dependent organ masses given in ICRP publication 89 (16). Existing software tools developed by Segars allow for the rapid scaling of human NURBS. One or more selected organs may be translated or rotated in any direction, scaled linearly in any direction (uniformly in 3 dimensions from the center by a fixed factor), and otherwise modified by the user. Slice-by-slice viewing of results is also afforded by this software tool. We started with Segars’ 2 adult models, scaled them to the ICRP 89 reference adult organ masses, and then continued to modify them to represent the smaller models by deforming the surfaces and then adjusting the organs to match the ICRP 89 models for each age. In general, the whole models were not simply scaled equally in all 3 dimensions; for example, the head size is much larger relative to the size of the whole body in young children than in adults. Overall scaling of the body was done to provide an appropriate body size and shape, and then organs were individually scaled. The ICRP 89–recommended organ mass values for adults and children of various reference ages (male and female newborn, 1-y-old, 5-y-old, 10-y-old, 15-y-old, and adult) were used as target masses for individual organs in each phantom. Mostly, organs were varied in all 3 dimensions equally, and then individual organs could be translated if needed to maintain internal consistency between the models. Age-dependent differences (e.g., head size in young children) were incorporated. When finalized, the modified models were rendered in voxel format with another software tool (using a resolution of 1.5 × 1.5 × 5 mm) and converted from binary to ASCII format. Hollow organs had both wall and content components; target masses from ICRP 89 were applied to both the wall and the content components. Organs were touching and even overlapping in some cases; Segars’ models have a priority scheme to handle this, which is invoked during voxelization. All original file sizes were 512 × 512 × 512; these were trimmed to a variable cubic format that eliminated excess rows and columns of voxels outside the phantom and were used with a GEANT4 (GEometry ANd Tracking)-based (18) Monte Carlo radiation transport program that we developed. Separate male and female models were made for the adult and pediatric models. Because the ICRP does not offer separate organ mass values for 10-y-old models and younger, the same model was used for both sexes, with the sex organs added as appropriate. For most organs, the difference between the NURBS reported and voxelized volumes were usually within 1%; for some small organs, however, the differences were sometimes greater. Absorbed fractions were converted to SAFs using the reference organ masses (Tables 1 and 2). Generally, 200,000 starting particles were used at each energy for each source organ. Uncertainties in calculated SAF values were generally less than 5% for most organs; for some small or distant organs or very low-energy particles, the uncertainties were necessarily higher. Reciprocity (19), averaging, and curve-smoothing methods were used to deal with noisier data. Many values at low energies were of high uncertainty or no hits were observed; these were reported as zeroes, and no attempts were made to extrapolate curves beyond the range of the reported data. For the pregnant female models, as noted by Shi et al. (14): “The Monte Carlo run employed approximately 1 × 107 histories to keep the statistical uncertainty to less than 10% for most of the target organs, except for photon energies below 30 keV and for some target organs that were too small or too far from the source organ. In the latter cases, uncertainties of 10% or more were observed. The results were discarded if the relative uncertainty was higher than 40% so we could include as many data points as possible. The cutoff energy for both electron and photon calculations in the EGS4-VLSI user code were set to be 10 keV.”
Selected Reference Organ Masses (g) for Pediatric/Adult Phantom Series
Selected Reference Organ Masses (g) for RPI Pregnant Female Phantom Series
Pregnant Female Models
A complete description of the methods used to develop SAFs for the pregnant female models developed at Rensselaer Polytechnic Institute (RPI) in collaboration with Vanderbilt has been published previously (13,14). Briefly, organ models were created by extracting the voxel information contained in the so-called boundary representations that consisted of polygonal meshes or NURBS. Individual organs and fetal models were integrated into the adult female whole body by careful adjustment of the organ shape and location to accommodate gestation-related anatomic changes in the body and to avoid organ overlap. For each organ, the volume and mass were specified manually according to reference values recommended in ICRP 89 (16). Once all the organs and total body weights had been adjusted, the surface models were voxelized for dose calculations using the EGS4 Monte Carlo code (20).
RESULTS
Figure 3 shows images of selected models from the adult and pediatric ICRP 89 phantom series, and Figure 4 shows renderings of the RPI pregnant female models. Figure 5 shows sample photon and electron SAFs for adult and pediatric models. Traditionally, all AFs for electron self-irradiation have been assumed to be 1.0 for an organ irradiating itself and 0.0 for other organs. Although this trend is generally true, when electron transport is performed there is some departure from this assumption at higher energies. This trend is more important for the smaller organs, particularly in the smaller phantoms.
Rendering of selected NURBS ICRP 89 phantoms: newborn female model (A), 5-y-old male model (B), 10-y-old female model (C), and 15-y-old male model (D).
Rendering of 3-dimensional organ and body surfaces for RPI pregnant female models: 3 mo (A), 6 mo (B), and 9 mo (C).
Sample SAF plots for ICRP 89 pediatric and adult models: male, photons stomach←liver (A), female, photons lungs←lungs (B), male, electrons kidneys←kidneys (C), and female, adrenals←kidneys (D).
DISCUSSION
The first generation of stylized anthropomorphic models has served the dosimetry community well for over 40 y for internal and external dose applications (5). Original body models used in dosimetry were spheric (21); this was mathematically convenient and provided conservative estimates of SAFs and doses but was not at all realistic, and there was no provision for organ cross-irradiation. The development of the stylized adult and pediatric phantom series by Cristy and Eckerman (1) (following the original design of Snyder et al. (22)) represented an important incremental change in methods for development of SAFs and a significant increase in model realism over a simple, uniform sphere. Monte Carlo methods were developed to estimate SAF values for the organs, and the science of dose calculations for internal and external simulated sources of radiation was greatly improved. However, the advent of image-based technologies that allowed the development of realistic body models for use in dose calculations brought an exciting opportunity to perform patient-individualized dose calculations in nuclear medicine therapy (23), as well as general dose assessment for internal or external sources. The models presented here represent the third generation of body models, with the organ masses based on those recommended by the ICRP for standardized individuals of various ages (16).
Realistic, voxel-based models derived from real human images have notable differences in organ proximity, compared with the traditional, stylized models, as was illustrated by Stabin and Yoriyaz (24). Organs are closer together, and many are in direct contact with each other (kidneys/adrenals, lungs/heart), whereas in stylized models, separation of organ spaces occurs because of the simplicity of the shapes used to model them. Changes in actual SAF values from the Cristy/Eckerman model series can be observed, but they are mostly small. This provides indirect validation of these results (because validation using physical phantoms is not possible with available technology); we also checked that the SAFs at low photon energies were approximately 1/(organ mass). Because the SAFs are similar to those of the previous generation of phantoms, the impact on calculated dose estimates for radiopharmaceuticals, other internal sources, or external sources should be minor but of interest for evaluation. Photon SAFs for lungs←liver for all 6 ages in the male phantoms demonstrate trends similar to those of the previous stylized phantoms, but the SAFs for the realistic phantoms are generally a bit higher because of the more realistic modeling of organ proximity (one lung actually rests on top of the liver). At very low energies, the differences are more pronounced, but these photons are typically not important to a complete dose assessment for a given radiopharmaceutical. In the adult female for lungs←liver, significant differences are noted at low energies, but at moderate to high energies the values are similar but still slightly higher. For stomach←spleen, the NURBS values are several orders of magnitude higher in the keV range and noticeably higher at all energies. In calculation of dose estimates for individual radiopharmaceuticals, doses to most organs are similar, but doses to certain organs (e.g., adrenals and pancreas) may be notably different. However, overall effects on effective dose values are small (to be shown in a subsequent publication).
It is useful to now explicitly treat electron SAFs. The OLINDA/EXM computer code (4) included an approximate correction for loss of electron energy in small objects at high electron energies, based on the AFs developed in unit-density spheres by Stabin and Konijnenberg (25). Explicit Monte Carlo transport provides a better estimation of this electron energy loss at high energies, which can be noted even for large organs such as kidneys (Fig. 5C) but more dramatically for small organs such as ovaries and testes. Deposition of electron energy on other organs is interesting. In most cases, this is simply a small contribution due to bremsstrahlung radiation, which is known to be an insignificant contribution to the overall dose in any internal dose problem (26). However, in cases in which organs are touching (e.g., kidneys/adrenals), β-energy is deposited in the first few millimeters of the tissue in the adjacent organ as well (Fig. 5D). Because these SAFs are for whole organs, inclusion of this component is not equivalent to averaging photon energy over the entire organ, which will be generally more uniform. But there is energy deposited, and we decided to leave this contribution in the reported SAFs.
Many other groups, including Petoussi-Henss et al. (6), Kramer et al. (28), Xu et al. (7), Zubal et al. (29), Lee et al. (30), and others, have developed standardized phantoms ultimately adopted by the ICRP (27). The ICRP 110 models also are based on the ICRP 89 organ weights but use manually segmented and adjusted organ representations from 2 human subjects. As such, the numeric values of SAFs from these models should be similar to those in the RADAR phantoms, because both model realistic organ geometries and use the same organ masses. Evaluation of several cases shows that this is the case (details to be provided in a later publication). We maintained the OLINDA/EXM 1.0 and 1.1 bone and marrow model (based on the modified Eckerman model) (31); we feel that this is an advantage in that it provides continuity between the codes in this area. We are aware that better models are under development (32), but they are not complete for all ages and sexes at present. We also maintained the traditional model for electron irradiation of hollow organ walls. The ICRP addressed the issue of the conservatism of assuming that the dose to the wall–content interface is one half that to the contents, but only for adults (33). We did explicitly model electron and photon transport using hollow-organ contents as the source, but we could not model the exact locations of sensitive cells in the walls, thickness of mucus layers, and other factors. The photon SAFs from the transport were used, but OLINDA/EXM will use the traditional approach for electrons until more complete data are available.
The RADAR phantom series, including the adult, pediatric, and pregnant female models, will be implemented in future versions of the OLINDA/EXM computer code (4), which has facilitated and standardized dose estimates for many years, building on the technical basis of the previous MIRDOSE software (3). New tables of dose estimates for radiopharmaceuticals in current use—organ doses and effective doses—will be generated with these new SAFs and will serve as an update to these dose estimates using this new generation of anthropomorphic phantoms. OLINDA/EXM still treats only average, whole-organ dose, based on an assumed uniform organ uptake. This approach is certainly appropriate for diagnostic applications, but for new drug applications and clinical evaluations this kind of dosimetry is routinely reported for therapeutic radiopharmaceuticals as well. More advanced methods have been demonstrated by several groups as research applications using 3-dimensional distributions of activity to produce nonhomogeneous dose distributions and dose–volume histograms for patient-individualized anatomic models (23,34). There is no present intention for OLINDA/EXM to treat this kind of detailed dosimetry; the role of this code is to provide average organ doses for standardized models.
CONCLUSION
A new generation of reference models for standardized internal and external dose calculations has been defined. The models will be implemented in standardized software for internal dose calculations and be used to produce new standardized dose estimates for radiopharmaceuticals and other applications.
DISCLOSURE STATEMENT
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
Work on this project was supported by grant NIH/NCI 1R42CA115122. No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Sep. 11, 2012.
- © 2012 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication March 20, 2012.
- Accepted for publication May 29, 2012.