library(sgeostat)
data(maas)
library(tripack)
maas.tri <- tri.mesh(maas$x, maas$y, duplicate="remove")
plot(maas.tri)
points(maas$x, maas$y, pch="x")
axis(1)
axis(2)
maas.vor <- voronoi.mosaic(maas$x, maas$y, duplicate="remove")
plot(maas.vor)
points(maas$x, maas$y, pch="x")