| smooth.wavelet {biwavelet} | R Documentation |
Smooth wavelet in both the time and scale domains
Description
The time smoothing uses a filter given by the absolute value of the wavelet function at each scale, normalized to have a total weight of unity, which is a Gaussian function for the Morlet wavelet. The scale smoothing is done with a boxcar function of width 0.6, which corresponds to the decorrelation scale of the Morlet wavelet.
Usage
smooth.wavelet(wave, dt, dj, scale)
Arguments
wave |
wavelet coefficients |
dt |
size of time steps |
dj |
number of octaves per scale |
scale |
wavelet scales |
Value
Returns the smoothed wavelet.
Note
This function is used internally for computing wavelet coherence. It is only appropriate for the morlet wavelet.
Author(s)
Tarik C. Gouhier (tarik.gouhier@gmail.com)
Code based on WTC MATLAB package written by Aslak Grinsted.
References
Torrence, C., and P. J. Webster. 1998. The annual cycle of persistence in the El Nino/Southern Oscillation. Quarterly Journal of the Royal Meteorological Society 124:1985-2004.
Examples
# Not run: smooth.wt1 <- smooth.wavelet(wave, dt, dj, scale)