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

In gstat: Spatial and Spatio-Temporal Geostatistical Modelling, Prediction and Simulation

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)

Author(s)

Edzer Pebesma

Examples

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

gstat documentation built on Sept. 11, 2024, 7:46 p.m.