X-ray CT image reconstruction via wavelet frame based regularization and Radon domain inpainting

Bin Dong, Jia Li, Zuowei Shen

Research output: Contribution to journalArticle

45 Citations (Scopus)

Abstract

X-ray computed tomography (CT) has been playing an important role in diagnostic of cancer and radiotherapy. However, high imaging dose added to healthy organs during CT scans is a serious clinical concern. Imaging dose in CT scans can be reduced by reducing the number of X-ray projections. In this paper, we consider 2D CT reconstructions using very small number of projections. Some regularization based reconstruction methods have already been proposed in the literature for such task, like the total variation (TV) based reconstruction (Sidky and Pan in Phys. Med. Biol. 53:4777, 2008; Sidky et al. in J. X-Ray Sci. Technol. 14(2):119-139, 2006; Jia et al. in Med. Phys. 37:1757, 2010; Choi et al. in Med. Phys. 37:5113, 2010) and balanced approach with wavelet frame based regularization (Jia et al. in Phys. Med. Biol. 56:3787-3807, 2011). For most of the existing methods, at least 40 projections is usually needed to get a satisfactory reconstruction. In order to keep radiation dose as minimal as possible, while increase the quality of the reconstructed images, one needs to enhance the resolution of the projected image in the Radon domain without increasing the total number of projections. The goal of this paper is to propose a CT reconstruction model with wavelet frame based regularization and Radon domain inpainting. The proposed model simultaneously reconstructs a high quality image and its corresponding high resolution measurements in Radon domain. In addition, we discovered that using the isotropic wavelet frame regularization proposed in Cai et al. (Image restorations: total variation, wavelet frames and beyond, 2011, preprint) is superior than using its anisotropic counterpart. Our proposed model, as well as other models presented in this paper, is solved rather efficiently by split Bregman algorithm (Goldstein and Osher in SIAM J. Imaging Sci. 2(2):323-343, 2009; Cai et al. in Multiscale Model. Simul. 8(2):337-369, 2009). Numerical simulations and comparisons will be presented at the end.

Original languageEnglish (US)
Pages (from-to)333-349
Number of pages17
JournalJournal of Scientific Computing
Volume54
Issue number2-3
DOIs
StatePublished - Feb 2013

Fingerprint

Inpainting
X-ray Tomography
Wavelet Frames
Radon
Computed Tomography
Image Reconstruction
Image reconstruction
Tomography
Regularization
X rays
Projection
Dose
Imaging
Total Variation
Imaging techniques
Dosimetry
Radiotherapy
Multiscale Model
Image Restoration
Numerical Comparisons

Keywords

  • Computed tomography
  • Radon domain inpainting
  • Split Bregman algorithm
  • Total variation
  • Wavelet frame

ASJC Scopus subject areas

  • Software
  • Computational Theory and Mathematics
  • Theoretical Computer Science
  • Engineering(all)

Cite this

X-ray CT image reconstruction via wavelet frame based regularization and Radon domain inpainting. / Dong, Bin; Li, Jia; Shen, Zuowei.

In: Journal of Scientific Computing, Vol. 54, No. 2-3, 02.2013, p. 333-349.

Research output: Contribution to journalArticle

@article{936c1853684741e28e65ca365a8a0a13,
title = "X-ray CT image reconstruction via wavelet frame based regularization and Radon domain inpainting",
abstract = "X-ray computed tomography (CT) has been playing an important role in diagnostic of cancer and radiotherapy. However, high imaging dose added to healthy organs during CT scans is a serious clinical concern. Imaging dose in CT scans can be reduced by reducing the number of X-ray projections. In this paper, we consider 2D CT reconstructions using very small number of projections. Some regularization based reconstruction methods have already been proposed in the literature for such task, like the total variation (TV) based reconstruction (Sidky and Pan in Phys. Med. Biol. 53:4777, 2008; Sidky et al. in J. X-Ray Sci. Technol. 14(2):119-139, 2006; Jia et al. in Med. Phys. 37:1757, 2010; Choi et al. in Med. Phys. 37:5113, 2010) and balanced approach with wavelet frame based regularization (Jia et al. in Phys. Med. Biol. 56:3787-3807, 2011). For most of the existing methods, at least 40 projections is usually needed to get a satisfactory reconstruction. In order to keep radiation dose as minimal as possible, while increase the quality of the reconstructed images, one needs to enhance the resolution of the projected image in the Radon domain without increasing the total number of projections. The goal of this paper is to propose a CT reconstruction model with wavelet frame based regularization and Radon domain inpainting. The proposed model simultaneously reconstructs a high quality image and its corresponding high resolution measurements in Radon domain. In addition, we discovered that using the isotropic wavelet frame regularization proposed in Cai et al. (Image restorations: total variation, wavelet frames and beyond, 2011, preprint) is superior than using its anisotropic counterpart. Our proposed model, as well as other models presented in this paper, is solved rather efficiently by split Bregman algorithm (Goldstein and Osher in SIAM J. Imaging Sci. 2(2):323-343, 2009; Cai et al. in Multiscale Model. Simul. 8(2):337-369, 2009). Numerical simulations and comparisons will be presented at the end.",
keywords = "Computed tomography, Radon domain inpainting, Split Bregman algorithm, Total variation, Wavelet frame",
author = "Bin Dong and Jia Li and Zuowei Shen",
year = "2013",
month = "2",
doi = "10.1007/s10915-012-9579-6",
language = "English (US)",
volume = "54",
pages = "333--349",
journal = "Journal of Scientific Computing",
issn = "0885-7474",
publisher = "Springer New York",
number = "2-3",

}

TY - JOUR

T1 - X-ray CT image reconstruction via wavelet frame based regularization and Radon domain inpainting

AU - Dong, Bin

AU - Li, Jia

AU - Shen, Zuowei

PY - 2013/2

Y1 - 2013/2

N2 - X-ray computed tomography (CT) has been playing an important role in diagnostic of cancer and radiotherapy. However, high imaging dose added to healthy organs during CT scans is a serious clinical concern. Imaging dose in CT scans can be reduced by reducing the number of X-ray projections. In this paper, we consider 2D CT reconstructions using very small number of projections. Some regularization based reconstruction methods have already been proposed in the literature for such task, like the total variation (TV) based reconstruction (Sidky and Pan in Phys. Med. Biol. 53:4777, 2008; Sidky et al. in J. X-Ray Sci. Technol. 14(2):119-139, 2006; Jia et al. in Med. Phys. 37:1757, 2010; Choi et al. in Med. Phys. 37:5113, 2010) and balanced approach with wavelet frame based regularization (Jia et al. in Phys. Med. Biol. 56:3787-3807, 2011). For most of the existing methods, at least 40 projections is usually needed to get a satisfactory reconstruction. In order to keep radiation dose as minimal as possible, while increase the quality of the reconstructed images, one needs to enhance the resolution of the projected image in the Radon domain without increasing the total number of projections. The goal of this paper is to propose a CT reconstruction model with wavelet frame based regularization and Radon domain inpainting. The proposed model simultaneously reconstructs a high quality image and its corresponding high resolution measurements in Radon domain. In addition, we discovered that using the isotropic wavelet frame regularization proposed in Cai et al. (Image restorations: total variation, wavelet frames and beyond, 2011, preprint) is superior than using its anisotropic counterpart. Our proposed model, as well as other models presented in this paper, is solved rather efficiently by split Bregman algorithm (Goldstein and Osher in SIAM J. Imaging Sci. 2(2):323-343, 2009; Cai et al. in Multiscale Model. Simul. 8(2):337-369, 2009). Numerical simulations and comparisons will be presented at the end.

AB - X-ray computed tomography (CT) has been playing an important role in diagnostic of cancer and radiotherapy. However, high imaging dose added to healthy organs during CT scans is a serious clinical concern. Imaging dose in CT scans can be reduced by reducing the number of X-ray projections. In this paper, we consider 2D CT reconstructions using very small number of projections. Some regularization based reconstruction methods have already been proposed in the literature for such task, like the total variation (TV) based reconstruction (Sidky and Pan in Phys. Med. Biol. 53:4777, 2008; Sidky et al. in J. X-Ray Sci. Technol. 14(2):119-139, 2006; Jia et al. in Med. Phys. 37:1757, 2010; Choi et al. in Med. Phys. 37:5113, 2010) and balanced approach with wavelet frame based regularization (Jia et al. in Phys. Med. Biol. 56:3787-3807, 2011). For most of the existing methods, at least 40 projections is usually needed to get a satisfactory reconstruction. In order to keep radiation dose as minimal as possible, while increase the quality of the reconstructed images, one needs to enhance the resolution of the projected image in the Radon domain without increasing the total number of projections. The goal of this paper is to propose a CT reconstruction model with wavelet frame based regularization and Radon domain inpainting. The proposed model simultaneously reconstructs a high quality image and its corresponding high resolution measurements in Radon domain. In addition, we discovered that using the isotropic wavelet frame regularization proposed in Cai et al. (Image restorations: total variation, wavelet frames and beyond, 2011, preprint) is superior than using its anisotropic counterpart. Our proposed model, as well as other models presented in this paper, is solved rather efficiently by split Bregman algorithm (Goldstein and Osher in SIAM J. Imaging Sci. 2(2):323-343, 2009; Cai et al. in Multiscale Model. Simul. 8(2):337-369, 2009). Numerical simulations and comparisons will be presented at the end.

KW - Computed tomography

KW - Radon domain inpainting

KW - Split Bregman algorithm

KW - Total variation

KW - Wavelet frame

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

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

U2 - 10.1007/s10915-012-9579-6

DO - 10.1007/s10915-012-9579-6

M3 - Article

AN - SCOPUS:84885445816

VL - 54

SP - 333

EP - 349

JO - Journal of Scientific Computing

JF - Journal of Scientific Computing

SN - 0885-7474

IS - 2-3

ER -