geospace-code / pymap3d

pure-Python (Numpy optional) 3D coordinate conversions for geospace ecef enu eci
https://geospace-code.github.io/pymap3d
BSD 2-Clause "Simplified" License
391 stars 85 forks source link

geodetic2enu produce different values from ecef2enu #49

Closed quonn77 closed 2 years ago

quonn77 commented 2 years ago

Invoking both geodetic2enu and ecef2enu with the same location and observable produce different output

import pymap3d as pm

lat1 = 45.977 lon1 = 7.658 h1 = 4531

lat0 = 46.017 lon0 = 7.750 h0 = 1673

x,y,z = pm.geodetic2ecef(lat1,lon1,h1)
e,n,u = pm.ecef2enu(x,y,z,lat0,lon0,h0)
e1,n1,u1 = pm.geodetic2enu(lat1,lon1,h1,lat0,lon0,h0)
print ("E=%f N=%f U=%f" % (e,n,u))
print ("E1=%f N1=%f U1=%f" % (e1,n1,u1))

Produce different output. I am expecting same result for enu the U output differs.

quonn77 commented 2 years ago

Same here...it was my mistake, I was printing wrong result because I have in the middle a geodetic2aer.....sorry again