This function extracts and returns a vector of mother wavelet coefficients, corresponding to a particular resolution level, from a wd wavelet decomposition object.
The pyramid of coefficients in a wavelet decomposition (returned from the wd function, say) are packed into a single vector in WaveThresh.
# S3 method for wd
accessD(wd, level, boundary=FALSE, aspect="Identity", ...)A vector containing the mother wavelet coefficients at the required resolution level (the coefficients might have been modified depending on the value of the aspect argument).
Wavelet decomposition object from which you wish to extract the mother wavelet coefficients.
The resolution level at which you wish to extract coefficients.
some methods of wavelet transform computation handle the boundaries by keeping some extra bookkeeping coefficients at either end of a resolution level. If this argument is TRUE then these bookkeeping coefficients are returned when the mother wavelets are returned. Otherwise, if FALSE, these coefficients are not returned.
The aspect argument permits the user to supply a function to modify the returned coefficients. The function is applied to the vector of coefficients before it is returned. This can be useful, say, with the complex DWT where you could supply aspect="Mod" if you wanted to return the modulus of the coefficients at a given resolution level. The default argument, "Identity", ensures that the coefficients are not modified before returning.
any other arguments
Version 3.5.3 Copyright Guy Nason 1994
G P Nason
The need for this function is a consequence of the pyramidal structure of Mallat's algorithm and the memory efficiency gain achieved by storing the pyramid as a linear vector. AccessD obtains information about where the smoothed data appears from the fl.dbase component of an wd object, in particular the array
fl.dbase$first.last.d
which gives a complete specification of index numbers and offsets for
wd.object$D.
Note that this function is a method for the generic function accessD.
Note also that this function only retrieves information from wd class objects. To insert coefficients into wd objects you have to use the putD function (or more precisely, the putD.wd method).
Mallat, S. G. (1989) A theory for multiresolution signal decomposition: the wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 11, 674--693.
Nason, G. P. and Silverman, B. W. (1994). The discrete wavelet transform in S. Journal of Computational and Graphical Statistics, 3, 163--191
wr, wd, accessD, filter.select, threshold
#
# Get the 4th resolution level of wavelet coefficients.
#
dat <- rnorm(128)
accessD(wd(dat), level=4)
Run the code above in your browser using DataLab