2015-07-21 50 views
1

我需要從轉換Gauss-Krüger座標爲WGS84座標系統(由Google使用系統)。我使用R來做到這一點,但我找不到解決此問題的簡單示例。Gauss-Krüger座標爲WGS84 in R

你能幫我嗎?

地址::HAUPTSTRASSE 62 70736費爾巴赫,德國

輸入(高斯 - 克呂格):3519358.50 5411371.00

輸出(WGS84):48.839580,9.262591

我的用於測試所述轉化數據

謝謝!

回答

3

我不能評論這個,所以我不得不在這裏發佈它作爲答案。但在我找到解決方案之前,輸出座標應該顛倒嗎?

解決方案:(改編自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 
+1

它的工作原理!!!!非常感謝。只有兩點意見:有必要加載庫rgdal。是的,座標是相反的 –

+0

是的,redal包是必需的,事實上,一旦你加載redal,sp會自動加載。我已經編輯了上面的答案。 – crwang