rgdal包lat / long - > UTM

时间:2011-11-02 12:45:21

标签: r gis

有人可以看出这有什么问题:

library(rgdal)

# Make a two-column matrix, col1 = long, col2 = lat
xy <- cbind(c(-107), c(26))
# Convert it to UTM coordinates (in units of meters)
project(xy, "+proj=utm +zone=51 ellps=WGS84")

据我所知,纬度大于84且小于80是无效的,但为什么我会这样做:

In project(xy, "+proj=utm +zone=51 ellps=WGS84") :
  1 projected point(s) not finite

谢谢!

1 个答案:

答案 0 :(得分:11)

经度使用正数(107,而不是-107)。负面是西半球