Dirk Bequé

Multi-source inverse geometry CT: a new system concept for x-ray computed tomography

Third-generation CT architectures are approaching fundamental limits. Spatial resolution is limited by the focal spot size and the detector cell size. Temporal resolution is limited by mechanical constraints on gantry rotation speed, and alternative geometries such as electron-beam CT and two-tube-two-detector CT come with severe tradeoffs in terms of image quality, dose-efficiency and complexity. Image noise is fundamentally linked to patient dose, and dose-efficiency is limited by finite detector efficiency and by limited spatio-temporal control over the X-ray flux. Finally, volumetric coverage is limited by detector size, scattered radiation, conebeam artifacts, Heel effect, and helical over-scan. We propose a new concept, multi-source inverse geometry CT, which allows CT to break through several of the above limitations. The proposed architecture has several advantages compared to third-generation CT: the detector is small and can have a high detection efficiency, the optical spot size is more consistent throughout the field-of-view, scatter is minimized even when eliminating the anti-scatter grid, the X-ray flux from each source can be modulated independently to achieve an optimal noise-dose tradeoff, and the geometry offers unlimited coverage without cone-beam artifacts. In this work we demonstrate the advantages of multi-source inverse geometry CT using computer simulations.

A Fourier-domain algorithm for total-variation regularized phase retrieval in differential X-ray phase contrast imaging

Phase retrieval in differential X-ray phase contrast imaging involves a one dimensional integration step. In the presence of noise, standard integration methods result in image blurring and streak artifacts. This work proposes a regularized integration method which takes the availability of two dimensional data as well as the integration-specific frequency-dependent noise amplification into account. In more detail, a Fourier-domain algorithm is developed comprising a frequency-dependent minimization of the total variation orthogonal to the direction of integration. For both simulated and experimental data, the novel method yielded strong artefact reduction without increased blurring superior to the results obtained by standard integration methods or regularization techniques in the image domain.

Energy dispersive X-ray diffraction spectral resolution considerations for security screening applications

Energy dispersive X-ray diffraction (EDXRD) is a very effective method for explosive and narcotic threat detection in baggage screening. The XRD profiles arise from the molecular interference when X-rays are coherently scattered by a substance. The accurate identification of the target material depends on the ability to detect and resolve the peaks present in the coherent scatter profiles. A high-energy resolution High Purity Germanium (HPGe) detector is therefore generally used in such type of systems. To evaluate the suitability of cost-effective room-temperature semiconductor detectors for next-generation baggage screening systems, an assessment of the minimal requirements for the system resolution is required. In this study a hybrid Monte Carlo code has been modified to account for the molecular interference function that gives rise to the coherent scatter signature. A model for a realistic response function for Cadmium Zinc Telluride (CZT) detectors is then used to convolve the spectral output. This simulation tool is then used to assess the system design features and their influence on spectral resolution.

Low-dose, phase-contrast mammography with high signal-to-noise ratio

Differential phase-contrast X-ray imaging using a Talbot-Lau interferometer has recently shown promising results for applications in medical imaging. However, reducing the applied radiation dose remains a major challenge. In this study, we consider the realization of a Talbot-Lau interferometer in a high Talbot order to increase the signal-to-noise ratio for low-dose applications. The quantitative performance of π and π/2 systems at high Talbot orders is analyzed through simulations, and the design energy and X-ray spectrum are optimized for mammography. It is found that operation even at very high Talbot orders is feasible and beneficial for image quality. As long as the X-ray spectrum is matched to the visibility spectrum, the SNR continuously increases with the Talbot order for π-systems. We find that the optimal X-ray spectra and design energies are almost independent of the Talbot order and that the overall imaging performance is robust against small variations in these parameters. Discontinuous spectra, such as that from molybdenum, are less robust because the characteristic lines may coincide with minima in the visibility spectra; however, they may offer slightly better performance. We verify this hypothesis by realizing a prototype system with a mean fringe visibility of above 40% at the seventh Talbot order. With this prototype, a proof-of-principle measurement of a freshly dissected breast at reasonable compression to 4 cm is conducted with a mean glandular dose of only 3 mGy but with a high SNR.

Iterative reconstruction for multi-source inverse geometry CT: a feasibility study

In a 3rd generation CT system, a single source projects the entire field of view (FOV) onto a large detector opposite to the source. In multi-source inverse geometry CT imaging, a multitude of sources sequentially project complementary parts of the FOV on a much smaller detector. These sources may be distributed in both the trans-axial and axial directions and jointly cover the entire FOV. Multi-source CT has several important advantages, including large axial coverage, improved dose-efficiency, and improved spatial resolution. One of the challenges of this concept is to ensure that no artifacts emerge in the reconstructed images where the sampling switches from one source to the next. This work studies iterative reconstruction for multi-source imaging and focuses on the appearance of such artifacts. For that purpose, phantom data are simulated using a realistic multi-source CT geometry, iteratively reconstructed and inspected for artifact content. More realistic experiments using rebinned clinical datasets (emulating a multi-source CT system) have also been performed. The results confirm the feasibility of artifact-free multi-source CT imaging in both full-scan and half-scan situations.

Fast one-dimensional wave-front propagation for x-ray differential phase-contrast imaging

Numerical wave-optical simulations of X-ray differential phase-contrast imaging using grating interferometry require the oversampling of gratings and object structures in the range of few micrometers. Consequently, fields of view of few millimeters already use large amounts of a computers main memory to store the propagating wave front, limiting the scope of the investigations to only small-scale problems. In this study, we apply an approximation to the Fresnel-Kirchhoff diffraction theory to overcome these restrictions by dividing the two-dimensional wave front up into 1D lines, which are processed separately. The approach enables simulations with samples of clinically relevant dimensions by significantly reducing the memory footprint and the execution time and, thus, allows the qualitative comparison of different setup configurations. We analyze advantages as well as limitations and present the simulation of a virtual mammography phantom of several centimeters of size.

Coronary artery motion estimation and compensation: A feasibility study

High temporal resolution and high spatial resolution are required to image the coronary arteries without motion artifacts. Several approaches have been pursued to achieve better temporal resolution including faster rotational speeds, and dual tube systems. In this paper, we present an alternative approach using motion estimation and compensation. The results demonstrate that the proposed methods can significantly reduce motion artifacts in coronary artery imaging.

Efficient algorithm for modeling keel-edge pinhole response

We describe an efficient method for computing the system response (sensitivity and point spread function) of a conventional gamma camera with multiple keel-edge pinhole collimators, taking into account geometric acceptance effects at the collimator as well as detector response. An efficient way of storing the resulting system matrix permits us to achieve very fast iterative reconstruction with reasonable memory requirements.

Geometric model, control & calibration of universal benchtop CT system

A universal benchtop CT system has recently been developed and installed at GE Global Research. The system consists of a source, detector and phantom table that are mounted on a set of high accuracy motion stages. These stages permit positioning and orienting the source, detector and phantom with respect to each other. The benchtop system is thereby able to simulate any realistic clinical CT system geometry. This work provides a geometrical model of this benchtop system and illustrates how the stages need to be controlled in order to simulate a particular CT system. This illustrates at the same time the geometrical capabilities of the system. The method and first results of the geometrical calibration of the system are also shown.

Fourier domain image fusion for differential X-ray phase-contrast breast imaging

X-Ray Phase-Contrast (XPC) imaging is a novel technology with a great potential for applications in clinical practice, with breast imaging being of special interest. This work introduces an intuitive methodology to combine and visualize relevant diagnostic features, present in the X-ray attenuation, phase shift and scattering information retrieved in XPC imaging, using a Fourier domain fusion algorithm. The method allows to present complementary information from the three acquired signals in one single image, minimizing the noise component and maintaining visual similarity to a conventional X-ray image, but with noticeable enhancement in diagnostic features, details and resolution. Radiologists experienced in mammography applied the image fusion method to XPC measurements of mastectomy samples and evaluated the feature content of each input and the fused image. This assessment validated that the combination of all the relevant diagnostic features, contained in the XPC images, was present in the fused image as well.