X-ray computed tomography (CT) is important for various applications, such as medical diagnosis, industrial nondestructive testing, and security inspection. However, traditional x-ray CT image values provide information only on the linear attenuation coefficients weighted by the spectrum of the incident x-ray beam and cannot accurately identify tumors and constituents. Moreover, this drawback persists in other advanced imaging modes, including molecular imaging with dynamic contrast-enhanced CT.1 As a state-of-the-art technique, spectral CT with multienergy-bin photon-counting detectors enables K-edge imaging. This method uses the K-absorption edge discontinuity of the attenuation coefficients of high-atomic-number elements, such as gadolinium and gold, whose K-edges lie within 50.2 and 80.7 keV, respectively. These K-edge features of contrast agents can be reconstructed and extracted for material-specific quantitative evaluation using data above and below the K-edge energy.2 K-edge imaging aided by gold nanoparticles (NPs) has been used to detect plaque inflammation degree in coronary artery diseases.3 However, the sensitivity of K-edge imaging depends on the scanned object size and its material composition due to beam hardening.4
Different from spectral CT or K-edge imaging, x-ray fluorescence computed tomography (XFCT) detects the characteristic photons from a material excited by initial x-rays. XFCT is a quantitative imaging modality that reconstructs the two-dimensional or three-dimensional distribution of nonradioactive high-atomic-number NPs within a specimen, such as gold, gadolinium, and iodine.56.–7 Unlike that in K-edge CT images, the NP contrast in XFCT images is completely independent from the type of surrounding tissue.8 Therefore, XFCT outperforms K-edge CT in imaging contrast agents at low concentrations. Bazalova et al.8 and Feng et al.9 proved that the contrast-to-noise ratio (CNR) of XFCT images is higher than that of transmission K-edge CT images for contrast concentrations below 0.4%. Moreover, XFCT imaging provides more depth localization than fluorescence molecular tomography,10,11 bioluminescence tomography,12 and x-ray luminescence CT with near-IR-emitting nanophosphor.13,14 This advantage is important for preclinical and potential clinical applications. Therefore, with the development of x-ray-excitable NPs, XFCT is a promising modality for molecular imaging in fundamental studies, drug development, and clinical experiments.
Early XFCT is performed on high-intensity monochromatic synchrotron radiation sources with large facilities.1518.104.22.168.–20 The pencil beam/single collimated beam from polychromatic x-ray tubes222.214.171.124.26.–27 is then used to improve the feasibility of using XFCT; however, this technique requires a long scanning time, e.g., 279 min28 and 180 h.29 One important reason for the low efficiency is that using a pencil beam setup requires the beam source to be shifted point by point while measuring a single slice of projection data. Similar processes of shifting on detectors are also required for a fan-beam setup with a single pixel detector.30 Jones and Cho31 demonstrated the simulation results of XFCT using a polychromatic cone-beam source with detector arrays and indicated that the total scanning time can be significantly reduced. Cong et al.32 proposed a fan-beam XFCT setup with detector arrays and evaluated its feasibility through simulation. Apart from the shifting, these XFCT setups rely on spectrum scanning to obtain the spectral information required by Compton scatter compensation. For photon-counting detectors, spectrum scanning generates an enormous waste of x-ray dose, which aggravates the problem.
In this paper, we report the first experimental results (to the best of our knowledge) obtained using a full-field fan-beam XFCT consisting of a conventional low-intensity x-ray tube, a multibin photon-counting detector array, and a tungsten pinhole collimator.
Methods and Experimental Setup
(2) X-ray fluorescence (XRF) photons are emitted isotropically from position while incident x-ray photons interact with contrast agents. Its intensity can be expressed as
(3) At last, some XRF photons emitted from position travel through the object and are recorded by photon-counting detectors
Then, the intensity map of generated XRF photons can be reconstructed using the iterative maximum-likelihood expectation maximization (ML-EM) algorithm as34
Figure 2(a) shows a photograph of our XFCT imaging setup. A low-brilliance x-ray beam (90 kV, 0.5 mA) was generated using a conventional Hamamatsu L12161-07 x-ray tube. The x-ray beam was collimated using lead plates in two directions to produce a fan beam. To reduce the total dose and scattering photons, a 0.3-mm Cu sheet was placed at the x-ray tube exit. A home-made polymethyl methacrylate (PMMA) phantom was set to rotate in precise steps using a computer-controlled motion stage, whereas the x-ray tube and the detector were stationary. The third-generation CT scanning mode was used with a rotational step of 1 deg over a full scan. At each angle view, XFCT data emitted from the whole phantom illuminated using the full-field fan-beam were collected by a linear photon-counting detector array (eV3500, eV PRODUCTS, Saxonburg, Pennsylvania). As shown in Fig. 2(b), a 5-mm tungsten pinhole collimator was placed before the detector. This x-ray detector used a 3-mm CdZnTe semiconductor, 256 pixels at a 0.5-mm pixel pitch, and five adjustable energy thresholds.33 The detector was placed at a 90-deg angle to the incident x-ray beam to minimize the number of scattering photons entering the detector. Furthermore, a 3-mm lead cover was used to shield the whole detector box from x-rays coming from outside the field of view. Moreover, the same x-ray detector was repositioned behind the phantom along the beam direction to provide transmission spectral CT measurements and correct the attenuation. A PMMA phantom containing water and gadolinium solution insertions was prepared, as shown in Fig. 2(c). The insertions consisted of 1%, 2%, 3%, and 4% (weight/volume) gadolinium solutions. The phantom dimension was 35 mm in diameter.
Figure 3 shows the sinograms (projection data, ) of the XFCT experiment. Per rotation, 360 projections were acquired at a frame duration of 5 s. This scenario resulted in a total measurement time of 30 min per slice and a dose equivalent of 900 mAs. The energy thresholds were set at 35, 40, 45, 50, 70, and 90 keV. The peak of gadolinium was located at 42.74 keV; hence, the data within the 40- to 45-keV energy bin were used for XFCT. Figure 3(a) shows the raw XFCT data. Although the detector was energy calibrated,35,36 it had obvious streak artifacts because of the inconsistencies among the different detector elements. Figure 3(b) shows the gain correction results based on Fig. 3(a), wherein most of the streak artifacts were removed. However, the background was strong due to the scattered photons. The state-of-the-art XFCT systems rely on spectrum analysis to distinguish a signal from scattering; single-pixel spectrometers serve this function but generate an enormous waste of acquisition time and x-ray dose. Inspired by a method designed for single-photon emission computed tomography systems,37 we propose to compensate for the Compton scattering based on energy bins rather than spectra. Considering that the scattering angle of each pixel is known, we can calculate the Compton scattering cross section based on the Klein–Nishina formula. With the incident tube spectrum obtained by simulation, we can determine the energy distribution of scattering photons. Based on this distribution, we can then estimate the scattering at 40 to 45 keV, which is mixed with the 42.74-keV XRF signal, by direct measurement of the scattering on nearby bins of 35 to 40 and 45 to 50 keV because XRF signals are theoretically absent. The scattering on the target energy bin can be calculated asFigure 3(c) shows the scattering correction results based on Fig. 3(b), wherein the background was effectively reduced.
Figure 4 shows the XFCT reconstruction images using the ML-EM algorithm. Figure 4(a) shows the image reconstructed from the raw data of Fig. 3(b). The image shows an obvious background beside the gadolinium solution regions caused by the scattered photons. Figure 4(b) shows the image reconstructed from the scattering correction data of Fig. 3(c) where the strong background in Fig. 4(a) was effectively removed. In the 40- to 45-keV energy bin, there are averagely 18.41 photons per pixel per angle in the range of the phantom, while after the scattering correction the number drops to 3.47. Notably, XRF photons are also attenuated by the object before entering the detector. Therefore, attenuation correction should also be considered.38 The attenuation map can be reconstructed using additional or simultaneous transmission spectral CT scan. Figure 4(c) shows the attenuation correction result based on Fig. 4(b) using the attenuation map reconstructed from the transmission spectral CT data in the 40- to 45-keV energy bin. With attenuation correction, the brightness of Fig. 4(c) is much higher than that of Fig. 4(b).
Average gadolinium concentrations calculated from the results of Fig. 4.
To quantify the XFCT reconstructions in Fig. 4, we calculated the average gadolinium concentration values in the regions-of-interest (ROI) marked by the red circles. The results are shown in Table 1. In the second row, the values in Fig. 4(a) are larger than the truth because the XFCT data in Fig. 4(a) contain many scattered photons that increased the calculated concentrations. After the scattering correction, the values in the third row were reduced and even smaller than the truth. As shown in the fourth row, attenuation correction can improve the reconstructed concentrations to correspond with the true values.
For comparison, spectral CT imaging was performed using the same hardware and scanning parameters, e.g., kVp and mAs. The difference was that the detector array was placed against the central x-ray beam, and the energy thresholds were changed to 40, 45, 50, 55, 73, and 90 keV because the K-edge of gadolinium was located at 50.229 keV. Figure 5(a) shows the decomposition results of gadolinium from the spectral CT reconstruction. Material decomposition magnified the inconsistencies among different detector elements and resulted in strong ring artifacts.
CNR comparisons between Figs. 5(a) and 5(b).
|CNR for Fig. 5(a)||0.5342||2.1075||2.7421||3.9733|
|CNR for Fig. 5(b)||1.1651||3.0917||4.6448||7.9373|
The images in Fig. 5 were evaluated as CNR by calculating the ratio of the difference between the mean value of each ROI (red circles) and background (PMMA, blue circles) and standard deviation of the background.
As shown in Table 2, the CNR of the reconstructed XFCT image is higher than that of the corresponding spectral CT image due to its strong background. Moreover, XFCT may have a higher sensitivity than spectral CT (including K-edge CT) in molecular contrast imaging. Therefore, XFCT could be superior to spectral CT in the presence of low concentrations of contrast agents.
A study of full-field fan-beam XFCT using conventional low-intensity x-ray tube and energy-sensitive photon-counting detectors is presented. Unlike other pencil beam studies, a fan beam produced using a polychromatic x-ray tube was used to stimulate emissions of XRF photons that were collected by passing through a tungsten pinhole collimator. Thus, a whole XFCT scan can be completed within 30 min using third-generation CT geometry, which is faster than other XFCT devices that use a pencil beam with a first-generation CT geometry. Experimental results showed that scattering and attenuation correction are indispensable for highly accurate XFCT reconstruction. Moreover, the CNR of the XFCT image is higher than that of the corresponding spectral or K-edge CT image because a strong background exists adjacent to the contrast agent regions in spectral or K-edge CT images. This full-field XFCT has high sensitivity in molecular contrast imaging and is a promising modality for biomedical imaging of molecular probes. A multipinhole collimator and large detector arrays may be used in future XFCT experiments to improve the efficiency of using XRF photons.
This work was partially supported by National Natural Science Foundation of China, Award Nos: 61571256, 81427803, and 11525521; and Beijing Municipal Science and Technology Commission, Award No. Z151100003915079.
Liang Li received his PhD from Tsinghua University in Beijing, China, in 2007. Currently, he is an associate professor in the Department of Engineering Physics, Tsinghua University. His research interests center on the mathematical and physical problems of x-ray imaging and its medical, industrial, and other applications.