使用R将纬度/长点的数据帧空间连接到多边形shapefil

Question

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

Answer 1

-G。经过一些小的改动后，Cocca的解决方案很有效。

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

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

再次感谢您的帮助！