Effect of insufficient temporal coherence on visibility contrast in X-ray grating interferometry

Empirical beam hardening and ring artifact correction for x-ray grating interferometry (EBHC-GI)

PURPOSE

Talbot-Lau grating interferometry enables the use of polychromatic x-ray sources, extending the range of potential applications amenable to phase contrast imaging. However, these sources introduce beam hardening effects not only from the samples but also from the gratings. As a result, grating inhomogeneities due to manufacturing imperfections can cause spectral nonuniformity artifacts when used with polychromatic sources. Consequently, the different energy dependencies of absorption, phase, and visibility contrasts impose challenges that so far have limited the achievable image quality. The purpose of this work was to develop and validate a correction strategy for grating-based x-ray imaging that accounts for beam hardening generated from both the imaged object and the gratings.

METHODS

The proposed two-variable polynomial expansion strategy was inspired by work performed to address beam hardening from a primary modulator. To account for the multicontrast nature of grating interferometry, this approach was extended to each contrast to obtain three sets of correction coefficients, which were determined empirically from a calibration scan. The method's feasibility was demonstrated using a tabletop Talbot-Lau grating interferometer micro-computed tomography (CT) system using CT acquisitions of a water sample and a silicon sample, representing low and high atomic number materials. Spectral artifacts such as cupping and ring artifacts were quantified using mean squared error (MSE) from the beam-hardening-free target image and standard deviation within a reconstructed image of the sample. Finally, the model developed using the water sample was applied to a fixated murine lung sample to demonstrate robustness for similar materials.

RESULTS

The water sample's absorption CT image was most impacted by spectral artifacts, but following correction to decrease ring artifacts, an 80% reduction in MSE and 57% reduction in standard deviation was observed. The silicon sample created severe artifacts in all contrasts, but following correction, MSE was reduced by 94% in absorption, 96% in phase, and 90% in visibility images. These improvements were due to the removal of ring artifacts for all contrasts and reduced cupping in absorption and phase images and reduced capping in visibility images. When the water calibration coefficients were applied to the lung sample, ring artifacts most prominent in the absorption contrast were eliminated.

CONCLUSIONS

The described method, which was developed to remove artifacts in absorption, phase, and normalized visibility micro-CT images due to beam hardening in the system gratings and imaged object, reduced the MSE by up to 96%. The method depends on calibrations that can be performed on any system and does not require detailed knowledge of the x-ray spectrum, detector energy response, grating attenuation properties and imperfections, or the geometry and composition of the imaged object.

High-speed rotating device for X-ray tomography with 10 ms temporal resolution

The temporal resolution of X-ray tomography, using a synchrotron radiation X-ray source, has been improved to millisecond order in recent years. However, the sample must be rotated at a speed of more than a few thousand revolutions per minute, which makes it difficult to control the environment around the sample. In this study, a high-speed rotation device has been developed, comprising two synchronized coaxial motors movable along the direction of the axis, which can stretch or compress the rotating sample. Using this device, tomograms of breaking rubber were successfully obtained at a temporal resolution of 10 ms.

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.

Wave optics simulation of grating-based X-ray phase-contrast imaging using 4D Mouse Whole Body (MOBY) phantom

PURPOSE

Demonstrate realistic simulation of grating-based x-ray phase-contrast imaging (GB-XPCI) using wave optics and the four-dimensional Mouse Whole Body (MOBY) phantom defined with non-uniform rational B-splines (NURBS).

METHODS

We use a full-wave approach, which uses wave optics for x-ray wave propagation from the source to the detector. This forward imaging model can be directly applied to NURBS-defined numerical phantoms such as MOBY. We assign the material properties (attenuation coefficient and electron density) of each model part using the data for adult human tissues. The Poisson noise is added to the simulated images based on the calculated photon flux at each pixel.

RESULTS

We simulate the intensity images of the MOBY phantom for eight different grating positions. From the simulated images, we calculate the absorption, differential phase, and normalized visibility contrast images. We also predict how the image quality is affected by different exposure times.

CONCLUSIONS

GB-XPCI can be simulated with the full-wave approach and a realistic numerical phantom defined with NURBS.

Modeling the Moiré fringe visibility of Talbot-Lau X-ray grating interferometry for single-frame multi-contrast imaging

Moiré fringe method in X-ray grating interferometry is characterized by its advantage to obtain multi-contrast data through single-frame imaging. However, how the visibility of the Moiré fringe is influenced by the system parameters, such as the misalignment angle, still lacks investigation, although it closely relates to the signal-to-noise ratio of the image data. In this paper, a simplified model of the Moiré fringe visibility is presented, which reveals that the Moiré fringe visibility can be expressed as the product of a misalignment-angle-related "sinc" function and a relatively independent factor. The following experimental results further suggest that the crosstalk between the detector pixels in the direction perpendicular to the Moiré fringe orientation is another main cause for the visibility reduction of the Moiré fringes.

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.

Hard X-ray imaging microscopy with self-imaging phenomenon

The self-imaging phenomenon referred to as the Talbot effect in the field of optics was discovered by H.F. Talbot in the 1830s, and is now widely used for imaging using not only visible light but also X-rays, electrons, neutrons, and matter waves. In this review, the author introduces the current progress being made in hard-X-ray imaging microscopy based on the self-imaging phenomenon. Hard-X-ray imaging microscopy is a promising technique for non-destructively visualizing internal structures in specimens with a spatial resolution up to a few tens of nanometers. The use of the self-imaging phenomenon makes it possible to realize highly sensitive phase-contrast X-ray imaging microscopes. These approaches have several advantages over conventional X-ray imaging microscopes, including the widely used Zernike X-ray phase-contrast microscopes, and can provide a powerful way of quantitative visualization with a high spatial resolution and a high sensitivity even for thick specimens.

New reconstruction method for few-view grating-based phase-contrast imaging via dictionary learning

Grating-based phase-contrast is a hot topic in recent years owing to its excellent imaging contrast capability on soft tissues. Although it is compatible with conventional X-ray tubes and applicable in many fields, long scanning time, and high radiation dose obstruct its wider use in clinical and medical fields, especially for computed tomography applications. In this study, we solve this challenge by reducing the projection views and compensating the loss of reconstruction quality through dual-dictionary learning algorithm. The algorithm is implemented in two steps. First, estimated high-quality absorption images are obtained from the first dual-quality dictionary learning, which uses the correspondence between high-quality images and low-quality ones reconstructed from highly under-sampled data. Then, the second absorption-phase dual-modality dictionary learning is adopted to yield both estimated phase and absorption images, resulting in complementary information for both modality images. Afterwards the absorption and phase images are gradually improved in iterative reconstructions. By using SSIM RMSE measurements and visual assessment for enlarged regions of interest, our proposed method can improve the resolution of these two modality images and recover smaller structures, as compared to conventional methods.