A weighted polynomial based material decomposition method for spectral x-ray CT imaging

Currently in photon counting based spectral x-ray computed tomography (CT) imaging, pre-reconstruction basis materials decomposition is an effective way to reconstruct densities of various materials. The iterative maximum-likelihood method requires precise spectrum information and is time-costly. In this paper, a novel non-iterative decomposition method based on polynomials is proposed for spectral CT, whose aim was to optimize the noise performance when there is more energy bins than the number of basis materials. Several subsets were taken from all the energy bins and conventional polynomials were established for each of them. The decomposition results from each polynomial were summed with pre-calculated weighting factors, which were designed to minimize the overall noises. Numerical studies showed that the decomposition noise of the proposed method was close to the Cramer-Rao lower bound under Poisson noises. Furthermore, experiments were carried out with an XCounter Filte X1 photon counting detector for two-material decomposition and three-material decomposition for validation.