This function is used to threshold the coefficients by group (or block).
BlockThresh(wc, j0, hatsigma, L, qmf, thresh = "hard")
wavelet coefficients.
coarsest decomposition scale.
estimator of noise variance.
Block size (n mod L must be 0).
Orthonormal quadrature mirror filter.
'hard' or 'soft'.
wcb
wavelet coefficient estimators.