使用R将经纬度的数据框空间连接到多边形shapefil

Question

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

Answer 1

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

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

1. 使用sp将SpatialPointsDataFrame数据帧转换为SpatialPointsDataFrame 。
2. 使用sp的proj4string函数重新投影SpatialPointsDataFrame以便它使用与我的多边形shapefile相同的投影。
3. 使用EcoSpatial的points.in.poly函数在空间上连接两个数据集。 这是关键步骤，因为它创建了一个新的SpatialPointsDataframe ，该SpatialPointsDataframe将相关的polgyon数据附加到每个xy点。
4. 使用as.data.frame将SpatialPointsDataframe强制SpatialPointsDataframe回常规R data.frame 。

再次感谢您的帮助！