Another nonlinear method Lucy-Rechardson algorithm[1]}, [2]} has been
successfully used on X-ray ring coded imaging in experiment[3]}.
And this algorithm has been most widely used for restoring images
with noise from counting statistics and the data approximate Poisson statistics.
In this program we first try to use the Lucy-Rechardson method to get reconstructed image.
The process of the imaging can be described as a spatial intensity distribution
\(g(x,y)\)  convolved by a neutron source spatial intensity distribution \(f(x,y)\)  with
a point spread function (PSF) \(h(x,y)\) . Add the noise \(n(x,y)\)  the detector
image can be expressed as
\(g(x,y)=f(x,y)\otimes h(x,y)+n(x,y)\) 
