Computational human phantom

From Wikipedia, the free encyclopedia

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-spline (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[edit]

The very first generation computational phantoms were developed to address the need 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. Phantoms such as the Shepp-Logan Phantom were used as models of a human head in the development and testing of image reconstruction algorithms.[3][4][5][6] 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 that resemble the human anatomy.

In general, stylized computational phantom is a mathematical representation of the human body which, 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[edit]

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

Phantoms derived from MIRD[edit]

"Family" phantom series[10]

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 the 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[edit]

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[edit]

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 computed tomography (CT) and magnetic resonance imaging (MRI) 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.[11]

Challenges for implementation[edit]

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.[1]

Basic development process of a voxel phantom[edit]

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[edit]

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.[12][13] 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[edit]

  • United States
    • Dr. Zubal and team at Yale University developed the VoxelMan phantom in 1994.[14] This original phantom was complete only from head-to-torso, and was designed specifically for improving nuclear medicine. Since its initial development it has been improved to include arms and legs to represent a complete human body, and a dedicated head has been completed which delineates the small internal sub-structures of the brain.[15]
    • 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).[16] 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.[17] Child computational phantoms had been severely underrepresented until this point. The team developed models ranging from newborn to mid-teens.
    • The U.S. Food and Drug Administration (FDA) developed voxel-based virtual family body phantoms [18] for X-ray dose study. Dr. Gu and Dr. Kyprianou improved the heart parts with male/female high-resolution computational heart phantoms (both voxel/mesh based) in 2011.[19] The main contribution is that the detail level of coronary arteries are visible in those phantoms.
  • Brazil
  • United Kingdom
    • The NORMAN phantom was developed by a team led by Dr. Dimbylow.[21] This was created by analyzing magnetic resonance 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.[22] The name of the phantom was ADELAIDE. This was the only teenage female phantom for a number of years.
  • Japan
  • Korea
    • Many computational phantoms have been created in Korea since 2004 by Drs. Lee and Kim.[25] 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.[26] 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.[27][28][29][30][31]

Recent developments[edit]

Statistical phantom[edit]

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.[32]

Boundary representation (third-generation) phantom[edit]

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.[33]

NURBS-based phantom[edit]

Surfaces of a non-uniform rational B-spline (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.[33] 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[edit]

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]


Real-life motion data (left) is acquired on a motion capture platform (center) and used to determine the posture of the CHAD phantom (right)[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 (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 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.[39]

From 2006 to 2009, researchers at the University of Florida designed a total of twelve "hybrid" male and female phantoms, representing newborn, 1-, 5-, 10-, and 15-year-old and adult male/females.[40][41][42] 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.[43] 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, researchers at Vanderbilt University, in collaboration with researchers from Duke University, developed a family of adult and pediatric phantoms by adapting the NURBS-based NCAT adult male and female phantoms.[43] ICRP-89 reference body and organ values were used to adjust NURBS surfaces.

In 2009 Cassola et al.[44] 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, researchers at RPI continued to create 5 more phantoms with different body mass index (BMI) ranging from 23 to 44 kg∙m-2.[45] These phantoms are used to study the correlation between BMI and organ doses resulting from CT and positron emission tomography (PET) examinations.

In 2011 researchers at Hanyang University, Korea, reported a polygon-surface reference Korean male phantom (PSRK-Man).[46] 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.

In 2012, researchers at RPI developed the Computational Human for Animated Dosimetry (CHAD) phantom, structured such that its posture could be adjusted in conjunction with data obtained using a motion capture system.[47] This phantom can be used to simulate the movement of a worker involved in an occupational of nuclear accident scenario, allowing researchers to gain an understanding of the impact of changing posture in the course of worker movement on radiation dose.

See also[edit]


  1. ^ a b 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. ^ Shepp, Larry; B. F. Logan (1974). "The Fourier Reconstruction of a Head Section". IEEE Transactions on Nuclear Science. Vol. NS-21.
  4. ^ Ellenberg, Jordan (February 22, 2010). "Fill in the Blanks: Using Math to Turn Lo-Res Datasets Into Hi-Res Samples". Wired. Retrieved 31 May 2013.
  5. ^ Müller, Jennifer L.; Siltanen, Samuli (2012-11-30). Linear and Nonlinear Inverse Problems with Practical Applications. SIAM. pp. 31–. ISBN 9781611972337. Retrieved 31 May 2013.
  6. ^ Koay, Cheng Guan; Joelle E. Sarlls; Evren Özarslan (2007). "Three-Dimensional Analytical Magnetic Resonance Imaging Phantom in the Fourier Domain" (PDF). Magn Reson Med. 58 (2): 430–436. doi:10.1002/mrm.21292. PMID 17616967. S2CID 16875383. Archived from the original (PDF) on 2013-02-16.
  7. ^ Report of the Task Group on Reference Man: ICRP Publication 23.
  8. ^ 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.
  9. ^ 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.
  10. ^ 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.
  11. ^ 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.
  12. ^ Gibbs, S. and Pujol, J. (1982). "A Monte Carlo method for patient dosimetry from diagnostic x-ray." Dentomaxillofac Radiol. 11, p. 25.
  13. ^ 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.
  14. ^ Zubal, I.G. et al. (1994). "Computerized three-dimensional segmented human anatomy." Med Phys, 21, p. 299.
  15. ^ LaRosa, Mary. "the Zubal Phantom". Retrieved 2015-08-17.
  16. ^ 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.
  17. ^ Lee, C. et al. (2006). "Whole-body voxel phantoms of paediatric patients—UF Series B." Phys Med Biol, 51, p. 4649.
  18. ^ A. Christ, W. Kainz, E. G. Hahn, K. Honegger, M. Zefferer, E. Neufeld, W. Rascher, R. Janka, W. Bautz, J. Chen, B. Kiefer, P. Schmitt, H. Hollenbach, J. Shen, M. Oberle, D. Szczerba, A. Kam, J. W. Guag, and N. Kuster. The virtual family development of surface based anatomical models of two adults and two children for dosimetric simulations. Phys. Med. Bio., 55(2):23–38, 2010.
  19. ^ S. Gu, R. Gupta and I. Kyprianou, Computational high-resolution heart phantoms for medical imaging and dosimetry simulations, Phys. Med. Biol., 56, no. 18, (2011): 5845-5864. [1]
  20. ^ 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.
  21. ^ 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.
  22. ^ 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.
  23. ^ 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.
  24. ^ 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.
  25. ^ 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.
  26. ^ 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.
  27. ^ 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 (3): 253–272. doi:10.1097/00004032-200403000-00003. PMID 14982227. S2CID 31201029.
  28. ^ Petoussi-Henss, N.; Zankl, M.; Fill, U.; Regulla, D. (2002). "The GSF family of voxel phantoms". Phys. Med. Biol. 47 (1): 89–106. Bibcode:2002PMB....47...89P. doi:10.1088/0031-9155/47/1/307. PMID 11814230. S2CID 250862634.
  29. ^ Zankl, Maria (2010). "The GSF voxel computational phantom family". Handbook of Anatomical Models for Radiation Dosimetry: 65–85.
  30. ^ 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 (2): 153–162. doi:10.1007/s004110100094. PMID 11484787. S2CID 29684856.
  31. ^ Zankl, M.; Wittmann, A. (2001). "The adult male voxel model "Golem" segmented from whole body CT patient data". Radiat. Environ. Biophys. 40 (2): 153–162. doi:10.1007/s004110100094. PMID 11484787. S2CID 29684856.
  32. ^ Mofrad F. B; et al. (2010). "Statistical construction of a Japanese male liver phantom for internal radionuclide dosimetry". Radiat Prot Dosimetry. 140 (2): 140–148. doi:10.1093/rpd/ncq164. PMID 20562118.
  33. ^ a b c 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 Polytechnic Institute
  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. ^ a b 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. ^ 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.
  41. ^ 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.
  42. ^ Lee C (2010). "The UF family of reference hybrid phantoms for computational radiation dosimetry". Phys. Med. Biol. 55 (2): 339–363. Bibcode:2010PMB....55..339L. doi:10.1088/0031-9155/55/2/002. PMC 2800036. PMID 20019401.
  43. ^ a b Stabin, M. et al. ICRP-89 based adult and pediatric phantom series, J NUCL MED MEETING ABSTRACTS, 49, 14, 2008.
  44. ^ Cassola V., Lima V., Kramer R., Khoury H. (2010). "FASH and MASH: female and male adult human phantoms based on polygon mesh surfaces: I. Development of the anatomy". Phys. Med. Biol. 55 (133): 133–162. Bibcode:2010PMB....55..133C. doi:10.1088/0031-9155/55/1/009. PMID 20009183. S2CID 8506045.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  45. ^ 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).
  46. ^ Kim CH (2010). "A polygon-surface reference Korean male phantom (PSRK-Man) and its direct implementation in Geant4 Monte Carlo simulation". Phys. Med. Biol. 56 (10): 3137–3161. doi:10.1088/0031-9155/56/10/016. PMID 21521906. S2CID 42132546.
  47. ^ J.A. Vazquez. (2012). The Need and Feasibility of a Dynamic Computational Human Phantom for Radiation Dosimetry Simulation Using Motion Capture Data (Master’s thesis). Troy, NY: Rensselaer Polytechnic Institute

External links[edit]