Bone-induced streak artifact suppression in sparse-view CT image reconstruction
© Jin et al.; licensee BioMed Central Ltd. 2012
Received: 10 February 2012
Accepted: 19 June 2012
Published: 2 August 2012
In sparse-view CT imaging, strong streak artifacts may appear around bony structures and they often compromise the image readability. Compressed sensing (CS) or total variation (TV) minimization-based image reconstruction method has reduced the streak artifacts to a great extent, but, sparse-view CT imaging still suffers from residual streak artifacts. We introduce a new bone-induced streak artifact reduction method in the CS-based image reconstruction.
We firstly identify the high-intensity bony regions from the image reconstructed by the filtered backprojection (FBP) method, and we calculate the sinogram stemming from the bony regions only. Then, we subtract the calculated sinogram, which stands for the bony regions, from the measured sinogram before performing the CS-based image reconstruction. The image reconstructed from the subtracted sinogram will stand for the soft tissues with little streak artifacts on it. To restore the original image intensity in the bony regions, we add the bony region image, which has been identified from the FBP image, to the soft tissue image to form a combined image. Then, we perform the CS-based image reconstruction again on the measured sinogram using the combined image as the initial condition of the iteration. For experimental validation of the proposed method, we take images of a contrast phantom and a rat using a micro-CT and we evaluate the reconstructed images based on two figures of merit, relative mean square error and total variation caused by the streak artifacts.
The images reconstructed by the proposed method have been found to have smaller streak artifacts than the ones reconstructed by the original CS-based method when visually inspected. The quantitative image evaluation studies have also shown that the proposed method outperforms the conventional CS-based method.
The proposed method can effectively suppress streak artifacts stemming from bony structures in sparse-view CT imaging.
Sparse-view CT is of great importance in clinical imaging for its potential to reduce the x-ray dose to the human subject and the scan time [1–3]. In sparse-view CT, less number of projection views than is required to satisfy the Nyquist sampling theorem is employed. Conventional filtered backprojection (FBP) based image reconstruction methods gives severe streak artifacts, sort of aliasing artifacts, in the images, which would hamper clinical utility of the sparse-view CT. Bony structures makes strongest streak artifacts, and physiological motions of the human subject, beam hardening, and photon starvation also make streak artifacts [4–7]. In sparse-view CT image reconstruction, iterative image reconstruction methods are usually employed since they outperform the conventional FBP methods in terms of signal-to-noise ratio (SNR) and streak artifacts [8–10]. Recent developments of compressed sensing (CS) or total variation (TV) minimization-based image reconstruction methods have reduced streak artifacts to the extent that sparse-view CT would be a plausible imaging modality for some clinical applications [11, 12]. Imaging guided radiation therapy (IGRT) using a cone-beam CT (CBCT) is one of the applications of great interests . It is now widely recognized that the CS-based image reconstruction can suppress streak artifacts to the unnoticeable level in the case of simple-structured-phantom imaging with the number of views as small as several tens [14–19]. But, in human imaging in which sparsity of the images is much lower than in the phantom imaging case, the CS-based image reconstruction methods often fail in suppressing streak artifacts. There have been a few reports on streak artifact suppression techniques in the CS-based image reconstruction. Leng et al. introduced a method to suppress respiration-induced streak artifacts in four-dimensional CBCT . They used a full-view image as a prior to suppress the streak artifacts in each respiratory phase image. They also proposed a method that a full-view image be used as a prior for the constraint in the CS-based image reconstruction from highly sparse-view projection data .
In this paper, we introduce a new sparse-view image reconstruction method to further reduce streak artifacts stemming from high-intensity objects like bony structures or metal implants. We incorporate bone segmentation into the CS-based image reconstruction to prevent streak artifact formation in the soft tissue regions. We have verified the proposed method using the projection data obtained from micro-CT scanning of a contrast phantom and a laboratory rat.
ART and CS
where A is the system matrix describing the forward projection in the CT scan . In ART, the above equation is solved in an iterative way that the difference between the projection data measured in the real scan and the projection data calculated from the estimated image is back-projected on to the image estimated at the previous iteration step. ART is known to have better performance than FBP in suppressing streak artifacts in sparse-view imaging. Many variants of ART with different iteration schemes have been proposed to improve the image quality and to reduce the computation time [22–25]. In this study, we use the ordered-subset simultaneous algebraic reconstruction technique (OS-SART)  for an ART solver.
f0 : = finit;
for k = 1: 1: K(main loop)
update fk by OS-SART from the projection data g;
for l = 1: 1: 10 (TV minimization loop)
compute the steepest decent direction d of TV;
The TV-minimization step has two control parameters, the maximum step size β in the steepest descent search, the reduction factor βred of the maximum step size after each iteration of the main loop. It is commonly known that the large step size of the steepest descent makes the image look smooth, and the small one makes the image look sharp [13, 14]. In this study, we empirically choose β and βred considering that too large β makes the image weak-contrasted whilst too small β makes the image very similar to the one reconstructed by ART .
Streak-artifact-suppressed CS image reconstruction (SAS-CS)
To reduce streak artifacts stemming from high-intensity structures like bones or metal implants, we combine CS-based image reconstruction approaches with the conventional FBP. Figure 1 shows the basic idea of the proposed method.
To identify the high-intensity region that makes strong streak artifacts, we first reconstruct an image using FBP, fFBP, from the acquired sinogram gacq. The resulting image may have streak artifacts around the high-intensity structures.
Step 2 : Extracting fbone from fFBP
From fFBP , we extract the high-intensity region, denoted as fbone, by applying a thresholding technique. We manually choose the global threshold Tbone by visual inspection of the image histogram.
Step 3 : Computing gbone by forward projecting fbone
From fbone, we compute forward projection of the high-intensity region to make the sinogram data gbone that accounts for the high-intensity region only.
Step 4 :
We subtract gbone from the measured sinogram, gacq, to exclude the components stemming from the high-intensity region.
Step 5 : ]
We use the subtracted sinogram gsoft, which account for the soft tissues only, for reconstruction of soft tissue images using the CS-based method. In this step of CS-based image reconstruction, we use a uniform image of zeroes as an initial guess of the CS-based image reconstruction.
Step 6 :
After reconstructing the soft tissue image fsoft via CS, we add the high-intensity region image fbone, which has been reconstructed by FBP, to the soft tissue image to get the composite image fsum.
Step 7 : .
To further refine the CT image, we perform the CS-based iterations again on the original sinogram with the initial guess of the CT image set to fsumobtained at the last step.
The CS-based image reconstruction in step 5 and 7 solves the constrained minimization problem defined in Eq. (2). Step 5 needs two inputs and three control parameters. The two inputs are the soft tissue sinogram, gsoft, and the initial guess of the reconstructed image which is all zeroed. The three control parameters are β and βred defined in the previous section, and K the maximum number of iterations of the main loop. In step 7, we perform the CS-based image reconstruction again using the same procedure as in step 5, but with the data inputs of gacq and fsum which has been obtained in step 6.
We have performed all the CT scans using the lab-built micro CT system described in our previous work . The micro CT system consists of a micro-focus x-ray source, a rotating object holder, a CMOS flat-panel detector. The micro-focus x-ray source (L8121-01, Hamamatsu, Japan) has a fixed tungsten anode having an angle of 25° against the electron beam and a 200 μm-thick beryllium exit window. The emitted x-ray beam has a span angle of 43°. The source has a variable focal spot size from 5 μm to 50 μm depending on the applied tube power. We have operated the micro-focus x-ray source in a continuous mode with a 1 mm-thick Al filter. We have used a commercially available flat-panel detector (C7942, Hamamatsu, Japan) as a 2D digital x-ray imager in the micro-CT system. The flat-panel detector consists of a 2240 × 2240 active matrix of transistors and photodiodes with a pixel pitch of 50 μm, and a CsI:Tl scintillator.
Image quality evaluation
where f is the sparse-view image reconstructed by the proposed method, fref the reference full-view image reconstructed by FBP, and fFBP the sparse-view image reconstructed by FBP, respectively. We calculate TV using Eq. (3). For the reference images, we use the 900-view images reconstucted by FBP which have little streak artifacts.
where f is the matrix form of the image vector f.
Means and standard deviations, RRMEs and SIs in the contrast phantom images
Contrast Phantom Case
Mean values ± standard deviations
0.1787 ± 0.0172
0.1784 ± 0.0033
0.1783 ± 0.0022
RRMEs and SIs in the rat abdomen images and rat pelvic floor images
Rat abdomen case
Rat pelvic floor case
In sparse-view imaging, reducing the extensive computation time of the CS-based image reconstruction is a great technical challenge for its application to clinical practice . Repetitive forward and backward projections account for most of the computations in the CS-based image reconstruction. Recent innovations in fast iterative image reconstructions based on graphic processing units (GPUs) have shown that sparse-view imaging may gain clinical applications in the near future [29, 30]. The computing cost of the CS-based image reconstruction is known to be higher than that of ART, which largely depends on the number of iterations to solve the minimization problem. Therefore, convergence speed of the CS-based image reconstruction is a crucial factor for its use in clinical practice. Recent development of a fast CS-based reconstruction algorithm based on the Barzilai-Borwein formulation has reduced the number of iterations to the extent that the CS-based image reconstruction could be used for real-time IGRT .
The computing cost of the proposed method, so called SAS-CS, has been found to be similar to that of the conventional CS-based reconstruction in that SAS-CS needs similar number of iterations. In addition to the repetitive forward and backward projections, SAS-CS needs additional non-iterative computations for the bone component subtraction from the measured sinogram. But, the computing cost for the bone segmentation is minimal as compared to that of the iterative computations. In fact, we have observed that SAS-CS slightly accelerates the convergence of the minimization. It seems that excluding the bone components from the measured projection data in the first iteration (step 5) of SAS-CS accounts for the convergence acceleration. In the second iteration (step 7) in which the bone components are also taken into account, the number of iterations similar to the one in step 5 suffice for further refinement of the reconstructed image in most cases. But we still need to speed up the computation for practical use of the proposed method. Recently developed fast algorithms, such as the adaptive-steepest-descent projection onto convex sets (ASD-POCS) [14, 19] or the Barzilai-Borwein formulation-based algorithm , may be used for our future studies to reduce the computation time.
In conclusion, the proposed method can suppress streak artifacts stemming from high-intensity objects in sparse-view CT imaging without significant increase of computing cost as compared to CS- or ART-based reconstructions. Experimental results obtained from the micro-CT imaging of a laboratory rat have demonstrated efficacy of the proposed method in suppressing bone-induced streak artifacts in sparse-view CT imaging.
This work was supported by the National Research Foundation (NRF) of Korea funded by the Korea government (MEST) (No. 2009-0078310).
- Bian JG, Siewerdsen JH, Han XA, Sidky EY, Prince JL, Pelizzari CA, Pan XC: Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT. Phys Med Biol 2010, 55: 6575–6599. 10.1088/0031-9155/55/22/001View ArticleGoogle Scholar
- Ritschl L, Bergner F, Fleischmann C, Kachelriess M: Improved total variation-based CT image reconstruction applied to clinical data. Phys Med Biol 2011, 56: 1545–1561. 10.1088/0031-9155/56/6/003View ArticleGoogle Scholar
- Li XL, Luo SQ, Talbert AJ, Eisner RL, DiBianca FA: A compressed sensing-based iterative algorithm for CT reconstruction and its possible application to phase contrast imaging. Biomed Eng Online 2011., 10: Google Scholar
- Brooks RA, Glover G, Talbert AJ, Eisner RL, DiBianca FA: Aliasing: a source of streaks in computed tomograms. J Comput Assist Tomo 1979, 3: 511–518. 10.1097/00004728-197908000-00014View ArticleGoogle Scholar
- Crawford CR, Kak AC: Aliasing artifacts in computerized tomography. Appl Optics 1979, 18: 3704–3711. 10.1364/AO.18.003704View ArticleGoogle Scholar
- Barrett JF, Keat N: Artifacts in CT: recognition and avoidance. Radiographics 2004, 24: 1679–1691. 10.1148/rg.246045065View ArticleGoogle Scholar
- Joseph PM, Schulz RA: View sampling requirements in fan beam computed-tomography. Med Phys 1980, 7: 692–702. 10.1118/1.594723View ArticleGoogle Scholar
- Bruyant PP, Sau J, Mallet JJ: Streak artifact reduction in filtered backprojection using a level line-Based interpolation method. J Nucl Med 2000, 41: 1913–1919.Google Scholar
- Wang G, Snyder DL, OSullivan JA, Vannier MW: Iterative deblurring for CT metal artifact reduction. IEEE T Med Imaging 1996, 15: 657–664. 10.1109/42.538943View ArticleGoogle Scholar
- De Man B, Nuyts J, Dupont P, Marchal G, Suetens P: Reduction of metal streak artifacts in x-ray computed tomography using a transmission maximum a posteriori algorithm. IEEE T Nucl Sci 2000, 47: 977–981. 10.1109/23.856534View ArticleGoogle Scholar
- Pan XC, Sidky EY, Vannier M: Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? Inverse Probl 2009, 25: 123009. 10.1088/0266-5611/25/12/123009MathSciNetView ArticleGoogle Scholar
- Tang J, Nett BE, Chen GH: Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms. Phys Med Biol 2009, 54: 5781–5804. 10.1088/0031-9155/54/19/008View ArticleGoogle Scholar
- Park JC, Song BY, Kim JS, Park SH, Kim HK, Liu ZW, Suh TS, Song WY: Fast compressed sensing-based CBCT reconstruction using Barzilai-Borwein formulation for application to on-line IGRT. Med Phys 2012, 39: 1207–1217. 10.1118/1.3679865View ArticleGoogle Scholar
- Sidky EY, Kao CM, Pan XH: Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT. J X-Ray Sci Technol 2006, 14: 119–139.Google Scholar
- Chen GH, Tang J, Leng SH: Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. Med Phys 2008, 35: 660–663. 10.1118/1.2836423View ArticleGoogle Scholar
- Yu HY, Wang G: A soft-threshold filtering approach for reconstruction from a limited number of projections. Phys Med Biol 2010, 55: 3905–3916. 10.1088/0031-9155/55/13/022View ArticleGoogle Scholar
- Song J, Liu QH, Johnson GA, Badea CT: Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT. Med Phys 2007, 34: 4476–4483. 10.1118/1.2795830View ArticleGoogle Scholar
- Zhang YH, Chan HP, Sahiner B, Wei J, Goodsitt MM, Hadjiiski LM, Ge J, Zhou CA: A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis. Med Phys 2006, 33: 3781–3795. 10.1118/1.2237543View ArticleGoogle Scholar
- Sidky EY, Pan XC: Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization. Phys Med Biol 2008, 53: 4777–4807. 10.1088/0031-9155/53/17/021View ArticleGoogle Scholar
- Leng S, Zambelli J, Tolakanahalli R, Nett B, Munro P, Star-Lack J, Paliwal B, Chena GH: Streaking artifacts reduction in four-dimensional cone-beam computed tomography. Med Phys 2008, 35: 4649–4659. 10.1118/1.2977736View ArticleGoogle Scholar
- Gordon R, Bender R, Herman GT: Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. J Theor Biol 1970, 29: 471–481. 10.1016/0022-5193(70)90109-8View ArticleGoogle Scholar
- Andersen AH, Kak AC: Simultaneous algebraic reconstruction technique (SART) - a superior implementation of the ART algorithm. Ultrasonic Imaging 1984, 6: 81–94. 10.1016/0161-7346(84)90008-7View ArticleGoogle Scholar
- Jiang M, Wang G: Convergence of the simultaneous algebraic reconstruction technique (SART). IEEE T Image Process 2003, 12: 957–961. 10.1109/TIP.2003.815295View ArticleGoogle Scholar
- Wang G, Jiang M: Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART). J X-Ray Sci Technol 2004, 12: 169–177.Google Scholar
- Gregor J, Benson T: Computational analysis and improvement of SIRT. IEEE T Med Imaging 2008, 27: 918–924.View ArticleGoogle Scholar
- Donoho DL: Compressed Sensing. IEEE T Inform Theory 2006, 52: 1289–1306.MathSciNetView ArticleGoogle Scholar
- Yu H, Wang G: Compressed sensing based interior tomography. Phys Med Biol 2009, 54: 2791–2805. 10.1088/0031-9155/54/9/014View ArticleGoogle Scholar
- Lee SC, Kim HK, Chun IK, Cho MH, Lee SY: A flat-panel detector based micro-CT system: performance evaluation for small-animal imaging. Phys Med Biol 2003, 48: 4173–4185. 10.1088/0031-9155/48/24/014View ArticleGoogle Scholar
- Xu F, Mueller K: Real-time 3D computed tomographic reconstruction using commodity graphics hardware. Phys Med Biol 2007, 52: 3405–3419. 10.1088/0031-9155/52/12/006View ArticleGoogle Scholar
- Jia X, Lou YF, Li RJ, Song WY, Jiang SB: GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation. Med Phys 2010, 37: 1757–1760. 10.1118/1.3371691View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.