Как создать многоугольники из точек с помощью пакетов R?

Question

У меня есть несколько наборов точек (для разных лет ~ 20)

Я хочу генерировать многоугольники для каждого набора точек, используя r пространственных пакетов.

Я знаю, что это может быть сделано с помощью ГИС, но, как я хочу что-то периодический процесс в R будет

полезно.

Answer 1

Вы не предоставили нам доступ к вашим данным, но вот пример для точек, представляющих города мира, используя подход, описанный Карсоном Фармером по телефону his blog. Надеюсь, это будет, чтобы вы начали ...

# Carson's Voronoi polygons function 
voronoipolygons <- function(x) { 
    require(deldir) 
    require(sp) 
    if (.hasSlot(x, 'coords')) { 
    crds <- x@coords 
    } else crds <- x 
    z <- deldir(crds[,1], crds[,2]) 
    w <- tile.list(z) 
    polys <- vector(mode='list', length=length(w)) 
    for (i in seq(along=polys)) { 
    pcrds <- cbind(w[[i]]$x, w[[i]]$y) 
    pcrds <- rbind(pcrds, pcrds[1,]) 
    polys[[i]] <- Polygons(list(Polygon(pcrds)), ID=as.character(i)) 
    } 
    SP <- SpatialPolygons(polys) 
    voronoi <- SpatialPolygonsDataFrame(SP, data=data.frame(x=crds[,1], 
    y=crds[,2], row.names=sapply(slot(SP, 'polygons'), 
    function(x) slot(x, 'ID')))) 
} 

Пример 1: Входной сигнал является SpatialPointsDataFrame:

# Read in a point shapefile to be converted to a Voronoi diagram 
library(rgdal) 
dsn <- system.file("vectors", package = "rgdal")[1] 
cities <- readOGR(dsn=dsn, layer="cities") 

v <- voronoipolygons(cities) 

plot(v) 

Пример 2: Ввод в векторы х, у координаты:

dat <- data.frame(x=runif(100), y=runif(100)) 
v2 <- voronoipolygons(dat) 
plot(v2) 

Answer 2

Тот же принцип, как показано на jbaums, но проще код:

library(dismo) 
library(rgdal) 
cities <- shapefile(file.path(system.file("vectors", package = "rgdal")[1], "cities")) 

v <- voronoi(cities) 
plot(v)