rtop (version 0.5-5)

rtopVariogram: create variogram for data with spatial support

Description

rtopVariogram will create binned variogram or cloud variogram of data with an areal support.

Usage

"rtopVariogram"(object, params = list(), ...) 
"rtopVariogram"(object, ...) 
"rtopVariogram"(object, formulaString, params=list(), cloud,  abins, dbins, ...) 
"rtopVariogram"(object, formulaString, params=list(), cloud,  abins, dbins, data.table = FALSE, ...)

Arguments

object

object of class rtop (see rtop-package) or a SpatialPolygonsDataFrame or SpatialPointsDataFrame with information about observations. If object is a SpatialPointsDataFrame, it must have a column with name area.

formulaString

formula that defines the dependent variable as a linear model of independent variables; suppose the dependent variable has name z, for ordinary and simple kriging use the formula z~1; for universal kriging, suppose z is linearly dependent on x and y, use the formula z~x+y. The formulaString defaults to "value~1" if value is a part of the data set. If not, the first column of the data set is used.

params

a set of parameters, used to modify the standard parameters for the rtop package, set in getRtopParams.

cloud

logical; if TRUE, calculate the semivariogram cloud, can be used to overrule the cloud parameter in params.

abins

possibility to set areal bins (not yet implemented)

dbins

possibility to set distance bins (not yet implemented)

data.table

an option to use data.table internally for the variogram computation for STSDF-objects

...

parameters to other functions called, e.g. gstat's variogram-function and to rtopVariogram.SpatialPointsDataFrame when the method is called with an object of a different class

Value

The function creates a variogram, either of type rtopVariogram or rtopVariogramCloud. This variogram is based on the variogram function from gstat, but has additional information about the spatial size or length of the observations. An rtop-object with the variogram added is returned if the function is called with an rtop-object as argument.For spatio-temporal objects (STSDF), the variogram is the spatially variogram, averaged for all time steps. There is a possibility to use data.table internally in this function, which can improve computation time for some cases.

See Also

rtop-package

Examples

## Not run: 
# library(rgdal)
# rpath = system.file("extdata",package="rtop")
# observations = readOGR(rpath,"observations")
# # Create a column with the specific runoff:
# observations$obs = observations$QSUMMER_OB/observations$AREASQKM
# 
# vario = rtopVariogram(observations, cloud = TRUE)
# ## End(Not run)