WaveThresh Help

InvBasis.wst


Invert a (packet-ordered) non-decimated wavelet transform using a given basis.

DESCRIPTION

Invert a (packet-ordered) non-decimated wavelet transform using a given basis.

USAGE

InvBasis.wst(wst, nv)

REQUIRED ARGUMENTS

wst
An object of class wst that contains coefficients of a packet ordered non-decimated wavelet transform (e.g. produced by the wst function) that you wish to invert.
nv
An object of class nv (a node vector) which specifies the particular basis in wst that you wish to use.

OPTIONAL ARGUMENTS

None.

VALUE

A vector containing the inverted reconstruction.

SIDE EFFECTS

None

DETAILS

The packet-ordered non-decimated wavelet transform computed by wst computes the coefficients of an input vector with respect to a library of all shifts of wavelet basis functions at all scales. Here "all shifts" means all integral shifts with respect to the finest scale coefficients, and "all scales" means all dyadic scales from 0 (the coarsest) to J-1 (the finest) where 2^J = n where n is the number of data points of the input vector. As such the packet-ordered non-decimated wavelet transform contains a library of all possible shifted wavelet bases.

Basis selection It is possible to select a particular basis and invert the representation with respect to the selected basis. This is what InvBasis.wst does. In WaveThresh such a basis is selected by creating a nv (node.vector) class object which identifies the basis. Then the function InvBasis takes the wavelet representation and the node.vector and inverts the representation with respect to the selected basis. The two functions MaNoVe and numtonv create a node.vector: the first by using a Coifman-Wickerhauser minimum entropy best-basis algorithm and the second by basis index.

Basis averaging. Rather than select a basis it is often useful to preserve information from all of the bases. For example, in curve estimation, after thresholding a wavelet representation the coefficients are coefficients of an estimate of the truth with respect to all of the shifted basis functions. Rather than select one of them we can average over all estimates. This sometimes gives a better curve estimate and can, for example, get rid of Gibbs effects. See Coifman and Donoho (1995) for more information on how to do curve estimation using the packet ordered non-decimated wavelet transform, thresholding and basis averaging.

For an example of denoising using the packet-ordered non-decimated wavelet transform and basis averaging see Johnstone and Silverman, 1997. The WaveThresh implementation of the basis averaging algorithm is to be found in Nason and Silverman, 1995

RELEASE

Version 3.6.0 Copyright Guy Nason 1995

SEE ALSO

AvBasis, InvBasis, MaNoVe, nv object, numtonv, wst, wst object.

EXAMPLES

#
# See example in the help to MaNoVe.wst
# for an example of InvBasis.wst
#