Narrow-Energy-Width CT Based on Multivoltage X-Ray Image Decomposition

Statistical iterative spectral CT imaging method based on blind separation of polychromatic projections

Spectral computed tomography (CT) can provide narrow-energy-width reconstructed images, thereby suppressing beam hardening artifacts and providing rich attenuation information for component characterization. We propose a statistical iterative spectral CT imaging method based on blind separation of polychromatic projections to improve the accuracy of narrow-energy-width image decomposition. For direct inversion in blind scenarios, we introduce the system matrix into the X-ray multispectral forward model to reduce indirect errors. A constrained optimization problem with edge-preserving regularization is established and decomposed into two sub-problems to be alternately solved. Experiments indicate that the novel algorithm obtains more accurate narrow-energy-width images than the state-of-the-art method.

Spectral CT imaging method based on blind separation of polychromatic projections with Poisson prior

The linear reconstruction of narrow-energy-width projections can suppress hardening artifacts in conventional computed tomography (CT). We develop a spectral CT blind separation algorithm for obtaining narrow-energy-width projections under a blind scenario where the incident spectra are unknown. The algorithm relies on an X-ray multispectral forward model. Based on the Poisson statistical properties of measurements, a constrained optimization problem is established and solved by a block coordinate descent algorithm that alternates between nonnegative matrix factorization and Gauss-Newton algorithm. Experiments indicate that the decomposed projections conform to the characteristics of narrow-energy-width projections. The new algorithm improves the accuracy of obtaining narrow-energy-width projections.