Material-separating regularizer for multi-energy x-ray tomography

Dual-energy X-ray tomography is considered in a context where the target under imaging consists of two distinct materials. The materials are assumed to be possibly intertwined space, but at any given location there only one material present. Further, energies chosen so that clear difference spectral dependence attenuation coefficients A novel regularizer presented for inverse problem reconstructing separate tomographic images combination things, (a) non-negativity constraint, and (b) penalty term containing inner product between images, promotes presence most pixel. preconditioned interior point method derived minimization regularization functional. Numerical tests with digital phantoms suggest new algorithm outperforms baseline method, Joint Total Variation regularization, terms correctly material-characterized pixels. While tested two-dimensional setting energies, approach readily generalizes three dimensions more number just needs match used imaging.