LAUSR

Standard positron emission tomography (PET) reconstruction techniques are based on maximum-likelihood (ML) optimization methods, such as the maximum-likelihood expectation-maximization (MLEM) algorithm and its variations. Most methodologies rely on a positivity constraint on the activity distribution image. Although this constraint is meaningful from a physical point of view, it can be a source of bias for low-count/high-background PET, which can compromise accurate quantification. Existing methods that allow for negative values in the estimated image usually utilize a modified log-likelihood, and therefore break the data statistics. In this paper, we propose to incorporate the positivity constraint on the projections only, by approximating the (penalized) log-likelihood function by an adequate sequence of objective functions that are easily maximized without constraint. This sequence is constructed such that there is hypo-convergence (a type of convergence that allows the convergence of the maximizers under some conditions) to the original log-likelihood, hence allowing us to achieve maximization with positivity constraint on the projections using simple settings. A complete proof of convergence under weak assumptions is given. We provide results of experiments on simulated data where we compare our methodology with the alternative direction method of multipliers (ADMM) method, showing that our algorithm converges to a maximizer, which stays in the desired feasibility set, with faster convergence than ADMM. We also show that this approach reduces the bias, as compared with MLEM images, in necrotic tumors—which are characterized by cold regions surrounded by hot structures—while reconstructing similar activity values in hot regions.