Gauss-Krüger在R中与WGS84坐标

时间:2015-07-21 14:33:29

标签: r coordinate-systems

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

你能帮我吗?

我测试转换的数据:

地址:Hauptstraße62 70736费尔巴赫,德国

输入(Gauss-Krüger):3519358.50 5411371.00

产出(WGS84):48.839580,9.266591

感谢!!!

1 个答案:

答案 0 :(得分:3)

我不能对此发表评论所以我不得不在此发布它作为答案。但在我找到解决方案之前,输出坐标应该反转吗?

解决方案:(改编自here

From: <sip:alice@127.0.0.11:5060>;tag=dqixl8jjv0w

输出:

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