This paper introduces the theoretical development of a numerical method, named NeAREst, for solving nonnegative 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 an...