Document Type

Conference Proceeding

Publication Date

2014

Journal or Book Title

AIP Conference Proceedings

Volume

1581

First Page

1848

Last Page

1855

DOI

10.1063/1.4865048

Abstract

We propose a method for reconstructing sparse images from polychromatic x-ray computed tomography (ct) measurements via mass attenuation coefficient discretization. The material of the inspected object and the incident spectrum are assumed to be unknown. We rewrite the Lambert-Beer’s law in terms of integral expressions of mass attenuation and discretize the resulting integrals. We then present a penalized constrained least-squares optimization approach forreconstructing the underlying object from log-domain measurements, where an active set approach is employed to estimate incident energy density parameters and the nonnegativity and sparsity of the image density map are imposed using negative-energy and smooth ℓ1-norm penalty terms. We propose a two-step scheme for refining the mass attenuation discretization grid by using higher sampling rate over the range with higher photon energy, and eliminating the discretization points that have little effect on accuracy of the forward projection model. This refinement allows us to successfully handle the characteristic lines (Dirac impulses) in the incident energy density spectrum. We compare the proposed method with the standard filtered backprojection, which ignores the polychromatic nature of the measurements and sparsity of theimage density map. Numerical simulations using both realistic simulated and real x-ray ct data are presented.

Comments

The following article appeared in AIP Conference Proceedings 1581 (2014): 1848 and may be found at doi:10.1063/1.4865048.

Rights

Copyright 2014 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics.

Copyright Owner

AIP Publishing LLC

Language

en

File Format

application/pdf

Share

COinS