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

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.


Introduction
Denoising is important in image processing.There are many techniques for denoising, such as using the wavelet transform [1][2][3][4][5][6][7][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][2][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 highfrequency 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 lowfrequency 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.

Principles of Wavelet Transform and SVD
2.1.Wavelet Transform and Threshold Denoising.For a twodimensional image, using the wavelet transform, one gets the four parts: the low-frequency image part and the three highfrequency 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 as where (: , ) 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 () ∈ [0, ), 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 as where  is universal threshold.

Introduction of SVD Denoising.
The SVD of matrix  ∈   1 × 2 ( 1 ≥  2 ) is defined as where  = ( 1 ,  2 , . . .,  1 ) ∈   1 × 2 and  = (V 1 , V 2 , . . ., V 2 ) ∈   1 × 2 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 rank() =  ( ≤  2 ). ∈   1 × 2 is the singular value matrix, with the diagonal elements (the nonzero singular values), as Then  can be expressed as by 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.

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.(e-1) the one-layer horizontal high-frequency image using the SVD denoising; (e-2) the one-layer vertical high-frequency image using the SVD denoising; (e-3) the one-layer diagonal high-frequency image with rotation of 45 degrees; (e-4) the denoised image of (e-3) using the SVD denoising; (e-5) the image of (e-4) with rotation restore; (e-6) the edge of the onelayer horizontal high-frequency image; (e-7) the edge of the one-layer vertical high-frequency image; (e-8) the edge of the one-layer diagonal high-frequency image; (e-9) the one-layer horizontal high-frequency image after the fusion; (e-10) the one-layer vertical high-frequency image after fusion; (e-11) the one-layer diagonal high-frequency image after fusion; (f-1) the denoised image using the method in this paper (db3 wavelet); (f-2) the denoising image using the method in this paper (sym3 wavelet); (g-1) the denoised image based on the traditional wavelet hard threshold (db3 wavelet); (g-2) the denoised image based on the traditional wavelet hard threshold (sym3 wavelet); (h-1) the denoised image based on the traditional wavelet soft threshold (db3 wavelet); (h-2) the denoised image based on the traditional wavelet soft threshold (sym3 wavelet); (i-1) the denoised image using the method in [3] for hard thresholding (db3 wavelet); (i-2) the denoised image using the method in [3] for hard thresholding (sym3 wavelet); (j-1) the denoised image using the method in [3] for soft thresholding (db3 wavelet); (j-2) the denoised image using the method in [3] for soft thresholding (sym3 wavelet); (k) the denoised image based on traditional SVD; (l) the denoised image using the method in this paper with the two-layer wavelet decomposition (db3 wavelet); (m) the denoised image using the method in this paper with the three-layer wavelet decomposition (db3 wavelet); (n) the PSNR results for different methods; (o) the PSNR results with 1∼3-layer wavelet decomposition (db3 wavelet); (p-1) the PSNR with three denoising methods ( 2 = 0.001, db3 wavelet); (p-2) the PSNR with three denoising methods ( 2 = 0.001, sym3 wavelet); (q-1) the PSNR with three denoising methods ( 2 = 0.005, db3 wavelet); (q-2) the PSNR with three denoising methods ( 2 = 0.005, sym3 wavelet); (r-1) the PSNR with three denoising methods ( 2 = 0.01, db3 wavelet); (r-2) the PSNR with three denoising methods ( 2 = 0.01, sym3 wavelet).

High-Frequency Edge
high-frequency information and using Roberts edge detection operator to detect the diagonal direction of the highfrequency information.Edge detection operators convolve with image through a series of windows with directions.These operators are 3×3 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.

Determination of the Reconstruction of Singular Values.
The number of the singular values used for the reconstruction has 1 ∼  options.The traditional approach is to choose the singular value threshold  ≤ √ 1  2 ; for Gauss white noise [4], the noise variance is  2 .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,  0 (, ) is the pixel gray value of the original image, () is the total number of singular values, and PSNR() is the PSNR of the reconstructed image using the total number of the  th singular values.

Simulations
In the simulations, two images are used, to verify the validity of the method; the size of both images is 512 × 512, as shown in Figures 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].

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].

Figure 1 :
Figure 1: Wavelet decomposition: the upper left corner is the smoothing approximation (low-frequency LL).HL: the lower left corner for the vertical part, LH: the upper right corner for the horizontal part, and HH: the lower right corner.
Extraction.Because the highfrequency 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