rasterVis (version 0.47)

Formula methods: Formula methods

Description

Formula methods

Usage

# S4 method for formula,Raster
xyplot(x, data, dirXY, maxpixels=1e5,
            alpha=0.05,
            xscale.components=xscale.raster, yscale.components=yscale.raster,
            par.settings=rasterTheme(),...)
# S4 method for formula,Raster
hexbinplot(x, data, dirXY,
            xscale.components=xscale.raster, yscale.components=yscale.raster,
            par.settings=rasterTheme(),...)

Arguments

x

A formula describing the variables to be related. It may include the layer names (which are internally converted to valid ones with make.names) and the x, y variables representing the coordinates of the Raster object. Besides, if dirXY is not missing, the variable dirXY can also be included in the formula.

data

A Raster object.

dirXY

A direction as a function of the coordinates (see xyLayer).

maxpixels

A numeric, for sampleRegular.

alpha

A numeric, transparency of the points.

xscale.components, yscale.components,par.settings

Customization of lattice. See xyplot for details.

Additional arguments for the xyplot and hexbinplot functions.

Examples

Run this code
# NOT RUN {
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
names(r)

xyplot(test~y, data=r, alpha=0.5)

# }
# NOT RUN {
##Solar irradiation data from CMSAF
##Data available from http://www.box.net/shared/rl51y1t9sldxk54ogd44

old <- getwd()
##change to your folder...
setwd('CMSAF')
listFich <- dir(pattern='2008')
stackSIS <- stack(listFich)SISmm <- SISmm*24 ##from irradiance (W/m2) to irradiation Wh/m2
setwd(old)

names(SISmm) <- month.abb

##Relation between the January & February versus July radiation for four
##differents longitude regions.
xyplot(Jan+Feb~Jul|cut(x, 4), data=SISmm, auto.key=list(space='right'))
##Faster with hexbinplot
hexbinplot(Jan~Jul|cut(x, 6), data=SISmm)
# }

Run the code above in your browser using DataCamp Workspace