Reducing negativity artifacts in emission tomography: Post-processing filtered backprojection solutions

The problem of negative artifacts in emission tomography reconstructions computed by filtered backprojection (FBP) is of practical concern particularly in low count studies. Statistical reconstruction methods based on maximum likelihood (ML) are automatically constrained to be non-negative but their excessive computational overhead has limited their use in operational settings. Motivated by the statistical character of the negativity artifact, the authors develop a simple post-processing technique that iteratively adjusts negative values by cancellation with positive values in a surrounding local neighborhood. The compute time of this approach is roughly equivalent to two applications of FBP. The approach was evaluated by numerical simulation in one- and two-dimensional settings. The studies compared smoothed versions of FBP, the post-processed FBP, and ML implemented by the expectation-maximization algorithm. The root mean square (RMS) error between the true and estimated source distribution was used to evaluated performance; in two dimensions, additional region-of-interest-based measures of reconstruction accuracy were also employed. In making comparisons between the different methods, the amount of smoothing applied to each reconstruction method was adapted to minimize the RMS error--this was found to be critical. After adjusting for this effect, the average RMS error for FBP was typically between 13% and 20% higher than ML. Similar results were found for the region-of-interest error