3

我需要在 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。

4

1 回答 1

4

在您的示例中,test不是data.frame一个向量,而是一个向量,其中的每个点都有多边形索引pts。尝试

pts$polygon_id = over(pts, polygon)

pts从一个SpatialPoints对象提升为一个SpatialPointsDataFrame

于 2015-09-12T11:14:04.437 回答