Reply by Fokko Beekhof●February 25, 20082008-02-25

Hello,
I'm trying to understand the Richardson-Lucy deconvolution algorithm,
primarily from the description here:
http://en.wikipedia.org/wiki/Richardson-Lucy_deconvolution
In this notation, u_i is a pixel in the latent image(unknown), c_i the
corresponding pixel in the observed image, and p_{i, j} is the
point-spread-function, i.e. the fraction of the pixel at j observed at i.
This is an iterative scheme, so an estimate for u is updated in every
iteration.
There are some things not entirely clear to me:
1) What should be the initial value for the latent image U ?
2) How does one detect convergence ?
Any help would be much appreciated!
F. Beekhof