# Image Denoising Based on Improved Wavelet Threshold Function for Wireless Camera Networks and Transmissions

Authors: Xiaoyu Wang; Xiaoxu Ou; Bo-Wei Chen; Mucheol Kim
Journal: International Journal of Distributed Sensor Networks
 (2015)

---

## Abstract

A new wavelet threshold denoising function and an improved threshold are proposed. It not only retains the advantages of hard and soft denoising functions but also overcomes the disadvantages of the continuity problem of hard denoising function and the constant deviation of the soft denoising function in the new method. In the case of the improved threshold conditions, the new threshold function has a better performance in outstanding image details. It can adapt to different images by joining an adjusting factor. Simulation results show that the new threshold function has a better ability of performing image details and a higher peak signal to noise ratio (PSNR).


---

## Body

## 1. Introduction

In the process of image transmission, the image is likely to be affected by noise in all aspects of transmission [1]. Due to these effects, the image will be interfered by noise. The visual effects of image are poor because of low image quality. At this point, the main noise in image transmission is the Gaussian noise. Fingerprint recognition, face recognition, and medical image are in great need of the image details, but a great interference to the image is brought by the image noise. These interferences are also promoting the development of the image denoising technology. Accuracy of minutiae extraction is very important in the automatic fingerprint identification technology based on matching minutiae [2]. Therefore, it is particularly important to preserve image details while removing the noise [3, 4]. The geometric method is the earliest and most traditional one for face detection and recognition. It usually needs to detect facial features, the relative position between the important facial features, and the relevant parameters of these features [5]; this method is badly in need of the detail information of the image. Image noise makes the image quality change so that it is difficult to extract the image feature. Image denoising can not only improve the practicability of the fingerprint, face recognition, and medical imaging but also enhance the visual effect of the image. The algorithm in this paper mainly aims at Gaussian noise.Image denoising algorithms can be divided into spatial domain algorithm and frequency domain algorithm [6–9]. Spatial domain algorithm is to process the image pixels directly. In the frequency domain algorithm, the image is transformed from time domain to frequency domain. Then the image is processed in frequency domain and transformed to that in time domain. Mean filtering, median filtering, and adaptive Wiener filtering are classical spatial domain algorithms. It also caused some problems while removing the noise in these methods, such as serious blur of the image, especially the edge and details of the image [10]. Fourier transform is in common use in frequency domain; however, the domain characteristics of it disappeared after Fourier transform. It can neither detect the time position and the degree of intensity after signal transformation nor describe the local properties of time domain of the image [11]. The emergence and improvement of wavelet transform theory solves this problem. So wavelet has been widely used in image processing [12, 13].Wavelet denoising algorithms can basically be divided into three categories: (1) modulus maxima algorithm based on wavelet transform [14]; (2) correlation algorithm [15]; (3) wavelet threshold algorithm [16, 17], including the hard threshold and soft threshold.In the process of image denoising, the selection of threshold and the construction of threshold function are the keys to denoising. A new threshold and a new construction of threshold function are proposed in this paper. The new threshold and the new threshold function are used to obtain a higher value of the peak signal of noise ratio while more image details are retained.


## 2. Wavelet Threshold Principle

### 2.1. Wavelet Transforms in Two Dimensions

An image can be thought of as a two-dimensional function; then the wavelet transform of an image is a two-dimensional wavelet transform. A two-dimensional scaling functionφ(x,y) and three two-dimensional wavelets ψH(x,y),  ψV(x,y), and ψD(x,y) are needed in two-dimensional wavelet transform. Each two-dimensional function above is the result of multiplying one-dimension function φ and ψ. Excluding the product of one-dimensional results, the remaining products produced a separable scaling function
(1)φ(x,y)=φ(x)φ(y)
and separable “directionally sensitive” wavelets
(2)ψHx,y=ψxφy,ψVx,y=φxψy,ψDx,y=ψxψy.
These wavelets measure variations of the function of an image, the intensity variations for images along the gray level of different directions. ψH measures variations along columns (e.g., horizontal edges), ψV responds to variations along rows (like vertical edges), and ψD corresponds to variations along diagonals [18].Given separable two-dimensional scaling and wavelet functions, then the two-dimensional wavelet transform becomes very simple. Firstly, we give a scaling function and a translated basis function:(3)φj,m,nx,y=2j/2φ2jx-m,2jy-n,ψj,m,nix,y=2j/2ψi2jx-m,2jy-n,i=H,V,D.
Here, i is the directional wavelets in above “directionally sensitive” wavelets (H represents the horizontal direction, V represents the vertical direction, and D represents the diagonal direction). j controls the width of φj,m,m(x,y), and the term 2j/2 controls the amplitude of φj,m,m(x,y).  m and n determine the position of φj,m,m(x,y) along x-axis and y-axis. The discrete wavelet transform for two-dimensional image f(x,y) of size M×N is
(4)Wφj0,m,n=1MN∑x=0M-1∑y=0N-1fx,yφj0,m,nx,y,Wψij,m,n=1MN∑x=0M-1∑y=0N-1f(x,y)ψj,m,ni(x,y)hhhhhhhhhhhhhhhhhhhi=H,V,D.j0 is the arbitrary initial scale of the scaling function and the coefficients of Wψij,m,n can be considered the approximation of f(x,y) at scale j0. The coefficients of Wψij,m,n add details of horizontal, vertical, and diagonal for j≥j0. In the Wψij,m,n function, we usually let j0=0 and select N=M=2J. Then we can know j=0,1,2,…,J-1 and m,n=0,1,2,…,2j-1. f(x,y) can be obtained via the inverse discrete wavelet transform while the Wφ and Wψi are given. Consider
(5)f(x,y)=1MN∑m∑nWφ(j0,m,n)φj0,m,n(x,y)+1MN∑i=H,V,D∑j=j0∞∑m∑nWψij,m,nψj,m,nix,y.
We can use Figure 1 [18] to represent two-dimensional wavelet transform and use Figure 2 [18] to represent its inverse transform.Figure 1
The discrete wavelet transform. Here,hψ(-n), hψ(-n), hψ(-m), and hφ(-m) are wavelet vectors. 2↓ represents downsampling by 2.Figure 2
The inverse discrete wavelet transform. Here,hψ(m), hφ(m), hψ(n), and hφ(n) are wavelet vectors. 2↑ signifies upsampling by 2.Four coefficient matrices can be obtained by the first wavelet transform of an image; one is the coefficient matrix of low frequency by the scaling function expansion and the other three coefficient matrices are generated by the translated basis functions. Wavelet threshold denoising is to modify these four coefficient matrices with a threshold and then use the four coefficient matrices after modification to reconstruct the image.


### 2.2. Wavelet Threshold Principle

Donoho and Johnstone proposed the theory of wavelet threshold shrinkage in 1992; the main theoretical basis of this theory is as follows [16, 19, 20]. The energy distribution of the signal in Besov space is mainly concentrated in the limited coefficients in wavelet domain, whereas the noise’s energy is distributed throughout the wavelet domain. Therefore, the transform coefficients of signal are greater than that of noise after wavelet decomposition. We can find a suitable threshold. When the amplitudes of wavelet coefficients are less than the preset threshold, it is considered that the amplitude is caused by the noise and the amplitude is set to zero; when the amplitudes of wavelet coefficients are greater than the threshold, it is considered that the amplitude is caused by signal and the amplitude is preserved or revised. The image is reconstructed by the processed wavelet coefficients in the steps above to achieve the purpose of image denoising. This algorithm is simple, has better visual effects, and can keep the edge information.Assume that the original image issij; the model of the noise image is
(6)fij=sij+nij,
where fij, sij, and nij represent the noise image, the original image, and the corrupting white Gaussian noise, respectively. i,j=1,2,3…N, N is an integer power of 2.The steps of wavelet threshold denoising are as follows from the analysis above.(1)
Select a suitable wavelet base and decomposition levels to decompose the imagefij, and the wavelet coefficient of the noise image is wij. Here, the input is fij and the output is wij.(2)
Select an appropriate threshold to process the high frequent wavelet coefficients of the noise image in each layer, and obtain the new coefficientsw^ij. Here, the input is wij and the output is w^ij.(3)
Reconstruct the image with the original low frequency wavelet coefficients andw^ij to get the denoising image. Here, the inputs are the original low frequency wavelet coefficients and w^ij. And the final output is the denoising image.


## 2.1. Wavelet Transforms in Two Dimensions

An image can be thought of as a two-dimensional function; then the wavelet transform of an image is a two-dimensional wavelet transform. A two-dimensional scaling functionφ(x,y) and three two-dimensional wavelets ψH(x,y),  ψV(x,y), and ψD(x,y) are needed in two-dimensional wavelet transform. Each two-dimensional function above is the result of multiplying one-dimension function φ and ψ. Excluding the product of one-dimensional results, the remaining products produced a separable scaling function
(1)φ(x,y)=φ(x)φ(y)
and separable “directionally sensitive” wavelets
(2)ψHx,y=ψxφy,ψVx,y=φxψy,ψDx,y=ψxψy.
These wavelets measure variations of the function of an image, the intensity variations for images along the gray level of different directions. ψH measures variations along columns (e.g., horizontal edges), ψV responds to variations along rows (like vertical edges), and ψD corresponds to variations along diagonals [18].Given separable two-dimensional scaling and wavelet functions, then the two-dimensional wavelet transform becomes very simple. Firstly, we give a scaling function and a translated basis function:(3)φj,m,nx,y=2j/2φ2jx-m,2jy-n,ψj,m,nix,y=2j/2ψi2jx-m,2jy-n,i=H,V,D.
Here, i is the directional wavelets in above “directionally sensitive” wavelets (H represents the horizontal direction, V represents the vertical direction, and D represents the diagonal direction). j controls the width of φj,m,m(x,y), and the term 2j/2 controls the amplitude of φj,m,m(x,y).  m and n determine the position of φj,m,m(x,y) along x-axis and y-axis. The discrete wavelet transform for two-dimensional image f(x,y) of size M×N is
(4)Wφj0,m,n=1MN∑x=0M-1∑y=0N-1fx,yφj0,m,nx,y,Wψij,m,n=1MN∑x=0M-1∑y=0N-1f(x,y)ψj,m,ni(x,y)hhhhhhhhhhhhhhhhhhhi=H,V,D.j0 is the arbitrary initial scale of the scaling function and the coefficients of Wψij,m,n can be considered the approximation of f(x,y) at scale j0. The coefficients of Wψij,m,n add details of horizontal, vertical, and diagonal for j≥j0. In the Wψij,m,n function, we usually let j0=0 and select N=M=2J. Then we can know j=0,1,2,…,J-1 and m,n=0,1,2,…,2j-1. f(x,y) can be obtained via the inverse discrete wavelet transform while the Wφ and Wψi are given. Consider
(5)f(x,y)=1MN∑m∑nWφ(j0,m,n)φj0,m,n(x,y)+1MN∑i=H,V,D∑j=j0∞∑m∑nWψij,m,nψj,m,nix,y.
We can use Figure 1 [18] to represent two-dimensional wavelet transform and use Figure 2 [18] to represent its inverse transform.Figure 1
The discrete wavelet transform. Here,hψ(-n), hψ(-n), hψ(-m), and hφ(-m) are wavelet vectors. 2↓ represents downsampling by 2.Figure 2
The inverse discrete wavelet transform. Here,hψ(m), hφ(m), hψ(n), and hφ(n) are wavelet vectors. 2↑ signifies upsampling by 2.Four coefficient matrices can be obtained by the first wavelet transform of an image; one is the coefficient matrix of low frequency by the scaling function expansion and the other three coefficient matrices are generated by the translated basis functions. Wavelet threshold denoising is to modify these four coefficient matrices with a threshold and then use the four coefficient matrices after modification to reconstruct the image.


## 2.2. Wavelet Threshold Principle

Donoho and Johnstone proposed the theory of wavelet threshold shrinkage in 1992; the main theoretical basis of this theory is as follows [16, 19, 20]. The energy distribution of the signal in Besov space is mainly concentrated in the limited coefficients in wavelet domain, whereas the noise’s energy is distributed throughout the wavelet domain. Therefore, the transform coefficients of signal are greater than that of noise after wavelet decomposition. We can find a suitable threshold. When the amplitudes of wavelet coefficients are less than the preset threshold, it is considered that the amplitude is caused by the noise and the amplitude is set to zero; when the amplitudes of wavelet coefficients are greater than the threshold, it is considered that the amplitude is caused by signal and the amplitude is preserved or revised. The image is reconstructed by the processed wavelet coefficients in the steps above to achieve the purpose of image denoising. This algorithm is simple, has better visual effects, and can keep the edge information.Assume that the original image issij; the model of the noise image is
(6)fij=sij+nij,
where fij, sij, and nij represent the noise image, the original image, and the corrupting white Gaussian noise, respectively. i,j=1,2,3…N, N is an integer power of 2.The steps of wavelet threshold denoising are as follows from the analysis above.(1)
Select a suitable wavelet base and decomposition levels to decompose the imagefij, and the wavelet coefficient of the noise image is wij. Here, the input is fij and the output is wij.(2)
Select an appropriate threshold to process the high frequent wavelet coefficients of the noise image in each layer, and obtain the new coefficientsw^ij. Here, the input is wij and the output is w^ij.(3)
Reconstruct the image with the original low frequency wavelet coefficients andw^ij to get the denoising image. Here, the inputs are the original low frequency wavelet coefficients and w^ij. And the final output is the denoising image.


## 3. Wavelet Threshold Denoising Function

### 3.1. Several Common Denoising Threshold Functions

(
1
) The hard threshold denoising function.
The function is(7)
w
^
i
j
=
w
i
j
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the hard threshold function.
(
2
) The soft threshold denoising function.
The function is(8)
w
^
i
j
=
sgn
⁡
w
i
j
w
i
j
-
λ
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the soft threshold function. sgn⁡(wij) represents the sign of wij.
(
3
) The improved soft threshold denoising function.
The function is(9)
w
^
i
j
=
sgn
⁡
w
i
j
w
i
j
-
α
λ
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the improved soft threshold denoising function. sgn⁡(wij) represents the sign of wij and α(0<α<1) is an adjustment factor of λ.
(
4
) The semisoft threshold denoising function.
The function is(10)w^ij=wij,wij>λ2sgn⁡wijλ2wij-λ1λ2-λ1,λ1≤wij≤λ20,wij<λ1,
where wij, sgn⁡(wij) represent the wavelet coefficients of the noisy image and the sign of wij, respectively, and w^ij is the result of the wavelet coefficients processed by the soft threshold function. λ1,λ2 are the thresholds of this function. The images of all the functions above are shown in Figure 3.Figure 3
The images of threshold denoising functions above.
(a)
The hard threshold denoising function
(b)
The soft threshold denoising function
(c)
The improved soft threshold denoising function
(d)
The semisoft threshold denoising function
(e)
The improved wavelet threshold functionThe soft threshold function and hard threshold function of Donoho have been widely used, but it remains inadequate. The edges of the image are blurred because the wavelet coefficients of the noisy image that are greater than the threshold are shrunk in the soft threshold denoising function [21]. The hard threshold denoising function can keep the local information of the image. However, the image will show the visual distortion caused by pseudo Gibbs phenomena because the hard threshold denoising function is not continuous at the threshold [22]. In addition, highlighting the details of the image is not particularly obvious in the hard threshold and soft threshold denoising function. A new threshold function will be introduced in the following paper.


### 3.2. The Improved Wavelet Threshold Function

The new function is(11)w^ij=sgn⁡wijwij-sinπ2λwijnλ,wij>λ0,wij>λ,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the new threshold function. sgn⁡(wij) represents the sign of wij. And n is an adjustment parameter.In the new threshold function, ifwij→λ, the value of w^ij→0. We can know that the new function is continuous at the threshold point so that the pseudo Gibbs phenomena can be avoided. If wij→∞, the value of w^ij→wij. It can be seen that when wij is larger, the difference between wij and w^ij becomes smaller gradually. It can solve the problem of fixed deviation brought by the soft threshold function, and the edge blur of the image can be degraded.  n can be adjusted depending on the different thresholds of the different images to make the new function adapt to more images. The image of the improved function is shown in Figure 3.


## 3.1. Several Common Denoising Threshold Functions

(
1
) The hard threshold denoising function.
The function is(7)
w
^
i
j
=
w
i
j
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the hard threshold function.
(
2
) The soft threshold denoising function.
The function is(8)
w
^
i
j
=
sgn
⁡
w
i
j
w
i
j
-
λ
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the soft threshold function. sgn⁡(wij) represents the sign of wij.
(
3
) The improved soft threshold denoising function.
The function is(9)
w
^
i
j
=
sgn
⁡
w
i
j
w
i
j
-
α
λ
,
w
i
j
>
λ
0
,
w
i
j
≤
λ
,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the improved soft threshold denoising function. sgn⁡(wij) represents the sign of wij and α(0<α<1) is an adjustment factor of λ.
(
4
) The semisoft threshold denoising function.
The function is(10)w^ij=wij,wij>λ2sgn⁡wijλ2wij-λ1λ2-λ1,λ1≤wij≤λ20,wij<λ1,
where wij, sgn⁡(wij) represent the wavelet coefficients of the noisy image and the sign of wij, respectively, and w^ij is the result of the wavelet coefficients processed by the soft threshold function. λ1,λ2 are the thresholds of this function. The images of all the functions above are shown in Figure 3.Figure 3
The images of threshold denoising functions above.
(a)
The hard threshold denoising function
(b)
The soft threshold denoising function
(c)
The improved soft threshold denoising function
(d)
The semisoft threshold denoising function
(e)
The improved wavelet threshold functionThe soft threshold function and hard threshold function of Donoho have been widely used, but it remains inadequate. The edges of the image are blurred because the wavelet coefficients of the noisy image that are greater than the threshold are shrunk in the soft threshold denoising function [21]. The hard threshold denoising function can keep the local information of the image. However, the image will show the visual distortion caused by pseudo Gibbs phenomena because the hard threshold denoising function is not continuous at the threshold [22]. In addition, highlighting the details of the image is not particularly obvious in the hard threshold and soft threshold denoising function. A new threshold function will be introduced in the following paper.


## 3.2. The Improved Wavelet Threshold Function

The new function is(11)w^ij=sgn⁡wijwij-sinπ2λwijnλ,wij>λ0,wij>λ,
where wij, λ represent the wavelet coefficients of the noisy image and the threshold, respectively, and w^ij is the result of the wavelet coefficients processed by the new threshold function. sgn⁡(wij) represents the sign of wij. And n is an adjustment parameter.In the new threshold function, ifwij→λ, the value of w^ij→0. We can know that the new function is continuous at the threshold point so that the pseudo Gibbs phenomena can be avoided. If wij→∞, the value of w^ij→wij. It can be seen that when wij is larger, the difference between wij and w^ij becomes smaller gradually. It can solve the problem of fixed deviation brought by the soft threshold function, and the edge blur of the image can be degraded.  n can be adjusted depending on the different thresholds of the different images to make the new function adapt to more images. The image of the improved function is shown in Figure 3.


## 4. Matlab Simulation Experiment and Results

### 4.1. Simulation Experiment

In order to show the superiority of the proposed algorithm, soft threshold and hard threshold denoising were selected as contrast. In this section, simulated experiments are performed with “kids.tif,” “mir.tif” (a medical image), and “105_7.tif” (an image in FVC2004 fingerprint database). This experiment has been performed using Matlab programming language (Matlab R2012b). The specific experimental steps are as follows.(1)
Add Gaussian white noise that the mean is 0 and the variance is 0.01 into the selected images.(2)
Decompose all the noisy images with wavelet base “sym4” to obtain the wavelet coefficients of each layer of all the images.(3)
Denoise the noisy images with the soft and hard threshold denoising method.(4)
Denoise the noisy images with the improved denoising method in this paper.(5)
Reconstruct image with the wavelet coefficients processed by soft, hard, and the new threshold denoising method, and compare the performance of the three methods.(6)
Use the peak signal to noise ratio (PSNR) as the evaluation criteria to evaluate each algorithm. Here, PSNR is defined as(12)PSNR=10lgfmax⁡21/MN∑i=0M-1∑j=0N-1fi,j-f0i,j2,where fmax⁡ is the maximum gray level of the image, M and N denote the size of the image, f(i,j) is the restored image, and f0(i,j) is the original image.


### 4.2. Experimental Results

The original image, noisy image, and the image after denoising processing are shown in Figure4.Figure 4
Images denoising. (a) Original image. (b) Image corrupted by Gaussian white noise that the mean is 0 and the variance is 0.01. (c) The denoising result of hard threshold. (d) The denoising result of soft threshold. (e) The denoising result of improved algorithm.
(a)
(b)
(c)
(d)
(e)Local enlarged images of denoising images are shown in Figure5.Figure 5
Local enlarged denoising images. (a) The hard threshold denoising algorithm. (b) The soft threshold denoising algorithm. (c) The improved wavelet threshold algorithm.
(a)
(b)
(c)We setn=0.02 for the improved threshold denoising function in the simulation experiment. PSNR of the different variance for these 3 images are shown in Tables 1, 2, and 3.Table 1
Comparison of PSNR (dB): kids.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
31.7121
27.6174
26.3008
Result of soft threshold
33.7362
29.2526
27.3461
Result of improved algorithm
34.2701
29.7518
28.0052Table 2
Comparison of PSNR (dB): mir.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
28.3733
25.6540
25.2654
Result of soft threshold
28.8702
26.6270
26.2286
Result of improved algorithm
29.9496
26.7387
25.9008Table 3
Comparison of PSNR (dB): 105_7.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
35.4000
34.6260
34.1692
Result of soft threshold
35.5536
35.1307
35.1261
Result of improved algorithm
36.4919
35.5582
35.1889


### 4.3. Selection of Threshold

Selecting the threshold is critical in the wavelet threshold denoising method. The suitability of it affects the denoising effect of image directly. In the VisuShrink algorithm proposed by Donoho, Universal threshold expression [19] is λ=δ2ln⁡(N), where N represents the length of the signal; δ=medianwij/0.6745, wij∈HH1, and HH1 is the diagonal coefficients of the first layer of wavelet decomposition for the noisy image.The threshold value used in this paper is a modification of [2]. The threshold expression in [2] is λ=δ2lb(m)/j, where δ, m, and j represent a standard deviation of the noisy image, the total number of wavelet coefficients of the image, and the level of decomposition. The noise amplitude is smaller and smaller while the signal amplitude is larger and larger with the increasing of decomposition levels of the noisy image. The threshold expression used in this paper is λj=δj2lb(m)/(j+2), where λj, m, and j represent the threshold of level j after wavelet decomposition of the noisy image, the total number of wavelet coefficients of the image,  and the level of decomposition. Here, σj=median(wij)/0.6745, where wij is the high frequency coefficients of the horizontal, vertical, and diagonal directions of layer j.


## 4.1. Simulation Experiment

In order to show the superiority of the proposed algorithm, soft threshold and hard threshold denoising were selected as contrast. In this section, simulated experiments are performed with “kids.tif,” “mir.tif” (a medical image), and “105_7.tif” (an image in FVC2004 fingerprint database). This experiment has been performed using Matlab programming language (Matlab R2012b). The specific experimental steps are as follows.(1)
Add Gaussian white noise that the mean is 0 and the variance is 0.01 into the selected images.(2)
Decompose all the noisy images with wavelet base “sym4” to obtain the wavelet coefficients of each layer of all the images.(3)
Denoise the noisy images with the soft and hard threshold denoising method.(4)
Denoise the noisy images with the improved denoising method in this paper.(5)
Reconstruct image with the wavelet coefficients processed by soft, hard, and the new threshold denoising method, and compare the performance of the three methods.(6)
Use the peak signal to noise ratio (PSNR) as the evaluation criteria to evaluate each algorithm. Here, PSNR is defined as(12)PSNR=10lgfmax⁡21/MN∑i=0M-1∑j=0N-1fi,j-f0i,j2,where fmax⁡ is the maximum gray level of the image, M and N denote the size of the image, f(i,j) is the restored image, and f0(i,j) is the original image.


## 4.2. Experimental Results

The original image, noisy image, and the image after denoising processing are shown in Figure4.Figure 4
Images denoising. (a) Original image. (b) Image corrupted by Gaussian white noise that the mean is 0 and the variance is 0.01. (c) The denoising result of hard threshold. (d) The denoising result of soft threshold. (e) The denoising result of improved algorithm.
(a)
(b)
(c)
(d)
(e)Local enlarged images of denoising images are shown in Figure5.Figure 5
Local enlarged denoising images. (a) The hard threshold denoising algorithm. (b) The soft threshold denoising algorithm. (c) The improved wavelet threshold algorithm.
(a)
(b)
(c)We setn=0.02 for the improved threshold denoising function in the simulation experiment. PSNR of the different variance for these 3 images are shown in Tables 1, 2, and 3.Table 1
Comparison of PSNR (dB): kids.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
31.7121
27.6174
26.3008
Result of soft threshold
33.7362
29.2526
27.3461
Result of improved algorithm
34.2701
29.7518
28.0052Table 2
Comparison of PSNR (dB): mir.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
28.3733
25.6540
25.2654
Result of soft threshold
28.8702
26.6270
26.2286
Result of improved algorithm
29.9496
26.7387
25.9008Table 3
Comparison of PSNR (dB): 105_7.tif.
Variance of Gaussian noise
0.01
0.03
0.05
Result of hard threshold
35.4000
34.6260
34.1692
Result of soft threshold
35.5536
35.1307
35.1261
Result of improved algorithm
36.4919
35.5582
35.1889


## 4.3. Selection of Threshold

Selecting the threshold is critical in the wavelet threshold denoising method. The suitability of it affects the denoising effect of image directly. In the VisuShrink algorithm proposed by Donoho, Universal threshold expression [19] is λ=δ2ln⁡(N), where N represents the length of the signal; δ=medianwij/0.6745, wij∈HH1, and HH1 is the diagonal coefficients of the first layer of wavelet decomposition for the noisy image.The threshold value used in this paper is a modification of [2]. The threshold expression in [2] is λ=δ2lb(m)/j, where δ, m, and j represent a standard deviation of the noisy image, the total number of wavelet coefficients of the image, and the level of decomposition. The noise amplitude is smaller and smaller while the signal amplitude is larger and larger with the increasing of decomposition levels of the noisy image. The threshold expression used in this paper is λj=δj2lb(m)/(j+2), where λj, m, and j represent the threshold of level j after wavelet decomposition of the noisy image, the total number of wavelet coefficients of the image,  and the level of decomposition. Here, σj=median(wij)/0.6745, where wij is the high frequency coefficients of the horizontal, vertical, and diagonal directions of layer j.


## 5. Conclusion

Wavelet transform has a good time-frequency characteristics. Because the frequency distributions of noise and image signal are different, wavelet transform for image denoising has obvious advantages. In this paper, we proposed a new threshold function to improve inadequacies of the soft threshold and hard threshold based on the original threshold function by analyzing the characteristics of wavelet threshold function. Under different noise conditions, the new threshold function can retain more image details and get a higher PSNR. Improved threshold denoising function has a good application prospect in the face and fingerprint denoising and medical imaging denoising.



---

*Source: 670216-2015-09-29.xml*