Synthesized Image Reconstruction for Post-Reconstruction Resolution Recovery

Resolution recovery (RR) techniques in positron emission tomography (PET) imaging aim to mitigate spatial resolution losses and related inaccuracies in quantification by using a model of the system's point spread function (PSF) during reconstruction or post-processing. However, including PSF modeling in fully 3-D image reconstruction is far from trivial as access to the scanner-specific forward and back-projectors is required, along with access to the 3-D sinogram data. Hence, post-reconstruction RR methods, such as the Richardson-Lucy (RL) algorithm, can be more practical. However, the RL method leads to relatively rapid noise amplification in early image iterations, giving inferior image quality compared to iterates obtained by placing the PSF model in the reconstruction algorithm. We propose a post-reconstruction RR method by synthesizing PET data by a forward projection of an initial real data reconstruction (such reconstructions are usually available via a scanner's standard reconstruction software). The synthetic PET data are then used to reconstruct an image, but crucially now including a modeled PSF within the system model used during reconstruction. Results from simulations and real data demonstrate the proposed method improves image quality compared to the RL algorithm, whilst avoiding the need for scanner-specific projectors and raw sinogram data (as required by standard PSF modeling within reconstruction).