Abstract

We consider a massively parallel implementation of Richardson-Lucy or maximum likelihood restoration with a spatially-variant point spread function (PSF). Richardson-Lucy iterates involve the computation of sums of the form: where O(x'_q) is the incident optical field estimate at discrete source location x'_q, I(x_q) is the measured discrete image at discrete field location x_q, and P(x_q, x'_q) is the discrete PSF – the probability that a photon from source region x'_q is incident on the detector at field region x_q. In general P is a function of source and field coordinates, and the computational burden of Eq. 1 is intractably large.

© 1995 Optical Society of America