Image Deconvolution under Poisson Noise
using SURE-LET Approach
Feng Xue, Jiaqi Liu, Gang Meng, Jing Yan and Min Zhao
National Key Laboratory of Science and Technology on Test Physics and
Numerical Mathematics, Beijing, 100076, P. R. China
ABSTRACT
We propose an image deconvolution algorithm when the data is contaminated by Poisson noise. By minimizing
Stein’s unbiased risk estimate (SURE), the SURE-LET method was firstly proposed to deal with Gaussian noise
corruption. Our key contribution is to demonstrate that the SURE-LET algorithm is also applicable for Poisson
noisy image and proposed an efficient algorithm.
The formulation of SURE requires knowledge of Gaussian noise variance. We experimentally found a simple
and direct link between the noise variance estimated by median absolute difference (MAD) method and the
optimal one that leads to the best deconvolution performance in terms of mean squared error (MSE). Extensive
experiments show that this optimal noise variance works satisfactorily for a wide range of natural images.
Keywords: Image deconvolution, SURE-LET, Poisson noise, noise variance, MAD
1. INTRODUCTION
Deconvolution is a longstanding problem in many areas of signal and image processing, e.g. biomedical imag-
ing
1
and astronomy.
2, 3
In presence of Poisson noise, several deconvolution methods have been proposed such
as Tikhonov-Miller inverse filter and Richardson-Lucy (RL) algorithms,
4, 5
and also been comprehensively re-
viewed.
2, 6
The RL has been extensively used in many applications, but as it tends to amplify the noise after a
few iterations, several extension have been proposed. For example, wavelet-regularized RL algorithm has been
proposed by several authors.
2, 7
Recently, it is increasingly recognized that most of those algorithms proposed for Gaussian noise are ap-
plicable and effective for Poisson noise.
8, 9
Both problems with different noise distributions share very similar
deconvolution methods, e.g. sparsity-promoting regularization,
8
total variation (TV) regularization and the
related optimization algorithms, e.g. ADMM and augmented Lagrangian algorithms.
10
Note that the recently proposed SURE-LET algorithm achieves the most state-of-the-art deconvolution per-
formance under Gaussian noise assumption.
11
By minimizing Stein’s unbiased risk estimate (SURE),
12, 13
SURE-
LET method linearly parametrizes the deconvolution process and simplifies the deconvolution as solution to a
low-order linear system of equations.
11
This paper is to develop a new algorithm to deal with Poisson noise,
based on the SURE-LET framework.
11
The key contribution of this work is to propose a simple method to
optimize the SURE as objective functional.
The paper is organized as follows. Section 2 is devoted to review the SURE-LET approach. In Section 3, we
propose a novel algorithm for Poisson noise within SURE-LET framework. Section 4 reports the experimental
results for discussion. Some concluding remarks are finally given in Section 5.
Throughout this paper, we use boldface lowercase letters, e.g. x ∈ R
N
, to denote N-dimensional real vectors,
where N is typically the number of pixels in an image. The n-th element of x is written as x
n
. The linear
transformations (matrices) R
N
→ R
M
are denoted by boldface uppercase letters, e.g. H ∈ R
M×N
. H
T
∈ R
N×M
denotes the transpose of matrix H.
This work was supported by the National Natural Science Foundation of China under Grant No. 61401013.
Further author information: (Send correspondence to Feng Xue)
Feng Xue: E-mail: fxue2012@gmail.com, Telephone: +86 10-88520038
Jiaqi Liu: ljq006@vip.sina.com; Gang Meng: mgang2012@yeah.net;
Jing Yan: jyan1998@126.com; Min Zhao: mzhao2012@yeah.net
AOPC 2015: Image Processing and Analysis, edited by Chunhua Shen, Weiping Yang, Honghai Liu,
Proc. of SPIE Vol. 9675, 96750B · © 2015 SPIE · CCC code: 0277-786X/15/$18 · doi: 10.1117/12.2197291
Proc. of SPIE Vol. 9675 96750B-1
Downloaded From: http://proceedings.spiedigitallibrary.org/ on 10/11/2015 Terms of Use: http://spiedigitallibrary.org/ss/TermsOfUse.aspx