Improved modeling of system response in list mode EM reconstruction of Compton scatter camera images
An improved List Mode EM method for reconstructing Compton scattering camera images has been developed. First, an approximate method for computation of the spatial variation in the detector sensitivity has been derived and validated by Monte Carlo computation. A technique for estimating the relative weight of system matrix coefficients for each gamma in the list has also been employed, as has a method for determining the relative probabilities of emission having some from pixels tallied in each list-mode back-projection. Finally, a technique has been developed for modeling the effects of Doppler broadening and finite detector energy resolution on the relative weights for pixels neighbor to those intersected by the back-projection, based on values for the FWHM of the spread in the cone angle computed by Monte Carlo. Memory issues typically associated with list mode reconstruction are circumvented by storing only a list of the pixels intersected by the back-projections, and computing the weights of the neighboring pixels at each iteration step. Simulated projection data has been generated for a representative Compton camera system (CSPRINT) for several source distributions and reconstructions performed. Reconstructions have also been performed for experimental data for distributed sources.