Search

Background: Dynamic chest radiography (DCR) is a recently developed functional x-ray imaging technique that detects pulmonary ventilation impairment as a decrease in changes in lung density during respiration. However, the diagnostic performance of DCR is uncertain owing to an insufficient number of clinical cases. One solution is virtual imaging trials (VITs), which is an emerging alternative method for efficiently evaluating medical imaging technology via computer simulation techniques., Purpose: This study aimed to estimate the typical threshold thickness of residual normal tissue below which the presence of emphysema may be detected by DCR via VITs using virtual patients with different physiques and a user-defined ground truth., Methods: Twenty extended cardiac-torso (XCAT) phantoms that exhibited changes in lung density during respiration were generated to simulate virtual patients. To simulate a locally collapsed lung, an air sphere was inserted into each lung regions in the phantom. The XCAT phantom was virtually projected using an x-ray simulator. The respiratory changes in pixel value (ΔPV) were measured on the projected air spheres (simulated lesions) to calculate the percentage of decrease (ΔPV%) relative to ΔPV exp-ins in the absence of an air sphere. The relationship between the amount of residual normal tissue and ΔPV% was fitted to a cubic approximation curve (hereafter, performance curve), and the threshold at which the ΔPV% began to decrease (normal-tissue thre ) was determined. The goodness of fit for each performance curve was evaluated according to the coefficient of determination (R 2 ) and the 95% confidence interval derived from the standard errors between the measured and theoretical values corresponding to each performance curve. The ΔPV% was also visualized as a color scaling to validate the results of the VITs in both virtual and clinical patients., Results: For each lung region in all body sizes, the ΔPV% decreased as the amount of residual normal tissue decreased and could be defined as a function of the amount of residual normal tissue in front of and behind the simulated lesions with high R 2 values. Meanwhile, the difference between the measured and theoretical values corresponding to each performance curve was only partially included in the 95% confidence interval. The normal-tissue thre values were 146.0, 179.5, and 170.9 mm for the upper, middle, and lower lungs, respectively, which were demonstrated in virtual patients and one real patient, where the value of the residual normal tissue was less than that of normal-tissue thre ; any reduction in the residual normal tissue was reflected as a reduced ΔPV and depicted as a reduced color intensity., Conclusions: The performance of DCR-based pulmonary impairment assessment depends on the amount of residual normal tissue in front of and behind the lesion rather than on the lesion size. The performance curve can be defined as a function of the amount of residual normal tissue in each lung region with a specific threshold of normal tissue remaining where lesions become detectable, shown as a decrease in ΔPV. The results of VITs are expected to accelerate future clinical trials for DCR-based pulmonary function assessment., (© 2024 The Author(s). Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Background: Photon-counting computed tomography (PCCT) has recently emerged into clinical use; however, its optimum imaging protocols and added benefits remains unknown in terms of providing more accurate lung density quantification compared to energy-integrating computed tomography (EICT) scanners., Purpose: To systematically assess the performance of a clinical PCCT scanner for lung density quantifications and compare it against EICT., Methods: This cross-sectional study involved a retrospective analysis of subjects scanned (August-December 2021) using a clinical PCCT system. The influence of altering reconstruction parameters was studied (reconstruction kernel, pixel size, slice thickness). A virtual CT dataset of anthropomorphic virtual subjects was acquired to demonstrate the correspondence of findings to clinical dataset, and to perform systematic imaging experiments, not possible using human subjects. The virtual subjects were imaged using a validated, scanner-specific CT simulator of a PCCT and two EICT (defined as EICT A and B) scanners. The images were evaluated using mean absolute error (MAE) of lung and emphysema density against their corresponding ground truth., Results: Clinical and virtual PCCT datasets showed similar trends, with sharper kernels and smaller voxel sizes increasing percentage of low-attenuation areas below -950 HU (LAA-950) by up to 15.7 ± 6.9% and 11.8 ± 5.5%, respectively. Under the conditions studied, higher doses, thinner slices, smaller pixel sizes, iterative reconstructions, and quantitative kernels with medium sharpness resulted in lower lung MAE values. While using these settings for PCCT, changes in the dose level (13 to 1.3 mGy), slice thickness (0.4 to 1.5 mm), pixel size (0.49 to 0.98 mm), reconstruction technique (70 keV-VMI to wFBP), and kernel (Qr48 to Qr60) increased lung MAE by 15.3 ± 2.0, 1.4 ± 0.6, 2.2 ± 0.3, 4.2 ± 0.8, and 9.1 ± 1.6 HU, respectively. At the optimum settings identified per scanner, PCCT images exhibited lower lung and emphysema MAE than those of EICT scanners (by 2.6 ± 1.0 and 9.6 ± 3.4 HU, compared to EICT A, and by 4.8 ± 0.8 and 7.4 ± 2.3 HU, compared to EICT B). The accuracy of lung density measurements was correlated with subjects' mean lung density (p < 0.05), measured by PCCT at optimum setting under the conditions studied., Conclusion: Photon-counting CT demonstrated superior performance in density quantifications, with its influences of imaging parameters in line with energy-integrating CT scanners. The technology offers improvement in lung quantifications, thus demonstrating potential toward more objective assessment of respiratory conditions., (© 2024 American Association of Physicists in Medicine.)

Background: As a leading cause of death, worldwide, cardiovascular disease is of great clinical importance. Among cardiovascular diseases, coronary artery disease (CAD) is a key contributor, and it is the attributed cause of death for 10% of all deaths annually. The prevalence of CAD is commensurate with the rise in new medical imaging technologies intended to aid in its diagnosis and treatment. The necessary clinical trials required to validate and optimize these technologies require a large cohort of carefully controlled patients, considerable time to complete, and can be prohibitively expensive. A safer, faster, less expensive alternative is using virtual imaging trials (VITs), utilizing virtual patients or phantoms combined with accurate computer models of imaging devices., Purpose: In this work, we develop realistic, physiologically-informed models for coronary plaques for application in cardiac imaging VITs., Methods: Histology images of plaques at micron-level resolution were used to train a deep convolutional generative adversarial network (DC-GAN) to create a library of anatomically variable plaque models with clinical anatomical realism. The stability of each plaque was evaluated by finite element analysis (FEA) in which plaque components and vessels were meshed as volumes, modeled as specialized tissues, and subjected to the range of normal coronary blood pressures. To demonstrate the utility of the plaque models, we combined them with the whole-body XCAT computational phantom to perform initial simulations comparing standard energy-integrating detector (EID) CT with photon-counting detector (PCD) CT., Results: Our results show the network is capable of generating realistic, anatomically variable plaques. Our simulation results provide an initial demonstration of the utility of the generated plaque models as targets to compare different imaging devices., Conclusions: Vast, realistic, and variable CAD pathologies can be generated to incorporate into computational phantoms for VITs. There they can serve as a known truth from which to optimize and evaluate cardiac imaging technologies quantitatively., (© 2024 American Association of Physicists in Medicine.)

Background: Digital anthropomorphic phantoms, such as the 4D extended cardiac-torso (XCAT) phantom, are actively used to develop, optimize, and evaluate a variety of imaging applications, allowing for realistic patient modeling and knowledge of ground truth. The XCAT phantom defines the activity and attenuation for a simulated patient, which includes a complete set of organs, muscle, bone, and soft tissue, while also accounting for cardiac and respiratory motion. However, the XCAT phantom does not currently include the lymphatic system, critical for evaluating medical imaging tasks such as sentinel node detection, node density measurement, and radiation dosimetry., Purpose: In this study, we aimed to develop a scalable lymphatic system in the XCAT phantom, to facilitate improved research of the lymphatic system in medical imaging. Using this scalable lymphatic system, we modeled the lymph node conglomerate pathology that is characteristically observed in primary mediastinal B-cell lymphoma (PMBCL). As an extended application, we evaluated positron emission tomography (PET) image quantification of metabolic tumor volume (MTV) and total lesion glycolysis (TLG) of these simulated lymphomas, though the phantoms may be applied to other imaging modalities and study design paradigms (e.g., image quality, detection)., Methods: A template model for the lymphatic system was developed based on anatomical data from the Visible Human Project of the National Library of Medicine. The segmented nodes and vessels were fit with non-uniform rational basis spline surfaces, and multichannel large deformation diffeomorphic metric mapping was used to propagate the template to different XCAT anatomies. To model conglomerates observed in PMBCL, lymph nodes were enlarged, converged within the mediastinum, and tracer concentration was increased. We used the phantoms as inputs to a PET simulation tool, which generated images using ordered subsets expectation maximization reconstruction with 2-8 mm Gaussian filters. Fixed thresholding (FT) and gradient segmentation were used to determine MTV and TLG. Percent bias (%Bias) and coefficient of variation (COV) were computed as measures of accuracy and precision, respectively, for each MTV and TLG measurement., Results: Using the methodology described above, we introduced a scalable lymphatic system in the XCAT phantom, which allows for the radioactivity and attenuation ground truth to be generated in 116 ± 2.5 s using a 2.3 GHz processor. Within the Rhinoceros interface, lymph node anatomy and function were modified to create a cohort of 10 phantoms with lymph node conglomerates. Using the lymphoma phantoms to evaluate PET quantification of MTV, mean %Bias values were -9.3%, -41.3%, and 20.9%, while COV values were 4.08%, 7.6%, and 3.4% using 25% FT, 40% FT, and gradient segmentations, respectively. Comparatively for TLG, mean %Bias values were -27.4%, -45.8%, and -16.0%, while COV values were 1.9%, 5.7%, and 1.4%, for the 25% FT, 40% FT, and gradient segmentations, respectively., Conclusions: In this work, we upgraded the XCAT phantom to include a lymphatic system, comprised of a network of 276 scalable lymph nodes and corresponding vessels. As an application, we created a cohort of phantoms with lymph node conglomerates to evaluate lymphoma quantification in PET imaging, which highlights an important application of this work., (© 2022 American Association of Physicists in Medicine.)

Purpose: We aim at developing a model-based algorithm that compensates for the effect of both pulse pileup (PP) and charge sharing (CS) and evaluates the performance using computer simulations., Methods: The proposed PCP algorithm for PP and CS compensation uses cascaded models for CS and PP we previously developed, maximizes Poisson log-likelihood, and uses an efficient three-step exhaustive search. For comparison, we also developed an LCP algorithm that combines models for a loss of counts (LCs) and CS. Two types of computer simulations, slab- and computed tomography (CT)-based, were performed to assess the performance of both PCP and LCP with 200 and 800 mA, (300 µm) 2  × 1.6-mm cadmium telluride detector, and a dead-time of 23 ns. A slab-based assessment used a pair of adipose and iodine with different thicknesses, attenuated X-rays, and assessed the bias and noise of the outputs from one detector pixel; a CT-based assessment simulated a chest/cardiac scan and a head-and-neck scan using 3D phantom and noisy cone-beam projections., Results: With the slab simulation, the PCP had little or no biases when the expected counts were sufficiently large, even though a probability of count loss (PCL) due to dead-time loss or PP was as high as 0.8. In contrast, the LCP had significant biases (>±2 cm of adipose) when the PCL was higher than 0.15. Biases were present with both PCP and LCP when the expected counts were less than 10-120 per datum, which was attributed to the fact that the maximum likelihood did not approach the asymptote. The noise of PCP was within 8% from the Cramér-Rao lower bounds for most cases when no significant bias was present. The two CT studies essentially agreed with the slab simulation study. PCP had little or no biases in the estimated basis line integrals, reconstructed basis density maps, and synthesized monoenergetic CT images. But the LCP had significant biases in basis line integrals when X-ray beams passed through lungs and near the body and neck contours, where the PCLs were above 0.15. As a consequence, basis density maps and monoenergetic CT images obtained by LCP had biases throughout the imaged space., Conclusion: We have developed the PCP algorithm that uses the PP-CS model. When the expected counts are more than 10-120 per datum, the PCP algorithm is statistically efficient and successfully compensates for the effect of the spectral distortion due to both PP and CS providing little or no biases in basis line integrals, basis density maps, and monoenergetic CT images regardless of count-rates. In contrast, the LCP algorithm, which models an LC due to pileup, produces severe biases when incident count-rates are high and the PCL is 0.15 or higher., (© 2022 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Purpose: The gold-standard method for estimation of patient-specific organ doses in digital tomosynthesis (DT) requires protocol-specific Monte Carlo (MC) simulations of radiation transport in anatomically accurate computational phantoms. Although accurate, MC simulations are computationally expensive, leading to a turnaround time in the order of core hours for simulating a single exam. This limits their clinical utility. The purpose of this study is to overcome this limitation by utilizing patient- and protocol-specific MC simulations to develop a comprehensive database of air-kerma-normalized organ dose coefficients for a virtual population of adult and pediatric patient models over an expanded set of exam protocols in DT for retrospective and prospective estimation of radiation dose in clinical tomosynthesis., Materials and Methods: A clinically representative virtual population of 14 patient models was used, with pediatric models (M and F) at ages 1, 5, 10, and 15 and adult patient models (M and F) with body mass index (BMIs) at 10th, 50th, and 90th percentiles of the US population. A graphics processing unit (GPU)-based MC simulation framework was used to simulate organ doses in the patient models, incorporating the scanner-specific configuration of a clinical DT system (VolumeRad, GE Healthcare, Waukesha, WI, USA) and an expanded set of exam protocols, including 21 distinct acquisition techniques for imaging a variety of anatomical regions (head and neck, thorax, spine, abdomen, and knee). Organ dose coefficients (h n ) were estimated by normalizing organ dose estimates to air kerma at 70 cm (X 70cm ) from the source in the scout view. The corresponding coefficients for projection radiography were approximated using organ doses estimated for the scout view. The organ dose coefficients were further used to compute air-kerma-normalized patient-specific effective dose coefficients (K n ) for all combinations of patients and protocols, and a comparative analysis examining the variation of radiation burden across sex, age, and exam protocols in DT, and with projection radiography was performed., Results: The database of organ dose coefficients (h n ) containing 294 distinct combinations of patients and exam protocols was developed and made publicly available. The values of K n were observed to produce estimates of effective dose in agreement with prior studies and consistent with magnitudes expected for pediatric and adult patients across the different exam protocols, with head and neck regions exhibiting relatively lower and thorax and C-spine (apsc, apcs) regions relatively higher magnitudes. The ratios (r = K n /K n ,rad ) quantifying the differences air-kerma-normalized patient-specific effective doses between DT and projection radiography were centered around 1.0 for all exam protocols, with the exception of protocols covering the knee region (pawk, patk)., Conclusions: This study developed a database of organ dose coefficients for a virtual population of 14 adult and pediatric XCAT patient models over a set of 21 exam protocols in DT. Using empirical measurements of air kerma in the clinic, these organ dose coefficients enable practical retrospective and prospective patient-specific radiation dosimetry. The computation of air-kerma-normalized patient-specific effective doses further enables the comparison of radiation burden to the patient populations between protocols and between imaging modalities (e.g., DT and projection radiography), as presented in this study., (© 2022 American Association of Physicists in Medicine.)

Purpose: Estimation of organ doses in digital tomosynthesis (DT) is challenging due to the lack of existing tools that accurately and flexibly model protocol- and view-specific collimations and motion trajectories of the source and detector for a variety of exam protocols, and the computational inefficiencies of conducting MC simulations. The purpose of this study was to overcome these limitations by developing and benchmarking a GPU-accelerated MC simulation framework compatible with patient-specific computational phantoms for individualized estimation of organ doses in DT., Materials and Methods: The framework for individualized estimation of dose in DT was developed as a two-step workflow: (1) a custom MATLAB code that accepts a patient-specific computational phantom and exam description (organ markers for defining the extremities of the anatomical region of interest, tube voltage, source-to-image distance, angular sweep range, number of projection views, and the pivot point to image distance - PPID) to compute the field of views (FOVs) for a clinical DT system, and (2) a MC tool (developed using MC-GPU) modeling the configuration of a clinical DT system to estimate organ doses based on the computed FOVs. Using this framework, we estimated organ doses for 28 radiosensitive organs in an adult reference patient model (M; 30 years) imaged using a commercial DT system (VolumeRad, GE Healthcare, Waukesha, WI). The estimates were benchmarked against values from a comparable organ dose estimation framework (reference dataset developed by the Advanced Laboratory for Radiation Dosimetry Studies at University of Florida) for a posterior-anterior chest exam. The resulting differences were quantified as percent relative errors and analyzed to identify any potential sources of bias and uncertainties. The timing performance (run duration in seconds) of the framework was also quantified for the same simulation to gauge the feasibility of the workflow for time-constrained clinical applications., Results: The organ dose estimates from the developed framework showed a close agreement with the reference dataset, with percent relative errors ranging from -6.9% to 5.0% and a mean absolute percent difference of 1.7% over all radiosensitive organs, with the exception of testes and eye lens, for which the percent relative errors were higher at -18.9% and -27.6%, respectively, due to their relative positioning outside the primary irradiation field, leading to fewer photons depositing energy and consequently higher errors in estimated organ doses. The run duration for the same simulation was 916.3 s, representing a substantial improvement in performance over existing nonparallelized MC tools., Conclusions: This study successfully developed and benchmarked a GPU-accelerated framework compatible with patient-specific anthropomorphic computational phantoms for accurate individualized estimation of organ doses in DT. By enabling patient-specific estimation of organ doses, this framework can aid clinicians and researchers by providing them with tools essential for tracking the radiation burden to patients for dose monitoring purposes and identifying the trends and relationships in organ doses for a patient population to optimize existing and develop new exam protocols., (© 2021 American Association of Physicists in Medicine.)

Purpose: Physicians utilize cerebral perfusion maps (e.g., cerebral blood flow, cerebral blood volume, transit time) to prescribe the plan of care for stroke patients. Variability in scanning techniques and post-processing software can result in differences between these perfusion maps. To determine which techniques are acceptable for clinical care, it is important to validate the accuracy and reproducibility of the perfusion maps. Validation using clinical data is challenging due to the lack of a gold standard to assess cerebral perfusion and the impracticality of scanning patients multiple times with different scanning techniques. In contrast, simulated data from a realistic digital phantom of the cerebral perfusion in acute stroke patients would enable studies to optimize and validate the scanning and post-processing techniques., Methods: We describe a complete framework to simulate CT perfusion studies for stroke assessment. We begin by expanding the XCAT brain phantom to enable spatially varying contrast agent dynamics and incorporate a realistic model of the dynamics in the cerebral vasculature derived from first principles. A dynamic CT simulator utilizes the time-concentration curves to define the contrast agent concentration in the object at each time point and generates CT perfusion images compatible with commercially available post-processing software. We also generate ground truth perfusion maps to which the maps generated by post-processing software can be compared., Results: We demonstrate a dynamic CT perfusion study of a simulated patient with an ischemic stroke and the resulting perfusion maps generated by post-processing software. We include a visual comparison between the computer-generated perfusion maps and the ground truth perfusion maps. The framework is highly tunable; users can modify the perfusion properties (e.g., occlusion location, CBF, CBV, and MTT), scanner specifications (e.g., focal spot size and detector configuration), scanning protocol (e.g., kVp and mAs), and reconstruction parameters (e.g., slice thickness and reconstruction filter)., Conclusions: This framework provides realistic test data with the underlying ground truth that enables a robust assessment of CT perfusion techniques and post-processing methods for stroke assessment., (© 2021 American Association of Physicists in Medicine.)

Purpose: Accurate preoperative assessment of tumor invasion/adhesion is crucial for planning appropriate operative procedures. Recent advances in digital radiography allow a motion analysis of lung tumors with dynamic chest radiography (DCR) with total exposure dose comparable to that of conventional chest radiography. The aim of this study was to investigate the feasibility of preoperative evaluation of pleural invasion/adhesion of lung tumors with DCR through a virtual clinical imaging study, using a four-dimensional (4D) extended cardiac-torso (XCAT) computational phantom., Methods: An XCAT phantom of an adult man (50th percentile in height and weight) with simulated respiratory and cardiac motions was generated to use as a virtual patient. To simulate lung tumors with and without pleural invasion, a 30-mm diameter tumor sphere was inserted into each lobe of the phantom. The virtual patient during respiration was virtually projected using an x-ray simulator in posteroanterior (PA) and oblique directions, and sequential bone suppression (BS) images were created. The measurement points (tumor, rib, and diaphragm) were automatically tracked on simulated images by a template matching technique. We calculated five quantitative metrics related to the movement distance and directions of the targeted tumor and evaluated whether DCR could distinguish between tumors with and without pleural invasion/adhesion., Results: Precise tracking of the targeted tumor was achieved on the simulated BS images without undue influence of rib shadows. There was a significant difference in all five quantitative metrics between the lung tumors with and without pleural invasion both on the oblique and PA projection views (P < 0.05). Quantitative metrics related to the movement distance were effective for tumors in the middle and lower lobes, while, those related to the movement directions were effective for tumors close to the frontal chest wall on the oblique projection view. The oblique views were useful for the evaluation of the space between the chest wall and a moving tumor., Conclusion: DCR could help distinguish between tumors with and without pleural invasion/adhesion based on the two-dimensional movement distance and direction using oblique and PA projection views. With anticipated improved image: processing to evaluate the respiratory displacement of lung tumors in the upper lobe or behind the heart, DCR holds promise for clinical assessment of tumor invasion/adhesion in the parietal pleura., (© 2021 American Association of Physicists in Medicine.)

Purpose: Many interventional procedures aim at changing soft tissue perfusion or blood flow. One problem at present is that soft tissue perfusion and its changes cannot be assessed in an interventional suite because cone-beam computed tomography is too slow (it takes 4-10 s per volume scan). In order to address the problem, we propose a novel method called IPEN for Intra-operative four-dimensional soft tissue PErfusion using a standard x-ray system with No gantry rotation., Methods: IPEN uses two input datasets: (a) the contours and locations of three-dimensional regions-of-interest (ROIs) such as arteries and sub-sections of cancerous lesions, and (b) a series of x-ray projection data obtained from an intra-arterial contrast injection to contrast enhancement to wash-out. IPEN then estimates a time-enhancement curve (TEC) for each ROI directly from projections without reconstructing cross-sectional images by maximizing the agreement between synthesized and measured projections with a temporal roughness penalty. When path lengths through ROIs are known for each x-ray beam, the ROI-specific enhancement can be accurately estimated from projections. Computer simulations are performed to assess the performance of the IPEN algorithm. Intra-arterial contrast-enhanced liver scans over 25 s were simulated using XCAT phantom version 2.0 with heterogeneous tissue textures and cancerous lesions. The following four sub-studies were performed: (a) The accuracy of the estimated TECs with overlapped lesions was evaluated at various noise (dose) levels with either homogeneous or heterogeneous lesion enhancement patterns; (b) the accuracy of IPEN with inaccurate ROI contours was assessed; (c) we investigated how overlapping ROIs and noise in projections affected the accuracy of the IPEN algorithm; and (d) the accuracy of the perfusion indices was assessed., Results: The TECs estimated by IPEN were sufficiently accurate at a reference dose level with the root-mean-square deviation (RMSD) of 0.0027 ± 0.0001 cm -1 or 13 ± 1 Hounsfield unit (mean ± standard deviation) for the homogeneous lesion enhancement and 0.0032 ± 0.0005 cm -1 for the heterogeneous enhancement (N = 20 each). The accuracy was degraded with decreasing doses: The RMSD with homogeneous enhancement was 0.0220 ± 0.0003 cm -1 for 20% of the reference dose level. Performing 3 × 3 pixel averaging on projection data improved the RMSDs to 0.0051 ± 0.0002 cm -1 for 20% dose. When the ROI contours were inaccurate, smaller ROI contours resulted in positive biases in TECs, whereas larger ROI contours produced negative biases. The bias remained small, within ± 0.0070 cm -1 , when the Sorenson-Dice coefficients (SDCs) were larger than 0.81. The RMSD of the TEC estimation was strongly associated with the condition of the problem, which can be empirically quantified using the condition number of a matrix A z that maps a vector of ROI enhancement values z to projection data and a weighted variance of projection data: a linear correlation coefficient (R) was 0.794 (P < 0.001). The perfusion index values computed from the estimated TECs agreed well with the true values (R ≥ 0.985, P < 0.0001)., Conclusion: The IPEN algorithm can estimate ROI-specific TECs with high accuracy especially when 3 × 3 pixel averaging is applied, even when lesion enhancement is heterogeneous, or ROI contours are inaccurate but the SDC is at least 0.81., (© 2020 American Association of Physicists in Medicine.)

Purpose: The purpose of this study was to simulate and validate organ doses from different computed tomography (CT) localizer radiograph geometries using Monte Carlo methods for a population of patients., Methods: A Monte Carlo method was developed to estimate organ doses from CT localizer radiographs using PENELOPE. The method was validated by comparing dosimetry estimates with measurements using an anthropomorphic phantom imbedded with thermoluminescent dosimeters (TLDs) scanned on a commercial CT system (Siemens SOMATOM Flash). The Monte Carlo simulation platform was then applied to conduct a population study with 57 adult computational phantoms (XCAT). In the population study, clinically relevant chest localizer protocols were simulated with the x-ray tube in anterior-posterior (AP), right lateral, and PA positions. Mean organ doses and associated standard deviations (in mGy) were then estimated for all simulations. The obtained organ doses were studied as a function of patient chest diameter. Organ doses for breast and lung were compared across different views and represented as a percentage of organ doses from rotational CT scans., Results: The validation study showed an agreement between the Monte Carlo and physical TLD measurements with a maximum percent difference of 15.5% and a mean difference of 3.5% across all organs. The XCAT population study showed that breast dose from AP localizers was the highest with a mean value of 0.24 mGy across patients, while the lung dose was relatively consistent across different localizer geometries. The organ dose estimates were found to vary across the patient population, partially explained by the changes in the patient chest diameter. The average effective dose was 0.18 mGy for AP, 0.09 mGy for lateral, and 0.08 mGy for PA localizer., Conclusion: A platform to estimate organ doses in CT localizer scans using Monte Carlo methods was implemented and validated based on comparison with physical dose measurements. The simulation platform was applied to a virtual patient population, where the localizer organ doses were found to range within 0.4%-8.6% of corresponding organ doses for a typical CT scan, 0.2%-3.3% of organ doses for a CT pulmonary angiography scan, and 1.1%-20.8% of organ doses for a low-dose lung cancer screening scan., (© 2019 American Association of Physicists in Medicine.)

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: 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.)

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.

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.

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: 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.

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.

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.

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: 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.

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.

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.

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.

Purpose: Radiation exposure from computed tomography (CT) to the public has increased the concern among radiation protection professionals. Being able to accurately assess the radiation dose patients receive during CT procedures is a crucial step in the management of CT dose. Currently, various computational anthropomorphic phantoms are used to assess radiation dose by different research groups. It is desirable to better understand how the dose results are affected by different choices of phantoms. In this study, the authors assessed the uncertainties in CT dose and risk estimation associated with different types of computational phantoms for a selected group of representative CT protocols., Methods: Routinely used CT examinations were categorized into ten body and three neurological examination categories. Organ doses, effective doses, risk indices, and conversion coefficients to effective dose and risk index (k and q factors, respectively) were estimated for these examinations for a clinical CT system (LightSpeed VCT, GE Healthcare). Four methods were used, each employing a different type of reference phantoms. The first and second methods employed a Monte Carlo program previously developed and validated in our laboratory. In the first method, the reference male and female extended cardiac-torso (XCAT) phantoms were used, which were initially created from the Visible Human data and later adjusted to match organ masses defined in ICRP publication 89. In the second method, the reference male and female phantoms described in ICRP publication 110 were used, which were initially developed from tomographic data of two patients and later modified to match ICRP 89 organ masses. The third method employed a commercial dosimetry spreadsheet (ImPACT group, London, England) with its own hermaphrodite stylized phantom. In the fourth method, another widely used dosimetry spreadsheet (CT-Expo, Medizinische Hochschule, Hannover, Germany) was employed together with its associated male and female stylized phantoms., Results: For fully irradiated organs, average coefficients of variation (COV) ranged from 0.07 to 0.22 across the four male phantoms and from 0.06 to 0.18 across the four female phantoms; for partially irradiated organs, average COV ranged from 0.13 to 0.30 across the four male phantoms and from 0.15 to 0.30 across the four female phantoms. Doses to the testes, breasts, and esophagus showed large variations between phantoms. COV for gender-averaged effective dose and k factor ranged from 0.03 to 0.23 and from 0.06 to 0.30, respectively. COV for male risk index and q factor ranged from 0.06 to 0.30 and from 0.05 to 0.36, respectively; COV for female risk index and q factor ranged from 0.06 to 0.49 and from 0.07 to 0.54, respectively., Conclusions: Despite closely matched organ mass, total body weight, and height, large differences in organ dose exist due to variation in organ location, spatial distribution, and dose approximation method. Dose differences for fully irradiated radiosensitive organs were much smaller than those for partially irradiated organs. Weighted dosimetry quantities including effective dose, male risk indices, k factors, and male q factors agreed well across phantoms. The female risk indices and q factors varied considerably across phantoms., (© 2012 American Association of Physicists in Medicine.)

Purpose: The authors previously introduced a methodology to generate a realistic three-dimensional (3D), high-resolution, computer-simulated breast phantom based on empirical data. One of the key components of such a phantom is that it provides a means to produce a realistic simulation of clinical breast compression. In the current study, they have evaluated a finite element (FE) model of compression and have demonstrated the effect of a variety of mechanical properties on the model using a dense mesh generated from empirical breast data. While several groups have demonstrated an effective compression simulation with lower density finite element meshes, the presented study offers a mesh density that is able to model the morphology of the inner breast structures more realistically than lower density meshes. This approach may prove beneficial for multimodality breast imaging research, since it provides a high level of anatomical detail throughout the simulation study., Methods: In this paper, the authors describe methods to improve the high-resolution performance of a FE compression model. In order to create the compressible breast phantom, dedicated breast CT data was segmented and a mesh was generated with 4-noded tetrahedral elements. Using an explicit FE solver to simulate breast compression, several properties were analyzed to evaluate their effect on the compression model including: mesh density, element type, density, and stiffness of various tissue types, friction between the skin and the compression plates, and breast density. Following compression, a simulated projection was generated to demonstrate the ability of the compressible breast phantom to produce realistic simulated mammographic images., Results: Small alterations in the properties of the breast model can change the final distribution of the tissue under compression by more than 1 cm; which ultimately results in different representations of the breast model in the simulated images. The model properties that impact displacement the most are mesh density, friction between the skin and the plates, and the relative stiffness of the different tissue types., Conclusions: The authors have developed a 3D, FE breast model that can yield high spatial resolution breast deformations under uniaxial compression for imaging research purposes and demonstrated that small changes in the mechanical properties can affect images generated using the phantom.

Purpose: Current methods for estimating and reporting radiation dose from CT examinations are largely patient-generic; the body size and hence dose variation from patient to patient is not reflected. Furthermore, the current protocol designs rely on dose as a surrogate for the risk of cancer incidence, neglecting the strong dependence of risk on age and gender. The purpose of this study was to develop a method for estimating patient-specific radiation dose and cancer risk from CT examinations., Methods: The study included two patients (a 5-week-old female patient and a 12-year-old male patient), who underwent 64-slice CT examinations (LightSpeed VCT, GE Healthcare) of the chest, abdomen, and pelvis at our institution in 2006. For each patient, a nonuniform rational B-spine (NURBS) based full-body computer model was created based on the patient's clinical CT data. Large organs and structures inside the image volume were individually segmented and modeled. Other organs were created by transforming an existing adult male or female full-body computer model (developed from visible human data) to match the framework defined by the segmented organs, referencing the organ volume and anthropometry data in ICRP Publication 89. A Monte Carlo program previously developed and validated for dose simulation on the LightSpeed VCT scanner was used to estimate patient-specific organ dose, from which effective dose and risks of cancer incidence were derived. Patient-specific organ dose and effective dose were compared with patient-generic CT dose quantities in current clinical use: the volume-weighted CT dose index (CTDIvol) and the effective dose derived from the dose-length product (DLP)., Results: The effective dose for the CT examination of the newborn patient (5.7 mSv) was higher but comparable to that for the CT examination of the teenager patient (4.9 mSv) due to the size-based clinical CT protocols at our institution, which employ lower scan techniques for smaller patients. However, the overall risk of cancer incidence attributable to the CT examination was much higher for the newborn (2.4 in 1000) than for the teenager (0.7 in 1000). For the two pediatric-aged patients in our study, CTDIvol underestimated dose to large organs in the scan coverage by 30%-48%. The effective dose derived from DLP using published conversion coefficients differed from that calculated using patient-specific organ dose values by -57% to 13%, when the tissue weighting factors of ICRP 60 were used, and by -63% to 28%, when the tissue weighting factors of ICRP 103 were used., Conclusions: It is possible to estimate patient-specific radiation dose and cancer risk from CT examinations by combining a validated Monte Carlo program with patient-specific anatomical models that are derived from the patients' clinical CT data and supplemented by transformed models of reference adults. With the construction of a large library of patient-specific computer models encompassing patients of all ages and weight percentiles, dose and risk can be estimated for any patient prior to or after a CT examination. Such information may aid in decisions for image utilization and can further guide the design and optimization of CT technologies and scan protocols.

Purpose: Radiation-dose awareness and optimization in CT can greatly benefit from a dose-reporting system that provides dose and risk estimates specific to each patient and each CT examination. As the first step toward patient-specific dose and risk estimation, this article aimed to develop a method for accurately assessing radiation dose from CT examinations., Methods: A Monte Carlo program was developed to model a CT system (LightSpeed VCT, GE Healthcare). The geometry of the system, the energy spectra of the x-ray source, the three-dimensional geometry of the bowtie filters, and the trajectories of source motions during axial and helical scans were explicitly modeled. To validate the accuracy of the program, a cylindrical phantom was built to enable dose measurements at seven different radial distances from its central axis. Simulated radial dose distributions in the cylindrical phantom were validated against ion chamber measurements for single axial scans at all combinations of tube potential and bowtie filter settings. The accuracy of the program was further validated using two anthropomorphic phantoms (a pediatric one-year-old phantom and an adult female phantom). Computer models of the two phantoms were created based on their CT data and were voxelized for input into the Monte Carlo program. Simulated dose at various organ locations was compared against measurements made with thermoluminescent dosimetry chips for both single axial and helical scans., Results: For the cylindrical phantom, simulations differed from measurements by -4.8% to 2.2%. For the two anthropomorphic phantoms, the discrepancies between simulations and measurements ranged between (-8.1%, 8.1%) and (-17.2%, 13.0%) for the single axial scans and the helical scans, respectively., Conclusions: The authors developed an accurate Monte Carlo program for assessing radiation dose from CT examinations. When combined with computer models of actual patients, the program can provide accurate dose estimates for specific patients.

Purpose: One issue with amplitude binning list-mode studies in SPECT for respiratory motion correction is that variation in the patient's respiratory pattern will result in binned motion states with little or no counts at various projection angles. The reduced counts result in limited-angle reconstruction artifacts which can impact the accuracy of the necessary motion estimation needed to correct the images. In this work, the authors investigate a method to overcome the effect of limited-angle reconstruction artifacts in SPECT when estimating respiratory motion., Methods: In the first pass of the reconstruction method, only the projection angles with significant counts in common between the binned respiratory states are used in order to better estimate the motion between them. After motion estimation, the estimates are used to correct for motion within iterative reconstruction using all of the acquired projection data., Results: Using simulated SPECT studies based on the NCAT phantom, the authors demonstrate the problem caused by having data available for only a limited number of angles when estimating motion and the utility of the proposed method in diminishing this error. For NCAT data sets with a clinically appropriate level of Poisson noise, the average registration error for motion with the proposed method was always less with the use of their algorithm, the reduction being statistically significant (p<0.05) in the majority of cases. The authors illustrate the ability of their method to correct the degradations caused by respiratory motion in short-axis slices and polar maps of the NCAT phantom for cases with 1 and 2 cm amplitudes of respiratory motion. In four cardiac-perfusion patients acquired on the same day, the authors demonstrate the large variability of the number of counts in the amplitude-binned projections. Finally, the authors demonstrate a visual improvement in the slices and polar maps of patient studies with the algorithm for respiratory motion correction., Conclusions: The authors' method shows promise in reducing errors in respiratory motion estimation despite the presence of limited-angle reconstruction effects due to irregularity in respiration. Improvements in image quality were observed in both simulated and clinical studies.

Purpose: The authors develop the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research., Methods: Highly detailed whole-body anatomies for the adult male and female were defined in the XCAT using nonuniform rational B-spline (NURBS) and subdivision surfaces based on segmentation of the Visible Male and Female anatomical datasets from the National Library of Medicine as well as patient datasets. Using the flexibility of these surfaces, the Visible Human anatomies were transformed to match body measurements and organ volumes for a 50th percentile (height and weight) male and female. The desired body measurements for the models were obtained using the PEOPLESIZE program that contains anthropometric dimensions categorized from 1st to the 99th percentile for US adults. The desired organ volumes were determined from ICRP Publication 89 [ICRP, "Basic anatomical and physiological data for use in radiological protection: reference values," ICRP Publication 89 (International Commission on Radiological Protection, New York, NY, 2002)]. The male and female anatomies serve as standard templates upon which anatomical variations may be modeled in the XCAT through user-defined parameters. Parametrized models for the cardiac and respiratory motions were also incorporated into the XCAT based on high-resolution cardiac- and respiratory-gated multislice CT data. To demonstrate the usefulness of the phantom, the authors show example simulation studies in PET, SPECT, and CT using publicly available simulation packages., Results: As demonstrated in the pilot studies, the 4D XCAT (which includes thousands of anatomical structures) can produce realistic imaging data when combined with accurate models of the imaging process. With the flexibility of the NURBS surface primitives, any number of different anatomies, cardiac or respiratory motions or patterns, and spatial resolutions can be simulated to perform imaging research., Conclusions: With the ability to produce realistic, predictive 3D and 4D imaging data from populations of normal and abnormal patients under various imaging parameters, the authors conclude that the XCAT provides an important tool in imaging research to evaluate and improve imaging devices and techniques. In the field of x-ray CT, the phantom may also provide the necessary foundation with which to optimize clinical CT applications in terms of image quality versus radiation dose, an area of research that is becoming more significant with the growing use of CT.

The initial process for creating a flexible three-dimensional computer-generated breast phantom based on empirical data is described. Dedicated breast computed-tomography data were processed to suppress noise and scatter artifacts in the reconstructed image set. An automated algorithm was developed to classify the breast into its primary components. A preliminary phantom defined using subdivision surfaces was generated from the segmented data. To demonstrate potential applications of the phantom, simulated mammographic image data were acquired of the phantom using a simplistic compression model and an analytic projection algorithm directly on the surface model. The simulated image was generated using a model for a polyenergetic cone-beam projection of the compressed phantom. The methods used to create the breast phantom generate resulting images that have a high level of tissue structure detail available and appear similar to actual mammograms. Fractal dimension measurements of simulated images of the phantom are comparatively similar to measurements from images of real human subjects. A realistic and geometrically defined breast phantom that can accurately simulate imaging data may have many applications in breast imaging research.

Estimating organ residence times is an essential part of patient-specific dosimetry for radioimmunotherapy (RIT). Quantitative imaging methods for RIT are often evaluated using a single physical or simulated phantom but are intended to be applied clinically where there is variability in patient anatomy, biodistribution, and biokinetics. To provide a more relevant evaluation, the authors have thus developed a population of phantoms with realistic variations in these factors and applied it to the evaluation of quantitative imaging methods both to find the best method and to demonstrate the effects of these variations. Using whole body scans and SPECT/CT images, organ shapes and time-activity curves of 111In ibritumomab tiuxetan were measured in dosimetrically important organs in seven patients undergoing a high dose therapy regimen. Based on these measurements, we created a 3D NURBS-based cardiac-torso (NCAT)-based phantom population. SPECT and planar data at realistic count levels were then simulated using previously validated Monte Carlo simulation tools. The projections from the population were used to evaluate the accuracy and variation in accuracy of residence time estimation methods that used a time series of SPECT and planar scans, Quantitative SPECT (QSPECT) reconstruction methods were used that compensated for attenuation, scatter, and the collimator-detector response. Planar images were processed with a conventional (CPlanar) method that used geometric mean attenuation and triple-energy window scatter compensation and a quantitative planar (QPlanar) processing method that used model-based compensation for image degrading effects. Residence times were estimated from activity estimates made at each of five time points. The authors also evaluated hybrid methods that used CPlanar or QPlanar time-activity curves rescaled to the activity estimated from a single QSPECT image. The methods were evaluated in terms of mean relative error and standard deviation of the relative error in the residence time estimates taken over the phantom population. The mean errors in the residence time estimates over all the organs were < 9.9% (pure QSPECT), < 13.2% (pure QPLanar), < 7.2% (hybrid QPlanar/QSPECT), < 19.2% (hybrid CPlanar/QSPECT), and 7%-159% (pure CPlanar). The standard deviations of the errors for all the organs over all the phantoms were < 9.9%, < 11.9%, < 10.8%, < 22.0%, and < 107.9% for the same methods, respectively. The processing methods differed both in terms of their average accuracy and the variation of the accuracy over the population of phantoms, thus demonstrating the importance of using a phantom population in evaluating quantitative imaging methods. Hybrid CPlanar/QSPECT provided improved accuracy compared to pure CPlanar and required the addition of only a single SPECT acquisition. The QPlanar or hybrid QPlanar/QSPECT methods had mean errors and standard deviations of errors that approached those of pure QSPECT while providing simplified image acquisition protocols, and thus may be more clinically practical.

The authors develop a unique CT simulation tool based on the 4D extended cardiac-torso (XCAT) phantom, a whole-body computer model of the human anatomy and physiology based on NURBS surfaces. Unlike current phantoms in CT based on simple mathematical primitives, the 4D XCAT provides an accurate representation of the complex human anatomy and has the advantage, due to its design, that its organ shapes can be changed to realistically model anatomical variations and patient motion. A disadvantage to the NURBS basis of the XCAT, however, is that the mathematical complexity of the surfaces makes the calculation of line integrals through the phantom difficult. They have to be calculated using iterative procedures; therefore, the calculation of CT projections is much slower than for simpler mathematical phantoms. To overcome this limitation, the authors used efficient ray tracing techniques from computer graphics, to develop a fast analytic projection algorithm to accurately calculate CT projections directly from the surface definition of the XCAT phantom given parameters defining the CT scanner and geometry. Using this tool, realistic high-resolution 3D and 4D projection images can be simulated and reconstructed from the XCAT within a reasonable amount of time. In comparison with other simulators with geometrically defined organs, the XCAT-based algorithm was found to be only three times slower in generating a projection data set of the same anatomical structures using a single 3.2 GHz processor. To overcome this decrease in speed would, therefore, only require running the projection algorithm in parallel over three processors. With the ever decreasing cost of computers and the rise of faster processors and multi-processor systems and clusters, this slowdown is basically inconsequential, especially given the vast improvement the XCAT offers in terms of realism and the ability to generate 3D and 4D data from anatomically diverse patients. As such, the authors conclude that the efficient XCAT-based CT simulator developed in this work will have applications in a broad range of CT imaging research.