Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms.

The relative performance of five fully 3D PET reconstruction algorithms is evaluated. The algorithms are a filtered backprojection (FBP) method and two variants each of the EM-ML and ART iterative methods. For each of the iterative methods, one variant makes use of voxels and the other makes use of 'blobs' (spherically symmetric functions smoothly decaying to zero at their boundaries) as basis functions in its discrete reconstruction model. The methods are evaluated from the point of view of the efficacy of the reconstructions produced by them for three typical medical tasks--estimation of the average activity inside specific regions of interest, detection of hot spots, and detection of cold spots. A free parameter is allowed in the description of each of the five algorithms; the parameters are determined by a training process during which a value of the free parameter is selected which (nearly) maximizes a technical figure of merit. Such training and the actual comparative evaluation is done by making use of randomly generated phantoms and their projection data. The methodology allows assignation of levels of statistical significance to claims of the relative superiority of one algorithm over another for a particular task. We find that using blobs as basis functions in the iterative algorithms is definitely advantageous over using voxels. This result has high statistical significance. (We also include a visual illustration of it.) Comparing FBP, EM-ML using blobs, and ART using blobs, we do not find a clear difference in the overall performance of the investigated variants of the methods. If anything, our results suggest that ART using blobs may be the most efficacious of the three.