Zhihua Qi

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.

Image reconstruction for fan-beam differential phase contrast computed tomography

Recently, x-ray differential phase contrast computed tomography (DPC-CT) has been experimentally implemented using a conventional tube combined with gratings. Images were reconstructed using a parallel-beam reconstruction formula. However, parallel-beam reconstruction formulae are not applicable when the parallel-beam approximation fails. In this paper, we present a new image reconstruction formula for fan-beam DPC-CT. There are several novel features of the new image reconstruction formula: (i) when the scanning angular range of data acquisition is larger than pi + gamma(m) (gamma(m) is the full fan angle), the entire field of view can be exactly reconstructed; (ii) when the scanning angular range is smaller than pi + gamma(m), a local region of interest (ROI) can be exactly reconstructed; (iii) it enables an exact reconstruction for a local ROI when the projection data are truncated at some view angles; (iv) it enlarges the imaging field of view when the detector is asymmetrically placed. In this last case, the data are truncated from every view angle. Numerical simulations have been conducted to validate the new reconstruction formula.

Performance studies of four-dimensional cone beam computed tomography

Four-dimensional cone beam computed tomography (4DCBCT) has been proposed to characterize the breathing motion of tumors before radiotherapy treatment. However, when the acquired cone beam projection data are retrospectively gated into several respiratory phases, the available data to reconstruct each phase is under-sampled and thus causes streaking artifacts in the reconstructed images. To solve the under-sampling problem and improve image quality in 4DCBCT, various methods have been developed. This paper presents performance studies of three different 4DCBCT methods based on different reconstruction algorithms. The aims of this paper are to study (1) the relationship between the accuracy of the extracted motion trajectories and the data acquisition time of a 4DCBCT scan and (2) the relationship between the accuracy of the extracted motion trajectories and the number of phase bins used to sort projection data. These aims will be applied to three different 4DCBCT methods: conventional filtered backprojection reconstruction (FBP), FBP with McKinnon-Bates correction (MB) and prior image constrained compressed sensing (PICCS) reconstruction. A hybrid phantom consisting of realistic chest anatomy and a moving elliptical object with known 3D motion trajectories was constructed by superimposing the analytical projection data of the moving object to the simulated projection data from a chest CT volume dataset. CBCT scans with gantry rotation times from 1 to 4 min were simulated, and the generated projection data were sorted into 5, 10 and 20 phase bins before different methods were used to reconstruct 4D images. The motion trajectories of the moving object were extracted using a fast free-form deformable registration algorithm. The root mean square errors (RMSE) of the extracted motion trajectories were evaluated for all simulated cases to quantitatively study the performance. The results demonstrate (1) longer acquisition times result in more accurate motion delineation for each method; (2) ten or more phase bins are necessary in 4DCBCT to ensure sufficient temporal resolution in tumor motion and (3) to achieve the same performance as FBP-4DCBCT with a 4 min data acquisition time, MB-4DCBCT and PICCS-4DCBCT need about 2- and 1 min data acquisition times, respectively.

Extraction of tumor motion trajectories using PICCS‐4DCBCT: A validation study

PURPOSE As a counterpart of 4DCT in the treatment planning stage of radiotherapy treatment, 4D cone beam computed tomography (4DCBCT) method has been proposed to verify tumor motion trajectories before radiation therapy treatment delivery. Besides 4DCBCT acquisition using slower gantry rotation speed or multiple rotations, a new method using the prior image constrained compressed sensing (PICCS) image reconstruction method and the standard 1-min data acquisition were proposed. In this paper, the PICCS-4DCBCT method was combined with deformable registration to validate its capability in motion trajectory extraction using physical phantom data, simulated human subject data from 4DCT and in vivo human subject data. METHODS Two methods were used to validate PICCS-4DCBCT for the purpose of respiratory motion delineation. The standard 1-min gantry rotation Cone Beam CT acquisition was used for both methods. In the first method, 4DCBCT projection data of a physical motion phantom were acquired using an on-board CBCT acquisition system (Varian Medical Systems, Palo Alto, CA). Using a deformable registration method, the object motion trajectories were extracted from both FBP and PICCS reconstructed 4DCBCT images, and compared against the programmed motion trajectories. In the second method, using a clinical 4DCT dataset, Cone Beam CT projections were simulated by forward projection. Using a deformable registration method, the tumor motion trajectories were extracted from the reconstructed 4DCT and PICCS-4DCBCT images. The performance of PICCS-4DCBCT is assessed against the 4DCT ground truth. The breathing period was varied in the simulation to study its effect on motion extraction. For both validation methods, the root mean square error (RMSE) and the maximum of the errors (MaxE) were used to quantify the accuracy of the extracted motion trajectories. After the validation, a clinical dataset was used to demonstrate the motion delineation capability of PICCS-4DCBCT for human subjects. RESULTS In both validation studies, the RMSEs of the extracted motion trajectories from PICCS-4DCBCT images are less than 0.7 mm, and their MaxEs are less than 1 mm, for all three directions. In comparison, FBP-4DCBCT shows considerably larger RMSEs in the physical phantom based validation. PICCS-4DCBCT also shows insensitivity to the breathing period in the 4DCT based validation. For the in vivo human subject study, high quality 3D motion trajectory of the tumor was obtained from PICCS-4DCBCT images and showed consistency with visual observation. CONCLUSIONS These results demonstrate accurate delineation of tumor motion trajectory can be achieved using PICCS-4DCBCT and the standard 1-min data acquisition.

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.