pyproj возвращает информацию для UTM

Я получаю x = [inf inf inf.....] и y = [inf inf inf...]. Кто-нибудь сталкивался с этой ошибкой раньше? Я пытаюсь преобразовать долготу и широту в значения 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)

python latitude-longitude utm pyproj
person Sushmitha Kudari    schedule 02.03.2020    source источник
comment
Можете ли вы дать ссылку на веб-сайт, где я могу скачать набор данных и прочитать об этом? Без него трудно тестировать.   -  person Rusty Robot    schedule 02.03.2020
comment
drive.google.com/open?id=1V9b-39NNh3wS52-ol0DjzA8RWsM7jfPP   -  person Sushmitha Kudari    schedule 02.03.2020

Ответы (1)


arrow_upward
0
arrow_downward

Я только что протестировал код с образцом позиции из Сан-Диего:

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

Вы можете проверить свои результаты с помощью онлайн-конвертера.

person Ionut Ticus    schedule 22.03.2020
comment
Вы пробовали фрагмент выше? Одним из возможных объяснений вашей ошибки может быть то, что вы поменяли местами lat с lon при чтении. - person Ionut Ticus; 27.03.2020