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

Question

我使用mapproj包将经/纬度集转换为地图投影（在本例中为lagrange）。 接下来，我用它来定义空间上的其他几个点。 我现在想将这些新点转换回纬度/经度，但是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，那么这可能是一个艰难的学习过程。

这篇文章详细介绍了将“长纬度”数据转换为其他预测的示例：

在rgdal R中使用spTransform重新投影空间点时出错