有人能看出这有什么问题吗：
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
Thanks!
使用正数表示经度（107，而不是 -107）。西半球为负
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
