[英]Converting latitude/longitude to UTM ~ zone projection
我在WGS84 Zone = 21 N下有一組緯度/經度坐標,我試圖重新投影為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
任何關於可能會導致錯誤的建議/指導將不勝感激。
謝謝。 阿南德
如果您的經緯儀說:
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
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.