R reprojecting map lat/long

Question

我用mapproj软件包将lat/long设置为地图投影（在本例中为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