Optimization of the Algorithm for the Implementation of Point Spread Function in the 3D-OSEM Reconstruction Algorithm Based on the List-Mode Micro PET Data

Abstract

Positron emission tomography (PET) is a popular research topic. People are becoming more interested in PET images as they become more widely available. However, the partial volume effect (PVE) in PET images remains one of the most influential factors causing the resolution of PET images to degrade. It is possible to reduce this PVE and achieve better image quality by measuring and modeling the point spread function (PSF) and then accounting for it inside the reconstruction algorithm. In this work, we examined the response characteristics of the MetisTM PET/CT system by acquiring 22Na point source at different locations in the field of view (FOV) of the scanner and reconstructing with small pixel size for images to obtain their radial, tangential, and axial full-width half maximum (FWHM). An image-based model of the PSF model was then obtained by fitting asymmetric two-dimensional Gaussians on the 22Na images. This PSF model determined by FWHM in three directions was integrated into a three-dimensional ordered subsets expectation maximization (3D-OSEM) algorithm based on a list-mode format to form a new PSF-OSEM algorithm. We used both algorithms to reconstruct point source, Derenzo phantom, and mouse PET images and performed qualitative and quantitative analyses. In the point source study, the PSF-OSEM algorithm reduced the FWHM of the point source PET image in three directions to about 0.67 mm, and in the phantom study, the PET image reconstructed by the PSF-OSEM algorithm had better visual effects. At the same time, the quantitative analysis results of the Derenzo phantom were better than the original 3D-OSEM algorithm. In the mouse experiment, the results of qualitative and quantitative analyses showed that the imaging quality of PSF-OSEM algorithm was better than that of 3D-OSEM algorithm. Our results show that adding the PSF model to the 3D-OSEM algorithm in the MetisTM PET/CT system helps to improve the resolution of the image and satisfy the qualitative and quantitative analysis criteria.