Robust Beam Hardening Artifacts Reduction for Computed Tomography Using Spectrum Modeling

The aim of this paper is to develop a fast and accurate beam hardening correction method by modeling physical interactions between X-ray photons and materials for computed tomography (CT) imaging. The nonlinear attenuation process of the X-ray projection is modeled by reprojecting a template image with the estimated polychromatic spectrum. By adding the scaled difference of the monochromatic reprojection data and the polychromatic reprojection to the raw projection data, the raw projection data is mapped into the corresponding monochromatic projection data, which is to reconstruct the beam hardening artifacts corrected images. The algorithm can also be implemented in image-domain which takes the uncorrected image volume as input when there is an adequate model of the spectrum. In this case, the scaled difference is reconstructed to yield a set of artifacts images that can be added directly to the uncorrected images. Numerical simulations, experimental phantom data, and animal data which are acquired on a modern diagnostic CT scanner (Discovery CT750 HD, GE Healthcare, WI, USA), and a modern C-Arm CT scanner (Artis Zee, Siemens Healthcare, Forchheim, Germany), respectively, are used to evaluate the proposed method. The results show the proposed method significantly reduced both cupping and streak artifacts, and successfully recovered the Housfield units accuracy. Extensive studies suggest the proposed model-based method successfully corrects the beam hardening artifacts. This paper is practically useful and is promising to be applied to commercial products.