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