# Image Denoising Using Singular Value Difference in the Wavelet Domain

Authors: Min Wang; Wei Yan; Shudao Zhou
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Singular value (SV) difference is the difference in the singular values between a noisy image and the original image; it varies regularly with noise intensity. This paper proposes an image denoising method using the singular value difference in the wavelet domain. First, the SV difference model is generated for different noise variances in the three directions of the wavelet transform and the noise variance of a new image is used to make the calculation by the diagonal part. Next, the single-level discrete 2-D wavelet transform is used to decompose each noisy image into its low-frequency and high-frequency parts. Then, singular value decomposition (SVD) is used to obtain the SVs of the three high-frequency parts. Finally, the three denoised high-frequency parts are reconstructed by SVD from the SV difference, and the final denoised image is obtained using the inverse wavelet transform. Experiments show the effectiveness of this method compared with relevant existing methods.


---

## Body

## 1. Introduction

Wavelet transform [1–6] and singular value decomposition (SVD) [7] have been widely used as transform domain methods [8] in image denoising. Wavelet transform can highlight the detailed information of an image, while SVD is generally used as a kind of nonlinear filter. Recently, several image denoising methods based on SVD have been proposed. In [9], a denoising algorithm based on adaptive SVD in the wavelet domain is proposed. An adaptive representation method is proposed [10] using the K-means and singular value decomposition (K-SVD), which uses a greedy algorithm to learn an overcomplete dictionary for image representation and denoising. An image denoising method based on K-SVD is proposed [11], which finds a few atoms from the dictionary having the best linear combination to represent each subblock. In [12], a denoising algorithm using adaptive SVD (ASVD) is proposed. A patch-based weighted-SVD denoising method is proposed [13] with feature retention. In [14], an image denoising method is proposed using higher order SVD (HOSVD). A denoising algorithm using SVD and wavelet transform is proposed [15], which enhances directional features. In [16], a denoising algorithm using SVD and Ridgelet transform is proposed; Ridgelet transform is slightly inferior in homogenous regions of nontextured images. A singular value decomposition (APBSVD) [17] that preserves the edge structure and avoids blurriness is proposed. In [18], a denoising method based on spatially adaptive iterative singular value thresholding (SAIST) is proposed. A new wavelet threshold determination method considering interscale correlation in signal denoising is proposed [19].Singular value (SV) difference, the difference in the singular values between a noisy image and the original image, varies regularly with noise variance. Therefore, we can use a standard SV difference to generate a new SV difference function using the different noise variance estimates. Then, the denoised image can be expressed as the SV of the noisy image minus the estimated SV difference. In this paper, an image denoising method is proposed using the SV difference in the wavelet domain. First, establish the SV difference model with different noise variances in the three directions after wavelet transform and estimate the noise variance of a new image. Then, use the single-level discrete 2-D wavelet transform to decompose each noisy image into a low-frequency and three high-frequency image parts, and use the SVD to obtain the SV of three high-frequency parts. Finally, reconstruct three denoised high-frequency parts by SVD with the SV difference and obtain the final denoised image by the inverse wavelet transform.The rest of this paper is organized as follows. In Section2, we briefly review the principles of wavelet transform and SVD denoising. In Section 3, we present the proposed denoising method in detail. In Section 4, we give some numerical experiments and performance analysis. Finally, Section 5 contains the concluding remarks.


## 2. Mathematical Preliminaries

### 2.1. Wavelet Transform and Variance Estimating

Using the wavelet transform, an image can be decomposed into its low-frequency and high-frequency parts. The low-frequency part represents the approximate energy of an image, while the three high-frequency parts represent the detailed information of an image, including the horizontal, vertical, and diagonal parts.New wavelet coefficients are estimated by thresholding the original wavelet coefficients. Common threshold functions include hard and soft threshold functions. These thresholds are determined by the varianceσof the noise, which is estimated as(1)σ=medianWj,10.6745,where Wj,1 is the median of the wavelet coefficients in the diagonal part.


### 2.2. SVD Denoising

Every two-dimensional imageAof size M×N(M≥N) can be decomposed into three matrices by SVD:(2)A=USVT,where U and V are the left and right singular matrices of A, with column vectors u⃑i and v⃑i, respectively. The rank of A is R(R≤N), and when the diagonal SV matrix S=diag⁡s1,s2,…,sR is a nonnegative matrix, the nonzero singular values can be arranged as s1≥s2≥⋯≥sR>0. These SVs reflect the energy distribution of the image, and si can be considered the representation coefficient. Therefore, image A can be expressed by ignoring SVs having value zero as follows:(3)A=∑i=1RsiuiviT.Wang et al. proposed an improved image denoising method based on wavelet and SVD transforms using the directional features [15]: use the SVD to filter the noise of the high-frequency parts with image rotations and the enhancement of the directional features; then rotate it back after filtering.


### 2.3. The Features of SV Difference

For a two-dimensional imageA, a noisy image B can be expressed as(4)B=A+N,where N represents the random noise.SV difference represents the difference in SVs between the noisy and original images, defined as(5)g=sB-sA,where sB and sA are the SVs of noisy image B and denoised image A, respectively.When noise is added to a clear image, it can be regarded as disturbance. This disturbance is bounded and is only related to the features of the noise, independent of the original image. Therefore, SV difference varies regularly with the noise variance. That is, if we know a certain SV difference with one variance and this SV difference has a regular relationship with variance, we can establish the SV difference model for an arbitrary variance. Moreover, for denoising, we can use the SV of the noisy image minus the SV difference to represent the SV of the denoised image. However, this denoising method based on SV difference applied to the entire image is suitable for flat images but very insensitive to detailed features. Therefore, we can consider applying SV difference on the detail-retained images after wavelet transform to remove the noise.To the original images, shown in Figures1(a) and 1(e), white Gaussian noise of various intensities are added, generating a series of SV differences of the three high-frequency parts. Figure 1 shows that the shape of these SV difference curves varies regularly with the intensity of the noise.Figure 1
SV difference curves of different images with different added noise: (a) first original image; (b) SV difference curves of (a) of the horizontal part; (c) SV difference curves of (a) of the vertical part; (d) SV difference curves of (a) of the diagonal part; (e) second original image; (f) SV difference curves of (e) of the horizontal part; (g) SV difference curves of (e) of the vertical part; (h) SV difference curves of (e) of the diagonal part.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)With the upper left corner of Figure1(e) as a reference, we extract the first image (of size 200×200) and then increase the length and width by 50 pixels, thus obtaining seven images to be considered as original images; we then add white Gaussian noise with variance of 0.02 to these original images. Figure 2 shows that the shape of the SV difference curves changes only with the size of the image, and it also varies regularly.Figure 2
Effects of different image sizes on SV difference: (a) SV difference curves of the horizontal part; (b) SV difference curves of the vertical part; (c) SV difference curves of the diagonal part.
(a)
(b)
(c)Figure3 contains seven SV difference curves, where the noisy images (with variance of 0.04) originate from Figure 1(e), extracting seven images (the size of 200×200) from upper left to lower right corner of each shift 50 pixels. Figure 3 shows that the SV difference does not change with the position of the image.Figure 3
SV difference curves at different positions in the image: (a) SV difference curves of the horizontal part; (b) SV difference curves of the vertical part; (c) SV difference curves of the diagonal part.
(a)
(b)
(c)


## 2.1. Wavelet Transform and Variance Estimating

Using the wavelet transform, an image can be decomposed into its low-frequency and high-frequency parts. The low-frequency part represents the approximate energy of an image, while the three high-frequency parts represent the detailed information of an image, including the horizontal, vertical, and diagonal parts.New wavelet coefficients are estimated by thresholding the original wavelet coefficients. Common threshold functions include hard and soft threshold functions. These thresholds are determined by the varianceσof the noise, which is estimated as(1)σ=medianWj,10.6745,where Wj,1 is the median of the wavelet coefficients in the diagonal part.


## 2.2. SVD Denoising

Every two-dimensional imageAof size M×N(M≥N) can be decomposed into three matrices by SVD:(2)A=USVT,where U and V are the left and right singular matrices of A, with column vectors u⃑i and v⃑i, respectively. The rank of A is R(R≤N), and when the diagonal SV matrix S=diag⁡s1,s2,…,sR is a nonnegative matrix, the nonzero singular values can be arranged as s1≥s2≥⋯≥sR>0. These SVs reflect the energy distribution of the image, and si can be considered the representation coefficient. Therefore, image A can be expressed by ignoring SVs having value zero as follows:(3)A=∑i=1RsiuiviT.Wang et al. proposed an improved image denoising method based on wavelet and SVD transforms using the directional features [15]: use the SVD to filter the noise of the high-frequency parts with image rotations and the enhancement of the directional features; then rotate it back after filtering.


## 2.3. The Features of SV Difference

For a two-dimensional imageA, a noisy image B can be expressed as(4)B=A+N,where N represents the random noise.SV difference represents the difference in SVs between the noisy and original images, defined as(5)g=sB-sA,where sB and sA are the SVs of noisy image B and denoised image A, respectively.When noise is added to a clear image, it can be regarded as disturbance. This disturbance is bounded and is only related to the features of the noise, independent of the original image. Therefore, SV difference varies regularly with the noise variance. That is, if we know a certain SV difference with one variance and this SV difference has a regular relationship with variance, we can establish the SV difference model for an arbitrary variance. Moreover, for denoising, we can use the SV of the noisy image minus the SV difference to represent the SV of the denoised image. However, this denoising method based on SV difference applied to the entire image is suitable for flat images but very insensitive to detailed features. Therefore, we can consider applying SV difference on the detail-retained images after wavelet transform to remove the noise.To the original images, shown in Figures1(a) and 1(e), white Gaussian noise of various intensities are added, generating a series of SV differences of the three high-frequency parts. Figure 1 shows that the shape of these SV difference curves varies regularly with the intensity of the noise.Figure 1
SV difference curves of different images with different added noise: (a) first original image; (b) SV difference curves of (a) of the horizontal part; (c) SV difference curves of (a) of the vertical part; (d) SV difference curves of (a) of the diagonal part; (e) second original image; (f) SV difference curves of (e) of the horizontal part; (g) SV difference curves of (e) of the vertical part; (h) SV difference curves of (e) of the diagonal part.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)With the upper left corner of Figure1(e) as a reference, we extract the first image (of size 200×200) and then increase the length and width by 50 pixels, thus obtaining seven images to be considered as original images; we then add white Gaussian noise with variance of 0.02 to these original images. Figure 2 shows that the shape of the SV difference curves changes only with the size of the image, and it also varies regularly.Figure 2
Effects of different image sizes on SV difference: (a) SV difference curves of the horizontal part; (b) SV difference curves of the vertical part; (c) SV difference curves of the diagonal part.
(a)
(b)
(c)Figure3 contains seven SV difference curves, where the noisy images (with variance of 0.04) originate from Figure 1(e), extracting seven images (the size of 200×200) from upper left to lower right corner of each shift 50 pixels. Figure 3 shows that the SV difference does not change with the position of the image.Figure 3
SV difference curves at different positions in the image: (a) SV difference curves of the horizontal part; (b) SV difference curves of the vertical part; (c) SV difference curves of the diagonal part.
(a)
(b)
(c)


## 3. Proposed Image Denoising Procedure

The steps of the proposed image denoising procedure using SV difference in the wavelet domain are given below (Figure4):(1)
Generate the SV difference model with different noise variance in the three high-frequency parts of the wavelet transform (Section3.1).(2)
Use the single-level discrete 2-D wavelet transform to decompose the new noisy image into its low-frequency part and three high-frequency parts, and then calculate the size and estimate the variance of the noise to get the new SV differences for this noisy image.(3)
Use SVD to get the noisy SV of the high-frequency parts, and subtract the above SV difference to obtain the denoised SV (Section3.2).(4)
Reconstruct the denoised high-frequency parts from the estimated denoised SVs and the corresponding eigenvectors by SVD.(5)
Finally, reconstruct the final denoised image by the inverse wavelet transform.Figure 4
Flow chart of this algorithm.


### 3.1. Establish the SV Difference Model

Figures1 and 2 show that the SV difference curve is sensitive to the size of the image and the added noise. We introduce two parameters related to the size of the image and the intensity of the added noise to adjust the SV difference function:(6)gnx=pnmqnσfnlmx,where n=H,V,D represents the horizontal, vertical, or diagonal direction, l and m are the size of the standard and processing images, respectively, p(m) and q(σ) represent the longitudinal shrinkage coefficients of different size and variance, respectively, and f(x) represents the SV difference function of the standard image, which can be obtained from the SV difference of the middle 200×200 pixels of the noise image (σ2=0.02) and the original image; x is the number of the SV.We use Figure1(a) as the test image; the fitting function of f(x), as shown in Figure 5, with 5-degree polynomial in the three directions can be estimated as(7)fHx=-1.5524×10-7x5+4.5650×10-5x4-4.9961×10-3x3+2.5727×10-1x2-11.6241x+587.4337fVx=-1.5199×10-7x5+4.3501×10-5x4-4.6762×10-3x3+2.4151×10-1x2-11.5083x+594.6933fDx=-1.6062×10-7x5+5.3118×10-5x4-6.5353×10-3x3+3.7854×10-1x2-16.2181x+676.9813.Figure 5
The standard SV difference functions in three directions: (a) SV difference curves of the horizontal part; (b) vertical SV difference curves of the vertical part; (c) diagonal SV difference curves.
(a)
(b)
(c)The longitudinal shrinkage coefficientsq(σ) associated with the variance are fitted by the maximum values in Figure 1, as shown in Figure 6. According to the shape of SV difference curves varies regularly with the variance of the noise of Figures 1(b)–1(d), the longitudinal shrinkage coefficients q(σ) in the three directions can be calculated by the fitting function of maximum values with 3-degree polynomial divided by the fitting function of standard variance (σ2=0.02)(8)qHσ=1513505σ3-377364σ2+39621σ+2381513505×0.023-377364×0.022+39621×0.02+238qVσ=939474σ3-268453σ2+34042σ+325939474×0.023-268453×0.022+34042×0.02+325qDσ=419329σ3-158053σ2+26564σ+592419329×0.023-158053×0.022+26564×0.02+592.Figure 6
Parameter functionsq(σ) in three directions: (a) maximum values curves of the horizontal part; (b) maximum values curves of the vertical part; (c) maximum values curves of the diagonal part.
(a)
(b)
(c)The longitudinal shrinkage coefficientsp(m) associated with the size are fitted by the maximum values in Figure 2, as shown in Figure 7. According to the shape of SV difference curves vary regularly with the size of the image of Figures 2(a)–2(c); the longitudinal shrinkage coefficients p(m) in the three directions can be calculated by the fitting function of maximum values with 3-degree polynomial divided by the fitting function of standard size (200):(9)pHm=2.1047×10-7×m3-1.0926×10-3×m2+1.8138×m+232.65682.1047×10-7×2003-1.0926×10-3×2002+1.8138×200+232.6568pVm=-4.7204×10-6×m3+3.2063×10-3×m2+0.6879×m+355.8800-4.7204×10-6×2003+3.2063×10-3×2002+0.6879×200+355.8800pDm=4.4546×10-6×m3-6.1128×10-3×m2+3.8470×m+92.75214.4546×10-6×2003-6.1128×10-3×2002+3.8470×200+92.7521.Figure 7
Parameter functionsp(m) in three directions: (a) maximum values curves of the horizontal part; (b) maximum values curves of the vertical part; (c) maximum values curves of the diagonal part.
(a)
(b)
(c)These coefficients of functions (6)–(9) are calculated from Figure 1(e) as the test image; because the noise is random, these coefficients will change slightly. Therefore, we can choose a random clear figure as the new testing standard figure to determine the coefficients of fitting function according to the above steps, or we can take the mean of the coefficients from multiple images as the final factor.


### 3.2. Estimate the Denoised SV

We get the approximate SV sequence of the original image through subtraction of the SV difference from the SVs of the noisy image:(10)snA=snB-gn,where snA and snB represent the SVs of the denoised and noisy images, respectively; gn represents the SV difference; n∈H,V,D represents the horizontal, vertical, or diagonal direction. We then reconstruct the denoised image by SVD.


## 3.1. Establish the SV Difference Model

Figures1 and 2 show that the SV difference curve is sensitive to the size of the image and the added noise. We introduce two parameters related to the size of the image and the intensity of the added noise to adjust the SV difference function:(6)gnx=pnmqnσfnlmx,where n=H,V,D represents the horizontal, vertical, or diagonal direction, l and m are the size of the standard and processing images, respectively, p(m) and q(σ) represent the longitudinal shrinkage coefficients of different size and variance, respectively, and f(x) represents the SV difference function of the standard image, which can be obtained from the SV difference of the middle 200×200 pixels of the noise image (σ2=0.02) and the original image; x is the number of the SV.We use Figure1(a) as the test image; the fitting function of f(x), as shown in Figure 5, with 5-degree polynomial in the three directions can be estimated as(7)fHx=-1.5524×10-7x5+4.5650×10-5x4-4.9961×10-3x3+2.5727×10-1x2-11.6241x+587.4337fVx=-1.5199×10-7x5+4.3501×10-5x4-4.6762×10-3x3+2.4151×10-1x2-11.5083x+594.6933fDx=-1.6062×10-7x5+5.3118×10-5x4-6.5353×10-3x3+3.7854×10-1x2-16.2181x+676.9813.Figure 5
The standard SV difference functions in three directions: (a) SV difference curves of the horizontal part; (b) vertical SV difference curves of the vertical part; (c) diagonal SV difference curves.
(a)
(b)
(c)The longitudinal shrinkage coefficientsq(σ) associated with the variance are fitted by the maximum values in Figure 1, as shown in Figure 6. According to the shape of SV difference curves varies regularly with the variance of the noise of Figures 1(b)–1(d), the longitudinal shrinkage coefficients q(σ) in the three directions can be calculated by the fitting function of maximum values with 3-degree polynomial divided by the fitting function of standard variance (σ2=0.02)(8)qHσ=1513505σ3-377364σ2+39621σ+2381513505×0.023-377364×0.022+39621×0.02+238qVσ=939474σ3-268453σ2+34042σ+325939474×0.023-268453×0.022+34042×0.02+325qDσ=419329σ3-158053σ2+26564σ+592419329×0.023-158053×0.022+26564×0.02+592.Figure 6
Parameter functionsq(σ) in three directions: (a) maximum values curves of the horizontal part; (b) maximum values curves of the vertical part; (c) maximum values curves of the diagonal part.
(a)
(b)
(c)The longitudinal shrinkage coefficientsp(m) associated with the size are fitted by the maximum values in Figure 2, as shown in Figure 7. According to the shape of SV difference curves vary regularly with the size of the image of Figures 2(a)–2(c); the longitudinal shrinkage coefficients p(m) in the three directions can be calculated by the fitting function of maximum values with 3-degree polynomial divided by the fitting function of standard size (200):(9)pHm=2.1047×10-7×m3-1.0926×10-3×m2+1.8138×m+232.65682.1047×10-7×2003-1.0926×10-3×2002+1.8138×200+232.6568pVm=-4.7204×10-6×m3+3.2063×10-3×m2+0.6879×m+355.8800-4.7204×10-6×2003+3.2063×10-3×2002+0.6879×200+355.8800pDm=4.4546×10-6×m3-6.1128×10-3×m2+3.8470×m+92.75214.4546×10-6×2003-6.1128×10-3×2002+3.8470×200+92.7521.Figure 7
Parameter functionsp(m) in three directions: (a) maximum values curves of the horizontal part; (b) maximum values curves of the vertical part; (c) maximum values curves of the diagonal part.
(a)
(b)
(c)These coefficients of functions (6)–(9) are calculated from Figure 1(e) as the test image; because the noise is random, these coefficients will change slightly. Therefore, we can choose a random clear figure as the new testing standard figure to determine the coefficients of fitting function according to the above steps, or we can take the mean of the coefficients from multiple images as the final factor.


## 3.2. Estimate the Denoised SV

We get the approximate SV sequence of the original image through subtraction of the SV difference from the SVs of the noisy image:(10)snA=snB-gn,where snA and snB represent the SVs of the denoised and noisy images, respectively; gn represents the SV difference; n∈H,V,D represents the horizontal, vertical, or diagonal direction. We then reconstruct the denoised image by SVD.


## 4. Simulations

To validate the proposed method, we conducted several experiments on two test images with added different white Gaussian noise; the results are presented in this section. The noise mean is 0, and the normalized variance is from 0.001 to 0.1. The size of both test images is512×512. Figures 8(a), 9(a), 11(a), 12(a), 13(a), and 14(a) are noisy images with the added noise of different variances.Figure 8
(a) The noisy image (σ2=0.001); (b) the low-frequency part of (a); (c) the horizontal part of (a); (d) the vertical part of (a); (e) the diagonal part of (a); (f) the SV of noisy horizontal part; (g) the SV of noisy vertical part; (h) the SV of noisy diagonal part; (i) the SV difference of the horizontal part; (j) the SV difference of the vertical part; (k) the SV difference of the diagonal part; (l) the denoised SV and original SV of the horizontal part; (m) the denoised SV and original SV of the vertical part; (n) the denoised SV and original SV of the diagonal part; (o) the reconstructed horizontal part; (p) the reconstructed vertical part; (q) the reconstructed diagonal part; (r) the denoised image using the method in this paper; (s) the denoised image based on the conventional wavelet hard threshold; (t) the denoised image based on the conventional wavelet soft threshold; (u) the denoised image based on the conventional SVD; (v) the denoised image using the method in [15]; (w) the denoised image using the method in [19] for hard thresholding; (x) the denoised image using the method in [19] for soft thresholding.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(q)
(r)
(s)
(t)
(u)
(v)
(w)
(x)Figure 9
(a) The noisy image (σ2=0.01); (b) the low-frequency part of (a); (c) the horizontal part of (a); (d) the vertical part of (a); (e) the diagonal part of (a); (f) the SV of noisy horizontal part; (g) the SV of noisy vertical part; (h) the SV of noisy diagonal part; (i) the SV difference of the horizontal part; (j) the SV difference of the vertical part; (k) the diagonal SV difference; (l) the denoised SV and original SV of the horizontal part; (m) denoised SV and original SV of the vertical part; (n) the diagonal denoised SV and original SV; (o) the reconstructed horizontal part; (p) the reconstructed vertical part; (q) the reconstructed diagonal part; (r) the denoised image using the method in this paper; (s) the denoised image based on the conventional wavelet hard threshold; (t) the denoised image based on the conventional wavelet soft threshold; (u) the denoised image based on the conventional SVD; (v) the denoised image using the method in [15]; (w) the denoised image using the method in [19] for hard thresholding; (x) the denoised image using the method in [19] for soft thresholding.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)
(p)
(q)
(r)
(s)
(t)
(u)
(v)
(w)
(x)For the one-layer wavelet decomposition with db3 wavelets, Figures8(r)–8(x) and 9(r)–9(x), 11(e)–11(k), and 12(b)–12(h) and 14(b)–14(h), respectively, illustrate the results of:(i)
the denoising method proposed in this paper,(ii)
the conventional wavelet hard threshold method,(iii)
the conventional wavelet soft threshold method,(iv)
the conventional SVD filtering method,(v)
the denoising method proposed in [15],(vi)
the denoising method proposed in [19] for hard and soft thresholding.We used peak signal to noise ratio (PSNR) and structural similarity index (SSIM) to evaluate the denoising effectiveness.Figures8–14 show that the denoising method proposed in this paper significantly outperforms PSNR and SSIM and provides a better denoising effect with greater detail preservation compared to those methods.Figure 10
(a) PSNR results with four methods of Figures8 and 9; (b) SSIM results with four methods of Figures 8 and 9; (c) PSNR results with three methods of Figures 8-9; (d) SSIM results with three methods of Figures 8 and 9.
(a)
(b)
(c)
(d)Figure 11
(a) The noisy image (σ2=0.01); (b) the horizontal denoised SV and original SV; (c) denoised SV and original SV of the vertical part; (d) the diagonal denoised SV and original SV; (e) the denoised image using the method in this paper; (f) the denoised image based on the conventional wavelet hard threshold; (g) the denoised image based on the conventional wavelet soft threshold; (h) the denoised image based on the conventional SVD; (i) the denoised image using the method in [15]; (j) the denoised image using the method in [19] for hard thresholding; (k) the denoised image using the method in [19] for soft thresholding; (l) PSNR results for different methods; (m) SSIM results for different methods.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)Figure 12
(a) The noisy image (σ2=0.1); (b) the denoised image using the method in this paper; (c) the denoised image based on the conventional wavelet hard threshold; (d) the denoised image based on the conventional wavelet soft threshold; (e) the denoised image based on the conventional SVD; (f) the denoised image using the method in [15]; (g) the denoised image using the method in [19] for hard thresholding; (h) the denoised image using the method in [19] for soft thresholding; (i) PSNR results for different methods; (j) SSIM results for different methods.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)Figure 13
(a) The noisy image (σ2=0.1); (b) the denoised image using the method in this paper; (c) the denoised image based on the conventional wavelet hard threshold; (d) the denoised image based on the conventional wavelet soft threshold; (e) the denoised image based on the conventional SVD; (f) the denoised image using the method in [15]; (g) the denoised image using the method in [19] for hard thresholding; (h) the denoised image using the method in [19] for soft thresholding; (i) PSNR results for different methods; (j) SSIM results for different methods.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)Figure 14
(a) The noisy image (σ2=0.01); (b) the denoised image using the method in this paper; (c) the denoised image based on the conventional wavelet hard threshold; (d) the denoised image based on the conventional wavelet soft threshold; (e) the denoised image based on the conventional SVD; (f) the denoised image using the method in [15]; (g) the denoised image using the method in [19] for hard thresholding; (h) the denoised image using the method in [19] for soft thresholding; (i) PSNR results for different methods; (j) SSIM results for different methods.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)


## 5. Conclusions

This paper provides an image denoising method using the singular value difference in the wavelet domain. Experimental results show that this proposed method is effective in denoising and has a better denoising performance compared with relevant existing methods. This method can be used to process images from remote sensing, from the medical domain, and other common color-image denoising problems. Because the longitudinal shrinkage coefficientsp(m), q(v), and f(x) are calculated by the polynomial fitting, the degree of polynomial selected and the testing value of variance and size all affect the coefficients of polynomial fitting function, which affects the final accuracy of image denoising. Therefore, further research still needs to be done on how to optimize the polynomial and choose the best standard image.



---

*Source: 1542509-2018-01-17.xml*