Comparison of two x-ray phase-contrast imaging methods with a microfocus source

In propagation-based X-ray phase-contrast (PB XPC) imaging, the measured image contains a mixture of absorption- and phase-contrast. To obtain separate images of the projected absorption and phase (i.e., refractive) properties of a sample, phase retrieval methods can be employed. It has been suggested that phase-retrieval can always improve image quality in PB XPC imaging. However, when objective (task-based) measures of image quality are employed, this is not necessarily true and phase retrieval can be detrimental. In this work, signal detection theory is utilized to quantify the performance of a Hotelling observer (HO) for detecting a known signal in a known background. Two cases are considered. In the first case, the HO acts directly on the measured intensity data. In the second case, the HO acts on either the retrieved phase or absorption image. We demonstrate that the performance of the HO is superior when acting on the measured intensity data. The loss of task-specific information induced by phase-retrieval is quantified by computing the efficiency of the HO as the ratio of the test statistic signal-to-noise ratio (SNR) for the two cases. The effect of the system geometry on this efficiency is systematically investigated. Our findings confirm that phase-retrieval can impair signal detection performance in XPC imaging.

X-ray phase contrast imaging (XPCI) shows promising potential for analyzing blood vessels affected by a cardiovascular disease and may become a valuable diagnostic tool in angiography. Identifying the most suitable XPCI method for providing optimal visibility at clinically acceptable doses is critical. This study presents a two-dimensional computational simulation in Geant4 to quantitatively compare two XPCI techniques—Inline XPCI and Edge Illumination Single Mask (EISM)—in visualizing atherosclerotic plaques within the projected thickness image of the lumen. The angiographic specimen analyzed was modeled as three concentric cylinders, with half of the inner cylinder simulating the lumen of the blood vessel and the other half simulating an atherosclerotic plaque, the intermediate cylinder simulating the vessel tissue, and the outer cylinder simulating the tissue around the blood vessel. The computational setup consisted of a W-anode polychromatic x-ray source and a photon counting detector, with six spectra and three sample-to-detector distances (magnifications) considered. To enable a direct comparison with Inline XPCI, a multi-material phase retrieval procedure was adapted to work within the EISM framework. Estimates of absorbed dose were also performed for each XPCI method. The results demonstrate that Inline XPCI achieves a superior contrast-to-noise ratio (CNR) as compared to EISM in both analyzed regions of the recovered image, across all spectra and magnifications, while also exhibiting lower absorbed doses. The difference in CNR between the two XPCI methods decreases as the voltage and magnification increases. While Inline XPCI shows superior contrast in the plaque region for all voltages and magnifications studied, EISM excels in recovering the projected thickness in the free-plaque region, achieving better contrast at low voltages for any magnification.

We have developed a totally incoherent, non-interferometric x-ray phase contrast imaging (XPCI) method. This is based on the edge illumination (EI) concept developed at the ELETTRA synchrotron in Italy in the late `90s. The method was subsequently adapted to the divergent beam generated by a conventional source, by replicating it for every detector line through suitable masks. The method was modelled both with the simplified ray-tracing and with the more rigorous wave-optics approach, and in both cases excellent agreement with the experimental results was found. The wave-optics model enabled assessing the methods' coherence requirements, showing that they are at least an order of magnitude more relaxed than in other methods, without this having negative consequences on the phase sensitivity. Our masks have large pitches (up to 50 times larger than in grating interferometry, for example), which allows for manufacturing through standard lithography, scalability, cost-effectiveness and easiness to align. When applied to a polychromatic and divergent beam generated by a conventional source, the method enables the detection of strong phase effects also with uncollimated, unapertured sources with focal spots of up to 100 μm, compatible with the state-of-the-art in mammography. When used at synchrotrons, it enables a contrast increase of orders of magnitude over other methods. Robust phase retrieval was proven for both coherent and incoherent sources, and additional advantages are compatibility with high x-ray energies and easy implementation of phase sensitivity in two directions simultaneously. This paper briefly summarizes these achievements and reviews some of the key results.

Propagation based X-ray phase-contrast imaging technique for the microstructure analysis of biological soft tissues

Sample radiation damage is one of the main drawbacks limiting applications of the x-ray phase-contrast imaging method. Recently, for x-ray grating-based phase contrast imaging, the reverse projection (RP) method has been proposed by Zhu et al. [Proc. Natl. Acad. Sci. U.S.A. 107, 13576-13581 (2010)]. Compared to the conventional phase stepping technique, the RP method allows a strong reduction of the exposure time and minimization of the errors induced by mechanical vibrations. However, so far, it has only been considered for parallel beam illumination, typical of a synchrotron radiation source. In this paper, a generalization of the RP method from parallel beam geometry to fan beam geometry was presented and discussed. Within parallel beam geometry, the RP method utilizes the conjugate characteristic between reverse projection images. Nevertheless, this characteristic is not directly satisfied for fan beam geometry. In this study, a phantom composed of known materials was constructed and the projection images of the phantom were calculated in a fan beam geometry. By considering single ray, the conjugate images of the projection images were derived from the projection dataset. After that, using the modified RP method the authors retrieved phase and absorption information from paired images. Extracted phase and absorption information of the phantom were in good agreement with theoretical values. Additionally, the slice reconstruction was performed and the results turned out to be in the authors' expectation. Theoretical calculations and numerical simulations confirm both feasibility and validity of the RP method under fan beam illumination. Because this method is simple, fast, and releases a relatively low dose, the authors believe that this research is very useful for the x-ray phase contrast imaging applications in clinical diagnosis, bioresearch, and industrial nondestructive testing.

We report a new preparative method for providing contrast through reduction in electron density that is uniquely suited for propagation-based differential x-ray phase contrast imaging. The method, which results in an air or fluid filled vasculature, makes possible visualization of the smallest microvessels, roughly down to 15 µm, in an excised murine liver, while preserving the tissue for subsequent histological workup. We show the utility of spatial frequency filtering for increasing the visibility of minute features characteristic of phase contrast imaging, and the capability of tomographic reconstruction to reveal microvessel structure and three-dimensional visualization of the sample. The effect of water evaporation from livers during x-ray imaging on the visibility of blood vessels is delineated. The deformed vascular tree in a cancerous murine liver is imaged.

Propagation-based imaging is an X-ray phase-contrast method that does not require any additional optical elements in the imaging system. Instead, phase-contrast effects are observed as measurable intensity variations after an incident X-ray wave field has propagated some distance between the object and plane of detection. Subsequent phase-retrieval processing of the recorded images is a necessary step in order to obtain quantitative information and can greatly increase the image signal-to-noise ratio, in particular in computed tomography. Most commonly this is done by approximating the object as a single material, which cannot always be justified. An example of this is the presence of bones and soft-tissues in biomedical imaging, which motivated the work presented here. We developed the combination of spectral X-ray imaging and propagation-based X-ray phasecontrast imaging to overcome the limitations of a single-material approximation. The addition of spectral data enables phase-retrieval on multi-material objects, as well as the use of processing methods well-known from conventional dual-energy imaging, such as material decomposition, virtual monochromatic images and effective atomic-number mapping. Here, we present first computed tomography results of spectral propagation-based phase-contrast. We discuss how the phase retrieval process can be interpreted as low- and high-pass filters. Lastly, we take a look at current and future technological challenges, and the potential application areas of spectral propagation-based phase-contrast imaging.

Objective. X-ray phase contrast imaging is a promising technique for future clinical diagnostic as it can provide enhanced contrast in soft tissues compared to traditional x-ray attenuation-contrast imaging. However, the strict requirements on the x-ray coherence and the precise alignment of optical elements limit its applications towards clinical use. To solve this problem, mesh-based x-ray phase contrast imaging method with one hexagonal mesh is proposed for easy alignment and better image visualization. Approach. The mesh produces structured illuminations and the detector captures its distortions to reconstruct the absorption, differential phase contrast (DPC) and dark-field (DF) images of the sample. In this work, we fabricated a hexagonal mesh to simultaneously retrieve DPC and DF signals in three different directions with single shot. A phase retrieval algorithm to obtain artifacts-free phase from DPC images with three different directions is put forward and false color dark-field image is also reconstructed with tri-directional images. Mesh-shifting method based on this hexagonal mesh modulator is also proposed to reconstruct images with better image quality at the expense of increased dose. Main results. In numerical simulations, the proposed hexagonal mesh outperforms the traditional square mesh in image evaluation metrics performance and false color visualization with the same radiation dose. The experimental results demonstrate its feasiblity in real imaging systems and its advantages in quantitive imaging and better visualization. The proposed hexagonal mesh is easy to fabricate and can be successfully applied to x-ray source with it spot size up to 300 μm. Significance. This work opens new possibilities for quantitative x-ray non-destructive imaging and may also be instructive for research fields such as x-ray structured illumination microscopy (SIM), x-ray spectral imaging and x-ray phase contrast and dark-field computed tomography (CT).

An x-ray microtomography station implemented at the X-ray Applied Crystallography Laboratory of the State University of Campinas is described. The station is based on a propagation based phase contrast imaging setup with a microfocus source and digital x-ray area detectors. Due to its simplicity, this setup is ideal for fast, high resolution imaging and microtomography of small biological specimens and materials research samples. It can also be coupled to gratings to use and develop new techniques as the harmonic spatial coherent imaging, which allow scattering contrast imaging. Details of the experimental setup, equipment, and software integration are described. Test microtomography for setup commissioning and characterization is shown. We conclude that phase contrast enhanced x-ray imaging and microtomography with resolution below 5 μm voxel size are possible and data sets as wide as 2000 × 2000 × 2000 voxels are obtained with this instrumentation.

Propagation-based X-ray phase-contrast imaging (PBI) is a powerful nondestructive imaging technique that can reveal the internal detailed structures in weakly absorbing samples. Extending PBI to CT (PBCT) enables high-resolution and high-contrast 3D visualization of microvasculature, which can be used for the understanding, diagnosis and therapy of diseases involving vasculopathy, such as cardiovascular disease, stroke and tumor. However, the long scan time for PBCT impedes its wider use in biomedical and preclinical microvascular studies. To address this issue, a novel CT reconstruction algorithm for PBCT is presented that aims at shortening the scan time for microvascular samples by reducing the number of projections while maintaining the high quality of reconstructed images. The proposed algorithm combines the filtered backprojection method into the iterative reconstruction framework, and a weighted guided image filtering approach (WGIF) is utilized to optimize the intermediate reconstructed images. Notably, the homogeneity assumption on the microvasculature sample is adopted as prior knowledge, and therefore, a prior image of microvasculature structures can be acquired by a k-means clustering approach. Then, the prior image is used as the guided image in the WGIF procedure to effectively suppress streaking artifacts and preserve microvasculature structures. To evaluate the effectiveness and capability of the proposed algorithm, simulation experiments on 3D microvasculature numerical phantom and real experiments with CT reconstruction on the microvasculature sample are performed. The results demonstrate that the proposed algorithm can, under noise-free and noisy conditions, significantly reduce the artifacts and effectively preserve the microvasculature structures on the reconstructed images and thus enables it to be used for clear and accurate 3D visualization of microvasculature from few-projection data. Therefore, for 3D visualization of microvasculature, the proposed algorithm can be considered an effective approach for reducing the scan time required by PBCT.

Visualization of temporomandibular joint articular cartilage using synchrotron-radiation X-ray phase-contrast imaging

The question as to whether the multiple least squares regression of X‐ray fluorescence analytical data used in an earlier paper and based on an equation derived for excitation by monochromatic radiation can be considered as valid with sufficient accuracy in an experimental system using polychromatic radiation as source of excitation was examined. The results of computer simulation of the process of X‐ray fluorescence show that although definite differences in the intensities of fluorescent radiation arise as a result of the substitution of polychromatic primary radiation for a monochromatic source, the flexibility of the particular equation used for regression eliminates such deviations by incorporating them in the experimental regression coefficients. The overall accuracy of the regression was found to be satisfactory by any practical criteria.

Phase-Contrast Tomography (PCT) is becoming an important technique for nondestructive, in-situ characterization of soft and hard condensed matter. This thesis sheds light on our progress in developing novel tomographic reconstruction algorithms in combination with in-situ experimental approaches that employ propagation-based PCT. These developments are aimed to improve the spatial and temporal resolutions of in-situ studies of the soft and hard condensed matter. The mechanical properties of materials such as steel are, to a large extent, defined by its microstructure. Propagation-based PCT (i.e. holotomography) can be used to visualize the microstructures associated with small variations in density or composition in the bulk of the specimen. This imaging technique has the advantage of being non-destructive and allows to carry out in-situ, time-resolved observations of miscrostructural changes during dynamic processes such as solid-state phase transformations. Sufficient spatial and temporal resolution can be achieved only when an extremely bright x-ray source is used. In order to test whether a novel table-top x-ray source, the MIRRORCLE-6X, can be used for implementation of a propagation-based PCT system, we have carried out a series of experiments. The characteristics of a propagation-based PCI system that employs the MIRRORCLE-6X were analyzed. In our further research we have focused on the experiments that could be performed at the European Synchrotron Radiation Facility (ESRF) in Grenoble, France. We reported on the investigation of the cementite microstructure in carbon steel. In this work we present the reconstruction of the three-dimensional morphology of cementite grains in the bulk of steel using a non-destructive imaging technique – x-ray Phase-Contrast Tomography (PCT). Complementary information about the crystalline structure of the ferrite grains surrounding the cementite is obtained using x-ray Diffraction-Contrast Tomography (DCT). The work was continued by a time-resolved x-ray PCT investigation of the austenite grain growth during the ferrite-to-austenite phase transformation in low-carbon steel. Experimental data allows to visualize the evolution of the three-dimensional morphology of the austenite-ferrite interfaces with a spatial resolution on the order of 1 ?m and a time resolution of less than 9 minutes. The in-situ measurements show that the theory for grain growth during diffusional solid-state phase transformations needs to be extended in order to describe the drastically lower mobility that we observed for flat interfaces compared to curved interfaces. We observe that some flat interfaces are completely immobile, which is not predicted by the state-of-the-art theories and showed - for the first time - that new computational model for grain growth are needed. Quantitative image reconstruction based on PCT data requires solving two inverse problems: phase retrieval and tomographic reconstruction. In most cases, the problem of phase retrieval is ill-posed and requires that some type of a prior knowledge about the reconstructed image (i.e. particular regularization method) is used. We present a novel algebraic approach suitable for phase retrieval using various linear models. In this approach Total Variation (TV) minimization is used for regularization of the linearized inverse problem by promoting the solution with a sparse gradient magnitude (i.e. piece-wise constant solution). In that case prior knowledge about the reconstructed image may allow to (partially) recover the unknown spatial frequencies that are undefined by the experimental data. The problem of tomographic reconstruction based on PCT data can be solved using a similar approach. We conclude our investigation of the image reconstruction algorithms by introducing several algebraic approaches suitable for tomographic reconstruction based on PCT data. In these approaches TV minimization is used to find a regularized solution to an underdetermined linear system based on a linearized representation of the PCT. When the density of the reconstructed object is piece-wise constant (or close to it) a virtually artifact-free solution can be computed for the tomographic problem. The proposed approach can also be used to radically improve the accuracy of the tomographic reconstruction based on incomplete data (e.g. small number of projections) or data with low signal to noise ratio.

Phase retrieval for X-ray differential phase contrast radiography with knowledge transfer learning from virtual differential absorption model

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.