# Blind Poissonian Image Deblurring Regularized by a Denoiser Constraint and Deep Image Prior

Authors: Yayuan Feng; Yu Shi; Dianjun Sun
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The denoising and deblurring of Poisson images are opposite inverse problems. Single image deblurring methods are sensitive to image noise. A single noise filter can effectively remove noise in advance, but it also damages blurred information. To simultaneously solve the denoising and deblurring of Poissonian images better, we learn the implicit deep image prior from a single degraded image and use the denoiser as a regularization term to constrain the latent clear image. Combined with the explicitL0 regularization prior of the image, the denoising and deblurring model of the Poisson image is established. Then, the split Bregman iteration strategy is used to optimize the point spread function estimation and latent clear image estimation. The experimental results demonstrate that the proposed method achieves good restoration results on a series of simulated and real blurred images with Poisson noise.


---

## Body

## 1. Introduction

The process of image acquisition in electron microscope imaging, astronomical imaging, and medical imaging is inevitably affected by environmental factors, which often cause the captured image to be disturbed by Poisson noise and blur degradation. The image degradation process can be modeled as follows:(1)y=PH∗x,where y is the observed degraded image, P represents the process in which the image is corrupted by Poisson noise, H represents the point spread function, and ∗ is the convolution operator. The problem of blind Poisson image deconvolution is to recover clear images and estimate the point spread function from the degraded image. The denoising and deblurring of Poisson images are opposite inverse problems. Single deblurring methods are sensitive to noise, and even a small amount of noise leads to biased point spread function estimation. Although the noise filter can effectively remove noise in advance, it also damages the blurred information and introduces more serious blur, which leads to amplified and biased point spread function estimation. To suppress Poisson noise and restore the blurred image simultaneously, we impose the denoiser as a regularization term to constrain the latent clear image, combined with an implicit deep image prior and an explicit L0 regularization prior of the image. Then, the energy functional regularization framework for the latent clear image and point spread function is modeled as follows:(2)Ex,H=minx,H1,Hx−ylogHx+Rx+KH,where 1,Hx−ylogHx is the data term, 1 represents a vector with all elements equal to 1, x is the latent clear image, which usually ensures that the clear image does not contain negative gray values, H is the point spread function, Rx is the regularization constraint term of the latent clear image, and KH is the regularization constraint term of a point spread function. Among them, regularization constraint terms Rx and KH play a vital role in the restoration of latent clear images and the estimation of point spread functions.


### 1.1. Relevant Work

Poisson image denoising and blurring is an image deconvolution problem. Image noise can be divided into three categories: additive noise, multiplicative noise, and Poisson noise. The logarithm of the multiplicative noise can be processed as additive noise of Gaussian distribution. Different from the first two kinds of noise, Poisson noise generally appears in the case of very small illuminance and amplification with high power electronic circuits and obeys a Poisson distribution. Using variance stabilizing transformation [1–5] and Anscombe transformation [6–9], the transformed data can be regarded as additive noise of Gaussian distribution. However, such methods may cause some data loss during the conversion process. Thus, based on these transformations, the denoising and deblurring results need to be further improved. In recent years, many scholars have proposed algorithm models for Poissonian image deblurring. Poisson log-likelihood is inseparable in Poisson image deblurring. Many scholars [10–15] transform this problem into an equivalent constrained optimization problem through the split Bregman method and then use the augmented Lagrange method to solve this constrained optimization problem. Setzer et al. [12] developed a Poisson denoising model based on minimizing the total variation regularization term (TV) and the I-divergence [16] as a similarity term to restore blurred images. The TV regularization term is sensitive to the regularization parameter setting; thus, Yan et al. [17] proposed an improved spatially adaptive total variation regularization algorithm, which can automatically balance the regularization strength between different regions and retain more edge information in restored images, to solve this problem. Most TV regularization-based restoration models can effectively reduce the noise in flat regions, but large staircase effects are also introduced in flat regions, and fine details are not preserved in complex structural regions, which limits the practical application of the TV regularization term. To solve this problem, Fang et al. [13] developed a Poisson image deconvolution method based on framelet regularization, which aims to adaptively capture multiscale edge structures in images. To make full use of the sparse information and nonlocal information of the image, Shi et al. [18] proposed a latent image estimation method based on nonlocal total variation and framelet regularization constraints, which achieved smooth denoising and deblurring results while maintaining details and edges.Generally, images corrupted by Poisson noise are accompanied by blur degradation. Therefore, the deblurring problem of Poisson images is not only a denoising problem but also a deblurring problem. In the process of image restoration, the prior information of degraded image plays a very important role [17–28]. In recent years, many scholars have used regularization priors to denoise or deblur images, such as dictionary learning [19], total variation [17, 18, 21, 22], nonlocal mean [26–28], and L0 regularization prior [20]. In the field of image denoising, a nonlocal mean regularization term was proposed by Buades et al. [26, 27] in 2005 to remove image noise. Almahdi and Hardie [28] developed a new recursive nonlocal mean denoising algorithm. In 2017, Romano et al. [29] provided a new idea for removing noise; they provided an alternative, stronger, and more flexible framework for removing additive white Gaussian noise by proposing an explicit image-adaptive Laplacian-based regularization function, which can use the selected denoiser to define the regularization term, called Regularization by Denoising (RED). The advantage of RED is that it can flexibly select the denoising engine and use existing denoising algorithms to define regularization terms. When using RED to restore degraded images, the authors mainly focus on images contaminated by additive white Gaussian noise. They tested and proved that RED can achieve good restoration results as a regularization term to remove additive white Gaussian noise. Different from the image corrupted by Gaussian noise, the image corrupted by Poisson noise causes low image contrast and the image details to be covered. To make full use of the texture region information in the image and preserve the details of the image to the greatest extent while removing Poisson noise, we select the nonlocal mean denoiser as the regularization term of Poisson image denoising in the proposed method in this paper. In the field of image deblurring, Tang et al. [21] proposed a nonblind image deblurring method by local and nonlocal total variation models. Gradient priors are commonly used for image deblurring. However, in practical applications, if only the gradient prior is used, degraded images cannot be restored well. In [20], Pan et al. proposed an image deblurring method based on the L0 regularization intensity and gradient prior. However, these deblurring methods all assume that noise is additive white Gaussian noise or impulse noise. To recover clear images from Poisson images effectively, we introduce the intensity prior of L0 regularization and the prior of RED as explicit prior regularization terms in the proposed method.With the success of deep networks in image processing, an increasing number of scholars have begun to use neural networks to solve the restoration problem of degraded images [22, 30–37]. Most scholars usually need the model to train a large number of data sample sets when using the deep CNN network to restore degraded images so that the model can learn the prior information of these data, thus restoring degraded images. Therefore, when using a neural network for image restoration, a large number of data sample sets need to be prepared first, and then the degraded image can be restored based on the trained model parameters. However, Ulyanov et al. [33] found in their research that the deep network structure is sufficient for capturing a large number of low-level image statistics before any learning. In other words, the only prior information when restoring degraded images is derived from the structure of the network. Based on this discovery, they proposed a method of learning implicit priors from a single degraded image by using a network, which is called DIP. The authors noted that an untrained CNN network was used to restore degraded images, while only a randomly initialized network was required. Inspired by this idea, Mataev et al. [32] proposed bringing in the concept of Regularization by Denoising (RED). They boost DIP by adding an explicit prior, which enriches the overall regularization effect. In the same year, Liu et al. [22] proposed combining implicit priors in DIP with a traditional TV regularization prior to improve the image quality of denoising or deblurring in DIP. However, these methods are all aimed at restoring degraded images contaminated by additive white Gaussian noise. It is still a challenge to denoise and deblur Poisson images using neural networks. At present, all blind Poisson image deconvolution methods have been proposed to use traditional explicit regularization methods to restore degraded images. Therefore, in this paper, we propose an algorithm for denoising and deblurring Poisson images by using neural networks. We combine an implicit regularization prior with two explicit regularization priors that are the prior of RED and the prior of L0 regularization based on intensity and use a neural network to restore degraded images from Poisson images. Experimental results demonstrate that, compared with traditional methods, Poisson images restored by deep networks can effectively suppress Poisson noise while preserving the real image edge details and detailed texture information of degraded images.


## 1.1. Relevant Work

Poisson image denoising and blurring is an image deconvolution problem. Image noise can be divided into three categories: additive noise, multiplicative noise, and Poisson noise. The logarithm of the multiplicative noise can be processed as additive noise of Gaussian distribution. Different from the first two kinds of noise, Poisson noise generally appears in the case of very small illuminance and amplification with high power electronic circuits and obeys a Poisson distribution. Using variance stabilizing transformation [1–5] and Anscombe transformation [6–9], the transformed data can be regarded as additive noise of Gaussian distribution. However, such methods may cause some data loss during the conversion process. Thus, based on these transformations, the denoising and deblurring results need to be further improved. In recent years, many scholars have proposed algorithm models for Poissonian image deblurring. Poisson log-likelihood is inseparable in Poisson image deblurring. Many scholars [10–15] transform this problem into an equivalent constrained optimization problem through the split Bregman method and then use the augmented Lagrange method to solve this constrained optimization problem. Setzer et al. [12] developed a Poisson denoising model based on minimizing the total variation regularization term (TV) and the I-divergence [16] as a similarity term to restore blurred images. The TV regularization term is sensitive to the regularization parameter setting; thus, Yan et al. [17] proposed an improved spatially adaptive total variation regularization algorithm, which can automatically balance the regularization strength between different regions and retain more edge information in restored images, to solve this problem. Most TV regularization-based restoration models can effectively reduce the noise in flat regions, but large staircase effects are also introduced in flat regions, and fine details are not preserved in complex structural regions, which limits the practical application of the TV regularization term. To solve this problem, Fang et al. [13] developed a Poisson image deconvolution method based on framelet regularization, which aims to adaptively capture multiscale edge structures in images. To make full use of the sparse information and nonlocal information of the image, Shi et al. [18] proposed a latent image estimation method based on nonlocal total variation and framelet regularization constraints, which achieved smooth denoising and deblurring results while maintaining details and edges.Generally, images corrupted by Poisson noise are accompanied by blur degradation. Therefore, the deblurring problem of Poisson images is not only a denoising problem but also a deblurring problem. In the process of image restoration, the prior information of degraded image plays a very important role [17–28]. In recent years, many scholars have used regularization priors to denoise or deblur images, such as dictionary learning [19], total variation [17, 18, 21, 22], nonlocal mean [26–28], and L0 regularization prior [20]. In the field of image denoising, a nonlocal mean regularization term was proposed by Buades et al. [26, 27] in 2005 to remove image noise. Almahdi and Hardie [28] developed a new recursive nonlocal mean denoising algorithm. In 2017, Romano et al. [29] provided a new idea for removing noise; they provided an alternative, stronger, and more flexible framework for removing additive white Gaussian noise by proposing an explicit image-adaptive Laplacian-based regularization function, which can use the selected denoiser to define the regularization term, called Regularization by Denoising (RED). The advantage of RED is that it can flexibly select the denoising engine and use existing denoising algorithms to define regularization terms. When using RED to restore degraded images, the authors mainly focus on images contaminated by additive white Gaussian noise. They tested and proved that RED can achieve good restoration results as a regularization term to remove additive white Gaussian noise. Different from the image corrupted by Gaussian noise, the image corrupted by Poisson noise causes low image contrast and the image details to be covered. To make full use of the texture region information in the image and preserve the details of the image to the greatest extent while removing Poisson noise, we select the nonlocal mean denoiser as the regularization term of Poisson image denoising in the proposed method in this paper. In the field of image deblurring, Tang et al. [21] proposed a nonblind image deblurring method by local and nonlocal total variation models. Gradient priors are commonly used for image deblurring. However, in practical applications, if only the gradient prior is used, degraded images cannot be restored well. In [20], Pan et al. proposed an image deblurring method based on the L0 regularization intensity and gradient prior. However, these deblurring methods all assume that noise is additive white Gaussian noise or impulse noise. To recover clear images from Poisson images effectively, we introduce the intensity prior of L0 regularization and the prior of RED as explicit prior regularization terms in the proposed method.With the success of deep networks in image processing, an increasing number of scholars have begun to use neural networks to solve the restoration problem of degraded images [22, 30–37]. Most scholars usually need the model to train a large number of data sample sets when using the deep CNN network to restore degraded images so that the model can learn the prior information of these data, thus restoring degraded images. Therefore, when using a neural network for image restoration, a large number of data sample sets need to be prepared first, and then the degraded image can be restored based on the trained model parameters. However, Ulyanov et al. [33] found in their research that the deep network structure is sufficient for capturing a large number of low-level image statistics before any learning. In other words, the only prior information when restoring degraded images is derived from the structure of the network. Based on this discovery, they proposed a method of learning implicit priors from a single degraded image by using a network, which is called DIP. The authors noted that an untrained CNN network was used to restore degraded images, while only a randomly initialized network was required. Inspired by this idea, Mataev et al. [32] proposed bringing in the concept of Regularization by Denoising (RED). They boost DIP by adding an explicit prior, which enriches the overall regularization effect. In the same year, Liu et al. [22] proposed combining implicit priors in DIP with a traditional TV regularization prior to improve the image quality of denoising or deblurring in DIP. However, these methods are all aimed at restoring degraded images contaminated by additive white Gaussian noise. It is still a challenge to denoise and deblur Poisson images using neural networks. At present, all blind Poisson image deconvolution methods have been proposed to use traditional explicit regularization methods to restore degraded images. Therefore, in this paper, we propose an algorithm for denoising and deblurring Poisson images by using neural networks. We combine an implicit regularization prior with two explicit regularization priors that are the prior of RED and the prior of L0 regularization based on intensity and use a neural network to restore degraded images from Poisson images. Experimental results demonstrate that, compared with traditional methods, Poisson images restored by deep networks can effectively suppress Poisson noise while preserving the real image edge details and detailed texture information of degraded images.


## 2. Proposed methodology

### 2.1. RED

The regularization term plays an important role in image deconvolution. Compared with other regularization terms, RED [29] has the advantage that existing denoising algorithms can be flexibly used to define regularization terms. The regularization term by denoising, which is an explicit image-adaptive Laplacian-based regularization function, uses the selected denoiser to define the regularization term. RED, as a regularization function, is given by the following formula:(3)ρx=12xTx−fx,where x is a degraded image, f· is the selected denoising engine, which is applied to degraded image x, and the selection of the denoising engine is flexible. x−fx is the denoising residual. We found in the experimental research that the nonlocal mean denoising device can make full use of the texture region information in the degraded image, which can remove the Poisson noise while maintaining the details of the image features to the greatest extent. Therefore, we use the nonlocal mean denoiser as the regularization term of Poisson image denoising in the proposed method in this paper.


### 2.2.L0 Regularization Term

Based on the observation of blurred images [20], the pixel intensity histogram of a blurred image is different from that of a clean image. In the experimental observation, it is found that the number of nonzero elements of pixel intensity in blurred images is large and the distribution is relatively dense. According to the pixel intensity property of blurred images, for blurred images, we have(4)Px=x0,where x0 represents the number of nonzero values in image x, and clear images and blurred images can be distinguished according to the criterion of the pixel value intensity distribution. Using this property of pixel value intensity distribution, clear images and blurred images can be distinguished in the process of image restoration, thus speeding up the convergence time of the algorithm. Therefore, the L0 regularization term based on the pixel intensity prior is adopted as the regularization prior term in the proposed method in this paper. The result analysis of convergence of L0 regularization is described in detail in Section 3.2.


### 2.3. DIP

In [33], the authors determined that the process of recovering degraded images in DIP is different from other CNN networks. It does not require a large amount of data-driven training but only requires the network to learn the implicit prior of degraded images in a randomly initialized neural network framework. The authors found that the neural network can learn the undamaged part of the image and then learn the damaged part of the image. For example, if the degraded image of white Gaussian noise is input to the network, the network will learn how to copy an image without noise first and then learn to copy the noise. Therefore, we can remove the network learning before the network learns to copy the noise so that we can obtain an image without noise. DIP defines the output of the randomly initialized neural network as(5)x=Fθz,where θ represents the parameters of the network, z is a fixed random vector, and Fθz represents the parameterized form of the network. Then, the objective function of DIP is as follows:(6)minθHFθz−x02,where H is a degraded operator and x0 is a degraded image. The network consists of linear convolution, upsampling, and nonlinear activation functions. Based on the given degraded image and the observation model, the model parameters can be approximated to the maximum likelihood through iteration.


### 2.4. Proposed Model

In equation (6), to restore degraded images through a deep network, data items for restoring additive white Gaussian noise images cannot solve the problem of Poisson image deblurring well. To obtain good denoising results while preserving the image details, we introduce equation (6) into the objective function 2 for Poissonian image deblurring and use RED as the denoising regularization term in the Poisson image deblurring model. We take the L0 norm of the latent clear image x as the regularization constraint in the Poisson image deblurring model. Then, the objective function 2 can be written as follows:(7)Ex,H=minx,H1,Hx−ylogHx+λ2xTx−fx+β2x0+τ2H22,where λ is the weight of RED and β is the weight of the L0 norm constrained regularization term.To solve the inseparable Poisson log-likelihood problem in equation (7), we use the split Bregman method to introduce auxiliary variable d1=Hx and transform equation (7) into a constrained problem. In combination with equation (5), if the output x=Fθz of the neural network is introduced into equation (7) as a constraint condition, equation (7) becomes the following equation:(8)mind1,x,H1,d1−ylogd1+λ2xTx−fx+β2x0+τ2H22s.t.d1=Hx,x=Fθz.


### 2.5. Optimal Estimation

To avoid the differentiation of the explicit denoising function, two auxiliary variables,u1 and u2, are introduced by using the split Bregman method. The existence of the L0 norm in the third term of equation (8) makes the solution somewhat difficult, the split Bregman method is used to introduce auxiliary variable v, and equation (8) is changed to the following equation:(9)mind1,x,θ,H1,d1−ylogd1+μ2d1−HFθz−u122+γ2x−Fθz−u222+λ2xTx−fx+α2x−v22+β2v0+τ2H22.According to the split Bregman method, seven unknown variablesd1, θ, x, u1, u2, v, and H are iteratively optimized and updated. Variables d1, u1, x, and H are fixed, and variable θ can be solved by the following formula:(10)minθμ2d1−HFθz−u122+γ2x−Fθz−u222.To solve equation (10), we use backpropagation to optimize update variable θ, where variables u1 and u2 are infinitely close to d1−HFθz and x−Fθz, respectively.Given fixed variablesd1, u1, and θ, the point spread function H can be solved by the following formula:(11)minHμ2d1−HFθz−u122+τ2H22.The update of variableHk+1 is obtained by the following equation:(12)Hk+1=μd1k−u1kTFθzk+1μFθzk+1Fθzk+1T+τ.Given fixed variablesθ, u1, and H, variable d1 can be solved by the following formula:(13)mind11,d1−ylogd1+μ2d1−HFθz−u122.The update of variabled1k+1, which can be solved by equation (13), is obtained by the following equation:(14)d1k+1=12μμHk+1Fθzk+1+μu1k−1+μHk+1Fθzk+1+μu1k−12+4μy,where d1k+1 is the value of d1 in the k+1th iteration and Fθzk+1 is the output of the neural network in the k+1th iteration.Given fixed variablesθ, u2, and v, variable x can be solved by the following formula:(15)minxγ2x−Fθz−u222+λ2xTx−fx+α2x−v22.We use gradient descent to updatexk+1:(16)xk+1=xk−cγxk−Fθzk+1−u2+λxk−fxk+αxk−v,where c should be selected to ensure a decrease.Given fixed variablesd1, H, and θ, auxiliary variable u1k+1 can be updated by the following formula:(17)u1k+1=u1k−d1k+1+Hk+1Fθzk+1.Given fixed variablesx and θ, auxiliary variable u2k+1 can be updated by the following formula:(18)u2k+1=u2k−xk+1+Fθzk+1.Given fixed variablex, auxiliary variable vk+1 can be updated by the following formula:(19)vk+1=xk+1,xk+12≥βα,0,otherwise.Algorithm1 proposed in this paper is as follows:Algorithm 1: Proposed method.
Input: degraded imageyInitialization:k = 0, H0=25, u10 = 0, u20 = 0, v0=0, d10=0, x0=y, randomly θ0Whileθk+1 is not converged doUpdateθk+1 by equation (10);UpdateHk+1 by equation (12);Updated1k+1 by equation (14);Updatexk+1 by equation (16);Updateu1k+1 by equation (17);Updateu2k+1 by equation (18);Updatevk+1 by equation (19);EndOutput: recovered clear imagex


## 2.1. RED

The regularization term plays an important role in image deconvolution. Compared with other regularization terms, RED [29] has the advantage that existing denoising algorithms can be flexibly used to define regularization terms. The regularization term by denoising, which is an explicit image-adaptive Laplacian-based regularization function, uses the selected denoiser to define the regularization term. RED, as a regularization function, is given by the following formula:(3)ρx=12xTx−fx,where x is a degraded image, f· is the selected denoising engine, which is applied to degraded image x, and the selection of the denoising engine is flexible. x−fx is the denoising residual. We found in the experimental research that the nonlocal mean denoising device can make full use of the texture region information in the degraded image, which can remove the Poisson noise while maintaining the details of the image features to the greatest extent. Therefore, we use the nonlocal mean denoiser as the regularization term of Poisson image denoising in the proposed method in this paper.


## 2.2.L0 Regularization Term

Based on the observation of blurred images [20], the pixel intensity histogram of a blurred image is different from that of a clean image. In the experimental observation, it is found that the number of nonzero elements of pixel intensity in blurred images is large and the distribution is relatively dense. According to the pixel intensity property of blurred images, for blurred images, we have(4)Px=x0,where x0 represents the number of nonzero values in image x, and clear images and blurred images can be distinguished according to the criterion of the pixel value intensity distribution. Using this property of pixel value intensity distribution, clear images and blurred images can be distinguished in the process of image restoration, thus speeding up the convergence time of the algorithm. Therefore, the L0 regularization term based on the pixel intensity prior is adopted as the regularization prior term in the proposed method in this paper. The result analysis of convergence of L0 regularization is described in detail in Section 3.2.


## 2.3. DIP

In [33], the authors determined that the process of recovering degraded images in DIP is different from other CNN networks. It does not require a large amount of data-driven training but only requires the network to learn the implicit prior of degraded images in a randomly initialized neural network framework. The authors found that the neural network can learn the undamaged part of the image and then learn the damaged part of the image. For example, if the degraded image of white Gaussian noise is input to the network, the network will learn how to copy an image without noise first and then learn to copy the noise. Therefore, we can remove the network learning before the network learns to copy the noise so that we can obtain an image without noise. DIP defines the output of the randomly initialized neural network as(5)x=Fθz,where θ represents the parameters of the network, z is a fixed random vector, and Fθz represents the parameterized form of the network. Then, the objective function of DIP is as follows:(6)minθHFθz−x02,where H is a degraded operator and x0 is a degraded image. The network consists of linear convolution, upsampling, and nonlinear activation functions. Based on the given degraded image and the observation model, the model parameters can be approximated to the maximum likelihood through iteration.


## 2.4. Proposed Model

In equation (6), to restore degraded images through a deep network, data items for restoring additive white Gaussian noise images cannot solve the problem of Poisson image deblurring well. To obtain good denoising results while preserving the image details, we introduce equation (6) into the objective function 2 for Poissonian image deblurring and use RED as the denoising regularization term in the Poisson image deblurring model. We take the L0 norm of the latent clear image x as the regularization constraint in the Poisson image deblurring model. Then, the objective function 2 can be written as follows:(7)Ex,H=minx,H1,Hx−ylogHx+λ2xTx−fx+β2x0+τ2H22,where λ is the weight of RED and β is the weight of the L0 norm constrained regularization term.To solve the inseparable Poisson log-likelihood problem in equation (7), we use the split Bregman method to introduce auxiliary variable d1=Hx and transform equation (7) into a constrained problem. In combination with equation (5), if the output x=Fθz of the neural network is introduced into equation (7) as a constraint condition, equation (7) becomes the following equation:(8)mind1,x,H1,d1−ylogd1+λ2xTx−fx+β2x0+τ2H22s.t.d1=Hx,x=Fθz.


## 2.5. Optimal Estimation

To avoid the differentiation of the explicit denoising function, two auxiliary variables,u1 and u2, are introduced by using the split Bregman method. The existence of the L0 norm in the third term of equation (8) makes the solution somewhat difficult, the split Bregman method is used to introduce auxiliary variable v, and equation (8) is changed to the following equation:(9)mind1,x,θ,H1,d1−ylogd1+μ2d1−HFθz−u122+γ2x−Fθz−u222+λ2xTx−fx+α2x−v22+β2v0+τ2H22.According to the split Bregman method, seven unknown variablesd1, θ, x, u1, u2, v, and H are iteratively optimized and updated. Variables d1, u1, x, and H are fixed, and variable θ can be solved by the following formula:(10)minθμ2d1−HFθz−u122+γ2x−Fθz−u222.To solve equation (10), we use backpropagation to optimize update variable θ, where variables u1 and u2 are infinitely close to d1−HFθz and x−Fθz, respectively.Given fixed variablesd1, u1, and θ, the point spread function H can be solved by the following formula:(11)minHμ2d1−HFθz−u122+τ2H22.The update of variableHk+1 is obtained by the following equation:(12)Hk+1=μd1k−u1kTFθzk+1μFθzk+1Fθzk+1T+τ.Given fixed variablesθ, u1, and H, variable d1 can be solved by the following formula:(13)mind11,d1−ylogd1+μ2d1−HFθz−u122.The update of variabled1k+1, which can be solved by equation (13), is obtained by the following equation:(14)d1k+1=12μμHk+1Fθzk+1+μu1k−1+μHk+1Fθzk+1+μu1k−12+4μy,where d1k+1 is the value of d1 in the k+1th iteration and Fθzk+1 is the output of the neural network in the k+1th iteration.Given fixed variablesθ, u2, and v, variable x can be solved by the following formula:(15)minxγ2x−Fθz−u222+λ2xTx−fx+α2x−v22.We use gradient descent to updatexk+1:(16)xk+1=xk−cγxk−Fθzk+1−u2+λxk−fxk+αxk−v,where c should be selected to ensure a decrease.Given fixed variablesd1, H, and θ, auxiliary variable u1k+1 can be updated by the following formula:(17)u1k+1=u1k−d1k+1+Hk+1Fθzk+1.Given fixed variablesx and θ, auxiliary variable u2k+1 can be updated by the following formula:(18)u2k+1=u2k−xk+1+Fθzk+1.Given fixed variablex, auxiliary variable vk+1 can be updated by the following formula:(19)vk+1=xk+1,xk+12≥βα,0,otherwise.Algorithm1 proposed in this paper is as follows:Algorithm 1: Proposed method.
Input: degraded imageyInitialization:k = 0, H0=25, u10 = 0, u20 = 0, v0=0, d10=0, x0=y, randomly θ0Whileθk+1 is not converged doUpdateθk+1 by equation (10);UpdateHk+1 by equation (12);Updated1k+1 by equation (14);Updatexk+1 by equation (16);Updateu1k+1 by equation (17);Updateu2k+1 by equation (18);Updatevk+1 by equation (19);EndOutput: recovered clear imagex


## 3. Experimental Results

In this section, we carry out experiments and image quality evaluation on simulated blurred images and real blurred images that are corrupted with Poisson noise and compare the experimental results with those of four methods: PID-Split [12], PIDSB-FA [13], PIDSB-NLFA [18], and DeepRED [32]. The experiment of the proposed method and the DeepRED method is implemented on a GTX 1080i GPU computer using the Python language in the Linux system. PID-Split, PIDSB-FA, and PIDSB-NLFA algorithms are tested by MATLAB 8a on a computer with an Intel(R) Core(TM) i5-7400 CPU and 8 GB RAM. In the section of parameter setting, we set the parameters of the PID-Split, PIDSB-FA, and PIDSB-NLFA methods according to the methods suggested in the original paper. For different Poisson degraded images, we set different parameters to obtain the best recovery results. The parameter settings in our proposed method are as follows: μ=0.01, γ=0.04, λ=0.02, α=0.01, β=0.002, τ=0.04. To evaluate the image quality, the Peak Signal-to-Noise Ratio (PSNR) and Visual Information Fidelity (VIF) are used in this section to evaluate the Poisson image recovered by each algorithm. Finally, the PSNR value of each algorithm is evaluated by using the “compare_psnr” function in the “skimage.measure” module in the Python language, and the VIF value of each algorithm is evaluated by using the “vifp_mscale” function in MATLAB. The higher the value of PSNR, the better the image recovery effect, the higher the value of VIF, and the higher the visual information fidelity of the image.


### 3.1. Simulated Poisson Image Denoising and Blurring

In this section, we compare the proposed method with the other four algorithms (PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED) on five images with simulated Poisson image denoising and blurring experiments. The five images are Lena (size = 256 × 256), Cameraman (size = 256 × 256), Butterfly (size = 256 × 256), Zebra (size = 584 × 387), and House (size = 256 × 256), as shown in Figure1.Figure 1
Original images. (a) Lena (size = 256 × 256). (b) Cameraman (size = 256 × 256). (c) Butterfly (size = 256 × 256). (d) Zebra (size = 584 × 387). (e) House (size = 256 × 256).
(a)(b)(c)(d)(e)In the simulated Poisson image experiment, we applied different levels of blur and noise to different images. The Gaussian blur kernel parameters we set for the Lena image of Figure1(a) were (size = 5, σ = 1.6), (size = 15, σ = 1.6), and (size = 25, σ = 2). Three different Gaussian blur kernels were used to convolve the remaining four images, and the Gaussian blur kernels were set to (size = 5, σ = 1.6), (size = 15, σ = 2), and (size = 25, σ = 1.6). Then, the “skimage.util.random_noise” function of the image processing algorithm set in Python was used to add three different Poisson noise levels to the five images, and the noise levels “Imax” were set to 3,000, 4,000, and 5,000. Particularly, since the Butterfly image with an image size of 256 × 256 in Figure 1(c) was insensitive to noise levels of 5,000 when the Gaussian blur kernel was (size = 15, σ = 2), we applied noise levels of 2,000, 3,000, and 4,000 to this image.In the experiment of deblurring simulated Poisson images with different noise levels and different blur kernel sizes, the image quality evaluation data recovered by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms and the proposed method in this paper are shown in Tables1 and 2. According to the image quality evaluation data in Tables 1 and 2, the PSNR value and VIF value of the proposed method in this paper were higher than those of the other four algorithms. We conducted a set of comparison experiments on the House image in Figure 1(e), and the results of the comparison experiments are shown in Figure 2. The House images have a large size of flat regions and a small amount of detailed information. In this group of comparative experiments, as shown in Figure 2(b), the House image was degraded by a Gaussian blur kernel of size 15 × 15 and standard deviation σ = 2 and polluted by Poisson noise (Imax = 3,000). As shown in Figure 2(c), although the image restored by the PID-Split algorithm suppressed noise well, it also caused staircase effects. The images recovered by the PIDSB-NLFA and PIDSB-FA algorithms are shown in Figures 2(d) and 2(e). Compared with Figure 2(c), the staircase effects were suppressed, and the noise in the flat regions was suppressed better. Figure 2(f) shows the result of restoration of the DeepRED, which retains more details but cannot remove Poisson noise well. The image restored by the proposed method in this paper is shown in Figure 2(g). Compared with the other four methods, it not only suppressed noise to the greatest extent but also had clearer details and the best visual perception.Table 1
PSNR (dB) values of the five methods.
ImageGaussian blurImaxDegraded imagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodLenaSize = 5,σ = 1.63,00024.8425.0625.2225.5725.8026.914,00024.8725.0725.2525.6225.8226.945,00024.9025.1125.2925.6825.8527.11Size = 15,σ = 1.63,00024.3625.0125.1225.2625.8126.334,00024.3825.0325.1625.2825.9226.475,00024.4025.0725.2025.2926.0926.57Size = 25,σ = 23,00022.9923.7024.4424.5524.5724.614,00023.0023.7224.4924.5724.5924.645,00023.0223.7524.5124.5924.6124.66CameramanSize = 5,σ = 1.63,00023.6325.3826.6826.7327.1227.464,00023.6625.5126.7326.7527.2028.055,00023.6725.6926.8526.9227.4028.14Size = 15,σ = 23,00021.7323.3924.0224.0924.2424.594,00021.7523.4324.1524.1724.3824.645,00021.7623.4824.2024.2324.4724.71Size = 25,σ = 1.63,00022.8324.3925.2925.3725.5126.294,00022.8424.4925.3225.4025.7426.305,00022.8824.6625.3525.4225.9226.50ButterflySize = 5,σ = 1.63,00022.8824.5025.2325.5225.5625.954,00022.9024.6325.2625.5525.6226.015,00022.9124.7825.2725.5825.6626.12Size = 15,σ = 22,00020.5721.9622.9923.2923.4723.983,00020.6022.0423.0323.4723.5924.054,00020.6322.2223.0723.5723.6924.13Size = 25,σ = 1.63,00021.8123.1723.9624.3324.4624.814,00021.8323.3323.9924.3624.5224.895,00021.8423.4824.0024.3924.5924.98ZebraSize = 5,σ = 1.63,00023.5824.8125.8926.5026.8727.174,00023.6124.8625.9626.5626.9627.245,00023.6424.9226.0226.6227.0527.30Size = 15,σ = 23,00021.6923.9525.0525.3226.3826.764,00021.7124.0325.1225.4326.5726.855,00021.7324.1125.2025.5426.6926.92Size = 25,σ = 1.63,00022.6824.4125.5126.0426.7827.224,00022.7124.4725.5726.1726.9027.355,00022.7324.5425.6526.2927.0327.42HouseSize = 5,σ = 1.63,00028.8130.0830.7130.4332.4432.714,00028.9330.1730.9130.4832.5432.885,00028.9930.2831.1030.6032.6133.18Size = 15,σ = 23,00026.1127.1227.6927.5228.3628.474,00026.1427.2027.7027.5728.3728.545,00026.1827.3327.7327.6228.3928.62Size = 25,σ = 1.63,00027.7128.7029.6129.4130.2930.874,00027.8428.7829.7329.4930.4731.085,00027.8528.8729.8129.5630.5631.15Table 2
VIF values of the five methods.
ImageGaussian blurImaxDegraded imagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodLenaSize = 5,σ = 1.63,0000.44820.43100.45130.51470.51820.57794,0000.45530.43270.45590.51220.51310.59255,0000.46030.43690.46140.51470.51360.5984Size = 15,σ = 1.63,0000.40260.43070.42980.49380.51780.55854,0000.40740.43380.43460.49420.52820.56835,0000.41160.43930.44040.49460.54100.5745Size = 25,σ = 23,0000.34840.39460.47610.49490.51300.52904,0000.35310.39730.48120.49620.51160.52975,0000.35660.40030.48420.49690.52900.5373CameramanSize = 5,σ = 1.63,0000.30860.40450.43320.43360.45130.45334,0000.31070.41330.44360.43950.45630.46755,0000.31400.42330.44500.44580.45730.4723Size = 15,σ = 23,0000.22550.31320.32600.38600.37430.39814,0000.22800.31660.33470.38670.38900.40175,0000.23010.31950.33720.38610.39500.4105Size = 25,σ = 1.63,0000.27050.37300.38540.38600.39960.42954,0000.27330.38050.38940.38760.40400.44075,0000.27600.39100.39030.38610.41850.4447ButterflySize = 5,σ = 1.63,0000.39400.47090.55350.55970.57680.58874,0000.39670.48150.55770.56140.59090.60285,0000.39820.49330.56020.56350.60740.6102Size = 15,σ = 22,0000.29250.38180.46220.47170.50180.51243,0000.29570.38770.46610.48490.52280.52844,0000.29910.40020.47130.49670.53350.5390Size = 25,σ = 1.63,0000.34710.41190.51170.52920.54780.56134,0000.34920.41890.51420.53250.54880.56305,0000.35070.42630.51630.53580.56600.5800ZebraSize = 5,σ = 1.63,0000.39130.40070.49050.51920.49130.52064,0000.39790.40380.49570.52780.51060.53025,0000.40200.4070.049920.51090.51790.5355Size = 15,σ = 23,0000.28320.33480.39820.40940.42950.44904,0000.28780.33870.40260.41550.43940.45905,0000.29110.34290.40710.42190.44480.4622Size = 25,σ = 1.63,0000.34020.36970.40660.47120.47630.49334,0000.34630.37310.41000.47890.48500.50165,0000.35040.37680.41360.48570.49590.5082HouseSize = 5,σ = 1.63,0000.38420.40280.44720.50400.50620.51564,0000.38980.40810.45630.50490.51030.51725,0000.39240.41270.46590.50470.51330.5320Size = 15,σ = 23,0000.29140.36200.4378043570.43340.45594,0000.29290.36880.43800.43730.43700.46805,0000.29460.37530.44040.43750.43890.4749Size = 25,σ = 1.63,0000.34290.37430.46530.47190.45220.49464,0000.34960.37830.47410.47240.46860.50375,0000.34950.38230.47520.46630.47030.5091Figure 2
Restoration of the House image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 15 × 15,σ = 2 and corrupted by Poisson noise with Imax = 3,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).Figure3 shows a set of comparative results for Zebra images in Figure 1(d). The size of Zebra image is 584 × 387. Figure 3(b) is an image degraded by a Gaussian blur kernel with a size of 25 × 25 and a standard deviation of 1.6 and Poisson noise with a size of 3,000. As can be seen from Figure 3(a), the zebra has a large size of flat regions and the grassland contains considerable detailed information. As seen in Figure 3(g), the proposed method in this paper not only suppressed noise but also preserved considerable detailed information. Figure 3(c) shows that the image restored by the PID-Split algorithm was too smooth and had staircase effects. Figures 3(d)–3(f) show the images recovered by the PIDSB-NLFA, PIDSB-FA, and DeepRED algorithms, respectively. Compared with Figure 3(c), they suppressed the staircase effects while retaining some details. Compared with the other four method, the proposed method in this paper not only suppressed the staircase effects but also preserved the details of the image to the greatest extent and was clearer.Figure 3
Restoration of the Zebra image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 25 × 25,σ = 1.6 and corrupted by Poisson noise with Imax = 3,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).In addition, the Butterfly image in Figure1(c) was degraded by the Gaussian blur kernel (size = 5 × 5, σ = 1.6) and then corrupted by Poisson noise (Imax = 5,000). As shown in Figure 4, compared with the images restored by the other four algorithms, the proposed method in this paper suppressed noise and restored clear images to the greatest extent. Figure 4(g) restored by our method looks natural and preserved the most detailed information of the images.Figure 4
Restoration of the Butterfly image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 5 × 5,σ = 1.6 and corrupted by Poisson noise with Imax = 5,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).


### 3.2. Convergence Analysis

The explicit regularization prior term adopts the combination of the intensity-basedL0 regularization term and RED. The purpose is to be able to use the pixel intensity distribution of the blurred image and the clear image to distinguish the restored image pixels while effectively filtering the image pixels and removing the Poisson noise of the degraded image. To verify the effectiveness of the L0 regularization prior, we used the method with the L0 norm regularization term and the method without the L0 norm to carry out simulation and comparison experiments on the House images of Figure 1(c) in Section 3.1. The experimental results are shown in Figure 5. Figure 5(a) shows the PSNR values of the image restored by the two algorithms. It can be clearly seen that the method with the L0 regularization term has the characteristic of fast convergence. With the increase in iteration times, the restored image reached a higher PSNR. Figure 5(b) shows the enlarged graph of the abscissa range (0∼500) in Figure 5(a). Figure 5(c) shows the loss of the objective function, and Figure 5(d) shows an enlarged graph with abscissa (0∼200) in Figure 5(c). It can be seen that the method with the L0 regularization term also achieved a smaller loss value.Figure 5
PSNR values and loss comparison. (a) PSNR. (b) PSNE (zoomed in). (c) Loss. (d) Loss (zoomed in).
(a)(b)(c)(d)


### 3.3. Denoising and Deblurring of Real Poisson Images

This section is an experimental comparison of denoising and deblurring of real Poisson images. Tests were carried out on four real astronomically degraded images. To evaluate the real Poisson image quality, we used the Blind/Referenceless Image Spatial Quality Evaluator (BRISQUE) to evaluate the real data. The smaller the value of BRISQUE, the better the image quality. The image quality evaluation data recovered by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms and the proposed method in this paper are shown in Table3.Table 3
BRISQUE values of the five methods.
ImagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodSaturn45.1343.3148.9641.4239.22Moon49.8143.3447.7539.2838.56Docking48.6549.1449.0541.9539.85Phoebe36.6132.9431.2630.4424.25Figure6(a) shows a satellite image, and Figure 6(b) shows an image restored by the PID-Split algorithm; the restored image is too smooth. Figures 6(c)–6(e) show the images recovered by the PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms, respectively. Although the staircase effects were suppressed to a large extent, ring effects exist at the edge. The image restored by the proposed method in this paper is shown in Figure 6(f). The restored image not only restored a large amount of detailed information and texture information but also suppressed the ring effects at the edge while suppressing the step effect. The image restored by the proposed method looks more natural.Figure 6
Restoration of real blurred image. (a) Saturn truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).The second real image is a moon image, as shown in Figure7(a). Figures 7(b)–7(e) show the images restored by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms, respectively. It can be seen from these three images that the images recovered by these methods had ringing effects of the edge, and the flat areas were too smooth. Figure 7(f) shows the image restored by the proposed method in this paper. It can be seen that the image restored by the proposed method in this paper retained more detailed information while suppressing noise and ringing effects.Figure 7
Restoration of real blurred image. (a) Moon truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).Figure8 shows a set of comparative results for the docking images of Shenzhou-9 and Tiangong-1. Restored images by RLTV, PID-Split, PIDSB-FA, DeepRED, and the proposed method are shown in Figures 8(b)–8(f), respectively. The image restored by the proposed method suppressed the ring effects at the edge and retained more detailed information. Figure 9(a) shows the Phoebe image. Comparing the restored images by the other four methods shown in Figures 9(b)–9(e), respectively, the image restored by the proposed method in this paper looks clear with fine details.Figure 8
Restoration of real blurred image. (a) The docking truth image of Shenzhou-9 and Tiangong-1. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).Figure 9
Restoration of real blurred image. (a) Phoebe truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).


## 3.1. Simulated Poisson Image Denoising and Blurring

In this section, we compare the proposed method with the other four algorithms (PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED) on five images with simulated Poisson image denoising and blurring experiments. The five images are Lena (size = 256 × 256), Cameraman (size = 256 × 256), Butterfly (size = 256 × 256), Zebra (size = 584 × 387), and House (size = 256 × 256), as shown in Figure1.Figure 1
Original images. (a) Lena (size = 256 × 256). (b) Cameraman (size = 256 × 256). (c) Butterfly (size = 256 × 256). (d) Zebra (size = 584 × 387). (e) House (size = 256 × 256).
(a)(b)(c)(d)(e)In the simulated Poisson image experiment, we applied different levels of blur and noise to different images. The Gaussian blur kernel parameters we set for the Lena image of Figure1(a) were (size = 5, σ = 1.6), (size = 15, σ = 1.6), and (size = 25, σ = 2). Three different Gaussian blur kernels were used to convolve the remaining four images, and the Gaussian blur kernels were set to (size = 5, σ = 1.6), (size = 15, σ = 2), and (size = 25, σ = 1.6). Then, the “skimage.util.random_noise” function of the image processing algorithm set in Python was used to add three different Poisson noise levels to the five images, and the noise levels “Imax” were set to 3,000, 4,000, and 5,000. Particularly, since the Butterfly image with an image size of 256 × 256 in Figure 1(c) was insensitive to noise levels of 5,000 when the Gaussian blur kernel was (size = 15, σ = 2), we applied noise levels of 2,000, 3,000, and 4,000 to this image.In the experiment of deblurring simulated Poisson images with different noise levels and different blur kernel sizes, the image quality evaluation data recovered by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms and the proposed method in this paper are shown in Tables1 and 2. According to the image quality evaluation data in Tables 1 and 2, the PSNR value and VIF value of the proposed method in this paper were higher than those of the other four algorithms. We conducted a set of comparison experiments on the House image in Figure 1(e), and the results of the comparison experiments are shown in Figure 2. The House images have a large size of flat regions and a small amount of detailed information. In this group of comparative experiments, as shown in Figure 2(b), the House image was degraded by a Gaussian blur kernel of size 15 × 15 and standard deviation σ = 2 and polluted by Poisson noise (Imax = 3,000). As shown in Figure 2(c), although the image restored by the PID-Split algorithm suppressed noise well, it also caused staircase effects. The images recovered by the PIDSB-NLFA and PIDSB-FA algorithms are shown in Figures 2(d) and 2(e). Compared with Figure 2(c), the staircase effects were suppressed, and the noise in the flat regions was suppressed better. Figure 2(f) shows the result of restoration of the DeepRED, which retains more details but cannot remove Poisson noise well. The image restored by the proposed method in this paper is shown in Figure 2(g). Compared with the other four methods, it not only suppressed noise to the greatest extent but also had clearer details and the best visual perception.Table 1
PSNR (dB) values of the five methods.
ImageGaussian blurImaxDegraded imagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodLenaSize = 5,σ = 1.63,00024.8425.0625.2225.5725.8026.914,00024.8725.0725.2525.6225.8226.945,00024.9025.1125.2925.6825.8527.11Size = 15,σ = 1.63,00024.3625.0125.1225.2625.8126.334,00024.3825.0325.1625.2825.9226.475,00024.4025.0725.2025.2926.0926.57Size = 25,σ = 23,00022.9923.7024.4424.5524.5724.614,00023.0023.7224.4924.5724.5924.645,00023.0223.7524.5124.5924.6124.66CameramanSize = 5,σ = 1.63,00023.6325.3826.6826.7327.1227.464,00023.6625.5126.7326.7527.2028.055,00023.6725.6926.8526.9227.4028.14Size = 15,σ = 23,00021.7323.3924.0224.0924.2424.594,00021.7523.4324.1524.1724.3824.645,00021.7623.4824.2024.2324.4724.71Size = 25,σ = 1.63,00022.8324.3925.2925.3725.5126.294,00022.8424.4925.3225.4025.7426.305,00022.8824.6625.3525.4225.9226.50ButterflySize = 5,σ = 1.63,00022.8824.5025.2325.5225.5625.954,00022.9024.6325.2625.5525.6226.015,00022.9124.7825.2725.5825.6626.12Size = 15,σ = 22,00020.5721.9622.9923.2923.4723.983,00020.6022.0423.0323.4723.5924.054,00020.6322.2223.0723.5723.6924.13Size = 25,σ = 1.63,00021.8123.1723.9624.3324.4624.814,00021.8323.3323.9924.3624.5224.895,00021.8423.4824.0024.3924.5924.98ZebraSize = 5,σ = 1.63,00023.5824.8125.8926.5026.8727.174,00023.6124.8625.9626.5626.9627.245,00023.6424.9226.0226.6227.0527.30Size = 15,σ = 23,00021.6923.9525.0525.3226.3826.764,00021.7124.0325.1225.4326.5726.855,00021.7324.1125.2025.5426.6926.92Size = 25,σ = 1.63,00022.6824.4125.5126.0426.7827.224,00022.7124.4725.5726.1726.9027.355,00022.7324.5425.6526.2927.0327.42HouseSize = 5,σ = 1.63,00028.8130.0830.7130.4332.4432.714,00028.9330.1730.9130.4832.5432.885,00028.9930.2831.1030.6032.6133.18Size = 15,σ = 23,00026.1127.1227.6927.5228.3628.474,00026.1427.2027.7027.5728.3728.545,00026.1827.3327.7327.6228.3928.62Size = 25,σ = 1.63,00027.7128.7029.6129.4130.2930.874,00027.8428.7829.7329.4930.4731.085,00027.8528.8729.8129.5630.5631.15Table 2
VIF values of the five methods.
ImageGaussian blurImaxDegraded imagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodLenaSize = 5,σ = 1.63,0000.44820.43100.45130.51470.51820.57794,0000.45530.43270.45590.51220.51310.59255,0000.46030.43690.46140.51470.51360.5984Size = 15,σ = 1.63,0000.40260.43070.42980.49380.51780.55854,0000.40740.43380.43460.49420.52820.56835,0000.41160.43930.44040.49460.54100.5745Size = 25,σ = 23,0000.34840.39460.47610.49490.51300.52904,0000.35310.39730.48120.49620.51160.52975,0000.35660.40030.48420.49690.52900.5373CameramanSize = 5,σ = 1.63,0000.30860.40450.43320.43360.45130.45334,0000.31070.41330.44360.43950.45630.46755,0000.31400.42330.44500.44580.45730.4723Size = 15,σ = 23,0000.22550.31320.32600.38600.37430.39814,0000.22800.31660.33470.38670.38900.40175,0000.23010.31950.33720.38610.39500.4105Size = 25,σ = 1.63,0000.27050.37300.38540.38600.39960.42954,0000.27330.38050.38940.38760.40400.44075,0000.27600.39100.39030.38610.41850.4447ButterflySize = 5,σ = 1.63,0000.39400.47090.55350.55970.57680.58874,0000.39670.48150.55770.56140.59090.60285,0000.39820.49330.56020.56350.60740.6102Size = 15,σ = 22,0000.29250.38180.46220.47170.50180.51243,0000.29570.38770.46610.48490.52280.52844,0000.29910.40020.47130.49670.53350.5390Size = 25,σ = 1.63,0000.34710.41190.51170.52920.54780.56134,0000.34920.41890.51420.53250.54880.56305,0000.35070.42630.51630.53580.56600.5800ZebraSize = 5,σ = 1.63,0000.39130.40070.49050.51920.49130.52064,0000.39790.40380.49570.52780.51060.53025,0000.40200.4070.049920.51090.51790.5355Size = 15,σ = 23,0000.28320.33480.39820.40940.42950.44904,0000.28780.33870.40260.41550.43940.45905,0000.29110.34290.40710.42190.44480.4622Size = 25,σ = 1.63,0000.34020.36970.40660.47120.47630.49334,0000.34630.37310.41000.47890.48500.50165,0000.35040.37680.41360.48570.49590.5082HouseSize = 5,σ = 1.63,0000.38420.40280.44720.50400.50620.51564,0000.38980.40810.45630.50490.51030.51725,0000.39240.41270.46590.50470.51330.5320Size = 15,σ = 23,0000.29140.36200.4378043570.43340.45594,0000.29290.36880.43800.43730.43700.46805,0000.29460.37530.44040.43750.43890.4749Size = 25,σ = 1.63,0000.34290.37430.46530.47190.45220.49464,0000.34960.37830.47410.47240.46860.50375,0000.34950.38230.47520.46630.47030.5091Figure 2
Restoration of the House image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 15 × 15,σ = 2 and corrupted by Poisson noise with Imax = 3,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).Figure3 shows a set of comparative results for Zebra images in Figure 1(d). The size of Zebra image is 584 × 387. Figure 3(b) is an image degraded by a Gaussian blur kernel with a size of 25 × 25 and a standard deviation of 1.6 and Poisson noise with a size of 3,000. As can be seen from Figure 3(a), the zebra has a large size of flat regions and the grassland contains considerable detailed information. As seen in Figure 3(g), the proposed method in this paper not only suppressed noise but also preserved considerable detailed information. Figure 3(c) shows that the image restored by the PID-Split algorithm was too smooth and had staircase effects. Figures 3(d)–3(f) show the images recovered by the PIDSB-NLFA, PIDSB-FA, and DeepRED algorithms, respectively. Compared with Figure 3(c), they suppressed the staircase effects while retaining some details. Compared with the other four method, the proposed method in this paper not only suppressed the staircase effects but also preserved the details of the image to the greatest extent and was clearer.Figure 3
Restoration of the Zebra image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 25 × 25,σ = 1.6 and corrupted by Poisson noise with Imax = 3,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).In addition, the Butterfly image in Figure1(c) was degraded by the Gaussian blur kernel (size = 5 × 5, σ = 1.6) and then corrupted by Poisson noise (Imax = 5,000). As shown in Figure 4, compared with the images restored by the other four algorithms, the proposed method in this paper suppressed noise and restored clear images to the greatest extent. Figure 4(g) restored by our method looks natural and preserved the most detailed information of the images.Figure 4
Restoration of the Butterfly image. (a) Ground truth image. (b) Degraded image blurred by Gaussian kernel of size 5 × 5,σ = 1.6 and corrupted by Poisson noise with Imax = 5,000. (c) Restored image by the PID-Split method. (d) Restored image by the PIDSB-FA method. (e) Restored image by the PIDSB-NLFA method. (f) Restored image by the DeepRED method. (g) Restored image by the proposed method. (h) Close-up views (close-up views of (a)–(g) image regions extracted from this example, respectively).


## 3.2. Convergence Analysis

The explicit regularization prior term adopts the combination of the intensity-basedL0 regularization term and RED. The purpose is to be able to use the pixel intensity distribution of the blurred image and the clear image to distinguish the restored image pixels while effectively filtering the image pixels and removing the Poisson noise of the degraded image. To verify the effectiveness of the L0 regularization prior, we used the method with the L0 norm regularization term and the method without the L0 norm to carry out simulation and comparison experiments on the House images of Figure 1(c) in Section 3.1. The experimental results are shown in Figure 5. Figure 5(a) shows the PSNR values of the image restored by the two algorithms. It can be clearly seen that the method with the L0 regularization term has the characteristic of fast convergence. With the increase in iteration times, the restored image reached a higher PSNR. Figure 5(b) shows the enlarged graph of the abscissa range (0∼500) in Figure 5(a). Figure 5(c) shows the loss of the objective function, and Figure 5(d) shows an enlarged graph with abscissa (0∼200) in Figure 5(c). It can be seen that the method with the L0 regularization term also achieved a smaller loss value.Figure 5
PSNR values and loss comparison. (a) PSNR. (b) PSNE (zoomed in). (c) Loss. (d) Loss (zoomed in).
(a)(b)(c)(d)


## 3.3. Denoising and Deblurring of Real Poisson Images

This section is an experimental comparison of denoising and deblurring of real Poisson images. Tests were carried out on four real astronomically degraded images. To evaluate the real Poisson image quality, we used the Blind/Referenceless Image Spatial Quality Evaluator (BRISQUE) to evaluate the real data. The smaller the value of BRISQUE, the better the image quality. The image quality evaluation data recovered by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms and the proposed method in this paper are shown in Table3.Table 3
BRISQUE values of the five methods.
ImagePID-SplitPIDSB-FAPIDSB-NLFADeepREDOur methodSaturn45.1343.3148.9641.4239.22Moon49.8143.3447.7539.2838.56Docking48.6549.1449.0541.9539.85Phoebe36.6132.9431.2630.4424.25Figure6(a) shows a satellite image, and Figure 6(b) shows an image restored by the PID-Split algorithm; the restored image is too smooth. Figures 6(c)–6(e) show the images recovered by the PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms, respectively. Although the staircase effects were suppressed to a large extent, ring effects exist at the edge. The image restored by the proposed method in this paper is shown in Figure 6(f). The restored image not only restored a large amount of detailed information and texture information but also suppressed the ring effects at the edge while suppressing the step effect. The image restored by the proposed method looks more natural.Figure 6
Restoration of real blurred image. (a) Saturn truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).The second real image is a moon image, as shown in Figure7(a). Figures 7(b)–7(e) show the images restored by the PID-Split, PIDSB-FA, PIDSB-NLFA, and DeepRED algorithms, respectively. It can be seen from these three images that the images recovered by these methods had ringing effects of the edge, and the flat areas were too smooth. Figure 7(f) shows the image restored by the proposed method in this paper. It can be seen that the image restored by the proposed method in this paper retained more detailed information while suppressing noise and ringing effects.Figure 7
Restoration of real blurred image. (a) Moon truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).Figure8 shows a set of comparative results for the docking images of Shenzhou-9 and Tiangong-1. Restored images by RLTV, PID-Split, PIDSB-FA, DeepRED, and the proposed method are shown in Figures 8(b)–8(f), respectively. The image restored by the proposed method suppressed the ring effects at the edge and retained more detailed information. Figure 9(a) shows the Phoebe image. Comparing the restored images by the other four methods shown in Figures 9(b)–9(e), respectively, the image restored by the proposed method in this paper looks clear with fine details.Figure 8
Restoration of real blurred image. (a) The docking truth image of Shenzhou-9 and Tiangong-1. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).Figure 9
Restoration of real blurred image. (a) Phoebe truth image. (b) Restored image by the PID-Split method. (c) Restored image by the PIDSB-FA method. (d) Restored image by the PIDSB-NLFA method. (e) Restored image by the DeepRED method. (f) Restored image by the proposed method. (g) Close-up views (close-up views of (a)–(f) image regions extracted from this example, respectively).


## 4. Conclusions

This paper proposed an algorithm for denoising and deblurring Poisson images by using neural networks. We combine an implicit regularization prior with two explicit regularization priors that are the prior of RED and the prior ofL0 regularization based on intensity and use a neural network to restore degraded images from Poisson images. In the proposed method, it has been proved that Poisson images restored by deep networks can effectively suppress Poisson noise. Meanwhile, the method with the L0 regularization term has the characteristic of fast convergence.



---

*Source: 9483521-2020-08-24.xml*