Abstract
An x-ray fluorescence computed tomography using a pinhole collimator (p-XFCT) for medical use depicts the 3-D distribution of imaging agent by collecting the x-ray fluorescence photons emitted from the imaging agent accumulated in an organ using a pinhole collimator and a 2-D detector. However, the number of the fluorescent x-ray photons are restricted by the pinhole, leading to acquisition of low S/N projections due to the photon starvation. Simultaneous acquisition of multiple projections using a multi-pinhole scheme overcomes the photon starvation. In this research, we clarified the measurement process and described the quantitative linear relationship between the quantities to be estimated and to be measured. We solved the multi-pinhole based XFCT (mp-XFCT) as a linear inverse problem regarding the described measurement process using ML-EM (maximum likelihood expectation-maximization) algorithm. We actually constructed an mp-XFCT imaging system at beamline AR-NE7A in KEK, and performed imaging experiment using a physical phantom to demonstrate the efficacy of mp-XFCT.
You have full access to this open access chapter, Download conference paper PDF
Similar content being viewed by others
Keywords
1 Introduction
In preclinical research, rodent models of human diseases yielded by genetic and tissue engineering are frequently used for elucidating their causes and developing their diagnosis and treatment. Nuclear imaging techniques such as positron emission tomography (PET) and single photon emission computed tomography (SPECT) are very important tools for in vivo observing functional information of the organs, i.e., the physiological and pathological characteristics, from the molecular viewpoint. However, they require radioactive imaging agents, leading to the measurement difficulty. Thus, a novel molecular imaging technique using non-radioactive imaging agents would be advantageous. Combining x-ray fluorescence analysis, which is one of the most sensitive physico-chemical analysis, with computed tomography, which provides the 3-D information, enables us to establish a low-invasive or nondestructive molecular imaging technique. So far, several types of the x-ray fluorescence computed tomography (XFCT) imaging scheme have been proposed [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. Among them a pinhole based XFCT (p-XFCT) outperforms the others in establishing both spatial resolution and fast scanning [13]. However, since the number of detected fluorescent photons are remarkably restricted by a pinhole, leading to acquisition of low S/N projections due to the photon starvation. Use of multiple pinholes simultaneously collects multiple projections, and then compensates the photon starvation effect [16, 17]. In this research, we clarified the measurement process of the multi-pinhole based XFCT (mp-XFCT) required for reconstruction. We actually constructed the mp-XFCT imaging system at beamline AR-NE7A in KEK, and performed an imaging experiment using a physical phantom to demonstrate the efficacy of mp-XFCT.
2 Method
2.1 Multi-Pinhole Based XFCT
Figure 1 shows the imaging geometry of mp-XFCT, which has seven pinholes, within the xyz-coordinate system. A monochromatic and parallel incident beam propagates along the x-axis, and a pinhole collimator is set so that its surface is perpendicular to the z-axis and the center of the pinhole lies on the z-axis. The surface of a two-dimensional detector is set perpendicular to the z-axis. The volumetric incident beam irradiates an object positioned in the vicinity of the origin, with the volumetric beam covering the object. Imaging agents in the object, such as I, Gd, and Au, are excited by the incident beam and isotropically emit fluorescent x-ray photons on de-excitation. Of all the fluorescent x-ray photons emitted from all of the imaging agents in the object, only the photons passing through the pinholes are acquired by the detector. Each view produces multiple projections corresponding to the individual pinholes. The projection acquisition is repeated while rotating the object around the y-axis. In the following, we assume that Kα1 fluorescence is dominant, although other fluorescents such as Kα2 and Kβ, can easily be incorporated into the subsequent expressions, despite a more complicated formulation.
2.2 Measurement Process for a Single Pinhole
Let us consider quantitatively how the fluorescent x-ray photons are detected at an element on the detective surface for a single pinhole. As shown in Fig. 1, consider a point Q inside the object on an incident x-ray flux, and a point S at the center of the detective element. We assume that point P is an intersection of the incident flux and the object surface, point Q is an intersection of the incident flux and the line connecting S to the pinhole center, and point R is an intersection of a line connecting Q to S and the object surface.
First, we will consider the measurement process when the above assumption holds. The process is divided into the following three subprocesses:
-
(a)
attenuation of the incident beam by the object from P to Q;
-
(b)
emission of fluorescent x-rays at Q; and
-
(c)
attenuation of the fluorescent x-rays by the object from Q to R, which then travel to S.
In subprocess (a), the density of incident x-ray photons having reached Q is given as I(a)(rQ) = I0 t exp(\( - \,\int {_{{\overline{PQ} }} } \;\mu_{a} \left( {E_{I} ,\,{\mathbf{r}}} \right)\,\;\text{d}l \)) [photons/mm2] owing to attenuation in line segment \( \overline{PQ} \), where rQ is the coordinate of Q, I0 [photons/mm2/s] is the number of incident photons per unit area and time, EI [keV], and μa(E, r) [1/mm] are the energy of incident x-ray and the three-dimensional map of the linear total attenuation coefficient inside the object at incident energy E, respectively, and \( t \)[s] is the exposure time.
In subprocess (b), the quantity of fluorescent x-rays emitted at Q is proportional to the product of the incident photon density I(a)(rQ) [photons/mm2] and the concentration of imaging element λ(rQ) [g/mm3], given as I(b)(rQ) = μph ω I(a)(rQ) λ(rQ) ΔrQ [photons], where μph [mm2/g] and ω are the photoelectric mass absorption coefficient of trace element and fluorescent yield, respectively, and ΔrQ [mm3] is an infinitesimal volume in the vicinity of Q.
In subprocess (c), of all the fluorescent x-ray fluxes emitted isotropically from Q, only the fluxes passing through the pinhole reach the detector, and the passing fluorescent fluxes are attenuated by the object. This attenuation is represented by the attenuation in the line segment \( \overline{QR} \), as the area of the pinhole is sufficiently small, although actually each flux is differently attenuated according to the propagation direction. Therefore, the fluorescent quantity emitted from Q and reaching S is given by I(c)(rS, rQ) = I(b)(rQ) (Ω(rS, rQ) /4π) exp(\( - \,\int {_{{\overline{QR} }} } \;\mu_{a} (E_{F} ,\,{\mathbf{r}}\text{)}\;\,\text{d}l \)) [photons], where rS is the coordinate of S, Ω(rS, rQ) [steradian] and EF [keV] are the solid angle subtended by Q through the pinhole toward the detective element whose center is S, and the energy of the fluorescent x-ray, respectively. For the subsequent discussion, we denote the fluorescent x-ray photons that are counted at the detective element at the coordinate whose center is rS as
where
and η is the quantum efficiency of the detector. If we specify the imaging geometry and the coordinates of S and Q, p(rS, rQ) can be known in advance, as the maps of μa(EI, r) and μa(EF, r) can be measured by ordinary computed tomography with monochromatic x-rays.
The above discussion concerns the contribution of fluorescent x-ray counts at S from Q, although the contributions from all regions around Q should actually be considered. Therefore, the total measured count at S is given as
where VS is a region in the object inside a cone subtended by S toward the pinhole.
For reconstruction, we discretize the object with voxels as shown in Fig. 2. Here we assume that in the projection acquisition process, the object is fixed to the xyz-coordinate system, and that the incident beam, the pinhole collimator, and the detective surface are rotated around the y-axis, with their positional relationship maintained, although in the actual measurement the object is rotated. We introduce the index j (= 1, 2, …, N) to identify the voxel in the lexicographic ordering corresponding to rQ, where N is the number of voxels. Additionally, we introduce the index i (= 1, 2, …, M) to identify the position of the individual detective elements throughout the whole projection acquisition, corresponding to rS, where the positions of detective elements throughout the whole projection acquisition are consecutively numbered, and M is the number of the positions. As a result, Eq. (3) is discretized to
where yi, pij, Vi and λj correspond to y(rS), p(rS, rQ)ΔrQ, VS and λ(rQ), respectively. yi and pij are known, and λj requires estimation.
2.3 Reconstruction for Multi-pinhole XFCT
Since use of a multi-pinhole collimator produces projections for the individual pinholes, Eq. (1) is rewritten as
where k is an index to identify the pinhole and L is the number of pinholes. As a result, we obtain a system of M × L linear equations. Then, we obtain L matrix equations yk = Pk λ (1 ≤ k ≤ L), where yk = (y k1 , y k2 , …, y k M )T, λ = (λ1, λ2, …, λN)T, and
In addition, we introduce a K (= L × M) -dimensional vector ψ and a K × N matrix Π, where
, and
As a result, a linear inverse problem,
should be solved, where vector ψ and matrix Π are known, and vector λ is to be estimated. For reconstruction, we straightforwardly applied the maximum likelihood expectation-maximization (ML-EM) algorithm to our liner inverse problem, ψ = Πλ [18, 19].
3 Experiment and Result
We constructed the mp-FXCT imaging system using the bending-magnet beamline AR-NE7A (6.5 GeV) at KEK. Since in this research we paid attention to K-shell fluorescence of I (Kα1: 28.3 keV), we by a double-crystal Bragg-Bragg monochromator using Si (111) single crystals selected the incident energy to 33.4 keV just above the K-edge energy of I (33.2 keV), so that we can generate as many fluorescent photons as possible. The flux rate was about 5.0 × 108 photons/mm2/s in front of the object. The cross section of the incident beam was collimated to 35 mm (horizontal) × 5 mm (vertical) by the slit.
A Pilatus 100 K with 487 × 195 elements (pixel size: 172 × 172 μm2) manufactured by DECTRIS Ltd. was used [15,16,17]. This detector was used because it offers a high signal-to-noise ratio thanks to no dark current, and it also offers a higher data transmission rate than other types of 2-D detectors, such as CCD-based detectors.
We fabricated a multi-pinhole collimator with 15 pinholes aligned in three rows containing five pinholes each. The collimator consisted of a 1-mm-thick PB plate and 15 W pinhole tips [16]. Each pinhole tip had a 0.1-mm-diameter hole. The pinholes were arranged so that the projections would not overlap.
The distance between the sample rotational axis and the collimator plane, and between the collimator plane and the detector surface were 27.8 mm, and 32.2 mm, respectively. The rotational stage and the detector were controlled by a PC.
We prepared a physical phantom which is a 10-mm-diameter PMMA made cylinder with six 2-mm-diameter channels filled with different concentrations of I solution (0.012, 0.025, 0.005, 0.1, 0.2, 0.3 I mg/ml), as shown in Fig. 2. The phantom was imaged with the exposure time per view of 2 min. 180 projections were obtained at a constant angular step of 2 degrees over 360 degrees.
We solved the inverse problem using the ordered subsets expectation maximization (OSEM) method which is an accelerated version of ML-EM and used widely in PET and SPECT reconstruction [18, 19]. The algorithm was implemented using C# on a PC with Windows 7, Intel Xeon CPU E5-2620, and 192 GB of memory. We set the number of subsets and of iterations were 3 and 5, respectively. The voxel size was 70 × 70 × 40; An edge of a voxel was 0.172 mm long, which is the same as the pixel dimension of the Pilatus detector.
Figure 3 depicts cross-sections acquired at the central levels of the 3-D reconstructed images (a) and their 3-D volume rendering images (b), where the left and right columns were reconstructed using a single pinhole and 15 pinholes, respectively. The back ground is noisier and the I-containing regions are less smooth in single-pinhole case than in the multi-pinhole case.
4 Conclusion
We described the imaging principle of multi-pinhole based x-ray fluorescence computed tomography (mp-XFCT) from the viewpoint of the measurement process. Then, we performed a phantom imaging experiment using an actual imaging system constructed in AR-NE7A and compare the result with the result in a single pinhole case to demonstrate that mp-XFCT improve an image quality by compensating the photon starvation effect caused by a pinhole.
In the present reconstruction process, scatter correction was not incorporated. In the experiment, the phantom size was relatively small, and then scatter components can be neglected. However, in a case of measuring a small animal much more scatter events will occur. Our future direction will be to devise the scatter correction and incorporate it with the reconstruction algorithm.
References
Cesareo, R., et al.: Nucl. Instrum. Meth. A 277, 559–672 (1989)
Hogan, J.P., et al.: IEEE Trans. Nucl. Sci. 38, 1721–1727 (1991)
Schroer, C.G.: Appl. Phys. Lett. 79, 1912–1914 (2001)
Glosio, B., et al.: J. Appl. Phys. 94, 145–156 (2003)
Rust, G.F., et al.: IEEE Trans. Nucl. Sci. 45, 75–88 (1998)
Yuasa, T., et al.: IEEE Trans. Nucl. Sci. 44, 54–62 (1997)
Yuasa, T., et al.: IEEE Trans. Nucl. Sci. 44, 1760–1769 (1997)
Huo, Q., et al.: Opt. Lett. 33, 2494–2496 (2008)
Huo, Q., et al.: X-Ray Spectrosc. 38, 439–445 (2009)
Takeda, T., et al.: J. Synchrotron Radiat. 16, 57–62 (2009)
Lwin, T.T., et al.: J. Synchrotron Radiat. 14, 158–162 (2007)
Lwin, T.T., et al.: J. Synchrotron Radiat. 15, 528–531 (2008)
Sunaguchi, N., et al.: Opt. Comm. 297, 210–215 (2013)
Bazalova-Carter, M., et al.: Radiother. Oncol. 15, S66–S57 (2015)
Sasaya, T., et al.: Sci. Rep. 7, 44143 (2017)
Sasaya, T., et al.: Sci. Rep. 7, 5742 (2017)
Sasaya, T., et al.: Nucl. Instrum. Meth. A 886, 71–76 (2018)
Lange, K., Garson, R.: J. Comput. Assist. Tomogr. 8, 306–316 (1984)
Shepp, L.A., Vardi, V.: IEEE Trans. Med. Imaging 1, 113–121 (1982)
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2018 Springer International Publishing AG, part of Springer Nature
About this paper
Cite this paper
Sasaya, T. et al. (2018). Multi-Pinhole Based X-Ray Fluorescence Computed Tomography: A Comparison with Single Pinhole Case. In: Mansouri, A., El Moataz, A., Nouboud, F., Mammass, D. (eds) Image and Signal Processing. ICISP 2018. Lecture Notes in Computer Science(), vol 10884. Springer, Cham. https://doi.org/10.1007/978-3-319-94211-7_23
Download citation
DOI: https://doi.org/10.1007/978-3-319-94211-7_23
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-94210-0
Online ISBN: 978-3-319-94211-7
eBook Packages: Computer ScienceComputer Science (R0)