- Open Access
Bone-induced streak artifact suppression in sparse-view CT image reconstruction
BioMedical Engineering OnLine volume 11, Article number: 44 (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
We use the algebraic reconstruction technique (ART) and CS as a platform for the sparse-view image reconstruction. ART is a minimum mean square error (MMSE) solver to find the image f that best matches the measured projection data g:
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.
The CS-based image reconstruction methods solve the following constrained optimization problem which has constraints of data fidelity and pixel positivity :
where Ψ is a sparsifying transform operator, and is the l 1 norm of an N-dimensional vector z. In this study, we use the discrete gradient transform for the sparsifying transform which has been widely used in the CS-based image reconstruction :
where i and j are the pixel indices in the x- and y-directions, respectively, and f is the 2D matrix form of f. The discrete gradient transform is often denoted as TV. We implement the CS-based image reconstruction algorithm using the OS-SART to enforce the data fidelity and the steepest descent method to minimize the TV in an alternating manner in the iteration. We summarize the algorithm for the CS-based image reconstruction in the following pseudo-code [19, 27].
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.
Step 1 :
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.
To validate the proposed method, we have performed CT scans of a contrast phantom and a sacrificed adult rat using the micro-CT. The contrast phantom consists of seven inserts six of which have physical densities similar to that of water and the rest of which has the bone-equivalent physical density. Figure 2 shows a schematic diagram of the contrast phantom along with the physical densities of the inserts. The seven inserts of 5 mm diameter were in a water bath made of an acryl cylinder of 40 mm diameter. We have made the inserts using the commercial electron density phantoms (Model 76-430, Nuclear Associates, NY, USA). We applied tube voltage and current of 40 kVp and 0.5 mA for the contrast phantom imaging, and 65 kVp and 0.34 mA for the rat imaging, respectively. To get reference images, we performed full-view scans with the number of views of 900 over 360 degrees. To get sparse-view projection data, we decimated the full-view projection data in the view direction.
Image quality evaluation
In order to evaluate the final image quality particularly in terms of streak artifact formation, we use two metrics, one for the total variations stemming from streak artifacts and the other for the relative errors of the final image with respect to the reference image. We define the normalized streak indicator (SI) using the total variation of the difference image  :
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.
To evaluate reconstruction errors as compared to the reference image, we use the relative root mean square error (RRME) defined by :
where f is the matrix form of the image vector f.
We have first reconstructed contrast phantom images the size of 512 × 512 using FBP from the 60-, and 900-view projection data acquired from the micro-CT scan. The images shown in Figure 3 have the pixel size of 85 × 85 μm2. As can be noticed from Figure 3b, the 900-view image to be used as a reference image shows little streak artifacts whilst the 60-view image shows strong streak artifacts stemming from the high-intensity insert placed at the center of the phantom. We have reconstructed the contrast phantom images using the three iterative reconstruction methods, ART (Figure 4a), CS (Figure 4b) and SAS-CS (Figure 4d). Figure 4c shows the soft-intensity image calculated in step 5 in SAS-CS. From the images, we can see that the ART image has strongest streak artifacts. In the CS image, the streak artifact has been remarkably reduced, but we can still see the residual streak artifacts. The SAS-CS image shows least level of steak artifacts among the three images. Figure 5 shows pixel-intensity profiles along the two lines shown in Figure 4a. Along the line 1 which runs the uniform background region, SAS-CS shows least level of fluctuation demonstrating its performance of streak-artifact suppression. Along the line 2 which runs over the inserts, CS and SAS-CS show similar level of fluctuations. Table 1 summarizes the RRMEs and SIs of the images reconstructed by the aforementioned three methods with respect to the reference image. Due to the strong streak artifacts, the ART image shows the biggest RRMEs whilst the SAS-CS image shows the least RRMEs. In terms of the normalized SI, the ART image also shows the worst performance and the SAS-CS image shows the best performance.
We have reconstructed rat abdomen images and pelvic floor images, with the matrix size of 512 × 512 and the pixel size of 120 × 120 μm2, using FBP from the 100-, and 900-view projection data. Figure 6 shows the rat abdomen images reconstructed by FBP. Here again, the 900-view FBP image is used as a reference image. The streak artifacts in the rat abdomen images are less conspicuous than they are in the contrast phantom images, but we can see clear streak artifacts from the 100-view images. Figure 7 shows the rat abdomen images reconstructed by the three iterative methods. Here again, the ART image (Figure 7a) shows the strongest streak artifacts and the SAS-CS image (Figure 7c) show less streak artifacts than the CS image (Figure 7b). Figure 8 shows the rat pelvic floor images reconstructed by FBP from the 100- and 900-view projection data. Due to the many bones on the pelvic floor, the 100-view FBP image shows strong streak artifacts. Figure 9 shows the pelvic floor images reconstructed by the iterative methods and the difference images taken from the reference image. Due to the high-intensity bones on the imaging plane, SAS-CS images (Figure 9c) also show residual streak artifacts which are, however, far less than the ART image (Figure 9a) and the CS image (Figure 9b). Table 2 summarizes the RRMEs and SIs of the rat pelvic images reconstructed by the aforementioned three methods with respect to the reference image. Due to the strong streak artifacts, the ART image shows the biggest RRME whilst the SAS-CS image shows the least RRMEs. In terms of the normalized SI, the ART image also shows the worst performance and the SAS-CS image shows the best performance.
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.
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/001
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/003
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:
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-00014
Crawford CR, Kak AC: Aliasing artifacts in computerized tomography. Appl Optics 1979, 18: 3704–3711. 10.1364/AO.18.003704
Barrett JF, Keat N: Artifacts in CT: recognition and avoidance. Radiographics 2004, 24: 1679–1691. 10.1148/rg.246045065
Joseph PM, Schulz RA: View sampling requirements in fan beam computed-tomography. Med Phys 1980, 7: 692–702. 10.1118/1.594723
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.
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.538943
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.856534
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/123009
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/008
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.3679865
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.
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.2836423
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/022
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.2795830
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.2237543
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/021
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.2977736
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-8
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-7
Jiang M, Wang G: Convergence of the simultaneous algebraic reconstruction technique (SART). IEEE T Image Process 2003, 12: 957–961. 10.1109/TIP.2003.815295
Wang G, Jiang M: Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART). J X-Ray Sci Technol 2004, 12: 169–177.
Gregor J, Benson T: Computational analysis and improvement of SIRT. IEEE T Med Imaging 2008, 27: 918–924.
Donoho DL: Compressed Sensing. IEEE T Inform Theory 2006, 52: 1289–1306.
Yu H, Wang G: Compressed sensing based interior tomography. Phys Med Biol 2009, 54: 2791–2805. 10.1088/0031-9155/54/9/014
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/014
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/006
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.3371691
This work was supported by the National Research Foundation (NRF) of Korea funded by the Korea government (MEST) (No. 2009-0078310).
The authors declare that they have no competing interests.
SO and OK conceived the study, implemented the proposed idea, evaluated the reconstructed images, and drafted the manuscript. JG and SY designed and performed the experiments, advised on the image quality evaluation. SY and OK supervised the study. All authors have read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Jin, S.O., Kim, J.G., Lee, S.Y. et al. Bone-induced streak artifact suppression in sparse-view CT image reconstruction. BioMed Eng OnLine 11, 44 (2012). https://doi.org/10.1186/1475-925X-11-44
- Streak artifact
- Sparse-view CT
- Compressed sensing
- Total variation