Sparse Poisson Noisy Image Deblurring

• Published in: IEEE transactions on image processing 2012-04, Vol.21 (4), p.1834-1846
• Main Authors: Carlavan, M., Blanc-Feraud, L.
• Format: Article
• Language: English
• Subjects: 3-D confocal microscopy deconvolution; Algorithms; Alternating-direction method (ADM); Computer Science; Computer Simulation; constrained minimization; Deconvolution; discrepancy principle; Image Enhancement - methods; Image Interpretation, Computer-Assisted - methods; Image Processing; Microscopy; Microscopy, Confocal - methods; Models, Statistical; Noise; Poisson Distribution; Poisson noise; regularizing parameter selection; Reproducibility of Results; Sensitivity and Specificity; Signal-To-Noise Ratio; Three dimensional displays; Wavelet transforms
• ISSN: 1057-7149; 1941-0042
• DOI: 10.1109/TIP.2011.2175934

Deblurring noisy Poisson images has recently been a subject of an increasing amount of works in many areas such as astronomy and biological imaging. In this paper, we focus on confocal microscopy, which is a very popular technique for 3-D imaging of biological living specimens that gives images with a very good resolution (several hundreds of nanometers), although degraded by both blur and Poisson noise. Deconvolution methods have been proposed to reduce these degradations, and in this paper, we focus on techniques that promote the introduction of an explicit prior on the solution. One difficulty of these techniques is to set the value of the parameter, which weights the tradeoff between the data term and the regularizing term. Only few works have been devoted to the research of an automatic selection of this regularizing parameter when considering Poisson noise; therefore, it is often set manually such that it gives the best visual results. We present here two recent methods to estimate this regularizing parameter, and we first propose an improvement of these estimators, which takes advantage of confocal images. Following these estimators, we secondly propose to express the problem of the deconvolution of Poisson noisy images as the minimization of a new constrained problem. The proposed constrained formulation is well suited to this application domain since it is directly expressed using the antilog likelihood of the Poisson distribution and therefore does not require any approximation. We show how to solve the unconstrained and constrained problems using the recent alternating-direction technique, and we present results on synthetic and real data using well-known priors, such as total variation and wavelet transforms. Among these wavelet transforms, we specially focus on the dual-tree complex wavelet transform and on the dictionary composed of curvelets and an undecimated wavelet transform.