A Sobel-TV Based Hybrid Model for Robust Image Denoising

The traditional Total-Variation algorithm has a good result to de-noise for noise image of small scale details, but it easily losses the details for the image with rich texture and tiny boundary. In order to solve this problem, this paper proposes a Sobel-TV model algorithm for image denoising. It uses TV model to de-noise and uses Sobel algorithm to control smoothness of image, which not only efficiently removes image noise but also simultaneously retail information, such as edge and texture. The experiments demonstrate that the proposed algorithm is simple, practical and generates better SNR, which is an important value to preprocess image.


Introduction
The actual image is inevitably disturbed by various noises while it may be introduced by the image formation process, image recording, image transmission, etc.So image denoising is a key component in image processing.The problem aims to remove noises and preserve edges and small scale structures such as textures, which is to enhance signal to noise ratio and improve image quality.Nowadays researchers have proposed many methods to solve this problem of image denoising, such as median filtering methods, smooth filtering methods, total variation (TV) [1] [2] diffusion methods, etc. Filtering methods tend to blur an image but edges in the image are not well preserved.Though TV De-noising method, [3] [4] can remove image noise and preserve image edges characteristics.This method will retain salt and pepper noise point as the edge of the image, which cannot remove salt and pepper noise [5]- [7].In this paper, we propose a Sobel-TV model algorithm for image denoising.The algorithm distills the image edge through using Sobel algorithm, then uses TV algorithm to remove image noise according to image edge information.The algorithm not only suppresses image noise, but also maintains the image edge and texture information.

Sobel-TV Model
The TV (total variation) method makes the image restoration problem transformed into a functional extremum problem through introducing the energy function, which is widely used in image denoising.The image denoising may be mathematically modeled as where 0 u is the noise image, u is the original image, n is noise, J is the Gaussian convolution operator.Set J I = (Unit matrix), and assuming that the image noise is Gaussian noise that the mean value µ is a zero and the variance is 2 σ , hence (1) simplifying gives: The image denosiing TV model to recover nosied image 0 u to original image u is ( ) ( ) In ( 3), ( ) In (4), ( ) ∇ is the gradient of image u.The TV model is a piecewise continuous function in the bounded variation space, hence, the model has good removal effect for the repeating patterns of small-scale details and noise, but it is likely to cause loss of detail for he rich texture image.
In (4) above, fidelity approximation ( ) keeps the original image features and reduces image distortion effect, diffusion term div u u remove noise of the image.When diffusion term removes noise of the image, it also make the whole image smooth, which losses the small feature in the image and blurres image.To solve the problem above, in (4) introduce a guide function ( ) , m x y , and its improved form as follows: , m x y has the feature as follows: the value of ( ) , m x y is the smaller in the larger gray value of the place; the value of ( ) , m x y is the bigger in the slower gray value of the place.Therefore, the function ( ) , m x y can selectively denoise or smooth image according to the gray-scale value changing in the each part of the image.From above, the function ( ) , m x y is expressed as follows: In ( 6), ( ) , g x y is the edge detection value in the coordinate point ( ) T is a preset threshold and its value may be changed according the demand of image processing, ( ) Through edge detection, 0 u can be transformed into ( ) , g x y .The math model of edge detection have first order differential operators and second order differential operators.Because the second-order differential operator computation is very complexity, we choice first-order differential operators, which is more flexibility.Firstorder differential operators that use to detect image edge have two popular methods: prewitt operator and sobel operator.The prewitt operator is very sensitive to noise and remove image noise, but it maintains worse edges in the image than sobel operator.Sobel operator weights the impact of the image pixel location, which can reduce the degree of edge blur, so it has better effect.
For digital image ( ) { } , f i j , the discrete form of the Sobel operator is defined as follows: In (7), A is vertical direction edge detection.
In ( 8), B is horizontal direction edge detection.From ( 7) and ( 8) above, we have , S i j is the image after the Sobel operator detects edge.Figure 1 displays sobel operator template.

Algorithm Implementation
Sobel-TV algorithm for image denoising is firstly using sobel operator to detect image edges, then using TV to remove image noise according to the image edge information.The specific algorithm is as follows: Step 1. Choose Sobel operator to detect original image ( ) Step 2. Using (12), remove image noise according to the image edge information. ( , div , , , , where ( ) is the result of ( ) λ is the parameter of the regular.In order to avoid the expression f f ∇ ∇ is not zero, we have introduced the parameter β by defining According to divergence, it arrives as follows.
where f x the first derivatives in the x direction; f y the first derivatives in the y direction; f xx the second derivatives in the x direction f yy the second derivatives in the y direction; f xy the second derivatives for f x on y.The boundary conditions is: The scale parameter formula is as follows: In this paper, the parameter is set as follows:

Numerical Experiments
We performed experiments on several images including the well-known Lean, Cameraman, Peppers and Barbara to show the performance of our proposed model.All numerical experiments were performed on 64-bit Windows 7 on a desktop with an Intel CPU at 3.0 GHz and 4 GB memory.In order to estimate the performance of different methods, we use Signal to Noise Ratio (SNR) as the image quality measure that is defined as, ( , mean SNR 10 log , , where u 0 is the original image, u * is the mean intensity value of u, and u is the restored image.We also adopt the structure similarity index measure (SSIM) to measure the similarity between two images, which is defined as where u µ is the average of u 0 , 2 σ is the variance of u 0 , * 0 u u σ is the covariance of u * and u 0 and c 1 c 2 are two constants to avoid instability.

Results of Comparison
In this subsection, we tested the algorithm on four different images.The noisy images are shown in    testing results demonstrate that the Sobel-TV algorithm for de-noising image is better than the other algorithm.
In Table 1 and Table 2, we list the SNR and SSIM values for mean filter, median filter and Sobel-TV.The best SNR and SSIM value for Sobel-TV are obtained.From Figures 3-6 and the numerical results in the Table 1, Table 2 it is evident that the proposed method is effective in both removing image noise and maintaining retail information.

Conclusion
In this paper, we have proposed Sobel-TV models that use Sobel-TV to remove the noise for digital images.The proposed algorithm is firstly using sobel operator to detect image edges, then using TV to remove image noise according to the image edge information, which not only removes noise but also preserves the image contours and texture detail information.We have also presented some promising evidence for our methods on a real world dataset.The result of experiment shows that the proposed algorithm is validating the efficacy.

Figure 2 (
Figure 2(a) and Figure 2(b)show the relationship between SNR or SSIM and the number of iterations.The results show that the SNR and SSIM is the maximum value when the number of iterations is 100.Therefore, the number of iterations in (12) chooses 100.

Figure 2 .Figure 3 .Figure 4 .Figure 5 .
Figure 2. (a) The relationship between SNR and the number of iterations; (b) The relationship between SSIM and the number of iterations.

Figure 6 .
Figure 6.The contract of de-noising for three methods to image Peppers.(a) Original image; (b) Noisy image; (c) Denoising image by mean filter; (d) De-noising image by median filter; (e) De-noising image by Sobel-TV.
and have Lipschitz continuous boundary; ∇ is the mode of the image gradient; the first one of the right of equation, called regex, is TV norm of image u, which relies on the edge of the image.Regular in the TV model has played a very important role, it allows that the image has non-continuous part, but does not allows that the image has oscillations part, so it can remove the noise of the image.The second of the right of equation is approximation term, which controls the different of the image u and the observation image 0

Table 1 .
Comparison of the SNR values.

Table 2 .
Comparison of the SSIM values.