histogram-methods: Histogram of Raster objects.

In rasterVis: Visualization methods for the raster package

Description

Draw histograms (with lattice) of Raster objects.

Usage

## S4 method for signature 'RasterLayer,missing'
histogram(x, data=NULL, maxpixels = 1e+05, nint=100, 
            xlab='', ylab='', main='', col='gray', ...)

## S4 method for signature 'RasterStackBrick,missing'
histogram(x, data=NULL, layers, FUN,
            maxpixels = 1e+05, nint=100,
            xlab='', ylab='', main='', col='gray',
            between=list(x=0.5, y=0.2),
            as.table=TRUE,
            xscale.components=xscale.raster,
            yscale.components=yscale.raster,
            scales=list(x=list(relation='free'),
              y=list(relation='free',
                draw=FALSE)),
            strip.names=c(TRUE, TRUE),
            par.settings=rasterTheme(),
            ...)

## S4 method for signature 'formula,Raster'
histogram(x, data, dirXY,
            maxpixels = 1e+05,
            xscale.components=xscale.raster,
            yscale.components=yscale.raster,
            strip=TRUE,
            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.
nint	Number of breaks for the histogram. See the documentation of histogram at lattice for details.
maxpixels	A numeric, for sampleRandom.
xlab, ylab, main, col	Arguments for histogram.
between, as.table, xscale.components, yscale.components, scales, strip, strip.names, par.settings	Graphical parameters of lattice. See xyplot for details.
...	Additional arguments for histogram

Author(s)

Oscar Perpiñán Lamigueiro

See Also

histogram,xscale.raster, yscale.raster, rasterTheme

Examples

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

## 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)
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

histogram(SISmm)
histogram(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')

histogram(~asinh(pop)|landClass, data=s,
            scales=list(relation='free'),
            strip=strip.custom(strip.levels=TRUE))

## End(Not run)

rasterVis documentation built on May 2, 2019, 6:49 p.m.