Abstract
Sparse angular sampling has been studied recently owing to its potential to decrease the radiation exposure from computed tomography (CT). In this study, we investigated the analytic reconstruction algorithm in sparse angular sampling using the sinogram interpolation method for improving image quality and computation speed. A prototype of the spectral CT system, which has a 64-pixel Cadmium Zinc Telluride (CZT)-based photon-counting detector, was used. The source-to-detector distance and the source-to-center of rotation distance were 1,200 and 1,015 mm, respectively. Two energy bins (23∼33 keV and 34∼44 keV) were set to obtain two reconstruction images. We used a PMMA phantom with height and radius of 50.0 mm and 17.5 mm, respectively. The phantom contained iodine, gadolinium, calcification, and lipid. The Feld-kamp-Davis-Kress (FDK) with the sinogram interpolation method and Maximum Likelihood Expectation Maximization (MLEM) algorithm were used to reconstruct the images. We evaluated the signal-to-noise ratio (SNR) of the materials. The SNRs of iodine, calcification, and liquid lipid were increased by 167.03%, 157.93%, and 41.77%, respectively, with the 23∼33 keV energy bin using the sinogram interpolation method. The SNRs of iodine, calcification, and liquid state lipid were also increased by 107.01%, 13.58%, and 27.39%, respectively, with the 34∼44 keV energy bin using the sinogram interpolation method. Although the FDK algorithm with the sinogram interpolation did not produce better results than the MLEM algorithm, it did result in comparable image quality to that of the MLEM algorithm. We believe that the sinogram interpolation method can be applied in various reconstruction studies using the analytic reconstruction algorithm. Therefore, the sinogram interpolation method can improve the image quality in sparse-angular sampling and be applied to CT applications.
REFERENCES
1.Kalke M., Siltanen S. Sinogram interpolation method for sparse-angle tomography. Appl Math. 5(3):423–441. 2014.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
2.Brooks RA., Weiss GH., Talbert AJ. A new approach to interpolation in computed tomography. J Comput Assisted Tomogr. 2(5):577–585. 1978.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
3.Kim D., Park S-J., Jo B., Kim H., Kim H-J. Investigation of sparse-angle view in cone beam computed tomography (CBCT) reconstruction algorithm using a sinogram interpolation method. IFMBE Proceedings. 2015, Canada, pp.74.
4.Siltanen S., Kolehmainen V., Järvenpää S, et al. Statistical inversion for medical x-ray tomography with few radiographs: I. General theory. Phys Med Biol. 48(10):1465–1490. 2003.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
5.Lahart MJ. Estimation of reconstructions in computed tomography. J Opt Soc Am. 71(10):1155–1161. 1981.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
6.Shikhaliev PM., Xu T., Molloi S. Photon counting CT: Concept and initial results. Med Phys. 32(2):427–436. 2005.
7.Lundqvist M., Cederstrom B., Chmill V., Danielsson M., Hasegawa B. Evaluation of a photon counting X-ray imaging system. IEEE Trans Nucl Sci. 48(4):1530–1536. 2001.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
8.Sarkar V., Shi C., Rassiah-Szegedi P, et al. The effect of a limited number of projections and reconstruction algorithms on the image quality of megavoltage digital tomosynthesis. J Appl Clin Med Phys. 10(3):155–172. 2009.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
9.Feldkamp LA., Davis LC., Kress JW. Practical cone-beam algorithm. J Opt Soc Am A. 1(6):612–619. 1984.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
10.Shepp LA., Vardi Y. Maximum likelihood reconstruction for emission tomography. IEEE Trans Med Imag. 1(2):113–122. 1982.
![crossref](/image/icon/bnr_ref_cross.gif)
![crossref](/image/icon/bnr_ref_cross.gif)
11.Dempster AP., Laird NM., Rubin DB. Maximum likelihood from incomplete data via the EM algorithm. J R Statist Soc B. 39(1):1–38. 1977.
12.https://en.wikipedia.org/wiki/Signal-to-noise_ratio_(imaging).
Fig. 4
Restored PMMA phantom sinogram. (a) Without interpolation in 36 views, (b) With interpolation in 36 angular views.
![pmp-27-105f4.tif](/upload/SynapseXML/2129pmp/thumb/pmp-27-105f4.gif)
Fig. 5
The reconstruction images of the PMMA phantom with the FDK and MLEM algorithms. (a∼c) Non-enhanced iodine image (23∼33 keV): (a) FDK without the sinogram interpolation method, (b) FDK with the sinogram interpolation method, (c) MLEM (10 iterations). (d∼e) Enhanced iodine image (34∼44 keV): (d) FDK without the sinogram interpolation method, (e) FDK with the sinogram interpolation method, (f) MLEM (10 iterations).
![pmp-27-105f5.tif](/upload/SynapseXML/2129pmp/thumb/pmp-27-105f5.gif)
Fig. 6
The SNRs in the reconstruction images by FDK with the sinogram interpolation method, FDK without the sinogram interpolation method, and MLEM (10 iterations) (23∼33 keV).
![pmp-27-105f6.tif](/upload/SynapseXML/2129pmp/thumb/pmp-27-105f6.gif)
Fig. 7
The SNRs in reconstruction images by FDK with the sinogram interpolation method, FDK without the sinogram interpolation method, and MLEM (10 iterations) (34∼44 keV).
![pmp-27-105f7.tif](/upload/SynapseXML/2129pmp/thumb/pmp-27-105f7.gif)
Fig. 8
Time consumption for reconstructing images by FDK with the sinogram interpolation method, FDK without the sinogram interpolation method, and MLEM (10 iterations).
![pmp-27-105f8.tif](/upload/SynapseXML/2129pmp/thumb/pmp-27-105f8.gif)
Table 1.
Detailed acquisition parameters of the spectral CT system.