【发布时间】:2020-08-25 14:13:27
【问题描述】:
这是从物理 SE question 交叉发布的,以防那里的人们不经常/很多使用 poliastro。请在此处回答之前检查那里,以防它已经被事件克服。
Python poliastro 包 documentation 指出,当从位置和速度矢量传播时,轨道的参考系将是“围绕吸引子的一个伪惯性系”,这意味着当地球是 ECI 时的一种形式吸引子。但是,我想确保它的纪元与IERS 指定的纪元相同。
我已经检查了类似post 的答案,发现缺乏清晰度。在这种情况下,我希望以下列形式使用 poliastro 的人提供具体答案:
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提交了一个拉取请求,它依赖astropy进行ECI 转换,其方式可能与poliastro 相同:github.com/geospace-code/pymap3d/pull/35#issue-477303148 -
Poliastro 在其文档中明确提到了天体:docs.poliastro.space/en/stable/api/safe/twobody/…
标签: python-3.x physics coordinate-systems satellite orbital-mechanics