pyproj为UTM返回inf

Question

我正在获取x= inf inf inf.....和y= inf inf inf...以前有没有人遇到过这个错误？我正在尝试将lon和lat转换为x，y值。我有以下代码，但我找不出哪些是不正确的。这是圣地亚哥的数据集。

import scipy.io;
from pyproj import Proj;

canyon = scipy.io.loadmat('INIT.MAT');
topo = canyon.get('siocanyon');

lat = tuple(topo[:, 0])
lon = tuple(topo[:, 1])
z = topo[:, 2]

myProj = Proj(proj='utm', zone=11, ellps='WGS84', preserve_units=False)
y, x = myProj(lon, lat)

Answer 1

我刚刚用圣地亚哥的一个示例位置测试了代码：

from pyproj import Proj

lat = 32.715687
lon = -117.161380
#z = 100

myProj = Proj(proj='utm', zone=11, ellps='WGS84', preserve_units=False)
x, y = myProj(lon, lat)
print(x, y)

给出

484876.3791751535 3619780.206742558

您可以使用online converter测试您的结果。