rasterVis (version 0.40)

densityplot-methods: Density plots for Raster objects.

Description

Draw kernel density plots (with lattice) of Raster objects.

Usage

"densityplot"(x, data=NULL, maxpixels = 1e+05, xlab='', ylab='', main='', col='black', ...)

"densityplot"(x, data=NULL, layers, FUN, maxpixels = 1e+05, xlab='', ylab='', main='', par.settings=rasterTheme(), ...)
"densityplot"(x, data, dirXY, maxpixels = 1e+05, xscale.components=xscale.raster, yscale.components=yscale.raster, auto.key = list(space = 'right'),  par.settings=rasterTheme(),...)

Arguments

x

A Raster* object or a formula.

data

NULL or a Raster object.

layers

A numeric or character which should indicate the layers to be displayed.

dirXY

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

FUN

A function to applied to the z slot of a RasterStackBrick object. The result of this function is used as the grouping variable of the plot.

maxpixels

A numeric, for sampleRandom.

xlab, ylab, main, col, xscale.components, yscale.components, par.settings, auto.key

Arguments for densityplot.

...

Additional arguments for densityplot

See Also

densityplot,xscale.raster, yscale.raster, rasterTheme

Examples

f <- system.file("external/test.grd", package="raster")
r <- raster(f)
densityplot(r)
s <- stack(r, r+500, r-500)
densityplot(s)

## Not run: 
# 
# old <- getwd()
# ##change to your folder...
# setwd('CMSAF')
# listFich <- dir(pattern='2008')
# stackSIS <- stack(listFich)
# stackSIS <- stackSIS*24 ##from irradiance (W/m2) to irradiation Wh/m2
# setwd(old)
# 
# idx <- seq(as.Date('2008-01-15'), as.Date('2008-12-15'), 'month')
# 
# SISmm <- setZ(stackSIS, idx)
# names(SISmm) <- month.abb
# 
# densityplot(SISmm)
# ##FUN applies to z if not NULL
# densityplot(SISmm, FUN=as.yearqtr)
# ## End(Not run)
## Not run: 
# ##http://neo.sci.gsfc.nasa.gov/Search.html?group=64
# pop <- raster('875430rgb-167772161.0.FLOAT.TIFF')
# pop[pop==99999] <- NA
# levelplot(pop, zscaleLog=10, par.settings=BTCTheme,
#           panel=panel.levelplot.raster, interpolate=TRUE)
# 
# ##http://neo.sci.gsfc.nasa.gov/Search.html?group=20
# landClass <- raster('241243rgb-167772161.0.TIFF')
# landClass[landClass==254] <- NA
# 
# 
# s <- stack(pop, landClass)
# names(s) <- c('pop', 'landClass')
# 
# densityplot(~asinh(pop)|landClass, data=s,
#             scales=list(relation='free'),
#             strip=strip.custom(strip.levels=TRUE))
# ## End(Not run)