Anisotropic X-Ray Dark-Field Tomography: A Continuous Model and its Discretization

Small-angle scattering tensor tomography algorithm for robust reconstruction of complex textures

The development of small-angle scattering tensor tomography has enabled the study of anisotropic nanostructures in a volume-resolved manner. It is of great value to have reconstruction methods that can handle many different nanostructural symmetries. For such a method to be employed by researchers from a wide range of backgrounds, it is crucial that its reliance on prior knowledge about the system is minimized, and that it is robust under various conditions. Here, a method is presented that employs band-limited spherical functions to enable the reconstruction of reciprocal-space maps of a wide variety of nanostructures. This method has been thoroughly tested and compared with existing methods in its ability to retrieve known reciprocal-space maps, as well as its robustness to changes in initial conditions, using both simulations and experimental data. It has also been evaluated for its computational performance. The anchoring of this method in a framework of integral geometry and linear algebra highlights its possibilities and limitations.

Fiber Orientation Estimation from X-ray Dark Field Images of Fiber Reinforced Polymers Using Constrained Spherical Deconvolution

The properties of fiber reinforced polymers are strongly related to the length and orientation of the fibers within the polymer matrix, the latter of which can be studied using X-ray computed tomography (XCT). Unfortunately, resolving individual fibers is challenging because they are small compared to the XCT voxel resolution and because of the low attenuation contrast between the fibers and the surrounding resin. To alleviate both problems, anisotropic dark field tomography via grating based interferometry (GBI) has been proposed. Here, the fiber orientations are extracted by applying a Funk-Radon transform (FRT) to the local scatter function. However, the FRT suffers from a low angular resolution, which complicates estimating fiber orientations for small fiber crossing angles. We propose constrained spherical deconvolution (CSD) as an alternative to the FRT to resolve fiber orientations. Instead of GBI, edge illumination phase contrast imaging is used because estimating fiber orientations with this technique has not yet been explored. Dark field images are generated by a Monte Carlo simulation framework. It is shown that the FRT cannot estimate the fiber orientation accurately for crossing angles smaller than 70∘, while CSD performs well down to a crossing angle of 50∘. In general, CSD outperforms the FRT in estimating fiber orientations.

Fast acquisition protocol for X-ray scattering tensor tomography

Microstructural information over an entire sample is important to understand the macroscopic behaviour of materials. X-ray scattering tensor tomography facilitates the investigation of the microstructural organisation in statistically large sample volumes. However, established acquisition protocols based on scanning small-angle X-ray scattering and X-ray grating interferometry inherently require long scan times even with highly brilliant X-ray sources. Recent developments in X-ray diffractive optics towards circular pattern arrays enable fast single-shot acquisition of the sample scattering properties with 2D omnidirectional sensitivity. X-ray scattering tensor tomography with the use of this circular grating array has been demonstrated. We propose here simple yet inherently rapid acquisition protocols for X-ray scattering tensor tomography leveraging on these new optical elements. Results from both simulation and experimental data, supported by a null space analysis, suggest that the proposed acquisition protocols are not only rapid but also corroborate that sufficient information for the accurate volumetric reconstruction of the scattering properties is provided. The proposed acquisition protocols will build the basis for rapid inspection and/or time-resolved tensor tomography of the microstructural organisation over an extended field of view.

Simulation study towards quantitative X-ray and neutron tensor tomography regarding the validity of linear approximations of dark-field anisotropy

Tensor tomography is fundamentally based on the assumption of a both anisotropic and linear contrast mechanism. While the X-ray or neutron dark-field contrast obtained with Talbot(-Lau) interferometers features the required anisotropy, a preceding detailed study of dark-field signal origination however found its specific orientation dependence to be a non-linear function of the underlying anisotropic mass distribution and its orientation, especially challenging the common assumption that dark-field signals are describable by a function over the unit sphere. Here, two approximative linear tensor models with reduced orientation dependence are investigated in a simulation study with regard to their applicability to grating based X-ray or neutron dark-field tensor tomography. By systematically simulating and reconstructing a large sample of isolated volume elements covering the full range of feasible anisotropies and orientations, direct correspondences are drawn between the respective tensors characterizing the physically based dark-field model used for signal synthesization and the mathematically motivated simplified models used for reconstruction. The anisotropy of freely rotating volume elements is thereby confirmed to be, for practical reconstruction purposes, approximable both as a function of the optical axis’ orientation or as a function of the interferometer’s grating orientation. The eigenvalues of the surrogate models’ tensors are found to exhibit fuzzy, yet almost linear relations to those of the synthesization model. Dominant orientations are found to be recoverable with a margin of error on the order of magnitude of 1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}∘. Although the input data must adequately address the full orientation dependence of dark-field anisotropy, the present results clearly support the general feasibility of quantitative X-ray dark-field tensor tomography within an inherent yet acceptable statistical margin of uncertainty.

Tomography of dark-field scatter including single-exposure Moiré fringe analysis with X-ray biprism interferometry-A simulation study

PURPOSE

In this work, we present tomographic simulations of a new hardware concept for X-ray phase-contrast interferometry wherein the phase gratings are replaced with an array of Fresnel biprisms, and Moiré fringe analysis is used instead of "phase stepping" popular with grating-based setups.

METHODS

Projections of a phantom consisting of four layers of parallel carbon microfibers is simulated using wave optics representation of X-ray electromagnetic waves. Simulated projections of a phantom with preferential scatter perpendicular to the direction of the fibers are performed to analyze the extraction of small-angle scatter from dark-field projections for the following: (1) biprism interferometry using Moiré fringe analysis; (2) grating interferometry using phase stepping with eight grating steps; and (3) grating interferometry using Moiré fringe analysis. Dark-field projections are modeled as projections of voxel intensities represented by a fixed finite vector basis set of scattering directions. An iterative MLEM algorithm is used to reconstruct, from simulated projection data, the coefficients of a fixed set of seven basis vectors at each voxel representing the small-angle scatter distribution.

RESULTS

Results of reconstructed vector coefficients are shown comparing the three simulated imaging configurations. The single-exposure Moiré fringe analysis shows not only an increase in noise because of one-eighth the number of projection samples but also is obtained with less dose and faster acquisition times. Furthermore, replacing grating interferometry with biprism interferometry provides better contrast-to-noise.

CONCLUSION

The simulations demonstrate the feasibility of the reconstruction of dark-field data acquired with a biprism interferometry system. With the potential of higher fringe visibility, biprism interferometry with Moiré fringe analysis might provide equal or better image quality to that of phase stepping methods with less imaging time and lower dose.

X-ray dark-field tomography reveals tooth cracks

Cracked tooth syndrome (CTS) is a common clinical finding for teeth, it affects about 5% of all adults each year. The finding of CTS is favored by several risk factors such as restorations, bruxism, occlusion habits, and age. Treatment options range, depending on the severity, from no treatment at all to tooth extraction. Early diagnosis of CTS is crucial for optimal treatment and symptom reduction. There is no standard procedure for an evidence-based diagnosis up to date. The diagnosis is a challenge by the fact that the symptoms, including pain and sensitivity to temperature stimuli, cannot be clearly linked to the disease. Commonly used visual inspection does not provide in-depth information and is limited by the resolution of human eyes. This can be overcome by magnifying optics or contrast enhancers, but the diagnosis will still strongly rely on the practicians experience. Other methods are symptom reproduction with percussions, thermal pulp tests or bite tests. Dental X-ray radiography, as well as computed tomography, rarely detect cracks as they are limited in resolution. Here, we investigate X-ray dark-field tomography (XDT) for the detection of tooth microcracks. XDT simultaneously detects X-ray small-angle scattering (SAXS) in addition to the attenuation, whereas it is most sensitive to the micrometer regime. Since SAXS originates from gradients in electron density, the signal is sensitive to the sample morphology. Microcracks create manifold interfaces which lead to a strong signal. Therefore, it is possible to detect structural changes originating from subpixel-sized structures without directly resolving them. Together with complementary attenuation information, which visualizes comparatively large cracks, cracks are detected on all length-scales for a whole tooth in a non-destructive way. Hence, this proof-of principle study on three ex-vivo teeth shows the potential of X-ray scattering for evidence-based detection of cracked teeth.

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.

The versatile X-ray beamline of the Munich Compact Light Source: design, instrumentation and applications

Inverse Compton scattering provides means to generate low-divergence partially coherent quasi-monochromatic, i.e. synchrotron-like, X-ray radiation on a laboratory scale. This enables the transfer of synchrotron techniques into university or industrial environments. Here, the Munich Compact Light Source is presented, which is such a compact synchrotron radiation facility based on an inverse Compton X-ray source (ICS). The recent improvements of the ICS are reported first and then the various experimental techniques which are most suited to the ICS installed at the Technical University of Munich are reviewed. For the latter, a multipurpose X-ray application beamline with two end-stations was designed. The beamline's design and geometry are presented in detail including the different set-ups as well as the available detector options. Application examples of the classes of experiments that can be performed are summarized afterwards. Among them are dynamic in vivo respiratory imaging, propagation-based phase-contrast imaging, grating-based phase-contrast imaging, X-ray microtomography, K-edge subtraction imaging and X-ray spectroscopy. Finally, plans to upgrade the beamline in order to enhance its capabilities are discussed.

Validation study of small-angle X-ray scattering tensor tomography

Small-angle scattering tensor tomography (SASTT) is a recently developed technique able to tomographically reconstruct the 3D reciprocal space from voxels within a bulk volume. SASTT extends the concept of X-ray computed tomography, which typically reconstructs scalar values, by reconstructing a tensor per voxel, which represents the local nanostructure 3D organization. In this study, the nanostructure orientation in a human trabecular-bone sample obtained by SASTT was validated by sectioning the sample and using 3D scanning small-angle X-ray scattering (3D sSAXS) to measure and analyze the orientation from single voxels within each thin section. Besides the presence of cutting artefacts from the slicing process, the nanostructure orientations obtained with the two independent methods were in good agreement, as quantified with the absolute value of the dot product calculated between the nanostructure main orientations obtained in each voxel. The average dot product per voxel over the full sample containing over 10 000 voxels was 0.84, and in six slices, in which fewer cutting artefacts were observed, the dot product increased to 0.91. In addition, SAXS tensor tomography not only yields orientation information but can also reconstruct the full 3D reciprocal-space map. It is shown that the measured anisotropic scattering for individual voxels was reproduced from the SASTT reconstruction in each voxel of the 3D sample. The scattering curves along different 3D directions are validated with data from single voxels, demonstrating SASTT's potential for a separate analysis of nanostructure orientation and structural information from the angle-dependent intensity distribution.

X-ray Fokker-Planck equation for paraxial imaging

The Fokker-Planck equation can be used in a partially-coherent imaging context to model the evolution of the intensity of a paraxial x-ray wave field with propagation. This forms a natural generalisation of the transport-of-intensity equation. The x-ray Fokker-Planck equation can simultaneously account for both propagation-based phase contrast, and the diffusive effects of sample-induced small-angle x-ray scattering, when forming an x-ray image of a thin sample. Two derivations are given for the Fokker-Planck equation associated with x-ray imaging, together with a Kramers-Moyal generalisation thereof. Both equations are underpinned by the concept of unresolved speckle due to unresolved sample micro-structure. These equations may be applied to the forward problem of modelling image formation in the presence of both coherent and diffusive energy transport. They may also be used to formulate associated inverse problems of retrieving the phase shifts due to a sample placed in an x-ray beam, together with the diffusive properties of the sample. The domain of applicability for the Fokker-Planck and Kramers-Moyal equations for paraxial imaging is at least as broad as that of the transport-of-intensity equation which they generalise, hence the technique is also expected to be useful for paraxial imaging using visible light, electrons and neutrons.

Brain Connectivity Exposed by Anisotropic X-ray Dark-field Tomography

To understand the interaction of different parts of the human brain it is essential to know how they are connected. Such connections are predominantly related to the brain's white matter, which forms the neuronal pathways, the axons. These axons, also referred to as nerve fibers, have a size on the micrometer scale and are therefore too small to be imaged by standard X-ray systems. In this paper, we use a grating interferometer and a method based on Anisotropic X-ray Dark-field Tomography (AXDT) with the goal to generate a three-dimensional tomographic reconstruction of these functional structures. A first preclinical survey shows that we successfully reconstruct the orientations of the brain fibers connectivity with our approach.

Design of Acquisition Schemes and Setup Geometry for Anisotropic X-ray Dark-Field Tomography (AXDT)

Anisotropic X-ray Dark-field Tomography (AXDT) is a new imaging technique for reconstructing the three-dimensional scattering profile within a sample using the dark-field signal measured in an X-ray grating interferometry setup. As in any tomographic measurement, the acquisition geometry plays a key role in the accurate reconstruction of the scattering information. More- over, the anisotropic nature of the dark-field signal poses additional challenges for designing the acquisition protocols. In this work, we present an efficient approach to measure scattering orientations spread over the unit sphere and prove its efficacy using the knowledge from conventional tomography. In addition, we conclude (using analytical and experimental results) that placing the gratings such that the grating bars make an angle of 45 degrees with respect to the vertical direction is the optimal setup configuration for AXDT.