Отношение широты / долготы к utm к широте / долготе - это очень некорректно, почему?

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

Теперь я преобразовываю его в широту / долготу и с помощью скрипта:

http: / /robotics.ai.uiuc.edu/~hyoon24/LatLongUTMconversion.py Я получаю результат, которого просто не может быть, lon до этого около 7, потом около 2.

zone, easting, northing = LLtoUTM(23, location.get_lat(), location.get_lon()) 

topUTM = northing + error
bottomUTM = northing - error
leftUTM = easting - error
rightUTM = easting + error
left, top = UTMtoLL(23, leftUTM, topUTM, zone)

Это ошибка в моем коде, или может быть скрипт ошибочен?

Итак, я попытался использовать pyproj, просто от широты / долготы до utm до широты / долготы, чтобы увидеть, что произойдет

>>> p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
>>> p

>>> x,y = p(47.9941214, 7.8509671)
>>> print x,y
5159550.36822 1114087.43925
>>> print p(x,y,inverse=True)
(47.971558538495991, 7.8546573140162605)

И здесь это не так далеко, как в сценарии сверху, но он все еще кажется достаточно неправильным, чтобы его нельзя было использовать. Почему? Что я могу сделать, чтобы получить более точные результаты?

РЕДАКТИРОВАТЬ:

Я запустил test (), и все тесты прошли успешно.

в файле epsg такого нет. Ближайшее, что я нашел, было следующее:

<32632> +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs <>

no tmerc. Также что мне нужно для передачи Towgs84 в качестве параметров? Те, что указаны выше?

11
задан Kara 18 March 2014 в 19:44
поделиться