將經度/緯度轉換爲UTM〜區域投影

我有一組緯度/長度座標在WGS84 Zone = 21 N下，我試圖重新投影到UTM格式，但我得到的數字真的很奇怪;將經度/緯度轉換爲UTM〜區域投影

library('rdgal') 
library('sp') 
plot_coord_SW <- data.frame(PLOT=c(1,2,3,8,9,10),east_proj=c(5.01964,5.022444,5.025106,5.016922,5.014864,5.011755),north_proj=c(-58.6252,-58.6247, -58.6248,-58.6203,-58.6203,-58.6211)) 

coordinates(plot_coord_SW) <- c("east_proj","north_proj") #converts to spatialpointsdataframe 
proj4string(plot_coord_SW) <- CRS("+proj=longlat +datum=WGS84") #sets the 
coordinate reference system=WGS84 
reproj_UTM <- spTransform(plot_coord_SW, CRS("+proj=utm +zone=21 ellps=WGS84")) 
      coordinates PLOT 
1 (3670898, -8215608) 1 
2 (3671045, -8215729) 2 
3 (3671124, -8215883) 3 
4 (3671325, -8215133) 8 
5 (3671255, -8215020) 9 
6 (3671066, -8214901) 10

我穿過與第一組座標[積1]與網上轉換的檢查，這是完全不同的，正確的基礎上，地圖投影：http://www.rcn.montana.edu/resources/converter.aspx

 coordinates PLOT

1（319820.9，555058.9） 1

任何意見/指導什麼錯誤可能會不勝感激。

謝謝。阿南德

來源

2017-08-25 Anand Roopsind

您的經緯度長是錯誤的，如果說

north_proj=c(-58.6252,-58.6247, -58.6248,-58.6203,-58.6203,-58.6211))

這意味着南; 如果通過east_proj你得到座標像你提到的那些交換north_proj：

> coordinates(reproj_UTM_32) 
    east_proj north_proj 
[1,] 319820.9 555058.9 
[2,] 319877.2 555368.8 
[3,] 319866.8 555663.2 
[4,] 320363.6 554757.0 
[5,] 320363.0 554529.4 
[6,] 320273.5 554185.8

來源

2017-08-25 22:34:34

啊....感謝埃利奧！這很簡單，非常感激 –

將經度/緯度轉換爲UTM〜區域投影

回答

相關問題