Noise amplification and ill-convergence of Richardson-Lucy deconvolution

Richardson-Lucy (RL) deconvolution optimizes the likelihood of the object estimate for an incoherent imaging system. It can offer an increase in contrast, but converges poorly, and shows enhancement of noise as the iteration progresses. We have discovered the underlying reason for this problematic convergence behaviour using a Cramér Rao Lower Bound (CRLB) analysis. An analytical expression for the CRLB diverges for spatial frequency components that approach the diffraction limit from below. The resulting mean noise variance per pixel diverges for large images. These results imply that a regular optimum of the likelihood does not exist, and that RL deconvolution is necessarily ill-convergent.