Evaluation of radiation dose to anthropomorphic paediatric models from positron-emitting labelled tracers

Evaluation of radiation dose to pediatric models from whole body PET/CT imaging

Positron emission tomography (PET)/computed tomography (CT) is a well-known modality for the diagnosis of various diseases in children and adult patients. On the other hand, younger patients are more radiosensitive than adults, so there are concerns about the level of ionizing radiation exposure in pediatric whole body PET/CT imaging. In this regard, comprehensive specific radiation dosimetry for whole body PET/CT imaging is highly desired for different ages, sizes, and shapes. Therefore, in this study, organ absorbed doses were evaluated for pediatric voxel models from 4 to 14 years old and compared with those of ICRP phantoms. Monte Carlo calculation was performed to evaluate S-value, absorbed dose, and effective dose from ¹⁸ F-FDG radiotracers and whole body CT scan for different computational models, including 4- to 14-year-old phantoms. The results showed that the S-value and, therefore, absorbed dose of ¹⁸ F-FDG strongly depended on the phantom anatomy. These variations were justified by the distance between source and target organs. Moreover, on average, the absorbed doses from whole body CT scans were 13.5 times lower than those from ¹⁸ F-FDG for all organs. According to the results, various anatomies and ages should be considered for accurate dose evaluation. These data can be used for specific risk assessment of the pediatric population in clinical nuclear imaging.

Development of an age-scalable 3D computational phantom in DICOM standard for late effects studies of childhood cancer survivors

Purpose

We previously developed an age-scalable 3D computational phantom that has been widely used for retrospective whole-body dose reconstructions of conventional two-dimensional historic radiation therapy (RT) treatments in late effects studies of childhood cancer survivors. This phantom is modeled in the FORTRAN programming language and is not readily applicable for dose reconstructions for survivors treated with contemporary RT whose treatment plans were designed using computed tomography images and complex treatment fields. The goal of this work was to adapt the current FORTRAN model of our age-scalable computational phantom into Digital Imaging and Communications in Medicine (DICOM) standard so that it can be used with any treatment planning system (TPS) to reconstruct contemporary RT. Additionally, we report a detailed description of the phantom's age-based scaling functions, information that was not previously published.

Method

We developed a Python script that adapts our phantom model from FORTRAN to DICOM. To validate the conversion, we compared geometric parameters for the phantom modeled in FORTRAN and DICOM scaled to ages 1 month, 6 months, 1, 2, 3, 5, 8, 10, 15, and 18 years. Specifically, we calculated the percent differences between the corner points and volume of each body region and the normalized mean square distance (NMSD) between each of the organs. In addition, we also calculated the percent difference between the heights of our DICOM age-scaled phantom and the heights (50th percentile) reported by the World Health Organization (WHO) and Centers for Disease Control and Prevention (CDC) for male and female children of the same ages. Additionally, we calculated the difference between the organ masses for our DICOM phantom and the organ masses for two reference phantoms (from International Comission on Radiation Protection (ICRP) 89 and the University of Florida/National Cancer Institute reference hybrid voxel phantoms) for ages newborn, 1, 5, 10, 15 and adult. Lastly, we conducted a feasibility study using our DICOM phantom for organ dose calculations in a commercial TPS. Specifically, we simulated a 6 MV photon right-sided flank field RT plan for our DICOM phantom scaled to age 3.9 years; treatment field parameters and age were typical of a Wilms tumor RT treatment in the Childhood Cancer Survivor Study. For comparison, the same treatment was simulated using our in-house dose calculation system with our FORTRAN phantom. The percent differences (between FORTRAN and DICOM) in mean dose and percent of volume receiving dose ⩾5 Gy were calculated for two organs at risk, liver and pancreas.

Results

The percent differences in corner points and the volumes of head, neck, and trunk body regions between our phantom modeled in FORTRAN and DICOM agreed within 3%. For all of the ages, the NMSDs were negliglible with a maximum NMSD of 7.80 × 10-2 mm for occiptital lobe of 1 month. The heights of our age-scaled phantom agreed with WHO/CDC data within 7% from infant to adult, and within 2% agreement for ages 5 years and older. We observed that organ masses in our phantom are less than the organ masses for other reference phantoms. Dose calculations done with our in-house calculation system (with FORTRAN phantom) and commercial TPS (with DICOM phantom) agreed within 7%.

Conclusion

We successfully adapted our phantom model from the FORTRAN language to DICOM standard and validated its geometric consistency. We also demonstrated that our phantom model is representative of population height data for infant to adult, but that the organ masses are smaller than in other reference phantoms and need further refinement. Our age-scalable computational phantom modeled in DICOM standard can be scaled to any age at RT and used within a commercial TPS to retrospectively reconstruct doses from contemporary RT in childhood cancer survivors.

Innovations in Computer Technologies Have Impacted Radiation Dosimetry Through Anatomically Realistic Phantoms and Fast Monte Carlo Simulations

Radiological physics principles have not changed in the past 60 y when computer technologies advanced exponentially. The research field of anatomical modeling for the purpose of radiation dose calculations has experienced an explosion in activity in the past two decades. Such an exciting advancement is due to the feasibility of creating three-dimensional geometric details of the human anatomy from tomographic imaging and of performing Monte Carlo radiation transport simulations on increasingly fast and cheap personal computers. The advent of a new type of high-performance computing hardware in recent years-graphics processing units-has made it feasible to carry out time-consuming Monte Carlo calculations at near real-time speeds. This paper introduces the history of three generations of computational human phantoms (the stylized medical internal radiation dosimetry-type phantoms, the voxelized tomographic phantoms, and the boundary representation deformable phantoms) and new development of the graphics processing unit-based Monte Carlo radiation dose calculations. Examples are given for research projects performed by my students in applying computational phantoms and a new Monte Carlo code, ARCHER, to problems in radiation protection, imaging, and radiotherapy. Finally, the paper discusses challenges and future opportunities for research.

Re-evaluation of pediatric 18F-FDG dosimetry: Cristy-Eckerman versus UF/NCI hybrid computational phantoms

Because of the concerns associated with radiation exposure at a young age, there is an increased interest in pediatric absorbed dose estimates for imaging agents. Almost all reported pediatric absorbed dose estimates, however, have been determined using adult pharmacokinetic data with radionuclide S values that take into account the anatomical differences between adults and children based upon the older Cristy-Eckerman (C-E) stylized phantoms. In this work, we use pediatric model-derived pharmacokinetics to compare absorbed dose and effective dose estimates for ¹⁸F-FDG in pediatric patients using S values generated from two different geometries of computational phantoms. Time-integrated activity coefficients of ¹⁸F-FDG in brain, lungs, heart wall, kidneys and liver, retrospectively, calculated from 35 pediatric patients at the Boston's Children Hospital were used. The absorbed dose calculation was performed in accordance with the Medical Internal Radiation Dose method using S values generated from the University of Florida/National Cancer Institute (UF/NCI) hybrid phantoms, as well as those from C-E stylized computational phantoms. The effective dose was computed using tissue-weighting factors from ICRP Publication 60 and ICRP Publication 103 for the C-E and UF/NCI, respectively. Substantial differences in the absorbed dose estimates between UF/NCI hybrid pediatric phantoms and the C-E stylized phantoms were found for the lungs, ovaries, red bone marrow and urinary bladder wall. Large discrepancies in the calculated dose values were observed in the bone marrow; ranging between  -26% to  +199%. The effective doses computed by the UF/NCI hybrid phantom S values were slightly different than those seen using the C-E stylized phantoms with percent differences of  -0.7%, 2.9% and  2.5% for a newborn, 1 year old and 5 year old, respectively. Differences in anatomical modeling features among computational phantoms used to perform Monte Carlo-based photon and electron transport simulations for ¹⁸F, and very likely for other radionuclides, impact internal organ dosimetry computations for pediatric nuclear medicine studies.

Effects of body habitus on internal radiation dose calculations using the 5-year-old anthropomorphic male models

Computational phantoms are commonly used in internal radiation dosimetry to assess the amount and distribution pattern of energy deposited in various parts of the human body from different internal radiation sources. Radiation dose assessments are commonly performed on predetermined reference computational phantoms while the argument for individualized patient-specific radiation dosimetry exists. This study aims to evaluate the influence of body habitus on internal dosimetry and to quantify the uncertainties in dose estimation correlated with the use of fixed reference models. The 5-year-old IT'IS male phantom was modified to match target anthropometric parameters, including body weight, body height and sitting height/stature ratio (SSR), determined from reference databases, thus enabling the creation of 125 5-year-old habitus-dependent male phantoms with 10th, 25th, 50th, 75th and 90th percentile body morphometries. We evaluated the absorbed fractions and the mean absorbed dose to the target region per unit cumulative activity in the source region (S-values) of F-18 in 46 source regions for the generated 125 anthropomorphic 5-year-old hybrid male phantoms using the Monte Carlo N-Particle eXtended general purpose Monte Carlo transport code and calculated the absorbed dose and effective dose of five ¹⁸F-labelled radiotracers for children of various habitus. For most organs, the S-value of F-18 presents stronger statistical correlations with body weight, standing height and sitting height than BMI and SSR. The self-absorbed fraction and self-absorbed S-values of F-18 and the absorbed dose and effective dose of ¹⁸F-labelled radiotracers present with the strongest statistical correlations with body weight. For ¹⁸F-Amino acids, ¹⁸F-Brain receptor substances, ¹⁸F-FDG, ¹⁸F-L-DOPA and ¹⁸F-FBPA, the mean absolute effective dose differences between phantoms of different habitus and fixed reference models are 11.4%, 11.3%, 10.8%, 13.3% and 11.4%, respectively. Total body weight, standing height and sitting height have considerable effects on human internal dosimetry. Radiation dose calculations for individual subjects using the most closely matched habitus-dependent computational phantom should be considered as an alternative to improve the accuracy of the estimates.

Development of computational pregnant female and fetus models and assessment of radiation dose from positron-emitting tracers

PURPOSE

Molecular imaging using PET and hybrid (PET/CT and PET/MR) modalities nowadays plays a pivotal role in the clinical setting for diagnosis and staging, treatment response monitoring, and radiation therapy treatment planning of a wide range of oncologic malignancies. The developing embryo/fetus presents a high sensitivity to ionizing radiation. Therefore, estimation of the radiation dose delivered to the embryo/fetus and pregnant patients from PET examinations to assess potential radiation risks is highly praised.

METHODS

We constructed eight embryo/fetus models at various gestation periods with 25 identified tissues according to reference data recommended by the ICRP publication 89 representing the anatomy of the developing embryo/fetus. The developed embryo/fetus models were integrated into realistic anthropomorphic computational phantoms of the pregnant female and used for estimating, using Monte Carlo calculations, S-values of common positron-emitting radionuclides, organ absorbed dose, and effective dose of a number of positron-emitting labeled radiotracers.

RESULTS

The absorbed dose is nonuniformly distributed in the fetus. The absorbed dose of the kidney and liver of the 8-week-old fetus are about 47.45 % and 44.76 % higher than the average absorbed dose of the fetal total body for all investigated radiotracers. For ¹⁸F-FDG, the fetal effective doses are 2.90E-02, 3.09E-02, 1.79E-02, 1.59E-02, 1.47E-02, 1.40E-02, 1.37E-02, and 1.27E-02 mSv/MBq at the 8th, 10th, 15th, 20th, 25th, 30th, 35th, and 38th weeks of gestation, respectively.

CONCLUSION

The developed pregnant female/fetus models matching the ICRP reference data can be exploited by dedicated software packages for internal and external dose calculations. The generated S-values will be useful to produce new standardized dose estimates to pregnant patients and embryo/fetus from a variety of positron-emitting labeled radiotracers.

Assessment of radiation dose in nuclear cardiovascular imaging using realistic computational models

PURPOSE

Nuclear cardiology plays an important role in clinical assessment and has enormous impact on the management of a variety of cardiovascular diseases. Pediatric patients at different age groups are exposed to a spectrum of radiation dose levels and associated cancer risks different from those of adults in diagnostic nuclear medicine procedures. Therefore, comprehensive radiation dosimetry evaluations for commonly used myocardial perfusion imaging (MPI) and viability radiotracers in target population (children and adults) at different age groups are highly desired.

METHODS

Using Monte Carlo calculations and biological effects of ionizing radiation VII model, we calculate the S-values for a number of radionuclides (Tl-201, Tc-99m, I-123, C-11, N-13, O-15, F-18, and Rb-82) and estimate the absorbed dose and effective dose for 12 MPI radiotracers in computational models including the newborn, 1-, 5-, 10-, 15-yr-old, and adult male and female computational phantoms.

RESULTS

For most organs, (201)Tl produces the highest absorbed dose whereas (82)Rb and (15)O-water produce the lowest absorbed dose. For the newborn baby and adult patient, the effective dose of (82)Rb is 48% and 77% lower than that of (99m)Tc-tetrofosmin (rest), respectively.

CONCLUSIONS

(82)Rb results in lower effective dose in adults compared to (99m)Tc-labeled tracers. However, this advantage is less apparent in children. The produced dosimetric databases for various radiotracers used in cardiovascular imaging, using new generation of computational models, can be used for risk-benefit assessment of a spectrum of patient population in clinical nuclear cardiology practice.

An exponential growth of computational phantom research in radiation protection, imaging, and radiotherapy: a review of the fifty-year history

Radiation dose calculation using models of the human anatomy has been a subject of great interest to radiation protection, medical imaging, and radiotherapy. However, early pioneers of this field did not foresee the exponential growth of research activity as observed today. This review article walks the reader through the history of the research and development in this field of study which started some 50 years ago. This review identifies a clear progression of computational phantom complexity which can be denoted by three distinct generations. The first generation of stylized phantoms, representing a grouping of less than dozen models, was initially developed in the 1960s at Oak Ridge National Laboratory to calculate internal doses from nuclear medicine procedures. Despite their anatomical simplicity, these computational phantoms were the best tools available at the time for internal/external dosimetry, image evaluation, and treatment dose evaluations. A second generation of a large number of voxelized phantoms arose rapidly in the late 1980s as a result of the increased availability of tomographic medical imaging and computers. Surprisingly, the last decade saw the emergence of the third generation of phantoms which are based on advanced geometries called boundary representation (BREP) in the form of Non-Uniform Rational B-Splines (NURBS) or polygonal meshes. This new class of phantoms now consists of over 287 models including those used for non-ionizing radiation applications. This review article aims to provide the reader with a general understanding of how the field of computational phantoms came about and the technical challenges it faced at different times. This goal is achieved by defining basic geometry modeling techniques and by analyzing selected phantoms in terms of geometrical features and dosimetric problems to be solved. The rich historical information is summarized in four tables that are aided by highlights in the text on how some of the most well-known phantoms were developed and used in practice. Some of the information covered in this review has not been previously reported, for example, the CAM and CAF phantoms developed in 1970s for space radiation applications. The author also clarifies confusion about 'population-average' prospective dosimetry needed for radiological protection under the current ICRP radiation protection system and 'individualized' retrospective dosimetry often performed for medical physics studies. To illustrate the impact of computational phantoms, a section of this article is devoted to examples from the author's own research group. Finally the author explains an unexpected finding during the course of preparing for this article that the phantoms from the past 50 years followed a pattern of exponential growth. The review ends on a brief discussion of future research needs (a supplementary file '3DPhantoms.pdf' to figure 15 is available for download that will allow a reader to interactively visualize the phantoms in 3D).