<< Chapter < Page Chapter >> Page >

Denoising algorithms based on the rwt

Soft-thresholding

In the traditional method of soft-thresholding, where the universal threshold is used, coefficients below a specified threshold are shrunk to zero while those above the threshold are shrunk by a factor of σ ^ ( 2 log ( N ) ) . On the orthogonal wavelet transforms, it has been shown to exhibit the following property:

Theorem 1 For a sequence of i.i.d. random variables z i N ( 0 , 1 ) , P ( max i = 1 . . N ( 2 log N ) ) 1 for N .

Bivariate shrinkage

Sendur and Selesnick proposed a bivariate shrinkage estimator by estimating the marginal variance of the wavelet coefficients via small neighborhoods as well as from the the corresponding neighborhoods of the parent coefficients. The developed method maintains the simplicity and intuition of soft-thresholding.

We can write

y k = w k + n k ,

where w k are the parent and child wavelet coefficients of the true, noise-free image and n k is the noise. We have for our variance, then, that

σ y 2 = σ k 2 + σ n 2 .

Noting that we will always be working with one coefficient at a time, we will suppress the k .

In , Sendur and Selesnick proposed a bivariate pdf for the wavelet coefficient w 1 and the parent w 2 to be

p w ( w ) = 3 2 π σ 2 exp ( - 3 σ ( w 1 2 + w 2 2 ) ,

where the marginal variance σ 2 is dependent upon the coefficient index k . They derived their MAP estimator to be

w ^ 1 = ( y 1 2 + y 2 2 - 3 σ n 2 σ ) + y 1 2 + y 2 2 y 1

To estimate the noise variance σ n 2 from the noisy wavelet coefficients, they used the median absolute deviance (MAD) estimator

σ ^ n 2 = m e d i a n ( | y i | ) 0 . 6745 , y i s u b b a n d H H ,

where the estimator uses the wavelet coeffiecients from the finest scale.

The marginal variance σ y 2 was estimated using neighborhoods around each wavelet coefficient as well as the corresponding neighborhood of the parent wavelet coefficient. For instance, for a 7x7 window, we take the neighborhood around y 1 , ( 4 , 4 ) to be the wavelet coefficients located in the square (1, 1), (1, 7), (7, 7), (7, 1) as well as the coefficients in the second level located in the same square; this square is denoted N ( k ) . The estimate used for σ y 2 is given by

σ ^ y 2 = 1 M y i N ( k ) y i 2 ,

where M is the size of the neighborhood N ( k ) . We can then estimate the standard deviation of the true wavelet coefficients through :

σ ^ = ( σ ^ y 2 - σ ^ n 2 ) + .

We then have the information we need to use equation .

Bls-gsm

Portilla, et. al. propose the BLS-GSM method for denoising digital images, which may be used with orthogonal and redundant wavelet transforms as well as with pyramidal schemes. They model neighborhoods of coefficients at adjacent positions and scales as the product of a Gaussian vector and a hidden positive scalar multiplier, so that the neighborhoods are defined similarly as in the BiShrink algorithm. The coefficient within each neighborhood around a reference coefficient of a subband are modeled with a Gaussian scale mixture (GSM) model. The chosen prior distribution is the Jeffrey's prior, p z ( z ) 1 z .

They assume the image has additive white Gaussian noise, although the algorithm also allows for nonwhite Gaussian noise. For a vector y corresponding to a neighborhood of N observed coefficients, we have

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, The art of the pfug. OpenStax CNX. Jun 05, 2013 Download for free at http://cnx.org/content/col10523/1.34
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'The art of the pfug' conversation and receive update notifications?

Ask