使用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 。

再次感謝您的幫助！