Interpretation and Utility of the Moments of Small-Angle X-Ray Scattering Distributions

On the equivalence of the X-ray scattering retrieval with beam tracking and analyser-based imaging using a synchrotron source

X-ray phase contrast imaging (XPCI) methods give access to contrast mechanisms that are based on the refractive properties of matter on top of the absorption coefficient in conventional x-ray imaging. Ultra small angle x-ray scattering (USAXS) is a phase contrast mechanism that arises due to multiple refraction events caused by physical features of a scale below the physical resolution of the used imaging system. USAXS contrast can therefore give insight into subresolution structural information, which is an ongoing research topic in the vast field of different XPCI techniques. In this study, we quantitatively compare the USAXS signal retrieved by the beam tracking XPCI technique with the gold standard of the analyzer based imaging XPCI technique using a synchrotron x-ray source. We find that, provided certain conditions are met, the two methods measure the same quantity.

On the quantification of sample microstructure using single-exposure x-ray dark-field imaging via a single-grid setup

The size of the smallest detectable sample feature in an x-ray imaging system is usually restricted by the spatial resolution of the system. This limitation can now be overcome using the diffusive dark-field signal, which is generated by unresolved phase effects or the ultra-small-angle x-ray scattering from unresolved sample microstructures. A quantitative measure of this dark-field signal can be useful in revealing the microstructure size or material for medical diagnosis, security screening and materials science. Recently, we derived a new method to quantify the diffusive dark-field signal in terms of a scattering angle using a single-exposure grid-based approach. In this manuscript, we look at the problem of quantifying the sample microstructure size from this single-exposure dark-field signal. We do this by quantifying the diffusive dark-field signal produced by 5 different sizes of polystyrene microspheres, ranging from 1.0 to 10.8 µm, to investigate how the strength of the extracted dark-field signal changes with the sample microstructure size, [Formula: see text]. We also explore the feasibility of performing single-exposure dark-field imaging with a simple equation for the optimal propagation distance, given microstructure with a specific size and thickness, and show consistency between this model and experimental data. Our theoretical model predicts that the dark-field scattering angle is inversely proportional to [Formula: see text], which is also consistent with our experimental data.

Differential phase contrast for quantitative imaging and spectro-microscopy at a nanoprobe beamline

The interaction of a focused X-ray beam with a sample in a scanning probe experiment can provide a variety of information about the interaction volume. In many scanning probe experiments X-ray fluorescence (XRF) is supplemented with measurements of the transmitted or scattered intensity using a pixelated detector. The automated extraction of different signals from an area pixelated detector is described, in particular the methodology for extracting differential phase contrast (DPC) is demonstrated and different processing methods are compared across a range of samples. The phase shift of the transmitted X-ray beam by the sample, extracted from DPC, is also compared with ptychography measurements to provide a qualitative and quantitative comparison. While ptychography produces a superior image, DPC can offer a simple, flexible method for phase contrast imaging which can provide fast results and feedback during an experiment; furthermore, for many science problems, such as registration of XRF in a lighter matrix, DPC can provide sufficient information to meet the experimental aims. As the DPC technique is a quantitative measurement, it can be expanded to spectroscopic studies and a demonstration of DPC for spectro-microscopy measurements is presented. Where ptychography can separate the absorption and phase shifts by the sample, quantitative interpretation of a DPC image or spectro-microscopy signal can only be performed directly when absorption is negligible or where the absorption contribution is known and the contributions can be fitted.

Ptychographic reconstruction with object initialization

X-ray ptychography is a cutting edge imaging technique providing ultra-high spatial resolutions. In ptychography, phase retrieval, i.e., the recovery of a complex valued signal from intensity-only measurements, is enabled by exploiting a redundancy of information contained in diffraction patterns measured with overlapping illuminations. For samples that are considerably larger than the probe we show that during the iteration the bulk information has to propagate from the sample edges to the center. This constitutes an inherent limitation of reconstruction speed for algorithms that use a flat initialization. Here, we experimentally demonstrate that a considerable improvement of computational speed can be achieved by utilizing a low resolution sample wavefront retrieved from measured diffraction patterns as object initialization. In addition, we show that this approach avoids phase artifacts associated with large phase gradients and may alleviate the requirements on phase structure within the probe. Object initialization is computationally fast, potentially beneficial for bulky sample and compatible with flat samples. Therefore, the presented approach is readily adaptable with established ptychographic reconstruction algorithms implying a wide spread use.

Quantifying the x-ray dark-field signal in single-grid imaging

X-ray dark-field imaging reveals the sample microstructure that is unresolved when using conventional methods of x-ray imaging. In this paper, we derive a new method to extract and quantify the x-ray dark-field signal collected using a single-grid imaging set-up, and relate the signal strength to the number of sample microstructures, N. This was achieved by modelling sample-induced changes to the shadow of the upstream grid, and fitting experimental data to this model. Our results suggested that the dark-field scattering angle from our spherical microstructures deviates slightly from the theoretical model of N, which was consistent with results from other experimental methods. We believe the approach outlined here can equip quantitative dark-field imaging of small samples, particularly in cases where only one sample exposure is possible, either due to sample movement or radiation dose limitations. Future directions include an extension into directional dark-field imaging.

Fluence adaptation for contrast-based dose optimization in x-ray phase-contrast imaging

PURPOSE

X-ray phase-contrast imaging (XPCI) can provide multiple contrasts with great potentials for clinical and industrial applications, including conventional attenuation, phase contrast, and dark field. Grating-based imaging (GBI) and edge-illumination (EI) are two promising types of XPCI as the conventional x-ray sources can be directly utilized. For the GBI and EI systems, the phase-stepping acquisition with multiple exposures at a constant fluence is usually adopted in the literature.This work, however, attempts to challenge such a constant fluence concept during the phase-stepping process and proposes a fluence adaptation mechanism for dose reduction.

METHOD

Given the importance of patient radiation dose for clinical applications, numerous studies have tried to reduce patient dose in XPCI by altering imaging system designs, data acquisition, and information retrieval. Recently, analytic multiorder moment analysis has been proposed to improve the computing efficiency. In these algorithms, multiple contrasts can be calculated by summing together the weighted phase-stepping curves (PSCs) with some kernel functions, which suggests us that the raw data at different steps have different contributions for the noise in retrieved contrasts. Therefore, it is possible to improve the noise performance by adjusting the fluence distribution during the phase-stepping process directly. Based on analytic retrieval formulas and the Gaussian noise model for detected signals, we derived an optimal adaptive fluence distribution, which is proportional to the absolute weighting kernel functions and the root of original sample PSCs acquired under the constant fluence. Considering that the original sample PSC might be unavailable, we proposed two practical forms for the GBI and EI systems, which are also able to reduce the contrast noise when comparing with the constant fluence distribution. Since the kernel functions are target contrast-dependent, our proposed fluence adaptation mechanism provides a way of realizing a contrast-based dose optimization while keeping the same noise level.

RESULTS

To validate our analyses, simulations and experiments are conducted for the GBI and EI systems. Simulated results demonstrate that the dose reduction ratio between our proposed fluence distributions and the typical constant one can be about 20% for the phase contrast, which is consistent with our theoretical predictions. Although the experimental noise reduction ratios are a little smaller than the theoretical ones, low-dose experiments observe better noise performance by our proposed method. Our simulated results also give out the effective ranges of the parameters of the PSCs, such as the visibility in the GBI, the standard deviation, and the mean value in the EI, providing a guidance for the use of our proposed approach in practice.

CONCLUSIONS

In this paper, we propose a fluence adaptation mechanism for contrast-based dose optimization in XPCI, which can be applied to the GBI and EI systems. Our proposed method explores a new direction for dose reduction, and may also be further extended to other types of XPCI systems and information retrieval algorithms.

Edge-illumination x-ray phase-contrast imaging

Although early demonstration dates back to the mid-sixties, x-ray phase-contrast imaging (XPCI) became hugely popular in the mid-90s, thanks to the advent of 3rd generation synchrotron facilities. Its ability to reveal object features that had so far been considered invisible to x-rays immediately suggested great potential for applications across the life and the physical sciences, and an increasing number of groups worldwide started experimenting with it. At that time, it looked like a synchrotron facility was strictly necessary to perform XPCI with some degree of efficiency-the only alternative being micro-focal sources, the limited flux of which imposed excessively long exposure times. However, new approaches emerged in the mid-00s that overcame this limitation, and allowed XPCI implementations with conventional, non-micro-focal x-ray sources. One of these approaches showing particular promise for 'real-world' applications is edge-illumination XPCI: this article describes the key steps in its evolution in the context of contemporary developments in XPCI research, and presents its current state-of-the-art, especially in terms of transition towards practical applications.

Fourier-based interpretation and noise analysis of the moments of small-angle x-ray scattering in grating-based x-ray phase contrast imaging

In grating-based x-ray phase contrast imaging, Fourier component analysis (FCA) is usually recognized as a gold standard to retrieve the contrasts including attenuation, phase and dark-field, since it is well-established on wave optics and is of high computational efficiency. Meanwhile, an alternative approach basing on the particle scattering theory is being developed and can provide similar contrasts with FCA by calculating multi-order moments of deconvolved small-angle x-ray scattering, so called as multi-order moment analysis (MMA). Although originated from quite different physics theories, the high consistency between the contrasts retrieved by FCA and MMA implies us that there may be some intrinsic connections between them, which has not been fully revealed to the best of our knowledge. In this work, we present a Fourier-based interpretation of MMA and conclude that the contrasts retrieved by MMA are actually the weighted compositions of Fourier coefficients, which means MMA delivers similar physical information as FCA. Based on the recognized cosine model, we also provide a truncated analytic MMA method, and its computational efficiency can be hundreds of times faster than the original deconvolution-based MMA method. Moreover, a noise analysis for our proposed truncated method is also conducted to further evaluate its performances. The results of numerical simulation and physical experiments support our analyses and conclusions.

Review and experimental verification of x-ray dark-field signal interpretations with respect to quantitative isotropic and anisotropic dark-field computed tomography

Talbot(-Lau) interferometric x-ray and neutron dark-field imaging has, over the past decade, gained substantial interest for its ability to provide insights into a sample's microstructure below the imaging resolution by means of ultra small angle scattering effects. Quantitative interpretations of such images depend on models of the signal origination process that relate the observable image contrast to underlying physical processes. A review of such models is given here and their relation to the wave optical derivations by Yashiro et al and Lynch et al as well as to small angle scattering is discussed. Fresnel scaling is introduced to explain the characteristic distance dependence observed in cone beam geometries. Moreover, a model describing the anisotropic signals of fibrous objects is derived. The Yashiro-Lynch model is experimentally verified both in radiographic and tomographic imaging in a monochromatic synchrotron setting, considering both the effects of material and positional dependence of the resulting dark-field contrast. The effect of varying sample-detector distance on the dark-field signal is shown to be non-negligible for tomographic imaging, yet can be largely compensated for by symmetric acquisition trajectories. The derived orientation dependence of the dark-field contrast of fibrous materials both with respect to variations in autocorrelation width and scattering cross section is experimentally validated using carbon fiber reinforced rods.

Grating-based spectral X-ray dark-field imaging for correlation with structural size measures

X-ray dark-field (XDF) imaging accesses information on the small-angle scattering properties of the sample. With grating interferometry, the measured scattering signal is related to the sample's autocorrelation function, which was previously demonstrated for simple samples, such as mono-dispersed microspheres for which the autocorrelation function is mathematically given. However, in potential clinical applications of XDF imaging, complex microstructures, such as lung parenchyma are under investigation. Their bahaviour in XDF imaging is not yet known and no mathematical description of the autocorrelation function is derived so far. In this work we demonstrate the previously established correlation of the XDF data of complex sample structures with their autocorrelation function to be impractical. Furthermore, we propose an applicable correlation between XDF and the sample's structural parameter on the basis of mean chord length, a medically-approved measure for alveolar structure, known to be affected by structural lung diseases. Our findings reveal a correlation between energy-dependent XDF imaging and the sample's mean chord length. By that, a connection between a medical measure for alveoli and XDF is achieved, which is particularly important regarding potential future XDF lung imaging applications for the assessment of alveoli size in diagnostic lung imaging.

Single-shot, omni-directional x-ray scattering imaging with a laboratory source and single-photon localization

Omni-directional, ultra-small-angle x-ray scattering imaging provides a method to measure the orientation of micro-structures without having to resolve them. In this letter, we use single-photon localization with the Timepix3 chip to demonstrate, to the best of our knowledge, the first laboratory-based implementation of single-shot, omni-directional x-ray scattering imaging using the beam-tracking technique. The setup allows a fast and accurate retrieval of the scattering signal using a simple absorption mask. We suggest that our new approach may enable faster laboratory-based tensor tomography and could be used for energy-resolved x-ray scattering imaging.

The trigonometric orthogonality of phase-stepping curves in grating-based x-ray phase-contrast imaging: Integral property and its implications for noise optimization

PURPOSE

Grating-based x-ray phase-contrast imaging (GPCI) is a promising technique for clinical applications as it can provide two newly emerging imaging modalities (differential phase-contrast and dark-field contrast) in addition to the conventional absorption contrast. As far, phase-stepping strategy is the most commonly used approach in GPCI to indirectly acquire differential phase-contrast and dark-field contrast. It is known that the obtained phase-stepping curves (PSCs) have the cosine property and the convolution property, leading to two types of information retrieval approaches in literature: the Fourier component analysis and the multi-order moment analysis. The purpose of this paper is to derive a new property of PSCs and apply the property to noise optimization for information retrieval.

METHODS

Based on the cosine expression of the flat PSC without the sample and the well-established convolution relationship between the flat PSC and the sample PSC, we reveal an important integral property of PSCs: the inner product of PSCs and an arbitrary function contains only zero-order and first-order components in the Fourier series. Furthermore, we apply the property to the direct multi-order moment analysis and propose a set of generalized forms including an optimal one in the presence of noise.

RESULTS

To validate the effectiveness of our analysis, we compared the simulated and real experiment results retrieved by the original direct multi-order moment analysis with the ones retrieved by our proposed noise-optimal form. A significant improvement of noise performance by our method is observed and the improvement ratio in differential phase-contrast is consistent with our theoretical calculation (39.2%).

CONCLUSIONS

In this paper, we reveal a new integral property of the acquired PSCs with and without samples in GPCI, which can be applied to information retrieval approaches like the direct multi-order moment analysis. Then we optimize these approaches to improve the noise performance, offering great potentials of dose reduction in practical applications.

Applying the Fokker-Planck equation to grating-based x-ray phase and dark-field imaging

X-ray imaging has conventionally relied upon attenuation to provide contrast. In recent years, two complementary modalities have been added; (a) phase contrast, which can capture low-density samples that are difficult to see using attenuation, and (b) dark-field x-ray imaging, which reveals the presence of sub-pixel sample structures. These three modalities can be accessed using a crystal analyser, a grating interferometer or by looking at a directly-resolved grid, grating or speckle pattern. Grating and grid-based methods extract a differential phase signal by measuring how far a feature in the illumination has been shifted transversely due to the presence of a sample. The dark-field signal is extracted by measuring how the visibility of the structured illumination is decreased, typically due to the presence of sub-pixel structures in a sample. The strength of the dark-field signal may depend on the grating period, the pixel size and the set-up distances, and additional dark-field signal contributions may be seen as a result of strong phase effects or other factors. In this paper we show that the finite-difference form of the Fokker-Planck equation can be applied to describe the drift (phase signal) and diffusion (dark-field signal) of the periodic or structured illumination used in phase contrast x-ray imaging with gratings, in order to better understand any cross-talk between attenuation, phase and dark-field x-ray signals. In future work, this mathematical description could be used as a basis for new approaches to the inverse problem of recovering both phase and dark-field information.

Quantitative characterization of ex vivo breast tissue via x-ray phase-contrast tomography

BACKGROUND

Grating-based X-ray phase-contrast imaging (GPCI) has received growing interests in recent years due to its high capability of visualizing soft tissue. Breast imaging is one of the most promising candidates for the first clinical application of this imaging modality.

OBJECTIVE

In this work, quantitative breast tissue characterization based on GPCI computed tomography (CT) is investigated with a laboratory X-ray tube through a comparison between attenuation-based CT images and phase-contrast CT images.

METHODS

The Hounsfield units (HU) scale was introduced to phase-contrast images due to its wide application in clinical medicine. In this work, instead of water, plastic cylinders composed of polyethylene terephthalate (PET) was treated as the calibration material. An alternative test-retest reliability (TRR) was presented to evaluate the repeatability of GPCI. Comparison between attenuation-based CT imaging and GPCI CT imaging was operated with the use of statistical analysis methods like histograms and receiver operating characteristic (ROC) curves.

RESULTS

The determined mean TRR related to cylinders is slightly larger in phase-contrast imaging (0.93) than that in attenuation-based imaging (0.89). With respect to distinguishing breast tissues, the AUC (area under curve) values of ROC curves of phase-contrast images are higher than that of attenuation-based images.

CONCLUSIONS

An ex vivo study of GPCI shows that it is a stable imaging modality for visualizing the breast tissue with good repeatability, and that it could be of potential for the diagnosis of breast cancer as well.