[英]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.