Joseph Zambelli

Radiation dose efficiency comparison between differential phase contrast CT and conventional absorption CT

PURPOSE The superior radiation dose efficiency of a newly implemented differential phase contrast CT imaging method compared to the conventional absorption CT method is demonstrated. METHODS A differential phase contrast CT imaging method has recently been implemented using conventional x-ray sources with a grating interferometer consisting of three gratings. This approach offers the possibility of simultaneous reconstruction of both attenuation contrast and phase contrast images from a single acquisition. This enables a direct comparison of radiation dose efficiency of both types of reconstructed images under identical conditions. Radiation dose efficiency was studied by measuring the change in contrast-to-noise ratio (CNR) with different exposure levels. A physical phantom of 28.5 mm diameter was constructed and used for measurement of CNR in both the absorption and phase contrast CT images, which were reconstructed from the same data set. RESULTS For three of the four materials studied, at any given exposure level, the CNR of the differential phase contrast CT images was superior to that of the corresponding absorption contrast CT images. The most dramatic improvement was noted in the contrast between PMMA and water, where the CNR was improved by a factor of approximately 5.5 in the differential phase contrast CT images. Additionally, the CNR of phase contrast CT is empirically shown to have the same square root dependence on exposure, as is the case for absorption contrast CT. CONCLUSIONS The differential phase contrast CT method provided higher CNR than conventional absorption CT at equivalent dose levels for most of the materials studied, and so may enable achievement of the same object visibility as conventional absorption CT methods at a lower exposure level.

Small-angle scattering computed tomography (SAS-CT) using a Talbot-Lau interferometer and a rotating anode x-ray tube: theory and experiments

X-ray differential phase contrast imaging methods, including projection imaging and the corresponding computed tomography (CT), have been implemented using a Talbot interferometer and either a synchrotron beam line or a low brilliance x-ray source generated by a stationary-anode x-ray tube. From small-angle scattering events which occur as an x-ray propagates through a medium, a signal intensity loss can be recorded and analyzed for an understanding of the micro-structures in an image object. This has been demonstrated using a Talbot-Lau interferometer and a stationary-anode x-ray tube. In this paper, theoretical principles and an experimental implementation of the corresponding CT imaging method are presented. First, a line integral is derived from analyzing the cross section of the small-angle scattering events. This method is referred to as small-angle scattering computed tomography (SAS-CT). Next, a Talbot-Lau interferometer and a rotating-anode x-ray tube were used to implement SAS-CT. A physical phantom and human breast tissue sample were used to demonstrate the reconstructed SAS-CT image volumes.

Quantitative imaging of electron density and effective atomic number using phase contrast CT.

Compared to single energy CT, which only provides information for x-ray linear attenuation coefficients, dual-energy CT is able to obtain both the electron density and effective atomic number for different materials in a quantitative way. In this study, as an alternative to dual-energy CT, a novel quantitative imaging method based on phase contrast CT is presented. Rather than requiring two projection data sets with different x-ray energy spectra, diffraction-grating-based phase contrast CT is capable of reconstructing images of both linear attenuation and refractive index decrement from the same projection data using a single x-ray energy spectra. From the two images, quantitative information of both the electron density and effective atomic number can be extracted. Two physical phantoms were constructed and used to validate the presented method. Experimental results demonstrate that (1) electron density can be accurately determined from refractive index decrement through a linear relationship, and (2) the effective atomic number can be explicitly derived from the ratio of the linear attenuation to refractive index decrement using a power function plus a constant. The presented method will provide insight into the technique of material separation and find its use in medical and industrial applications.

Multicontrast x-ray computed tomography imaging using Talbot-Lau interferometry without phase stepping

PURPOSE The purpose of this work is to demonstrate that multicontrast computed tomography (CT) imaging can be performed using a Talbot-Lau interferometer without phase stepping, thus allowing for an acquisition scheme like that used for standard absorption CT. METHODS Rather than using phase stepping to extract refraction, small-angle scattering (SAS), and absorption signals, the two gratings of a Talbot-Lau interferometer were rotated slightly to generate a moiré pattern on the detector. A Fourier analysis of the moiré pattern was performed to obtain separate projection images of each of the three contrast signals, all from the same single-shot of x-ray exposure. After the signals were extracted from the detector data for all view angles, image reconstruction was performed to obtain absorption, refraction, and SAS CT images. A physical phantom was scanned to validate the proposed data acquisition method. The results were compared with a phantom scan using the standard phase stepping approach. RESULTS The reconstruction of each contrast mechanism produced the expected results. Signal levels and contrasts match those obtained using the phase stepping technique. CONCLUSIONS Absorption, refraction, and SAS CT imaging can be achieved using the Talbot-Lau interferometer without the additional overhead of long scan time and phase stepping.

Time-Resolved Interventional Cardiac C-arm Cone-Beam CT: An Application of the PICCS Algorithm

Time-resolved cardiac imaging is particularly interesting in the interventional setting since it would provide both image guidance for accurate procedural planning and cardiac functional evaluations directly in the operating room. Imaging the heart in vivo using a slowly rotating C-arm system is extremely challenging due to the limitations of the data acquisition system and the high temporal resolution required to avoid motion artifacts. In this paper, a data acquisition scheme and an image reconstruction method are proposed to achieve time-resolved cardiac cone-beam computed tomography imaging with isotropic spatial resolution and high temporal resolution using a slowly rotating C-arm system. The data are acquired within 14 s using a single gantry rotation with a short scan angular range. The enabling image reconstruction method is the prior image constrained compressed sensing (PICCS) algorithm. The prior image is reconstructed from data acquired over all cardiac phases. Each cardiac phase is then reconstructed from the retrospectively gated cardiac data using the PICCS algorithm. To validate the method, several studies were performed. Both numerical simulations using a hybrid motion phantom with static background anatomy as well as physical phantom studies have been used to demonstrate that the proposed method enables accurate reconstruction of image objects with a high isotropic spatial resolution. A canine animal model scanned in vivo was used to further validate the method.

Scaling law for noise variance and spatial resolution in differential phase contrast computed tomography

PURPOSE The noise variance versus spatial resolution relationship in differential phase contrast (DPC) projection imaging and computed tomography (CT) are derived and compared to conventional absorption-based x-ray projection imaging and CT. METHODS The scaling law for DPC-CT is theoretically derived and subsequently validated with phantom results from an experimental Talbot-Lau interferometer system. RESULTS For the DPC imaging method, the noise variance in the differential projection images follows the same inverse-square law with spatial resolution as in conventional absorption-based x-ray imaging projections. However, both in theory and experimental results, in DPC-CT the noise variance scales with spatial resolution following an inverse linear relationship with fixed slice thickness. CONCLUSIONS The scaling law in DPC-CT implies a lesser noise, and therefore dose, penalty for moving to higher spatial resolutions when compared to conventional absorption-based CT in order to maintain the same contrast-to-noise ratio.

Design and development of C-arm based cone-beam CT for image-guided interventions : Initial results

X-ray cone-beam computed tomography (CBCT) is of importance in image-guided intervention (IGI) and image-guided radiation therapy (IGRT). In this paper, we present a cone-beam CT data acquisition system using a GE INNOVA 4100 (GE Healthcare Technologies, Waukesha, Wisconsin) clinical system. This new cone-beam data acquisition mode was developed for research purposes without interfering with any clinical function of the system. It provides us a basic imaging pipeline for more advanced cone-beam data acquisition methods. It also provides us a platform to study and overcome the limiting factors such as cone-beam artifacts and limiting low contrast resolution in current C-arm based cone-beam CT systems. A geometrical calibration method was developed to experimentally determine parameters of the scanning geometry to correct the image reconstruction for geometric non-idealities. Extensive phantom studies and some small animal studies have been conducted to evaluate the performance of our cone-beam CT data acquisition system.

Interior tomography in x-ray differential phase contrast CT imaging.

Differential phase contrast computed tomography (DPC-CT) is an x-ray imaging method that uses the wave properties of imaging photons as the contrast mechanism. It has been demonstrated that DPC images can be obtained using a conventional x-ray tube and a Talbot-Lau-type interferometer. Due to the limited size of the gratings, current data acquisition systems only offer a limited field of view, and thus are prone to data truncation. As a result, the reconstructed DPC-CT image may suffer from image artifacts and increased inaccuracy in the reconstructed image values. In this paper, we demonstrate that a small region of interest (ROI) within a large object can be accurately and stably reconstructed using fully truncated projection datasets provided that a priori information on electron density is known for a small region inside the ROI. The method reconstructs an image iteratively to satisfy a group of physical conditions by using a projection onto convex set (POCS) approach. In this work, this POCS algorithm is validated using both numerical simulations and physical phantom experimental data. In both cases, the root mean square error is reduced by an order of magnitude with respect to the truncated analytic reconstructions. Truncation artifacts observed in the latter reconstructions are eliminated using the POCS algorithm.

Fundamental relationship between the noise properties of grating-based differential phase contrast CT and absorption CT: Theoretical framework using a cascaded system model and experimental validation

PURPOSE Using a grating interferometer, a conventional x-ray cone beam computed tomography (CT) data acquisition system can be used to simultaneously generate both conventional absorption CT (ACT) and differential phase contrast CT (DPC-CT) images from a single data acquisition. Since the two CT images were extracted from the same set of x-ray projections, it is expected that intrinsic relationships exist between the noise properties of the two contrast mechanisms. The purpose of this paper is to investigate these relationships. METHODS First, a theoretical framework was developed using a cascaded system model analysis to investigate the relationship between the noise power spectra (NPS) of DPC-CT and ACT. Based on the derived analytical expressions of the NPS, the relationship between the spatial-frequency-dependent noise equivalent quanta (NEQ) of DPC-CT and ACT was derived. From these fundamental relationships, the NPS and NEQ of the DPC-CT system can be derived from the corresponding ACT system or vice versa. To validate these theoretical relationships, a benchtop cone beam DPC-CT/ACT system was used to experimentally measure the modulation transfer function (MTF) and NPS of both DPC-CT and ACT. The measured three-dimensional (3D) MTF and NPS were then combined to generate the corresponding 3D NEQ. RESULTS Two fundamental relationships have been theoretically derived and experimentally validated for the NPS and NEQ of DPC-CT and ACT: (1) the 3D NPS of DPC-CT is quantitatively related to the corresponding 3D NPS of ACT by an inplane-only spatial-frequency-dependent factor 1∕f (2), the ratio of window functions applied to DPC-CT and ACT, and a numerical factor C(g) determined by the geometry and efficiency of the grating interferometer. Note that the frequency-dependent factor is independent of the frequency component f(z) perpendicular to the axial plane. (2) The 3D NEQ of DPC-CT is related to the corresponding 3D NEQ of ACT by an f (2) scaling factor and numerical factors that depend on both the attenuation and refraction properties of the image object, as well as C(g) and the MTF of the grating interferometer. CONCLUSIONS The performance of a DPC-CT system is intrinsically related to the corresponding ACT system. As long as the NPS and NEQ of an ACT system is known, the corresponding NPS and NEQ of the DPC-CT system can be readily estimated using additional characteristics of the grating interferometer.

Grating-based phase contrast tomosynthesis imaging: Proof-of-concept experimental studies

PURPOSE This paper concerns the feasibility of x-ray differential phase contrast (DPC) tomosynthesis imaging using a grating-based DPC benchtop experimental system, which is equipped with a commercial digital flat-panel detector and a medical-grade rotating-anode x-ray tube. An extensive system characterization was performed to quantify its imaging performance. METHODS The major components of the benchtop system include a diagnostic x-ray tube with a 1.0 mm nominal focal spot size, a flat-panel detector with 96 μm pixel pitch, a sample stage that rotates within a limited angular span of ± 30°, and a Talbot-Lau interferometer with three x-ray gratings. A total of 21 projection views acquired with 3° increments were used to reconstruct three sets of tomosynthetic image volumes, including the conventional absorption contrast tomosynthesis image volume (AC-tomo) reconstructed using the filtered-backprojection (FBP) algorithm with the ramp kernel, the phase contrast tomosynthesis image volume (PC-tomo) reconstructed using FBP with a Hilbert kernel, and the differential phase contrast tomosynthesis image volume (DPC-tomo) reconstructed using the shift-and-add algorithm. Three inhouse physical phantoms containing tissue-surrogate materials were used to characterize the signal linearity, the signal difference-to-noise ratio (SDNR), the three-dimensional noise power spectrum (3D NPS), and the through-plane artifact spread function (ASF). RESULTS While DPC-tomo highlights edges and interfaces in the image object, PC-tomo removes the differential nature of the DPC projection data and its pixel values are linearly related to the decrement of the real part of the x-ray refractive index. The SDNR values of polyoxymethylene in water and polystyrene in oil are 1.5 and 1.0, respectively, in AC-tomo, and the values were improved to 3.0 and 2.0, respectively, in PC-tomo. PC-tomo and AC-tomo demonstrate equivalent ASF, but their noise characteristics quantified by the 3D NPS were found to be different due to the difference in the tomosynthesis image reconstruction algorithms. CONCLUSIONS It is feasible to simultaneously generate x-ray differential phase contrast, phase contrast, and absorption contrast tomosynthesis images using a grating-based data acquisition setup. The method shows promise in improving the visibility of several low-density materials and therefore merits further investigation.