gstat (version 2.0-3)

vgmArea: point-point, point-area or area-area semivariance

Description

Compute point-point, point-area or area-area variogram values from point model

Usage

vgmArea(x, y = x, vgm, ndiscr = 16, verbose = FALSE, covariance = TRUE)

Arguments

x

object of class SpatialPoints or SpatialPolygons

y

object of class SpatialPoints or SpatialPolygons

vgm

variogram model, see vgm

ndiscr

number of points to discretize an area, using spsample

verbose

give progress bar

covariance

logical; compute covariances, rather than semivariances?

Value

semivariance or covariance matrix of dimension length(x) x lenght(y)

Examples

# NOT RUN {
library(sp)
demo(meuse, ask = FALSE, echo = FALSE)
vgmArea(meuse[1:5,], vgm = vgm(1, "Exp", 1000)) # point-point
vgmArea(meuse[1:5,], meuse.area, vgm = vgm(1, "Exp", 1000)) # point-area
# }