Single photon emission computed tomography (SPECT) is performed with multiple emission energies. For quantitative or qualitative SPECT, the image formation process for emissions at different energy ranges is modeled (44, 46, 48, 50) separately. Different scatter, different attenuation, and/or different collimator-detector response models corresponding to different energy ranges are used in the reconstruction.