A method of iteratively reconstructing images from data acquired by converging collimation SPECT, comprising the steps of dividing a reconstruction space into multiple iso-resolution zones that are oriented non-parallel to the detector, and correcting for attenuation by applying a same attenuation value to all voxels in the same iso-resolution zone.