Voxel-based Phantom Research Articles

Voxel-Based Frog Phantom for Internal Dose Evaluation

A voxel-based frog phantom has been developed for radiation protection of the environment. The voxel-based frog phantom was applied to evaluating self-absorbed fractions (self-AFs), which are defined as the fraction of energy emitted by a radiation source that is absorbed within the source organ. The self-AFs were evaluated for both photons and electrons in the spleen, kidneys, and liver using Monte Carlo simulations. Furthermore, self-S values (mGy/MBq.s) for 18F and 90Y in the organs were calculated using the results of the self-AFs. Consequently, the voxel-based frog phantom was found to be useful for the organ dose evaluations, which have not been proposed by the International Commission on Radiological Protection (ICRP). It was also confirmed that the self-AFs and self-S values are largely dependent on the mass of the source organ.

Calculation of Dose Contributions of Electron and Charged Heavy Particles inside Phantoms Irradiated by Monoenergetic Neutron

The radiation-transport code PHITS with an event generator mode has been applied to analyze energy depositions of electrons and charged heavy particles in two spherical phantoms and a voxel-based mouse phantom upon neutron irradiation. The calculations using the spherical phantoms quantitatively clarified the type and energy of charged particles which are released through interactions of neutrons with the phantom elements and contribute to the radiation dose. The relative contribution of electrons increased with an increase in the size of the phantom and with a decrease in the energy of the incident neutrons. Calculations with the voxel-based mouse phantom for 2.0-MeV neutron irradiation revealed that the doses to different locations inside the body are uniform, and that the energy is mainly deposited by recoil protons. The present study has demonstrated that analysis using PHITS can yield dose distributions that are accurate enough for RBE evaluation.

Effect of Radiotherapy Volume and Dose on Secondary Cancer Risk in Stage I Testicular Seminoma

To estimate and compare the secondary cancer risk (SCR) due to para-aortic (PA), dogleg field (DLF), or extensive field (EF) radiotherapy (RT) at different dose levels for Stage I testicular seminoma. The organ equivalent dose concept with a linear, plateau, and linear-exponential dose-response model was applied to the dose distributions to estimate the SCR. The dose distributions were calculated in a voxel-based anthropomorphic phantom. Three different three-dimensional plans were computed: PA, DLF, and EF. The plans were calculated with 6-MV photons and two opposed fields, using 20 Gy in 10 fractions. The estimated cumulative SCR for a 75-year-old patient treated with PA-RT at age 35 was 23.3% (linear model), 20.9% (plateau model), and 20.8% (linear-exponential model) compared with 19.8% for the general population. Dependent on the model, PA-RT compared with DLF-RT reduced the SCR by 48-63% or 64-69% when normalized to EF-RT. For PA-RT, the linear dose-response model predicted a decrease of 45% in the SCR, using 20 Gy instead of 30 Gy; the linear-exponential dose-response model predicted no change in SCR. Our model suggested that the SCR after PA-RT for Stage I testicular seminoma is reduced by approximately one-half to two-thirds compared with DLF-RT, independent of the dose-response model. The SCR is expected to be equal or lower with 20 Gy than with 30 Gy. In the absence of mature patient data, the organ equivalent dose concept offers the best potential method of estimating the SCR when discussing treatment options with patients.

Novel approach to stationary transmission scanning using Compton scattered radiation

Transmission scanning-based estimation of the attenuation map plays a crucial role in quantitative radionuclide imaging. X-ray computed tomography (CT) reconstructs directly the attenuation coefficients map from data transmitted through the object. This paper proposes an alternative route for reconstructing the object attenuation map by exploiting Compton scatter of transmitted radiation from an externally placed radionuclide source. In contrast to conventional procedures, data acquisition is realized as a series of images parameterized by the Compton scattering angle and registered on a stationary gamma camera operating without spatial displacement. Numerical simulation results using realistic voxel-based phantoms are presented to illustrate the efficiency of this new transmission scanning approach for attenuation map reconstruction. The encouraging results presented in this paper may suggest the possibility of proposing a new concept for emission/transmission imaging using scattered radiation, which has many advantages compared to conventional technologies.

Voxel-based models representing the male and female ICRP reference adult--the skeleton

For the forthcoming update of organ dose conversion coefficients, the International Commission on Radiological Protection (ICRP) will use voxel-based computational phantoms due to their improved anatomical realism compared with the class of mathematical or stylized phantoms used previously. According to the ICRP philosophy, these phantoms should be representative of the male and female reference adults with respect to their external dimensions, their organ topology and their organ masses. To meet these requirements, reference models of an adult male and adult female have been constructed at the GSF, based on existing voxel models segmented from tomographic images of two individuals whose body height and weight closely resemble the ICRP Publication 89 reference values. The skeleton is a highly complex structure of the body, composed of cortical bone, trabecular bone, red and yellow bone marrow and endosteum ('bone surfaces' in their older terminology). The skeleton of the reference phantoms consists of 19 individually segmented bones and bone groups. Sub-division of these bones into the above-mentioned constituents would be necessary in order to allow a direct calculation of dose to red bone marrow and endosteum. However, the dimensions of the trabeculae, the cavities containing bone marrow and the endosteum layer lining these cavities are clearly smaller than the resolution of a normal CT scan and, thus, these volumes could not be segmented in the tomographic images. As an attempt to represent the gross spatial distribution of these regions as realistically as possible at the given voxel resolution, 48 individual organ identification numbers were assigned to various parts of the skeleton: every segmented bone was subdivided into an outer shell of cortical bone and a spongious core; in the shafts of the long bones, a medullary cavity was additionally segmented. Using the data from ICRP Publication 89 on elemental tissue composition, from ICRU Report 46 on material mass densities, and from ICRP Publication 70 on the distribution of the red bone marrow among and marrow cellularity in individual bones, individual elemental compositions for these segmented bone regions were derived. Thus, most of the relevant source and target regions of the skeleton were provided. Dose calculations using these regions will be based on fluence-to-dose response functions that are multiplied with the particle fluence inside specific bone regions to give the dose quantities of interest to the target tissues.

SEECAL utilizing voxel-based SAFs

The computer program SEECAL, written by Cristy and Eckerman from the Oak Ridge National Laboratory, USA, calculates specific effective energies (also known as S factors in the MIRD terminology) for an adult male, an adult female and five paediatric ages. Its dosimetric methodology is that of the ICRP. Among other parameters, SEECAL requires input data on specific absorbed fractions (SAF) and utilises those derived from the MIRD-type stylised anthropomorphic phantoms. SEECAL has been used worldwide for dose estimations concerning occupational or public exposures due to radionuclides incorporated into the body and has formed the basis for programs developed by other laboratories to calculate, for example, dose to the patients undergoing nuclear medicine procedures. The revised version of SEECAL is at the moment limited to adults and utilises the photon SAFs derived with Monte Carlo methods for the new reference male and female voxel-based phantoms to be adopted by the ICRP.

Mouse S-Factors Based on Monte Carlo Simulations in the Anatomical Realistic Moby Phantom for Internal Dosimetry

Biokinetic and dosimetry studies in small animals often precede clinical radionuclide therapies. As in human studies, a reliable evaluation of therapeutic efficacy is essential and must be based on accurate dosimetry, which must be based on a realistic dosimetry model. The aim of this study was to evaluate the differences in the results when using a more anatomic realistic mouse phantom, as compared to previously mathematically described phantoms, based mainly on ellipsoids and cylinders. The difference in results from the two Monte Carlo codes, EGS4 and MCNPX 2.6a, was also evaluated. An anatomical correct mouse phantom (Moby) was developed by Segars et al. for the evaluation and optimization of the in vivo imaging of mice. The Moby phantom is based on surfaces, which allows for an easy and flexible definition of organ sizes. It includes respiratory movements and a beating heart. It also allows for a redefinition of the location of several internal organs. The execution of the Moby program generates a three-dimensional voxel-based phantom of a specified size, which was modified and used as input for Monte Carlo simulations of absorbed fractions and S-factors. The radiation transport was simulated both with the EGS4 system and the MCNPX 2.6a code. Calculations were done for the radionuclides (18)F, (124)I, (131)I, (111)In, (177)Lu, and (90)Y. S-factors were calculated using in-house-developed IDL programs and compared with results from previously published models. The comparison of S-factors obtained by the Moby model and mathematical phantoms showed that these, in many cases, were within the same range, whereas for some organs, they were underestimated in the mathematical phantoms. The results were closer to the more anatomically realistic phantom than to the mathematical phantoms, with some exceptions. When investing differences between MCNPX 2.6a and EGS4 using the Moby phantom, results indicated some differences in absorbed fractions for electrons. This reason may be owing to differences in the codes regarding the theory for which electron transport are simulated. It is possible to calculate S-factors that are specific for small animals, such as mice. The Moby phantom is useful as a dosimetry model because it is anatomically realistic, but still very flexible, with 35 accurately segmented regions.

Organ and effective doses in pediatric patients undergoing helical multislice computed tomography examination

As multidetector computed tomography (CT) serves as an increasingly frequent diagnostic modality, radiation risks to patients became a greater concern, especially for children due to their inherently higher radiosensitivity to stochastic radiation damage. Current dose evaluation protocols include the computed tomography dose index (CTDI) or point detector measurements using anthropomorphic phantoms that do not sufficiently provide accurate information of the organ-averaged absorbed dose and corresponding effective dose to pediatric patients. In this study, organ and effective doses to pediatric patients under helical multislice computed tomography (MSCT) examinations were evaluated using an extensive series of anthropomorphic computational phantoms and Monte Carlo radiation transport simulations. Ten pediatric phantoms, five stylized (equation-based) ORNL phantoms (newborn, 1-year, 5-year, 10-year, and 15-year) and five tomographic (voxel-based) UF phantoms (9-month male, 4-year female, 8-year female, 11-year male, and 14-year male) were implemented into MCNPX for simulation, where a source subroutine was written to explicitly simulate the helical motion of the CT x-ray source and the fan beam angle and collimator width. Ionization chamber measurements were performed and used to normalize the Monte Carlo simulation results. On average, for the same tube current setting, a tube potential of 100 kVp resulted in effective doses that were 105% higher than seen at 80 kVp, and 210% higher at 120 kVp regardless of phantom type. Overall, the ORNL phantom series was shown to yield values of effective dose that were reasonably consistent with those of the gender-specific UF phantom series for CT examinations of the head, pelvis, and torso. However, the ORNL phantoms consistently overestimated values of the effective dose as seen in the UF phantom for MSCT scans of the chest, and underestimated values of the effective dose for abdominal CT scans. These discrepancies increased with increasing kVp. Finally, absorbed doses to select radiation sensitive organs such as the gonads, red bone marrow, colon, and thyroid were evaluated and compared between phantom types. Specific anatomical problems identified in the stylized phantoms included excessive pelvic shielding of the ovaries in the female phantoms, enhanced red bone marrow dose to the arms and rib cage for chest exams, an unrealistic and constant torso thickness resulting in excessive x-ray attenuation in the regions of the abdominal organs, and incorrect positioning of the thyroid within the stylized phantom neck resulting in insufficient shielding by clavicles and scapulae for lateral beam angles. To ensure more accurate estimates of organ absorbed dose in multislice CT, it is recommended that voxel-based phantoms, potentially tailored to individual body morphometry, be utilized in any future prospective epidemiological studies of medically exposed children.

Calculation of the photonuclear yield using an anthropomorphic phantom

The photoneutron production from different elements in major body organs using an anthropomorphic voxel-based phantom irradiated with 15, 18 and 25 MV photon beam spectra from Medical Linear Accelerator were estimated by the mean of MCNPX Monte Carlo code. The soft tissue composition recommended by ICRU44 was used in the simulations. Other elements present, less than 0.005% by weight collectively were neglected. The results of this initial simulation will provide the basis for further simulation to estimate the changes in the elemental concentration in the irradiated volume.

Monte Carlo simulations of absorbed dose in a mouse phantom from 18-fluorine compounds.

The purpose of this study was to calculate internal absorbed dose distribution in mice from pre-clinical small animal PET imaging procedures with fluorine-18 labeled compounds (18FDG, 18FLT, and fluoride ion). The GATE Monte Carlo software and a realistic, voxel-based mouse phantom that included a subcutaneous tumor were used to perform simulations. Discretized time-activity curves obtained from dynamic in vivo studies with each of the compounds were used to set the activity concentration in the simulations. For 18FDG, a realistic range of uptake ratios was considered for the heart and tumor. For each simulated time frame, the biodistribution of the radionuclide in the phantom was considered constant, and a sufficient number of decays were simulated to achieve low statistical uncertainty. Absorbed dose, which was scaled to take into account radioactive decay, integration with time, and changes in biological distribution was reported in mGy per MBq of administered activity for several organs and uptake scenarios. The mean absorbed dose ranged from a few mGy/MBq to hundreds of mGy/MBq. Major organs receive an absorbed dose in a range for which biological effects have been reported. The effects on a given investigation are hard to predict; however, investigators should be aware of potential perturbations especially when the studied organ receives high absorbed dose and when longitudinal imaging protocols are considered.

A high-resolution anthropomorphic voxel-based tomographic phantom for proton therapy of the eye

Proton therapy is increasingly used in medical treatments for cancer patients due to the sharp dose conformity offered by the characteristic Bragg peak. Proton beam interactions with the eye will be simulated using the MCNPX Monte Carlo code and available nuclear cross-section data to calculate the dose distribution in the eye gel and surrounding organs. A high-resolution eye model will be employed using a 3D geometrical voxel-based anthropomorphic head phantom obtained from the Visible Human Project (female data). Manual segmentation of the eye, carried out by the Medical Physics group at the University of Surrey resulted in 15 identified structures. This work emphasizes the use of a realistic phantom for accurately predicting dose deposition by protons.

Age-dependent organ and effective dose coefficients for external photons: a comparison of stylized and voxel-based paediatric phantoms

This present study investigates the anatomical realism of conventional stylized models of children by comparing organ dose conversion coefficients for the ORNL paediatric phantom series with those determined in the UF (University of Florida) voxel paediatric phantoms. The latter includes whole-body models of a 9 month male, 4 year female, 8 year female, 11 year male and a 14 year male. Of these phantoms, the 1 year, 5 year and 10 year ORNL phantoms, and 9 month male, 4 year female and 11 year male UF voxel phantoms were selected for side-by-side comparisons under idealized external photon irradiation. Organ absorbed dose per unit air kerma (Gy/Gy) for various radiosensitive organs and tissues were calculated for monoenergetic photons over the energy range of 15 keV to 10 MeV and for six irradiation geometries: anterior–posterior (AP), posterior–anterior (PA), right lateral (RLAT), left lateral (LLAT), rotational (ROT) and isotropic (ISO). Differences in organ dose conversion coefficients for the gonads, bone marrow, colon, lung and stomach, to which prominent tissue weighting factors are assigned, were depicted and analysed. Two major causes of observed differences were suggested: differences in organ shape and position and the differences in tissue shielding by overlying tissue regions within the phantoms. Significant discrepancies caused by anatomical differences between the two types of phantoms are also reported for several organs, and in particular, the thyroid and urinary bladder. The results of this study suggest that the paediatric series of ORNL phantoms also have less realistic internal organ and body anatomy and that dose conversion coefficients from these stylized phantoms should be re-evaluated using paediatric voxel phantoms.

Comparison between effective doses for voxel-based and stylized exposure models from photon and electron irradiation

For the last two decades, the organ and tissue equivalent dose as well as effective dose conversion coefficients recommended by the International Commission on Radiological Protection (ICRP) have been determined with exposure models based on stylized MIRD5-type phantoms representing the human body with its radiosensitive organs and tissues according to the ICRP Reference Man released in Publication No. 23, on Monte Carlo codes sometimes simulating rather simplified radiation physics and on tissue compositions from different sources. Meanwhile the International Commission on Radiation Units and Measurements (ICRU) has published reference data for human tissue compositions in Publication No. 44, and the ICRP has released a new report on anatomical reference data in Publication No. 89. As a consequence many of the components of the traditional stylized exposure models used to determine the effective dose in the past have to be replaced: Monte Carlo codes, human phantoms and tissue compositions. This paper presents results of comprehensive investigations on the dosimetric consequences to be expected from the replacement of the traditional stylized exposure models by the voxel-based exposure models. Calculations have been performed with the EGS4 Monte Carlo code for external and internal exposures to photons and electrons with the stylized, gender-specific MIRD5-type phantoms ADAM and EVA on the one hand and with the recently developed tomographic or voxel-based phantoms MAX and FAX on the other hand for a variety of exposure conditions. Ratios of effective doses for the voxel-based and the stylized exposure models will be presented for external and internal exposures to photons and electrons as a function of the energy and the geometry of the radiation field. The data indicate that for the exposure conditions considered in these investigations the effective dose may change between +60% and −50% after the replacement of the traditional exposure models by the voxel-based exposure models.

Characterization and correction of temporal artifacts in CT.

Image artifacts, caused by a temporally delayed response of a computed tomography (CT) detector, were investigated. To study its consequences, a computer model for a standard third generation CT scanner was devised and simulations were carried out. Resulting image artifacts were studied for various voxel-based and mathematical phantoms using three time constants, which characterize the delayed signal responses of different detector configurations. Furthermore, a theory was developed to understand temporal artifacts in reconstructed images. A filter function was derived, which compensates for temporal artifacts. For a given phantom, simulations and theory demonstrate that artifacts scale with the time constant, characteristic for the detector, and the angular speed of the scanner.

OEDIPE: A Personalized Dosimetric Tool Associating Voxel-Based Models with MCNPX

A new tool, named OEDIPE (a French acronym that stands for "Tool for Personalized Internal Dose Assessment") was developed to carry out personalized internal dosimetry calculations for nuclear medicine (for both therapeutic and diagnostic procedures) and for radiation safety (in the case of internal contamination). It was developed under the PV-Wave visual data analysis system by the Institute of Radioprotection and Nuclear Safety (IRSN) in collaboration with the French Institute of Health and Medical Research (INSERM). This software creates anthropomorphic voxel-based phantoms from computed tomography (CT) and magnetic resonance imaging (MRI) patient images through the use of a friendly graphical user interface (GUI). Several tools have been built-in to allow for image segmentation. Source data, including VOI localization and cumulated activities, are assessed by single photon emission computed tomography (SPECT) images, and the source may be specified in any number of organs either as a point source or a homogeneously distributed source. It is also possible to choose the dosimetric parameters required for the study (mean organ dose or a dose distribution). Phantom, source, and dosimetric parameters are automatically written into a file. That file is then processed by the Monte Carlo code MCNPX (LANL) to perform the actual dose calculation. OEDIPE can compute either the absorbed dose in each organ (in a few minutes), or the absorbed dose in each voxel of the phantom (i.e. the spatial dose distribution at a tissue level) in a few hours or more. OEDIPE automatically reads the MCNPX output file and processes results to give a list of absorbed doses in each organ or a plot of isodose curves superimposed onto the phantom. Because of the long calculation times required to compute an absorbed dose within an entire whole-body phantom at a spatial resolution of a few millimeters, modifications were made to reduce computational times to reasonable values. To illustrate this tool, results of a dosimetric study of technetium-99m labeling of a bone-scanning agent are presented. OEDIPE is a tool that can be used for patient-specific dosimetry--for example, in targeted radiotherapy--by taking into account the individual patient anatomy, including tumors.

MAX meets ADAM: a dosimetric comparison between a voxel-based and a mathematical model for external exposure to photons

The International Commission on Radiological Protection intends to revise the organ and tissue equivalent dose conversion coefficients published in various reports. For this purpose the mathematical human medical internal radiation dose (MIRD) phantoms, actually in use, have to be replaced by recently developed voxel-based phantoms. This study investigates the dosimetric consequences, especially with respect to the effective male dose, if not only a MIRD phantom is replaced by a voxel phantom, but also if the tissue compositions and the radiation transport codes are changed. This task will be resolved by systematically replacing in the mathematical ADAM/GSF exposure model, first the radiation transport code, then the tissue composition and finally the phantom anatomy, in order to arrive at the voxel-based MAX/EGS4 exposure model. The results show that the combined effect of these replacements can decrease the effective male dose by up to 25% for external exposures to photons for incident energies above 30 keV for different field geometries, mainly because of increased shielding by a heterogeneous skeleton and by the overlying adipose and muscle tissue, and also because of the positions internal organs have in a realistically designed human body compared to their positions in the mathematically constructed phantom.

Photon specific absorbed fractions calculated in the trunk of an adult male voxel-based phantom.

A new approach for calculating internal dose estimates was developed through the use of a more realistic computational model of the human body. The study demonstrates the capability of building a patient-specific phantom with voxel-based data for the simulation of radiation transport and energy deposition using Monte Carlo methods such as the MCNP-4B code. MCNP-4B was used to calculate absorbed fractions for photons in a voxel-based phantom, and values were compared to reference values from traditional phantoms used for many years. Results obtained in general agreed well with previous values, but considerable differences were found in some cases due to two major causes; differences in the organ masses between the phantoms and the occurrence of organ overlap in the voxel-based phantom (which is not well modeled in the mathematical phantoms). These new techniques offer promise of developing a new generation of more realistic phantoms for internal, as well as external, dose assessment. The principal area of implementation in internal dose assessment should be the development of patient-specific dose estimates in nuclear medicine therapy, such as radioimmunotherapy (RIT). However, as new voxel-based phantoms for different individuals can be developed, they may also be used with the techniques developed here to derive new absorbed fractions and replace the traditional values usedfor other applications in internal and external dose assessment, which have been based on mathematical constructs that are not always very representative of real human organs.

Absorbed dose estimates to structures of the brain and head using a high-resolution voxel-based head phantom.

The purpose of this article is to demonstrate the viability of using a high-resolution 3-D head phantom in Monte Carlo N-Particle (MCNP) for boron neutron capture therapy (BNCT) structure dosimetry. This work describes a high-resolution voxel-based model of a human head and its use for calculating absorbed doses to the structures of the brain. The Zubal head phantom is a 3-D model of a human head that can be displayed and manipulated on a computer. Several changes were made to the original head phantom which now contains over 29 critical structures of the brain and head. The modified phantom is a 85 x 109 x 120 lattice of voxels, where each voxel is 2.2 x 2.2 x 1.4 mm3. This model was translated into MCNP lattice format. As a proof of principle study, two MCNP absorbed dose calculations were made (left and right lateral irradiations) using a uniformly distributed neutron disk source with an 1/E energy spectrum. Additionally, the results of these two calculations were combined to estimate the absorbed doses from a bilateral irradiation. Radiobiologically equivalent (RBE) doses were calculated for all structures and were normalized to 12.8 Gy-Eq. For a left lateral irradiation, the left motor cortex receives the limiting RBE dose. For a bilateral irradiation, the insula cortices receive the limiting dose. Among the nonencephalic structures, the parotid glands receive RBE doses that were within 15% of the limiting dose.

Absorbed fractions in a voxel-based phantom calculated with the MCNP-4B code.

A new approach for calculating internal dose estimates was developed through the use of a more realistic computational model of the human body. The present technique shows the capability to build a patient-specific phantom with tomography data (a voxel-based phantom) for the simulation of radiation transport and energy deposition using Monte Carlo methods such as in the MCNP-4B code. MCNP-4B absorbed fractions for photons in the mathematical phantom of Snyder et al. agreed well with reference values. Results obtained through radiation transport simulation in the voxel-based phantom, in general, agreed well with reference values. Considerable discrepancies, however, were found in some cases due to two major causes: differences in the organ masses between the phantoms and the occurrence of organ overlap in the voxel-based phantom, which is not considered in the mathematical phantom.

Computerized three-dimensional segmented human anatomy.

Manual segmentation of 129 x-ray CT transverse slices of a living male human has been done and a computerized 3-dimensional volume array modeling all major internal structures of the body has been created. Each voxel of the volume contains a index number designating it as belonging to a given organ or internal structure. The original x-ray CT images were reconstructed in a 512 x 512 matrix with a resolution of 1 mm in the x,y plane. The z-axis resolution is 1 cm from neck to midthigh and 0.5 cm from neck to crown of the head. This volume array represents a high resolution model of the human anatomy and can serve as a voxel-based anthropomorphic phantom suitable for many computer-based modeling and simulation calculations.