Richardson-Lucy deconvolution ?

Started by Fokko Beekhof February 25, 2008

I'm trying to understand the Richardson-Lucy deconvolution algorithm,
primarily from the description here:

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

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