Small-angle scatter tomography with a photon-counting detector array

Method of sparse-view coded-aperture x-ray diffraction tomography

Objective.X-ray diffraction (XRD) has been considered as a valuable diagnostic technology providing material specific 'finger-print' information i.e. XRD pattern to distinguish different biological tissues. XRD tomography (XRDT) further obtains spatial-resolved XRD pattern distribution, which has become a frontier biological sample inspection method. Currently, XRD computed tomography (XRD-CT) featured by the conventional CT scan mode with rotation has the best spatial resolution among various XRDT methods, but its scan process takes hours. Meanwhile, snapshot XRDT methods such as coded-aperture XRDT (CA-XRDT) aim at direct imaging without scan movements. With compressed-sensing acquisition applied, CA-XRDT significantly shortens data acquisition time. However, the snapshot acquisition results in a significant drop in spatial resolution. Hence, we need an advanced XRDT method that significantly accelerates XRD-CT acquisition and still maintains an acceptable imaging accuracy for biological sample inspection.Approach.Inspired by the high spatial resolution of XRD-CT from rotational scan and the fast compressed-sensing acquisition in snapshot CA-XRDT (SnapshotCA-XRDT), we proposed a new XRDT imaging method: sparse-view rotational CA-XRDT (RotationCA-XRDT). It takes SnapshotCA-XRDT as a preliminary depth-resolved XRDT method, and combines rotational scan to significantly improve the spatial resolution. A model-based iterative reconstruction (MBIR) method is adopted for RotationCA-XRDT. Moreover, we suggest a refined system model calculation for the RotationCA-XRDT MBIR which is a key factor to improve reconstruction image quality.Main results.We conducted our experimental validation based on Monte-Carlo simulation for a breast sample. The results show that the proposed RotationCA-XRDT method succeeded in producing good images for detecting 2 mm square carcinoma with a 15-view scan. The spatial resolution is significantly improved from current SnapshotCA-XRDT methods. With our refined system model, MBIR can obtain high quality images with little artifacts.Significance.In this work, we proposed a new high spatial resolution XRDT method combining coded-aperture compressed-sensing acquisition and sparse-view scan. The proposed RotationCA-XRDT method obtained significantly better image resolution than current SnapshotCA-XRDT methods in the field. It is of great potential for biological sample XRDT inspection. The proposed RotationCA-XRDT is the fastest millimetre-resolution XRDT method in the field which reduces the scan time from hours to minutes.

Depth resolved pencil beam radiography using AI - a proof of principle study

AIMS

Clinical radiographic imaging is seated upon the principle of differential keV photon transmission through an object. At clinical x-ray energies the scattering of photons causes signal noise and is utilized solely for transmission measurements. However, scatter - particularly Compton scatter, is characterizable. In this work we hypothesized that modern radiation sources and detectors paired with deep learning techniques can use scattered photon information constructively to resolve superimposed attenuators in planar x-ray imaging.

METHODS

We simulated a monoenergetic x-ray imaging system consisting of a pencil beam x-ray source directed at an imaging target positioned in front of a high spatial- and energy-resolution detector array. The setup maximizes information capture of transmitted photons by measuring off-axis scatter location and energy. The signal was analyzed by a convolutional neural network, and a description of scattering material along the axis of the beam was derived. The system was virtually designed/tested using Monte Carlo processing of simple phantoms consisting of 10 pseudo-randomly stacked air/bone/water materials, and the network was trained by solving a classification problem.

RESULTS

From our simulations we were able to resolve traversed material depth information to a high degree, within our simple imaging task. The average accuracy of the material identification along the beam was 0.91±0.01, with slightly higher accuracy towards the entrance/exit peripheral surfaces of the object. The average sensitivity and specificity was 0.91 and 0.95, respectively.

CONCLUSIONS

Our work provides proof of principle that deep learning techniques can be used to analyze scattered photon patterns which can constructively contribute to the information content in radiography, here used to infer depth information in a traditional 2D planar setup. This principle, and our results, demonstrate that the information in Compton scattered photons may provide a basis for further development. The work was limited by simple testing scenarios and without yet integrating complexities or optimizations. The ability to scale performance to the clinic remains unexplored and requires further study.

Reciprocal-FDK reconstruction for x-ray diffraction computed tomography

Objective. X-ray diffraction (XRD) technology uses x-ray small-angle scattering signal for material analysis, which is highly sensitive to material inter-molecular structure. To meet the high spatial resolution requirement in applications such as medical imaging, XRD computed tomography (XRDCT) has been proposed to provide XRD intensity with improved spatial resolution from point-wise XRD scan. In XRDCT, 2D spatial tomography corresponds to a 3D reconstruction problem with the third dimension being the XRD spectrum dimension, i.e. the momentum transfer dimension. Current works in the field have studied reconstruction methods for either angular-dispersive XRDCT or energy-dispersive XRDCT for small samples. The approximations used are only suitable for regions near the XRDCT iso-center. A new XRDCT reconstruction method is needed for more general imaging applications. Approach. We derive a new FDK-type reconstruction method (Reciprocal-FDK) for XRDCT without limitation on object size. By introducing a set of reciprocal variables, the XRDCT model is transformed into a classical cone-parallel CT model, which is an extension of a circular-trajectory cone-beam CT model, after which the FDK method is applied for XRDCT reconstruction. Main results. Both analytical simulation and Monte Carlo simulation experiments are conducted to validate the XRDCT reconstruction method. The results show that when compared to existing analytical reconstruction methods, there are improvements in the proposed Reciprocal-FDK method with regard to relative structure reconstruction and XRD pattern peak reconstruction. Since cone-parallel CT does not satisfy the data completeness condition, cone-angle effect affects the reconstruction accuracy of XRDCT. The property of cone-angle effect in XRDCT is also analyzed with ablation studies. Significance. We propose a general analytical reconstruction method for XRDCT without constraint on object size. Reciprocal-FDK provides a complete derivation and theoretical support for XRDCT reconstruction by analogy to the well-studied cone-parallel CT model. In addition, the intrinsic problem with the XRDCT data model and the corresponding reconstruction error are discussed for the first time.

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.

Novel Detection Scheme for X-ray Small-Angle Scattering

X-ray imaging techniques, including x-ray radiography and computed tomography, have been in use for decades and proven effective and indispensable in diagnosis and therapy due to their fine resolution and fast acquisition speed. However, the innate disadvantage of x-ray is the poor soft tissue contrast. Small-angle scattering signals were shown to provide unique information about the abnormality of soft tissues that is complementary to the traditional attenuation image. Currently, there is no effective small-angle scattering detection system. In this paper, we propose a new "collimation" design dedicated to capture a small-angle scattering radiographic image directly, which carries critical pathological information for differentiation between normal and abnormal tissues. Our design consists of two interlaced gratings so that both the primary flux and Compton scattering photons are effectively blocked to leave the apertures mainly open to small-angle scattering photons. Theoretical analysis and Monte Carlo simulations demonstrate that small-angle scattering radiography is feasible with our proposed technology.

Development and assessment of a multi-beam continuous-phantom-motion x-ray scatter projection imaging system

X-ray image formation using scattered radiation can yield a superior contrast-to-noise ratio compared to conventional transmission x-ray imaging. A barrier to practical implementation of scatter imaging systems has been slow image acquisition. We have developed a projection imaging system which uses five monoenergetic pencil beams in combination with continuous phantom motion to achieve acquisition times that are practical for medical and security applications. The system was configured at the Canadian Light Source synchrotron and consists of a primary collimator, motorized stages for phantom translation, a flat-panel x-ray detector for measuring scattered x rays, and photodiodes for simultaneously measuring transmitted x rays. Image generation requires several corrections to raw data artifacts arising from the nature of the detector, x-ray source, and acquisition procedure. We developed a novel correction for pixel location inaccuracy arising from continuous phantom motion. A five-beam system had nearly five times faster acquisition than a single-beam system. Continuous motion acquisition was approximately 30 times faster than step-and-shoot acquisition. The total acquisition time for a 9 cm × 5 cm phantom with 8425 pixels was just over 2 min. Image quality was also assessed, in part to determine its relation to acquisition speed. The width of sharp material boundaries was found to be at a minimum equal to the pencil beam width (1.75 mm) and to have an additional width equal to the product of the phantom translation speed and the acquisition time per pixel (up to 1.0 mm in our experiments). Contrast-detail performance was independent of acquisition speed, depending only on phantom entrance x-ray fluence. Pixel signal-to-noise ratio measurements indicate that detector readout noise is important for the scatter data, even for phantom air kerma as high as 30 mGy. Images could be improved with a detector having lower readout noise and higher sensitivity. Its spatial resolution could be moderate. We confirmed that for the same range of λ^-1 sin(θ/2), where λ is the x-ray wavelength and θ is the scattering angle, scatter images acquired using different beam energies (33-70 keV) had nearly identical contrast.

Three-dimensional reciprocal space x-ray coherent scattering tomography of two-dimensional object

PURPOSE

X-ray coherent scattering tomography is a powerful tool in discriminating biological tissues and bio-compatible materials. Conventional x-ray scattering tomography framework can only resolve isotropic scattering profile under the assumption that the material is amorphous or in powder form, which is not true especially for biological samples with orientation-dependent structure. Previous tomography schemes based on x-ray coherent scattering failed to preserve the scattering pattern from samples with preferred orientations, or required elaborated data acquisition scheme, which could limit its application in practical settings. Here, we demonstrate a simple imaging modality to preserve the anisotropic scattering signal in three-dimensional reciprocal (momentum transfer) space of a two-dimensional sample layer.

METHODS

By incorporating detector movement along the direction of x-ray beam, combined with a tomographic data acquisition scheme, we match the five dimensions of the measurements with the five dimensions (three in momentum transfer domain, and two in spatial domain) of the object. We employed a collimated pencil beam of a table-top copper-anode x-ray tube, along with a panel detector to investigate the feasibility of our method.

RESULTS

We have demonstrated x-ray coherent scattering tomographic imaging at a spatial resolution ~2 mm and momentum transfer resolution 0.01 Å^-1 for the rotation-invariant scattering direction. For any arbitrary, non-rotation-invariant direction, the same spatial and momentum transfer resolution can be achieved based on the spatial information from the rotation-invariant direction. The reconstructed scattering profile of each pixel from the experiment is consistent with the x-ray diffraction profile of each material. The three-dimensional scattering pattern recovered from the measurement reveals the partially ordered molecular structure of Teflon wrap in our sample.

CONCLUSIONS

We extend the applicability of conventional x-ray coherent scattering tomography to the reconstruction of two-dimensional samples with anisotropic scattering profile by introducing additional degree of freedom on the detector. The presented method has the potential to achieve low-cost, high-specificity material discrimination based on x-ray coherent scattering.