Detail-preserving PET reconstruction with sparse image representation and anatomical priors

Abstract

Positron emission tomography (PET) reconstruction is an illposed inverse problem which typically involves fitting a high-dimensional forward model of the imaging process to noisy, and sometimes undersampled photon emission data. To improve the image quality, prior information derived from anatomical images of the same subject has been previously used in the penalised maximum likelihood (PML) method to regularise the model complexity and selectively smooth the image on a voxel basis in PET reconstruction. In this work, we propose a novel perspective of incorporating the prior information by exploring the sparse property of natural images. Instead of a regular voxel grid, the sparse image representation jointly determined by the prior image and the PET data is used in reconstruction to leverage between the image details and smoothness, and this prior is integrated into the PET forward model and has a closed-form expectation maximisation (EM) solution. Simulations show that the proposed approach achieves improved bias versus variance trade-off and higher contrast recovery than the current state-of-the-art methods, and preserves the image details better. Application to clinical PET data shows promising results.