Computational human phantom

Computational human phantoms are models of the human body used in computerized analysis. Since the 1960s, the radiological science community has developed and applied these models for ionizing radiation dosimetry studies. These models have become increasingly accurate with respect to the internal structure of the human body. As computing evolved, so did the phantoms. Graduating from phantoms based on simple quadratic equations to voxelized phantoms, which were based on actual medical images of the human body, was a major step. The newest models are based on more advanced mathematics, such as Non-Uniform Rational B-Splines (NURBS) and polygon meshes, which allow for 4-D phantoms where simulations can take place not only 3-dimensional space but in time as well. Phantoms have been developed for a wide variety of humans, from children to adolescents to adults, male and female, as well as pregnant women. With such a variety of phantoms, many kinds of simulations can be run, from dose received from medical imaging procedures to nuclear medicine. Over the years, the results of these simulations have created an assortment of standards that have been adopted in the International Commission on Radiological Protection (ICRP) recommendations.^[1]

Stylized (First-Generation) Computational Phantoms

The very first generation computational phantoms were developed for the needs to better assess organ doses from internally deposited radioactive materials in workers and patients. Until the late 1950s, the ICRP still used very simple models.^[2] In these calculations, each organ of the body was assumed to be represented as a sphere with an “effective radius.” The radionuclide of interest was assumed to be located at the center of the sphere and the “effective absorbed energy” was calculated for each organ. However, scientists attempted to model individual organs of the body and ultimately the entire human body in a realistic manner, the efforts of which led to stylized anthropomorphic phantoms—those resemble the human anatomy.

In general, stylized computational phantom is a mathematical representation of the human body, when coupled with a Monte Carlo radiation transport computer code, can be used to track the radiation interactions and energy deposition in the body. The feature of stylized computational phantom is finely tuned by adjusting individual parameters of the mathematical equations, which describes the volume, position, and shape of individual organs. Stylized computational phantom has a long history of development through the 1960s to 1980s.

MIRD Phantom

The MIRD phantom^[3] was developed by developed by Fisher and Snyder at Oak Ridge National Laboratory (ORNL) in 1960s with 22 internal organs and more than 100 sub-regions.^[4][5] It is the first anthropomorphic phantom representing a hermaphrodite adult for internal dosimetry.

Phantoms derived from MIRD

“Family” phantom series^[6]

Based on MIRD phantom, many derivations of phantoms were developed for the following decades. The major types of phantom include: stylized “Family” phantom series developed in 1980s by Cristy and Eckerman; “ADAM and EVA” developed by GSF, Germany; CAM (Computerized Anatomical Man) phantom developed by NASA unknown by the mainstream radiation protection dosimetry community, etc.

Limitation on stylized phantom

Although many efforts were undertaken to diversify and extend its applications in radiation protection, radiation therapy, and medical imaging, one cannot overcome its inborn limitation. The representation of internal organs in this mathematical phantom was crude, by capturing only the most general description of the position and geometry of each organ. With the powerful computer and tomographic imaging technologies became available in the late 1980s, the history launched a new era of voxel phantoms.

Voxel (Second-Generation) Phantoms

VIP-Man phantom developed by Dr. Xu and team at Rensselaer Polytechnic Institute in Troy, NY.^[7]

The stylized phantoms provided only basic information with a large degree of error. More accurate methods of simulating a human body were necessary to advance. To allow further research, the computer technology had to become more powerful and more readily available. This did not occur until the 1980s. The real breakthrough occurred when CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) devices could generate highly accurate images of internal organs in three dimensions and in digital format. Researchers discovered that they could take that diagnostic data and transform it into a voxel (volumetric pixel) format, essentially re-creating the human body in digital form in 3D. Today there are over 38 human phantoms in voxel format, for many different uses.^[8]

Challenges for implementation

Two major issues with development of the reference phantoms are difficulty in obtaining useful images and handling the large amount of data created from these images. CT scans give the human body a large dose of ionizing radiation – something the computational phantom was designed to circumvent in the first place. MRI images take a long time to process. Furthermore, most scans of a single subject cover only a small portion of the body, whereas a full scan series is needed for useful data. Handling this data is also difficult. While the newer computers had hard drives large enough to store the data, the memory requirements for processing the images to the desired voxel size were often too steep.^[9]

Basic development process of a voxel phantom

While there have been many voxel phantoms developed, they have all followed a similar path to completion. First, they must obtain the raw data, from CT scans, MRI imaging, or direct imaging through photography. Second, the components of the body must be segmented, or identified and separated from the rest. Third, the density of each component must be identified, along with the composition of each. Lastly, the data must be unified into a single 3D structure so it may be used for analysis.

Early developments

The earliest work on voxelized phantoms occurred independently at about the same time by Dr. Gibbs, of Vanderbilt University, and Dr. Zankl at the National Research Center for Environment and Health (GSF) in Germany.^[10][11] This occurred about 1982. Dr. Gibb’s work started with X-ray images, not CT or MRI images, for the reconstruction of a human phantom which was used for medical dose simulations. M. Zankl and team did use CT imaging to create 12 phantoms, ranging from BABY to VISIBLE HUMAN.

Advancements in voxel phantom design by country

• United States
  • Dr. Zubal and team at Yale University developed the VoxelMan phantom in 1994.^[12] This phantom was complete only from head-to-torso, and was designed specifically for improving nuclear medicine.
  • In 2000, Dr. George Xu and two students at Rensselaer Polytechnic Institute (RPI) created the VIP-Man phantom from data retrieved from the National Library of Medicine’s (NLM) Visible Human Project (VHP).^[13] This phantom was the most complex model to date, with over 3.7 billion voxels. This model was used in many studies concerning health physics and medical physics.
  • Dr. Bolch and team at the University of Florida created a set of pediatric phantoms from 2002 to 2006.^[14] Child computational phantoms had been severely underrepresented until this point. The team developed models ranging from newborn to mid-teens.

• Brazil
  • Dr. Kramer in Brazil improved on the data from the Yale Zubal team in an attempt to create a phantom similar to ICRP requirements, and created the MAX phantom.^[15]

• United Kingdom
  • The NORMAN phantom was developed by a team led by Dr. Dimbylow.^[16] This was created by analyzing MR images of a human male in 1996. In 2005, the team created a female phantom.

• Australia
  • At Flinders University, Dr. Caon and team created a torso phantom to simulate a teenage girl in 1999.^[17] The name of the phantom was ADELAIDE. This was the only teenage female phantom for a number of years.

• Japan
  • The first Asian phantom was developed by Dr. Saito and team at the Japanese Atomic Energy Research Institute (JAERI) in 2001.^[18] This was primarily used for radiation dosimetry studies.
  • Another group, led by Dr. Nagaoka at the National Institute of Information and Communications Technology (NIICT), created a male and female phantom around the same time period as the JAERI group.^[19] These were created from MR images.

• Korea
  • Many computational phantoms have been created in Korea since 2004 by Drs. Lee and Kim.^[20] Both male and female phantoms have been created. The High-Definition Reference Korean (HDRK) was created by color pictures of a cadaver, similar to the construction of the VIP-Man phantom from RPI.

• China
  • In the mid-2000s, the Chinese government authorized the creation of their own version of the VHP.^[21] The data was used by Dr. Zhang and team at the China Institute for Radiation Protection to create the CNMAN phantom, the most accurate computational phantom to date.

• Germany
  • M. Zankl and co-workers used CT images to construct a variety of individual voxel phantoms, including three paediatric ones and a woman in the 24th week of pregnancy.^{[22] [23] [24] [25] [26]}

Recent developments

Since voxel phantoms have been superseded by mesh- and NURBS-based phantoms, there has been limited development in the past few years. Many of the advancements have come through manipulations of current phantoms to correct for height and weight variations of the population.^[27]

In 2009, the International Commission on Radiological Protection (ICRP) and the International Commission on Radiation Units and Measurements (ICRU) adopted voxel-based adult male and female reference computational phantoms for the computation of organ dose coefficients ^[28] .

Statistical phantom

A computational framework was presented, based on statistical shape modelling, for construction of race-specific organ models for internal radionuclide dosimetry and other nuclear-medicine applications. The proposed technique used to create the race-specific statistical phantom maintains anatomic realism and provides the statistical parameters for application to radionuclide dosimetry.^[29]

Boundary Representation (BREP) phantom

4-D BREP Phantom used to model a breathing human torso^[30]

Boundary representation (BREP) phantoms are computational human models that contain exterior and interior anatomical features of a human body using boundary representation method. In the realm of health and medical physics they are primarily used for ionizing radiation dosimetry.

In the development of computational human phantoms, of particular interest is the concept of a “deformable” phantom whose geometry can be conveniently transformed to fit particular physical organ shapes, volumes, or body postures. Design of this type of phantom is realized by Non-Uniform Rational B-Spline (NURBS) method or polygonal mesh method, which are usually collectively called BREP methods. Compared to the voxel phantoms, BREP phantoms are better suited for geometry deformation and adjustment, because a larger set of computerized operations are available, such as extrusion, chamfering, blending, drafting, shelling and tweaking. A major advantage of BREP phantoms is their ability to morph into an existing reference phantom or into the anatomy of a real worker or patient, which makes individual-specific dose calculation possible.^[31]

NURBS-based phantom

Surfaces of a NURBS-based phantom are defined by NURBS equations which are formulated by a set of control points. The shape and volume of a NURBS surface vary with the coordinates of control points. This feature is useful in designing a time-dependent 4D human body modeling.^[32] An example is given by NCAT phantoms by Segars et al, which is used to simulate cardiac and respiratory motions with more realistic modeling of the cardiac system.

Polygonal mesh-based phantom

A polygonal mesh is composed of a set of vertices, edges, and faces that specify the shape of a polyhedral object in 3D space. The surfaces of the phantom are defined by a large amount of polygonal meshes, most commonly triangles. The polygonal mesh has three remarkable advantages in developing whole-body phantoms. Firstly, mesh surfaces depicting human anatomy can be conveniently obtained from real patient images or commercial human anatomy mesh models. Secondly, the polygonal mesh-based phantom has considerable flexibility in adjusting and fine-tuning its geometry, allowing the simulation of very complex anatomies. Thirdly, many commercial computer aided design (CAD) software, such as Rhinoceros, AutoCAD, Visualization Toolkit (VTK), provide built-in functions able to rapidly convert polygonal mesh into NURBS.^[33]

Development

Pregnant (9 months) adult female phantom^[34]

Segars was the precursor of applying NURBS to phantom design. In 2001 his doctoral thesis described the method of developing a dynamic NURBS-based cardiac-torso (NCAT) phantom in detail. The phantom has a 4D beating heart model which was derived from 4D tagged magnetic resonance image (MRI) data. The remaining organs in the torso of the phantom were designed based on the Visible Human Project CT data set and were composed of 3D NURBS surfaces. Respiratory motion was also incorporated into this phantom.

In 2005, Xu et al. at Rensselaer Polytechnic Institute (RPI) used the 3D VIP-Man phantom to simulate respiratory motions by adopting the gated respiratory motion data of the NCAT phantom.^[35] The 4D VIP-Man Chest phantom was used to study external-beam treatment planning for a lung cancer patient.^[36] In 2007, Xu’s research group reported creation of a series of polygon-based phantoms representing a pregnant woman and her fetus at the end of 3, 6, and 9 month gestations (RPI Pregnant Females).^[37] The mesh data were initially obtained from separate anatomical information sources including a non-pregnant female, a 7-month pregnant woman CT data set, and a mesh model of the fetus. In 2008, two triangular mesh-based phantoms were created, named as RPI Deformable Adult Male and Female (RPI-AM, RPI-FM).^[38][39] The anatomic parameters of the phantoms were made consistent with two datasets: the mass and density of internal organs originated from ICRP-23 and ICRP-89, and the whole-body height and weight percentile data were obtained from the National Health and Nutrition Examination Survey (NHANES 1999-2002). Later on, to study the relationship between breast size and lung dosimetry, a new group of phantoms were produced by altering the breast geometry of RPI-AF.^[40]

From 2006 to 2009, Bolch et al. at University of Florida (UF) designed a total of twelve “hybrid” male and female phantoms, representing newborn, 1, 5, 10, 15-year old and adult male/females.^[41][42][43] The phantoms are addressed as “hybrid” because most organs and tissues were modeled by NURBS surfaces whereas the skeleton, brain and extra-thoracic airways were modeled by polygonal surfaces.^[44] Anatomic parameters of the phantoms were adjusted to match 4 reference datasets, i.e., standard anthropometric data, reference organ masses from ICRP Publication 89, reference elemental compositions provided in ICRP 89 as well as ICRU Report 46, and reference data on the alimentary tract organs given in ICRP Publications 89 and 100.

In 2008, Stabin et al. at Vanderbilt University, in collaboration with Segars from Duke University, developed a family of adult and pediatric phantoms by adapting the NURBS-based NCAT adult male and female phantoms.^[45] ICRP-89 reference body and organ values were used to adjust NURBS surfaces.

In 2009 Cassola et al^[46] at the Federal University of Pernambuco, Brazil, developed a pair of polygonal mesh-based phantoms in standing posture, FASH (Female Adult meSH) and MASH (Male Adult meSH). The methodology is very similar but not entirely identical to the one implemented in the designing of RPI-AM and RPI-FM.

In 2010, based on existing RPI-AM, Xu et al. continued to create 5 more phantoms with different Body Mass Index (BMI) ranging from 23 to 44 kg∙m-2.^[47] These phantoms are used to study the correlation between BMI and organ doses resulting from CT and PET examinations.

In 2011 Kim et al. at Hanyang University, Korea, reported a polygon-surface reference Korean male phantom (PSRK-Man).^[48] This phantom was constructed by converting the Visible Korean Human-Man (VKH-man) into a polygonal mesh-based phantom. The height, weight, geometry of organs and tissues were adjusted to match the Reference Korean data. Without voxelization the PSRK-man could be directly implemented in Geant4 Monte Carlo simulation using a built-in function, but the computation time was 70~150 times longer than that required by High Definition Reference Korean-Man (HDRK-Man), a voxelized phantom derived also from VKH-man.

External Links

• List of computational human phantoms
• Consortium of Computational Human Phantoms (CCHP)
• Rensselaer Radiation Measurement and Dosimetry Group
• Helmholtz Zentrum München, Department of Radiation Sciences, Research Unit Medical Radiation Physics and Diagnostics

References

1. ^ Xu, X.G.; Eckerman, K.F. Handbook of Anatomical Models for Radiation Dosimetry. Taylor & Francis, 2010. ISBN 978-1-4200-5979-3.
2. ^ ICRP. Report of Committee II on Permissible Dose for Internal Radiation International Commission on Radiological Protection (Oxford: Pergamon Press), 1959.
3. ^ Report of the Task Group on Reference Man: ICRP Publication 23.
4. ^ Fisher, H.L.J. and Snyder, W.S. “Variation of dose delivered by 137Cs as a function of body size from infancy to adulthood.” ORNL-4007 (Oak Ridge, TN: Oak Ridge National Laboratory), P. 221, 1966.
5. ^ Fisher, H.L.J. and Snyder, W.S. “Distribution of dose delivered in the body size from a source of gamma rays distributed uniformly in an organ”, ORNL-4168 (Oak Ridge, TN: Oak Ridge National Laboratory), p. 245, 1967.
6. ^ Kramer, R. et al. All about FAX: A female adult voXel phantom for Monte Carlo calculation in radiation protection dosimetry, Phys Med Biol, 49, 5203, 2004.
7. ^ Photo courtesy of Dr. George Xu, Rensselaer Polyechnic Institute http://www.rpi.edu/dept/radsafe/public_html/project.htm
8. ^ Zaidi, H. and Xu, X.G. (2007). “Computational anthropomorphic models of the human anatomy: The path to realistic Monte Carlo modeling in radiological sciences”, Annu Rev Biomed Eng. 9,p. 471.
9. ^ Xu, X.G.; Eckerman, K.F. Handbook of Anatomical Models for Radiation Dosimetry. Taylor & Francis, 2010. ISBN 978-1-4200-5979-3.
10. ^ Gibbs, S. and Pujol, J. (1982). “A Monte Carlo method for patient dosimetry from diagnostic x-ray.” Dentomaxillofac Radiol. 11, p. 25.
11. ^ Zankl, M. et al. (1988). “The construction of computer tomographic phantoms and their application in radiology and radiation protection.” Radiat Environ Biophys, 27, p. 153.
12. ^ Zubal, I.G. et al. (1994). “Computerized three-dimensional segmented human anatomy.” Med Phys, 21, p. 299.
13. ^ Xu, X.G., Chao, T.C., and Bozkurt, A. (2000) “VIP-Man: An image-based whole-body adult male model constructed from color photographs of the Visible Human Project for multi-particle Monte Carlo calculations.” Health Phys, 78, p. 476.
14. ^ Lee, C. et al. (2006). “Whole-body voxel phantoms of paediatric patients—UF Series B.” Phys Med Biol, 51, p. 4649.
15. ^ Kramer, R. et al. (2003). “All about MAX: A male adult voxel phantom for Monte Carlo calculations in radiation protection dosimetry.” Phys Med Biol, 48, p. 1239.
16. ^ Dimbylow, P.J. (1996). “The development of realistic voxel phantoms for electromagnetic field dosimetry,” in Proceedings of the Workshop on Voxel Phantom Development, Chilton, U.K.
17. ^ Caon, M., Bibbo, G., and Pattison, J. (1999). “An EGS4-ready tomographic computational model of a 14-year-old female torso for calculating organ doses from CT examinations.” Phys Med Biol, 44, p. 2213.
18. ^ Saito, K. et al. (2001). “Construction of a computed tomographic phantom for a Japanese male adult and dose calculation system.” Radiat Environ Biophys, 40, p. 69.
19. ^ Nagaoka, T. et al. (2004). “Development of realistic high-resolution whole-body voxel models of Japanese adult males and females of average height and weight, and application of models to radio-frequency electromagnetic-field dosimetry.” Phys Med Biol, 49, p. 1.
20. ^ Kim, C.H. et al. (2008). “HDRK-Man: A whole-body voxel model based on high-resolution color slice images of a Korean adult male cadaver.” Phys Med Biol, 53, p. 4093.
21. ^ Zhang, B.Q. et al. (2007). “CNMAN: A Chinese adult male voxel phantom constructed from color photographs of a visible anatomical data set.” Radiat Prot Dosim, 124, p. 130.
22. ^ Fill, U.; Zankl, M., Petoussi-Henss, N., Siebert, M., Regulla, D. (2004). "Adult female voxel models of different stature and photon conversion coefficients for radiation protection". Health Physics 86: 253–272. 
23. ^ Petoussi-Henss, N.; Zankl, M., Fill, U., Regulla, D. (2002). "The GSF family of voxel phantoms". Phys. Med. Biol. 47: 89–106. 
24. ^ Zankl, Maria (2010). "The GSF voxel computational phantom family. Handbook of anatomical models for radiation dosimetry". Handbook of anatomical models for radiation dosimetry: 65–85. 
25. ^ Zankl, M.; Veit, R., Williams, G., Schneider, K., Fendel, H., Petoussi, N., Drexler, G (2001). "The construction of computer tomographic phantoms and their application in radiology and radiation protection". Radiat. Environ. Biophys. 40: 153–162. 
26. ^ Zankl, M.; Wittmann, A. (2001). "The adult male voxel model "Golem" segmented from whole body CT patient data". Radiat. Environ. Biophys. 40: 153–162. 
27. ^ Kramer, G.H. et al. (2010). “Tools for Creating and Manipulating Voxel Phantoms.” Health Physics. 98(3), pp 542-548. doi: 10.1097/HP.0b013e3181c34ced.
28. ^ ICRP (2009). Adult reference computational phantoms. ICRP Publication 110. Oxford, UK. 
29. ^ Mofrad, F. B et al. (2010). “Statistical construction of a Japanese male liver phantom for internal radionuclide dosimetry.” Radiat Prot Dosimetry. 140(2), pp 140-148. doi: 10.1093/rpd/ncq164.
30. ^ Photo courtesy of Dr. George Xu, Rensselaer Polyechnic Institute http://www.rpi.edu/dept/radsafe/public_html/project.htm
31. ^ Na, YH. et al. Deformable adult human phantoms for radiation protection dosimetry: anthropometric data representing size distributions of adult worker populations and software algorithms, Phys Med Biol, 55, 3789, 2010.
32. ^ Na, YH. et al. Deformable adult human phantoms for radiation protection dosimetry: anthropometric data representing size distributions of adult worker populations and software algorithms, Phys Med Biol, 55, 3789, 2010.
33. ^ Na, YH. et al. Deformable adult human phantoms for radiation protection dosimetry: anthropometric data representing size distributions of adult worker populations and software algorithms, Phys Med Biol, 55, 3789, 2010.
34. ^ Photo courtesy of Dr. George Xu, Rensselaer Polyechnic Institute http://www.rpi.edu/dept/radsafe/public_html/project.htm
35. ^ Xu, X.G. and Shi, C. Preliminary development of a 4D anatomical model for Monte Carlo simulations, Monte Carlo 2005 Topical Meeting. The Monte Carlo Method:Versatility Unbounded In a Dynamic Computing World, Chattanooga, TN, April 17–21, 2005 2005.
36. ^ Zhang, J. et al. Development of a geometry-based respiratory motion-simulating patient model for radiation treatment dosimetry, Journal of Applied Clinical Medical Physics, 9, 16, 2008.
37. ^ Xu, X.G. et al. A boundary-representation method for designing whole-body radiation dosimetry models: pregnant females at the ends of three gestational periods--RPI-P3, -P6 and -P9, Phys Med Biol, 52, 7023, 2007.
38. ^ Xu, X.G, Zhang, J.Y., and Na, Y.H. Preliminary Data for Mesh-Based Deformable Phantom Development: Is it Possible to Design Person-Specific Phantoms On-demand. The International Conference on Radiation Shielding-11, April 14–17, 2008.
39. ^ Hegenbart,L. et al. A Monte Carlo study of lung counting efficiency for female workers of different breast sizes using deformable phantoms Phys. Med. Biol. 53, 5527, 2008.
40. ^ Hegenbart,L. et al. A Monte Carlo study of lung counting efficiency for female workers of different breast sizes using deformable phantoms Phys. Med. Biol. 53, 5527, 2008.
41. ^ Lee, C. et al. Hybrid computational phantoms of the male and female newborn patient: NURBS-based whole-body models, Phys Med Biol, 52, 3309, 2007.
42. ^ Lee, C. et al. Hybrid computational phantoms of the 15-year male and female adolescent: Applications to CT organ dosimetry for patients of variable morphometry, Medical Physics, 35, 2366, 2008.
43. ^ Lee, C. "The UF family of reference hybrid phantoms for computational radiation dosimetry" Phys. Med. Biol. 55, no.2 (2010): 339–363.
44. ^ Stabin, M. et al. ICRP-89 based adult and pediatric phantom series, J NUCL MED MEETING ABSTRACTS, 49, 14, 2008.
45. ^ Stabin, M. et al. ICRP-89 based adult and pediatric phantom series, J NUCL MED MEETING ABSTRACTS, 49, 14, 2008.
46. ^ , VF. "FASH and MASH: female and male adult human phantoms based on polygon mesh surfaces: I. Development of the anatomy." Phys. Med. Biol. 55, no. 133 (2010): 133–162.
47. ^ A. Ding, M. Mille, P. F. Caracappa, X. G. Xu, “Impact of body size of obese patients on PET/CT dose estimates: Monte Carlo calculations using a set of BMI-adjustable phantoms,” 53rd Annual Meeting of the American Association of Physicists in Medicine, Vancouver, British Columbia, Canada, July 31-August 4 (2011).
48. ^ Kim, CH. "A polygon-surface reference Korean male phantom (PSRK-Man) and its direct implementation in Geant4 Monte Carlo simulation" Phys. Med. Biol. 56, no.10 (2010): 3137–3161.