Advanced
A Study on Wavelet-based Image Denoising Using a Modified Adaptive Thresholding Method
A Study on Wavelet-based Image Denoising Using a Modified Adaptive Thresholding Method
Journal of information and communication convergence engineering. 2012. Mar, 10(1): 45-52
Copyright ©2012, The Korean Institute of Information and Commucation Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : December 28, 2011
  • Accepted : January 24, 2012
  • Published : March 31, 2012
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Gao Yinyu

Abstract
Thedenoising of a natural image corrupted by Gaussian noise is a long established problem in signal or image processing. Today the research is focus on the wavelet domain, especially using the wavelet threshold method. In this paper, a wavelet-based image denoising modified adaptive thresholding method is proposed. The proposed method computes thethreshold adaptively based on the scale level and adaptively estimates wavelet coefficients by using a modified thresholding function that considers the dependency between the parent coefficient and child coefficient and the soft thresholding function at different scales. Experimental results show that the proposed method provides high peak signal-to-noise ratio results and preserves the detailed information of the original image well, resulting in a superior quality image.
Keywords
I . INTRODUCTION
The denoising of a natural image corrupted by additive white Gaussian noise (AWGN) is a classical problem in the signal processing community. The corruption of an image by noise is common during its acquisition or transmission. The aim of denoising is to remove the noise while keeping the signal featuresas much as possible. Traditional algorithms, such as the standard median (SM) filter and mean filter (MF) perform image denoising in the spatial domain [1 - 4 , 11] . In recent years, wavelet transform-based image denoising algorithms have shown remarkable success.
In 1992, Donoho and Johnstone [5] presented a method named wavelet shrinkage, and showed its obvious efficiency for signal denoising and inverse problem solving. In this method, a discrete wavelet transform (DWT) is performed on the noisy signal first. Then with a preset threshold, coefficients with a magnitude smaller than the threshold are set to zero while those with a larger magnitude are kept and used to estimate the noiseless coefficients. Finally, an inverse discrete wavelet transform (IDWT) reconstructs the signal from the estimated coefficients. Later, Donoho [7] proposed a sample thresholding rule that sets all the coefficients smaller than the universal threshold to zero and shrinks the rest of the coefficientsby the threshold (soft thresholding) or leaves them with out change (hard thresholding). Donoho and Johnstone [6] also presented a thresholding method using Stein's risk estimator called Sure Shrink. Sendur and Selesnick [8 , 9] proposed BiShrink for denoising. The BiShrinkage function indicates that the estimated wavelet coefficients dependon the parent coefficients.
In this paper, a modified wavelet based image denoising thresholding method is proposed. A modified thresholding function considering the inter scale dependency between the parent coefficient and child coefficient and conventional soft thresholding function are adaptively used to estimate coefficients in different scales. Then the peak signal to noise ratio (PSNR) is used as a standard evaluation of the results of denoising the images; experimental results show that our proposed method gives a higher PSNR and with fewer visual artifacts compared with other methods.
II. INTERSCALE DEPENDENCY
Fig. 1 illustrates the sub-band regions of the two dimensional (2D) critically sampled wavelet transformat level 3, which can be label as HH k , HL k , LH k , and LL J , where scale k = 1, 2, 3, ..., J [8] . J takes the same value as the level for decomposition.
PPT Slide
Lager Image
Sub-band regions of the critically sampled wavelet transform.
The wavelet coefficients of the k th sub-band will be the parent of the coefficients of the k -1 th sub-band. Thisway, P(S) can be denoted as the parent sub-band of the sub-band S. For example, if S is HH 1 , then P(S) is HH 2 , and if S is HL 2 , then P(S) is HL 3 . According to the inter scale dependency principle, each parent coefficient represents a lower frequency component than its children.
III. TRADITIONAL WAVELET SHRINKAGE
Wavelet thresholding is a simple non-linear technique that operates on one wavelet coefficient at a time. The most popular thresholding methods are soft thresholding and hard thresholding. They are given by (1) and (2), respectively. Soft thresholding is a continuous function while hard thresholding is a discontinuous function which causes some artifacts in a denoised image. Therefore, the soft thresholding method is preferred over hard thresholding [4 - 7] .
PPT Slide
Lager Image
PPT Slide
Lager Image
Where λ is the threshold and ŵ and w are the modified and noisy wavelet coefficients, respectively, and ( i, j ) is the location of wavelet coefficient.
- A. VisuShrink
VisuShrink is thresholding by applying the universal threshold ( λuniv ) that was proposed by Donoho and Johnstone [5 , 6] . It is given by
PPT Slide
Lager Image
Where σn is noise variance and M is the size of the image. In practice, the noise variance is estimated by the median estimator [7] using (4).
PPT Slide
Lager Image
- B. SureShrink
Another well-known method for threshold selection is SureShrink, which is based on the method of Stein’s unbiased risk estimator (SURE) [6] . It is a sub-band adaptive threshold technique. A separate threshold is computed for each sub-band based on Stein’s unbiased risk estimator as follows:
PPT Slide
Lager Image
Where #(·) means the number satisfying the condition and “ ˄ ” means the lesser value of the two parameters.
Then the SureShrink threshold λsure is defined as
PPT Slide
Lager Image
Where λsure denotes the value that minimizes (5).
SureShrink is alocal thresholding method having different thresholds for different levels. SureShrink follows the soft thresholding rule.
- C. BiShrink
In [8] is proposed a Laplace bivariate model for the wavelet coefficients of the useful component of the input image and a bivariate Gaussian model for the wavelet coefficients of the noise. The model captures the inter scale dependence between a wavelet coefficient and its parent. Then Sendur and Selesnick [8 , 9] proposed BiShrink for denoising according to this model. BiShrink is a simple non-linear shrinkage function for wavelet denoising, which generalizes the soft thresholding approach of Donoho and Johnstone [5] and it depends on both the coefficient and its parent, yielding improved results for wavelet-based image denoising.
Let w 2 represent the parent of w 1 . The parent is located at the same geometrical coordinates as the child, but at the successive scale. The problem is formulated in the wavelet domain as y 1 = w 1 + n 1 and y 2 = w 2 + n 2 to take into account the statistical dependencies between a coefficient and its parent. y 1 and y 2 are noisy observations of w 1 and w 2 , and n 1 and n 2 are noise samples [8 - 10] . Then y is
PPT Slide
Lager Image
Where w = ( w 1 , w 2 ), y = ( y 1 , y 2 ) , and n = ( n 1 , n 2 ).
The BiShrink formula is given by
PPT Slide
Lager Image
Here σ is the marginal variance for each wavelet coefficient. ( g ) + is defined as
PPT Slide
Lager Image
IV. PROPOSED METHOD
In this paper, an efficient wavelet based image denoising thresholding method is proposed. The proposed method computes the threshold adaptively based on scale level and uses a modified thresholding function which considers the inter scale dependency between parent coefficient and child coefficient and traditional soft thresholding function to estimate the coefficients located in different scales.
For scales k = 1, 2, 3, ..., J -1, the proposed threshold and thresholding function as follows.
Firstly, proposed method computes the threshold value for this scale using the marginal variance and the neighborhood variance. The marginal variance for the s th coefficient will be estimated using neighboring coefficients in the region N(s) [8 , 9] . Here N(s) is defined as all coefficients within a square-shaped window that is centered at the s th coefficient as illustrated in Fig. 2 .
PPT Slide
Lager Image
Illustration of neighborhood N(s).
Our method also tries to estimate the marginal variance σ 2 for the s th wavelet coefficient. In order to confirm the value of marginal variance, our method firstly calculates the variance of neighborhood σ m 2 .
PPT Slide
Lager Image
Where N is the size of the neighborhood N(s) . Then, the marginal variance can be estimated as (11).
PPT Slide
Lager Image
Using (10) and (11), the parent coefficient marginal variance σ p and child coefficient marginal variance σ c can be calculated when using the same window size for the parent and current sub-bands.
Then, the threshold for scales k = 1, 2, 3, ..., J -1, is defined as
PPT Slide
Lager Image
Here L stands for the level chosen for decomposition.
Secondly, a modified thresholding method is proposed to estimate these noisy coefficients. In order to improve the denoising results, our method tries to reduce the fixed-bias, but it cannot be reduced to zero. The expression of the proposed thresholding function is as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Fig. 3 shows the plot of our modified thresholding function. Our results clearly show that the estimated value should depend on the parent value.
In order to show clearly our proposed thresholding function, the relationship between the noisy child and estimated coefficient by a 2D figure when the noisy parent is 0 is plotted in Fig. 4 .
For the coarsest scale k = J , an adaptive sub-band dependent threshold is used for image denoising via wavelet soft thresholding. The adaptive threshold λJ can be estimated as
PPT Slide
Lager Image
Here, σ s is the noise-free signal standard deviation ateach high sub-band(HH J , HL J , LH J ) for scale J , and σ s is computed by an approximate maximum likelihood (ML) estimator.
PPT Slide
Lager Image
Where Ns 2 is the size of each sub-band.
PPT Slide
Lager Image
Proposed thresholding function.
PPT Slide
Lager Image
Illustration of relationship between noisy child and estimated coefficient.
Next, using the soft thresholding function expressed as (1) to estimate the coefficients at scale k = J , the estimated coefficient for scale J is then defined as ŵ j .
Thus, the finally estimated coefficient ŵ can be summarized as follows:
PPT Slide
Lager Image
V. EXPERIMENTAL RESULT ANALYSES
To verify the validity of the proposed method, it is compared with the Mean filter [12], VisuShrink that chooses the universal soft thresholding denoising technique, SureShrink, and BiShrink. The wavelet transform used in all the wavelet domain methods is Db8 in 3-level ( L ) decompositions. For MF and the proposed method, a 5×5 window size ( N ( s )) is used to process the images. And the 512×512 grayscale images Lena and Boat, shown in Fig. 5 , are used for experimentation.
PPT Slide
Lager Image
The original test images. (a) Lena, (b) boat.
The performance of the estimators is measured by using the classical PSNR.
PPT Slide
Lager Image
Where M is the size of the image. O is the original image, Y is the restored image, and G is the gray level of the image.
Figs. 6 and 7 show the visual results of the simulation for
Denoisingresults (PSNRs) of different methods for Lenaand boat.PSNR: peak signal to noise ratio.
PPT Slide
Lager Image
Denoisingresults (PSNRs) of different methods for Lenaand boat. PSNR: peak signal to noise ratio.
subjective visual quality comparison, where the standard deviation of the noise added to the original images is 30. In order to compare the results more clearly, the distribution map of the absolute difference between the original image and the processed image is denoted as Ad and shown in each figure. In this paper, we choose the interval value as 20 to plot the distribution map. The absolute difference can be formulated as
PPT Slide
Lager Image
Figs. 6 a and 7 aare the noisy images that was corrupted by AWGN with σ =30. Figs. 6 b-f and 7 b-f show the results of the restoration by VisuShrink, SureShrink, Mean filter, BiShrink, and the proposed method, respectively. As shown in Figs. 6 and 7 , the distribution map of absolute difference demonstrates the results especially well.
The image that VisuShrink produces is too smooth while the SureShrink and MF methods leave more noise remaining. BiShrink and our method are considerably superior to VisuShrink and SureShrink. However, BiShrink produces more disturbing artifacts than our method. Our method also retains more textural details of the Lena and Boat images.
A performance comparison of the proposed method in terms of PSNR with other reported methods is shown in Table 1. These methods were tested using different noise levels of 15, 20, 25, and 30. It is obvious that the proposed method almost consistently produces the highest PSNRs for the two test images in all noise levels.
For instance, when the test image Lena was corrupted by AWGN with a standard deviation σ =30, the PSNR values after being processed by VisuShrink, SureShrink, Mean filter, and BiShrink were 25.70 dB, 25.49 dB, 27.79 dB, and 27.92 dB, respectively. However, the PSNR value of the proposed method was 28.98 dB.
PPT Slide
Lager Image
Comparing the results of various denoising methods, for Lena corrupted by noise σ=30. (a) Noisy image, (b) VisuShrink, (c) SureShrink, (d) Mean filter, (e) BiShrink, (f) proposed method.
PPT Slide
Lager Image
Comparing the results of various denoising methods, for boat corrupted by noise σ=30. (a) Noisy image, (b) VisuShrink, (c) SureShrink, (d) Mean filter, (e) BiShrink, (f) proposed method.
VI. CONCLUSIONS
In this paper, an efficient wavelet based image denoising adaptive thresholding method is proposed. The proposed method selects the threshold adaptively based on scale levels so that it leads to improved results. We also proposed a modified thresholding function that considered the dependency between the parent coefficient and child coefficient.
Performance of the proposed algorithm is measured by applying the PSNR and the proposed method provides high PSNR results. Our method also shows the good quality of restoring images that are compared with the VisuShrink, SureShrink, BiShrink, and Mean filter. Furthermore, our proposed algorithm can retain edges and details, while eliminating most of the noise.
References
Nasri M , Nezamabadi-Pour H 2009 "Image denoising in the wavelet domain using a new adaptive thresholding function" Neurocomputing 72 (4-6) 1012 - 1025    DOI : 10.1016/j.neucom.2008.04.016
Portilla J , Strela V , Wainwright M. J , Simoncelli E. P 2003 "Image denoising using scale mixtures of Gaussians in the wavelet domain" IEEE Transactions on Image Processing 12 (11) 1338 - 1351    DOI : 10.1109/TIP.2003.818640
Luisier F , Blu T , Unser M 2007 "A new SURE approach to image denoising: inter scale or thonormal wavelet thresholding" IEEE Transactions on Image Processing 16 (3) 593 - 606    DOI : 10.1109/TIP.2007.891064
Gonzalez R. C , Woods R. E 2002 Digital image processing 2nd ed Upper Saddle River Prentice Hall
Donoho D. L , Johnstone I. M 1992 "Ideal spatial adaptation by wavelet shrinkage" Department of StatisticsTechnical Report Stanford University
Donoho D. L , Johnstone I. M 1995 "Adapting to unknown smoothness via wavelet shrinkage" Journal of the American Statistical Association 90 (432) 1200 - 1224
Donoho D. L 1995 "De-noising by soft-thresholding" IEEE Transactions on Information Theory 41 (3) 613 - 627    DOI : 10.1109/18.382009
Sendur L , Selesnick I. W 2002 "Bivariate shrinkage functions for wavelet-based denoising exploiting interscale dependency" IEEE Transactions on Signal Processing 50 (11) 2744 - 2756    DOI : 10.1109/TSP.2002.804091
Sendur L , Selesnick I. W 2002 "Bivariate shrinkage with local variance estimation" IEEE Signal Processing Letters 9 (12) 438 - 441    DOI : 10.1109/LSP.2002.806054
Gao Y 2011 "Image denoising using adaptive threshold method in wavelet domain" International Journal of Maritime Information and Communication Sciences 9 (6) 763 - 768
Gao Y 2011 "Restoration of images contaminated by mixed Gaussian and impulse noise using a complex method" International Journal of Maritime Information and Communication Sciences 9 (3) 336 - 340