Monte Carlo model of a polychromatic laboratory based edge illumination x-ray phase contrast system

Edge illumination x-ray phase contrast simulations using the CAD-ASTRA toolbox

Edge illumination x-ray phase contrast imaging (XPCI) provides increased contrast for low absorbing materials compared to attenuation images and sheds light on the material microstructure through dark field contrast. To apply XPCI in areas such as non-destructive testing and inline inspection, where scanned samples are increasingly compared to simulated reference images, accurate and efficient simulation software is required. However, currently available simulators rely on expensive Monte Carlo techniques or wave-optics frameworks, resulting in long simulation times. Furthermore, these simulators are often not optimized to work with computer-aided design (CAD) models, a common and memory-efficient method to represent manufactured objects, hindering their integration in an inspection pipeline. In this work, we address these shortcomings by introducing an edge illumination XPCI simulation framework built upon the recently developed CAD-ASTRA toolbox. CAD-ASTRA allows for the efficient simulation of x-ray projections from CAD models through GPU-accelerated ray tracing and supports ray refraction in a geometric optics framework. The edge illumination implementation is validated and its performance is benchmarked against GATE, a state-of-the-art Monte Carlo simulator, revealing a simulation speed increase of up to three orders of magnitude, while maintaining high accuracy in the resulting images.

Grating designs for cone beam edge illumination X-ray phase contrast imaging: a simulation study

Edge illumination is an emerging X-ray phase contrast imaging technique providing attenuation, phase and dark field contrast. Despite the successful transition from synchrotron to lab sources, the cone beam geometry of lab systems limits the effectiveness of using conventional planar gratings. The non-parallel incidence of X-rays introduces shadowing effects, worsening with increasing cone angle. To overcome this limitation, several alternative grating designs can be considered. In this paper, the effectiveness of three alternative designs is compared to conventional gratings using numerical simulations. Improvements in flux and contrast are discussed, taking into account practical considerations concerning the implementation of the designs.

Virtual grating approach for Monte Carlo simulations of edge illumination-based x-ray phase contrast imaging

The design of new x-ray phase contrast imaging setups often relies on Monte Carlo simulations for prospective parameter studies. Monte Carlo simulations are known to be accurate but time consuming, leading to long simulation times, especially when many parameter variations are required. This is certainly the case for imaging methods relying on absorbing masks or gratings, with various tunable properties, such as pitch, aperture size, and thickness. In this work, we present the virtual grating approach to overcome this limitation. By replacing the gratings in the simulation with virtual gratings, the parameters of the gratings can be changed after the simulation, thereby significantly reducing the overall simulation time. The method is validated by comparison to explicit grating simulations, followed by representative demonstration cases.

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.

Semi-classical Monte Carlo algorithm for the simulation of X-ray grating interferometry

Traditional simulation techniques such as wave optics methods and Monte Carlo (MC) particle transport cannot model both interference and inelastic scattering phenomena within one framework. Based on the rules of quantum mechanics to calculate probabilities, we propose a new semi-classical MC algorithm for efficient and simultaneous modeling of scattering and interference processes. The similarities to MC particle transport allow the implementation as a flexible c++ object oriented extension of EGSnrc-a well-established MC toolkit. In addition to previously proposed Huygens principle based transport through optics components, new variance reduction techniques for the transport through gratings are presented as transport options to achieve the required improvement in speed and memory costs necessary for an efficient exploration (system design-dose estimations) of the medical implementation of X-ray grating interferometry (GI), an emerging imaging technique currently subject of tremendous efforts towards clinical translation. The feasibility of simulation of interference effects is confirmed in four academic cases and an experimental table-top GI setup. Comparison with conventional MC transport show that deposited energy features of EGSnrc are conserved.

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.

Large field of view, fast and low dose multimodal phase-contrast imaging at high x-ray energy

X-ray phase contrast imaging (XPCI) is an innovative imaging technique which extends the contrast capabilities of 'conventional' absorption based x-ray systems. However, so far all XPCI implementations have suffered from one or more of the following limitations: low x-ray energies, small field of view (FOV) and long acquisition times. Those limitations relegated XPCI to a 'research-only' technique with an uncertain future in terms of large scale, high impact applications. We recently succeeded in designing, realizing and testing an XPCI system, which achieves significant steps toward simultaneously overcoming these limitations. Our system combines, for the first time, large FOV, high energy and fast scanning. Importantly, it is capable of providing high image quality at low x-ray doses, compatible with or even below those currently used in medical imaging. This extends the use of XPCI to areas which were unpractical or even inaccessible to previous XPCI solutions. We expect this will enable a long overdue translation into application fields such as security screening, industrial inspections and large FOV medical radiography - all with the inherent advantages of the XPCI multimodality.