Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through Radon inversion initialization

Wojciech Zbijewski, Freek J. Beekman

Research output: Contribution to journalArticle

Abstract

Statistical reconstruction (SR) methods provide a general and flexible framework for obtaining tomographic images from projections. For several applications SR has been shown to outperform analytical algorithms in terms of resolution -noise trade-off achieved in the reconstructions. A disadvantage of SR is the long computational time required to obtain the reconstructions, in particular when large data sets characteristic for x-ray computer tomography (CT) are involved. As was shown recently, by combining statistical methods with block iterative acceleration schemes [e.g., like in the ordered subsets convex (OSC) algorithm], the reconstruction time for x-ray CT applications can be reduced by about two orders of magnitude. There are, however, some factors lengthening the reconstruction process that hamper both accelerated and standard statistical algorithms to similar degree. In this simulation study based on monoenergetic and scatter-free projection data, we demonstrate that one of these factors is the extremely high number of iterations needed to remove artifacts that can appear around high-contrast structures. We also show (using the OSC method) that these artifacts can be adequately suppressed if statistical reconstruction is initialized with images generated by means of Radon inversion algorithms like filtered back projection (FBP). This allows the reconstruction time to be shortened by even as much as one order of magnitude. Although the initialization of the statistical algorithm with FBP image introduces some additional noise into the first iteration of OSC reconstruction, the resolution -noise trade-off and the contrast-to-noise ratio of final images are not markedly compromised.

Original languageEnglish (US)
Pages (from-to)62-69
Number of pages8
JournalMedical Physics
Volume31
Issue number1
DOIs
StatePublished - Jan 2004
Externally publishedYes

Fingerprint

Radon
Artifacts
Tomography
X-Rays
Noise

Keywords

  • FBP
  • Imaging
  • Iterative reconstruction
  • Maximum likelihood
  • Ordered subsets
  • Tomography

ASJC Scopus subject areas

  • Biophysics

Cite this

@article{ab39395226294e0e9ac0602511e822c3,
title = "Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through Radon inversion initialization",
abstract = "Statistical reconstruction (SR) methods provide a general and flexible framework for obtaining tomographic images from projections. For several applications SR has been shown to outperform analytical algorithms in terms of resolution -noise trade-off achieved in the reconstructions. A disadvantage of SR is the long computational time required to obtain the reconstructions, in particular when large data sets characteristic for x-ray computer tomography (CT) are involved. As was shown recently, by combining statistical methods with block iterative acceleration schemes [e.g., like in the ordered subsets convex (OSC) algorithm], the reconstruction time for x-ray CT applications can be reduced by about two orders of magnitude. There are, however, some factors lengthening the reconstruction process that hamper both accelerated and standard statistical algorithms to similar degree. In this simulation study based on monoenergetic and scatter-free projection data, we demonstrate that one of these factors is the extremely high number of iterations needed to remove artifacts that can appear around high-contrast structures. We also show (using the OSC method) that these artifacts can be adequately suppressed if statistical reconstruction is initialized with images generated by means of Radon inversion algorithms like filtered back projection (FBP). This allows the reconstruction time to be shortened by even as much as one order of magnitude. Although the initialization of the statistical algorithm with FBP image introduces some additional noise into the first iteration of OSC reconstruction, the resolution -noise trade-off and the contrast-to-noise ratio of final images are not markedly compromised.",
keywords = "FBP, Imaging, Iterative reconstruction, Maximum likelihood, Ordered subsets, Tomography",
author = "Wojciech Zbijewski and Beekman, {Freek J.}",
year = "2004",
month = "1",
doi = "10.1118/1.1631091",
language = "English (US)",
volume = "31",
pages = "62--69",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "1",

}

TY - JOUR

T1 - Suppression of intensity transition artifacts in statistical x-ray computer tomography reconstruction through Radon inversion initialization

AU - Zbijewski, Wojciech

AU - Beekman, Freek J.

PY - 2004/1

Y1 - 2004/1

N2 - Statistical reconstruction (SR) methods provide a general and flexible framework for obtaining tomographic images from projections. For several applications SR has been shown to outperform analytical algorithms in terms of resolution -noise trade-off achieved in the reconstructions. A disadvantage of SR is the long computational time required to obtain the reconstructions, in particular when large data sets characteristic for x-ray computer tomography (CT) are involved. As was shown recently, by combining statistical methods with block iterative acceleration schemes [e.g., like in the ordered subsets convex (OSC) algorithm], the reconstruction time for x-ray CT applications can be reduced by about two orders of magnitude. There are, however, some factors lengthening the reconstruction process that hamper both accelerated and standard statistical algorithms to similar degree. In this simulation study based on monoenergetic and scatter-free projection data, we demonstrate that one of these factors is the extremely high number of iterations needed to remove artifacts that can appear around high-contrast structures. We also show (using the OSC method) that these artifacts can be adequately suppressed if statistical reconstruction is initialized with images generated by means of Radon inversion algorithms like filtered back projection (FBP). This allows the reconstruction time to be shortened by even as much as one order of magnitude. Although the initialization of the statistical algorithm with FBP image introduces some additional noise into the first iteration of OSC reconstruction, the resolution -noise trade-off and the contrast-to-noise ratio of final images are not markedly compromised.

AB - Statistical reconstruction (SR) methods provide a general and flexible framework for obtaining tomographic images from projections. For several applications SR has been shown to outperform analytical algorithms in terms of resolution -noise trade-off achieved in the reconstructions. A disadvantage of SR is the long computational time required to obtain the reconstructions, in particular when large data sets characteristic for x-ray computer tomography (CT) are involved. As was shown recently, by combining statistical methods with block iterative acceleration schemes [e.g., like in the ordered subsets convex (OSC) algorithm], the reconstruction time for x-ray CT applications can be reduced by about two orders of magnitude. There are, however, some factors lengthening the reconstruction process that hamper both accelerated and standard statistical algorithms to similar degree. In this simulation study based on monoenergetic and scatter-free projection data, we demonstrate that one of these factors is the extremely high number of iterations needed to remove artifacts that can appear around high-contrast structures. We also show (using the OSC method) that these artifacts can be adequately suppressed if statistical reconstruction is initialized with images generated by means of Radon inversion algorithms like filtered back projection (FBP). This allows the reconstruction time to be shortened by even as much as one order of magnitude. Although the initialization of the statistical algorithm with FBP image introduces some additional noise into the first iteration of OSC reconstruction, the resolution -noise trade-off and the contrast-to-noise ratio of final images are not markedly compromised.

KW - FBP

KW - Imaging

KW - Iterative reconstruction

KW - Maximum likelihood

KW - Ordered subsets

KW - Tomography

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

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

U2 - 10.1118/1.1631091

DO - 10.1118/1.1631091

M3 - Article

VL - 31

SP - 62

EP - 69

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 1

ER -