Theoretical Evaluation of the Detectability of Random Lesions in Bayesian Emission Reconstruction

Detecting cancerous lesion is an important task in positron emission tomography (PET). Bayesian methods based on the maximum a posteriori principle (also called penalized maximum likelihood methods) have been developed to deal with the low signal to noise ratio in the emission data. Similar to the filter cut-off frequency in the filtered backprojection method, the prior parameters in Bayesian reconstruction control the resolution and noise trade-off and hence affect detectability of lesions in reconstructed images. Bayesian reconstructions are difficult to analyze because the resolution and noise properties are nonlinear and object-dependent. Most research has been based on Monte Carlo simulations, which are very time consuming. Building on the recent progress on the theoretical analysis of image properties of statistical reconstructions and the development of numerical observers, here we develop a theoretical approach for fast computation of lesion detectability in Bayesian reconstruction. The results can be used to choose the optimum hyperparameter for the maximum lesion detectability. New in this work is the use of theoretical expressions that explicitly model the statistical variation of the lesion and background without assuming that the object variation is (locally) stationary. The theoretical results are validated using Monte Carlo simulations. The comparisons show good agreement between the theoretical predications and the Monte Carlo results

Analysis of Observer Performance in Unknown-Location Tasks for Tomographic Image Reconstruction

Our goal is to optimize regularized image reconstruction for emission tomography with respect to lesion detectability in the reconstructed images. We consider model observers whose decision variable is the maximum value of a local test statistic within a search area. Previous approaches have used simulations to evaluate the performance of such observers. We propose an alternative approach, where approximations of tail probabilities for the maximum of correlated Gaussian random fields facilitate analytical evaluation of detection performance. We illustrate how these approximations, which are reasonably accurate at low probability of false alarm operating points, can be used to optimize regularization with respect to lesion detectability.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85914/1/Fessler33.pd

Analysis of Observer Performance in Detecting Signals with Location Uncertainty for Regularized Tomographic Image Reconstruction

Our goal is to optimize regularized image reconstruction methods for emission tomography with respect to the task of detecting small lesions in the reconstructed images. To reflect medical practice realistically, we consider the location of the lesion to be unknown. This location uncertainty significantly complicates the mathematical analysis of model observer performance. We consider model observers whose decisions are based on finding the maximum value of a local test statistic over all possible locations. Khurd and Gindi (SPIE 2004) and Qi and Huesman (SPIE 2004) described analytical approximations of the moments of the local test statistics and used Monte Carlo simulations to evaluate the localization performance of such "maximum observers". We propose here an alternative approach, where tail probability approximations developed by Adler (AAP 2000) facilitate analytical evaluation of the detection performance of these observers. We illustrate how these approximations can be used to evaluate the probability of detection (for low probability of false alarm operating points) for the maximum channelized hotelling observer. Using our analyses, one can rank and optimize image reconstruction methods without requiring time-consuming Monte Carlo simulations.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85960/1/Fessler205.pd

Fast Predictions of Variance Images for Fan-Beam Transmission Tomography With Quadratic Regularization

Accurate predictions of image variances can be useful for reconstruction algorithm analysis and for the design of regularization methods. Computing the predicted variance at every pixel using matrix-based approximations is impractical. Even most recently adopted methods that are based on local discrete Fourier approximations are impractical since they would require a forward and backprojection and two fast Fourier transform (FFT) calculations for every pixel, particularly for shift-variant systems like fan-beam tomography. This paper describes new "analytical" approaches to predicting the approximate variance maps of 2-D images that are reconstructed by penalized-likelihood estimation with quadratic regularization in fan-beam geometries. The simplest of the proposed analytical approaches requires computation equivalent to one backprojection and some summations, so it is computationally practical even for the data sizes in X-ray computed tomography (CT). Simulation results show that it gives accurate predictions of the variance maps. The parallel-beam geometry is a simple special case of the fan-beam analysis. The analysis is also applicable to 2-D positron emission tomography (PET).Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86007/1/Fessler37.pd