Nonlinear statistical iterative reconstruction for propagation-based phase-contrast tomography

Image quality and scan time optimisation for in situ phase contrast x-ray tomography of the intervertebral disc

In-line phase contrast synchrotron tomography combined with in situ mechanical loading enables the characterisation of soft tissue micromechanics via digital volume correlation (DVC) within whole organs. Optimising scan time is important for reducing radiation dose from multiple scans and to limit sample movement during acquisition. Also, although contrasted edges provided by in-line phase contrast tomography of soft tissues are useful for DVC, the effect of phase contrast imaging on its accuracy has yet to be investigated. Due to limited time at synchrotron facilities, scan parameters are often decided during imaging and their effect on DVC accuracy is not fully understood. Here, we used previously published data of intervertebral disc phase contrast tomography to evaluate the influence of i) fibrous image texture, ii) number of projections, iii) tomographic reconstruction method, and iv) phase contrast propagation distance on DVC results. A greater understanding of how image texture influences optimal DVC tracking was obtained by visualising objective function mapping, enabling tracking inaccuracies to be identified. When reducing the number of projections, DVC was minimally affected by image high frequency noise but with a compromise in accuracy. Iterative reconstruction methods improved image signal-to-noise and consequently significantly lowered DVC displacement uncertainty. Propagation distance was shown to affect DVC accuracy. Consistent DVC results were achieved within a propagation distance range which provided contrast to the smallest scale features, where; too short a distance provided insufficient features to track, whereas too long led to edge effect inconsistencies, particularly at greater deformations. Although limited to a single sample type and image setup, this study provides general guidelines for future investigations when optimising image quality and scan times for in situ phase contrast x-ray tomography of fibrous connective tissues.

X-ray computed tomography using partially coherent Fresnel diffraction with application to an optical fiber

A reconstruction algorithm for partially coherent x-ray computed tomography (XCT) including Fresnel diffraction is developed and applied to an optical fiber. The algorithm is applicable to a high-resolution tube-based laboratory-scale x-ray tomography instrument. The computing time is only a few times longer than the projective counterpart. The algorithm is used to reconstruct, with projections and diffraction, a tilt series acquired at the micrometer scale of a graded-index optical fiber using maximum likelihood and a Bayesian method based on the work of Bouman and Sauer. The inclusion of Fresnel diffraction removes some reconstruction artifacts and use of a Bayesian prior probability distribution removes others, resulting in a substantially more accurate reconstruction.

Material Decomposition Using Spectral Propagation-Based Phase-Contrast X-Ray Imaging

Material decomposition in X-ray imaging uses the energy-dependence of attenuation to digitally decompose an object into specific constituent materials, generally at the cost of enhanced image noise. Propagation-based X-ray phase-contrast imaging is a developing technique that can be used to reduce image noise, in particular from weakly attenuating objects. In this paper, we combine spectral phase-contrast imaging with material decomposition to both better visualize weakly attenuating features and separate them from overlying objects in radiography. We derive an algorithm that performs both tasks simultaneously and verify it against numerical simulations and experimental measurements of ideal two-component samples composed of pure aluminum and poly(methyl methacrylate). Additionally, we showcase first imaging results of a rabbit kitten's lung. The attenuation signal of a thorax, in particular, is dominated by the strongly attenuating bones of the ribcage. Combined with the weak soft tissue signal, this makes it difficult to visualize the fine anatomical structures across the whole lung. In all cases, clean material decomposition was achieved, without residual phase-contrast effects, from which we generate an un-obstructed image of the lung, free of bones. Spectral propagation-based phase-contrast imaging has the potential to be a valuable tool, not only in future lung research, but also in other systems for which phase-contrast imaging in combination with material decomposition proves to be advantageous.

Nanoscale imaging of shale fragments with coherent X-ray diffraction

Despite the abundance of shales in the Earth's crust and their industrial and environmental importance, their microscale physical properties are poorly understood, owing to the presence of many structurally related mineral phases and a porous network structure spanning several length scales. Here, the use of coherent X-ray diffraction imaging (CXDI) to study the internal structure of microscopic shale fragments is demonstrated. Simultaneous wide-angle X-ray diffraction (WAXD) measurement facilitated the study of the mineralogy of the shale microparticles. It was possible to identify pyrite nanocrystals as inclusions in the quartz-clay matrix and the volume of closed unconnected pores was estimated. The combined CXDI-WAXD analysis enabled the establishment of a correlation between sample morphology and crystallite shape and size. The results highlight the potential of the combined CXDI-WAXD approach as an upcoming imaging modality for 3D nanoscale studies of shales and other geological formations via serial measurements of microscopic fragments.

Model-Based Iterative Reconstruction for Propagation-Based Phase-Contrast X-Ray CT including Models for the Source and the Detector

Propagation-based phase-contrast X-ray computed tomography is a valuable tool for high-resolution visualization of biological samples, giving distinct improvements in terms of contrast and dose requirements compared to conventional attenuation-based computed tomography. Due to its ease of implementation and advances in laboratory X-ray sources, this imaging technique is increasingly being transferred from synchrotron facilities to laboratory environments. This however poses additional challenges, such as the limited spatial coherence and flux of laboratory sources, resulting in worse resolution and higher noise levels. Here we extend a previously developed iterative reconstruction algorithm for this imaging technique to include models for the reduced spatial coherence and the signal spreading of efficient scintillator-based detectors directly into the physical forward model. Furthermore, we employ a noise model which accounts for the full covariance statistics of the image formation process. In addition, we extend the model describing the interference effects such that it now matches the formalism of the widely-used single-material phase-retrieval algorithm, which is based on the sample-homogeneity assumption. We perform a simulation study as well as an experimental study at a laboratory inverse Compton source and compare our approach to the conventional analytical approaches. We find that the modeling of the source and the detector inside the physical forward model can tremendously improve the resolution at matched noise levels and that the modeling of the covariance statistics reduces overshoots (i.e. incorrect increase / decrease in sample properties) at the sample edges significantly.

High resolution 3D visualization of the spinal cord in a post-mortem murine model

A crucial issue in the development of therapies to treat pathologies of the central nervous system is represented by the availability of non-invasive methods to study the three-dimensional morphology of spinal cord, with a resolution able to characterize its complex vascular and neuronal organization. X-ray phase contrast micro-tomography enables a high-quality, 3D visualization of both the vascular and neuronal network simultaneously without the need of contrast agents, destructive sample preparations or sectioning. Until now, high resolution investigations of the post-mortem spinal cord in murine models have mostly been performed in spinal cords removed from the spinal canal. We present here post-mortem phase contrast micro-tomography images reconstructed using advanced computational tools to obtain high-resolution and high-contrast 3D images of the fixed spinal cord without removing the bones and preserving the richness of micro-details available when measuring exposed spinal cords. We believe that it represents a significant step toward the in-vivo application.

Metric-guided regularisation parameter selection for statistical iterative reconstruction in computed tomography

As iterative reconstruction in Computed Tomography (CT) is an ill-posed problem, additional prior information has to be used to get a physically meaningful result (close to ground truth if available). However, the amount of influence of the regularisation prior is crucial to the outcome of the reconstruction. Therefore, we propose a scheme for tuning the strength of the prior via a certain image metric. In this work, the parameter is tuned for minimal histogram entropy in selected regions of the reconstruction as histogram entropy is a very basic approach to characterise the information content of data. We performed a sweep over different regularisation parameters showing that the histogram entropy is a suitable metric as it is well behaved over a wide range of parameters. The parameter determination is a feedback loop approach we applied to numerically simulated FORBILD phantom data and verified with an experimental measurement of a micro-CT device. The outcome is evaluated visually and quantitatively by means of root mean squared error (RMSE) and structural similarity (SSIM) for the simulation and visually for the measured sample (no ground truth available). The final reconstructed images exhibit noise-suppressed iterative reconstruction. For both datasets, the optimisation is robust where its initial value is concerned. The parameter tuning approach shows that the proposed metric-driven feedback loop is a promising tool for finding a suitable regularisation parameter in statistical iterative reconstruction.

X-ray imaging of a water bear offers a new look at tardigrade internal anatomy

BACKGROUND

Tardigrades (water bears) are microscopic invertebrates of which the anatomy has been well studied using traditional techniques, but a comprehensive three-dimensional reconstruction has never been performed. In order to close this gap, we employed X-ray computed tomography (CT), a technique that is becoming increasingly popular in zoology for producing high-resolution, three-dimensional (3D) scans of whole specimens. While CT has long been used to scan larger samples, its use in some microscopic animals can be problematic, as they are often too small for conventional CT yet too large for high-resolution, optics-based soft X-ray microscopy. This size gap continues to be narrowed with advancements in technology, with high-resolution imaging now being possible using both large synchrotron devices and, more recently, laboratory-based instruments.

RESULTS

Here we use a recently developed prototype lab-based nano-computed tomography device to image a 152 μm-long tardigrade at high resolution (200-270 nm pixel size). The resulting dataset allowed us to visualize the anatomy of the tardigrade in 3D and analyze the spatial relationships of the internal structures. Segmentation of the major structures of the body enabled the direct measurement of their respective volumes. Furthermore, we segmented every storage cell individually and quantified their volume distribution. We compare our measurements to those from published studies in which other techniques were used.

CONCLUSIONS

The data presented herein demonstrate the utility of CT imaging as a powerful supplementary tool for studies of tardigrade anatomy, especially for quantitative volume measurements. This nanoCT study represents the smallest complete animal ever imaged using CT, and offers new 3D insights into the spatial relationships of the internal organs of water bears.

Propagation-based phase-contrast tomography of a guinea pig inner ear with cochlear implant using a model-based iterative reconstruction algorithm

Propagation-based phase-contrast computed tomography has become a valuable tool for visualization of three-dimensional biological samples, due to its high contrast between materials with similar attenuation properties. However, one of the most-widely used phase-retrieval algorithms imposes a homogeneity assumption onto the sample, which leads to artifacts for numerous applications where this assumption is violated. Prominent examples are biological samples with highly-absorbing implants. Using synchrotron radiation, we demonstrate by the example of a guinea pig inner ear with a cochlear implant electrode, how a recently developed model-based iterative algorithm for propagation-based phase-contrast computed tomography yields distinct benefits for such a task. We find that the model-based approach improves the overall image quality, removes the detrimental influence of the implant and accurately visualizes the cochlea.