Achromatic approach to phase-based multi-modal imaging with conventional X-ray sources

T staging esophageal tumors with x rays

With histopathology results typically taking several days, the ability to stage tumors during interventions could provide a step change in various cancer interventions. X-ray technology has advanced significantly in recent years with the introduction of phase-based imaging methods. These have been adapted for use in standard labs rather than specialized facilities such as synchrotrons, and approaches that enable fast 3D scans with conventional x-ray sources have been developed. This opens the possibility to produce 3D images with enhanced soft tissue contrast at a level of detail comparable to histopathology, in times sufficiently short to be compatible with use during surgical interventions. In this paper we discuss the application of one such approach to human esophagi obtained from esophagectomy interventions. We demonstrate that the image quality is sufficiently high to enable tumor T staging based on the x-ray datasets alone. Alongside detection of involved margins with potentially life-saving implications, staging tumors intra-operatively has the potential to change patient pathways, facilitating optimization of therapeutic interventions during the procedure itself. Besides a prospective intra-operative use, the availability of high-quality 3D images of entire esophageal tumors can support histopathological characterization, from enabling "right slice first time" approaches to understanding the histopathology in the full 3D context of the surrounding tumor environment.

Volumetric High-Resolution X-Ray Phase-Contrast Virtual Histology of Breast Specimens With a Compact Laboratory System

The assessment of margin involvement is a fundamental task in breast conserving surgery to prevent recurrences and reoperations. It is usually performed through histology, which makes the process time consuming and can prevent the complete volumetric analysis of large specimens. X-ray phase contrast tomography combines high resolution, sufficient penetration depth and high soft tissue contrast, and can therefore provide a potential solution to this problem. In this work, we used a high-resolution implementation of the edge illumination X-ray phase contrast tomography based on "pixel-skipping" X-ray masks and sample dithering, to provide high definition virtual slices of breast specimens. The scanner was originally designed for intra-operative applications in which short scanning times were prioritised over spatial resolution; however, thanks to the versatility of edge illumination, high-resolution capabilities can be obtained with the same system simply by swapping x-ray masks without this imposing a reduction in the available field of view. This makes possible an improved visibility of fine tissue strands, enabling a direct comparison of selected CT slices with histology, and providing a tool to identify suspect features in large specimens before slicing. Combined with our previous results on fast specimen scanning, this works paves the way for the design of a multi-resolution EI scanner providing intra-operative capabilities as well as serving as a digital pathology system.

The effect of a variable focal spot size on the contrast channels retrieved in edge-illumination X-ray phase contrast imaging

Multi-modal X-ray imaging allows the extraction of phase and dark-field (or “Ultra-small Angle Scatter”) images alongside conventional attenuation ones. Recently, scan-based systems using conventional sources that can simultaneously output the above three images on relatively large-size objects have been developed by various groups. One limitation is the need for some degree of spatial coherence, achieved either through the use of microfocal sources, or by placing an absorption grating in front of an extended source. Both these solutions limit the amount of flux available for imaging, with the latter also leading to a more complex setup with additional alignment requirements. Edge-illumination partly overcomes this as it was proven to work with focal spots of up to 100 micron. While high-flux, 100 micron focal spot sources do exist, their comparatively large footprint and high cost can be obstacles to widespread translation. A simple solution consists in placing a single slit in front of a large focal spot source. We used a tunable slit to study the system performance at various effective focal spot sizes, by extracting transmission, phase and dark-field images of the same specimens for a range of slit widths. We show that consistent, repeatable results are obtained for varying X-ray statistics and effective focal spot sizes. As the slit width is increased, the expected reduction in the raw differential phase peaks is observed, compensated for in the retrieval process by a broadened sensitivity function. This leads to the same values being correctly retrieved, but with a slightly larger error bar i.e. a reduction in phase sensitivity. Concurrently, a slight increase in the dark-field signal is also observed.

X-ray computed tomography for 3D plant imaging

X-ray computed tomography (CT) is a valuable tool for 3D imaging of plant tissues and organs. Applications include the study of plant development and organ morphogenesis, as well as modeling of transport processes in plants. Some challenges remain, however, including attaining higher contrast for easier quantification, increasing the resolution for imaging subcellular features, and decreasing image acquisition and processing time for high-throughput phenotyping. In addition, phase contrast, multispectral, dark-field, soft X-ray, and time-resolved imaging are emerging. At the same time, a large amount of 3D image data are becoming available, posing challenges for data management. We review recent advances in the area of X-ray CT for plant imaging, and describe opportunities for using such images for studying transport processes in plants.

Optimization of multipoint phase retrieval in edge illumination X-ray imaging: A theoretical and experimental analysis

PURPOSE

In this work, an analytical model describing the noise in the retrieved three contrast channels, transmission, refraction, and ultra small-angle scattering, obtained with edge illumination X-ray phase-based imaging system is presented and compared to experimental data.

METHODS

In EI, images acquired at different displacements of the presample mask (i.e., different illumination levels referred to as points on the "illumination curve"), followed by pixel-wise curve fitting, are exploited to quantitatively retrieve the three contrast channels. Therefore, the noise in the final image will depend on the error associated with the fitting process. We use a model based on the derivation of the standard error on fitted parameters, which relies on the calculation of the covariance matrix, to estimate the noise and the cross-channel correlation as a function of the position of the sampling points. In particular, we investigated the most common cases of 3 and 5 sampling points. In addition, simulations have been used to better understand the role of the integration time for each sampling point. Finally, the model is validated by comparison with the experimental data acquired with an edge illumination setup based on a tungsten rotating anode X-ray source and a photon counting detector.

RESULTS

We found a good match between the predictions of the model and the experimental data. In particular, for the investigated cases, an arrangement of the sampling points leading to minimum noise and cross-channel correlation can be found. Simulations revealed that, given a fixed overall scanning time, its distribution into the smallest possible number of sampling points needed for phase retrieval leads to minimum noise thanks to higher statistics per point.

CONCLUSIONS

This work presents an analytical model describing the noise in the various contrast channels retrieved in edge illumination as a function of the illumination curve sampling. In particular, an optimal sampling scheme leading to minimum noise has been determined for the case where 3 or 5 sampling points are used, which represent two of the most common acquisition schemes. In addition, the correlation between noise in the different channels and the role of the number of points and exposure time have been also investigated. In general, our results suggest a series of procedures that should be followed in order to optimize the experimental acquisitions.

Comparison of microstructural imaging of the root canal isthmus using propagation-based X-ray phase-contrast and absorption micro-computed tomography

Clear and complete microstructural imaging of the root canal isthmus is an important part of pathological investigations in research and clinical practice. X-ray micro-computed tomography (μCT) is a widely used non-destructive imaging technique, which allows for distortion-free three-dimensional (3D) visualisation. While absorption μCT typically has poor contrast resolution for observing the root canal isthmus, especially for weak-absorbing tissues, propagation-based X-ray phase-contrast imaging (PBI) is a powerful imaging method, which in its combination with μCT (PB-PCμCT) enables high-resolution and high-contrast microstructural imaging of the weak-absorbing tissues in samples. To investigate the feasibility and ability of PB-PCμCT in microstructural imaging of the root canal isthmus, conventional absorption μCT and PB-PCμCT experiments were performed. The two-dimensional (2D) and 3D comparison results demonstrated that, compared to absorption μCT, PB-PCμCT has the ability to image the root canal isthmus more clearly and completely, and thus, it has great potential to serve as a valuable tool for biomedical and preclinical studies on the root canal isthmus.

A Preliminary Investigation into the Use of Edge Illumination X-ray Phase Contrast Micro-CT for Preclinical Imaging

Purpose

To enable a preliminary assessment of the suitability of edge illumination (EI) x-ray phase contrast (XPC) micro x-ray computed tomography (micro-CT) to preclinical imaging. Specifically, to understand how different acquisition schemes and their combination with dedicated data processing affect contrast-to-noise ratio (CNR) and spatial resolution, while providing control over scan time and radiation dose delivery.

Procedures

Deceased mice (n = 3) were scanned with an EI XPC micro-CT setup operated under different settings, leading to scan times between 18 h and 13 min. For the shortest scan, the entrance dose was measured with a calibrated PTW 23344 ion chamber. Different data processing methods were applied, retrieving either separate attenuation and phase images, or hybrid (combined attenuation and phase) images. A quantitative comparison was performed based on CNR and spatial resolution measurements for a soft tissue interface.

Results

All phase-based images have led to a higher CNR for the considered soft tissue interface than the attenuation image, independent of scan time. The best relative CNR (a sixfold increase) was observed in one of the hybrid images. Spatial resolution was found to be connected to scan time, with a resolution of approximately 20 μm and 60 μm achieved for the longest and shortest scans, respectively. An entrance dose of approximately 300 mGy was estimated for the scan performed within 13 min.

Conclusions

Despite their preliminary nature, our results suggest that EI XPC bears potential for enhancing the utility of preclinical micro-CT, and, pending further research and development, could ultimately become a valuable technique in this field.

Quantitative X-ray phase contrast computed tomography with grating interferometry : Biomedical applications of quantitative X-ray grating-based phase contrast computed tomography

The ability of biomedical imaging data to be of quantitative nature is getting increasingly important with the ongoing developments in data science. In contrast to conventional attenuation-based X-ray imaging, grating-based phase contrast computed tomography (GBPC-CT) is a phase contrast micro-CT imaging technique that can provide high soft tissue contrast at high spatial resolution. While there is a variety of different phase contrast imaging techniques, GBPC-CT can be applied with laboratory X-ray sources and enables quantitative determination of electron density and effective atomic number. In this review article, we present quantitative GBPC-CT with the focus on biomedical applications.

Fast, non-iterative algorithm for quantitative integration of X-ray differential phase-contrast images

X-ray phase contrast imaging is gaining importance as an imaging tool. However, it is common for X-ray phase detection techniques to be sensitive to the derivatives of the phase. Therefore, the integration of differential phase images is a fundamental step both to access quantitative pixel content and for further analysis such as segmentation. The integration of noisy data leads to artefacts with a severe impact on image quality and on its quantitative content. In this work, an integration method based on the Wiener filter is presented and tested using simulated and real data obtained with the edge illumination differential X-ray phase imaging method. The method is shown to provide high image quality while preserving the quantitative pixel content of the integrated image. In addition, it requires a short computational time making it suitable for large datasets.

Evaluation of a compact multicontrast and multiresolution X-ray phase contrast edge illumination system for small animal imaging

PURPOSE

In this work the performance of a compact multiresolution and multicontrast x-ray phase system based on edge illumination is investigated. It has been designed for small animal imaging and with a limited footprint for ease of deployment in laboratories.

METHODS

The presented edge illumination system is based on a compact microfocus tungsten x-ray source combined with a flat panel detector. The source has a maximum output of 10 W when the minimum spot size of about 15 μm is used. The system has an overall length of 70 cm. A new double sample mask design, obtained by arranging both skipped and nonskipped configurations on the same structure, provides dual resolution capability. To test the system, we carried out computed tomography (CT) scans of a plastic phantom with different source settings using both single-image and multi-image acquisition schemes at different spatial resolutions. In addition, CT scans of an ex-vivo mouse specimen were acquired at the best identified working conditions to demonstrate the application of the presented system to small animal imaging.

RESULTS

We found this system delivers good image quality, allowing for an efficient material separation and improving detail visibility in small animals thanks to the higher signal-to-noise ratio (SNR) of phase contrast with respect to conventional attenuation contrast. The system offers high versatility in terms of spatial resolution thanks to the double sample mask design integrated into a single scanner. The availability of both multi- and single-image acquisition schemes coupled with their dedicated retrieval algorithms, allows different working modes which can be selected based on user preference. Multi-image acquisition provides quantitative separation of the real and imaginary part of the refractive index, however, it requires a long scanning time. On the other hand, the single image approach delivers the best material separation and image quality at all the investigated source settings with a shorter scanning time but at the cost of quantitativeness. Finally, we also observed that the single image approach combined with a high-power x-ray source may result in a fast acquisition protocol compatible with in-vivo imaging.

Simplified retrieval method for Edge Illumination X-ray phase contrast imaging allowing multi-modal imaging with fewer input frames

We present data from an implementation of Edge Illumination (EI) that uses a detector aperture designed for increasing dynamic range, suitable for clinically relevant X-ray energies and demonstrated here using synchrotron radiation. By utilising a sufficiently large crosstalk between pixels, this implementation enables single-scan imaging for phase and absorption, and double-scan for phase, absorption and dark field imaging. The presence of the detector mask enables a direct comparison between conventional EI and beam tracking (BT), which we conduct through Monte Carlo and analytical modelling in the case of a single-scan being used for the retrieval of all three contrasts. In the present case, where the X-ray beam width is comparable to the pixel size, we provide an analysis on best-positioning of the beam on the detector for accurate signal retrieval. Further, we demonstrate an application of this method by distinguishing different concentrations of microbubbles via their dark field signals at high energy using an EI system.

Improved iterative tomographic reconstruction for x-ray imaging with edge-illumination

Iterative tomographic reconstruction has been established as a viable alternative for data analysis in phase-sensitive x-ray imaging based on the edge-illumination principle. However, previously published approaches did not account for drifts of optical elements during a scan, which can lead to artefacts. Up to now, the strategy to reduce these artefacts was to acquire additional intermediate flat field images, which were used to correct the sinograms. Here, we expand the theoretical model to take these effects into account and demonstrate a significant reduction of (ring)-artefacts in the final reconstructions, while allowing for a significant reduction of scan time and dose. We further improve the model by including the capability to reconstruct combined absorption and phase contrast slices, which we experimentally demonstrate to deliver improved contrast to noise ratios compared to previously employed single shot approaches.

Advancement of X-ray radiography using microfocus X-ray source in conjunction with amplitude grating and SOI pixel detector, SOPHIAS

We show how to improve microfocus X-ray radiography by using the SOPHIAS silicon-on-insulator pixel detector in conjunction with an amplitude grating. Single-exposure multi-energy absorption and differential phase contrast imaging was performed using the single amplitude grating method. The sensitivity in differential phase contrast imaging in a two-pixel-pitch setup was enhanced by 39% in comparison with the previously reported method [Rev. Sci. Instrum. 81, 113702 (2010).] by analyzing charge-sharing effects. Small-angle-scattering imaging was also possible in the two-pixel-pitch setup by counting the number of X-ray photons passing the pixel boundaries.

Recent advances in edge illumination x-ray phase-contrast tomography

Edge illumination (EI) is an x-ray phase-contrast imaging technique, exploiting sensitivity to x-ray refraction to visualize features, which are often not detected by conventional absorption-based radiography. The method does not require a high degree of spatial coherence and is achromatic and, therefore, can be implemented with both synchrotron radiation and commercial x-ray tubes. Using different retrieval algorithms, information about an object's attenuation, refraction, and scattering properties can be obtained. In recent years, a theoretical framework has been developed that enables EI computed tomography (CT) and, hence, three-dimensional imaging. This review provides a summary of these advances, covering the development of different image acquisition schemes, retrieval approaches, and applications. These developments constitute an integral part in the transformation of EI CT into a widely spread imaging tool for use in a range of fields.

Hybrid framework for feasible modeling of an edge illumination X-ray phase-contrast imaging system at a human scale

Here, we present a hybrid approach for simulating an edge illumination X-ray phase-contrast imaging (EIXPCi) set-up using graphics processor units (GPU) with a high degree of accuracy. In this study, the applicability of pixel, mesh and non-uniform rational B-splines (NURBS) objects to carry out realistic maps of X-ray phase-contrast distribution at a human scale is accounted for by using numerical anthropomorphic phantoms and a very fast and robust simulation framework which integrates total interaction probabilities along selected X-ray paths. We exploit the mathematical and algorithmic properties of NURBS and describe how to represent human scale phantoms in an edge illumination X-ray phase-contrast model. The presented implementation allows the modeling of a variety of physical interactions of x-rays with different mathematically described objects and the recording of quantities, e.g. path integrals, interaction sites and deposited energies. Furthermore, our efficient, scalable and optimized hybrid Monte Carlo and ray-tracing projector can be used in iterative reconstruction algorithms on multi GPU heterogeneous systems. The preliminary results of our innovative approach show the fine performance of an edge illumination X-ray phase-contrast medical imaging system on various human-like soft tissues with noticeably reduced computation time. Our approach to the EIXPCi modeling confirms that building a true imaging system at a human scale should be possible and the simulations presented here aim at its future development.

Interpretation and Utility of the Moments of Small-Angle X-Ray Scattering Distributions

Small angle x-ray scattering has been proven to be a valuable method for accessing structural information below the spatial resolution limit implied by direct imaging. Here, we theoretically derive the relation that links the subpixel differential phase signal provided by the sample to the moments of scattering distributions accessible by refraction sensitive x-ray imaging techniques. As an important special case we explain the scatter or dark-field contrast in terms of the sample's phase signal. Further, we establish that, for binary phase objects, the nth moment scales with the difference of the refractive index decrement to the power of n. Finally, we experimentally demonstrate the utility of the moments by quantitatively determining the particle sizes of a range of powders with a laboratory-based setup.

X-ray Phase-Contrast Radiography and Tomography with a Multiaperture Analyzer

We present a multiaperture analyzer setup for performing x-ray phase contrast imaging in planar and three-dimensional modalities. The method is based on strongly structuring the x-ray beam with an amplitude modulator, before it reaches the sample, and on a multiaperture analyzing element before detection. A multislice representation of the sample is used to establish a quantitative relation between projection images and the corresponding three-dimensional distributions, leading to successful tomographic reconstruction. Sample absorption, phase, and scattering are retrieved from the measurement of five intensity projections. The method is tested on custom-built phantoms with synchrotron radiation: sample absorption and phase can be reliably retrieved also in combination with strong scatterers, simultaneously attaining high sensitivity and dynamic range.

Edge illumination X-ray phase tomography of multi-material samples using a single-image phase retrieval algorithm

In this paper we present a single-image phase retrieval algorithm for multi-material samples, developed for the edge illumination (EI) X-ray phase contrast imaging method. The theoretical derivation is provided, along with any assumptions made. The algorithm is evaluated quantitatively using both simulated and experimental results from a computed tomography (CT) scan using the EI laboratory implementation. Qualitative CT results are provided for a biological sample containing both bone and soft-tissue. Using a single EI image per projection and knowledge of the complex refractive index, the algorithm can accurately retrieve the interface between a given pair of materials. A composite CT slice can be created by splicing together multiple CT reconstructions, each retrieved for a different pair of materials.

A method for high-energy, low-dose mammography using edge illumination x-ray phase-contrast imaging

Since the breast is one of the most radiosensitive organs, mammography is arguably the area where lowering radiation dose is of the uttermost importance. Phase-based x-ray imaging methods can provide opportunities in this sense, since they do not require x-rays to be stopped in tissue for image contrast to be generated. Therefore, x-ray energy can be considerably increased compared to those usually exploited by conventional mammography. In this article we show how a novel, optimized approach can lead to considerable dose reductions. This was achieved by matching the edge-illumination phase method, which reaches very high angular sensitivity also at high x-ray energies, to an appropriate image processing algorithm and to a virtually noise-free detection technology capable of reaching almost 100% efficiency at the same energies. Importantly, while proof-of-concept was obtained at a synchrotron, the method has potential for a translation to conventional sources.

Small angle x-ray scattering with edge-illumination

Sensitivity to sub-pixel sample features has been demonstrated as a valuable capability of phase contrast x-ray imaging. Here, we report on a method to obtain angular-resolved small angle x-ray scattering distributions with edge-illumination- based imaging utilizing incoherent illumination from an x-ray tube. Our approach provides both the three established image modalities (absorption, differential phase and scatter strength), plus a number of additional contrasts related to unresolved sample features. The complementarity of these contrasts is experimentally validated by using different materials in powder form. As a significant application example we show that the extended complementary contrasts could allow the diagnosis of pulmonary emphysema in a murine model. In support of this, we demonstrate that the properties of the retrieved scattering distributions are consistent with the expectation of increased feature sizes related to pulmonary emphysema. Combined with the simplicity of implementation of edge-illumination, these findings suggest a high potential for exploiting extended sub-pixel contrasts in the diagnosis of lung diseases and beyond.

Asymmetric masks for laboratory-based X-ray phase-contrast imaging with edge illumination

We report on an asymmetric mask concept that enables X-ray phase-contrast imaging without requiring any movement in the system during data acquisition. The method is compatible with laboratory equipment, namely a commercial detector and a rotating anode tube. The only motion required is that of the object under investigation which is scanned through the imaging system. Two proof-of-principle optical elements were designed, fabricated and experimentally tested. Quantitative measurements on samples of known shape and composition were compared to theory with good agreement. The method is capable of measuring the attenuation, refraction and (ultra-small-angle) X-ray scattering, does not have coherence requirements and naturally adapts to all those situations in which the X-ray image is obtained by scanning a sample through the imaging system.

X-ray absorption, phase and dark-field tomography through a beam tracking approach

We present a development of the beam-tracking approach that allows its implementation in computed tomography. One absorbing mask placed before the sample and a high resolution detector are used to track variations in the beam intensity distribution caused by the sample. Absorption, refraction, and dark-field are retrieved through a multi-Gaussian interpolation of the beam. Standard filtered back projection is used to reconstruct three dimensional maps of the real and imaginary part of the refractive index, and of the dark-field signal. While the method is here demonstrated using synchrotron radiation, its low coherence requirements suggest a possible implementation with laboratory sources.