Я пробовал следующее, ввод: данные широты / долготы затем я рассчитаю прямоугольник вокруг него, скажем, на 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 в качестве параметров? Те, что указаны выше?