Gauss-Krüger 坐标到 R 中的 WGS84答案

【问题标题】：Gauss-Krüger coordinates to WGS84 in RGauss-Krüger 坐标到 R 中的 WGS84
【发布时间】：2015-07-21 14:33:28
【问题描述】：

我需要将 Gauss-Krüger 坐标转换为 WGS84 坐标（Google 系统使用）。我正在使用 R 来执行此操作，但我找不到解决此问题的简单示例。

你能帮帮我吗？

我用于测试转换的数据：

地址：Hauptstraße 62 70736 费尔巴赫，德国

输入（高斯-克鲁格）：3519358.50 5411371.00

输出 (WGS84)：48.839580、9.262591

谢谢！！！

【问题讨论】：

标签： r coordinate-systems

【解决方案1】：

我无法对此发表评论，因此我不得不在此处将其发布为答案。但在我得到解决方案之前，输出坐标是否应该反转？

解决方案：（改编自here）

library(rgdal)
# Creating data
GK <- data.frame(cbind("X_GK"=3519358.50,"Y_GK"=5411371.00))

#Spatial Information, Coordinates Transform
coordinates(GK) <- c("X_GK", "Y_GK")
proj4string(GK) <- CRS("+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs") # Defining Gauss Krüger
GK.WGS84 <- spTransform(GK, CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")) # tranforming to WGS84 longlat

GK.WGS84

输出：

> GK.WGS84
SpatialPoints:
         X_GK     Y_GK
[1,] 9.262686 48.83949
Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

【讨论】：

它有效！！！！非常感谢。只有两个 cmets：需要加载库 rgdal。是的，坐标被颠倒了
是的，redal包是必须的，其实一旦加载redal，sp就会自动加载。我已经编辑了上面的答案。