本文介绍了rgdal软件包lat/long-> UTM的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
有人可以看到这是怎么回事:
Can someone see what's wrong with this:
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都是无效的,但为什么会出现:
I understand that Latitude greater than 84 and smaller than 80 are invalid but why do I get:
In project(xy, "+proj=utm +zone=51 ellps=WGS84") :
1 projected point(s) not finite
谢谢!
推荐答案
使用正数表示经度(107,而不是-107).对西半球是负面的
use positive number for longitude (107, instead of -107). negative is for western hemisphere
这篇关于rgdal软件包lat/long-> UTM的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!