X-ray computed tomography (CT) is a crucial tool for non-invasive medical diagnosis that uses differences in materials' attenuation coefficients to generate contrast and provide 3D information. Grating-based dark-field-contrast X-ray imaging is an innovative technique that utilizes small-angle scattering to generate additional co-registered images with additional microstructural information. While it is already possible to perform human chest dark-field radiography, it is assumed that its diagnostic value increases when performed in a tomographic setup. However, the susceptibility of Talbot-Lau interferometers to mechanical vibrations coupled with a need to minimize data acquisition times has hindered its application in clinical routines and the combination of X-ray dark-field imaging and large field-of-view (FOV) tomography in the past. In this work, we propose a processing pipeline to address this issue in a human-sized clinical dark-field CT prototype. We present the corrective measures that are applied in the employed processing and reconstruction algorithms to mitigate the effects of vibrations and deformations of the interferometer gratings. This is achieved by identifying spatially and temporally variable vibrations in air reference scans. By translating the found correlations to the sample scan, we can identify and mitigate relevant fluctuation modes for scans with arbitrary sample sizes. This approach effectively eliminates the requirement for sample-free detector area, while still distinctly separating fluctuation and sample information. As a result, samples of arbitrary dimensions can be reconstructed without being affected by vibration artifacts. To demonstrate the viability of the technique for human-scale objects, we present reconstructions of an anthropomorphic thorax phantom.

X-ray dark-field imaging enables a spatially-resolved visualization of ultra-small-angle X-ray scattering. Using phantom measurements, we demonstrate that a material's effective dark-field signal may be reduced by modification of the visibility spectrum by other dark-field-active objects in the beam. This is the dark-field equivalent of conventional beam-hardening, and is distinct from related, known effects, where the dark-field signal is modified by attenuation or phase shifts. We present a theoretical model for this group of effects and verify it by comparison to the measurements. These findings have significant implications for the interpretation of dark-field signal strength in polychromatic measurements.

Grating-based X-ray phase-contrast and in particular dark-field radiography are promising new imaging modalities for medical applications. Currently, the potential advantage of dark-field imaging in early-stage diagnosis of pulmonary diseases in humans is being investigated. These studies make use of a comparatively large scanning interferometer at short acquisition times, which comes at the expense of a significantly reduced mechanical stability as compared to tabletop laboratory setups. Vibrations create random fluctuations of the grating alignment, causing artifacts in the resulting images. Here, we describe a novel maximum likelihood method for estimating this motion, thereby preventing these artifacts. It is tailored to scanning setups and does not require any sample-free areas. Unlike any previously described method, it accounts for motion in between as well as during exposures.

Dark-field radiography, a new X-ray imaging method, has recently been applied to human chest imaging for the first time. It employs conventional X-ray devices in combination with a Talbot-Lau interferometer with a large field of view, providing both attenuation and dark-field radiographs. It is well known that sample scatter creates artifacts in both modalities. Here, we demonstrate that also X-ray scatter generated by the interferometer as well as detector crosstalk create artifacts in the dark-field radiographs, in addition to the expected loss of spatial resolution. We propose deconvolution-based correction methods for the induced artifacts. The kernel for detector crosstalk is measured and fitted to a model, while the kernel for scatter from the analyzer grating is calculated by a Monte-Carlo simulation. To correct for scatter from the sample, we adapt an algorithm used for scatter correction in conventional radiography. We validate the obtained corrections with a water phantom. Finally, we show the impact of detector crosstalk, scatter from the analyzer grating and scatter from the sample and their successful correction on dark-field images of a human thorax.

Grating-based phase- and dark-field-contrast X-ray imaging is a novel technology that aims to extend conventional attenuation-based X-ray imaging by unlocking two additional contrast modalities. The so called phase-contrast and dark-field channels provide enhanced soft tissue contrast and additional microstructural information. Accessing this additional information comes at the expense of a more intricate measurement setup and necessitates sophisticated data processing. A big challenge for translating grating-based dark-field computed tomography to medical applications lies in minimizing the data acquisition time. While a continuously moving detector is ideal, it prohibits conventional phase stepping techniques that require multiple projections under the same angle with different grating positions. One solution to this problem is the so-called sliding window processing approach that is compatible with continuous data acquisition. However, conventional sliding window techniques lead to crosstalk-artifacts between the three image channels, if the projection of the sample moves too fast on the detector within a processing window. In this work we introduce a new interpretation of the phase retrieval problem for continuous acquisitions as a demodulation problem. In this interpretation, we identify the origin of the crosstalk-artifacts as partially overlapping modulation side bands. Furthermore, we present three algorithmic extensions that improve the conventional sliding-window-based phase retrieval and mitigate crosstalk-artifacts. The presented algorithms are tested in a simulation study and on experimental data from a human-scale dark-field CT prototype. In both cases they achieve a substantial reduction of the occurring crosstalk-artifacts.

X-ray computed tomography (CT) is an invaluable imaging technique for non-invasive medical diagnosis. However, for soft tissue in the human body the difference in attenuation is inherently small. Grating-based X-ray phase-contrast is a relatively novel imaging method which detects additional interaction mechanisms between photons and matter, namely refraction and small-angle scattering, to generate additional images with different contrast. The experimental setup involves a Talbot-Lau interferometer whose susceptibility to mechanical vibrations hindered acquisition schemes suitable for clinical routine in the past. We present a processing pipeline to identify spatially and temporally variable fluctuations occurring in an interferometer installed on a continuously rotating clinical CT gantry. The correlations of the vibrations in the modular grating setup are exploited to identify a small number of relevant fluctuation modes, allowing for a sample reconstruction free of vibration artifacts.

X-ray computed tomography (CT) is an important non-destructive imaging technique, particularly in clinical diagnostics. Even with the latest innovations like dual-energy and photon-counting CT, the image contrast is solely generated from attenuation in the tissue. An extension - fully compatible with these novelties - is dark-field CT, which retrieves an additional, so-called dark-field contrast. Unlike the attenuation channel, the dark-field channel is sensitive to tissue microstructure and porosity below the resolution of the imaging system, which allows additional insights into the health of the lung tissue or the structure of calcifications. The potential clinical value has been demonstrated in several preclinical studies and recently also in radiography patient studies. Just recently the first dark-field CT for the human body was established at the Technical University of Munich and in this paper, we discuss the performance of this prototype. We evaluate the interferometer components and the imposed challenges that the integration into the CT gantry brings by comparing the results to simulations and measurements at a laboratory setup. The influence of the clinical X-ray source on the Talbot-Lau interferometer and the impact of vibrations, which are immanent on the clinical CT gantry, are analyzed in detail to reveal their characteristic frequencies and origin. A beam hardening correction is introduced as an important step to adapt to the poly-chromatic spectrum and make quantitative dark-field imaging possible. We close with an analysis of the image resolution and the applied patient dose, and conclude that the performance is sufficient to suggest initial patient studies using the presented dark-field CT system.

Computed tomography (CT) as an important clinical diagnostics method can profit from extension with dark-field imaging, as it is currently restricted to X-rays' attenuation contrast only. Dark-field imaging allows access to more tissue properties, such as micro-structural texture or porosity. The up-scaling process to clinical scale is complex because several design constraints must be considered. The two most important ones are that the finest grating is limited by current manufacturing technology to a [Formula: see text] period and that the interferometer should fit into the CT gantry with minimal modifications only. In this work we discuss why an inverse interferometer and a triangular G1 profile are advantageous and make a compact and sensitive interferometer implementation feasible. Our evaluation of the triangular grating profile reveals a deviation in the interference pattern compared to standard grating profiles, which must be considered in the subsequent data processing. An analysis of the grating orientation demonstrates that currently only a vertical layout can be combined with cylindrical bending of the gratings. We also provide an in-depth discussion, including a new simulation approach, of the impact of the extended X-ray source spot which can lead to large performance loss and present supporting experimental results. This analysis reveals a vastly increased sensitivity to geometry and grating period deviations, which must be considered early in the system design process.

This work presents an approach for image reconstruction in clinical low-dose tomography that combines principles from sparse signal processing with ideas from deep learning. First, we describe sparse signal representation in terms of dictionaries from a statistical perspective and interpret dictionary learning as a process of aligning the distribution that arises from a generative model with the empirical distribution of true signals. As a result, we can see that sparse coding with learned dictionaries resembles a specific variational autoencoder, where the encoder is a sparse coding algorithm and the decoder is a linear function. Next, we show that dictionary learning can also benefit from computational advancements introduced in the context of deep learning, such as parallelism and stochastic optimization. Finally, we show that regularization by dictionaries achieves competitive performance in computed tomography reconstruction compared to state-of-the-art model-based and data-driven approaches, while being unsupervised with respect to tomographic data.

Grating-based differential phase-contrast imaging has proven to be feasible with conventional X-ray sources. The polychromatic spectrum generally limits the performance of the interferometer but benefit can be gained with an energy-sensitive detector. In the presented work, we employ the energy-discrimination capability to correct for phase-wrapping artefacts. We propose to use the phase shifts, which are measured in distinct energy bins, to estimate the optimal phase shift in the sense of maximum likelihood. We demonstrate that our method is able to correct for phase-wrapping artefacts, to improve the contrast-to-noise ratio and to reduce beam hardening due to the modelled energy dependency. The method is evaluated on experimental data which are measured with a laboratory Talbot-Lau interferometer equipped with a conventional polychromatic X-ray source and an energy-sensitive photon-counting pixel detector. Our work shows, that spectral imaging is an important step to move differential phase-contrast imaging closer to pre-clinical and clinical applications, where phase wrapping is particularly problematic.

Photon-counting detector technology has enabled the first experimental investigations of energy-resolved computed tomography (CT) imaging and the potential use for K-edge imaging. However, limitations in regards to detecter technology have been imposing a limit to effective count rates. As a consequence, this has resulted in high noise levels in the obtained images given scan time limitations in CT imaging applications. It has been well recognized in the area of low-dose imaging with conventional CT that iterative image reconstruction provides a superior signal to noise ratio compared to traditional filtered backprojection techniques. Furthermore, iterative reconstruction methods also allow for incorporation of a roughness penalty function in order to make a trade-off between noise and spatial resolution in the reconstructed images. In this work, we investigate statistically-principled iterative image reconstruction from material-decomposed sinograms in spectral CT. The proposed reconstruction algorithm seeks to minimize a penalized likelihood-based cost functional, where the parameters of the likelihood function are estimated by computing the Fisher information matrix associated with the material decomposition step. The performance of the proposed reconstruction method is quantitatively investigated by use of computer-simulated and experimental phantom data. The potential for improved K-edge imaging is also demonstrated in an animal experiment.

The mathematical analysis of exact filtered back-projection algorithms is strictly related to Radon inversion. We show how filter-lines can be defined for the helical trajectory, which serve for the extraction of contributions of particular kinds of Radon-planes. Due to the Fourier-slice theorem, Radon-planes with few intersections with the helix are associated with low-frequency contributions to transversal slices. This insight leads to different applications of the new method. The application presented here enables the incorporation of an arbitrary amount of redundant data in an approximate way. This means that the back-projection is not restricted to an n-Pi interval. A detailed mathematical analysis, in which we demonstrate how the defined filter-lines work, concludes this paper.

In this paper, we formulate a reconstruction algorithm for an n-Pi acquisition, where n can be any positive odd integer. The algorithm is a generalization of the method presented in (Bontus et al. 2003). It is based on the results obtained by Katsevich (2004). For the algorithm, different sets of filter-lines have to be defined. We describe the variation of these lines along the detector in some detail, before we discuss, how the method gives all Radon-plane contributions the correct weighting. The different sets of filter-lines are all contained within the n-Pi window, such that a practical realization is possible. Reconstruction results, which we present in the final section, show convincing image quality.

A new class of acquisition schemes for helical cone-beam computed tomography (CB-CT) scanning is introduced, and their effect on the reconstruction methods is analyzed. These acquisition schemes are based on a new detector shape that is bounded by the helix. It will be shown that the data acquired with these schemes are compatible with exact reconstruction methods, and the adaptation of exact reconstruction algorithms to the new acquisition geometry is described. At the same time, the so-called PI-sufficiency condition is fulfilled. Moreover, a good fit to the acquisition requirements of the various medical applications of cone-beam CT is achieved. In contrast to other helical cone-beam acquisition and reconstruction methods, the n-PI-method introduced in this publication allows for variable pitches of the acquisition helix. This additional feature will introduce a higher flexibility into the acquisition protocols of future medical cone-beam scanners. An approximative n-PI-filtered backprojection (n-PI-FBP) reconstruction method is presented and verified. It yields convincing image quality.