Solving non-negative linear inverse problems with the NeAREst method

This paper introduces the theoretical development of a numerical method, named NeAREst, for solving non-negative linear inverse problems, which arise often from physical or probabilistic models, especially, in image estimation with limited and indirect measurements. The Richardson-Lucy (RL) iteration is omnipresent in conventional methods that are based on probabilistic assumptions, arguments and techniques. Without resorting to probabilistic assumptions, NeAREst retains many appealing properties of the RL iteration by utilizing it as the substrate process and provides much needed mechanisms for acceleration as well as for selection of a target solution when many admissible ones exist.