Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 391212, 6 pages

http://dx.doi.org/10.1155/2013/391212

Research Article

A Novel CT Imaging System with Adjacent Double X-Ray Sources

Mou An and Yaoqin Xie

Key Laboratory of Health Informatics, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, 1068 Xueyuan Avenue, Shenzhen 518055, China

Received 28 June 2013; Revised 14 September 2013; Accepted 19 October 2013

Academic Editor: Tianye Niu

Copyright © 2013 Mou An and Yaoqin Xie. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Current computed tomography (CT) scanners rotate fast to reduce motion artifact. X-ray tube must work in a high power to make the image clear under short exposure time. However, the life span of such a tube may be shortened. In this paper, we propose a novel double sources CT imaging system, which puts two of the same X-ray sources closely with each other. The system is different from current dual source CT with orthogonal X-ray sources. In our system, each projection is taken twice by these two sources to enhance the exposure value and then recovered to a single source projection for image reconstruction. The proposed system can work like normal single source CT system, while halving down the working power for each tube.

1. Introduction

Current computed tomography (CT) scanners acquire multiple projection images (~1000 frames) [1, 2]. A more powerful X-ray tube could emit more X-ray photons during one exposure frame and then cost less time to ensure constant exposure dose. Twice the power gives the possibility of increasing the scanning speed. However, the cathode filament electric current must be doubled to ensure emitting twice the number of electrons, which is harmful to the lifetime of the cathode filament and the positive plate. Zhang et al. [3–6] proposed a multiplexing radiography technology based on carbon nanotube field emission, which is an effective method to gather X-ray radiation from more than one tube to enhance the exposure rate in one frame. In this case, the X-ray sources must be modulated easily. Otherwise, the maximum tube voltage of these sources is limited by current technology.

In this paper, a novel CT imaging system is developed by placing two of the same X-ray tubes close to each other, along with a common detector. The system is different from current dual source CT system, where two sets of tube detectors are placed orthogonally [7–9]. In this system, the detector detects X-rays coming from two sources simultaneously. Two sources working together can make the exposure rate the same as a single source with half the exposure time.

In the following sections, first we describe the structure of adjacent double X-ray sources and then propose a method to separate the overlapped projections for image reconstruction. For demonstration, fan-beam data acquisition and image reconstruction are presented, but the method can be easily extended to cone beam image reconstruction. Lastly, we discuss the potential benefits and limitations of the proposed method.

2. Materials and Methods

2.1. CT System with Adjacent Double X-Ray Sources

The structure of traditional CT scanner system is illustrated in Figure 1(a). It can be seen that one X-ray source located in the focal point of the arc detector and synchronously rotates in anticlockwise direction along with the arc detector [10–12]. The rotation axis is the midpoint of the radius of curvature. Figure 1(b) gives us the proposed adjacent double X-ray sources system structure. In this system, one arc detector faces two X-ray sources. Different from Figure 1(a), there are two of the same X-ray sources that synchronously rotate with the arc detector. The rotation axis is no more the midpoint of the radius of curvature but the center of the arc detector.

(a)

(b)

Figure 1: (a) Current CT scanner system using a single-beam X-ray source and arc detector, whose focus is located in the source. (b) The proposed adjacent double X-ray sources system, in which the focus of the arc detector is located in the rotating center.

In this illustration about proposed double X-ray sources scanning system, sources are labeled as S1 and S2. As shown in Figure 1(b), θ is the angle between two sources relative to the rotation center. Two sources synchronously rotate in anticlockwise direction along with the arc detector. If we assume that the rotation speed is ω and the time between two consecutive sampling events is δ t, then . In the processing of projections acquisition, the position of S1 at one sampling event will be replaced by S2 in the next sampling event. And, as shown in Figure 1(b), the object must be covered by X-rays coming from both of the two sources. Both sources are turned on during the whole scanning procedure.

2.2. Projection Separating

Just like the Hadamard multiplexing radiography (HMR) method [5], we assume the Nth projection under two sources has the form and the Nth projection under one source has the form . The total projections number is . Let be the original projections serials; the overlapped projections serials are related to the original serials by a linear transform:

The -matrices consist of only 1 and 0, which correspond to the existence or not of the single source. The definite form is as follows:

where the row number of matrix means the acquisition sequence number in one scanning circle. And the column number of each row means the position of the sources.

As an example for the matrix of order , the convolution process that occurs in the experiment can be expressed succinctly in a matrix notation as

The number of rows of the S-matrices equals M, that is, the projections number, and must meet the condition . For a big enough M, the arc angle of detector satisfies the relationship , which means that after taking projections, the current position of projections received by the rotated detector is completely separated from that received on the original position, as shown in Figure 2.

Figure 2: Rotation of the adjacent double X-ray sources system. represents the two sources at the Nth projection angle in one scanning circle. After the sources detector pairs rotate degrees, the Lth projection is completely separated from the first projection .

For simplicity and clarity, in Figure 3, we changed the arc detector to flat. In this illustration, we assume that , , , and and are the projections about the same angle but the detector drifts. The vertical dashed line in Figure 3 cuts into two parts. It is required that the object must be covered by X-rays coming from both of the two sources; as a result, the right side of in is zero. On the other hand, the extension line part of in is zero too. The positions of the projected images in and totally coincide. As a result, it can be absolutely expressed by ; then we can subtract just by . In the same way, in equals ; we can also subtract just by . We do this thing until the projections drift totally outside the vertical dashed line; then we can recover . The general expression is performed as

Here, if , means the sources rotate more than one circle and equals . L means that after the sources rotate degrees, the projection is totally separated from the original projection .

Figure 3: Overlapped projections with the drifted detector. In two adjacent overlapped projections, there exist two original projected signals that come from the same angle source and that are accepted by the drifted detector.

2.3. Image Reconstruction

The flow chart of the general imaging and data processing procedure for this adjacent double X-ray sources system is shown in Figure 4. After the overlapped projections are separated, the traditional fan-beam reconstruction algorithm is applied for image reconstruction [13–15].

Figure 4: Flow chart of the imaging procedure of the adjacent double X-ray sources system. The incident X-ray comes from two sources transmitted through the object and formed an overlapped X-ray projection recorded by the arc X-ray detector. After a complete set of overlapped projections were acquired, the proposed algorithm was applied to the projection data to recover the original projection images.

By means of variable substitution, we get the FBP method for the fan-beam reconstruction [15]

where and is the projection accepted by arc detector D1 on the rotated degree , as illustrated in Figure 5. The opening angle relative to the source S between X-ray passed through the reconstructed point and the detector’s central line (dotted line) is . If we replace the detector D1 by D2, the received image by D2 is and . Then the FBP method (5) can be rewritten as

where .

Figure 5: Illustration of the two kinds of arc detectors. is the center of detector and is the center of detector .

3. Results and Discussion

To demonstrate the feasibility of the proposed system, the Shepp-Logan phantom was applied to simulate double sources exposure. The opening angle of the two sources to the rotation center was . The radius of the circular orbit was 60 cm; that is,  cm. The maximum fan-subtending angle was . The projections were evenly acquired in with a sampling angle of , and was equally sampled in with a total of 512 X-ray sums. Figure 6 gives the noise-free sinogram of the projections under one X-ray source, overlapped and separated under double X-ray sources. Because neither of the two X-ray sources was on the central line of the arc detector, after being separated, the projections (c) were a little vertically drifted compared with (a). During the reconstruction, this should be taken into account.

(a)

(b)

(c)

Figure 6: Sinogram of projections: from left to right it is under single X-ray source, overlapped and separated under double X-ray sources.

The Shepp-Logan phantom is shown in Figure 7(a). for comparison, the reconstructed images using (5) for single X-ray source and using (6) for double X-ray sources are shown in Figures 7(b) and 7(c). Each of them is 512 × 512 pixels and displayed based on the conventional greyscale window setup of [0.1, 0.3] [16]. The image reconstructed from the overlapped projections seems to have the same quality as the image reconstructed from the projections under single X-ray source. Mean square error (MSE) of Figure 7(b) is 0.039 and for Figure 7(c) is 0.040, respectively. Quantitative agreement between the reconstructed image in Figure 7(c) and the true phantom in Figure 7(a) can be seen in Figure 8.

(a)

(b)

(c)

Figure 7: (a) Shepp-Logan phantom. (b) Reconstructed image from the projections under single X-ray source. (c) Reconstructed image from the signals separated from the overlapped projections under two X-ray sources. Display window: [0.1 0.3].

(a)

(b)

Figure 8: Profile comparison between reconstructed (solid red line) and true image (dashed blue line) on the horizontal and vertical central line.

In the simulation above, we present preliminary studies that demonstrate image reconstruction noise-free with the overlapped projections under adjacent double X-ray sources. In this imaging system, two X-ray sources must be the same, that is, the machine model, radiation field of photons. In practice, radiation field could be adjusted by the working voltage and current. Our results are provided about the fan-beam reconstruction, but the method is suitable for cone beam CT; just make sure that the two adjacent X-ray sources are placed parallel to the rotating direction and that the opening angle of the two sources to rotation center agree with the acquisition angle. In normal CT scanning, we take hundreds of projections during one circle. The angle rotated between two sampling is so small that two X-ray sources must be placed close to each other. We could place two identical cathode filament in one X-ray tube. However, sparse projections reconstruction acquires fewer projections than conventional methods; the sampling angle is relatively big, which facilitates placing two X-ray tubes together.

As compared to traditional CT system, the proposed adjacent double X-ray sources system can squeeze double X-ray photons during one unit exposure time; thus, it can easily enhance the illumination for some situations that need a high exposure rate such as fat patients. Also, high exposure rate during one unit exposure time gives the ability to halve down the projections acquisition time to reduce the motion artifact. On the other hand, if fast exposure is not required, two sources that work in a half power can give enough X-ray photons, which can prolong the lifetime of the tubes.

The main advantage of such a configuration is the power, since the two tubes can operate simultaneously. However, several issues should be tackled before it can be used in clinic. First, scatter is a big issue that can influence the quality of the reconstructed images. Figure 9 is the schematic diagram of the antiscatter collimators. Each lead inclined to one focus, that is, the source. In our proposed system, two sources are placed in cross-section perpendicular to the rotating axis and parallel to the rotating direction. This kind of antiscatter collimators could not be used. But we can take out the leads on one direction; for example, we take out the leads parallel to -axis. Then the rest of the lead inclined to line across the two sources. This is a compromised way to reduce part of the scatter. Another big issue is the CNR after the preprocessing step in (4). According to the result in [5], the separated projections under three sources are still clear enough.

Figure 9: Antiscatter collimators.

4. Conclusions

In this paper, we developed a novel structure of CT system with adjacent double X-ray sources. In the proposed system, two X-ray sources work simultaneously. Therefore, the tube can work longer under the same exposure dose as traditional single source or it can give a high exposure rate for some situations such as fat patients. Because the sources are not necessary to be modulated and the reconstruction method is based on the FBP algorithm, the system can be realized by current X-ray tube technology and CT structure.

Acknowledgments

This work was supported in part by Grants from the National Natural Science Foundation of China (NSFC: 81171402), the NSFC Joint Research Fund for Overseas Research Chinese, the Hong Kong and Macao Young Scholars (30928030), the National Basic Research Program 973 (2010CB732606) from the Ministry of Science and Technology of China, and the Guangdong Innovative Research Team Program (no. 2011S013) of China.

References

J. Hsieh, Computed Tomography: Principles, Design, Artifacts, and Recent Advances, SPIE Press, 2003.

E. Seeram, Computed Tomography, W.B. Saunders, Philadelphia, Pa, USA, 2001.

J. Zhang, G. Yang, Y. Z. Lee, S. Chang, J. P. Lu, and O. Zhou, “Multiplexing radiography using a carbon nanotube based x-ray source,” Applied Physics Letters, vol. 89, no. 6, Article ID 064106, 2006.

J. Zhang, G. Yang, Y. Lee, S. Chang, J. P. Lu, and O. Zhou, “Multiplexing radiography based on carbon nanotube field emission x-ray technology,” in Medical Imaging 2007: Physics of Medical Imaging, vol. 6510 of Proceedings of the SPIE, San Diego, Calif, USA, February 2007.

J. Zhang, G. Yang, S. Chang, J. P. Lu, and O. Zhou, “Hadamard multiplexing radiography based on carbon nanotube field emission multi-pixel x-ray technology,” in Medical Imaging 2008: Physics of Medical Imaging, vol. 6913 of Proceedings of the SPIE, San Diego, Calif, USA, February 2008.

G. Cao, J. Zhang, O. Zhou, and J. Lu, “Temporal multiplexing radiography for dynamic x-ray imaging,” Review of Scientific Instruments, vol. 80, no. 9, Article ID 093902, 2009.

M. Petersilka, H. Bruder, B. Krauss, K. Stierstorfer, and T. G. Flohr, “Technical principles of dual source CT,” European Journal of Radiology, vol. 68, no. 3, pp. 362–368, 2008.

S. Nicolaou, A. Eftekhari, T. Sedlic et al., “The utilization of dual source CT in imaging of polytrauma,” European Journal of Radiology, vol. 68, no. 3, pp. 398–408, 2008.

T. G. Flohr, C. H. McCollough, H. Bruder, et al., “First performance evaluation of a dual-source CT (DSCT) system,” European Radiology, vol. 16, no. 2, pp. 256–268, 2006.

G. T. Gullberg, “The reconstruction of fan-beam data by filtering the back-projection,” Computer Graphics and Image Processing, vol. 10, no. 1, pp. 30–47, 1979.

F. Natterer, “Sampling in fan beam tomography,” SIAM Journal on Applied Mathematics, vol. 53, no. 2, pp. 358–380, 1993.

G.-H. Chen, R. Tokalkanahalli, T. Zhuang, B. E. Nett, and J. Hsieh, “Development and evaluation of an exact fan-beam reconstruction algorithm using an equal weighting scheme via locally compensated filtered backprojection (LCFBP),” Medical Physics, vol. 33, no. 2, pp. 475–481, 2006.

B. K. P. Horn, “Fan-beam reconstruction methods,” Proceedings of the IEEE, vol. 67, no. 12, pp. 1616–1623, 1979.

G. Besson, “CT fan-beam parametrizations leading to shift-invariant filtering,” Inverse Problems, vol. 12, no. 6, pp. 815–833, 1996.

G. Zeng, Medical Image Reconstruction, Higher Education Press, 2010.

T. Niu and L. Zhu, “Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: phantom studies,” Medical Physics, vol. 39, no. 7, article 4588, 2012.