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.

Efficient iterative reconstruction with beam shape compensation for THz computed tomography.

Terahertz (THz) computed tomography is an emerging nondestructive and non-ionizing imaging method. Most THz reconstruction methods rely on the Radon transform, originating from x-ray imaging, in which x rays propagate in straight lines. However, a THz beam has a finite width, and ignoring its shape results in blurred reconstructed images. Moreover, accounting for the THz beam model in a straightforward way in an iterative reconstruction method results in extreme demands in memory and in slow convergence. In this paper, we propose an efficient iterative reconstruction that incorporates the THz beam shape, while avoiding the above disadvantages. Both simulation and real experiments show that our approach results in improved resolution recovery in the reconstructed image. Furthermore, we propose a suitable preconditioner to improve the convergence speed of our reconstruction.

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.

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.

poly-DART: A discrete algebraic reconstruction technique for polychromatic X-ray CT.

The discrete algebraic reconstruction technique (DART) is a tomographic method to reconstruct images from X-ray projections in which prior knowledge on the number of object materials is exploited. In monochromatic X-ray CT (e.g., synchrotron), DART has been shown to lead to high-quality reconstructions, even with a low number of projections or a limited scanning view. However, most X-ray sources are polychromatic, leading to beam hardening effects, which significantly degrade the performance of DART. In this work, we propose a new discrete tomography algorithm, poly-DART, that exploits sparsity in the attenuation values using DART and simultaneously accounts for the polychromatic nature of the X-ray source. The results show that poly-DART leads to a vastly improved segmentation on polychromatic data obtained from Monte Carlo simulations as well as on experimental data, compared to DART.

Parametric Reconstruction of Glass Fiber-reinforced Polymer Composites from X-ray Projection Data-A Simulation Study.

We present a new approach to estimate geometry parameters of glass fibers in glass fiber-reinforced polymers from simulated X-ray micro-computed tomography scans. Traditionally, these parameters are estimated using a multi-step procedure including image reconstruction, pre-processing, segmentation and analysis of features of interest. Each step in this chain introduces errors that propagate through the pipeline and impair the accuracy of the estimated parameters. In the approach presented in this paper, we reconstruct volumes from a low number of projection angles using an iterative reconstruction technique and then estimate position, direction and length of the contained fibers incorporating a priori knowledge about their shape, modeled as a geometric representation, which is then optimized. Using simulation experiments, we show that our method can estimate those representations even in presence of noisy data and only very few projection angles available.

MoVIT: a tomographic reconstruction framework for 4D-CT.

4D computed tomography (4D-CT) aims to visualise the temporal dynamics of a 3D sample with a sufficiently high temporal and spatial resolution. Successive time frames are typically obtained by sequential scanning, followed by independent reconstruction of each 3D dataset. Such an approach requires a large number of projections for each scan to obtain images with sufficient quality (in terms of artefacts and SNR). Hence, there is a clear trade-off between the rotation speed of the gantry (i.e. time resolution) and the quality of the reconstructed images. In this paper, the MotionVector-based Iterative Technique (MoVIT) is introduced which reconstructs a particular time frame by including the projections of neighbouring time frames as well. It is shown that such a strategy improves the trade-off between the rotation speed and the SNR. The framework is tested on both numerical simulations and on 4D X-ray CT datasets of polyurethane foam under compression. Results show that reconstructions obtained with MoVIT have a significantly higher SNR compared to the SNR of conventional 4D reconstructions.

Fast and flexible X-ray tomography using the ASTRA toolbox.

Object reconstruction from a series of projection images, such as in computed tomography (CT), is a popular tool in many different application fields. Existing commercial software typically provides sufficiently accurate and convenient-to-use reconstruction tools to the end-user. However, in applications where a non-standard acquisition protocol is used, or where advanced reconstruction methods are required, the standard software tools often are incapable of computing accurate reconstruction images. This article introduces the ASTRA Toolbox. Aimed at researchers across multiple tomographic application fields, the ASTRA Toolbox provides a highly efficient and highly flexible open source set of tools for tomographic projection and reconstruction. The main features of the ASTRA Toolbox are discussed and several use cases are presented.

Measuring Lattice Strain in Three Dimensions through Electron Microscopy.

The three-dimensional (3D) atomic structure of nanomaterials, including strain, is crucial to understand their properties. Here, we investigate lattice strain in Au nanodecahedra using electron tomography. Although different electron tomography techniques enabled 3D characterizations of nanostructures at the atomic level, a reliable determination of lattice strain is not straightforward. We therefore propose a novel model-based approach from which atomic coordinates are measured. Our findings demonstrate the importance of investigating lattice strain in 3D.

Dynamic intensity normalization using eigen flat fields in X-ray imaging.

In X-ray imaging, it is common practice to normalize the acquired projection data with averaged flat fields taken prior to the scan. Unfortunately, due to source instabilities, vibrating beamline components such as the monochromator, time varying detector properties, or other confounding factors, flat fields are often far from stationary, resulting in significant systematic errors in intensity normalization. In this work, a simple and efficient method is proposed to account for dynamically varying flat fields. Through principal component analysis of a set of flat fields, eigen flat fields are computed. A linear combination of the most important eigen flat fields is then used to individually normalize each X-ray projection. Experiments show that the proposed dynamic flat field correction leads to a substantial reduction of systematic errors in projection intensity normalization compared to conventional flat field correction.

Fiber Orientation Estimation from X-ray Dark Field Images of Fiber Reinforced Polymers Using Constrained Spherical Deconvolution.

The properties of fiber reinforced polymers are strongly related to the length and orientation of the fibers within the polymer matrix, the latter of which can be studied using X-ray computed tomography (XCT). Unfortunately, resolving individual fibers is challenging because they are small compared to the XCT voxel resolution and because of the low attenuation contrast between the fibers and the surrounding resin. To alleviate both problems, anisotropic dark field tomography via grating based interferometry (GBI) has been proposed. Here, the fiber orientations are extracted by applying a Funk-Radon transform (FRT) to the local scatter function. However, the FRT suffers from a low angular resolution, which complicates estimating fiber orientations for small fiber crossing angles. We propose constrained spherical deconvolution (CSD) as an alternative to the FRT to resolve fiber orientations. Instead of GBI, edge illumination phase contrast imaging is used because estimating fiber orientations with this technique has not yet been explored. Dark field images are generated by a Monte Carlo simulation framework. It is shown that the FRT cannot estimate the fiber orientation accurately for crossing angles smaller than 70∘, while CSD performs well down to a crossing angle of 50∘. In general, CSD outperforms the FRT in estimating fiber orientations.

Mesh modeling of system geometry and anatomy phantoms for realistic GATE simulations and their inclusion in SPECT reconstruction.

Objective.Monte-Carlo simulation studies have been essential for advancing various developments in single photon emission computed tomography (SPECT) imaging, such as system design and accurate image reconstruction. Among the simulation software available, Geant4 application for tomographic emission (GATE) is one of the most used simulation toolkits in nuclear medicine, which allows building systems and attenuation phantom geometries based on the combination of idealized volumes. However, these idealized volumes are inadequate for modeling free-form shape components of such geometries. Recent GATE versions alleviate these major limitations by allowing users to import triangulated surface meshes.Approach.In this study, we describe our mesh-based simulations of a next-generation multi-pinhole SPECT system dedicated to clinical brain imaging, called AdaptiSPECT-C. To simulate realistic imaging data, we incorporated in our simulation the XCAT phantom, which provides an advanced anatomical description of the human body. An additional challenge with the AdaptiSPECT-C geometry is that the default voxelized XCAT attenuation phantom was not usable in our simulation due to intersection of objects of dissimilar materials caused by overlap of the air containing regions of the XCAT beyond the surface of the phantom and the components of the imaging system.Main results.We validated our mesh-based modeling against the one constructed by idealized volumes for a simplified single vertex configuration of AdaptiSPECT-C through simulated projection data of123I-activity distributions. We resolved the overlap conflict by creating and incorporating a mesh-based attenuation phantom following a volume hierarchy. We then evaluated our reconstructions with attenuation and scatter correction for projections obtained from simulation consisting of mesh-based modeling of the system and the attenuation phantom for brain imaging. Our approach demonstrated similar performance as the reference scheme simulated in air for uniform and clinical-like123I-IMP brain perfusion source distributions.Significance.This work enables the simulation of complex SPECT acquisitions and reconstructions for emulating realistic imaging data close to those of actual patients.

Automatic landmark detection and mapping for 2D/3D registration with BoneNet.

The 3D musculoskeletal motion of animals is of interest for various biological studies and can be derived from X-ray fluoroscopy acquisitions by means of image matching or manual landmark annotation and mapping. While the image matching method requires a robust similarity measure (intensity-based) or an expensive computation (tomographic reconstruction-based), the manual annotation method depends on the experience of operators. In this paper, we tackle these challenges by a strategic approach that consists of two building blocks: an automated 3D landmark extraction technique and a deep neural network for 2D landmarks detection. For 3D landmark extraction, we propose a technique based on the shortest voxel coordinate variance to extract the 3D landmarks from the 3D tomographic reconstruction of an object. For 2D landmark detection, we propose a customized ResNet18-based neural network, BoneNet, to automatically detect geometrical landmarks on X-ray fluoroscopy images. With a deeper network architecture in comparison to the original ResNet18 model, BoneNet can extract and propagate feature vectors for accurate 2D landmark inference. The 3D poses of the animal are then reconstructed by aligning the extracted 2D landmarks from X-ray radiographs and the corresponding 3D landmarks in a 3D object reference model. Our proposed method is validated on X-ray images, simulated from a real piglet hindlimb 3D computed tomography scan and does not require manual annotation of landmark positions. The simulation results show that BoneNet is able to accurately detect the 2D landmarks in simulated, noisy 2D X-ray images, resulting in promising rigid and articulated parameter estimations.

Adjoint image warping using multivariate splines with application to four-dimensional computed tomography.

PURPOSE: Adjoint image warping is an important tool to solve image reconstruction problems that warp the unknown image in the forward model. This includes four-dimensional computed tomography (4D-CT) models in which images are compared against recorded projection images of various time frames using image warping as a model of the motion. The inversion of these models requires the adjoint of image warping, which up to now has been substituted by approximations. We introduce an efficient implementation of the exact adjoints of multivariate spline based image warping, and compare it against previously used alternatives. METHODS: Using symbolic computer algebra, we computed a list of 64 polynomials that allow us to compute a matrix representation of trivariate cubic image warping. By combining an on-the-fly computation of this matrix with a parallelized implementation of columnwise matrix multiplication, we obtained an efficient, low memory implementation of the adjoint action of 3D cubic image warping. We used this operator in the solution of a previously proposed 4D-CT reconstruction model in which the image of a single subscan was compared against projection data of multiple subscans by warping and then projecting the image. We compared the properties of our exact adjoint with those of approximate adjoints by warping along inverted motion. RESULTS: Our method requires halve the memory to store motion between subscans, compared to methods that need to compute and store an approximate inverse of the motion. It also avoids the computation time to invert the motion and the tunable parameter of the number of iterations used to perform this inversion. Yet, a similar and often better reconstruction quality was obtained in comparison with these more expensive methods, especially when the motion is large. When compared against a simpler method that is similar to ours in computational demands, our method achieves a higher reconstruction quality in general. CONCLUSIONS: Our implementation of the exact adjoint of cubic image warping improves efficiency and provides accurate reconstructions.

Geometry Calibration of a Modular Stereo Cone-Beam X-ray CT System.

Compared to single source systems, stereo X-ray CT systems allow acquiring projection data within a reduced amount of time, for an extended field-of-view, or for dual X-ray energies. To exploit the benefit of a dual X-ray system, its acquisition geometry needs to be calibrated. Unfortunately, in modular stereo X-ray CT setups, geometry misalignment occurs each time the setup is changed, which calls for an efficient calibration procedure. Although many studies have been dealing with geometry calibration of an X-ray CT system, little research targets the calibration of a dual cone-beam X-ray CT system. In this work, we present a phantom-based calibration procedure to accurately estimate the geometry of a stereo cone-beam X-ray CT system. With simulated as well as real experiments, it is shown that the calibration procedure can be used to accurately estimate the geometry of a modular stereo X-ray CT system thereby reducing the misalignment artifacts in the reconstruction volumes.

Extended imaging volume in cone-beam x-ray tomography using the weighted simultaneous iterative reconstruction technique.

An issue in computerized x-ray tomography is the limited size of available detectors relative to objects of interest. A solution was provided in the past two decades by positioning the detector in a lateral offset position, increasing the effective field of view (FOV) and thus the diameter of the reconstructed volume. However, this introduced artifacts in the obtained reconstructions, caused by projection truncation and data redundancy. These issues can be addressed by incorporating an additional data weighting step in the reconstruction algorithms, known as redundancy weighting. In this work, we present an implementation of redundancy weighting in the widely-used simultaneous iterative reconstruction technique (SIRT), yielding the weighted SIRT (W-SIRT) method. The new technique is validated using geometric phantoms and a rabbit specimen, by performing both simulation studies as well as physical experiments. The experiments are carried out in a highly flexible stereoscopic x-ray system equipped with x-ray image intensifiers (XRIIs). The simulations showed that higher values of contrast-to-noise ratio could be obtained using the W-SIRT approach as compared to a weighted implementation of the simultaneous algebraic reconstruction technique (SART). The convergence rate of the W-SIRT was accelerated by including a relaxation parameter in the W-SIRT algorithm, creating the aW-SIRT algorithm. This allowed to obtain the same results as the W-SIRT algorithm, but at half the number of iterations, yielding a much shorter computation time. The aW-SIRT algorithm has proven to perform well for both large as well as small regions of overlap, outperforming the pre-convolutional Feldkamp-David-Kress algorithm for small overlap regions (or large detector offsets). The experiments confirmed the results of the simulations. Using the aW-SIRT algorithm, the effective FOV was increased by >75%, only limited by experimental constraints. Although an XRII is used in this work, the method readily applies to flat-panel detectors as well.

Fast simulation of yttrium-90 bremsstrahlung photons with GATE.

PURPOSE: Multiple investigators have recently reported the use of yttrium-90 (90Y) bremsstrahlung single photon emission computed tomography (SPECT) imaging for the dosimetry of targeted radionuclide therapies. Because Monte Carlo (MC) simulations are useful for studying SPECT imaging, this study investigates the MC simulation of 90Y bremsstrahlung photons in SPECT. To overcome the computationally expensive simulation of electrons, the authors propose a fast way to simulate the emission of 90Y bremsstrahlung photons based on prerecorded bremsstrahlung photon probability density functions (PDFs). METHODS: The accuracy of bremsstrahlung photon simulation is evaluated in two steps. First, the validity of the fast bremsstrahlung photon generator is checked. To that end, fast and analog simulations of photons emitted from a 90Y point source in a water phantom are compared. The same setup is then used to verify the accuracy of the bremsstrahlung photon simulations, comparing the results obtained with PDFs generated from both simulated and measured data to measurements. In both cases, the energy spectra and point spread functions of the photons detected in a scintillation camera are used. RESULTS: Results show that the fast simulation method is responsible for a 5% overestimation of the low-energy fluence (below 75 keV) of the bremsstrahlung photons detected using a scintillation camera. The spatial distribution of the detected photons is, however, accurately reproduced with the fast method and a computational acceleration of approximately 17-fold is achieved. When measured PDFs are used in the simulations, the simulated energy spectrum of photons emitted from a point source of 90Y in a water phantom and detected in a scintillation camera closely approximates the measured spectrum. The PSF of the photons imaged in the 50-300 keV energy window is also accurately estimated with a 12.4% underestimation of the full width at half maximum and 4.5% underestimation of the full width at tenth maximum. CONCLUSIONS: Despite its limited accuracy, the fast bremsstrahlung photon generator is well suited for the simulation of bremsstrahlung photons emitted in large homogeneous organs, such as the liver, and detected in a scintillation camera. The computational acceleration makes it very useful for future investigations of 90Y bremsstrahlung SPECT imaging.

Fast and memory-efficient Monte Carlo-based image reconstruction for whole-body PET.

PURPOSE: Several studies have shown the benefit of an accurate system modeling using Monte Carlo techniques. For state-of-the-art whole-body positron emission tomography (PET) scanners, Monte Carlo-based image reconstruction is associated with a significant computational cost to calculate the system matrix as well as a large memory capacity to store it. In this article, the authors present a simulation-reconstruction framework to solve these problems on the Philips Gemini GS PET scanner. METHODS: A fast, realistic system matrix simulation module was developed using egs_pet, which is an efficient PET simulation code based on EGSnrc. The generated system matrix was then used in a rotator-based ordered subset expectation maximization (OS-EM) algorithm, which exploits the rotational symmetry of a cylindrical PET scanner. The system matrix was further compressed by using sparse storage techniques. RESULTS: The system matrix simulation took five days on 50 cores of Xeon 2.66 GHz, resulting in a system matrix of 2.01 GB. The entire system matrix could be stored in the main memory of a standard personal computer. The image quality in terms of contrast-noise trade-offs was considerably improved compared to a standard OS-EM algorithm. The image quality was also compared to the clinical software on the scanner using routine parameter settings. The contrast recovery coefficient of small hot spheres and cold spheres was significantly improved. CONCLUSIONS: The results indicated that the proposed framework could be used for this PET scanner with improved image quality. This method could also be applied to other state-of-the-art whole-body PET scanners and preclinical PET scanners with a similar shape.

A Machine Learning Approach to Growth Direction Finding for Automated Planting of Bulbous Plants.

In agricultural robotics, a unique challenge exists in the automated planting of bulbous plants: the estimation of the bulb's growth direction. To date, no existing work addresses this challenge. Therefore, we propose the first robotic vision framework for the estimation of a plant bulb's growth direction. The framework takes as input three x-ray images of the bulb and extracts shape, edge, and texture features from each image. These features are then fed into a machine learning regression algorithm in order to predict the 2D projection of the bulb's growth direction. Using the x-ray system's geometry, these 2D estimates are then mapped to the 3D world coordinate space, where a filtering on the estimate's variance is used to determine whether the estimate is reliable. We applied our algorithm on 27,200 x-ray simulations from T. Apeldoorn bulbs on a standard desktop workstation. Results indicate that our machine learning framework is fast enough to meet industry standards (<0.1 seconds per bulb) while providing acceptable accuracy (e.g. error < 30° in 98.40% of cases using an artificial 3-layer neural network). The high success rates of the proposed framework indicate that it is worthwhile to proceed with the development and testing of a physical prototype of a robotic bulb planting system.