Paper
21 May 1999 Fast computation of the covariance of MAP reconstructions of PET images
Jinyi Qi, Richard M. Leahy
Author Affiliations +
Abstract
We develop an approximate theoretical formula for fast computation of the covariance of PET images reconstructed using maximum a posteriori (MAP) estimation. The results assume a Poisson likelihood for the data and a quadratic prior on the image. The covariance for each voxel is computed using 2D FFTs and is a function of a single data dependent parameter. This parameter is computed using a modified backprojection. For a small region of interest (ROI), the correlation can be assumed to be locally stationary so that computation of the variance of an ROI can be performed very rapidly. Previous approximate formulae for the variance of MAP estimators have performed poorly in areas of low activity since they do not account for the non- negativity constraints that are routinely used in MAP algorithms. Here a `truncated Gaussian' model is used to compensate for the effect of the non-negativity constraints. Accuracy of the theoretical expressions is evaluated using both Monte Carlo simulations and a multiple-frame 15O- water brain study. The Monte Carlo studies show that the truncated Gaussian model is effective in compensating for the effect of the non-negativity constraint. These results also show good agreement between Monte Carlo covariances and the theoretical approximations. The 15O-water brain study further confirms the accuracy of the theoretical approximations.
© (1999) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Jinyi Qi and Richard M. Leahy "Fast computation of the covariance of MAP reconstructions of PET images", Proc. SPIE 3661, Medical Imaging 1999: Image Processing, (21 May 1999); https://doi.org/10.1117/12.348589
Lens.org Logo
CITATIONS
Cited by 13 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Monte Carlo methods

Positron emission tomography

Brain

Brain mapping

Fourier transforms

Tin

Data modeling

Back to Top