Search

The purpose of this study was to develop detailed and realistic models of the cortical and trabecular bones in the spine, ribs, and sternum and incorporate them into the library of virtual human phantoms (XCAT). Cortical bone was modeled by 3D morphological erosion of XCAT homogenously defined bones with an average thickness measured from the CT dataset upon which each individual XCAT phantom was based. The trabecular texture was modeled using a power law synthesis algorithm where the parameters were tuned using high-resolution anatomical images of the Human Visible Female. The synthesized bone textures were added into the XCAT phantoms. To qualitatively evaluate the improved realism of the bone modeling, CT simulations of the XCAT phantoms were acquired with and without the textured bone modeling. The 3D power spectrum of the anatomical images exhibited a power law behavior (R 2 = 0.84), as expected in fractal and porous textures. The proposed texture synthesis algorithm was able to synthesize textures emulating real anatomical images, with the simulated CT images with the prototyped bones were more realistic than those simulated with the original XCAT models. Incorporating intra-organ structures, the "textured" phantoms are envisioned to be used to conduct virtual clinical trials in the context of medical imaging in cases where the actual trials are infeasible due to the lack of ground truth, cost, or potential risks to the patients.

The 4D extended cardiac-torso (XCAT) phantom has provided a valuable tool to study the effects of anatomy and motion on medical images, especially cardiac motion. One limitation of the XCAT was that it did not have a physiological basis which to realistically simulate variations in cardiac function. In this work, we incorporate into the XCAT anatomy the four-chamber FE Living Heart Model (LHM) developed by the Living Heart Project (LHP). The LHM represents the state of the art in cardiac FE simulation because of its ability to accurately replicate the biomechanical motion of the entire heart and its variations. We create a new series of 4D phantoms capable of simulating patients with varying body sizes and shapes; cardiac positions, orientations, and dynamics. While extendable to other imaging modalities and technologies, our goal is to use the FE-enhanced XCAT models to investigate the optimal use of computed tomography (CT) for the evaluation of coronary artery disease (CAD). With the ability to simulate realistic, predictive, patient quality 4D imaging data, the enhanced XCAT models will enable optimization studies to identify the most promising systems or system parameters for further clinical validation.

The purpose of this paper was to extend the extended cardiac-torso (XCAT) series of computational phantoms to include a detailed lung architecture including airways and pulmonary vasculature. Eleven XCAT phantoms of varying anatomy were used in this paper. The lung lobes and initial branches of the airways, pulmonary arteries, and veins were previously defined in each XCAT model. These models were extended from the initial branches of the airways and vessels to the level of terminal branches using an anatomically-based volume-filling branching algorithm. This algorithm grew the airway and vasculature branches separately and iteratively without intersecting each other using cylindrical models with diameters estimated by order-based anatomical measurements. Geometrical features of the extended branches were compared with the literature anatomy values to quantitatively evaluate the models. These features include branching angle, length to diameter ratio, daughter to parent diameter ratio, asymmetrical branching pattern, diameter, and length ratios. The XCAT phantoms were then used to simulate CT images to qualitatively compare them with the original phantom images. The proposed growth model produced 46369 ± 12521 airways, 44737 ± 11773 arteries, and 39819 ± 9988 veins to the XCAT phantoms. Furthermore, the growth model was shown to produce asymmetrical airway, artery, and vein networks with geometrical attributes close to morphometry and model based studies. The simulated CT images of the phantoms were judged to be more realistic, including more airways and pulmonary vessels compared with the original phantoms. Future work will seek to add a heterogeneous parenchymal background into the XCAT lungs to make the phantoms even more representative of human anatomy, paving the way towards the use of XCAT models as a tool to virtually evaluate the current and emerging medical imaging technologies.

The four-dimensional (4-D) eXtended CArdiac-Torso (XCAT) series of phantoms was developed to provide accurate computerized models of the human anatomy and physiology. The XCAT series encompasses a vast population of phantoms of varying ages from newborn to adult, each including parameterized models for the cardiac and respiratory motions. With great flexibility in the XCAT's design, any number of body sizes, different anatomies, cardiac or respiratory motions or patterns, patient positions and orientations, and spatial resolutions can be simulated. As such, the XCAT phantoms are gaining a wide use in biomedical imaging research. There they can provide a virtual patient base from which to quantitatively evaluate and improve imaging instrumentation, data acquisition, techniques, and image reconstruction and processing methods which can lead to improved image quality and more accurate clinical diagnoses. The phantoms have also found great use in radiation dosimetry, radiation therapy, medical device design, and even the security and defense industry. This review paper highlights some specific areas in which the XCAT phantoms have found use within biomedical imaging and other fields. From these examples, we illustrate the increasingly important role that computerized phantoms and computer simulation are playing in the research community.

The purpose of this study was to develop a dynamic physical cardiac phantom with a realistic coronary plaque to investigate stenosis measurement accuracy under clinically relevant heart-rates. The coronary plaque model (5 mm diameter, 50% stenosis, and 32 mm long) was designed and 3D-printed with tissue equivalent materials (calcified plaque with iodine-enhanced lumen). Realistic cardiac motion was modeled by converting computational cardiac motion vectors into compression and rotation profiles executed by a commercial base cardiac phantom. The phantom was imaged on a dual-source CT system applying a retrospective gated coronary CT angiography (CCTA) protocol using synthesized motion-synchronized electrocardiogram (ECG) waveforms. Multiplanar reformatted images were reconstructed along vessel centerlines. Enhanced lumens were segmented by five independent operators. On average, stenosis measurement accuracy was 0.9% positively biased for the motion-free condition. Average measurement accuracy monotonically decreased from 0.9% positive bias for the motion-free condition to 18.5% negative bias at 90 beats per minute. Contrast-to-noise ratio, lumen circularity, and segmentation conformity also decreased monotonically with increasing heart-rate. These results demonstrate successful implementation of a base cardiac phantom with a 3D-printed coronary plaque model, relevant motion profile, and coordinated ECG waveform. They further show the utility of the model to ascertain metrics of CCTA accuracy and image quality under realistic plaque, motion, and acquisition conditions.

Purpose: The limited number of 3D patient-based breast phantoms available could be augmented by synthetic breast phantoms in order to facilitate virtual clinical trials (VCTs) using model observers for breast imaging optimization and evaluation., Methods: These synthetic breast phantoms were developed using Principal Component Analysis (PCA) to reduce the number of dimensions needed to describe a training set of images. PCA decomposed a training set of M breast CT volumes (with millions of voxels each) into an M-1-dimensional space of eigenvectors, which we call eigenbreasts. Each of the training breast phantoms was compactly represented by the mean image plus a weighted sum of eigenbreasts. The distribution of weights observed from training was then sampled to create new synthesized breast phantoms., Results: The resulting synthesized breast phantoms demonstrated a high degree of realism, as supported by an observer study. Two out of three experienced physicist observers were unable to distinguish between the synthesized breast phantoms and the patient-based phantoms. The fibroglandular density and noise power law exponent of the synthesized breast phantoms agreed well with the training data., Conclusions: Our method extends our series of digital breast phantoms based on breast CT data, providing the capability to generate new, statistically varying ensembles consisting of tens of thousands of virtual subjects. This work represents an important step toward conducting future virtual trials for task-based assessment of breast imaging, where it is vital to have a large ensemble of realistic phantoms for statistical power as well as clinical relevance., (© 2017 American Association of Physicists in Medicine.)

Purpose: To develop a 4-dimensional computerized magnetic resonance imaging phantom with image textures extracted from real patient scans for liver motion studies., Methods: The proposed phantom was developed based on the current version of 4-dimensional extended cardiac-torso computerized phantom and a clinical magnetic resonance scan. Initially, the extended cardiac-torso phantom was voxelized in abdominal-chest region at the end of exhalation phase. Structures/tissues were classified into 4 categories: (1) Seven key textured organs, including liver, gallbladder, spleen, stomach, heart, kidneys, and pancreas, were mapped from a clinical T1-weighted liver magnetic resonance scan using deformable registration. (2) Large textured soft tissue volumes were simulated via an iterative pattern generation method using the same magnetic resonance scan. (3) Lung and intestine structures were generated by assigning uniform intensity with proper noise modeling. (4) Bony structures were generated by assigning the magnetic resonance values. A spherical hypointensity tumor was inserted into the liver. Other respiratory phases of the 4-dimensional phantom were generated using the backward deformation vector fields exported by the extended cardiac-torso program, except that bony structures were generated separately for each phase. A weighted image filtering process was utilized to improve the overall tissue smoothness at each phase., Results: Three 4-dimensional series with different motion amplitudes were generated. The developed motion phantom produced good illustrations of abdominal-chest region with anatomical structures in key organs and texture patterns in large soft tissue volumes. In a standard series, the tumor volume was measured as 13.90 ± 0.11 cm 3 in a respiratory cycle and the tumor's maximum center-of-mass shift was 2.95 cm/1.84 cm on superior-inferior/anterior-posterior directions. The organ motion during the respiratory cycle was well rendered. The developed motion phantom has the flexibility of motion pattern variation, organ geometry variation, and tumor modeling variation., Conclusions: A 4-D computerized phantom was developed and could be used to produce image series with synthetic magnetic resonance textures for magnetic resonance imaging research of liver motion.

This study aimed to estimate the organ dose reduction potential for organ-dose-based tube current modulated (ODM) thoracic computed tomography (CT) with a wide dose reduction arc. Twenty-one computational anthropomorphic phantoms (XCAT) were used to create a virtual patient population with clinical anatomic variations. The phantoms were created based on patient images with normal anatomy (age range: 27 to 66 years, weight range: 52.0 to 105.8 kg). For each phantom, two breast tissue compositions were simulated: [Formula: see text] and [Formula: see text] (glandular-to-adipose ratio). A validated Monte Carlo program (PENELOPE, Universitat de Barcelona, Spain) was used to estimate the organ dose for standard tube current modulation (TCM) (SmartmA, GE Healthcare) and ODM (GE Healthcare) for a commercial CT scanner (Revolution, GE Healthcare) using a typical clinical thoracic CT protocol. Both organ dose and [Formula: see text]-to-organ dose conversion coefficients ([Formula: see text] factors) were compared between TCM and ODM. ODM significantly reduced all radiosensitive organ doses ([Formula: see text]). The breast dose was reduced by [Formula: see text]. For [Formula: see text] factors, organs in the anterior region (e.g., thyroid and stomach) exhibited substantial decreases, and the medial, distributed, and posterior region saw either an increase of less than 5% or no significant change. ODM significantly reduced organ doses especially for radiosensitive superficial anterior organs such as the breasts.

The purpose of this study was to investigate relationships between patient attributes and organ dose for a population of computational phantoms for 20 tomosynthesis and radiography protocols. Organ dose was estimated from 54 adult computational phantoms (age: 18 to 78 years, weight 52 to 117 kg) using a validated Monte-Carlo simulation (PENELOPE) of a system capable of performing tomosynthesis and radiography. The geometry and field of view for each exam were modeled to match clinical protocols. For each protocol, the energy deposited in each organ was estimated by the simulations, converted to dose units, and then normalized by exposure in air. Dose to radiosensitive organs was studied as a function of average patient thickness in the region of interest and as a function of body mass index. For tomosynthesis, organ doses were also studied as a function of x-ray tube position. This work developed comprehensive information for organ dose dependencies across a range of tomosynthesis and radiography protocols. The results showed a protocol-dependent exponential decrease with an increasing patient size. There was a variability in organ dose across the patient population, which should be incorporated in the metrology of organ dose. The results can be used to prospectively and retrospectively estimate organ dose for tomosynthesis and radiography.

ECG gated cardiac PET imaging measures functional parameters such as left ventricle (LV) ejection fraction (EF), providing diagnostic and prognostic information for management of patients with coronary artery disease (CAD). Respiratory motion degrades spatial resolution and affects the accuracy in measuring the LV volumes for EF calculation. The goal of this study is to systematically investigate the effect of respiratory motion correction on the estimation of end-diastolic volume (EDV), end-systolic volume (ESV), and EF, especially on the separation of normal and abnormal EFs. We developed a respiratory motion incorporated 4D PET image reconstruction technique which uses all gated-frame data to acquire a motion-suppressed image. Using the standard XCAT phantom and two individual-specific volunteer XCAT phantoms, we simulated dual-gated myocardial perfusion imaging data for normally and abnormally beating hearts. With and without respiratory motion correction, we measured the EDV, ESV, and EF from the cardiac-gated reconstructed images. For all the phantoms, the estimated volumes increased and the biases significantly reduced with motion correction compared with those without. Furthermore, the improvement of ESV measurement in the abnormally beating heart led to better separation of normal and abnormal EFs. The simulation study demonstrated the significant effect of respiratory motion correction on cardiac imaging data with motion amplitude as small as 0.7 cm. The larger the motion amplitude the more improvement respiratory motion correction brought about on the EF measurement. Using data-driven respiratory gating, we also demonstrated the effect of respiratory motion correction on estimating the above functional parameters from list mode patient data. Respiratory motion correction has been shown to improve the accuracy of EF measurement in clinical cardiac PET imaging.

Purpose To develop a method to incorporate the propagation of contrast material into computational anthropomorphic phantoms for estimation of organ dose at computed tomography (CT). Materials and Methods A patient-specific physiologically based pharmacokinetic (PBPK) model of the human cardiovascular system was incorporated into 58 extended cardiac-torso (XCAT) patient phantoms. The PBPK model comprised compartmental models of vessels and organs unique to each XCAT model. For typical injection protocols, the dynamics of the contrast material in the body were described according to a series of patient-specific iodine mass-balance differential equations, the solutions to which provided the contrast material concentration time curves for each compartment. Each organ was assigned to a corresponding time-varying iodinated contrast agent to create the contrast material-enhanced five-dimensional XCAT models, in which the fifth dimension represents the dynamics of contrast material. To validate the accuracy of the models, simulated aortic and hepatic contrast-enhancement results throughout the models were compared with previously published clinical data by using the percentage of discrepancy in the mean, time to 90% peak, peak value, and slope of enhancement in a paired t test at the 95% significance level. Results The PBPK model allowed effective prediction of the time-varying concentration curves of various contrast material administrations in each organ for different patient models. The contrast-enhancement results were in agreement with results of previously published clinical data, with mean percentage, time to 90% peak, peak value, and slope of less than 10% (P > .74), 4%, 7%, and 14% for uniphasic and 12% (P > .56), 4%, 12%, and 14% for biphasic injection protocols, respectively. The exception was hepatic enhancement results calculated for a uniphasic injection protocol for which the discrepancy was less than 25%. Conclusion A technique to model the propagation of contrast material in XCAT human models was developed. The models with added contrast material propagation can be applied to simulate contrast-enhanced CT examinations. © RSNA, 2017 Online supplemental material is available for this article.

Purpose To estimate the radiation dose as a result of contrast medium administration in a typical abdominal computed tomographic (CT) examination across a library of contrast material-enhanced computational patient models. Materials and Methods In part II of this study, first, the technique described in part I of this study was applied to enhance the extended cardiac-torso models with patient-specific iodine-time profiles reflecting the administration of contrast material. Second, the patient models were deployed to assess the patient-specific organ dose as a function of time in a typical abdominal CT examination using Monte Carlo simulation. In this hypothesis-generating study, organ dose refers to the total energy deposited in the unit mass of the tissue inclusive of iodine. Third, a study was performed as a strategy to anticipate the biologically relevant dose (absorbed dose to tissue) in highly perfused organs such as the liver and kidney. The time-varying organ-dose increment values relative to those for unenhanced CT examinations were reported. Results The results from the patient models subjected to the injection protocol indicated up to a total 53%, 30%, 35%, 54%, 27%, 18%, 17%, and 24% increase in radiation dose delivered to the heart, spleen, liver, kidneys, stomach, colon, small intestine, and pancreas, respectively. The biologically relevant dose increase with respect to the dose at an unenhanced CT examination was in the range of 0%-18% increase for the liver and 27% for the kidney across 58 patient models. Conclusion The administration of contrast medium increases the total radiation dose. However, radiation dose, while relevant to be included in estimating the risk associated with contrast-enhanced CT, may still not fully characterize the total biologic effects. Therefore, given the fact that many CT diagnostic decisions would be impossible without the use of iodine, this study suggests the need to consider the effect of iodinated contrast material on the organ doses to patients undergoing CT studies when designing CT protocols. © RSNA, 2017 Online supplemental material is available for this article.

Background: The estimation of organ doses and effective doses for children receiving CT examinations is of high interest. Newer, more realistic anthropomorphic body models can provide information on individual organ doses and improved estimates of effective dose., Materials and Methods: Previously developed body models representing 50th-percentile individuals at reference ages (newborn, 1, 5, 10 and 15 years) were modified to represent 10th, 25th, 75th and 90th height percentiles for both genders and an expanded range of ages (3, 8 and 13 years). We calculated doses for 80 pediatric reference phantoms from simulated chest-abdomen-pelvis exams on a model of a Philips Brilliance 64 CT scanner. Individual organ and effective doses were normalized to dose-length product (DLP) and fit as a function of body diameter., Results: We calculated organ and effective doses for 80 reference phantoms and plotted them against body diameter. The data were well fit with an exponential function. We found DLP-normalized organ dose to correlate strongly with body diameter (R 2 >0.95 for most organs). Similarly, we found a very strong correlation with body diameter for DLP-normalized effective dose (R 2 >0.99). Our results were compared to other studies and we found average agreement of approximately 10%., Conclusion: We provide organ and effective doses for a total of 80 reference phantoms representing normal-stature children ranging in age and body size. This information will be valuable in replacing the types of vendor-reported doses available. These data will also permit the recording and tracking of individual patient doses. Moreover, this comprehensive dose database will facilitate patient matching and the ability to predict patient-individualized dose prior to examination.

In Monte Carlo simulation of CT dose, many input parameters are required (e.g. bowtie filter properties and scan start/end location). Our goal was to examine the uncertainties in patient dose when input parameters were inaccurate. Using a validated Monte Carlo program, organ dose from a chest CT scan was simulated for an average-size female phantom using a reference set of input parameter values (treated as the truth). Additional simulations were performed in which errors were purposely introduced into the input parameter values. The effects on four dose quantities were analyzed: organ dose (mGy/mAs), effective dose (mSv/mAs), CTDI vol -normalized organ dose (unitless), and DLP-normalized effective dose (mSv/mGy · cm). At 120 kVp, when spectral half value layer deviated from its true value by  ±1.0 mm Al, the four dose quantities had errors of 18%, 7%, 14% and 2%, respectively. None of the dose quantities were affected significantly by errors in photon path length through the graphite section of the bowtie filter; path length error as large as 5 mm produced dose errors of  ⩽2%. In contrast, error of this magnitude in the aluminum section produced dose errors of  ⩽14%. At a total collimation of 38.4 mm, when radiation beam width deviated from its true value by  ±  3 mm, dose errors were  ⩽7%. Errors in tube starting angle had little impact on effective dose (errors  ⩽  1%); however, they produced organ dose errors as high as 66%. When the assumed scan length was longer by 4 cm than the truth, organ dose errors were up to 137%. The corresponding error was 24% for effective dose, but only 3% for DLP-normalized effective dose. Lastly, when the scan isocenter deviated from the patient's anatomical center by 5 cm, organ and effective dose errors were up 18% and 8%, respectively.

Purpose: Amplitude-based respiratory gating is known to capture the extent of respiratory motion (RM) accurately but results in residual motion in the presence of respiratory hysteresis. In our previous study, we proposed and developed a novel approach to account for respiratory hysteresis by applying the Bouc-Wen (BW) model of hysteresis to external surrogate signals of anterior/posterior motion of the abdomen and chest with respiration. In this work, using simulated and clinical SPECT myocardial perfusion imaging (MPI) studies, we investigate the effects of respiratory hysteresis and evaluate the benefit of correcting it using the proposed BW model in comparison with the abdomen signal typically employed clinically., Methods: The MRI navigator data acquired in free-breathing human volunteers were used in the specially modified 4D NCAT phantoms to allow simulating three types of respiratory patterns: monotonic, mild hysteresis, and strong hysteresis with normal myocardial uptake, and perfusion defects in the anterior, lateral, inferior, and septal locations of the mid-ventricular wall. Clinical scans were performed using a Tc-99m sestamibi MPI protocol while recording respiratory signals from thoracic and abdomen regions using a visual tracking system (VTS). The performance of the correction using the respiratory signals was assessed through polar map analysis in phantom and 10 clinical studies selected on the basis of having substantial RM., Results: In phantom studies, simulations illustrating normal myocardial uptake showed significant differences (P < 0.001) in the uniformity of the polar maps between the RM uncorrected and corrected. No significant differences were seen in the polar map uniformity across the RM corrections. Studies simulating perfusion defects showed significantly decreased errors (P < 0.001) in defect severity and extent for the RM corrected compared to the uncorrected. Only for the strong hysteretic pattern, there was a significant difference (P < 0.001) among the RM corrections. The errors in defect severity and extent for the RM correction using abdomen signal were significantly higher compared to that of the BW (severity = -4.0%, P < 0.001; extent = -65.4%, P < 0.01) and chest (severity = -4.1%, P < 0.001; extent = -52.5%, P < 0.01) signals. In clinical studies, the quantitative analysis of the polar maps demonstrated qualitative and quantitative but not statistically significant differences (P = 0.73) between the correction methods that used the BW signal and the abdominal signal., Conclusions: This study shows that hysteresis in respiration affects the extent of residual motion left in the RM-binned data, which can impact wall uniformity and the visualization of defects. Thus, there appears to be the potential for improved accuracy in reconstruction in the presence of hysteretic RM with the BW model method providing a possible step in the direction of improvement., (© 2016 American Association of Physicists in Medicine.)

Purpose: This study aimed to investigate the breast dose reduction potential of a breast-positioning (BP) technique for thoracic CT examinations with organ-based tube current modulation (OTCM)., Methods: This study included 13 female anthropomorphic computational phantoms (XCAT, age range: 27-65 y.o., weight range: 52-105.8 kg). Each phantom was modified to simulate three breast sizes in standard supine geometry. The modeled breasts were then morphed to emulate BP that constrained the majority of the breast tissue inside the 120° anterior tube current (mA) reduction zone. The OTCM mA value was modeled using a ray-tracing program, which reduced the mA to 20% in the anterior region with a corresponding increase to the posterior region. The organ doses were estimated by a validated Monte Carlo program for a typical clinical CT system (SOMATOM Definition Flash, Siemens Healthcare). The simulated organ doses and organ doses normalized by CTDI vol were used to compare three CT protocols: attenuation-based tube current modulation (ATCM), OTCM, and OTCM with BP (OTCM BP )., Results: On average, compared to ATCM, OTCM reduced breast dose by 19.3 ± 4.5%, whereas OTCM BP reduced breast dose by 38.6 ± 8.1% (an additional 23.8 ± 9.4%). The dose saving of OTCM BP was more significant for larger breasts (on average 33, 38, and 44% reduction for 0.5, 1, and 2 kg breasts, respectively). Compared to ATCM, OTCM BP also reduced thymus and heart dose by 15.1 ± 7.4% and 15.9 ± 6.2% respectively., Conclusions: In thoracic CT examinations, OTCM with a breast-positioning technique can markedly reduce unnecessary exposure to radiosensitive organs in anterior chest wall, specifically breast tissue. The breast dose reduction is more notable for women with larger breasts., (© 2016 American Association of Physicists in Medicine.)

Currently, computed tomography (CT) dosimetry relies on surrogates for dose, such as CT dose index and size-specific dose estimates, rather than dose per se . Organ dose is considered as the gold standard for radiation dosimetry. However, organ dose estimation requires precise knowledge of organ locations. Regional imparted energy and dose can also be used to quantify radiation burden and are beneficial because they do not require knowledge of organ size or location. This work investigated an automated technique to retrospectively estimate the imparted energy from tube current-modulated (TCM) CT exams across 13 protocols. Monte Carlo simulations of various head and body TCM CT examinations across various tube potentials and TCM strengths were performed on 58 adult computational extended cardiac-torso phantoms to develop relationships between scanned mass and imparted energy normalized by dose length product. Results from the Monte Carlo simulations indicate that normalized imparted energy increases with increasing both scanned mass and tube potential, but it is relatively unaffected by the strength of the TCM. The automated algorithm was tested on 40 clinical datasets with a 98% success rate.

This study aims to characterize the effect of background tissue density and heterogeneity on the detection of irregular masses in breast tomosynthesis, while demonstrating the capability of the sophisticated tools that can be used in the design, implementation, and performance analysis of virtual clinical trials (VCTs). Twenty breast phantoms from the extended cardiac-torso (XCAT) family, generated based on dedicated breast computed tomography of human subjects, were used to extract a total of 2173 volumes of interest (VOIs) from simulated tomosynthesis images. Five different lesions, modeled after human subject tomosynthesis images, were embedded in the breasts and combined with the lesion absent condition yielded a total of [Formula: see text] VOIs. Effects of background tissue density and heterogeneity on the detection of the lesions were studied by implementing a composite hypothesis signal detection paradigm with location known exactly, lesion known exactly or statistically, and background known statistically. Using the area under the receiver operating characteristic curve, detection performance deteriorated as density was increased, yielding findings consistent with clinical studies. A human observer study was performed on a subset of the simulated tomosynthesis images, confirming the detection performance trends with respect to density and serving as a validation of the implemented detector. Performance of the implemented detector varied substantially across the 20 breasts. Furthermore, background tissue density and heterogeneity affected the log-likelihood ratio test statistic differently under lesion absent and lesion present conditions. Therefore, considering background tissue variability in tissue models can change the outcomes of a VCT and is hence of crucial importance. The XCAT breast phantoms have the potential to address this concern by offering realistic modeling of background tissue variability based on a wide range of human subjects, comprising various breast shapes, sizes, and densities.

Estimating organ dose for clinical patients requires accurate modeling of the patient anatomy and the dose field of the CT exam. The modeling of patient anatomy can be achieved using a library of representative computational phantoms (Samei et al 2014 Pediatr. Radiol. 44 460-7). The modeling of the dose field can be challenging for CT exams performed with a tube current modulation (TCM) technique. The purpose of this work was to effectively model the dose field for TCM exams using a convolution-based method. A framework was further proposed for prospective and retrospective organ dose estimation in clinical practice. The study included 60 adult patients (age range: 18-70 years, weight range: 60-180 kg). Patient-specific computational phantoms were generated based on patient CT image datasets. A previously validated Monte Carlo simulation program was used to model a clinical CT scanner (SOMATOM Definition Flash, Siemens Healthcare, Forchheim, Germany). A practical strategy was developed to achieve real-time organ dose estimation for a given clinical patient. CTDIvol-normalized organ dose coefficients ([Formula: see text]) under constant tube current were estimated and modeled as a function of patient size. Each clinical patient in the library was optimally matched to another computational phantom to obtain a representation of organ location/distribution. The patient organ distribution was convolved with a dose distribution profile to generate [Formula: see text] values that quantified the regional dose field for each organ. The organ dose was estimated by multiplying [Formula: see text] with the organ dose coefficients ([Formula: see text]). To validate the accuracy of this dose estimation technique, the organ dose of the original clinical patient was estimated using Monte Carlo program with TCM profiles explicitly modeled. The discrepancy between the estimated organ dose and dose simulated using TCM Monte Carlo program was quantified. We further compared the convolution-based organ dose estimation method with two other strategies with different approaches of quantifying the irradiation field. The proposed convolution-based estimation method showed good accuracy with the organ dose simulated using the TCM Monte Carlo simulation. The average percentage error (normalized by CTDIvol) was generally within 10% across all organs and modulation profiles, except for organs located in the pelvic and shoulder regions. This study developed an improved method that accurately quantifies the irradiation field under TCM scans. The results suggested that organ dose could be estimated in real-time both prospectively (with the localizer information only) and retrospectively (with acquired CT data).

Purpose: The authors are developing a series of computational breast phantoms based on breast CT data for imaging research. In this work, the authors develop a program that will allow a user to alter the phantoms to simulate the effect of gravity and compression of the breast (craniocaudal or mediolateral oblique) making the phantoms applicable to multimodality imaging., Methods: This application utilizes a template finite-element (FE) breast model that can be applied to their presegmented voxelized breast phantoms. The FE model is automatically fit to the geometry of a given breast phantom, and the material properties of each element are set based on the segmented voxels contained within the element. The loading and boundary conditions, which include gravity, are then assigned based on a user-defined position and compression. The effect of applying these loads to the breast is computed using a multistage contact analysis in FEBio, a freely available and well-validated FE software package specifically designed for biomedical applications. The resulting deformation of the breast is then applied to a boundary mesh representation of the phantom that can be used for simulating medical images. An efficient script performs the above actions seamlessly. The user only needs to specify which voxelized breast phantom to use, the compressed thickness, and orientation of the breast., Results: The authors utilized their FE application to simulate compressed states of the breast indicative of mammography and tomosynthesis. Gravity and compression were simulated on example phantoms and used to generate mammograms in the craniocaudal or mediolateral oblique views. The simulated mammograms show a high degree of realism illustrating the utility of the FE method in simulating imaging data of repositioned and compressed breasts., Conclusions: The breast phantoms and the compression software can become a useful resource to the breast imaging research community. These phantoms can then be used to evaluate and compare imaging modalities that involve different positioning and compression of the breast.

Cubic Hermite hexahedral finite element meshes have some well-known advantages over linear tetrahedral finite element meshes in biomechanical and anatomic modeling using isogeometric analysis. These include faster convergence rates as well as the ability to easily model rule-based anatomic features such as cardiac fiber directions. However, it is not possible to create closed complex objects with only regular nodes; these objects require the presence of extraordinary nodes (nodes with 3 or >= 5 adjacent elements in 2D) in the mesh. The presence of extraordinary nodes requires new constraints on the derivatives of adjacent elements to maintain continuity. We have developed a new method that uses an ensemble coordinate frame at the nodes and a local-to-global mapping to maintain continuity. In this paper, we make use of this mapping to create cubic Hermite models of the human ventricles and a four-chamber heart. We also extend the methods to the finite element equations to perform biomechanics simulations using these meshes. The new methods are validated using simple test models and applied to anatomically accurate ventricular meshes with valve annuli to simulate complete cardiac cycle simulations.

The objective of this investigation was to determine the effectiveness of three motion reducing strategies in diminishing the degrading impact of respiratory motion on the detection of small solitary pulmonary nodules (SPN) in single photon emission computed tomographic (SPECT) imaging in comparison to a standard clinical acquisition and the ideal case of imaging in the absence of respiratory motion. To do this non-uniform rational B-spline cardiac-torso (NCAT) phantoms based on human-volunteer CT studies were generated spanning the respiratory cycle for a normal background distribution of Tc-99m NeoTect. Similarly, spherical phantoms of 1.0 cm diameter were generated to model small SPN for each of 150 uniquely located sites within the lungs whose respiratory motion was based on the motion of normal structures in the volunteer CT studies. The SIMIND Monte Carlo program was used to produce SPECT projection data from these. Normal and single-lesion containing SPECT projection sets with a clinically realistic Poisson noise level were created for the cases of: 1) the end-expiration (EE) frame with all counts, 2) respiration-averaged motion with all counts, 3) one-fourth of the 32 frames centered around EE (Quarter-Binning), 4) one-half of the 32 frames centered around EE (Half-Binning), and 5) eight temporally binned frames spanning the respiratory cycle. Each of the sets of combined projection data were reconstructed with RBI-EM with system spatial-resolution compensation (RC). Based on the known motion for each of the 150 different lesions, the reconstructed volumes of respiratory bins were shifted so as to superimpose the locations of the SPN onto that in the first bin (Reconstruct and Shift). Five human-observers performed localization receiver operating characteristics (LROC) studies of SPN detection. The observer results were analyzed for statistical significance differences in SPN detection accuracy among the three correction strategies, the standard acquisition, and the ideal case of the absence of respiratory motion. Our human-observer LROC determined that Quarter-Binning and Half-Binning strategies resulted in SPN detection accuracy statistically significantly below (P < 0.05) that of standard clinical acquisition, whereas the Reconstruct and Shift strategy resulted in a detection accuracy not statistically significantly different from that of the ideal case. This investigation demonstrates that tumor detection based on acquisitions associated with less than all the counts which could potentially be employed may result in poorer detection despite limiting the motion of the lesion. The Reconstruct and Shift method results in tumor detection that is equivalent to ideal motion correction.

The objectives of this investigation were to model the respiratory motion of solitary pulmonary nodules (SPN) and then use this model to determine the impact of respiratory motion on the localization and detection of small SPN in SPECT imaging for four reconstruction strategies. The respiratory motion of SPN was based on that of normal anatomic structures in the lungs determined from breath-held CT images of a volunteer acquired at two different stages of respiration. End-expiration (EE) and time-averaged (Frame Av) non-uniform-B-spline cardiac torso (NCAT) digital-anthropomorphic phantoms were created using this information for respiratory motion within the lungs. SPN were represented as 1 cm diameter spheres which underwent linear motion during respiration between the EE and end-inspiration (EI) time points. The SIMIND Monte Carlo program was used to produce SPECT projection data simulating Tc-99m depreotide (NeoTect) imaging. The projections were reconstructed using 1) no correction (NC), 2) attenuation correction (AC), 3) resolution compensation (RC), and 4) attenuation correction, scatter correction, and resolution compensation (AC&#95;SC&#95;RC). A human-observer localization receiver operating characteristics (LROC) study was then performed to determine the difference in localization and detection accuracy with and without the presence of respiratory motion. The LROC comparison determined that respiratory motion degrades tumor detection for all four reconstruction strategies, thus correction for SPN motion would be expected to improve detection accuracy. The inclusion of RC in reconstruction improved detection accuracy for both EE and Frame Av over NC and AC. Also the magnitude of the impact of motion was least for AC&#95;SC&#95;RC.

Purpose: To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data., Methods: A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms., Results: After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms., Conclusions: This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

Purpose: In CT imaging, a desirable quality assurance dose quantity should account for the dose variability across scan parameters and scanner models. Recently, AAPM Task Group 111 proposed to use equilibrium dose-pitch product (Dˆeq), in place of CT dose index (CTDI100), for scan modes involving table translation. The purpose of this work is to investigate whether this new concept better accounts for the tube voltage (kVp) dependence of organ dose than the conventional CTDI100., Methods: Three extended cardiac-torso patient models were included in this study. They represented normal-weight, overweight, and obese patients with abdomen-pelvis diameters ranging between 23 and 36 cm and body mass indices ranging between 20 and 42. A Monte Carlo program developed and validated for a 128-slice CT system was used to simulate organ dose for abdomen-pelvis scans at five tube voltages (70, 80, 100, 120, 140 kVp) with a pitch of 0.8 and a collimation of 38.4 mm. The same Monte Carlo program was also used to obtain CTDI100 and Dˆeq as well as their volume-averaged values, CTDIvol and Dˆeq,vol., Results: With other scan parameters kept constant, organ dose itself depended strongly on kVp. For the normal-weight patient model, the coefficient of variation (COV) across the five kVp values ranged between 72% and 75% for nine organs (liver, spleen, stomach, pancreas, kidneys, colon, small intestine, bladder, and ovaries) fully encompassed by the primary radiation beam. The COV generally increased with increasing patient size, ranging between 72%-77% and 76%-81% for the overweight and obese patient models, respectively. One-way analysis of variance for the effect of kVp was highly significant for all patient models (p<2×10(-26)). When organ dose was normalized by CTDIvol, the COV across kVps reduced to 5%-11%, 6%-15%, 12%-22% for the normal-weight, overweight, and obese patients, respectively. The effect of kVp was still highly significant (p=0.0001-0.004). When organ dose was normalized by Dˆeq,vol, the COV further reduced to 1%-8%, 3%-11%, 9%-19%, respectively. The effect of kVp was still significant for the obese patient model (p=0.004), but no longer significant for the normal-weight and overweight patient models (p=0.4 and 0.09, respectively). Finally, if organ dose conversion factors obtained at 120 kVp were used to approximate the values at 70 and 80 kVp, the resulting errors in the estimated organ dose were significantly reduced when the conversion factors were based on Dˆeq,vol instead of CTDIvol., Conclusions: In adult abdomen-pelvis CT, equilibrium dose-pitch product better accounts for the kVp dependence of organ dose than CTDI100.

Previously, the authors developed a series of eight realistic digital mouse and rat whole body phantoms based on NURBS technology to facilitate internal and external dose calculations in various species of rodents. In this paper, two body phantoms of adult beagles are described based on voxel images converted to NURBS models. Specific absorbed fractions for activity in 24 organs are presented in these models. CT images were acquired of an adult male and female beagle. The images were segmented, and the organs and structures were modeled using NURBS surfaces and polygon meshes. Each model was voxelized at a resolution of 0.75 × 0.75 × 2 mm. The voxel versions were implemented in GEANT4 radiation transport codes to calculate specific absorbed fractions (SAFs) using internal photon and electron sources. Photon and electron SAFs were then calculated for relevant organs in both models. The SAFs for photons and electrons were compatible with results observed by others. Absorbed fractions for electrons for organ self-irradiation were significantly less than 1.0 at energies above 0.5 MeV, as expected for many of these small-sized organs, and measurable cross irradiation was observed for many organ pairs for high-energy electrons (as would be emitted by nuclides like 32P, 90Y, or 188Re). The SAFs were used with standardized decay data to develop dose factors (DFs) for radiation dose calculations using the RADAR Method. These two new realistic models of male and female beagle dogs will be useful in radiation dosimetry calculations for external or internal simulated sources.

Purpose: We previously developed a set of highly detailed 4D reference pediatric extended cardiac-torso (XCAT) phantoms at ages of newborn, 1, 5, 10, and 15 yr with organ and tissue masses matched to ICRP Publication 89 values. In this work, we extended this reference set to a series of 64 pediatric phantoms of varying age and height and body mass percentiles representative of the public at large. The models will provide a library of pediatric phantoms for optimizing pediatric imaging protocols., Methods: High resolution positron emission tomography-computed tomography data obtained from the Duke University database were reviewed by a practicing experienced radiologist for anatomic regularity. The CT portion of the data was then segmented with manual and semiautomatic methods to form a target model defined using nonuniform rational B-spline surfaces. A multichannel large deformation diffeomorphic metric mapping algorithm was used to calculate the transform from the best age matching pediatric XCAT reference phantom to the patient target. The transform was used to complete the target, filling in the nonsegmented structures and defining models for the cardiac and respiratory motions. The complete phantoms, consisting of thousands of structures, were then manually inspected for anatomical accuracy. The mass for each major tissue was calculated and compared to linearly interpolated ICRP values for different ages., Results: Sixty four new pediatric phantoms were created in this manner. Each model contains the same level of detail as the original XCAT reference phantoms and also includes parameterized models for the cardiac and respiratory motions. For the phantoms that were 10 yr old and younger, we included both sets of reproductive organs. This gave them the capability to simulate both male and female anatomy. With this, the population can be expanded to 92. Wide anatomical variation was clearly seen amongst the phantom models, both in organ shape and size, even for models of the same age and sex. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities., Conclusions: This work provides a large cohort of highly detailed pediatric phantoms with 4D capabilities of varying age, height, and body mass. The population of phantoms will provide a vital tool with which to optimize 3D and 4D pediatric imaging devices and techniques in terms of image quality and radiation-absorbed dose.

Purpose: We investigate the differences without/with respiratory motion correction in apparent imaging agent localization induced in reconstructed emission images when the attenuation maps used for attenuation correction (from CT) are misaligned with the patient anatomy during emission imaging due to differences in respiratory state., Methods: We investigated use of attenuation maps acquired at different states of a 2 cm amplitude respiratory cycle (at end-expiration, at end-inspiration, the center map, the average transmission map, and a large breath-hold beyond range of respiration during emission imaging) to correct for attenuation in MLEM reconstruction for several anatomical variants of the NCAT phantom which included both with and without non-rigid motion between heart and sub-diaphragmatic regions (such as liver, kidneys etc). We tested these cases with and without emission motion correction and attenuation map alignment/non-alignment., Results: For the NCAT default male anatomy the false count-reduction due to breathing was largely removed upon emission motion correction for the large majority of the cases. Exceptions (for the default male) were for the cases when using the large-breathhold end-inspiration map (TI&#95;EXT), when we used the end-expiration (TE) map, and to a smaller extent, the end-inspiration map (TI). However moving the attenuation maps rigidly to align the heart region, reduced the remaining count-reduction artifacts. For the female patient count-reduction remained post motion correction using rigid map-alignment due to the breast soft-tissue misalignment. Quantitatively, after the transmission (rigid) alignment correction, the polar-map 17-segment RMS error with respect to the reference (motion-less case) reduced by 46.5% on average for the extreme breathhold case. The reductions were 40.8% for end-expiration map and 31.9% for end-inspiration cases on the average, comparable to the semi-ideal case where each state uses its own attenuation map for correction., Conclusions: Two main conclusions are that even rigid emission motion correction to rigidly align the heart region to the attenuation map helps in average cases to reduce the count-reduction artifacts and secondly, within the limits of the study (ex. rigid correction) when there is lung tissue inferior to the heart as with the NCAT phantom employed in this study endexpiration maps (TE) might best be avoided as they may create more artifacts than the end-inspiration (TI) maps.

X-ray phase-contrast imaging (XPCI) can dramatically improve soft tissue contrast in X-ray medical imaging. Despite worldwide efforts to develop novel XPCI systems, a numerical framework to rigorously predict the performance of a clinical XPCI system at a human scale is not yet available. We have developed such a tool by combining a numerical anthropomorphic phantom defined with non-uniform rational B-splines (NURBS) and a wave optics-based simulator that can accurately capture the phase-contrast signal from a human-scaled numerical phantom. Using a synchrotron-based, high-performance XPCI system, we provide qualitative comparison between simulated and experimental images. Our tool can be used to simulate the performance of XPCI on various disease entities and compare proposed XPCI systems in an unbiased manner.

Purpose: Physical phantoms are essential for the development, optimization, and evaluation of x-ray breast imaging systems. Recognizing the major effect of anatomy on image quality and clinical performance, such phantoms should ideally reflect the three-dimensional structure of the human breast. Currently, there is no commercially available three-dimensional physical breast phantom that is anthropomorphic. The authors present the development of a new suite of physical breast phantoms based on human data., Methods: The phantoms were designed to match the extended cardiac-torso virtual breast phantoms that were based on dedicated breast computed tomography images of human subjects. The phantoms were fabricated by high-resolution multimaterial additive manufacturing (3D printing) technology. The glandular equivalency of the photopolymer materials was measured relative to breast tissue-equivalent plastic materials. Based on the current state-of-the-art in the technology and available materials, two variations were fabricated. The first was a dual-material phantom, the Doublet. Fibroglandular tissue and skin were represented by the most radiographically dense material available; adipose tissue was represented by the least radiographically dense material. The second variation, the Singlet, was fabricated with a single material to represent fibroglandular tissue and skin. It was subsequently filled with adipose-equivalent materials including oil, beeswax, and permanent urethane-based polymer. Simulated microcalcification clusters were further included in the phantoms via crushed eggshells. The phantoms were imaged and characterized visually and quantitatively., Results: The mammographic projections and tomosynthesis reconstructed images of the fabricated phantoms yielded realistic breast background. The mammograms of the phantoms demonstrated close correlation with simulated mammographic projection images of the corresponding virtual phantoms. Furthermore, power-law descriptions of the phantom images were in general agreement with real human images. The Singlet approach offered more realistic contrast as compared to the Doublet approach, but at the expense of air bubbles and air pockets that formed during the filling process., Conclusions: The presented physical breast phantoms and their matching virtual breast phantoms offer realistic breast anatomy, patient variability, and ease of use, making them a potential candidate for performing both system quality control testing and virtual clinical trials.

Purpose: Computed tomography (CT) has been widely used worldwide as a tool for medical diagnosis and imaging. However, despite its significant clinical benefits, CT radiation dose at the population level has become a subject of public attention and concern. In this light, optimizing radiation dose has become a core responsibility for the CT community. As a fundamental step to manage and optimize dose, it may be beneficial to have accurate and prospective knowledge about the radiation dose for an individual patient. In this study, the authors developed a framework to prospectively estimate organ dose for chest and abdominopelvic CT exams under tube current modulation (TCM)., Methods: The organ dose is mainly dependent on two key factors: patient anatomy and irradiation field. A prediction process was developed to accurately model both factors. To model the anatomical diversity and complexity in the patient population, the authors used a previously developed library of computational phantoms with broad distributions of sizes, ages, and genders. A selected clinical patient, represented by a computational phantom in the study, was optimally matched with another computational phantom in the library to obtain a representation of the patient's anatomy. To model the irradiation field, a previously validated Monte Carlo program was used to model CT scanner systems. The tube current profiles were modeled using a ray-tracing program as previously reported that theoretically emulated the variability of modulation profiles from major CT machine manufacturers Li et al., [Phys. Med. Biol. 59, 4525-4548 (2014)]. The prediction of organ dose was achieved using the following process: (1) CTDIvol-normalized-organ dose coefficients (horgan) for fixed tube current were first estimated as the prediction basis for the computational phantoms; (2) each computation phantom, regarded as a clinical patient, was optimally matched with one computational phantom in the library; (3) to account for the effect of the TCM scheme, a weighted organ-specific CTDIvol [denoted as CTDIvol organ,weighted] was computed for each organ based on the TCM profile and the anatomy of the "matched" phantom; (4) the organ dose was predicted by multiplying the weighted organ-specific CTDIvol with the organ dose coefficients (horgan). To quantify the prediction accuracy, each predicted organ dose was compared with the corresponding organ dose simulated from the Monte Carlo program with the TCM profile explicitly modeled., Results: The predicted organ dose showed good agreements with the simulated organ dose across all organs and modulation profiles. The average percentage error in organ dose estimation was generally within 20% across all organs and modulation profiles, except for organs located in the pelvic and shoulder regions. For an average CTDIvol of a CT exam of 10 mGy, the average error at full modulation strength (α = 1) across all organs was 0.91 mGy for chest exams, and 0.82 mGy for abdominopelvic exams., Conclusions: This study developed a quantitative model to predict organ dose for clinical chest and abdominopelvic scans. Such information may aid in the design of optimized CT protocols in relation to a targeted level of image quality.

Purpose: With increased focus on radiation dose from medical imaging, prospective radiation dose estimates are becoming increasingly desired. Using available populations of adult and pediatric patient phantoms, radiation dose calculations can be catalogued and prospectively applied to individual patients that best match certain anatomical characteristics. In doing so, the knowledge of organ size and location is a required element. Here, the authors develop a predictive model of organ locations and volumes based on an analysis of adult and pediatric computed tomography (CT) data., Methods: Fifty eight adult and 69 pediatric CT datasets were segmented and utilized in the study. The maximum and minimum points of the organs were recorded with respect to the axial distance from the tip of the sacrum. The axial width, midpoint, and volume of each organ were calculated. Linear correlations between these three organ parameters and patient age, BMI, weight, and height were determined., Results: No statistically significant correlations were found in adult patients between the axial width, midpoint, and volume of the organs versus the patient age or BMI. Slight, positive linear trends were found for organ midpoint versus patient weight (max r(2) = 0.382, mean r(2) = 0.236). Similar trends were found for organ midpoint versus height (max r(2) = 0.439, mean r(2) = 0.200) and for organ volume versus height (max r(2) = 0.410, mean r(2) = 0.153). Gaussian fits performed on probability density functions of the adult organs resulted in r(2)-values ranging from 0.96 to 0.996. The pediatric patients showed much stronger correlations overall. Strong correlations were observed between organ axial midpoint versus age, height, and weight (max r(2) = 0.842, mean r(2) = 0.790; max r(2) = 0.949, mean r(2) = 0.894; and max r(2) = 0.870, mean r(2) = 0.847, respectively). Moderate linear correlations were also observed for organ axial width versus height (max r(2) = 0.772, mean r(2) = 0.562) and for organ volume versus height (max r(2) = 0.781, mean r(2) = 0.601)., Conclusions: Adult patients exhibited small variations in organ volume and location with respect to height and weight, but no meaningful correlation existed between these parameters and age or BMI. Once adulthood is reached, organ morphology and positioning seem to remain static. However, clear trends are evident between pediatric organ locations versus age, height, and weight. Such information can be incorporated into a matching methodology that may provide the highest probability of representing the anatomy of a patient undergoing a clinical exam to prospectively estimate the radiation dose.

Among the various metrics to quantify CT radiation dose, organ dose is generally regarded as one of the best to reflect patient radiation burden. Organ dose is dependent on two main factors, namely patient anatomy and irradiation field. An accurate estimation of organ dose requires detailed modeling of both factors. The modeling of patient anatomy needs to reflect the anatomical diversity and complexity across the population so that the attributes of a given clinical patient can be properly accounted for. The modeling of the irradiation field needs to accurately reflect the CT system condition, especially the tube current modulation (TCM) technique. We present an atlas-based method to model patient anatomy via a library of computational phantoms with representative ages, sizes and genders. A clinical patient is matched with a corresponding computational phantom to obtain a representation of patient anatomy. The irradiation field of the CT system is modeled using a validated Monte Carlo simulation program. The tube current modulation profiles are simulated using a manufacturer-generalizable ray-tracing algorithm. Combining the patient model, Monte Carlo results, and TCM profile, organ doses are obtained by multiplying organ dose values from a fixed mA scan (normalized to CTDIvol-normalized, denoted as h organ ) and an adjustment factor that reflects the specific irradiation of each organ. The accuracy of the proposed method was quantified by simulating clinical abdominopelvic examinations of 58 patients. The predicted organ doses showed good agreement with simulated organ dose across all organs and modulation schemes. For an average CTDIvol of a CT exam of 10 mGy, the absolute median error across all organs was 0.64 mGy (-0.21 and 0.97 for 25th and 75th percentiles, respectively). The percentage differences were within 15%. The study demonstrates that it is feasible to estimate organ doses in clinical CT examinations for protocols without and with tube current modulation. The methodology can be used for both prospective and retrospective estimation of organ dose.

The radiation dose involved in any medical imaging modality that uses ionizing radiation needs to be well understood by the medical physics and clinical community. This is especially true of screening modalities. Digital breast tomosynthesis (DBT) has recently been introduced into the clinic and is being used for screening for breast cancer in the general population. Therefore, it is important that the medical physics community have the required information to be able to understand, estimate, and communicate the radiation dose levels involved in breast tomosynthesis imaging. For this purpose, the American Association of Physicists in Medicine Task Group 223 on Dosimetry in Tomosynthesis Imaging has prepared this report that discusses dosimetry in breast imaging in general, and describes a methodology and provides the data necessary to estimate mean breast glandular dose from a tomosynthesis acquisition. In an effort to maximize familiarity with the procedures and data provided in this Report, the methodology to perform the dose estimation in DBT is based as much as possible on that used in mammography dose estimation.

Body CT scans are routinely performed using tube-current-modulation (TCM) technology. There is notable variability across CT manufacturers in terms of how TCM technology is implemented. Some manufacturers aim to provide uniform image noise across body regions and patient sizes, whereas others aim to provide lower noise for smaller patients. The purpose of this study was to conduct a theoretical investigation to understand how manufacturer-dependent TCM scheme affects organ dose, and to develop a generic approach for assessing organ dose across TCM schemes. The adult reference female extended cardiac-torso (XCAT) phantom was used for this study. A ray-tracing method was developed to calculate the attenuation of the phantom for a given projection angle based on phantom anatomy, CT system geometry, x-ray energy spectrum, and bowtie filter filtration. The tube current (mA) for a given projection angle was then calculated as a log-linear function of the attenuation along that projection. The slope of this function, termed modulation control strength, α, was varied from 0 to 1 to emulate the variability in TCM technology. Using a validated Monte Carlo program, organ dose was simulated for five α values (α = 0, 0.25, 0.5, 0.75, and 1) in the absence and presence of a realistic system mA limit. Organ dose was further normalized by volume-weighted CT dose index (CTDIvol) to obtain conversion factors (h factors) that are relatively independent of system specifics and scan parameters. For both chest and abdomen-pelvis scans and for 24 radiosensitive organs, organ dose conversion factors varied with α, following second-order polynomial equations. This result suggested the need for α-specific organ dose conversion factors (i.e., conversion factors specific to the modulation scheme used). On the other hand, across the full range of α values, organ dose in a TCM scan could be derived from the conversion factors established for a fixed-mA scan (hFIXED). This was possible by multiplying hFIXED by a revised definition of CTDIvol that accounts for two factors: (a) the tube currents at the location of an organ and (b) the variation in organ volume along the longitudinal direction. This α-generic approach represents an approximation. The error associated with this approximation was evaluated using the α-specific organ dose (i.e., the organ dose obtained by using α-specific mA profiles as inputs into the Monte Carlo simulation) as the reference standard. When the mA profiles were constrained by a realistic system limit, this α-generic approach had errors of less than ~20% for the full range of α values. This was the case for 24 radiosensitive organs in both chest and abdomen-pelvis CT scans with the exception of thyroid in the chest scan and bladder in the abdomen-pelvis scan. For these two organs, the errors were less than ~40%. The results of this theoretical study suggested that knowing the mA modulation profile and the fixed-mA conversion factors, organ dose may be estimated for a TCM scan independent of the specific modulation scheme applied.

Background: Computer simulations are important for validating novel image acquisition and reconstruction strategies. In cardiovascular magnetic resonance (CMR), numerical simulations need to combine anatomical information and the effects of cardiac and/or respiratory motion. To this end, a framework for realistic CMR simulations is proposed and its use for image reconstruction from undersampled data is demonstrated., Methods: The extended Cardiac-Torso (XCAT) anatomical phantom framework with various motion options was used as a basis for the numerical phantoms. Different tissue, dynamic contrast and signal models, multiple receiver coils and noise are simulated. Arbitrary trajectories and undersampled acquisition can be selected. The utility of the framework is demonstrated for accelerated cine and first-pass myocardial perfusion imaging using k-t PCA and k-t SPARSE., Results: MRXCAT phantoms allow for realistic simulation of CMR including optional cardiac and respiratory motion. Example reconstructions from simulated undersampled k-t parallel imaging demonstrate the feasibility of simulated acquisition and reconstruction using the presented framework. Myocardial blood flow assessment from simulated myocardial perfusion images highlights the suitability of MRXCAT for quantitative post-processing simulation., Conclusion: The proposed MRXCAT phantom framework enables versatile and realistic simulations of CMR including breathhold and free-breathing acquisitions.

Purpose: This study aimed to provide a comprehensive patient-specific organ dose estimation across a multiplicity of computed tomography (CT) examination protocols., Methods: A validated Monte Carlo program was employed to model a common CT system (LightSpeed VCT, GE Healthcare). The organ and effective doses were estimated from 13 commonly used body and neurological CT examination. The dose estimation was performed on 58 adult computational extended cardiac-torso phantoms (35 male, 23 female, mean age 51.5 years, mean weight 80.2 kg). The organ dose normalized by CTDIvol (h factor) and effective dose normalized by the dose length product (DLP) (k factor) were calculated from the results. A mathematical model was derived for the correlation between the h and k factors with the patient size across the protocols. Based on this mathematical model, a dose estimation iPhone operating system application was designed and developed to be used as a tool to estimate dose to the patients for a variety of routinely used CT examinations., Results: The organ dose results across all the protocols showed an exponential decrease with patient body size. The correlation was generally strong for the organs which were fully or partially located inside the scan coverage (Pearson sample correlation coefficient (r) of 0.49). The correlation was weaker for organs outside the scan coverage for which distance between the organ and the irradiation area was a stronger predictor of dose to the organ. For body protocols, the effective dose before and after normalization by DLP decreased exponentially with increasing patient's body diameter (r > 0.85). The exponential relationship between effective dose and patient's body diameter was significantly weaker for neurological protocols (r < 0.41), where the trunk length was a slightly stronger predictor of effective dose (0.15 < r < 0.46)., Conclusions: While the most accurate estimation of a patient dose requires specific modeling of the patient anatomy, a first order approximation of organ and effective doses from routine CT scan protocols can be reasonably estimated using size specific factors. Estimation accuracy is generally poor for organ outside the scan range and for neurological protocols. The dose calculator designed in this study can be used to conveniently estimate and report the dose values for a patient across a multiplicity of CT scan protocols.

Mammography is currently the most widely utilized tool for detection and diagnosis of breast cancer. However, in women with dense breast tissue, tissue overlap may obscure lesions. Digital breast tomosynthesis can reduce tissue overlap. Furthermore, imaging with contrast enhancement can provide additional functional information about lesions, such as morphology and kinetics, which in turn may improve lesion identification and characterization. The performance of these imaging techniques is strongly dependent on the structural composition of the breast, which varies significantly among patients. Therefore, imaging system and imaging technique optimization should take patient variability into consideration. Furthermore, optimization of imaging techniques that employ contrast agents should include the temporally varying breast composition with respect to the contrast agent uptake kinetics. To these ends, we have developed a suite of 4-D virtual breast phantoms, which are incorporated with the kinetics of contrast agent propagation in different tissues and can realistically model normal breast parenchyma as well as benign and malignant lesions. This development presents a new approach in performing simulation studies using truly anthropomorphic models. To demonstrate the utility of the proposed 4-D phantoms, we present a simplified example study to compare the performance of 14 imaging paradigms qualitatively and quantitatively.

Purpose: Understanding the radiation dose to a patient is essential when considering the use of an ionizing diagnostic imaging test for clinical diagnosis and screening. Using Monte Carlo simulations, the authors estimated the three-dimensional organ-dose distribution from neutron and gamma irradiation of the male liver, female liver, and female breasts for neutron- and gamma-stimulated spectroscopic imaging., Methods: Monte Carlo simulations were developed using the Geant4 GATE application and a voxelized XCAT human phantom. A male and a female whole body XCAT phantom was voxelized into 256 × 256 × 600 voxels (3.125 × 3.125 × 3.125 mm(3)). A monoenergetic rectangular beam of 5.0 MeV neutrons or 7.0 MeV photons was made incident on a 2 cm thick slice of the phantom. The beam was rotated at eight different angles around the phantom ranging from 0° to 180°. Absorbed dose was calculated for each individual organ in the body and dose volume histograms were computed to analyze the absolute and relative doses in each organ., Results: The neutron irradiations of the liver showed the highest organ dose absorption in the liver, with appreciably lower doses in other proximal organs. The dose distribution within the irradiated slice exhibited substantial attenuation with increasing depth along the beam path, attenuating to ~15% of the maximum value at the beam exit side. The gamma irradiation of the liver imparted the highest organ dose to the stomach wall. The dose distribution from the gammas showed a region of dose buildup at the beam entrance, followed by a relatively uniform dose distribution to all of the deep tissue structures, attenuating to ~75% of the maximum value at the beam exit side. For the breast scans, both the neutron and gamma irradiation registered maximum organ doses in the breasts, with all other organs receiving less than 1% of the breast dose. Effective doses ranged from 0.22 to 0.37 mSv for the neutron scans and 41 to 66 mSv for the gamma scans., Conclusions: Neutron and gamma irradiation of a primary target organ was found to impart the majority of the total dose to the primary target organ (and other large organs) within the beam plane and considerably lower dose to proximal organs outside of the beam. These results also indicate that despite the use of a highly scattering particle such as a neutron, the dose from neutron stimulated emission computed tomography scans is on par with other clinical imaging techniques such as x-ray computed tomography (x-ray CT). Given the high nonuniformity in the dose across an organ during the neutron scan, care must be taken when computing average doses from neutron irradiations. The effective doses from neutron scanning were found to be comparable to x-ray CT. Further technique modifications are needed to reduce the effective dose levels from the gamma scans.

Purpose: The authors previously developed an adult population of 4D extended cardiac-torso (XCAT) phantoms for multimodality imaging research. In this work, the authors develop a reference set of 4D pediatric XCAT phantoms consisting of male and female anatomies at ages of newborn, 1, 5, 10, and 15 years. These models will serve as the foundation from which the authors will create a vast population of pediatric phantoms for optimizing pediatric CT imaging protocols., Methods: Each phantom was based on a unique set of CT data from a normal patient obtained from the Duke University database. The datasets were selected to best match the reference values for height and weight for the different ages and genders according to ICRP Publication 89. The major organs and structures were segmented from the CT data and used to create an initial pediatric model defined using nonuniform rational B-spline surfaces. The CT data covered the entire torso and part of the head. To complete the body, the authors manually added on the top of the head and the arms and legs using scaled versions of the XCAT adult models or additional models created from cadaver data. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from a template XCAT phantom (male or female 50th percentile adult) to the target pediatric model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. The masses of the organs in each phantom were matched to the reference values given in ICRP Publication 89. The new reference models were checked for anatomical accuracy via visual inspection., Results: The authors created a set of ten pediatric reference phantoms that have the same level of detail and functionality as the original XCAT phantom adults. Each consists of thousands of anatomical structures and includes parameterized models for the cardiac and respiratory motions. Based on patient data, the phantoms capture the anatomic variations of childhood, such as the development of bone in the skull, pelvis, and long bones, and the growth of the vertebrae and organs. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities., Conclusions: The development of patient-derived pediatric computational phantoms is useful in providing variable anatomies for simulation. Future work will expand this ten-phantom base to a host of pediatric phantoms representative of the public at large. This can provide a means to evaluate and improve pediatric imaging devices and to optimize CT protocols in terms of image quality and radiation dose.

Purpose: To estimate organ dose from pediatric chest and abdominopelvic computed tomography (CT) examinations and evaluate the dependency of organ dose coefficients on patient size and CT scanner models., Materials and Methods: The institutional review board approved this HIPAA-compliant study and did not require informed patient consent. A validated Monte Carlo program was used to perform simulations in 42 pediatric patient models (age range, 0-16 years; weight range, 2-80 kg; 24 boys, 18 girls). Multidetector CT scanners were modeled on those from two commercial manufacturers (LightSpeed VCT, GE Healthcare, Waukesha, Wis; SOMATOM Definition Flash, Siemens Healthcare, Forchheim, Germany). Organ doses were estimated for each patient model for routine chest and abdominopelvic examinations and were normalized by volume CT dose index (CTDI(vol)). The relationships between CTDI(vol)-normalized organ dose coefficients and average patient diameters were evaluated across scanner models., Results: For organs within the image coverage, CTDI(vol)-normalized organ dose coefficients largely showed a strong exponential relationship with the average patient diameter (R(2) > 0.9). The average percentage differences between the two scanner models were generally within 10%. For distributed organs and organs on the periphery of or outside the image coverage, the differences were generally larger (average, 3%-32%) mainly because of the effect of overranging., Conclusion: It is feasible to estimate patient-specific organ dose for a given examination with the knowledge of patient size and the CTDI(vol). These CTDI(vol)-normalized organ dose coefficients enable one to readily estimate patient-specific organ dose for pediatric patients in clinical settings. This dose information, and, as appropriate, attendant risk estimations, can provide more substantive information for the individual patient for both clinical and research applications and can yield more expansive information on dose profiles across patient populations within a practice., (© RSNA, 2013.)

Purpose: Given the radiation concerns inherent to the x-ray modalities, accurately estimating the radiation doses that patients receive during different imaging modalities is crucial. This study estimated organ doses, effective doses, and risk indices for the three clinical chest x-ray imaging techniques (chest radiography, tomosynthesis, and CT) using 59 anatomically variable voxelized phantoms and Monte Carlo simulation methods., Methods: A total of 59 computational anthropomorphic male and female extended cardiac-torso (XCAT) adult phantoms were used in this study. Organ doses and effective doses were estimated for a clinical radiography system with the capability of conducting chest radiography and tomosynthesis (Definium 8000, VolumeRAD, GE Healthcare) and a clinical CT system (LightSpeed VCT, GE Healthcare). A Monte Carlo dose simulation program (PENELOPE, version 2006, Universitat de Barcelona, Spain) was used to mimic these two clinical systems. The Duke University (Durham, NC) technique charts were used to determine the clinical techniques for the radiographic modalities. An exponential relationship between CTDIvol and patient diameter was used to determine the absolute dose values for CT. The simulations of the two clinical systems compute organ and tissue doses, which were then used to calculate effective dose and risk index. The calculation of the two dose metrics used the tissue weighting factors from ICRP Publication 103 and BEIR VII report., Results: The average effective dose of the chest posteroanterior examination was found to be 0.04 mSv, which was 1.3% that of the chest CT examination. The average effective dose of the chest tomosynthesis examination was found to be about ten times that of the chest posteroanterior examination and about 12% that of the chest CT examination. With increasing patient average chest diameter, both the effective dose and risk index for CT increased considerably in an exponential fashion, while these two dose metrics only increased slightly for radiographic modalities and for chest tomosynthesis. Effective and organ doses normalized to mAs all illustrated an exponential decrease with increasing patient size. As a surface organ, breast doses had less correlation with body size than that of lungs or liver., Conclusions: Patient body size has a much greater impact on radiation dose of chest CT examinations than chest radiography and tomosynthesis. The size of a patient should be considered when choosing the best thoracic imaging modality.

Recent studies have shown the feasibility of estimating patient dose from a CT exam using CTDI(vol)-normalized-organ dose (denoted as h), DLP-normalized-effective dose (denoted as k), and DLP-normalized-risk index (denoted as q). However, previous studies were limited to a small number of phantom models. The purpose of this work was to provide dose coefficients (h, k, and q) across a large number of computational models covering a broad range of patient anatomy, age, size percentile, and gender. The study consisted of 100 patient computer models (age range, 0 to 78 y.o.; weight range, 2-180 kg) including 42 pediatric models (age range, 0 to 16 y.o.; weight range, 2-80 kg) and 58 adult models (age range, 18 to 78 y.o.; weight range, 57-180 kg). Multi-detector array CT scanners from two commercial manufacturers (LightSpeed VCT, GE Healthcare; SOMATOM Definition Flash, Siemens Healthcare) were included. A previously-validated Monte Carlo program was used to simulate organ dose for each patient model and each scanner, from which h, k, and q were derived. The relationships between h, k, and q and patient characteristics (size, age, and gender) were ascertained. The differences in conversion coefficients across the scanners were further characterized. CTDI(vol)-normalized-organ dose (h) showed an exponential decrease with increasing patient size. For organs within the image coverage, the average differences of h across scanners were less than 15%. That value increased to 29% for organs on the periphery or outside the image coverage, and to 8% for distributed organs, respectively. The DLP-normalized-effective dose (k) decreased exponentially with increasing patient size. For a given gender, the DLP-normalized-risk index (q) showed an exponential decrease with both increasing patient size and patient age. The average differences in k and q across scanners were 8% and 10%, respectively. This study demonstrated that the knowledge of patient information and CTDIvol/DLP values may be used to estimate organ dose, effective dose, and risk index in abdominopelvic CT based on the coefficients derived from a large population of pediatric and adult patients.

Purpose: To develop a technique to estimate onboard 4D-CBCT using prior information and limited-angle projections for potential 4D target verification of lung radiotherapy., Methods: Each phase of onboard 4D-CBCT is considered as a deformation from one selected phase (prior volume) of the planning 4D-CT. The deformation field maps (DFMs) are solved using a motion modeling and free-form deformation (MM-FD) technique. In the MM-FD technique, the DFMs are estimated using a motion model which is extracted from planning 4D-CT based on principal component analysis (PCA). The motion model parameters are optimized by matching the digitally reconstructed radiographs of the deformed volumes to the limited-angle onboard projections (data fidelity constraint). Afterward, the estimated DFMs are fine-tuned using a FD model based on data fidelity constraint and deformation energy minimization. The 4D digital extended-cardiac-torso phantom was used to evaluate the MM-FD technique. A lung patient with a 30 mm diameter lesion was simulated with various anatomical and respirational changes from planning 4D-CT to onboard volume, including changes of respiration amplitude, lesion size and lesion average-position, and phase shift between lesion and body respiratory cycle. The lesions were contoured in both the estimated and "ground-truth" onboard 4D-CBCT for comparison. 3D volume percentage-difference (VPD) and center-of-mass shift (COMS) were calculated to evaluate the estimation accuracy of three techniques: MM-FD, MM-only, and FD-only. Different onboard projection acquisition scenarios and projection noise levels were simulated to investigate their effects on the estimation accuracy., Results: For all simulated patient and projection acquisition scenarios, the mean VPD (±S.D.)∕COMS (±S.D.) between lesions in prior images and "ground-truth" onboard images were 136.11% (±42.76%)∕15.5 mm (±3.9 mm). Using orthogonal-view 15°-each scan angle, the mean VPD∕COMS between the lesion in estimated and "ground-truth" onboard images for MM-only, FD-only, and MM-FD techniques were 60.10% (±27.17%)∕4.9 mm (±3.0 mm), 96.07% (±31.48%)∕12.1 mm (±3.9 mm) and 11.45% (±9.37%)∕1.3 mm (±1.3 mm), respectively. For orthogonal-view 30°-each scan angle, the corresponding results were 59.16% (±26.66%)∕4.9 mm (±3.0 mm), 75.98% (±27.21%)∕9.9 mm (±4.0 mm), and 5.22% (±2.12%)∕0.5 mm (±0.4 mm). For single-view scan angles of 3°, 30°, and 60°, the results for MM-FD technique were 32.77% (±17.87%)∕3.2 mm (±2.2 mm), 24.57% (±18.18%)∕2.9 mm (±2.0 mm), and 10.48% (±9.50%)∕1.1 mm (±1.3 mm), respectively. For projection angular-sampling-intervals of 0.6°, 1.2°, and 2.5° with the orthogonal-view 30°-each scan angle, the MM-FD technique generated similar VPD (maximum deviation 2.91%) and COMS (maximum deviation 0.6 mm), while sparser sampling yielded larger VPD∕COMS. With equal number of projections, the estimation results using scattered 360° scan angle were slightly better than those using orthogonal-view 30°-each scan angle. The estimation accuracy of MM-FD technique declined as noise level increased., Conclusions: The MM-FD technique substantially improves the estimation accuracy for onboard 4D-CBCT using prior planning 4D-CT and limited-angle projections, compared to the MM-only and FD-only techniques. It can potentially be used for the inter/intrafractional 4D-localization verification.

A novel technique based on Fourier transform theory has been developed that directly extracts respiratory information from projections without the use of external surrogates. While the feasibility has been demonstrated with three patients, a more extensive validation is necessary. Therefore, the purpose of this work is to investigate the effects of a variety of respiratory and anatomical scenarios on the performance of the technique with the 4D digital extended cardiac torso phantom. FT-phase and FT-magnitude methods were each applied to identify peak-inspiration projections and quantitatively compared to the gold standard of visual identification. Both methods proved to be robust across the studied scenarios with average differences in respiratory phase <10% and percentage of projections assigned within 10% of the gold standard >90%, when incorporating minor modifications to region-of-interest (ROI) selection and/or low-frequency location for select cases of DA and lung percentage in the field of view of the projection. Nevertheless, in the instance where one method initially faltered, the other method prevailed and successfully identified peak-inspiration projections. This is promising because it suggests that the two methods provide complementary information to each other. To ensure appropriate clinical adaptation of markerless, self-sorted four-dimensional cone-beam CT (4D-CBCT), perhaps an optimal integration of the two methods can be developed.

High-order cubic Hermite finite elements have been valuable in modeling cardiac geometry, fiber orientations, biomechanics, and electrophysiology, but their use in solving three-dimensional problems has been limited to ventricular models with simple topologies. Here, we utilized a subdivision surface scheme and derived a generalization of the "local-to-global" derivative mapping scheme of cubic Hermite finite elements to construct bicubic and tricubic Hermite models of the human atria with extraordinary vertices from computed tomography images of a patient with atrial fibrillation. To an accuracy of 0.6 mm, we were able to capture the left atrial geometry with only 142 bicubic Hermite finite elements, and the right atrial geometry with only 90. The left and right atrial bicubic Hermite meshes were G1 continuous everywhere except in the one-neighborhood of extraordinary vertices, where the mean dot products of normals at adjacent elements were 0.928 and 0.925. We also constructed two biatrial tricubic Hermite models and defined fiber orientation fields in agreement with diagrammatic data from the literature using only 42 angle parameters. The meshes all have good quality metrics, uniform element sizes, and elements with aspect ratios near unity, and are shared with the public. These new methods will allow for more compact and efficient patient-specific models of human atrial and whole heart physiology., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Purpose: The authors previously developed the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. The XCAT consisted of highly detailed whole-body models for the standard male and female adult, including the cardiac and respiratory motions. In this work, the authors extend the XCAT beyond these reference anatomies by developing a series of anatomically variable 4D XCAT adult phantoms for imaging research, the first library of 4D computational phantoms., Methods: The initial anatomy of each phantom was based on chest-abdomen-pelvis computed tomography data from normal patients obtained from the Duke University database. The major organs and structures for each phantom were segmented from the corresponding data and defined using nonuniform rational B-spline surfaces. To complete the body, the authors manually added on the head, arms, and legs using the original XCAT adult male and female anatomies. The structures were scaled to best match the age and anatomy of the patient. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from the template XCAT phantom (male or female) to the target patient model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. Each new phantom was refined by checking for anatomical accuracy via inspection of the models., Results: Using these methods, the authors created a series of computerized phantoms with thousands of anatomical structures and modeling cardiac and respiratory motions. The database consists of 58 (35 male and 23 female) anatomically variable phantoms in total. Like the original XCAT, these phantoms can be combined with existing simulation packages to simulate realistic imaging data. Each new phantom contains parameterized models for the anatomy and the cardiac and respiratory motions and can, therefore, serve as a jumping point from which to create an unlimited number of 3D and 4D variations for imaging research., Conclusions: A population of phantoms that includes a range of anatomical variations representative of the public at large is needed to more closely mimic a clinical study or trial. The series of anatomically variable phantoms developed in this work provide a valuable resource for investigating 3D and 4D imaging devices and the effects of anatomy and motion in imaging. Combined with Monte Carlo simulation programs, the phantoms also provide a valuable tool to investigate patient-specific dose and image quality, and optimization for adults undergoing imaging procedures.

Purpose: The authors previously reported on a three-dimensional computer-generated breast phantom, based on empirical human image data, including a realistic finite-element based compression model that was capable of simulating multimodality imaging data. The computerized breast phantoms are a hybrid of two phantom generation techniques, combining empirical breast CT (bCT) data with flexible computer graphics techniques. However, to date, these phantoms have been based on single human subjects. In this paper, the authors report on a new method to generate multiple phantoms, simulating additional subjects from the limited set of original dedicated breast CT data. The authors developed an image morphing technique to construct new phantoms by gradually transitioning between two human subject datasets, with the potential to generate hundreds of additional pseudoindependent phantoms from the limited bCT cases. The authors conducted a preliminary subjective assessment with a limited number of observers (n = 4) to illustrate how realistic the simulated images generated with the pseudoindependent phantoms appeared., Methods: Several mesh-based geometric transformations were developed to generate distorted breast datasets from the original human subject data. Segmented bCT data from two different human subjects were used as the "base" and "target" for morphing. Several combinations of transformations were applied to morph between the "base' and "target" datasets such as changing the breast shape, rotating the glandular data, and changing the distribution of the glandular tissue. Following the morphing, regions of skin and fat were assigned to the morphed dataset in order to appropriately assign mechanical properties during the compression simulation. The resulting morphed breast was compressed using a finite element algorithm and simulated mammograms were generated using techniques described previously. Sixty-two simulated mammograms, generated from morphing three human subject datasets, were used in a preliminary observer evaluation where four board certified breast radiologists with varying amounts of experience ranked the level of realism (from 1 = "fake" to 10 = "real") of the simulated images., Results: The morphing technique was able to successfully generate new and unique morphed datasets from the original human subject data. The radiologists evaluated the realism of simulated mammograms generated from the morphed and unmorphed human subject datasets and scored the realism with an average ranking of 5.87 ± 1.99, confirming that overall the phantom image datasets appeared more "real" than "fake." Moreover, there was not a significant difference (p > 0.1) between the realism of the unmorphed datasets (6.0 ± 1.95) compared to the morphed datasets (5.86 ± 1.99). Three of the four observers had overall average rankings of 6.89 ± 0.89, 6.9 ± 1.24, 6.76 ± 1.22, whereas the fourth observer ranked them noticeably lower at 2.94 ± 0.7., Conclusions: This work presents a technique that can be used to generate a suite of realistic computerized breast phantoms from a limited number of human subjects. This suite of flexible breast phantoms can be used for multimodality imaging research to provide a known truth while concurrently producing realistic simulated imaging data.

Unlabelled: A new generation of reference computational phantoms, based on image-based models tied to the reference masses defined by the International Commission on Radiological Protection (ICRP) for dose calculations, is presented., Methods: Anatomic models based on nonuniform rational b-spline modeling techniques were used to define reference male and female adults, 15-y-olds, 10-y-olds, 5-y-olds, 1-y-olds, newborns, and pregnant women at 3 stages of gestation, using the defined reference organ masses in ICRP publication 89. Absorbed fractions and specific absorbed fractions for internal emitters were derived using standard Monte Carlo radiation transport simulation codes., Results: Differences were notable between many pairs of organs in specific absorbed fractions because of the improved realism of the models, with adjacent organs usually closer and sometimes touching. Final estimates of absorbed dose for radiopharmaceuticals, for example, were only slightly different overall, as many of the differences were small and most pronounced at low radiation energies. Some new important organs were defined (salivary glands, prostate, eyes, and esophagus), and the identity of a few gastrointestinal tract organs changed., Conclusion: A new generation of reference models for standardized internal and external dose calculations has been defined. The models will be implemented in standardized software for internal dose calculations and be used to produce new standardized dose estimates for radiopharmaceuticals and other applications.

Purpose: In computed tomography (CT), organ dose, effective dose, and risk index can be estimated from volume-weighted CT dose index (CTDI(vol)) or dose-length product (DLP) using conversion coefficients. Studies have investigated how these coefficients vary across scanner models, scan parameters, and patient size. However, their variability across CT protocols has not been systematically studied. Furthermore, earlier studies of the effect of patient size have not included obese individuals, which currently represent more than one-third of U.S. adults. The purpose of this study was to assess the effects of protocol and obesity on dose and risk conversion coefficients in adult body CT., Methods: Whole-body computational phantoms were created from clinical CT images of six adult patients (three males, three females), representing normal-weight patients and patients of three obesity classes. Body CT protocols at our institution were selected and categorized into ten examination categories based on anatomical region examined. A validated Monte Carlo program was used to estimate organ dose. Organ dose estimates were normalized by CTDI(vol) and size-specific dose estimate (SSDE) to obtain organ dose conversion coefficients (denoted as h and h(ss) factors, respectively). Assuming each phantom to be 20, 40, and 60 years old, effective dose and risk index were calculated and normalized by DLP to obtain effective dose and risk index conversion coefficients (denoted as k and q factors, respectively). Coefficient of variation was used to quantify the variability of each conversion coefficient across examination categories. The effect of obesity was assessed by comparing each obese phantom with the normal-weight phantom of the same gender., Results: For a given organ, the variability of h factor across examination categories that encompassed the entire organ volume was generally within 15%. However, k factor varied more across examination categories (15%-27%). For all three ages, the variability of q factor was small for male (<10%), but large for female phantoms (21%-43%). Relative to the normal-weight phantoms, the reduction in h factor (an average across fully encompassed organs) was 17%-42%, 17%-40%, and 51%-63% for obese-class-I, obese-class-II, and obese-class-III phantoms, respectively. h(ss) factor was not independent of patient diameter and generally decreased with increasing obesity. Relative to the normal-weight phantoms, the reduction in k factor was 12%-40%, 14%-46%, and 44%-59% for obese-class-I, obese-class-II, and obese-class-III phantoms, respectively. The respective reduction in q factor was 11%-36%, 17%-42%, and 48%-59% at 20 years of age and similar at other ages., Conclusions: In adult body CT, dose to an organ fully encompassed by the primary radiation beam can be estimated from CTDI(vol) using a protocol-independent conversion coefficient. However, fully encompassed organs only account for 50% ± 19% of k factor and 46% ± 24% of q factor. Dose received by partially encompassed organs is also substantial. To estimate effective dose and risk index from DLP, it is necessary to use conversion coefficients specific to the anatomical region examined. Obesity has a significant effect on dose and risk conversion coefficients, which cannot be predicted using body diameter alone. SSDE-normalized organ dose is not independent of diameter. SSDE itself generally overestimates organ dose for obese patients.