Estimating Precision in Functional Images

Abstract

Functional imaging of biologic parameters like in vivo tissue metabolism is made possible by Positron Emission Tomography (PET). Many techniques have been suggested for extracting such images from dynamic time-course sequences of reconstructed PET scans. Quantitating the precision of these estimates is important for drawing inferences on the biologic parameters. Analytic variance formulas are not immediate owing to the nonlinear methods used in extraction. The usual resampling approach is infeasible because each image reconstruction in PET is a computationally demanding solution to a high-dimensional linear inverse problem. We suggest an alternative simulation approach that approximates the distribution of reconstructed PET scans and performs a parametric bootstrap in the imaging domain. Results on a simplified model chosen to match the characteristics of PET reconstruction are very encouraging. Mixture analysis is used to estimate functional images; however, the suggested approach is general enough to extend to other techniques or imaging methods.