A proximal iteration for deconvolving Poisson noisy images using sparse representations

Abstract

We propose an image deconvolution algorithm when the data is contaminated by Poisson noise. The image to restore is assumed to be sparsely represented in a dictionary of waveforms such as the wavelet or curvelet transforms. Our key contributions are: First, we handle the Poisson noise properly by using the Anscombe variance stabilizing transform leading to a {\it non-linear} degradation equation with additive Gaussian noise. Second, the deconvolution problem is formulated as the minimization of a convex functional with a data-fidelity term reflecting the noise properties, and a non-smooth sparsity-promoting penalties over the image representation coefficients (e.g. ${\ell }_1$-norm). Third, a fast iterative backward-forward splitting algorithm is proposed to solve the minimization problem. We derive existence and uniqueness conditions of the solution, and establish convergence of the iterative algorithm. Finally, a GCV-based model selection procedure is proposed to objectively select the regularization parameter. Experimental results are carried out to show the striking benefits gained from taking into account the Poisson statistics of the noise. These results also suggest that using sparse-domain regularization may be tractable in many deconvolution applications with Poisson noise such as astronomy and microscopy.

References in corpus (2)

• On the "degrees of freedom" of the lasso
• Poisson inverse problems

Cited by in corpus (28)

• Sparse Regularization via Convex Analysis
• Restoration of Poissonian Images Using Alternating Direction Optimization
• Generalized Forward-Backward Splitting
• This is SPIRAL-TAP: Sparse Poisson Intensity Reconstruction ALgorithms - Theory and Practice
• Compressed sensing performance bounds under Poisson noise
• Proximal Splitting Methods in Signal Processing
• D${}^3$PO - Denoising, Deconvolving, and Decomposing Photon Observations
• Imaging Atomic-Scale Chemistry from Fused Multi-Modal Electron Microscopy
• Plug-and-Play Quantum Adaptive Denoiser for Deconvolving Poisson Noisy Images
• Projected Nesterov's Proximal-Gradient Algorithm for Sparse Signal Reconstruction with a Convex Constraint
• Non-parametric PSF estimation from celestial transit solar images using blind deconvolution
• Deconvolution of confocal microscopy images using proximal iteration and sparse representations
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Epigraphical splitting for solving constrained convex formulations of inverse problems with proximal tools
• Sparse Signal Recovery under Poisson Statistics
• Proximity Operators of Discrete Information Divergences
• Parallel Proximal Algorithm for Image Restoration Using Hybrid Regularization -- Extended Version
• Hierarchical Convex Optimization by the Hybrid Steepest Descent Method with Proximal Splitting Operators -- Enhancements of SVM and Lasso
• Deconvolution of Poissonian Images Using Variable Splitting and Augmented Lagrangian Optimization
• Deconvolution under Poisson noise using exact data fidelity and synthesis or analysis sparsity priors
• Reconstruction Error Bounds for Compressed Sensing under Poisson or Poisson-Gaussian Noise Using Variance Stabilization Transforms
• Performance Bounds for LASSO under Multiplicative Noise: Applications to Pooled RT-PCR Testing
• Class-specific Poisson denoising by patch-based importance sampling
• Reconstruction Error Bounds for Compressed Sensing under Poisson Noise using the Square Root of the Jensen-Shannon Divergence
• Nested iterative algorithms for convex constrained image recovery problems
• Properties on n-dimensional convolution for image deconvolution
• Multiview Attenuation Estimation and Correction
• Linear inverse problems with noise: primal and primal-dual splitting