R中多边形上的覆盖点

Question

我需要在 100,000 个多边形上叠加 40,000 个点。多边形是使用以下方法创建的：

polygon <- gBuffer(pc, width=500, byid=TRUE)

叠加层是使用以下方法创建的：

test <- over(pts, polygon)

因此，我有一个数据框，其中包含落在每个多边形（均值）中的点的属性。我的问题是是否可以将多边形的 gps 坐标（质心）cbind() 到我的数据框？

我查看了 google 和 stackoverflow，但没有找到关于该主题的任何内容。解决方案可能很简单，但我看不到 :)

任何帮助将不胜感激！

谢谢 马丁

编辑： 那么如果我想选择多边形呢？ 有（一个或多个）点。我用过

test <- over(polygons, pts, returnList = TRUE)

如何将多边形的坐标绑定到数据框列表（每个对应一个多边形）。我需要信息。作为列表，因为我需要计算行驶距离。 谢谢 M.

Answer 1

在您的示例中，test 不是data.frame，而是具有pts 中每个点的多边形索引的向量。试试

pts$polygon_id = over(pts, polygon)

这会将pts 从SpatialPoints 对象提升为SpatialPointsDataFrame。