奇怪的坐标单位答案

【问题标题】：Strange coordinate units奇怪的坐标单位
【发布时间】：2018-06-26 15:58:03
【问题描述】：

我收到了一个需要对其进行分析的 .gdb 文件。这是一个 ESRI 文件，但我无权访问任何 ESRI 产品。我使用“rgdal”包在 R 中打开它，但它显示的坐标似乎非常偏离。

例如，(-9288065, 4604652) 应该映射到肯塔基州。我最初的想法是缺少小数，所以应该是 -92.88 和 -46.04，但这映射到肯塔基州以北的某个地方，这是错误的。

这些坐标单位是不是很眼熟？我如何将它们转换为纬度和经度？

【问题讨论】：

您可能会在专门针对 GIS 问题的 Stack Exchange 站点获得更好的答案：Geographic Information Systems。
那些坐标不是纬度/经度，但可能是某种 NAD83 州平面，我猜是肯塔基州。您需要知道它们是什么类型的投影，以及它们的单位是英尺还是米。只是试了几下，我无法接近，所以它们可能不是任何国家飞机。询问发送给您的人它在哪个投影中。
看看在你的空间对象上调用sp::proj4string。这将告诉您投影。我不认为它缺少小数；相反，它看起来可能是以英尺为单位的投影，可能会出现在 proj4 字符串中。您还可以通过 CRS 代码在spatialreference.org 获取有关预测的信息
@camille 结果如下： "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 + units=m +nadgrids=@null +no_defs"，我不知道怎么解释（我没有GIS背景）。
"+proj=merc" 表示它是墨卡托投影，"+units=m" 表示它以米为单位。查看有关更改投影的教程here。你知道你想要什么投影吗？

标签： r gis

【解决方案1】：

在@camille 的帮助下解决了。

对于未来的用户：

library(rgdal)
library(raster)

myFeatureClass<-readOGR('your_data_here')

#checks the projection type
sp::proj4string(myFeatureClass)

#convert projection to longlat 
final_data <- as.data.frame(spTransform(myFeatureClass, CRS("+proj=longlat +datum=WGS84")))

【讨论】：