An Adaptive Image Denoising Model Based on Tikhonov and TV Regularizations

To avoid the staircase artifacts, an adaptive image denoising model is proposed by the weighted combination of Tikhonov regularization and total variation regularization. In our model, Tikhonov regularization and total variation regularization can be adaptively selected based on the gradient information of the image. When the pixels belong to the smooth regions, Tikhonov regularization is adopted, which can eliminate the staircase artifacts. When the pixels locate at the edges, total variation regularization is selected, which can preserve the edges. We employ the split Bregman method to solve our model. Experimental results demonstrate that our model can obtain better performance than those of other models.


Introduction
Images are inevitably contaminated by noise during acquisition and transmission and so image denoising is one of the most fundamental tasks in image processing and computer vision [1].Image denoising aims at preserving edges and finescale structural information while removing the noise in the image.During the past two decades, a great deal of research [2][3][4] has gone into excogitating models for removing noise while preserving edges and other fine-scale image details.Recently, partial differential equation (PDE) and variationalbased methods have gained popularity for image denoising [5,6].Among them, the best known one is TV model proposed by Rudin et al. [7].Although this second-order PDE can succeed in preserving edges, it suffers from the staircase effect in the smooth regions.To prevent the staircasing, a modified model proposed by Blomgren et al. [8] was to automatically adapt the gradient-based exponent to fit the data term.It has the following form: min Gradient-based exponent has value 1 near the edges and is close to 2 in the smooth region, so that near the edges it behaves exactly like the ROF model and behaves like Tikhonov regularization method in the smooth regions.However, this model is nonconvex and difficult to solve.Another approach to overcome the limitation of TV is to introduce higher-order derivatives into the energy function [9,10].One of the most classical fourth-order PDEs is introduced by You and Kaveh (Y-K), which measured the oscillations in the noisy image according to the two different functions [11].But, a major challenge higher-order PDEs face is numerical computation.At the same time, another challenge is to blur the edges in the denoising.In addition, they can also produce speckle noise in the recovered image.Recently, some hybrid models were proposed which combined the second-order partial differential equation models and the fourth-order PDEs [12,13].Gholami and Hosseini [14] proposed a novel method to combine total variation regularization and Tikhonov regularization to reconstruct piecewise smooth signals.
Inspired by their work, we propose an adaptive image denoising model by the weighted combination of Tikhonov regularization and total variation regularization.Unlike the model [14] which firstly divided the image into piecewise constant component and smooth component, our model directly uses the gradient information to define the weight 2 Advances in Multimedia functions which judge whether the pixels locate at the edges or the smooth regions.When the pixels belong to the smooth regions, Tikhonov regularization is adopted, which can eliminate the staircase artifacts.When the pixels locate at the edges, total variation regularization is selected, which can preserve the edges.Our model does not solve the highorder PDE and is also a convex combination.All numerical experiments verify the efficiency of our proposed model.
The remainder of this paper is organized as follows.In Section 2, we briefly review Tikhonov regularization and TV regularization and then discuss their advantages and limitations.In Section 3, we give the proposed model and numerical implementation in detail.The experimental results are given in Section 4. Finally, this paper is concluded in Section 5.

Related Works and Discussion
In this paper, the observed noisy image  0 (Ω ⊂  2 → ) can be defined as  0 =  + ; here,  is the original image, and  is an additive Gaussian white noise with the mean 0 and variance  2 .Mathematically speaking, the denoising problem is an ill-posed inverse problem.One of the most known techniques to solve this problem is regularization method and energy minimization as follows: Here, the first term () is a regularization term, while the second term is fidelity term which ensures that the denoised image is close to the observed image, and  > 0 is the regularization parameter which balances the regularization term and the fidelity one.Much research has been put forward to remove the noises from the noisy image by either exploiting a particular regularization technique [15,16] or taking the advantages of a set of regularizations using a combination of them [14,17,18].In the image denoising, Tikhonov regularization and total variation regularization are two kinds of regularization techniques and they are widely used in the past three decades.
2.1.Tikhonov Regularization.Tikhonov regularization method is one of the earliest used regularizations, in which the cost function contains an  2 -norm regularization term of the magnitude of the image gradient [19].The model based on the Tikhonov regularization method is defined as follows: where ‖ ⋅ ‖ 2 is  2 norm of the bounded variation.We know that minimizing formula (3) yields the Euler-Lagrange equation Seen from (4), Tikhonov regularization technique is actually isotropic diffusion.Hence, while removing the noise in the observed image  0 , prominent structures like edges and jumps in the denoised image  are blurred.However, the first-and second-order derivative operators are only used in Tikhonov regularization, so it has computational advantages [20].
2.2.Total Variation Regularization.In addition, there are other regularization techniques which preserve the edge information in the recovered image [21].Among these edgepreserving regularization methods [7], the best known one is the total variation regularization.Rudin et al. proposed using the  1 norm of the gradient of the image, instead of the  2 norm, so the formulation for the variational model (ROF) is as follows: subject to Here Ω is a bounded domain of the image with Lipschitz boundary.BV(Ω) is the space of functions with bounded variation.|Ω| represents the area of the image.By the gradient descent method, we get the Euler-Lagrange equation of TV denoising model from ( 5) as follows: When |∇| = √|∇  | 2 + |∇  | 2 , formula ( 5) is the isotropic TV of .∇   and ∇   are the horizontal and the vertical gradient operators, respectively.The isotropic TV regularization performs well in preserving edges while eliminating noise.However, the fine-scale details such as textures are often smoothed out with noise in the image.Esedoǵlu and Osher [21] proposed anisotropic total variation (ATV) regularization as |∇| = |∇  | + |∇  |.The uniqueness and existence of solutions of ATV are shown in [22].However, the ATV is still TV-type regularization which is not good at preserving fine-scale details and may cause staircase artifacts in the smooth regions.these two regularizations and take the energy function of the image for image denoising as follows:

The Proposed Method
Here the functions () are defined as where  is constant and  is defined as the image gradient modulus value.The function () ∈ [0, 1] decreases monotonically from 1 to 0. When the pixels belong to edges, the value of  is large, so () → 1.When the pixels belong to flat regions,  is small, so () → 0. From the proposed model, we can see that when () → 0, which means the pixel locates in the edge, TV regularization term is adopted; when () → 1, which means the pixel belongs to the smooth region, Tikhonov regularization is used.So we can see from our model (8) the following.
(1) In the region containing edges, when () is close to 0, the new model adopts the total variation (TV) model, so the model is good at preserving the edges of the image.
(2) In the smooth areas of the image, () is close to 1, and the new model will highlight Tikhonov regularization model, so the model can eliminate the staircase effect.
(3) In the flat areas of image containing less image features and much noise, Tikhonov regularization model and total variation model work together.

Numerical Implementation.
To solve (8), we use split Bregman method proposed by Goldstein and Osher [23].Split Bregman is a flexible and effective tool for solving where () = (|∇|). is used to replace (|∇|).We can see that the first and second terms are not directly crossed.Second, by introducing the variable , the constrained minimization problem (10) can be converted to the following unconstrained problem: Advances in Multimedia where the penalty parameter  is a position constant.Using split Bregman method, the unconstrained problem ( 11) can be solved by the following iterative equations: To solve the first problem in ( 12), we let The partial derivatives of ( 13) are showed as follows: Therefore the Euler-Lagrange equation of ( 14) is According to the gradient descent method, we can derive We also derive the Euler-Lagrange equation of the second problem in (12) as follows: This is a linear equation of , so additional operator split (AOS) iteration is used to solve (17).
In sum, the new model can be implemented as follows.
Algorithm 1.The proposed algorithm for image denoising is as follows.
Step 1. Input the observed image  0 .
Until the stop condition is satisfied.

Experiments
To evaluate the performance of the proposed model, in this section, we have compared our model with other models in terms of the visual quality of denoised images, peak signal to noise ratio (PSNR), and structural similarity index measure (SSIM).PSNR is defined by with where |Ω| denotes the image size,  is the recovered image, and  0 is the noise-free original image.MSE represents the noise variance in the recovered image.PSNR is used to measure the noise removal effect and the larger value indicates the smoother recovered image.However, PSNR sometimes disaccords with human visual judgments.So, SSIM (structural similarity) is used to assess the noise removal quality because the SSIM criterion is closer to the human vision system [24].SSIM is defined by where   and  2  are the mean and variance of , respectively,   0 is the covariance of  and  0 , and  1 and  2 are two constants to avoid instability.SSIM measures the structure similarity between the original image and the recovered image.A stopping condition for all experiments is defined as follows: +1 − Here  is a given positive number.We set  = 10 − 3 in the paper.The other denoising parameters in the experiments are set as follows: the time step Δ = 0.2 and grid step size ℎ = 1,  = 0.02, and  = 1 to 25.We find that  has a deep influence on the denoising results.All experiments were implemented by MATLAB R2009a and performed on 64-bit Windows 7 on the desk with an Intel CPU of 1.7 GHz and 4 GB memory.is the noisy image with Gaussian noise whose standard deviation is 30 and mean value is zero.Figures 1(c)-1(h) show the denoising results of nonlocal means model [2], Y-K model [11], Perona-Malik model [4], total variation (TV) method [7], modified TV (MTV) [25], and our proposed model.Figure 2 shows the locally enlarged images of Figure 1.From Figures 1 and 2, we can observe that our model yields better results in the visual sense than other models in the visual sense, since it eliminates the staircase effect of TV model and preserves edges as well as TV model.In other words, our model can take advantage of both the filter based on Tikhonov regularization and the one based on total variation regularization.
Table 1 shows the PSNR, SSIM, and computational times of Figure 1.From Table 1, we also see that our model also outperforms other models in PSNR and SSIM.The experiments of our model are solved iteratively like TV, PM, and MTV.Seen from Table 1, the proposed method uses less computational time than Y-K and NLM, but the computational time of our proposed method is higher than those of TV, PM, and MTV.

Experiments with Natural Images.
To verify the performance of our model comprehensively, we have also implemented the quantitative and qualitative evaluation on the natural images from http://decsai.ugr.es/cvg/dbimagenes/.Figures 3, 4, and 5 exhibit the denoising results on the noisy images "lena, " "peppers, " and "house" with  = 30.The figures show that the denoising results of our model look better than other models.The results of TV model and modified TV cause the staircase effect.Blurring edges come out in Perona-Malik model.The result of Y-K has speckles.Figure 6 shows the partial enlarged images of Figure 3. Tables 2 and 3 show the PSNR and SSIM of our experiments, from which we can see that the proposed algorithm in all noise power has the highest PSNR and SSIM among the six methods.In addition, we show that the experimental results on the "lena" image with  = 10, 20, 30, 40, and 50 in Figure 7 and Tables 4 and 5 give the quantitative experimental results of PSNR and SSIM, respectively.These experimental results demonstrate that our model has better performance.

Comparison with BM3D.
Block-matching and 3D filtering (BM3D) [3] is the state-of-the-art denoising method, which enhances the sparsity by grouping similar 2D image blocks into 3D groups based on the fact that an image has a locally sparse representation in transform-domain.We compare the proposed model with BM3D whose code is downloaded from http://www.cs.tut.fi/∼foi/GCF-BM3D/.As shown in Figure 8 and Table 6, we find that the proposed model is not necessarily better than BM3D in removing the noises from the images.But, our model is still novel considering that it is a significant improvement on the other models based on PDE.

Conclusion
To eliminate the so-called staircase effect in many known image denoising models and speckles in the high-order PDE, we propose an improved model for image denoising with a weighted combination of Tikhonov regularization and total variation regularization.Numerical experiments show that our proposed model can recover the images from the noisy images, has the highest PSNR and SSIM among the six methods and can preserve important structures, such as edges and corners.

Figure 1 :Figure 2 :
Figure 1: Result of the synthetic image by the five models: (a) original image, (b) noisy image, (c) result by NLM, (d) result by Y-K model, (e) result by PM, (f) result by TV, (g) result by MTV, and (h) result by our model.

Figure 3 :Figure 4 :
Figure 3: Results of the lena image by the five models: (a) original image, (b) noisy image, (c) result by NLM, (d) result by Y-K model, (e) result by PM, (f) result by TV, (g) result by MTV, and (h) result by our model.

Figure 5 :Figure 6 : 2 + 𝜆 2 (
Figure 5: Results of the house image by the five models: (a) original image, (b) noisy image, (c) result by NLM, (d) result by Y-K model, (e) result by PM, (f) result by TV, (g) result by MTV, and (h) result by our model.

Figure 1 (
a) is the original image and Figure 1(b)

Figure 7 :
Figure 7: PSNR and SSIM values of six models for denoising lena image.

Figure 8 :
Figure 8: Results of our model and BM3D: (a) original image, (b) noisy image with the mean value of 0 and standard deviation of  = 30, (c) results by our model, and (d) results by BM3D.

Table 1 :
PSNR and SSIM comparisons with different methods for the synthetic image (256 × 256).

Table 4 :
PSNR comparison with different methods for the "lena" image.

Table 5 :
SSIM comparison with different methods for the "lena" image.

Table 6 :
PSNR and SSIM comparisons between our model and BM3D.