rgdal 包纬度/经度 -> UTM

Question

有人可以看到这有什么问题吗:

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

谢谢!

Answer 1

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