papyri
{ } Raw JSON

bundles / skimage 0.26.1rc0.dev0+git20260530.b607368ff / skimage / restoration / _denoise / denoise_wavelet

function

skimage.restoration._denoise:denoise_wavelet

source: /dev/scikit-image/src/skimage/restoration/_denoise.py :850

Signature

def   denoise_wavelet ( image ,  sigma = None ,  wavelet = db1 ,  mode = soft ,  wavelet_levels = None ,  convert2ycbcr = False ,  method = BayesShrink ,  rescale_sigma = True ,  * ,  channel_axis = None )

Summary

Perform wavelet denoising on an image.

Parameters

image : ndarray (M[, N[, ...P]][, C]) of ints, uints or floats

Input data to be denoised. image can be of any numeric type, but it is cast into an ndarray of floats for the computation of the denoised image.

sigma : float or list, optional

The noise standard deviation used when computing the wavelet detail coefficient threshold(s). When None (default), the noise standard deviation is estimated via the method in [2].

wavelet : str, optional

The type of wavelet to perform and can be any of the options pywt.wavelist outputs. The default is 'db1'. For example, wavelet can be any of {'db2', 'haar', 'sym9'} and many more.

mode : {'soft', 'hard'}, optional

An optional argument to choose the type of denoising performed. It noted that choosing soft thresholding given additive noise finds the best approximation of the original image.

wavelet_levels : int or None, optional

The number of wavelet decomposition levels to use. The default is three less than the maximum number of possible decomposition levels.

convert2ycbcr : bool, optional

If True and channel_axis is set, do the wavelet denoising in the YCbCr colorspace instead of the RGB color space. This typically results in better performance for RGB images.

method : {'BayesShrink', 'VisuShrink'}, optional

Thresholding method to be used. The currently supported methods are "BayesShrink" [1] and "VisuShrink" [2]. Defaults to "BayesShrink".

rescale_sigma : bool, optional

If False, no rescaling of the user-provided sigma will be performed. The default of True rescales sigma appropriately if the image is rescaled internally.

channel_axis : int or None, optional

If None, the image is assumed to be grayscale (single-channel). Otherwise, this parameter indicates which axis of the array corresponds to channels.

Returns

out : ndarray

Denoised image.

Notes

The wavelet domain is a sparse representation of the image, and can be thought of similarly to the frequency domain of the Fourier transform. Sparse representations have most values zero or near-zero and truly random noise is (usually) represented by many small values in the wavelet domain. Setting all values below some threshold to 0 reduces the noise in the image, but larger thresholds also decrease the detail present in the image.

If the input is 3D, this function performs wavelet denoising on each color plane separately.

Many wavelet coefficient thresholding approaches have been proposed. By default, denoise_wavelet applies BayesShrink, which is an adaptive thresholding method that computes separate thresholds for each wavelet sub-band as described in [1].

If method == "VisuShrink", a single "universal threshold" is applied to all wavelet detail coefficients as described in [2]. This threshold is designed to remove all Gaussian noise at a given sigma with high probability, but tends to produce images that appear overly smooth.

Although any of the wavelets from PyWavelets can be selected, the thresholding methods assume an orthogonal wavelet transform and may not choose the threshold appropriately for biorthogonal wavelets. Orthogonal wavelets are desirable because white noise in the input remains white noise in the subbands. Biorthogonal wavelets lead to colored noise in the subbands. Additionally, the orthogonal wavelets in PyWavelets are orthonormal so that noise variance in the subbands remains identical to the noise variance of the input. Example orthogonal wavelets are the Daubechies (e.g. 'db2') or symmlet (e.g. 'sym2') families.

Examples

from skimage import color, data
✓
img = img_as_float(data.astronaut())
img = color.rgb2gray(img)
⚠
rng = np.random.default_rng()
✓
img += 0.1 * rng.standard_normal(img.shape)
img = np.clip(img, 0, 1)
denoised_img = denoise_wavelet(img, sigma=0.1, rescale_sigma=True)
⚠

Aliases

  • skimage.restoration.denoise_wavelet