plot.Wimage {fields}R Documentation

Plots 2-d wavelet coefficents by level and type


Produces a panel of images using the split.screen tools that organize the wavelet coefficients from a multiresolution by location, resolution level and type.


plot.Wimage(x, cut.min, graphics.reset = TRUE, common.range = FALSE, 
    color.table = tim.colors(128), Nlevel = NULL, with.lines = FALSE, 
    omd.width = 0.2, ...)


x matrix of coefficients
cut.min defines minimum level
graphics.reset If TRUE will reset device to orignal settings including closing the split.screen mode.
common.range If TRUE image plots will be on common color scale.
color.table Color table to be used for image plots. Default is Tim Hoar's favorite.
Nlevel Number of levels to plot. Default is to plot all.
with.lines If TRUE will add white outlines of pixels. Default is FALSE to prevent from tiny pixels from turning white!
omd.width Fration of device surface width devoted to the vertical color legend strips. Default usually accomodates most axis labels.
... Other Graphical parameters to be passed to the par function.


As with most complicated graphical figures you basically get what it draws although it should be easy to modify this function for customization. The split.screen set of graphical functions are used to divide up the plotting real estate into subplots of different sizes. The user can experiment with different outer margin space (omd.width) Use cex.axis argument in the call to change the size of the numerals in the color strip.

By setting grahics.reset to FALSE the function returns a matrix giving the split screen ids to reference each of the individual plots. Use the screen function to move to given plot and use a high level plotting function to overlay information. (See example below.) Use close.screen( all=TRUE) to turn off split.screen mode for subsequent and normal plotting.


Doug Nychka


old.par<- par(no.readonly=TRUE) # these functions may leave the device in 
# with some funny defaults

#multiresolution of John Lennon 
Wtransform.image( lennon, cut.min=16)->look 
plot.Wimage( look, cut.min=16)

# adding information 
plot.Wimage( look, cut.min=16, Nlevel=3, graphics.reset=FALSE)-> plot.layout
# plot.layout here is a 4X3 matrix with the screen numbers

# move to the smooth coefficients plot

 screen( plot.layout[1,1])
 plot( c(.5,16.5), c( .5,16.5), type="n",axes=FALSE) 
 points( 8,8, cex=2, pch="+")

# NOTE: just points( 8,8) will not work here. This has to do with split.screen not 
# reseting the plotting pars correctly. 

 box(col=6, lwd=2)# just for fun 

close.screen( all=TRUE)
par( old.par) # reset  to old settings 

[Package fields version 3.3.1 Index]