Application of an expectation maximization method to the reconstruction of X-ray-tube spectra from transmission data

Abstract

An expectation maximization method is applied to the reconstruction of X-ray tube spectra from transmission measurements in the energy range 7–40keV. A semiconductor single-photon counting detector, ionization chambers and a scintillator-based detector are used for the experimental measurement of the transmission. The number of iterations required to reach an approximate solution is estimated on the basis of the measurement error, according to the discrepancy principle. The effectiveness of the stopping rule is studied on simulated data and validated with experiments. The quality of the reconstruction depends on the information available on the source itself and the possibility to add this knowledge to the solution process is investigated. The method can produce good approximations provided that the amount of noise in the data can be estimated.