Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography

Jing Wang, Tianfang Li, Hongbing Lu, Zhengrong Liang

Research output: Contribution to journalArticle

317 Citations (Scopus)

Abstract

Reconstructing low-dose X-ray computed tomography (CT) images is a noise problem. This work investigated a penalized weighted least-squares (PWLS) approach to address this problem in two dimensions, where the WLS considers first-and second-order noise moments and the penalty models signal spatial correlations. Three different implementations were studied for the PWLS minimization. One utilizes a Markov random field (MRF) Gibbs functional to consider spatial correlations among nearby detector bins and projection views in sinogram space and minimizes the PWLS cost function by iterative Gauss-Seidel algorithm. Another employs Karhunen-Loève (KL) transform to de-correlate data signals among nearby views and minimizes the PWLS adaptively to each KL component by analytical calculation, where the spatial correlation among nearby bins is modeled - by the same Gibbs functional. The third one models the spatial correlations among image pixels in image domain also by a MRF Gibbs functional and minimizes the PWLS by iterative successive over-relaxation algorithm. In these three implementations, a quadratic functional regularization was chosen for the MRF model. Phantom experiments showed a comparable performance of these three PWLS-based methods in terms of suppressing noise-induced streak artifacts and preserving resolution in the reconstructed images. Computer simulations concurred with the phantom experiments in terms of noise-resolution tradeoff and detectability in low contrast environment. The KL-PWLS implementation may have the advantage in terms of computation for high-resolution dynamic low-dose CT imaging.

Original languageEnglish (US)
Pages (from-to)1272-1283
Number of pages12
JournalIEEE Transactions on Medical Imaging
Volume25
Issue number10
DOIs
StatePublished - Oct 2006

Fingerprint

X Ray Computed Tomography
Computer-Assisted Image Processing
Noise abatement
Least-Squares Analysis
Image reconstruction
Dosimetry
Tomography
Noise
Bins
X rays
Love
Cost functions
Pixels
Experiments
Detectors
Imaging techniques
Computer simulation
Computer Simulation
Artifacts
Costs and Cost Analysis

Keywords

  • Low-dose X-ray computed tomography
  • Noise reduction
  • Penalized weighted least-squares

ASJC Scopus subject areas

  • Biomedical Engineering
  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology
  • Electrical and Electronic Engineering
  • Computer Science Applications
  • Computational Theory and Mathematics

Cite this

Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography. / Wang, Jing; Li, Tianfang; Lu, Hongbing; Liang, Zhengrong.

In: IEEE Transactions on Medical Imaging, Vol. 25, No. 10, 10.2006, p. 1272-1283.

Research output: Contribution to journalArticle

@article{6d3a25b97edb4062b92568b41caf9576,
title = "Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography",
abstract = "Reconstructing low-dose X-ray computed tomography (CT) images is a noise problem. This work investigated a penalized weighted least-squares (PWLS) approach to address this problem in two dimensions, where the WLS considers first-and second-order noise moments and the penalty models signal spatial correlations. Three different implementations were studied for the PWLS minimization. One utilizes a Markov random field (MRF) Gibbs functional to consider spatial correlations among nearby detector bins and projection views in sinogram space and minimizes the PWLS cost function by iterative Gauss-Seidel algorithm. Another employs Karhunen-Lo{\`e}ve (KL) transform to de-correlate data signals among nearby views and minimizes the PWLS adaptively to each KL component by analytical calculation, where the spatial correlation among nearby bins is modeled - by the same Gibbs functional. The third one models the spatial correlations among image pixels in image domain also by a MRF Gibbs functional and minimizes the PWLS by iterative successive over-relaxation algorithm. In these three implementations, a quadratic functional regularization was chosen for the MRF model. Phantom experiments showed a comparable performance of these three PWLS-based methods in terms of suppressing noise-induced streak artifacts and preserving resolution in the reconstructed images. Computer simulations concurred with the phantom experiments in terms of noise-resolution tradeoff and detectability in low contrast environment. The KL-PWLS implementation may have the advantage in terms of computation for high-resolution dynamic low-dose CT imaging.",
keywords = "Low-dose X-ray computed tomography, Noise reduction, Penalized weighted least-squares",
author = "Jing Wang and Tianfang Li and Hongbing Lu and Zhengrong Liang",
year = "2006",
month = "10",
doi = "10.1109/TMI.2006.882141",
language = "English (US)",
volume = "25",
pages = "1272--1283",
journal = "IEEE Transactions on Medical Imaging",
issn = "0278-0062",
publisher = "Institute of Electrical and Electronics Engineers Inc.",
number = "10",

}

TY - JOUR

T1 - Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography

AU - Wang, Jing

AU - Li, Tianfang

AU - Lu, Hongbing

AU - Liang, Zhengrong

PY - 2006/10

Y1 - 2006/10

N2 - Reconstructing low-dose X-ray computed tomography (CT) images is a noise problem. This work investigated a penalized weighted least-squares (PWLS) approach to address this problem in two dimensions, where the WLS considers first-and second-order noise moments and the penalty models signal spatial correlations. Three different implementations were studied for the PWLS minimization. One utilizes a Markov random field (MRF) Gibbs functional to consider spatial correlations among nearby detector bins and projection views in sinogram space and minimizes the PWLS cost function by iterative Gauss-Seidel algorithm. Another employs Karhunen-Loève (KL) transform to de-correlate data signals among nearby views and minimizes the PWLS adaptively to each KL component by analytical calculation, where the spatial correlation among nearby bins is modeled - by the same Gibbs functional. The third one models the spatial correlations among image pixels in image domain also by a MRF Gibbs functional and minimizes the PWLS by iterative successive over-relaxation algorithm. In these three implementations, a quadratic functional regularization was chosen for the MRF model. Phantom experiments showed a comparable performance of these three PWLS-based methods in terms of suppressing noise-induced streak artifacts and preserving resolution in the reconstructed images. Computer simulations concurred with the phantom experiments in terms of noise-resolution tradeoff and detectability in low contrast environment. The KL-PWLS implementation may have the advantage in terms of computation for high-resolution dynamic low-dose CT imaging.

AB - Reconstructing low-dose X-ray computed tomography (CT) images is a noise problem. This work investigated a penalized weighted least-squares (PWLS) approach to address this problem in two dimensions, where the WLS considers first-and second-order noise moments and the penalty models signal spatial correlations. Three different implementations were studied for the PWLS minimization. One utilizes a Markov random field (MRF) Gibbs functional to consider spatial correlations among nearby detector bins and projection views in sinogram space and minimizes the PWLS cost function by iterative Gauss-Seidel algorithm. Another employs Karhunen-Loève (KL) transform to de-correlate data signals among nearby views and minimizes the PWLS adaptively to each KL component by analytical calculation, where the spatial correlation among nearby bins is modeled - by the same Gibbs functional. The third one models the spatial correlations among image pixels in image domain also by a MRF Gibbs functional and minimizes the PWLS by iterative successive over-relaxation algorithm. In these three implementations, a quadratic functional regularization was chosen for the MRF model. Phantom experiments showed a comparable performance of these three PWLS-based methods in terms of suppressing noise-induced streak artifacts and preserving resolution in the reconstructed images. Computer simulations concurred with the phantom experiments in terms of noise-resolution tradeoff and detectability in low contrast environment. The KL-PWLS implementation may have the advantage in terms of computation for high-resolution dynamic low-dose CT imaging.

KW - Low-dose X-ray computed tomography

KW - Noise reduction

KW - Penalized weighted least-squares

UR - http://www.scopus.com/inward/record.url?scp=33947260155&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33947260155&partnerID=8YFLogxK

U2 - 10.1109/TMI.2006.882141

DO - 10.1109/TMI.2006.882141

M3 - Article

C2 - 17024831

AN - SCOPUS:33947260155

VL - 25

SP - 1272

EP - 1283

JO - IEEE Transactions on Medical Imaging

JF - IEEE Transactions on Medical Imaging

SN - 0278-0062

IS - 10

ER -