gavo.stc.geo module

Coordinate systems for positions on earth.

class gavo.stc.geo.WGS84[source]

Bases: object

the WGS84 reference system.

GM = 398600500000000.0
J2 = 0.00108263
a = 6378137.0
f = 0.0033528106647474805
omega = 7.292115e-05
gavo.stc.geo.geocToGeod(long, phip, rho=1, refSys=<class 'gavo.stc.geo.WGS84'>)[source]

returns geodetic coordinates long, phi for geocentric coordinates.

refSys defaults is the reference system the geodetic coordinates are expressed in.

This will not work at the poles – patches welcome.

See Astron. Almanac, Appendix K; we go for the iterative solution discussed there.

gavo.stc.geo.geodToGeoc(long, phi, height, refSys=<class 'gavo.stc.geo.WGS84'>)[source]

returns geocentric coordinates lambda, phi’, rho for geodetic coordinates.

refSys defaults is the reference system the geodetic coordinates are expressed in.

height is in meter, long, phi in degrees.

See Astron. Almanac, Appendix K.