地理空间数据：形状文件的奇怪经度/纬度答案

【问题标题】：Geospatial data: weird longitude/latitude to the shapefile地理空间数据：形状文件的奇怪经度/纬度
【发布时间】：2016-09-15 21:14:01
【问题描述】：

我刚刚开始使用地理空间数据，并且有一个与经度/纬度相关的问题。我下载了法语部门的 shapefile (http://professionnels.ign.fr/geofla) 并使用 R 读取数据：

library(rgdal)
library(maptools)
map <- readOGR(dsn="...",layer="DEPARTEMENT")
map_f <- fortify(map)
head(map_f)

如果我查看经纬度，它是类似 long = 885661.8、lat = 6679942 的数字。相比之下，巴黎的经纬度组合是 2.352222、48.85661 (geocode("Paris"))。所以我真的不明白为什么这些数字的大小如此不同？有人可以帮助我或提供一些好的参考吗？

提前致谢！

【问题讨论】：

欢迎来到！请提供一个可重现的示例，您提供的链接是法语网站，为什么不将实际链接发布到 shapefile...
感谢您的回复。这个链接应该可以做到：wxs-telechargement.ign.fr/oikr5jryiph0iwhw36053ptm/…。在此先感谢:)
无论您发布什么链接都要求我接受一些我不愿意做的认证/验证，请尝试将 shapefile 保存到保管箱并发布保管箱链接……总之，我怀疑你会通过提供粗略的链接获得有关 SO 的帮助...
@CyrusMohammadian 是正确的，我们确实要求在 R 标签中提供可重复的示例。话虽如此，他们通常不喜欢保管箱链接，因为它们经常断开，但我想有总比没有好。我没有收到带有该链接的认证请求，但它确实给了我一个压缩文件，这些文件存在很高的恶意软件风险。
纬度/经度并不奇怪，它们只是简单的投影。您能否提供map_f 的输出，以便我们可以看到分配给数据的投影。您可能只需要重新投影即可获得熟悉的纬度/经度。一个很好的参考是the difference between projected and geographic coordinate systems。

【解决方案1】：

坐标看起来是投影的（可能是 UTM），而不是经纬度。我认为这些字段只是贴错了标签。我不懂法语，但在快速浏览您链接的网站时，这引起了我的注意：

大都会：(RGF 93) 投影 Lambert93

En outre-mer : (系统法律)- Projections UTM

两者都不是地理坐标系。请追踪元数据并检查坐标系。您还可以检查proj4string(map) 以查看是否定义了坐标系。

【讨论】：

您好，该命令给了我： "+proj=utm +zone=20 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0" 。我不确定所有这些意味着什么，所以我想我有一些阅读要做
您的数据位于投影坐标系中：通用横轴墨卡托、20 区和 WGS84 基准面。坐标单位以米为单位，表示东向和北向。
试试spTransform(map, CRS("+proj=longlat +datum=WGS84"))