Resolution analysis of a volumetric coded aperture X-ray diffraction imaging system.

BACKGROUND: A coded aperture X-ray diffraction (XRD) imaging system can measure the X-ray diffraction form factor from an object in three dimensions -X, Y and Z (depth), broadening the potential application of this technology. However, to optimize XRD systems for specific applications, it is critical to understand how to predict and quantify system performance for each use case. OBJECTIVE: The purpose of this work is to present and validate 3D spatial resolution models for XRD imaging systems with a detector-side coded aperture. METHODS: A fan beam coded aperture XRD system was used to scan 3D printed resolution phantoms placed at various locations throughout the system's field of view. The multiplexed scatter data were reconstructed using a model-based iterative reconstruction algorithm, and the resulting volumetric images were evaluated using multiple resolution criteria to compare against the known phantom resolution. We considered the full width at half max and Sparrow criterion as measures of the resolution and compared our results against analytical resolution models from the literature as well as a new theory for predicting the system resolution based on geometric arguments. RESULTS: We show that our experimental measurements are bounded by the multitude of theoretical resolution predictions, which accurately predict the observed trends and order of magnitude of the spatial and form factor resolutions. However, we find that the expected and observed resolution can vary by approximately a factor of two depending on the choice of metric and model considered. We observe depth resolutions of 7-16 mm and transverse resolutions of 0.6-2 mm for objects throughout the field of view. Furthermore, we observe tradeoffs between the spatial resolution and XRD form factor resolution as a function of sample location. CONCLUSION: The theories evaluated in this study provide a useful framework for estimating the 3D spatial resolution of a detector side coded aperture XRD imaging system. The assumptions and simplifications required by these theories can impact the overall accuracy of describing a particular system, but they also can add to the generalizability of their predictions. Furthermore, understanding the implications of the assumptions behind each theory can help predict performance, as shown by our data's placement between the conservative and idealized theories, and better guide future systems for optimized designs.

Evaluating county-level lung cancer incidence from environmental radiation exposure, PM2.5, and other exposures with regression and machine learning models.

Characterizing the interplay between exposures shaping the human exposome is vital for uncovering the etiology of complex diseases. For example, cancer risk is modified by a range of multifactorial external environmental exposures. Environmental, socioeconomic, and lifestyle factors all shape lung cancer risk. However, epidemiological studies of radon aimed at identifying populations at high risk for lung cancer often fail to consider multiple exposures simultaneously. For example, moderating factors, such as PM2.5, may affect the transport of radon progeny to lung tissue. This ecological analysis leveraged a population-level dataset from the National Cancer Institute's Surveillance, Epidemiology, and End-Results data (2013-17) to simultaneously investigate the effect of multiple sources of low-dose radiation (gross [Formula: see text] activity and indoor radon) and PM2.5 on lung cancer incidence rates in the USA. County-level factors (environmental, sociodemographic, lifestyle) were controlled for, and Poisson regression and random forest models were used to assess the association between radon exposure and lung and bronchus cancer incidence rates. Tree-based machine learning (ML) method perform better than traditional regression: Poisson regression: 6.29/7.13 (mean absolute percentage error, MAPE), 12.70/12.77 (root mean square error, RMSE); Poisson random forest regression: 1.22/1.16 (MAPE), 8.01/8.15 (RMSE). The effect of PM2.5 increased with the concentration of environmental radon, thereby confirming findings from previous studies that investigated the possible synergistic effect of radon and PM2.5 on health outcomes. In summary, the results demonstrated (1) a need to consider multiple environmental exposures when assessing radon exposure's association with lung cancer risk, thereby highlighting (1) the importance of an exposomics framework and (2) that employing ML models may capture the complex interplay between environmental exposures and health, as in the case of indoor radon exposure and lung cancer incidence.

Medical Physics Leadership Academy Journal Club (Leadership Club) program: Two-year review in building a community of leaders.

The American Association of Physicists in Medicine began the Medical Physics Leadership Academy Journal Club in the fall of 2020. The initiative was launched to provide a forum for medical physicists to learn about leadership topics using published material, discuss and reflect on the material, and consider incorporating the discussed skills into their professional practice. This report presents the framework for the MPLA Journal Club program, describes the lessons learned over the last 2 years, summarizes the data collected from attendees, and highlights the roadmap for the program moving forward.

Exploring the Spatial Patterning of Sociodemographic Disparities in Extreme Heat Exposure at Multiple Scales Across the Conterminous United States.

Climate change has led to an increase in heat-related morbidity and mortality. The impact of heat on health is unequally distributed amongst different socioeconomic and demographic groups. We use high-resolution daily air temperature-based heat wave intensity (HWI) and neighborhood-scale sociodemographic information from the conterminous United States to evaluate the spatial patterning of extreme heat exposure disparities. Assuming differences in spatial patterns at national, regional, and local scales; we assess disparities in heat exposure across race, housing characteristics, and poverty level. Our findings indicate small differences in HWI based on these factors at the national level, with the magnitude and direction of the differences varying by region. The starkest differences are present over the Northeast and Midwest, where primarily Black neighborhoods are exposed to higher HWI than predominantly White areas. At the local level, we find the largest difference by socioeconomic status. We also find that residents of nontraditional housing are more vulnerable to heat exposure. Previous studies have either evaluated such disparities for specific cities and/or used a satellite-based land surface temperature, which, although correlated with air temperature, does not provide the true measure of heat exposure. This study is the first of its kind to incorporate high-resolution gridded air temperature-based heat exposure in the evaluation of sociodemographic disparities at a national scale. The analysis suggests the unequal distribution of heat wave intensities across communities-with higher heat exposures characterizing areas with high proportions of minorities, low socioeconomic status, and homes in need of retrofitting to combat climate change.

Dose coefficients for organ dosimetry in tomosynthesis imaging of adults and pediatrics across diverse protocols.

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 (hn ) were estimated by normalizing organ dose estimates to air kerma at 70 cm (X70cm ) 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 (Kn ) 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 (hn ) containing 294 distinct combinations of patients and exam protocols was developed and made publicly available. The values of Kn 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 = Kn /Kn ,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.

A GPU-accelerated framework for individualized estimation of organ doses in digital tomosynthesis.

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.

Application of machine learning classifiers to X-ray diffraction imaging with medically relevant phantoms.

PURPOSE: Recent studies have demonstrated the ability to rapidly produce large field of view X-ray diffraction (XRD) images, which provide rich new data relevant to the understanding and analysis of disease. However, work has only just begun on developing algorithms that maximize the performance toward decision-making and diagnostic tasks. In this study, we present the implementation of and comparison between rules-based and machine learning (ML) classifiers on XRD images of medically relevant phantoms to explore the potential for increased classification performance. METHODS: Medically relevant phantoms were utilized to provide well-characterized ground-truths for comparing classifier performance. Water and polylactic acid (PLA) plastic were used as surrogates for cancerous and healthy tissue, respectively, and phantoms were created with varying levels of spatial complexity and biologically relevant features for quantitative testing of classifier performance. Our previously developed X-ray scanner was used to acquire co-registered X-ray transmission and diffraction images of the phantoms. For classification algorithms, we explored and compared two rules-based classifiers (cross-correlation, or matched-filter, and linear least-squares unmixing) and two ML classifiers (support vector machines and shallow neural networks). Reference XRD spectra (measured by a commercial diffractometer) were provided to the rules-based algorithms, while 60% of the measured XRD pixels were used for training of the ML algorithms. The area under the receiver operating characteristic curve (AUC) was used as a comparative metric between the classification algorithms, along with the accuracy performance at the midpoint threshold for each classifier. RESULTS: The AUC values for material classification were 0.994 (cross-correlation [CC]), 0.994 (least-squares [LS]), 0.995 (support vector machine [SVM]), and 0.999 (shallow neural network [SNN]). Setting the classification threshold to the midpoint for each classifier resulted in accuracy values of CC = 96.48%, LS = 96.48%, SVM = 97.36%, and SNN = 98.94%. If only considering pixels ±3 mm from water-PLA boundaries (where partial volume effects could occur due to imaging resolution limits), the classification accuracies were CC = 89.32%, LS = 89.32%, SVM = 92.03%, and SNN = 96.79%, demonstrating an even larger improvement produced by the machine-learned algorithms in spatial regions critical for imaging tasks. Classification by transmission data alone produced an AUC of 0.773 and accuracy of 85.45%, well below the performance levels of any of the classifiers applied to XRD image data. CONCLUSIONS: We demonstrated that ML-based classifiers outperformed rules-based approaches in terms of overall classification accuracy and improved the spatially resolved classification performance on XRD images of medical phantoms. In particular, the ML algorithms demonstrated considerably improved performance whenever multiple materials existed in a single voxel. The quantitative performance gains demonstrate an avenue to extract and harness XRD imaging data to improve material analysis for research, industrial, and clinical applications.

X-ray fan beam coded aperture transmission and diffraction imaging for fast material analysis.

X-ray transmission imaging has been used in a variety of applications for high-resolution measurements based on shape and density. Similarly, X-ray diffraction (XRD) imaging has been used widely for molecular structure-based identification of materials. Combining these X-ray methods has the potential to provide high-resolution material identification, exceeding the capabilities of either modality alone. However, XRD imaging methods have been limited in application by their long measurement times and poor spatial resolution, which has generally precluded combined, rapid measurements of X-ray transmission and diffraction. In this work, we present a novel X-ray fan beam coded aperture transmission and diffraction imaging system, developed using commercially available components, for rapid and accurate non-destructive imaging of industrial and biomedical specimens. The imaging system uses a 160 kV Bremsstrahlung X-ray source while achieving a spatial resolution of ≈ 1 × 1 mm2 and a spectral accuracy of > 95% with only 15 s exposures per 150 mm fan beam slice. Applications of this technology are reported in geological imaging, pharmaceutical inspection, and medical diagnosis. The performance of the imaging system indicates improved material differentiation relative to transmission imaging alone at scan times suitable for a variety of industrial and biomedical applications.

A GPU-accelerated framework for rapid estimation of scanner-specific scatter in CT for virtual imaging trials.

Virtual imaging trials (VITs), defined as the process of conducting clinical imaging trials using computer simulations, offer a time- and cost-effective alternative to traditional imaging trials for CT. The clinical potential of VITs hinges on the realism of simulations modeling the image acquisition process, where the accurate scanner-specific simulation of scatter in a time-feasible manner poses a particular challenge. To meet this need, this study proposes, develops, and validates a rapid scatter estimation framework, based on GPU-accelerated Monte Carlo (MC) simulations and denoising methods, for estimating scatter in single source, dual-source, and photon-counting CT. A CT simulator was developed to incorporate parametric models for an anti-scatter grid and a curved energy integrating detector with an energy-dependent response. The scatter estimates from the simulator were validated using physical measurements acquired on a clinical CT system using the standard single-blocker method. The MC simulator was further extended to incorporate a pre-validated model for a PCD and an additional source-detector pair to model cross scatter in dual-source configurations. To estimate scatter with desirable levels of statistical noise using a manageable computational load, two denoising methods using a (1) convolutional neural network and an (2) optimized Gaussian filter were further deployed. The viability of this framework for clinical VITs was assessed by integrating it with a scanner-specific ray-tracer program to simulate images for an image quality (Mercury) and an anthropomorphic phantom (XCAT). The simulated scatter-to-primary ratios agreed with physical measurements within 4.4% ± 10.8% across all projection angles and kVs. The differences of ∼121 HU between images with and without scatter, signifying the importance of scatter for simulating clinical images. The denoising methods preserved the magnitudes and trends observed in the reference scatter distributions, with an averaged rRMSE value of 0.91 and 0.97 for the two methods, respectively. The execution time of ∼30 s for simulating scatter in a single projection with a desirable level of statistical noise indicates a major improvement in performance, making our tool an eligible candidate for conducting extensive VITs spanning multiple patients and scan protocols.

Comparison of 12 surrogates to characterize CT radiation risk across a clinical population.

OBJECTIVES: Quantifying radiation burden is essential for justification, optimization, and personalization of CT procedures and can be characterized by a variety of risk surrogates inducing different radiological risk reflections. This study compared how twelve such metrics can characterize risk across patient populations. METHODS: This study included 1394 CT examinations (abdominopelvic and chest). Organ doses were calculated using Monte Carlo methods. The following risk surrogates were considered: volume computed tomography dose index (CTDIvol), dose-length product (DLP), size-specific dose estimate (SSDE), DLP-based effective dose (EDk ), dose to a defining organ (ODD), effective dose and risk index based on organ doses (EDOD, RI), and risk index for a 20-year-old patient (RIrp). The last three metrics were also calculated for a reference ICRP-110 model (ODD,0, ED0, and RI0). Lastly, motivated by the ICRP, an adjusted-effective dose was calculated as [Formula: see text]. A linear regression was applied to assess each metric's dependency on RI. The results were characterized in terms of risk sensitivity index (RSI) and risk differentiability index (RDI). RESULTS: The analysis reported significant differences between the metrics with EDr showing the best concordance with RI in terms of RSI and RDI. Across all metrics and protocols, RSI ranged between 0.37 (SSDE) and 1.29 (RI0); RDI ranged between 0.39 (EDk) and 0.01 (EDr) cancers × 103patients × 100 mGy. CONCLUSION: Different risk surrogates lead to different population risk characterizations. EDr exhibited a close characterization of population risk, also showing the best differentiability. Care should be exercised in drawing risk predictions from unrepresentative risk metrics applied to a population. KEY POINTS: • Radiation risk characterization in CT populations is strongly affected by the surrogate used to describe it. • Different risk surrogates can lead to different characterization of population risk. • Healthcare professionals should exercise care in ascribing an implicit risk to factors that do not closely reflect risk.

Simulation based evaluation of a fan beam coded aperture x-ray diffraction imaging system for biospecimen analysis.

X-ray diffraction (XRD) imaging yields spatially resolved, material-specific information, which can aid medical diagnosis and inform treatment. In this work we used simulations to analyze the utility of fan beam coded aperture XRD imaging for fast, high-resolution scatter imaging of biospecimens for tissue assessment. To evaluate the proposed system's utility in a specific task, we employed a deterministic model to produce simulated data from biologically realistic breast tissue phantoms and model-based reconstruction to recover a spatial map of the XRD signatures throughout the phantoms. We found an XRD spatial resolution of ≈1 mm with a mean reconstructed spectral accuracy of 0.98 ± 0.01 for a simulated 1 × 150 mm2 fan beam operating at 160 kVp, 10 mA, and 4.5 s exposures. A classifier for cancer detection was developed utilizing cross-correlation of XRD spectra against a spectral library, with a receiver operating characteristic curve with an area under the curve value of 0.972. Our results indicated a potential diagnostic modality that could aid in tasks ranging from analysis of ex-vivo pathology biospecimens to intraoperative cancer margin assessment, motivating future work to develop an experimental system while enabling the development of improved algorithms for imaging and tissue analysis-based classification performance.

Patient-Informed Organ Dose Estimation in Clinical CT: Implementation and Effective Dose Assessment in 1048 Clinical Patients.

OBJECTIVE. The purpose of this study is to comprehensively implement a patient-informed organ dose monitoring framework for clinical CT and compare the effective dose (ED) according to the patient-informed organ dose with ED according to the dose-length product (DLP) in 1048 patients. MATERIALS AND METHODS. Organ doses for a given examination are computed by matching the topogram to a computational phantom from a library of anthropomorphic phantoms and scaling the fixed tube current dose coefficients by the examination volume CT dose index (CTDIvol) and the tube-current modulation using a previously validated convolution-based technique. In this study, the library was expanded to 58 adult, 56 pediatric, five pregnant, and 12 International Commission on Radiological Protection (ICRP) reference models, and the technique was extended to include multiple protocols, a bias correction, and uncertainty estimates. The method was implemented in a clinical monitoring system to estimate organ dose and organ dose-based ED for 647 abdomen-pelvis and 401 chest examinations, which were compared with DLP-based ED using a t test. RESULTS. For the majority of the organs, the maximum errors in organ dose estimation were 18% and 8%, averaged across all protocols, without and with bias correction, respectively. For the patient examinations, DLP-based ED was significantly different from organ dose-based ED by as much as 190.9% and 234.7% for chest and abdomen-pelvis scans, respectively (mean, 9.0% and 24.3%). The differences were statistically significant (p < .001) and exhibited overestimation for larger-sized patients and underestimation for smaller-sized patients. CONCLUSION. A patient-informed organ dose estimation framework was comprehensively implemented applicable to clinical imaging of adult, pediatric, and pregnant patients. Compared with organ dose-based ED, DLP-based ED may overestimate effective dose for larger-sized patients and underestimate it for smaller-sized patients.

Virtual clinical trial for quantifying the effects of beam collimation and pitch on image quality in computed tomography.

Purpose: To utilize a virtual clinical trial (VCT) construct to investigate the effects of beam collimation and pitch on image quality (IQ) in computed tomography (CT) under different respiratory and cardiac motion rates. Approach: A computational human model [extended cardiac-torso (XCAT) phantom] with added lung lesions was used to simulate seven different rates of cardiac and respiratory motions. A validated CT simulator (DukeSim) was used in this study. A supplemental validation was done to ensure the accuracy of DukeSim across different pitches and beam collimations. Each XCAT phantom was imaged using the CT simulator at multiple pitches (0.5 to 1.5) and beam collimations (19.2 to 57.6 mm) at a constant dose level. The images were compared against the ground truth using three task-generic IQ metrics in the lungs. Additionally, the bias and variability in radiomics (morphological) feature measurements were quantified for task-specific lung lesion quantification across the studied imaging conditions. Results: All task-generic metrics degraded by 1.6% to 13.3% with increasing pitch. When imaged with motion, increasing pitch reduced motion artifacts. The IQ slightly degraded (1.3%) with changes in the studied beam collimations. Patient motion exhibited negative effects (within 7%) on the IQ. Among all features across all imaging conditions studies, compactness2 and elongation showed the largest ( - 26.5 % , 7.8%) and smallest ( - 0.8 % , 2.7%) relative bias and variability. The radiomics results were robust across the motion profiles studied. Conclusions: While high pitch and large beam collimations can negatively affect the quality of CT images, they are desirable for fast imaging. Further, our results showed no major adverse effects in morphology quantification of lung lesions with the increase in pitch or beam collimation. VCTs, such as the one demonstrated in this study, represent a viable methodology for experiments in CT.

A real-time Monte Carlo tool for individualized dose estimations in clinical CT.

The increasing awareness of the adverse effects associated with radiation exposure in computed tomography (CT) has necessesitated the quantification of dose delivered to patients for better risk assessment in the clinic. The current methods for dose quantification used in the clinic are approximations, lacking realistic models for the irradiation conditions utilized in the scan and the anatomy of the patient being imaged, which limits their relevance for a particular patient. The established gold-standard technique for individualized dose quantification uses Monte Carlo (MC) simulations to obtain patient-specific estimates of organ dose in anatomically realistic computational phantoms to provide patient-specific estimates of organ dose. Although accurate, MC simulations are computationally expensive, which limits their utility for time-constrained applications in the clinic. To overcome these shortcomings, a real-time GPU-based MC tool based on FDA's MC-GPU framework was developed for patient and scanner-specific dosimetry in the clinic. The tool was validated against (1) AAPM's TG-195 reference datasets and (2) physical measurements of dose acquired using TLD chips in adult and pediatric anthropomorphic phantoms. To demonstrate its utility towards providing individualized dose estimates, it was integrated with an automatic segmentation method for generating patient-specific models, which were then used to estimate patient- and scanner-specific organ doses for a select population of 50 adult patients using a clinically relevant CT protocol. The organ dose estimates were compared to corresponding dose estimates from a previously validated MC method based on Penelope. The dose estimates from our MC tool agreed within 5% for all organs (except thyroid) tabulated by TG-195 and within 10% for all TLD locations in the adult and pediactric phantoms, across all validation cases. Compared against Penelope, the organ dose estimates agreed within 3% on average for all organs in the patient population study. The average run duration for each patient was estimated at 23.79 s, representing a significant speedup (~700×) over existing non-parallelized MC methods. The accuracy of dose estimates combined with a significant improvement in execution times suggests a feasible solution utilizing the proposed MC tool for real-time individualized dosimetry in the clinic.

Modeling "Textured" Bones in Virtual Human Phantoms.

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

Organ doses from CT localizer radiographs: Development, validation, and application of a Monte Carlo estimation technique.

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.

Development of a scanner-specific simulation framework for photon-counting computed tomography.

The aim of this study was to develop and validate a simulation platform that generates photon-counting CT images of voxelized phantoms with detailed modeling of manufacturer-specific components including the geometry and physics of the x-ray source, source filtrations, anti-scatter grids, and photon-counting detectors. The simulator generates projection images accounting for both primary and scattered photons using a computational phantom, scanner configuration, and imaging settings. Beam hardening artifacts are corrected using a spectrum and threshold dependent water correction algorithm. Physical and computational versions of a clinical phantom (ACR) were used for validation purposes. The physical phantom was imaged using a research prototype photon-counting CT (Siemens Healthcare) with standard (macro) mode, at four dose levels and with two energy thresholds. The computational phantom was imaged with the developed simulator with the same parameters and settings used in the actual acquisition. Images from both the real and simulated acquisitions were reconstructed using a reconstruction software (FreeCT). Primary image quality metrics such as noise magnitude, noise ratio, noise correlation coefficients, noise power spectrum, CT number, in-plane modulation transfer function, and slice sensitivity profiles were extracted from both real and simulated data and compared. The simulator was further evaluated for imaging contrast materials (bismuth, iodine, and gadolinium) at three concentration levels and six energy thresholds. Qualitatively, the simulated images showed similar appearance to the real ones. Quantitatively, the average relative error in image quality measurements were all less than 4% across all the measurements. The developed simulator will enable systematic optimization and evaluation of the emerging photon-counting computed tomography technology.

DukeSim: A Realistic, Rapid, and Scanner-Specific Simulation Framework in Computed Tomography.

The purpose of this study was to develop a CT simulation platform that is: 1) compatible with voxel-based computational phantoms; 2) capable of modeling the geometry and physics of commercial CT scanners; and 3) computationally efficient. Such a simulation platform is designed to enable the virtual evaluation and optimization of CT protocols and parameters for achieving a targeted image quality while reducing radiation dose. Given a voxelized computational phantom and a parameter file describing the desired scanner and protocol, the developed platform DukeSim calculates projection images using a combination of ray-tracing and Monte Carlo techniques. DukeSim includes detailed models for the detector quantum efficiency, quantum and electronic noise, detector crosstalk, subsampling of the detector and focal spot areas, focal spot wobbling, and the bowtie filter. DukeSim was accelerated using GPU computing. The platform was validated using physical and computational versions of a phantom (Mercury phantom). Clinical and simulated CT scans of the phantom were acquired at multiple dose levels using a commercial CT scanner (Somatom Definition Flash; Siemens Healthcare). The real and simulated images were compared in terms of image contrast, noise magnitude, noise texture, and spatial resolution. The relative error between the clinical and simulated images was less than 1.4%, 0.5%, 2.6%, and 3%, for image contrast, noise magnitude, noise texture, and spatial resolution, respectively, demonstrating the high realism of DukeSim. The runtime, dependent on the imaging task and the hardware, was approximately 2-3 minutes per rotation in our study using a computer with 4 GPUs. DukeSim, when combined with realistic human phantoms, provides the necessary toolset with which to perform large-scale and realistic virtual clinical trials in a patient and scanner-specific manner.

X-ray diffraction tomography with limited projection information.

X-ray diffraction tomography (XDT) records the spatially-resolved X-ray diffraction profile of an extended object. Compared to conventional transmission-based tomography, XDT displays high intrinsic contrast among materials of similar electron density and improves the accuracy in material identification thanks to the molecular structural information carried by diffracted photons. However, due to the weak diffraction signal, a tomographic scan covering the entire object typically requires a synchrotron facility to make the acquisition time more manageable. Imaging applications in medical and industrial settings usually do not require the examination of the entire object. Therefore, a diffraction tomography modality covering only the region of interest (ROI) and subsequent image reconstruction techniques with truncated projections are highly desirable. Here we propose a table-top diffraction tomography system that can resolve the spatially-variant diffraction form factor from internal regions within extended samples. We demonstrate that the interior reconstruction maintains the material contrast while reducing the imaging time by 6 folds. The presented method could accelerate the acquisition of XDT and be applied in portable imaging applications with a reduced radiation dose.