r - 使用 R 将 lat/long 点的数据框空间连接到多边形 shapefil

Question

我希望编写一个 R 程序来读取纬度/经度点的数据帧和 13 个多边形的 shapefile，然后识别每个纬度/经度点位于哪个多边形内。是否有一个 R 包可以让我进行这种类型的空间连接？我一直在使用 rgdal 来读取多边形 shapefile，但我没有看到用于与我的数据框中的 lat/long 列进行空间连接的选项。谢谢！

Answer 1

G. Cocca 的解决方案经过一些小的修改后效果很好。

R最终，解决方案需要使用'slibrary(sp)和的四个步骤library(EcoSpatial)：

1. 将 lat/long 点数据帧转换为SpatialPointsDataFrameusing sp.
2. 使用sp'sproj4string函数重新投影，SpatialPointsDataFrame以便它使用与我的多边形 shapefile 相同的投影。
3. 使用EcoSpatial'spoints.in.poly函数在空间上连接两个数据集。这是关键步骤，因为它创建了一个新的SpatialPointsDataframe，将相关的 polgyon 数据附加到每个 xy 点。
4. 用来as.data.frame强制SpatialPointsDataframe后面变成正则R data.frame。

再次感谢您的帮助！