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测试结果。