terra (version 1.9-46)

voronoi: Voronoi diagram and Delaunay triangles

Description

Get a Voronoi diagram or Delaunay triangles for points, or the nodes of lines or polygons

Usage

# S4 method for SpatVector
voronoi(x, bnd=NULL, tolerance=0, as.lines=FALSE, deldir=FALSE) 
# S4 method for SpatVector
delaunay(x, tolerance=0, as.lines=FALSE, constrained=FALSE)

Value

SpatVector

Arguments

x

SpatVector

bnd

SpatVector to set the outer boundary of the voronoi diagram

tolerance

numeric >= 0. Snapping tolerance applied to the input vertices before computing the diagram (0 is no snapping). Expressed in the units of the CRS of x, i.e. meters or feet for projected data and degrees for lon/lat data.

as.lines

logical. If TRUE, lines are returned without the outer boundary

constrained

logical. If TRUE, a constrained delaunay triangulation is returned

deldir

logical. If TRUE, the deldir is used instead of the GEOS C++ library method. It has been reported that deldir does not choke on very large data sets

Examples

wkt <- c("MULTIPOLYGON ( ((40 40, 20 45, 45 30, 40 40)), 
  ((20 35, 10 30, 10 10, 30 5, 45 20, 20 35),(30 20, 20 15, 20 25, 30 20)))", 
  "POLYGON ((0 -5, 10 0, 10 -10, 0 -5))")
x <- vect(wkt)
v <- voronoi(x)
v

d <- delaunay(x)
d

plot(v, lwd=2, col=rainbow(15))
lines(x, col="gray", lwd=2)
points(x)