Noise amplification and ill-convergence of Richardson-Lucy deconvolution
Y. Liu (TU Delft - ImPhys/Stallinga group)
Spozmai Panezai (TU Delft - ImPhys/Rieger group)
Y. Wang (TU Delft - ImPhys/Stallinga group)
S Stallinga (TU Delft - ImPhys/Stallinga group)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
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.