WaveThresh
Help
put.wd3D
Put wavelet coefficient array into
a 3D wavelet object
DESCRIPTION
This function put an array of wavelet coefficients,
corresponding to a particular resolution level into a
wd wavelet decomposition object.
The pyramid of coefficients in a wavelet decomposition (returned from the
wd3D function, say) are packed into a single array in
WaveThresh3.
USAGE
putD.wd3D(x, v)
REQUIRED ARGUMENTS
- x
- 3D Wavelet decomposition object into which you
wish to insert the wavelet coefficients.
-
v
- This argument is a list with the following components:
a
-
A 3-dimensional array with each dimension of length equal to
two to the power of
lev which is the level at which
you wish to insert the coefficients into x.
lev
- The level at which you wish to insert the
coefficients into
x.
block
- A character string indicating which coefficient block you
wish to insert the coefficients into. This can be one of
GGG, GGH, GHG, GHH, HGG, HGH, HHG. Additionally this can be
HHH when the
lev argument above is zero.
OPTIONAL ARGUMENTS
None.
VALUE
A new wd3D object is returned with the coefficients
at level lev in block given by block are replaced
by the contents of a, if a is of the correct
dimensions!
DETAILS
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 an
array.
Note that this function is a method for the
generic function putD.
RELEASE
Version 3.9.6 Copyright Guy Nason 1997
SEE ALSO
accessD,
accessD.wd3D,
print.wd3D,
putD,
putDwd3Dcheck,
summary.wd3D,
threshold.wd3D,
wd3D,
wd3D object,
wr3D.
EXAMPLES
#
# Generate some test data
#
a <- array(rnorm(8*8*8), dim=c(8,8,8))
#
# Perform the 3D DWT
#
awd3D <- wd3D(a)
#
# Replace the second level coefficients by uniform random variables
# in block GGG (for some reason)
#
#
newsubarray <- list(a = array(runif(4*4*4), dim=c(4,4,4)), lev=2, block="GGG")
awd3D <- putD(awd3D, v=newsubarray)