所以我打算确定一组数百个地理坐标中的哪些点位于形状多边形中。我能够将所有坐标(纬度和经度)转换为SpatialPointsDataFrame,我也可以在R中读取和绘制shapefile。 我的shp文件和点的图表是:
从stackoverflow中的类似问题我试图在sp包中使用函数point.in.polygon以及over()函数。 多边形:polygon1 要点:spdf
pt.in.poly <- sp::over( spdf, polygon1, fn = NULL)
或
new_shape <- point.in.poly(spdf, polygon1)
但在这两种情况下,我得到的结果是:
Error: identicalCRS(x, y) is not TRUE
如果它是一个坐标系统,我该如何解决?或者我做错了什么?