plotMODWT2D: Plot Two-Dimensional Wavelet Decomposition on "sampSurf"...
In ssWavelets: Wavelet Functionality for Package 'sampSurf'

Description Usage Arguments Details Value Author(s) See Also Examples

This routine will plot any of the individual “ssMODWT” decompositions below (waveType) at any level either as the raw wavelet coefficients or various wavelet variances (see Details for some restrictions). Covariance objects, “ssCovMODWT”, are also supported.

plotMODWT2D(ssMODWT, 
            level = 1, 
            waveType = c("ISO", "LH", "HL", "HH", "LL"), 
            decompType = c("modwt", "mra"), 
            type = c("raw", "var"), 
            isoSmooth = TRUE, 
            showPlot = TRUE, 
            showIZs = TRUE, 
            addLattice = FALSE, 
            col = NA, 
            boxCol = "gray", 
            title = NA, 
            ...
           )

`ssMODWT`	An object of class “`ssMODWT`” for raw or variance display, or an object of class “`ssCovMODWT`” for covariance display.
`level`	The decomposition level desired, j = 1,…,J0. Note that for `waveType = 'ISO'` and `type = 'var'`, it is legal to specify `j = 0`; please see the details below for an explanation.
`waveType`	‘ISO’ is for isotropic; the next three are the anisotropic decompositions: ‘LH’ horizontal; ‘HL’ vertical; ‘HH’ diagonal. Finally, the ‘LL’ is the smooth at the highest scale (level J0). Please see the vignette for more details.
`decompType`	Either “MODWT” or “MRA” decomposition may be displayed for `type = 'raw'`. Only the former is available for variance or covariance displays.
`type`	‘raw’ corresponds to the raw wavelet coefficients. ‘var’ is for a plot of the variance.
`isoSmooth`	`TRUE`: include ‘LLJ’ in the result for isotropic at level J (i.e., J0); `FALSE`: do not include the smooth component.
`showPlot`	`TRUE`: display the plot; `FALSE`: no display.
`showIZs`	`TRUE`: display the inclusion zones for the “Stem” objects; `FALSE`: no display.
`addLattice`	`TRUE`: create a lattice version of the graph (see details); `FALSE`: no lattice version.
`col`	A vector palette of colors for the surface; if this is `NA`, the default, then `palMODWT` is used.
`boxCol`	Outline color for perimeter box.
`title`	Add a title to the plot if desired. The default (`NA`) is a title (and subtitles) constructed from the arguments passed. Otherwise, one can pass their own title or `title =` “” for no title.
`...`	Gobbled.

This routine has a fair bit of functionality and it is the basis for the plotLevel2D function as well. The vignette has a number of examples that can be reviewed for more information.

Note that not all combinations of argument options are available. For example, variance is calculated only for decompType = “MODWT” as noted above. Also, level = 0 is only meaningful at this point for waveType = “ISO” and type = “var”, which presupposes decompType = “MODWT”. Also, only type = “var” is allowed when the ssMODWT argument passed is an object of class “ssCovMODWT”. This latter restriction is because there the covariance is composed of two “ssMODWT” objects; therefore, type = ‘raw’ would be ambiguous. One can always plot the raw wavelets for the individual “ssMODWT” objects that were used to create the covariance object.

With addLattice = TRUE, the actual graphical object is created using levelplot in the rasterVis package. Therefore, the latter must be available to use this option. The “lattice” graph is returned and may be plotted as usual.

A list with...

`r:`	The final “`RasterLayer`” object.
`plt:`	The “lattice” plot object if desired; `NULL` otherwise.

Jeffrey H. Gove

plotLevel2D

#
# 1. creates a sampSurf object with horizontal point sampling
# 2. creates a J_0 = 3-level MODWT decomposition object
# 3. a basic plot...
#
tr = Tract(c(x = 64, y = 64), cellSize = 1) #square tract ~0.5ha
btr = bufferedTract(10, tr)
ag3m = angleGauge(3)        #metric BAF
sshps = sampSurf(10, btr, iZone = 'horizontalPointIZ', angleGauge = ag3m,
        topDiam = c(0,0), startSeed = 123)
modwt.hps = ssMODWT(sshps, J = 3)
## Not run: 
r1 = plotMODWT2D(modwt.hps, waveType = 'ISO', type = 'var', level = 1)

## End(Not run)