R将地图重新​​投影到纬度/经度

Question

我使用mapproj包将lat / long集转换为map投影（在本例中为lagrange）。接下来我用它来定义空间上的其他几个点。我现在想将这些新点转换回lat / long，但是mapproj不提供这样的功能，我也没有在其他地方找到它。

cord<-data.frame(x, y)

lagrange<-mapproject(cord$x, cord$y, projection="lagrange")
#Some  
#clustering
#here

#Here are cluster centers, but in the form of lagrange projection. I want a lat/long cord of them
cluster.centers

如何在R中完成？

Answer 1

我阅读了包文档，我不相信作者希望他的一组投影可以用于绘图之外。

在R中使用地理空间数据的“标准”方法是通过 sp 对象的集合。转换sp对象的方法是通过 spTransform 。对象和方法都可以在 sp 包中找到。应该有一些教程可供选择，如果你不熟悉GIS和R，这可能是一个陡峭的学习曲线。

本文详细介绍了将“长拉”数据转换为其他投影的示例：

Error when re-projecting spatial points using spTransform in rgdal R