Background
The combination of magnetic resonance imaging (MRI)-anatomical information with positron emission tomography (PET) image reconstruction has been shown to improve PET image quality in terms of spatial resolution and image noise, especially in brain PET imaging. There are different approaches to combine MRI and PET available, being the use of a Bayesian framework the most extended. Generally, the strength of the prior is controlled by a hyperparameter that needs to be tuned depending on the acquired statistics/counts and the desired image quality in the resulting PET image. However, comparisons between methods is scant, and it is not clear how sensitive they are to the different levels of statistics that can be measured in a PET scan.
Methods
In this work we employed maximum a posteriori (MAP) reconstruction with MRI information to guide the PET reconstruction, and evaluated the performance of several prior models and optimization methods with a fixed and adaptive hyperparameter. Different simulated scenarios, and measured data using different radiotracers at different levels of statistics were employed for the evaluation. Comparisons in image quality and quanti cation between methods were performed in different brain cortical and sub-cortical regions.
Results
Simulated data showed that an adaptive hyperparameter consistently outperformed a fixed hyperparameter for every image reconstruction algorithm implemented. The best performance was achieved with a model combining the Bowsher prior weighted by similarity coefficients based on joint entropy between PET and MRI. One-step-late (OSL) and preconditioned gradient ascent (PGA) optimization methods performed similarly at any level of statistics and number of iterations, so long as the hyperparameter was adaptive.
Conclusions
Results with simulated and measured data agreed that MAP reconstruction outperformed OSEM reconstruction, especially at low level of statistics, without any need of tuning. High resolution and low noise images were obtained using MAP reconstruction for 5{30 minutes scan times, showing negligible image quality difference for different radiotracers.