R - rgdal 패키지로 UTM 좌표를 변환 할 때의 문제

Question

큰 UTM 좌표 (영역 33)를 lat-long (WGS84) 좌표로 변환하는 작업이 주어졌습니다. 위치는 노르웨이 국경 근처에 위치해야합니다. 나는 rgdal 패키지를 사용하여 R의 좌표를 변환하는 시도

#extracting some examples from data set 

lat <- c(7790281, 7726438, 7266202, 7259480, 7271802) 
long <- c(1053817, 1054025, 451754, 475228, 462235) 
df <- data.frame(lat, long) 

#converting from UTM33 to WGS84 

library(rgdal) 
tmp <- data.frame(coords.x = df$lat, coords.y = df$long) 
coordinates(tmp) <-c("coords.x","coords.y") 
proj4string(tmp) <- CRS("+proj=utm +zone=33 ellps=WGS84") #UTM zone 33 
CRS.new <- CRS("+init=epsg:4326") # WGS84 
coords <- spTransform(tmp,CRS.new) 

#making a dataframe out of coordinates 
coords <- data.frame(lat=coords@coords[,1], long=coords@coords[,2]) 

이 그러나

 lat  long 
1 69.19936 5.840162 
2 68.92219 5.857922 
3 66.55129 2.585804 
4 66.52160 2.721701 
5 66.58242 2.644578 

을 산출, 온라인 변환 도구 (http://www.rcn.montana.edu/resources/converter.aspx) 내가 얻을로 좌표를 변환 할 때 다음과 같은 서로 다른 결과

lat.corr long.corr 
1 69.63112 29.37223 
2 69.07562 29.00514 
3 65.51455 13.95675 
4 65.45688 14.46554 
5 65.56618 14.18178 

두지도를 동일한지도에 플로팅하면 rgdal로 변환 한 점이 잘못 배치 된 것을 분명히 알 수 있습니다. 만약

Map

은 문제가 무엇인지 이해하려고 노력하지만, 고려 및 연구의 몇 시간 후, 나는 아직도 그것을 찾을 수 없습니다. 문제가 무엇인지 찾아내는 데 도움이 될 것입니다!

편집 : UTM 영역 경계를 가로 지르는 Easting 좌표와 관련이있을 수 있습니다. rgdal이 문제를 해결할 수 있습니까?

Answer 1

나는 그것이 있어야

tmp <- data.frame(coords.x = df$lat, coords.y = df$long) 

당신에게 혼합 위도와 경도를 추측 :

tmp <- data.frame(coords.x = df$lon, coords.y = df$lat) 

그리고 당신은 얻을 :

> (coords <- data.frame(lon=coords@coords[,1], lat=coords@coords[,2])) 
    lon  lat 
1 29.37172 69.63114 
2 29.00472 69.07564 
3 13.95675 65.51455 
4 14.46554 65.45688 
5 14.18178 65.56618