search

ros-geographic-info / geographic_info

ROS packages for geographic information
http://ros.org/wiki/geographic_info
66 stars 62 forks source link

Bad UTM Transformation #24

Closed FletcherFT closed 3 years ago

FletcherFT commented 6 years ago

Comparing the geodesy and pygeodesy modules:

from pygeodesy import utm as pyutm
from geodesy import utm
latlon=(-41.0,147.0)
pyutm.toUtm(latlon[0],latlon[1]) # [Z:55, H:S, E:500000, N:5461243]
utm.fromLatLong(latlon[0],latlon[1]).toPoint() # x: 500000.0 y: -4538757.06191 z: 0.0

Easting seems to be fine, Northing is the problem:

latlon=(-40.0,146.0)
pyutm.toUtm(latlon[0],latlon[1]) # [Z:55, H:S, E:414640, N:5571764]
utm.fromLatLong(latlon[0],latlon[1]).toPoint() # x: 414639.538157 y: -4428236.06463 z: 0.0

Confirmed that the pygeodesy module is giving correct UTMs according to geoplaner. Am I missing an important aspect of the geodesy.py namespace?