我不能評論這個,所以我不得不在這裏發佈它作爲答案。但在我找到解決方案之前,輸出座標應該顛倒嗎?
解決方案:(改編自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
它的工作原理!!!!非常感謝。只有兩點意見:有必要加載庫rgdal。是的,座標是相反的 –
是的,redal包是必需的,事實上,一旦你加載redal,sp會自動加載。我已經編輯了上面的答案。 – crwang