Research Article  Open Access
Restoring Poissonian Images by a Combined FirstOrder and SecondOrder Variation Approach
Abstract
The restoration of blurred images corrupted by Poisson noise is an important topic in imaging science. The problem has recently received considerable attention in recent years. In this paper, we propose a combined firstorder and secondorder variation model to restore blurred images corrupted by Poisson noise. Our model can substantially reduce the staircase effect, while preserving edges in the restored images, since it combines advantages of the firstorder and secondorder total variation. We study the issues of existence and uniqueness of a minimizer for this variational model. Moreover, we employ a gradient descent method to solve the associated EulerLagrange equation. Numerical results demonstrate the validity and efficiency of the proposed method for Poisson noise removal problem.
1. Introduction
Image restoration problem has been widely studied in the areas of image processing. The goal of image restoration is to reconstruct an approximation of an original image from a blurred and noisy one [1–4]. Over the last decades, most of the literature deal with the additive noise and multiplicative noise model [5–11]. Given a blurred and noisy image or , where is the original image, represents the noise, and the blurring operator is a point spread function (PSF); the restoration problem involving the additive noise or multiplicative noise is to recover from the observed image . The additive noise and multiplicative noise models have been extensively studied. We refer the reader to [12–19] and references therein for a review of various methods. In fact, there are other types of noise such as Poisson noise. In this paper, we consider the problem of seeking the approximations of original images from blurred images corrupted by Poisson noise. The restoration of blurred images corrupted by Poisson noise is an important task in various applications such as astronomical imaging, electronic microscopy, singleparticle emissioncomputed tomography (SPECT), and positron emission tomography (PET) [20–24]. In a Poisson noise model, the recorded image is a realization of a random variable . A statistical model for is given by where is the expected value of the background which is assumed to be constant. The given data of the problem are the imaging operator , the expected value of the background , and the detected image . It is known that the major difficulty for finding the original image of (1) is the models and algorithms developed for deconvolution of additive noise, or multiplicative noise images cannot be directly applied to the restoration of Poissonian images. Therefore, it is important to devise efficient and reliable algorithms for restoring blurred images corrupted by Poisson noise. The problem of restoration of Poissonian images has received considerable attention in recent years.
We note that, in a discrete setting, the intensity in the pixel is a random variable that follows a Poisson law of mean where the structured matrix related to the boundary conditions is called a blurring matrix. Considering that these random variables are independent, the probability distribution of , for given , , and , can be written as The use of Stirlings formula leads to the wellknown Csiszár [25] divergence measure between and . Indeed, by taking the negative logarithm of the likelihood function, we obtain This measure generalizes the KullbackLeibler divergence or crossentropy measure to accommodate functions whose integrals are not constant, as they would be if they were probability distributions. Dropping the terms independent of in , we obtain Therefore, the maximum likelihood estimator of the original image is the minimizer with respect to of the negativelog Poisson likelihood functional .
As is well known, the structured matrix has many singular values of different orders of magnitude close to the origin. This makes the minimizer of (4) very sensitive to the noise in the righthand side . Thus a regularization method may be employed to compute the approximate solution that is less sensitive to noise than the naive solution. In general, regularization methods formulate the image restoration problem as a minimization problem of the form: where is a regularization function and is the regularization parameter that controls the balance between and . We note that the regularization parameter is an important quantity which controls the properties of the regularized solution, and should therefore be chosen with care. Throughout the years a variety of parameter choice strategies such as the discrepancy principle, the curve, and generalized crossvalidation (GCV) have been proposed for Poisson noise removal problem [26–28].
It is shown in [29] that the choice of the regularization function is related to the features of the images to be restored. The use of different regularization functions is an active research area. Probably one of the most popular regularization functions is the Tikhonov regularization function , where is an auxiliary operator chosen among the identity and loworder differential operators; see [30] for more details. In [31], Landi and Piccolomini proposed a quasiNewton projection method for deblurring Poissoncorrupted images after formulating the image restoration problem as a nonnegatively constrained minimization problem where the objective function is the sum of the KullbackLeibler divergence, used to express fidelity to the data in the presence of Poisson noise, and of a Tikhonov regularization term. Their numerical results show the potential of the method both in terms of relative error reduction and computational efficiency. In [32, 33], the authors developed a cost functional which incorporates the statistics of the noise in the image data and Tikhonov regularization to induce stability and employed an efficient hybrid gradient projectionreduced Newton (GPRNCG) method for the problem of restoration of Poissonian images. Finally, they gave a comparison between the RichardsonLucy [34, 35] iterative regularization method and GPRNCG algorithm for the restoration of blurred images corrupted by Poisson noise. By the measures of performance used in this comparison, GPRNCG was more efficient than RichardsonLucy iteration.
It is known that Tikhonov regularization estimate is similar to lowpass filtering; therefore, it produces a smoothing effect on the restored images; that is, it penalizes edges, which is not a good approximation of the original image if it contains edges. To overcome this shortcoming, Rudin et al. [5] proposed a total variation (TV) based regularization technique, which preserves the edge information in the restored image. In the case of TV regularization, the estimated solution is obtained by minimizing the object function: where is the discrete TV regularization term. It is shown in [36] that the TV penalty function acts more like a model selection device by identifying a small number of critical points at which is allowed to jump. The number of these selected jump points is controlled by . Especially, if is piecewise constant with a finite number of jump discontinuities, then gives the sum of magnitudes of the jumps.
In [21], Bardsley and Luttman considered the TV regularization for poissonian images and showed that the problem of computing the minimizer of the resulting TVpenalized Poisson likelihood functional is well posed. They then proved that as the errors in the data and in the blurring operator tend to zero, the resulting minimizer converges to the minimizer of the exact likelihood function. Finally, the practical effectiveness of the approach is demonstrated on synthetically generated data, and a nonnegatively constrained, projected quasiNewton method was introduced. In [22], Figueiredo and BioucasDias proposed an approach to deconvolving Poissonian images with the TV function, which is based upon an alternating direction optimization method. Using standard convex analysis tools, they presented sufficient conditions for existence and uniqueness of solutions of these optimization problems. In [23], Setzer et al. considered the restoration of blurred images corrupted by Poisson noise by minimizing an energy functional consisting of the KullbackLeibler divergence as similarity term and the TV regularization term. Their minimizing algorithm uses alternating split Bregman techniques which can be reinterpreted as DouglasRachford splitting applied to the dual problem. In contrast to recently developed iterative algorithms, their algorithm contains no inner iterations and produces nonnegative images. The high efficiency of their algorithm in comparison to other recently developed algorithms is demonstrated by artificial and realworld numerical examples. We refer the reader to [24, 37] and references therein for other efficient methods solving the restoration of blurred images corrupted by Poisson noise with the TV function.
Although the total variation regularization is extremely popular in a variety of applications, it also gives rise to some undesired effects. Both from a theoretical and experimental point of view, it has been shown that the TV norm transforms smooth signal into piecewise constants, the socalled staircase effect. Staircase solutions fail to satisfy the evaluation of visual quality, and they can develop false edges that do not exist in the true image. To attenuate the staircase effect, there is a growing interest in the literature for replacing the TV norm by a highorder TV norm. The motivation behind such attempt is to restore potentially a wider range of images, which comprise more than merely piecewise constant regions. The majority of the highorder norms involve secondorder differential operators because piecewisevanishing secondorder derivatives lead to piecewiselinear solutions that better fit smooth intensity changes; see [38] for more details.
Secondorder regularization schemes have been considered so far in the literature mainly for dealing with the staircase effect while preserving the edge information in the restored image. There are two classes of secondorder regularization methods for image restoration problems. The first class combines a secondorder regularizer with the TV norm. For example, a technique in [39, 40] combining the TV filter with a fourthorder PDE filter was proposed to preserve edges and at the same time avoid the staircase effect in smooth regions for noise removal. In [41], Papafitsoros and Schönlieb considered a highorder model involving convex functions of the TV and the TV of the first derivatives for image restoration problems and used the split Bregman method to solve numerically the corresponding discretized problem. Chan et al. proposed a secondorder model to substantially reduce the staircase effect, while preserving sharp jump discontinuities for image restoration in [42]. The second class employs a secondorder regularizer in a standalone way. The secondorder regularizer was first proposed in [43] for additive noise removal in medical magnetic resonance images. The LLT model presented in [43] can overcome the staircase effect that occurs with the TV norm filter and better preserve the fine details in nonblocky images. In order to test its practical potential the authors have applied their method to a wide range of real images, including structural and functional MRI data sets. The main strength of their method was the ability to process signals with a smooth change in the intensity value. Numerical results have shown that, compared with some related partial differential equation (PDE) models, the LLT model is rather robust in removing noise and handling edges. In [44, 45], the authors reconsidered the fourthorder PDE model for additive noise removal and employed the dual algorithm of Chambolle for solving the highorder problems.
In order to get restored images with edge preserved and simultaneously smooth region reconstructed for Poisson noise removal, it is natural to utilize a combined firstorder and secondorder total variation technique. In this paper, we consider to modify the total variation model by adding a highorder functional for restoring blurred images corrupted by Poisson noise. Our model can substantially reduce the staircase effect, while preserving edges in the restored images, since it combines advantages of the firstorder and secondorder total variation. We study the issues of existence and uniqueness of a minimizer for this variational model. Moreover, we employ a gradient descent method to solve the associated EulerLagrange equation. Some numerical results demonstrate the validity and efficiency of the proposed method for Poisson noise removal problem.
The organization of this paper is outlined as follows. In the next section, we present a combined firstorder and secondorder variation model to restore blurred images corrupted by Poisson noise and study the issues of existence and uniqueness of a minimizer for this variational model. A gradient descent method to solve the associated EulerLagrange equation is presented in Section 3. Some numerical experiments are given to illustrate the performance of the proposed algorithm in Section 4. We give concluding remarks in the last section.
2. A Combined FirstOrder and SecondOrder Variation Model
In the Bayesian approach, a prior probability density for is also specified, and the posterior density: given by Bayes’s law, is maximized with respect to . In (7), defines the prior probability distribution of the measurements , which are given by measurements, and the likelihood function is the conditional probability of observing a fixed for a variable . The maximizer of is called the maximum a posteriori (MAP) estimate. Since does not depend on , it is not needed in the computation of the MAP estimate and thus can be ignored. Taking the natural logarithm of (7) and dropping the terms independent of , we obtain that maximizing (7) is equivalent to minimizing with respect to , where . It is clear that corresponds to the regularization term in the classical penalized likelihood approach to regularization. However in the Bayesian setting, is the probability density known as the prior from which the unknown is assumed to arise. Thus the prior knowledge regarding the characteristics of can be formulated in the form of a probability density , and this yields to a natural, and statistically rigorous, motivation for the regularization method [27].
Standard Tikhonov regularization corresponds to the following choice of the prior: This corresponds to the assumption that the prior for is a zeromean Gaussian random variable with covariance matrix , which has the effect of penalizing reconstructions with large norm. For norm of the gradient regularization, the penalty has the similar form where is the gradient operator. The use of this regularization function has the effect of penalizing reconstructions that are not smooth. As an alternative to Tikhonov regularization including standard Tikhonov regularization and gradient regularization for Poisson noise removal problem, total variation regularization is another regularization technique that allows for the presence of sharp edges in the resulting reconstruction. In the case of total variation regularization, we have where .
Optimization techniques such as alternation direction method of multipliers (ADMM) and alternating split Bregman algorithm were proposed for the restoration of blurred images corrupted by Poisson noise by minimizing an energy functional consisting of the KullbackLeibler divergence as the similarity term and the TV regularization term; see [21–24] for more details. In the other hand, some Newtonrelated gradient methods were proposed for deblurring Poissoncorrupted images by formulating the image restoration problem as a nonnegatively constrained minimization problem where the objective function is the sum of the KullbackLeibler divergence, used to express fidelity to the data in the presence of Poisson noise, and of a Tikhonov regularization term; see [31–33] for more details.
To attenuate the staircase effect in the TV regularization, we consider a combined firstorder and secondorder variation model to restore blurred images corrupted by Poisson noise in this work. In the proposed model, the prior probability density is of the form: where . It leads to present the following functional for restoring blurred images corrupted with Poisson noise: where the parameter is used to control the balance between the edges and the smooth surface and is the regularization parameter which measures the tradeoff between the fidelity term and the regularized term. The parameter can be computed by using the information of the edges and smooth regions of the resulting image obtained by smoothing the observed image with the lowpass filters such as the median filter and the Gauss filter.
The functional in (13) is defined on the set of such that ; in particular, must be positive almost everywhere. Some basic notations and properties concerning the space and space can be found in [39, 40]. Motivated by [46], we have the following result to show the existence and uniqueness of the minimizer for the model (13).
Theorem 1. Let be a bounded, open subset of with Lipschitz boundary. Assume that is a positive bounded function and is positive definite. Then for such that has a unique minimizer for the model (13).
Proof. First, a straightforward computation can show that is bounded below by . Since the penalty function , we obtain that is bounded below. Thus, we can choose a sequence such that . Hence, and are bounded. Using Jensen’s inequality, we have
which indicates that is bounded. This and the boundedness of yields that the sequence is bounded in . Following the compactness of in , we deduce that there exists such that a subsequence converges to pointwise almost everywhere in . By the lower semicontinuity of the space, we obtain that . Since is bounded below, we get that by using Fatou’s lemma. Therefore, and minimizes for all such that .
Since is positive definite and is positive, it immediately follows that is strictly convex. Obviously, is a convex function. We conclude that is strictly convex, as the sum of a convex function and of a strictly convex function. Therefore, the minimizer is unique.
3. Computational Method
In this section, our aim is to propose a timemarching gradient descent method for computing the minimizer of the model (13). From [5, 39, 40] we know that minimizing (13) for a given constant yields the associate EulerLagrange equation: with the boundary conditions where denotes the unit outernormal vector of .
Using the gradient descent method, we are able to derive the associated heat flow for the model (13): In this work, we use the finite difference scheme to discretize (17); see [5, 39, 40] for more details. Denoting the step space by , we employ the following discretization used in the implementations; see Table 1.
In the discretization, the notation and is near 0. Denoting the time step by , we get the following explicit computational scheme:
Now we discuss the choice of the weighting parameter in (18). Due to the strengths and weaknesses of the firstorder and secondorder variation approach, it is desirable that the weighting parameter along edges and in flat regions, emphasizing the restoration properties for the firstorder total variation. To emphasize the restoration properties of the secondorder total variation in smooth regions, we want . Specially, the resulting algorithm of (18) is just the TV regularization method for Poisson noise removal problem when , while the resulting algorithm is the highorder TV regularization method when . Usually, we may compute the parameter by using the information of the edges and smooth regions of the resulting image obtained by smoothing the observed image with the lowpass filters such as the median filter and the Gauss filter. We have carried out some numerical experiments. We observe that the fixed obtained by computing the observed image can give very good results. In this work, in order to detect edges and smoothing regions much better, we adopt the method for updating as proposed in [40]. The results obtained by carrying out various numerical examples show that the updating procedure behaves better for our model than the fixed . It is because, as the iteration proceeds, the edges and smoothing regions of recovered image are closer to the original image, then the parameter computed by the updating scheme can be better suitable for restoration. So we employ the method in [40] for updating in our numerical experiments. More precisely, suppose that is the th iterative solution, and we update the parameter as follows: where , and it means that only the absolute largest jumps are unaffected of the highorder regularization. As reported in [40], for large and small values of the parameter is closer to , and for intermediate values of the parameter approaches 0, which means that the highorder filter dominates the computation and the staircase effect is suppressed. Since only small jumps should be suppressed with the highorder regularization, it is a good choice for .
We are now in a position to describe the timemarching gradient descent algorithm (Algorithm 1) for restoring blurred images corrupted by Poisson noise.
4. Numerical Experiments
In this section, we present some numerical results to illustrate the performance of the proposed model for Poisson noise removal problem. We compare our model with the one proposed in [21] (TV method), the one proposed in [31] (Tikhonov method), and the one proposed in [47] (HTV method). All computations of the present paper were carried out in Matlab 7.10. The results were obtained by running the Matlab codes on an Intel(R) Core(TM) i3 CPU (2.27 GHz, 2.27 GHz) computer with RAM of 2048 M.
The quality of the restoration results with different methods is compared quantitatively by using the signaltonoise ratio (SNR), the relative error (ReErr), and structural SIMilarity index (SSIM). They are defined as follows: where and are the ideal image and the restored image, respectively, and where and are averages of and , respectively, and are the variance of and , respectively, and is the covariance of and . The positive constants and can be thought of as stabilizing constants for nearzero denominator values. The SSIM is a wellknown quality metric used to measure the similarity between two images. The SSIM method was developed by Wang et al. [48] and is based on three specific statistical measures that are much closer to how the human eye perceives differences between images. In the following experiments, we will use SSIM map to reveal areas of high or low similarity between two images, the whiter SSIM map, and the closer between the two images.
In order to give a good comparison, we employ the timemarching gradient descent algorithm for the four models. In the timemarching gradient descent algorithm, the parameter is introduced to avoid divisions by zero. The initial guess is chosen to be the degraded image in all tests. We set the step size as 0.1 in order to obtain a stable iterative procedure. We terminate the iterations for the methods and accept as the computed approximation of the ideal image as soon as the maximum number of allowed outer iterations has been carried out or the relative differences between consecutive iterates satisfy In addition, we carry out many experiments with different regularization parameters in the models, andthe one with the best results is presented in this work. In this way, we have a fair comparison since we compare the restorations for four different methods.
In the first experiment, we consider the “Cameraman” image of size . The degraded image in Figure 1(b) is obtained by performing the blurring operation psfGauss proposed in [49] on the original image in Figure 1(a) with the background and adding the Poisson noise. Note that there is no parameter associated with the Poisson noise, but the noise magnitude depends on the absolute image intensities. The amount of noise in a region of the image increases with the intensity of the image there. The restored images by all four methods are shown in Figures 1(c)–1(f). From these figures, compared with the tikhonov regularization method, the TV regularization method, and the high TV regularization method, the proposed approach yields better results in image restoration since it avoids the staircase effect of the general TV methods, and the edges are preserved as well as the general TV methods. In Figures 2(c)–2(f), we have enlarged some details of the four restored images. As it is seen in the zoomed parts, the proposed method outperforms the other three methods. In Figures 2(c)–2(f), we also present the SSIM maps of the restored images in Figures 2(g)–2(j). Note that the SSIM maps of the restored image by the proposed method are whiter than the other methods; that is, our method can get better results. In Table 2, we compare their restoration results in SNRs and ReErrs. We observe from Table 2 that both the SNR and ReErr values of the restored images by the proposed method are better than those by the other three methods.

(a) True image
(b) Degraded image
(c) Tikhonov
(d) TV
(e) HTV
(f) Hybrid TV
(a) True
(b) Degraded
(c) Tikhonov
(d) TV
(e) HTV
(f) Hybrid TV
(g) Tikhonov
(h) TV
(i) HTV
(j) Hybrid TV
Moffat blur is considered in the second example. It is known that the PSF of an astronomical telescope is often modeled by the Moffat function. The “Lena” image shown in Figure 3(a) is degraded by the Moffat blur psfMoffat proposed in [50] with and the Poisson noise. The information of restored images by the four methods is displayed in Figures 3(c)–3(f). From the visual quality of restored images, the proposed regularization method is better than the other three methods. In Figures 4(c)–4(f), we display the zoom parts of the restored images (the part shown as the white rectangle in Figure 3(a)). From Figures 4(c)–4(f), we see that the restored image obtained by our method has more details than those by the other methods. The comparison of SSIM maps shown in Figures 4(g)–4(j) also proves that our method can get better results. We report the SNR and RelErr values in Table 2. From the table, we know that our method behaves much better.
(a) True image
(b) Degraded image
(c) Tikhonov
(d) TV
(e) HTV
(f) Hybrid TV
(a) True
(b) Degraded
(c) Tikhonov
(d) TV
(e) HTV
(f) Hybrid TV
(g) Tikhonov
(h) TV
(i) HTV
(j) Hybrid TV
In the third example, the original astronomical object shown in Figure 5(a) is degraded by the given PSF in Figure 5(b) and the Poisson noise. The degraded image is displayed in Figure 5(c). The relative error between the noisy image and the original image is 0.3150.
(a) True image
(b) PSF
(c) Degraded image
(d) Theta
(e) Tikhonov
(f) TV
(g) HTV
(h) Hybrid TV
In Figures 5(e)–5(h), we show the restored images by the four different methods. From these figures, we observe that the restored image by our model is more better than the other models in terms of the staircase effect and edge preservation. In addition, we plot the choice of in Figure 5(d). From the figure, we see that only small jumps are suppressed with the secondorder regularization. The SNRs and relative errors are reported in Table 2. The SNR by our method is higher than those by the other methods, and the relative error by our method is smaller than those by the other methods. From these figures and Table 2, we observe that the proposed method outperforms the tikhonov regularization method, the TV regularization method, and the high TV regularization method.
5. Conclusions
In this paper, we propose a new variational model to restore blurred images corrupted by Poisson noise. Based on the good feature of highorder functional, we propose a model by adding an extra highorder functional term in the total variation model. Our model combines advantages of the firstorder and secondorder total variation. It can substantially reduce the staircase effect, while preserving edges in the restored images. The issues of existence and uniqueness of a minimizer for this variational model is discussed. At last, we employ a gradient descent method to solve the associated EulerLagrange equation. The numerical experiments show that the proposed method outperforms some existing restoration methods in terms of the SNR, relative error, and SSIM map for Poisson noise removal problem. The comparisons between the reconstructed images obtained by the four methods show that the proposed one can alleviate the staircase effect significantly while preserve edges.
Conflict of Interests
All of the coauthors do not have a direct financial relation with the trademarks mentioned in our paper that might lead to a conflict of interests for any of the coauthors.
Acknowledgment
L. Jiang and J. Huang are supported by NSFC (10871034), X.G. Lv and J. Liu are supported by NSFC (60973015 and 61170311), Sichuan Province Sci. & Tech. Research Project (2011JY0002 and 12ZC1802), and Chinese Universities Specialized Research Fund for the Doctoral Program (20110185110020).
References
 H. Andrew and B. Hunt, Digital Image Restoration,, Prentice Hall, Englewood Cliffs, NJ, USA, 1977.
 M. Bertero and P. Boccacci, Introduction to Inverse Problems in Imaging, Institute of Physics Publishing, London, UK, 1998. View at: Publisher Site  Zentralblatt MATH  MathSciNet
 T. F. Chan and J. Shen, Image Processing and Analysis: Variational, PDE, Wavelet, and Stochastic Methods, SIAM, Philadelphia, Pa, USA, 2005. View at: Publisher Site  Zentralblatt MATH  MathSciNet
 M. Banham and A. Katsaggelos, “Digital image restoration,” IEEE Signal Processing Magazine, vol. 14, no. 2, pp. 24–41, 1997. View at: Google Scholar
 L. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Physical D, vol. 60, no. 1–4, pp. 259–268, 1992. View at: Publisher Site  Google Scholar
 A. Chambolle and P. L. Lions, “Image recovery via total variation minimization and related problems,” Numerische Mathematik, vol. 76, no. 2, pp. 167–188, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Rudin, P. L. Lions, and S. Osher, “Multiplicative denoising and deblurring: theory and algorithms,” in Geometric Level Set Methods in Imaging, Vision, and Graphics, S. Osher and N. Paragios, Eds., pp. 103–119, Springer, New York, NY, USA, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 G. Aubert and J. F. Aujol, “A variational approach to removing multiplicative noise,” SIAM Journal on Applied Mathematics, vol. 68, no. 4, pp. 925–946, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. M. Huang, M. K. Ng, and Y. W. Wen, “A new total variation method for multiplicative noise removal,” SIAM Journal on Imaging Sciences, vol. 2, no. 1, pp. 20–40, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 B. Chen, J. L. Cai, W. S. Chen, and Y. Li, “A multiplicative noise removal approach based on partial differential equation model,” Mathematical Problems in Engineering, vol. 2012, Article ID 242043, 14 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Correia, M. Carbillet, P. Boccacci, M. Bertero, and L. Fini, “Restoration of interferometric images—I. The software package AIRY,” Astronomy and Astrophysics, vol. 387, no. 2, pp. 733–743, 2002. View at: Publisher Site  Google Scholar
 A. Chambolle, “An algorithm for total variation minimization and applications,” Journal of Mathematical Imaging and Vision, vol. 20, no. 12, pp. 89–97, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 T. Chan, S. Esedoglu, and F. Park, “Image decomposition combining staircase reduction and texture extraction,” Journal of Visual Communication and Image Representation, vol. 18, no. 6, pp. 464–486, 2007. View at: Publisher Site  Google Scholar
 S. Becker, J. Bobin, and E. J. Candès, “NESTA: a fast and accurate firstorder method for sparse recovery,” SIAM Journal on Imaging Sciences, vol. 4, no. 1, pp. 1–39, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Beck and M. Teboulle, “Fast gradientbased algorithms for constrained total variation image denoising and deblurring problems,” IEEE Transactions on Image Processing, vol. 18, no. 11, pp. 2419–2434, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Y. L. Wang, J. F. Yang, W. T. Yin, and Y. Zhang, “A new alternating minimization algorithm for total variation image reconstruction,” SIAM Journal on Imaging Sciences, vol. 1, no. 3, pp. 248–272, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. BioucasDias and M. A. T. Figueiredo, “Multiplicative noise removal using variable splitting and constrained optimization,” IEEE Transactions on Image Processing, vol. 19, no. 7, pp. 1720–1730, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S. Yun and H. Woo, “A new multiplicative denoising variational model based on $m$th root transformation,” IEEE Transactions on Image Processing, vol. 21, no. 5, pp. 2523–2533, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 H. Woo and S. Yun, “Alternating minimization algorithm for speckle reduction with a shifting technique,” IEEE Transactions on Image Processing, vol. 21, no. 4, pp. 1701–1714, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 M. Bertero, P. Boccacci, G. Desiderà, and G. Vicidomini, “Image deblurring with Poisson data: from cells to galaxies,” Inverse Problems, vol. 25, no. 12, Article ID 123006, 26 pages, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. Bardsley and A. Luttman, “Total variationpenalized Poisson likelihood estimation for illposed problems,” Advances in Computational Mathematics, vol. 31, no. 13, pp. 35–59, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. T. Figueiredo and J. M. BioucasDias, “Restoration of Poissonian images using alternating direction optimization,” IEEE Transactions on Image Processing, vol. 19, no. 12, pp. 3133–3145, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S. Setzer, G. Steidl, and T. Teuber, “Deblurring Poissonian images by split Bregman techniques,” Journal of Visual Communication and Image Representation, vol. 21, no. 3, pp. 193–199, 2010. View at: Publisher Site  Google Scholar
 D. Q. Chen and L. Z. Cheng, “Deconvolving Poissonian images by a novel hybrid variational model,” Journal of Visual Communication and Image Representation, vol. 22, no. 7, pp. 643–652, 2011. View at: Publisher Site  Google Scholar
 I. Csiszár, “Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems,” The Annals of Statistics, vol. 19, no. 4, pp. 2032–2066, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. Bardsley and J. Goldes, “Techniques for regularization parameter and hyperparameter selection in PET and SPECT imaging,” Inverse Problems in Science and Engineering, vol. 19, no. 2, pp. 267–280, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. Bardsley and J. Goldes, “Regularization parameter selection methods for illposed Poisson maximum likelihood estimation,” Inverse Problems, vol. 25, no. 9, Article ID 095005, 18 pages, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Q. Chen and L. Z. Cheng, “Spatially adapted regularization parameter selection based on the local discrepancy function for Poissonian image deblurring,” Inverse Problems, vol. 28, no. 1, Article ID 015004, 24 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Landi and E. L. Piccolomini, “NPTool: a Matlab software for nonnegative image restoration with Newton projection methods,” Numerical Algorithms, 2012. View at: Publisher Site  Google Scholar
 A. Tikhonov and V. Arsenin, Solution of IllPoised Problems, Winston, Washington, DC, USA, 1977.
 G. Landi and E. L. Piccolomini, “An improved Newton projection method for nonnegative deblurring of Poissoncorrupted images with Tikhonov regularization,” Numerical Algorithms, vol. 60, no. 1, pp. 169–188, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. Bardsley and C. R. Vogel, “A nonnegatively constrained convex programming method for image reconstruction,” SIAM Journal on Scientific Computing, vol. 25, no. 4, pp. 1326–1343, 2003/04. View at: Publisher Site  Google Scholar  MathSciNet
 J. M. Bardsley and N. Laobeul, “Tikhonov regularized Poisson likelihood estimation: theoretical justification and a computational method,” Inverse Problems in Science and Engineering, vol. 16, no. 2, pp. 199–215, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 W. H. Richardson, “Bayesanbased iterative methods of image restoration,” Journal of the Optical Society of America, vol. 62, no. 1, pp. 55–59, 1972. View at: Publisher Site  Google Scholar
 L. B. Lucy, “An iterative technique for the rectification of observed images,” The Astronomical Journal, vol. 79, pp. 745–754, 1974. View at: Publisher Site  Google Scholar
 V. Agarwal, A. V. Gribok, and M. A. Abidi, “Image restoration using ${L}_{1}$ norm penalty function,” Inverse Problems in Science and Engineering, vol. 15, no. 8, pp. 785–809, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 G. Landi and E. L. Piccolomini, “An efficient method for nonnegatively constrained total variationbased denoising of medical images corrupted by Poisson noise,” Computerized Medical Imaging and Graphics, vol. 36, no. 1, pp. 38–46, 2012. View at: Publisher Site  Google Scholar
 S. Lefkimmiatis, A. Bourquard, and M. Unser, “Hessianbased norm regularization for image restoration with biomedical applications,” IEEE Transactions on Image Processing, vol. 21, no. 3, pp. 983–995, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 F. Li, C. M. Shen, J. S. Fan, and C. L. Shen, “Image restoration combining a total variational filter and a fourthorder filter,” Journal of Visual Communication and Image Representation, vol. 18, no. 4, pp. 322–330, 2007. View at: Publisher Site  Google Scholar
 M. Lysker and X. C. Tai, “Iterative image restoration combining total varition minimization and a secondorder functional,” International Journal of Computer Vision, vol. 66, no. 1, pp. 5–18, 2006. View at: Publisher Site  Google Scholar
 K. Papafitsoros and C. B. Schönlieb, “A combined first and second order variational approach for image reconstruction,” http://arxiv.org/abs/1202.6341. View at: Google Scholar
 T. Chan, A. Marquina, and P. Mulet, “Highorder total variationbased image restoration,” SIAM Journal on Scientific Computing, vol. 22, no. 2, pp. 503–516, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Lysaker, A. Lundervold, and X. C. Tai, “Noise removal using fourthorder partial differential equation with applications to medical magnetic resonance images in space and time,” IEEE Transactions on Image Processing, vol. 12, no. 12, pp. 1579–1590, 2003. View at: Google Scholar
 G. Steidl, “A note on the dual treatment of higherorder regularization functionals,” Computing, vol. 76, no. 12, pp. 135–148, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Z. Chen, J. P. Song, and X. C. Tai, “A dual algorithm for minimization of the LLT model,” Advances in Computational Mathematics, vol. 31, no. 1–3, pp. 115–130, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Le, R. Chartrand, and T. J. Asaki, “A variational approach to reconstructing images corrupted by Poisson noise,” Journal of Mathematical Imaging and Vision, vol. 27, no. 3, pp. 257–263, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 W. F. Zhou and Q. G. Li, “Poisson noise removal scheme based on fourthorder PDE by alternating minimization algorithm,” Abstract and Applied Analysis, vol. 2012, Article ID 965281, 14 pages, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Wang, A. C. Bovik, H. R. Sheikh, and E. P. Simoncelli, “Image quality assessment: from error visibility to structural similarity,” IEEE Transactions on Image Processing, vol. 13, no. 4, pp. 600–612, 2004. View at: Publisher Site  Google Scholar
 J. G. Nagy, K. Palmer, and L. Perrone, “Iterative methods for image deblurring: a Matlab objectoriented approach,” Numerical Algorithms, vol. 36, no. 1, pp. 73–93, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. C. Hansen, J. G. Nagy, and D. P. O'Leary, Deblurring Images: Matrices, Spectra, and Filtering, vol. 3 of Fundamentals of Algorithms, SIAM, Philadelphia, Pa, USA, 2006. View at: Publisher Site  MathSciNet
Copyright
Copyright © 2013 Le Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.