Как получить широту / долготу из x,y декартовой координаты tmerc проекции?

Я пытаюсь преобразовать данные из декартовых координат OpenDrive в значения широты / долготы. Для pt = [6.714150516498e+05, 5.434880530093e+06] на немецком шоссе A9, я хочу знать, что такое координаты широты / долготы.

Я пытался решить это таким образом, но я получаю широту / долготу координат пустыни в Нигерии.

from pyproj import Proj, transform
inProj = Proj("+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs")
outProj = Proj(init='epsg:4326', preserve_units=True)

w = [6.714150516498e+05, 5.434880530093e+06]

lat, lon = transform(inProj, outProj, w[0], w[1])
print(lat, lon)

1 ответ

Решение

преобразование возвращается (x,y), Это означает: x = lon, y = lat,

>>> print transform(inProj, outProj, w[0], w[1])[::-1]
(49.04294637744738, 11.345684678187824)

Который где-то в Германии рядом с трассой А9.

Другие вопросы по тегам