A reconstruction algorithm for coherent scatter computed tomography based on filtered back-projection

Fast X-ray diffraction (XRD) tomography for enhanced identification of materials

X-ray computed tomography (CT) is a commercially established modality for imaging large objects like passenger luggage. CT can provide the density and the effective atomic number, which is not always sufficient to identify threats like explosives and narcotics, since they can have a similar composition to benign plastics, glass, or light metals. In these cases, X-ray diffraction (XRD) may be better suited to distinguish the threats. Unfortunately, the diffracted photon flux is typically much weaker than the transmitted one. Measurement of quality XRD data is therefore slower compared to CT, which is an economic challenge for potential customers like airports. In this article we numerically analyze a novel low-cost scanner design which captures CT and XRD signals simultaneously, and uses the least possible collimation to maximize the flux. To simulate a realistic instrument, we propose a forward model that includes the resolution-limiting effects of the polychromatic spectrum, the detector, and all the finite-size geometric factors. We then show how to reconstruct XRD patterns from a large phantom with multiple diffracting objects. We include a reasonable amount of photon counting noise (Poisson statistics), as well as measurement bias (incoherent scattering). Our XRD reconstruction adds material-specific information, albeit at a low resolution, to the already existing CT image, thus improving threat detection. Our theoretical model is implemented in GPU (Graphics Processing Unit) accelerated software which can be used to further optimize scanner designs for applications in security, healthcare, and manufacturing quality control.

Compressed sensing two-dimensional Bragg scatter imaging

Here we introduce a new reconstruction technique for two-dimensional Bragg scattering tomography (BST), based on the Radon transform models of Webber and Miller [Inverse Probl. Imaging15, 683 (2021).10.3934/ipi.2021010]. Our method uses a combination of ideas from multibang control and microlocal analysis to construct an objective function which can regularize the BST artifacts; specifically the boundary artifacts due to sharp cutoff in sinogram space (as observed in [arXiv preprint, arXiv:2007.00208 (2020)]), and artifacts arising from approximations made in constructing the model used for inversion. We then test our algorithm in a variety of Monte Carlo (MC) simulated examples of practical interest in airport baggage screening and threat detection. The data used in our studies is generated with a novel Monte-Carlo code presented here. The model, which is available from the authors upon request, captures both the Bragg scatter effects described by BST as well as beam attenuation and Compton scatter.

Accuracy assessment and characterization of x-ray coded aperture coherent scatter spectral imaging for breast cancer classification

Although transmission-based x-ray imaging is the most commonly used imaging approach for breast cancer detection, it exhibits false negative rates higher than 15%. To improve cancer detection accuracy, x-ray coherent scatter computed tomography (CSCT) has been explored to potentially detect cancer with greater consistency. However, the 10-min scan duration of CSCT limits its possible clinical applications. The coded aperture coherent scatter spectral imaging (CACSSI) technique has been shown to reduce scan time through enabling single-angle imaging while providing high detection accuracy. Here, we use Monte Carlo simulations to test analytical optimization studies of the CACSSI technique, specifically for detecting cancer in ex vivo breast samples. An anthropomorphic breast tissue phantom was modeled, a CACSSI imaging system was virtually simulated to image the phantom, a diagnostic voxel classification algorithm was applied to all reconstructed voxels in the phantom, and receiver-operator characteristics analysis of the voxel classification was used to evaluate and characterize the imaging system for a range of parameters that have been optimized in a prior analytical study. The results indicate that CACSSI is able to identify the distribution of cancerous and healthy tissues (i.e., fibroglandular, adipose, or a mix of the two) in tissue samples with a cancerous voxel identification area-under-the-curve of 0.94 through a scan lasting less than 10 s per slice. These results show that coded aperture scatter imaging has the potential to provide scatter images that automatically differentiate cancerous and healthy tissue within ex vivo samples. Furthermore, the results indicate potential CACSSI imaging system configurations for implementation in subsequent imaging development studies.

Direct three-dimensional coherently scattered x-ray microtomography

PURPOSE

It has been shown that coherently scattered x rays can be used to discriminate and identify specific components in a mixture of low atomic weight materials. The authors demonstrated a new method of doing coherently scattered x-ray tomography with a thin sheet of x ray.

METHODS

A collimated x-ray fan-beam, a parallel polycapillary collimator, and a phantom consisting of several biocompatible materials of low attenuation-based contrast were used to investigate the feasibility of the method. Because of the particular experimental setup, only the phantom translation perpendicular to the x-ray beam is needed and, thus, there is no need of Radon-type tomographic reconstruction, except for the correction of the attenuation to the primary and scattered x rays, which was performed by using a conventional attenuation-based tomographic image data set. The coherent scatter image contrast changes with momentum transfer among component materials in the specimen were investigated with multiple x-ray sources with narrow bandwidth spectra generated with anode and filter combinations of Cu/Ni (8 keV), Mo/Zr (18 keV), and Ag/Pd (22 keV) and at multiple scatter angles by orienting the detector and polycapillary collimator at different angles to the illuminating x ray.

RESULTS

The contrast among different materials changes with the x-ray source energy and the angle at which the image was measured. The coherent scatter profiles obtained from the coherent scatter images are consistent with the published results.

CONCLUSIONS

This method can be used to directly generate the three-dimensional coherent scatter images of small animal, biopsies, or other small objects with low atomic weight biological or similar synthetic materials with low attenuation contrast. With equipment optimized, submillimeter spatial resolution may be achieved.

High speed laser tomography system

A high speed laser tomography system was developed capable of acquiring three-dimensional (3D) images of optically thin clouds of moving micron-sized particles. It operates by parallel-shifting an illuminating laser sheet with a pair of galvanometer-driven mirrors and synchronously recording two-dimensional (2D) images of thin slices of the imaged volume. The maximum scanning speed achieved was 120,000 slices/s, sequences of 24 volume scans (up to 256 slices each) have been obtained. The 2D slices were stacked to form 3D images of the volume, then the positions of the particles were identified and followed in the consecutive scans. The system was used to image a complex plasma with particles moving at speeds up to cm/s.