densityplot-methods {rasterVis} | R Documentation |
Draw kernel density plots (with lattice
) of Raster objects.
## S4 method for signature 'RasterLayer,missing' densityplot(x, data=NULL, maxpixels = 1e+05, xlab='', ylab='', main='', col='black', ...) ## S4 method for signature 'RasterStackBrick,missing' densityplot(x, data=NULL, layers, FUN, maxpixels = 1e+05, xlab='', ylab='', main='', par.settings=rasterTheme(), ...) ## S4 method for signature 'formula,Raster' densityplot(x, data, dirXY, maxpixels = 1e+05, xscale.components=xscale.raster, yscale.components=yscale.raster, auto.key = list(space = 'right'), par.settings=rasterTheme(),...)
x |
A |
data |
|
layers |
A numeric or character which should indicate the layers to be displayed. |
dirXY |
A direction as a function of the coordinates (see
|
FUN |
A function to applied to the |
maxpixels |
A numeric, for |
xlab, ylab, main, col, xscale.components, yscale.components,
par.settings, auto.key |
Arguments for |
... |
Additional arguments for |
Oscar Perpiñán Lamigueiro
densityplot
,xscale.raster
,
yscale.raster
, rasterTheme
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)