Generalization of median root prior reconstruction

Penalized iterative algorithms for image reconstruction in emission tomography contain conditions on which kind of images are accepted as solutions. The penalty term has commonly been a function of pairwise pixel differences in the activity in a local neighborhood, such that smooth images are favored. Attempts to ensure better edge and detail preservation involve difficult tailoring of parameter values or the penalty function itself. The previously introduced median root prior (MRP) favors locally monotonic images. MRP preserves sharp edges while reducing locally nonmonotonic noise at the same time. Quantitative properties of MRP are good, because differences in the neighboring pixel values are not penalized as such. The median is used as an estimate for a penalty reference, against which the pixel value is compared when setting the penalty. In order to generalize the class of MRP-type of priors, the standard median was replaced by other order statistic operations, the L and finite-impluse-response median hybrid (FMH) filters. They allow for smoother appearance as they apply linear weighting together with robust nonlinear operations. The images reconstructed using the new MRP-L and MRP-FMH priors are visually more conventional. Good quantitative properties of MRP are not significantly altered by the new priors.