Research Article | Open Access

Min Wang, Zhen Li, Xiangjun Duan, Wei Li, "An Image Denoising Method with Enhancement of the Directional Features Based on Wavelet and SVD Transforms", *Mathematical Problems in Engineering*, vol. 2015, Article ID 469350, 9 pages, 2015. https://doi.org/10.1155/2015/469350

# An Image Denoising Method with Enhancement of the Directional Features Based on Wavelet and SVD Transforms

**Academic Editor:**Yann Favennec

#### Abstract

This paper proposes an image denoising method, using the wavelet transform and the singular value decomposition (SVD), with the enhancement of the directional features. First, use the single-level discrete 2D wavelet transform to decompose the noised image into the low-frequency image part and the high-frequency parts (the horizontal, vertical, and diagonal parts), with the edge extracted and retained to avoid edge loss. Then, use the SVD to filter the noise of the high-frequency parts with image rotations and the enhancement of the directional features: to filter the diagonal part, one needs first to rotate it 45 degrees and rotate it back after filtering. Finally, reconstruct the image from the low-frequency part and the filtered high-frequency parts by the inverse wavelet transform to get the final denoising image. Experiments show the effectiveness of this method, compared with relevant methods.

#### 1. Introduction

Denoising is important in image processing. There are many techniques for denoising, such as using the wavelet transform [1–8] and singular value decomposition (SVD) filters [9]. Wavelet transform has multiresolution properties, and the noise of the image is mainly in the high-frequency parts [1–3]; the wavelet threshold method is most widely used for noise filtering, which is effective, but it also removes some useful edge information.

Singular value decomposition (SVD) is generally used as a kind of nonlinear filter, with good numerical robustness [9]. SVD can handle different types of noise, without the prior knowledge about the noise [9]. The directional features of SVD can be used to filter the high-frequency parts of an image, if with proper handling.

This paper proposes a method for image denoising, which utilizes the properties of the wavelet transform and the SVD, with consideration of enhancement of the directional features in the original image.

The method is described as follows. First, use the single-level discrete 2D wavelet transform to decompose the noised image into the low-frequency image part and the high-frequency parts (the horizontal, vertical, and diagonal parts), with the edge extracted and retained to avoid edge loss. Then use the SVD to filter the noise of the high-frequency parts with image rotations and the enhancement of the directional features: to filter the diagonal part, one needs first to rotate it 45 degrees and rotate it back after filtering. Finally, reconstruct the image from the low-frequency part and the filtered high-frequency parts, by the inverse wavelet transform, to get the final denoising image.

The experiments show that this method is effective in denoising, while it retains the original edge details, compared with the peak signal-to-noise ratio (PSNR) using the relevant methods, and it has better performance than the recent improved wavelet threshold denoising method in [3].

The reminder of the paper is arranged as follows. Section 2 is the principles of wavelet transform and SVD. Section 3 is the denoising method. Section 4 is the numerical experiments. Section 5 is the conclusion.

#### 2. Principles of Wavelet Transform and SVD

##### 2.1. Wavelet Transform and Threshold Denoising

For a two-dimensional image, using the wavelet transform, one gets the four parts: the low-frequency image part and the three high-frequency parts (the horizontal, vertical, and diagonal parts), as shown in Figure 1.

Common threshold functions for denoising are as follows: (1) hard threshold functions and (2) soft threshold functions. He et al. proposed an improved wavelet threshold denoising method [3]: a new index is proposed to measure the interscale correlation of wavelet coefficients aswhere represents the wavelet coefficients in location point . is an index measuring the interscale correlation of wavelet coefficients. Because the universal threshold often has the risk of overkilling the useful information, therefore , the threshold needs to be shrunk to retain location point , and when , the threshold of location point remains invariant. A new threshold determination method considering interscale correlation is presented aswhere is universal threshold.

##### 2.2. Introduction of SVD Denoising

The SVD of matrix is defined aswhere and are the left singular matrix and the right singular matrix of , the column vector of and is the left singular vectors and the right singular vector of . The rank of is . is the singular value matrix, with the diagonal elements (the nonzero singular values), as . Then can be expressed asby removing zero singular value of . The singular value reflects the energy distribution of the matrix. We choose the appropriate number of or singular value threshold to reconstruct

For an image with only vertical lines, the image signal energy after SVD is concentrated in a few larger singular values. This is the directional characteristics of SVD.

#### 3. Image Denoising in This Paper

The algorithm of this paper is as follows (Figure 2):(1)use the single-level or multilevel discrete 2D wavelet transform to decompose the noised image into the low-frequency image part and the high-frequency parts (the horizontal, vertical, and diagonal parts);(2)use the edge detection operator to extract the edges (Section 3.1);(3)use the SVD to filter the noise of the high-frequency parts, with the singular values determined by the PSNR of the image (Section 3.2);(4)use the SVD to filter the noise of the high-frequency parts with image rotations and the enhancement of the directional features: to filter the diagonal part, one needs to first rotate it 45 degrees and rotate it back after filtering (one illustration is shown in Figure 3);(5)finally, reconstruct the image from the low-frequency part and the filtered high-frequency parts, using the inverse wavelet transform by the inverse wavelet transform, to get the final denoising image.

**(a) Original image**

**(b) Image (a) with noise**

**(c) Rotating the image (b)**

**(d) SVD denoising image**

##### 3.1. High-Frequency Edge Extraction

Because the high-frequency parts correspond to the edge portion of the original image, we can use the single-direction edge detection operator to extract the edge details: using Sobel edge detection operator to detect the horizontal and vertical direction of the high-frequency information and using Roberts edge detection operator to detect the diagonal direction of the high-frequency information. Edge detection operators convolve with image through a series of windows with directions. These operators are matrices. The following are the Sobel operators for edge detection in the horizontal and vertical directions, respectively:and the following are the Roberts operators for edge detection in 45° and 135° directions, respectively:

Therefore, each operator can be used to match the two mutually perpendicular directions that are needed for different edge direction.

##### 3.2. Determination of the Reconstruction of Singular Values

The number of the singular values used for the reconstruction has options. The traditional approach is to choose the singular value threshold ; for Gauss white noise [4], the noise variance is . This formula gives the upper limit value only, but not the optimal threshold. This paper uses the greatest gradient of the PSNR as the measure criterion:where is the pixel gray value of the denoising image, is the pixel gray value of the original image, is the total number of singular values, and is the PSNR of the reconstructed image using the total number of the singular values.

#### 4. Simulations

In the simulations, two images are used, to verify the validity of the method; the size of both images is , as shown in Figures 4(a) and 5(a). Figures 4(c)–4(m) are the denoised images for Figure 4(a).

For the one-layer wavelet decomposition, Figures 4(f-1)–4(k), respectively, illustrate the results of the following:(i)the image denoising method in this paper (with db3 and sym3 wavelets in the wavelet transform, resp.);(ii)the traditional wavelet hard threshold method (with db3 and sym3 wavelets in the wavelet transform, resp.);(iii)the traditional wavelet soft threshold method (with db3 and sym3 wavelets in the wavelet transform, resp.);(iv)the hard thresholding method proposed in [3] (with db3 and sym3 wavelets in the wavelet transform, resp.);(v)the soft thresholding method proposed in [3] (with db3 and sym3 wavelets in the wavelet transform, resp.);(vi)the SVD filtering method.

The results of the image denoising method in this paper, using the 2~3-layer wavelet decomposition (with db3 wavelet), are obtained in Figures 4(l) and 4(m).

For Figure 5(a), the results are shown in Figures 5(c)–5(i), using the method of this paper; as a comparison, Figures 5(j) and 5(k) are the results using the traditional wavelet hard threshold and soft threshold methods and the method proposed in the literature [3].

#### 5. Conclusions

This paper provides an image denoising method, using the wavelet transform and the singular value decomposition (SVD), with the enhancement of the directional features. The experiments show that this method is effective in denoising, while it retains the original edge details, compared with the results using the relevant methods, and it also has a better performance than the recent wavelet threshold denoising method in [3].

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### References

- M. K. Mihçak, I. Kozintsev, K. Ramchandran, and P. Moulin, “Low-complexity image denoising based on statistical modeling of wavelet coefficients,”
*IEEE Signal Processing Letters*, vol. 6, no. 12, pp. 300–303, 1999. View at: Publisher Site | Google Scholar - J.-L. Starck, E. J. Candès, and D. L. Donoho, “The curvelet transform for image denoising,”
*IEEE Transactions on Image Processing*, vol. 11, no. 6, pp. 670–684, 2002. View at: Publisher Site | Google Scholar | MathSciNet - C. He, J. Xing, J. Li, Q. Yang, and R. Wang, “A new wavelet threshold determination method considering interscale correlation in signal denoising,”
*Mathematical Problems in Engineering*, vol. 2015, Article ID 280251, 9 pages, 2015. View at: Publisher Site | Google Scholar - R. Yang and M. Ren, “Wavelet denoising using principal component analysis,”
*Expert Systems with Applications*, vol. 38, no. 1, pp. 1073–1076, 2011. View at: Publisher Site | Google Scholar - A. Fathi and A. R. Naghsh-Nilchi, “Efficient image denoising method based on a new adaptive wavelet packet thresholding function,”
*IEEE Transactions on Image Processing*, vol. 21, no. 9, pp. 3981–3990, 2012. View at: Publisher Site | Google Scholar | MathSciNet - L.-K. Shark and C. Yu, “Denoising by optimal fuzzy thresholding in wavelet domain,”
*Electronics Letters*, vol. 36, no. 6, pp. 581–582, 2000. View at: Publisher Site | Google Scholar - T.-H. Yi, H.-N. Li, and X.-Y. Zhao, “Noise smoothing for structural vibration test signals using an improved wavelet thresholding technique,”
*Sensors*, vol. 12, no. 8, pp. 11205–11220, 2012. View at: Publisher Site | Google Scholar - F. Autin, J.-M. Freyermuth, and R. von Sachs, “Combining thresholding rules: a new way to improve the performance of wavelet estimators,”
*Journal of Nonparametric Statistics*, vol. 24, no. 4, pp. 905–922, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - N. Sae-Bae and S. Udomhunsakul, “Adaptive block-based singular value decomposition filtering,” in
*Proceedings of the Computer Graphics, Imaging and Visualisation (CGIV '07)*, pp. 298–303, IEEE, Bangkok, Thailand, August 2007. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2015 Min Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.