从R中的lat long列表创建多边形

我有一个包含Id、lat和long的超过500k行的列表，如下所示：

Id  Latitude    Longitude   
1   7896    -50.33766   -22.23764   
2   7896    -50.33768   -22.23767   
3   7896    -50.33768   -22.23768   
4   7896    -50.33770   -22.23775

我需要根据Id创建一个多边形几何体，每个几何体生成一个不同的多边形

如何使用此数据创建SpatialPolygonDataFrame

有人能帮我吗

---

谢谢。

考虑到数据的

Id

列中有几个唯一的值，您可以对每个值应用一个函数来创建

空间多边形数据帧的列表。然后，您可以使用
do.call
查找列表。但是，这也应该适用于一个唯一的Id

library(sp)

# recreating your data

data = structure(list(Id = c(7896, 7896, 7896, 7896), 
                      Latitude = c(-50.33766, -50.33768, -50.33768, -50.3377),
                      Longitude = c(-22.23764, -22.23767, -22.23768, -22.23775)), 
                 class = "data.frame", row.names = c(NA, -4L))

# Applying a function to each unique value in the Id column of data

poly_list <- sapply(unique(data$Id), function(x){
  # create a polygon for each unique Id
  poly = Polygon(data[data$Id == x, 3:2]) # or Polygon(data[data$Id == x, c('Longitude', 'Latitude')])
  # create a spatial polygon from the polygon
  polys = SpatialPolygons(list(Polygons(list(poly), ID = x)))
  # convert the spatial polygons to spatial polygons dataframe
  as(polys, 'SpatialPolygonsDataFrame')
})

# rbind the list of spatial polygons dataframe

poly_list <- do.call(rbind, poly_list)

# visualize 

plot(poly_list)

库（sp）
#重新创建数据
数据=结构（列表Id=c（7896789678967896），
纬度=c（-50.33766，-50.33768，-50.33768，-50.3377），
经度=c（-22.23764，-22.23767，-22.23768，-22.23775）），
class=“data.frame”，row.names=c（NA，-4L））
#将函数应用于数据Id列中的每个唯一值
多边形列表