问题描述
这是来自物理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)
解决方法
在GitHub / poliastro上存在一些与poliastro如何处理参考框架有关的未解决问题。 polipoli似乎正在尝试使用地心天参考系统(GCRS)作为其以地球为中心的惯性参考系。通过查看poliastro master分支中的代码,以及查看与参考框架相关的许多未解决问题,我将它们称为伪GCRS框架。
根据定义,GCRS是相对论正确的参照系。 Poliastro完全是牛顿式的。它忽略了光速问题。它忽略了地球时间(TT)和重心动态时间(TDB)之间的差异。它忽略了测地线进动。它忽略了很多事情。结合以上内容,便是我撰写poliastro将伪GCRS框架用于以地球为中心的惯性框架的原因。