Inclusion of a GaAs detector model in the Photon Counting Toolkit software for the study of breast imaging systems

We present an upgraded version of the Photon Counting Toolkit (PcTK), a freely available by request MATLAB tool for the simulation of semiconductor-based photon counting detectors (PCD), which has been extended and validated to account for gallium arsenide (GaAs)-based PCD(s). The modified PcTK version was validated by performing simulations and acquiring experimental data for three different cases. The LAMBDA 60 K module planar detector (X-Spectrum GmbH, Germany) based on the Medipix3 ASIC technology was used in all cases. This detector has a 500-μm thick GaAs sensor and a 256 × 256-pixel array with 55 μm pixel size. The first validation was a comparison between simulated and measured spectra from a 109Cd radionuclide source. In the second validation study, experimental measurements and simulations of mammography spectra were generated to observe the performance of the GaAs version of the PcTK with polychromatic radiation used in conventional x-ray imaging systems. The third validation study used single event analysis to validate the spatio-energetic model of the extended PcTK version. Overall, the software provided a good agreement between simulated and experimental data, validating the accuracy of the GaAs model. The software could be an attractive tool for accurate simulation of breast imaging modalities relying on photon counting detectors and therefore could assist in their characterization and optimization.

Classification of breast microcalcifications with GaAs photon-counting spectral mammography using an inverse problem approach

The purpose of this study was to investigate the use of a Gallium Arsenide (GaAs) photon-counting spectral mammography system to differentiate between Type I and Type II calcifications. Type I calcifications, consisting of calcium oxalate dihydrate (CO) or weddellite compounds are more often associated with benign lesions in the breast, and Type II calcifications containing hydroxyapatite (HA) are associated with both benign and malignant lesions in the breast. To be able to differentiate between these two calcification types, it is necessary to be able to estimate the full spectrum of the x-ray beam transmitted through the breast. We propose a novel method for estimating the energy-dependent x-ray transmission fraction of a beam using a photon counting detector with a limited number of energy bins. Using the estimated x-ray transmission through microcalcifications, it was observed that calcification type can be accurately estimated with machine learning. The study was carried out on a custom-built laboratory benchtop system using the SANTIS 0804 GaAs detector prototype system from DECTRIS Ltd with two energy thresholds enabled. Four energy thresholds detector was simulated by taking two separate acquisitions in which two energy thresholds were enabled for each acquisition and set at (12 keV, 21 keV) and then (29 keV, 36 keV). Measurements were performed using BR3D (CIRS, Norfolk, VA) breast imaging phantoms mimicking 100% adipose and 100% glandular tissues swirled together in an approximate 50/50 ratio by weight with the addition of in-house-developed synthetic microcalcifications. First, an inverse problem-based approach was used to estimate the full energy x-ray transmission fraction factor using known basis transmission factors from varying thicknesses of aluminum and polymethyl methacrylate (PMMA). Second, the classification of Type I and Type II calcifications was performed using the estimated energy-dependent transmission fraction factors for the pixels containing calcifications. The results were analyzed using receiver operating characteristic (ROC) analysis and demonstrated good discrimination performance with the area under the ROC curve greater than 84%. They indicated that GaAs photon-counting spectral mammography has potential use as a non-invasive method for discrimination between Type I and Type II calcifications. Results from this study suggested that GaAs-based spectral mammography could serve as a non-invasive measure for ruling out malignancy of calcifications found in the breast. Additional studies in more clinically realistic conditions involving breast tissues samples with smaller microcalcification specks should be performed to further explore the feasibility of this approach.

Theoretical comparison and optimization of cadmium telluride and gallium arsenide photon-counting detectors for contrast-enhanced spectral mammography

Purpose

Most photon-counting detectors (PCDs) being developed use cadmium telluride (CdTe), which has nonoptimal characteristic x-ray emission with energies in the range used for breast imaging. New PCD using a gallium arsenide (GaAs) has been developed. Since GaAs has characteristic x-rays lower in energy than those of CdTe, it is hypothesized that this new PCD might be beneficial for spectral x-ray breast imaging.

Approach

We performed simulations using realistic mammography x-ray spectra with both CdTe and GaAs PCDs. Five different experiments were conducted, each comparing the performance of CdTe and GaAs: (1) sensitivity of iodine quantification to charge cloud size and electronic noise, (2) effect of photon spectrum on iodine quantification, (3) effect of varying the number of energy bins, (4) a dose analysis to assess any possible dose reduction from using either detector, and (5) spectral performance of ideal CdTe and GaAs PCDs. For each study, 3 sets of 5000 noise realizations were used to calculate the Cramer-Rao lower bound (CRLB) of iodine quantification.

Results

For all spectra studied, GaAs gave a lower CRLB for iodine quantification, with 10 of the 12 spectra showing a statistically significant difference ( p ≤ 0.05 ). The photon energy spectrum that optimized iodine detection for both detector materials was the 40 kVp beam with 2-mm Al filtration, which produced CRLBs of 0.282    cm 2 and 0.257    cm 2 for CdTe and GaAs, respectively, when using five energy bins.

Conclusion

GaAs is a promising detector material for contrast-enhanced spectral mammography that offers better spectral performance than CdTe.

An empirical method for geometric calibration of a photon counting detector-based cone beam CT system

BACKGROUND

Geometric calibration is essential in developing a reliable computed tomography (CT) system. It involves estimating the geometry under which the angular projections are acquired. Geometric calibration of cone beam CTs employing small area detectors, such as currently available photon counting detectors (PCDs), is challenging when using traditional-based methods due to detectors' limited areas.

OBJECTIVE

This study presented an empirical method for the geometric calibration of small area PCD-based cone beam CT systems.

METHODS

Unlike the traditional methods, we developed an iterative optimization procedure to determine geometric parameters using the reconstructed images of small metal ball bearings (BBs) embedded in a custom-built phantom. An objective function incorporating the sphericities and symmetries of the embedded BBs was defined to assess performance of the reconstruction algorithm with the given initial estimated set of geometric parameters. The optimal parameter values were those which minimized the objective function. The TIGRE toolbox was employed for fast tomographic reconstruction. To evaluate the proposed method, computer simulations were carried out using various numbers of spheres placed in various locations. Furthermore, efficacy of the method was experimentally assessed using a custom-made benchtop PCD-based cone beam CT.

RESULTS

Computer simulations validated the accuracy and reproducibility of the proposed method. The precise estimation of the geometric parameters of the benchtop revealed high-quality imaging in CT reconstruction of a breast phantom. Within the phantom, the cylindrical holes, fibers, and speck groups were imaged in high fidelity. The CNR analysis further revealed the quantitative improvements of the reconstruction performed with the estimated parameters using the proposed method.

CONCLUSION

Apart from the computational cost, we concluded that the method was easy to implement and robust.

Fabrication of microcalcifications for insertion into phantoms used to evaluate x-ray breast imaging systems

Physical breast phantoms can be used to evaluate x-ray imaging systems such as mammography, digital breast tomosynthesis and dedicated breast computed tomography (bCT). These phantoms typically attempt to mimic x-ray attenuation properties of adipose and fibroglandular tissues within the breast. In order to use these phantoms for task-based objective assessment of image quality, relevant diagnostic features should be modeled within the phantom, such as mass lesions and/or microcalcifications. Evaluating imaging system performance in detecting microcalcifications is of particular interest due to its' clinical significance. Many previously-developed phantoms have used materials that model microcalcifications using unrealistic chemical composition, which do not accurately portray their desired x-ray attenuation and scatter properties. We report here on a new method for developing real microcalcification simulants that can be embedded in breast phantoms. This was achieved in several steps, including cross-linking hydroxyapatite and calcium oxalate powders with a binder called polyvinylpyrrolidone (PVP), and mechanical compression. The fabricated microcalcifications were evaluated by measuring their x-ray attenuation and scatter properties using x-ray spectroscopy and x-ray diffraction systems, respectively, and were demonstrated with x-ray mammography and bCT images. Results suggest that using these microcalcification models will make breast phantoms more realistic for use in evaluating task-based detection performance of the abovementioned breast imaging techniques, and bode well for extending their use to spectral imaging and x-ray coherent scatter computed tomography.

Erratum: Characterization of a GaAs photon-counting detector for mammography

The erratum corrects an error in Fig. 6 of the originally published article.

Exploring CNN potential in discriminating benign and malignant calcifications in conventional and dual-energy FFDM: simulations and experimental observations

Purpose: Deep convolutional neural networks (CNN) have demonstrated impressive success in various image classification tasks. We investigated the use of CNNs to distinguish between benign and malignant microcalcifications, using either conventional or dual-energy mammography x-ray images. The two kinds of calcifications, known as type-I (calcium oxalate crystals) and type-II (calcium phosphate aggregations), have different attenuation properties in the mammographic energy range. However, variations in microcalcification shape, size, and density as well as compressed breast thickness and breast tissue background make this a challenging discrimination task for the human visual system. Approach: Simulations (conventional and dual-energy mammography) and phantom experiments (conventional mammography only) were conducted using the range of breast thicknesses and randomly shaped microcalcifications. The off-the-shelf Resnet-18 CNN was trained on the regions of interest with calcification clusters of the two kinds. Results: Both Monte Carlo simulations and experimental phantom data suggest that deep neural networks can be trained to separate the two classes of calcifications with high accuracy, using dual-energy mammograms. Conclusions: Our work shows the encouraging results of using the CNNs for non-invasive testing for type-I and type-II microcalcifications and may stimulate further research in this area with expanding presence of the novel breast imaging modalities like dual-energy mammography or systems using photon-counting detectors.

Characterization of a GaAs photon-counting detector for mammography

Purpose: The purpose of this study was to evaluate the potential of a prototype gallium arsenide (GaAs) photon-counting detector (PCD) for imaging of the breast. Approach: First, the contrast-to-noise ratio (CNR) using different aluminum/poly(methyl methacrylate) (PMMA) phantoms of different thicknesses were measured. Second, microcalcification detection accuracy using a receiver operating characteristic study with three observers reading an ensemble of images was measured. Finally, the feasibility of using a GaAs system with two energy bins for contrast-enhanced mammography was investigated. Results: For the first two studies, the GaAs detector was compared with a commercial mammography system. The CNR was estimated by imaging 18-, 36-, and 110 - μ m -thick aluminum targets placed on top of 6 cm of PMMA plates and was found to be similar or better over a range of exposures. We observed a similar performance of detecting microcalcifications with the GaAs detector over a range of clinically applicable dose levels with a small increase at lower dose levels. The results also showed that contrast-enhanced spectral mammography using a GaAs PCD is feasible and beneficial. Conclusions: Results from this study suggest that performance with this new detector seems either slightly improved or equivalent to a commercial mammography system that used an energy-integrated detector. No attempt at optimizing exposure techniques for the GaAs detector was performed. Further research is needed to determine optimal acquisition parameters for the GaAs detector and to develop more sophisticated material decomposition algorithms that promise to provide improved quantitative estimates of iodine uptake.

Assessment of task-based performance from five clinical DBT systems using an anthropomorphic breast phantom

PURPOSE

Digital breast tomosynthesis (DBT) is a limited-angle tomographic breast imaging modality that can be used for breast cancer screening in conjunction with full-field digital mammography (FFDM) or synthetic mammography (SM). Currently, there are five commercial DBT systems that have been approved by the U.S. FDA for breast cancer screening, all varying greatly in design and imaging protocol. Because the systems are different in technical specifications, there is a need for a quantitative approach for assessing them. In this study, the DBT systems are assessed using a novel methodology with an inkjet-printed anthropomorphic phantom and four alternative forced choice (4AFC) study scheme.

METHOD

A breast phantom was fabricated using inkjet printing and parchment paper. The phantom contained 5-mm spiculated masses fabricated with potassium iodide (KI)-doped ink and microcalcifications (MCs) made with calcium hydroxyapatite. Images of the phantom were acquired on all five systems with DBT, FFDM, and SM modalities where available using beam settings under automatic exposure control. A 4AFC study was conducted to assess reader performance with each signal under each modality. Statistical analysis was performed on the data to determine proportion correct (PC), standard deviations, and levels of significance.

RESULTS

For masses, overall detection was highest with DBT. The difference in PC was statistically significant between DBT and SM for most systems. A relationship was observed between increasing PC and greater gantry span. For MCs, performance was highest with DBT and FFDM compared to SM. The difference between PC of DBT and PC of SM was statistically significant for all manufacturers.

CONCLUSIONS

This methodology represents a novel approach for evaluating systems. This study is the first of its kind to use an inkjet-printed anthropomorphic phantom with realistic signals to assess performance of clinical DBT imaging systems.

Label-free X-ray estimation of brain amyloid burden

Amyloid plaque deposits in the brain are indicative of Alzheimer’s and other diseases. Measurements of brain amyloid burden in small animals require laborious post-mortem histological analysis or resource-intensive, contrast-enhanced imaging techniques. We describe a label-free method based on spectral small-angle X-ray scattering with a polychromatic beam for in vivo estimation of brain amyloid burden. Our findings comparing 5XFAD versus wild-type mice correlate well with histology, showing promise for a fast and practical in vivo technique.

Feasibility of a label-free X-ray method to estimate brain amyloid load in small animals

BACKGROUND

Amyloid plaque in the brain is associated with a wide range of neurodegenerative diseases such as Alzheimer's and Parkinson's and defined as aggregates of amyloid fibrils rich in β-sheet structures.

NEW METHOD

We report a label-free method based on small-angle X-ray scattering (SAXS) to estimate amyloid load in an intact mouse head with skull. The method is based on recording and analyzing the X rays elastically scattered from the β-sheets of amyloid plaques in a mouse head at angles smaller than 10° and energies between 30 and 45 keV. The method is demonstrated by acquiring the spectral SAXS data of an amyloid model and an excised head from a wild-type mouse for 600 s.

RESULTS

We captured the distinct scattering peaks of the amyloid plaques at momentum transfer (q) of 6 and 13 nm^-1 associated with β-sheet structure. We first show linear correlation between the mass fraction of the amyloid target and the area under the peak (AUP) of the scattering curve. We report results for estimating amyloid load in a fixed mouse head by recovering the characteristic scattering signal from the amyloid target situated at various locations. The coefficient of variation in the amyloid load estimate is found to be less than 10%.

COMPARISON WITH EXISTING METHODS

There are no previously described label-free X-ray methods for the estimation of amyloid load in an intact head.

CONCLUSIONS

We demonstrated the feasibility of a label-free method based on SAXS to potentially estimate brain amyloid in small animals.

CT metal artifact reduction algorithms: Toward a framework for objective performance assessment

Purpose

Although several metal artifact reduction (MAR) algorithms for computed tomography (CT) scanning are commercially available, no quantitative, rigorous, and reproducible method exists for assessing their performance. The lack of assessment methods poses a challenge to regulators, consumers, and industry. We explored a phantom‐based framework for assessing an important aspect of MAR performance: how applying MAR in the presence of metal affects model observer performance at a low‐contrast detectability (LCD) task This work is, to our knowledge, the first model observer–based framework for the evaluation of MAR algorithms in the published literature.

Methods

We designed a numerical head phantom with metal implants. In order to incorporate an element of randomness, the phantom included a rotatable inset with an inhomogeneous background. We generated simulated projection data for the phantom. We applied two variants of a simple MAR algorithm, sinogram inpainting, to the projection data, that we reconstructed using filtered backprojection. To assess how MAR affected observer performance, we examined the detectability of a signal at the center of a region of interest (ROI) by a channelized Hotelling observer (CHO). As a figure of merit, we used the area under the ROC curve (AUC).

Results

We used simulation to test our framework on two variants of the MAR technique of sinogram inpainting. We found that our method was able to resolve the difference in two different MAR algorithms’ effect on LCD task performance, as well as the difference in task performances when MAR was applied, vs not.

Conclusion

We laid out a phantom‐based framework for objective assessment of how MAR impacts low‐contrast detectability, that we tested on two MAR algorithms. Our results demonstrate the importance of testing MAR performance over a range of object and imaging parameters, since applying MAR does not always improve the quality of an image for a given diagnostic task. Our framework is an initial step toward developing a more comprehensive objective assessment method for MAR, which would require developing additional phantoms and methods specific to various clinical applications of MAR, and increasing study efficiency.

Identification of amyloid plaques in the brain using an x-ray photon-counting strip detector

Brain aggregates of β amyloid (βA) protein plaques have been widely recognized as associated with many neurodegenerative diseases, and their identification can assist in the early diagnosis of Alzheimer’s disease. We investigate the feasibility of using a spectral x-ray coherent scatter system with a silicon strip photon-counting detector for identifying brain βA protein plaques. This approach is based on differences in the structure of amyloid, white and grey matter in the brain. We simulated an energy- and angular-dispersive X-ray diffraction system with an x-ray pencil beam and Silicon strip sensor, energy-resolving detectors. The polychromatic beam is geometrically focused toward a region of interest in the brain. First, the open-source MC-GPU code for Monte Carlo transport was modified to accommodate the detector model. Second, brain phantoms with and without βA were simulated to assess the method and determine the radiation dose required to obtain acceptable statistical power. For βA targets of 3, 4 and 5 mm sizes in a 15-cm brain model, the required incident exposure was about 0.44 mR from a 60 kVp tungsten spectrum and 3.5 mm of added aluminum filtration. The results suggest that the proposed x-ray coherent scatter technique enables the use of high energy x-ray spectra and therefore has the potential to be used for accurate in vivo detection and quantification of βA in the brain within acceptable radiation dose levels.

Classification of breast microcalcifications using dual-energy mammography

The potential of dual-energy mammography for microcalcification classification was investigated with simulation and phantom studies. Classification of type I/II calcifications was performed using the tissue attenuation ratio as a performance metric. The simulation and phantom studies were carried out using breast phantoms of 50% fibroglandular and 50% adipose tissue composition and thicknessess ranging from 3 to 6 cm. The phantoms included models of microcalcifications ranging in size between 200 and 900    μ m . The simulation study was carried out with fixed MGD of 1.5 mGy using various low- and high-kVp spectra, aluminum filtration thicknesses, and exposure distribution ratios to predict an optimized imaging protocol for the phantom study. Attenuation ratio values were calculated for microcalcification signals of different types at two different voltage settings. ROC analysis showed that classification performance as indicated by the area under the ROC curve was always greater than 0.95 for 1.5 mGy deposited mean glandular dose. This study provides encouraging first results in classifying malignant and benign microcalcifications based solely on dual-energy mammography images.

Reproducing two-dimensional mammograms with three-dimensional printed phantoms

Mammography is currently the standard imaging modality used to screen women for breast abnormalities, and, as a result, it is a tool of great importance for the early detection of breast cancer. Physical phantoms are commonly used as surrogates of breast tissue to evaluate some aspects of the performance of mammography systems. However, most phantoms do not reproduce the anatomic heterogeneity of real breasts. New fabrication technologies, such as three-dimensional (3-D) printing, have created the opportunity to build more complex, anatomically realistic breast phantoms that could potentially assist in the evaluation of mammography systems. The reproducibility and relative low cost of 3-D printed objects might also enable the development of collections of representative patient models that could be used to assess the effect of anatomical variability on system performance, hence making bench testing studies a step closer to clinical trials. The primary objective of this work is to present a simple, easily reproducible methodology to design and print 3-D objects that replicate the attenuation profile observed in real two-dimensional mammograms. The secondary objective is to evaluate the capabilities and limitations of the competing 3-D printing technologies and characterize the x-ray properties of the different materials they use. Printable phantoms can be created using the open-source code introduced, which processes a raw mammography image to estimate the amount of x-ray attenuation at each pixel, and outputs a triangle mesh object that encodes the observed attenuation map. The conversion from the observed pixel gray value to a column of printed material with equivalent attenuation requires certain assumptions and knowledge of multiple imaging system parameters, such as x-ray energy spectrum, source-to-object distance, compressed breast thickness, and average breast material attenuation. To validate the proposed methodology, x-ray projections of printed phantoms were acquired with a clinical mammography system. The quality of the printing process was evaluated by comparing the mammograms of the printed phantoms and the original mammograms used to create the phantoms. The structural similarity index and the root-mean-square error were used as objective metrics to compare the two images. A detailed description of the software, a characterization of the printed materials using x-ray spectroscopy, and an evaluation of the realism of the sample printed phantoms are presented.

Feasibility of estimating volumetric breast density from mammographic x-ray spectra using a cadmium telluride photon-counting detector

PURPOSE

Mammographic density of glandular breast tissue has a masking effect that can reduce lesion detection accuracy and is also a strong risk factor for breast cancer. Therefore, accurate quantitative estimation of breast density is clinically important. In this study, we investigate experimentally the feasibility of quantifying volumetric breast density with spectral mammography using a CdTe-based photon-counting detector.

METHODS

To demonstrate proof-of-principle, this study was carried out using the single pixel Amptek XR-100T-CdTe detector. The total number of x rays recorded by the detector from a single pencil-beam projection through 50%/50% of adipose/glandular mass fraction-equivalent phantoms was measured. Material decomposition assuming two, four, and eight energy bins was then applied to characterize the inspected phantom into adipose and glandular using log-likelihood estimation, taking into account the polychromatic source, the detector response function, and the energy-dependent attenuation.

RESULTS

Measurement tests were carried out for different doses, kVp settings, and different breast sizes. For dose of 1 mGy and above, the percent relative root mean square (RMS) errors of the estimated breast density was measured below 7% for all three phantom studies. It was also observed that some decrease in RMS errors was achieved using eight energy bins. For 3 and 4 cm thick phantoms, performance at 40 and 45 kVp showed similar performance. However, it was observed that 45 kVp showed better performance for a phantom thickness of 6 cm at low dose levels due to increased statistical variation at lower photon count levels with 40 kVp.

CONCLUSION

The results of the current study suggest that photon-counting spectral mammography systems using CdTe detectors have the potential to be used for accurate quantification of volumetric breast density on a pixel-to-pixel basis, with an RMS error of less than 7%.

Small-angle X-ray scattering characteristics of mouse brain: Planar imaging measurements and tomographic imaging simulations

Small-angle x-ray scattering (SAXS) imaging can differentiate tissue types based on their nanoscale molecular structure. However, characterization of the coherent scattering cross-section profile of relevant tissues is needed to optimally design SAXS imaging techniques for a variety of biomedical applications. Reported measured nervous tissue x-ray scattering cross sections under a synchrotron source have had limited agreement. We report a set of x-ray cross-section measurements obtained from planar SAXS imaging of 1 mm thick mouse brain (APP/PS1 wild-type) coronal slices using an 8 keV laboratory x-ray source. Two characteristic peaks were found at 0.96 and 1.60 nm⁻¹ attributed to myelin. The peak intensities varied by location in the slice. We found that regions of gray matter, white matter, and corpus callosum could be segmented by their increasing intensities of myelin peaks respectively. Measured small-angle x-ray scattering cross sections were then used to define brain tissue scattering properties in a GPU-accelerated Monte Carlo simulation of SAXS computed tomography (CT) using a higher monochromatic x-ray energy (20 keV) to study design trade-offs for noninvasive in vivo SAXS imaging on a small-animal head including radiation dose, signal-to-noise ratio (SNR), and the effect of skull presence on the previous two metrics. Simulation results show the estimated total dose to the mouse head for a single SAXS-CT slice was 149.4 mGy. The pixel SNR was approximately 30.8 for white matter material whether or not a skull was present. In this early-stage proof-of-principle work, we have demonstrated our brain cross-section data and simulation tools can be used to assess optimal instrument parameters for dedicated small-animal SAXS-CT prototypes.

Non-invasive classification of breast microcalcifications using x-ray coherent scatter computed tomography

We investigate the use of energy dispersive x-ray coherent scatter computed tomography (ED-CSCT) as a non-invasive diagnostic method to differentiate between type I and type II breast calcifications. This approach is sensitive to the differences of composition and internal crystal structure of different types of microcalcifications. The study is carried out by simulating a CSCT system with a scanning pencil beam, considering a polychromatic x-ray source and an energy-resolving photon counting detector. In a first step, the multidimensional angle and energy distributed CSCT data is reduced to the projection-space distributions of only a few components, corresponding to the expected target composition: adipose, glandular tissue, weddellite (calcium oxalate) for type I calcifications, and hydroxyapatite for type II calcifications. The maximum-likelihood estimation of scatter components algorithm used, operating in the projection space, takes into account the polychromatic source, the detector response function and the energy dependent attenuation. In the second step, component images are reconstructed from the corresponding estimated component projections using filtered backprojection. In a preliminary step the coherent scatter differential cross sections for hydroxyapatite and weddellite minerals were determined experimentally. The classification of type I or II calcifications is done using the relative contrasts of their components as the criterion. Simulation tests were carried out for different doses and energy resolutions for multiple realizations. The results were analyzed using relative/receiver operating characteristic methodology and show good discrimination ability at medium and higher doses. The noninvasive CSCT technique shows potential to further improve the breast diagnostic accuracy and reduce the number of breast biopsies.

Maximum-likelihood estimation of scatter components algorithm for x-ray coherent scatter computed tomography of the breast

Coherent scatter computed tomography (CSCT) is a reconstructive x-ray imaging technique that yields the spatially resolved coherent-scatter cross section of the investigated object revealing structural information of tissue under investigation. In the original CSCT proposals the reconstruction of images from coherently scattered x-rays is done at each scattering angle separately using analytic reconstruction. In this work we develop a maximum likelihood estimation of scatter components algorithm (ML-ESCA) that iteratively reconstructs images using a few material component basis functions from coherent scatter projection data. The proposed algorithm combines the measured scatter data at different angles into one reconstruction equation with only a few component images. Also, it accounts for data acquisition statistics and physics, modeling effects such as polychromatic energy spectrum and detector response function. We test the algorithm with simulated projection data obtained with a pencil beam setup using a new version of MC-GPU code, a Graphical Processing Unit version of PENELOPE Monte Carlo particle transport simulation code, that incorporates an improved model of x-ray coherent scattering using experimentally measured molecular interference functions. The results obtained for breast imaging phantoms using adipose and glandular tissue cross sections show that the new algorithm can separate imaging data into basic adipose and water components at radiation doses comparable with Breast Computed Tomography. Simulation results also show the potential for imaging microcalcifications. Overall, the component images obtained with ML-ESCA algorithm have a less noisy appearance than the images obtained with the conventional filtered back projection algorithm for each individual scattering angle. An optimization study for x-ray energy range selection for breast CSCT is also presented.

Monte Carlo X-ray transport simulation of small-angle X-ray scattering instruments using measured sample cross sections

Small-angle X-ray scattering (SAXS) has recently been proposed as a novel noninvasivein vivomolecular imaging technique to characterize molecular interactions deep within the body using high-contrast probes. This article describes a detailed Monte Carlo X-ray transport simulation technique that utilizes user-provided cross sections to describe X-ray interaction in virtual samples and explore SAXS instrument design choices. The accuracy of the simulation code is validated with sample material cross sections derived from analytical models and empirical measurements of a homogeneous spherical gold nanoparticle (GNP) monomer, a dimer and heterogeneous mixtures of the two in aqueous solution. Analytical and measured scattering profiles from these samples were converted to cross sections using an absolute water standard. Our Monte Carlo estimates of the fraction of dimers from analytically derived and empirically derived cross sections are strongly correlated, with less than 1.5 and 16% error, respectively, to the expected concentration of monomer and dimer species. In addition, a variety of monoenergetic X-ray beams were simulated to investigate coherent scatteringversusradiation dose for a range of sample sizes. For GNP spheres in aqueous solution, the energy range that produces the most coherent scattering at the detector per deposited energy was between 31 and 49 keV for a sample thickness of 1 mm to 10 cm. The method described here for the detailed simulation of SAXS using measured and modeled cross sections will enable instrumentation optimization forin vivomolecular imaging applications.