The Python poliastro
package documentation states that, when propagating from position and velocity vectors, the orbit's reference frame will be "one pseudo-inertial frame around the attractor," which in the case of earth-orbiting objects, means some form of ECI. However, I want to make sure I am matching up coordinates because I am using another bit of code to convert a position on Earth to ECI using the IERS.
I have already checked answers to a similar post and found the clarity to be lacking. I want a specific answer in this case from someone who has used poliastro
in the following form:
from astropy import units
from poliastro.bodies import Earth
from poliastro.twobody import Orbit
from poliastro.twobody.propagation import kepler
ephemerisPosition = (ephemerisPosition * units.earthRad).to(units.meter)
ephemerisVelocity = (ephemerisVelocity * units.earthRad / units.day).to(units.meter / units/second)
seconds = ephememerisStart + secondsOffset
for i, time in enumerate(seconds):
ss = Orbit.from_vectors(Earth, ephemerisPosition, ephemerisVelocity).propagate(time * units.second, method=kepler)
pymap3d
which relies onastropy
for its coordinate conversions, probably the same waypoliastro
does: github.com/geospace-code/pymap3d/pull/35#issue-477303148 $\endgroup$poliastro
mentionsastropy
explicitly in its documentation: docs.poliastro.space/en/stable/api/safe/twobody/… $\endgroup$