11

I would like to be more specific about what I want to do. I get coordinates in ECI and I need to get the latitude and longitude from this. How can I do? I was searching but I could'nt find anything about it. Thanks again.

(I'm doing a small program in java that shows the position of a satellite in a given time. So, I used the NORAD SGP algorithm, and I have the position (x,y,z) and velocity(Vx,Vy,Vz). But the coordinates system used by this algorithm is the ECI, according what I read. Now I need to draw the satellite in a map, but I can't convert this coordinates to some system that could help me. I think if I can convert it to longitude and latitude it would be easy to draw. Could you help me? how can I do it? What is the best option(UTM,etc)? Thanks.)

voodoo14
  • 515
  • 2
  • 6
  • 10

2 Answers2

9

What you want to do is called ECI/ECEF (cartessian) to Geodetic (lat/lon) conversion. This conversion is the most complex of all the geodetic conversions as the closed form solution is complicated. See page 34 of Stevens and Lewis, Aircraft Control and Simulation for a discussion of the coordinate systems: http://books.google.com/books/about/Aircraft_control_and_simulation.html?id=T0Ux6av4btIC

ECI to geodetic is a two step process:

The first step is the easiest in that you need to convert ECI (earth centered inertial) to ECEF (earth centered/earth fixed).

The second step is to convert ECEF to geodetic. You can read about solving this via Newton-Ralphson here: http://en.wikipedia.org/wiki/Geodetic_system

However, if I remember correctly, Newton-Raphson becomes unstable around the poles. The closed form solutions are much more complicated. I have successfully implemented Zhu's method. The advantage of the closed form solution is no iterations and there are no singularities (technically there are singularities but not above the earth). The reference: J. Zhu. Conversion of earth-centered earth-fixed coordinates to geodetic coordinates. Technical Report IEEE Log NO. T-AES/30/3/1666, IEEE, December 1993.

TreyA
  • 3,339
  • 2
  • 20
  • 26
  • 1
    Any chance you can expand on converting ECI to ECEF? – MattSayar Mar 06 '14 at 00:56
  • 2
    ECI and ECEF share the same origin and z-axis. They differ by the Earth's rotation rate so the transformation is just a 2-D cos sin; -sin cos type transform. – TreyA Mar 06 '14 at 03:22
0

The PyMap3D library for Python has the eci2ecef and the ecef2geodetic functions. They interface with Astropy by default and are quite good, thanks to the hard work of Michael Hirsch (SciVision - he accepts donations).

You need, of course, to know the times at which the ECI coordinates were observed; without that you cannot convert them to ECEF. I recommend not using the simpler eci2geodetic function unless high precision and accuracy are not required; it will be somewhat quicker because it does not account for Earth's nutation, etc.).

brethvoice
  • 350
  • 1
  • 4
  • 14
  • P.S. If you want to geek out a bit, check out Michael's PhD dissertation: https://zenodo.org/record/3239516#.YTJ3BfdJF9B – brethvoice Sep 03 '21 at 19:35