Combining Monte Carlo methods with coherent wave optics for the simulation of phase-sensitive X-ray imaging

A hybrid simulation method towards the gamma ray phase contrast imaging for metallic material

A high efficiency simulation method for propagation-based phase-contrast imaging, called directional macro-wavefront (DMWF), is developed with the aim of simulating high-energy phase-contrast imaging. This method takes both Monte Carlo and wave optical propagation into consideration. Traditional wave-optics-based simulation methods for phase-contrast imaging encounter unacceptable computational complexity when high-energy radiation is used. In contrast, this method effectively addresses this issue by using macro-wavefront integration. Several simulation examples using typical parameters of inverse Compton scattering sources are presented to illustrate the excellent energy adaptability and efficiency of the DMWF method. This method provides a more efficient approach for phase-contrast imaging simulations, which will drive the advancement of high-energy phase-contrast imaging.

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.

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.

Evaluation of simulators for x-ray speckle-based phase contrast imaging

X-ray phase contrast imaging (PCI) denotes a group of highly sensitive imaging techniques that permits imaging at scales ranging from nanoscopic to the medical. Recently introduced, speckle-based imaging has seen a rapid development because of its experimental simplicity and its capability to retrieve the refraction, the scattering and the absorption of a sample using a conventional x-ray set-up. Precise simulation would permit to optimise the imaging setups for different applications, but until now works on simulation of x-ray speckle-based PCI have been very few. In this work we evaluate different simulation codes, based on Monte-Carlo, analytical ray-tracing and wave-optics Fresnel propagation. The simulation results are compared to both synchrotron and conventional imaging experiments to permits their validation. We obtain a strong similarity between simulated and experimental data. We discuss the validity and applicability of each approach.

Imaging in thick samples, a phased Monte Carlo radiation transfer algorithm

SIGNIFICANCE

Optical microscopy is characterized by the ability to get high resolution, below 1  μm, high contrast, functional and quantitative images. The use of shaped illumination, such as with lightsheet microscopy, has led to greater three-dimensional isotropic resolution with low phototoxicity. However, in most complex samples and tissues, optical imaging is limited by scattering. Many solutions to this issue have been proposed, from using passive approaches such as Bessel beam illumination to active methods incorporating aberration correction, but making fair comparisons between different approaches has proven to be challenging.

AIM

We present a phase-encoded Monte Carlo radiation transfer algorithm (φMC) capable of comparing the merits of different illumination strategies or predicting the performance of an individual approach.

APPROACH

We show that φMC is capable of modeling interference phenomena such as Gaussian or Bessel beams and compare the model with experiment.

RESULTS

Using this verified model, we show that, for a sample with homogeneously distributed scatterers, there is no inherent advantage to illuminating a sample with a conical wave (Bessel beam) instead of a spherical wave (Gaussian beam), except for maintaining a greater depth of focus.

CONCLUSION

φMC is adaptable to any illumination geometry, sample property, or beam type (such as fractal or layered scatterer distribution) and as such provides a powerful predictive tool for optical imaging in thick samples.

In Silico Phase-Contrast X-Ray Imaging of Anthropomorphic Voxel-Based Phantoms

Propagation-based phase-contrast X-ray imaging is an emerging technique that can improve dose efficiency in clinical imaging. In silico tools are key to understanding the fundamental imaging mechanisms and develop new applications. Here, due to the coherent nature of the phase-contrast effects, tools based on wave propagation (WP) are preferred over Monte Carlo (MC) based methods. WP simulations require very high wave-front sampling which typically limits simulations to small idealized objects. Virtual anthropomorphic voxel-based phantoms are typically provided with a resolution lower than imposed sampling requirements and, thus, cannot be directly translated for use in WP simulations. In the present paper we propose a general strategy to enable the use of these phantoms for WP simulations. The strategy is based on upsampling in the 3D domain followed by projection resulting in high-resolution maps of the projected thickness for each phantom material. These maps can then be efficiently used for simulations of Fresnel diffraction to generate in silico phase-contrast X-ray images. We demonstrate the strategy on an anthropomorphic breast phantom to simulate propagation-based phase-contrast mammography using a laboratory micro-focus X-ray source.

X-ray phase contrast simulation for grating-based interferometry using GATE

The overall importance of x-ray phase contrast (XPC) imaging has grown substantially in the last decades, in particular with the recent advent of compact lab-based XPC systems. For optimizing the experimental XPC setup, as well as benchmarking and testing new acquisition and reconstruction techniques, Monte Carlo (MC) simulations are a valuable tool. GATE, an open source application layer on top of the Geant4 simulation software, is a versatile MC tool primarily intended for various types of medical imaging simulations. To our knowledge, however, there is no GATE-based academic simulation software available for XPC imaging. In this paper, we extend the GATE framework with new physics-based tools for accurate XPC simulations. Our approach combines Monte Carlo simulations in GATE for modelling the x-ray interactions in the sample with subsequent numerical wave propagation, starting from the GATE output.

Wave optics simulation of grating-based X-ray phase-contrast imaging using 4D Mouse Whole Body (MOBY) phantom

PURPOSE

Demonstrate realistic simulation of grating-based x-ray phase-contrast imaging (GB-XPCI) using wave optics and the four-dimensional Mouse Whole Body (MOBY) phantom defined with non-uniform rational B-splines (NURBS).

METHODS

We use a full-wave approach, which uses wave optics for x-ray wave propagation from the source to the detector. This forward imaging model can be directly applied to NURBS-defined numerical phantoms such as MOBY. We assign the material properties (attenuation coefficient and electron density) of each model part using the data for adult human tissues. The Poisson noise is added to the simulated images based on the calculated photon flux at each pixel.

RESULTS

We simulate the intensity images of the MOBY phantom for eight different grating positions. From the simulated images, we calculate the absorption, differential phase, and normalized visibility contrast images. We also predict how the image quality is affected by different exposure times.

CONCLUSIONS

GB-XPCI can be simulated with the full-wave approach and a realistic numerical phantom defined with NURBS.

Towards Monte Carlo simulation of X-ray phase contrast using GATE

We describe the first developments towards a Monte Carlo X-ray phase contrast imaging simulator for the medical imaging and radiotherapy simulation software GATE. Phase contrast imaging is an imaging modality taking advantage of the phase shift of X-rays. This modality produces images with a higher sensitivity than conventional, attenuation based imaging. As the first developments towards Monte Carlo phase contrast simulation, we implemented a Monte Carlo process for the refraction and total reflection of X-rays, as well as an analytical wave optics approach for generating Fresnel diffraction patterns. The implementation is validated against data acquired using a laboratory X-ray tomography system. The overall agreement between the simulations and the data is encouraging, which motivates further development of Monte Carlo based simulation of X-ray phase contrast imaging. These developments have been released in GATE version 8.2.

Visibility simulation of realistic grating interferometers including grating geometries and energy spectra

The performance of X-ray and neutron grating interferometers is characterised by their visibility, which is a measure for the maximum achievable contrast. In this study we show how the real grating geometry in a grating interferometer with three gratings impacts the interference and self projection that leads to visibility in the first place. We quantify the individual contributions of wavelength distributions and grating shapes in terms of visibility reduction by determining the absolute as well as relative effect of each contribution. The understanding of the impact of changed geometry and wavelength distributions on the interference of neutrons/X-rays allows us to present the first fully quantitative model of a grating interferometer setup. We demonstrate the capability of the simulation framework by building a model of the neutron grating interferometer at the ICON beamline and directly comparing simulated and measured visibility values. The general nature of the model makes it possible to extend it to any given grating interferometer for both X-rays and neutrons.

A generalized quantitative interpretation of dark-field contrast for highly concentrated microsphere suspensions

In X-ray grating interferometry, dark-field contrast arises due to partial extinction of the detected interference fringes. This is also called visibility reduction and is attributed to small-angle scattering from unresolved structures in the imaged object. In recent years, analytical quantitative frameworks of dark-field contrast have been developed for highly diluted monodisperse microsphere suspensions with maximum 6% volume fraction. These frameworks assume that scattering particles are separated by large enough distances, which make any interparticle scattering interference negligible. In this paper, we start from the small-angle scattering intensity equation and, by linking Fourier and real-space, we introduce the structure factor and thus extend the analytical and experimental quantitative interpretation of dark-field contrast, for a range of suspensions with volume fractions reaching 40%. The structure factor accounts for interparticle scattering interference. Without introducing any additional fitting parameters, we successfully predict the experimental values measured at the TOMCAT beamline, Swiss Light Source. Finally, we apply this theoretical framework to an experiment probing a range of system correlation lengths by acquiring dark-field images at different energies. This proposed method has the potential to be applied in single-shot-mode using a polychromatic X-ray tube setup and a single-photon-counting energy-resolving detector.

Effects of unresolvable edges in grating-based X-ray differential phase imaging

We investigated effects of unresolvable sharp edges on images obtained in a grating-based X-ray differential phase imaging technique. Results of numerical calculations for monochromatic X-rays show that an unresolvable sharp edge generates not only differential-phase contrast but also visibility contrast. The latter shows that the visibility contrast has another major origin other than ultra-small-angle X-ray scattering (USAXS) from randomly distributed unresolvable microstructures, which has been considered the main origin for the contrast. The effects were experimentally confirmed using a synchrotron X-ray source.

Inclusion of coherence in Monte Carlo models for simulation of x-ray phase contrast imaging

Interest in phase contrast imaging methods based on electromagnetic wave coherence has increased significantly recently, particularly at X-ray energies. This is giving rise to a demand for effective simulation methods. Coherent imaging approaches are usually based on wave optics, which require significant computational resources, particularly for producing 2D images. Monte Carlo (MC) methods, used to track individual particles/photons for particle physics, are not considered appropriate for describing coherence effects. Previous preliminary work has evaluated the possibility of incorporating coherence in Monte Carlo codes. However, in this paper, we present the implementation of refraction in a model that is based on time of flight calculations and the Huygens-Fresnel principle, which allow reproducing the formation of phase contrast images in partially and fully coherent experimental conditions. The model is implemented in the FLUKA Monte Carlo code and X-ray phase contrast imaging simulations are compared with experiments and wave optics calculations.