satellite - PyEphem 中地球卫星的矩形坐标

Question

我已经搜索了大量的 PyEphem 网页，下载了代码，并调查了对象，但我没有找到一种方法来获得某个名义惯性坐标系中地球卫星的地心直角坐标。我错过了一些明显的东西吗？我正在寻找“传统”的 x、y、z 和 x-dot、y-dot、z-dot。非常感谢您的任何建议。

Answer 1

好问题!
libastro内部的卫星定位程序使用 x y z内部坐标，但 libastro在 Python 有机会看到它们之前将它们扔掉。我本可以尝试修补 libastro然后 PyEphem 使数据可用，但是在研究您的问题时，我发现最近做了大量工作来更新 SGP4 卫星跟踪算法并为其提供测试套件:

http://www.celestrak.com/publications/AIAA/2006-6753/

在这个发现的激励下，我做了一些更好的事情:我花了周末创建了一个新的纯 Python sgp4基于此引用实现的包!

http://pypi.python.org/pypi/sgp4

很快我将围绕这个卫星预测引擎构建一个名为 Skyfield 的完整天文库，但目前您应该能够直接从 Python 包索引(上面的链接)下载它并调用它原始以获得您的传统坐标:

from sgp4.earth_gravity import wgs72
from sgp4.io import twoline2rv

line1 = ('1 00005U 58002B   00179.78495062  '
         '.00000023  00000-0  28098-4 0  4753')
line2 = ('2 00005  34.2682 348.7242 1859667 '
         '331.7664  19.3264 10.82419157413667')

satellite = twoline2rv(line1, line2, wgs72)
position, velocity = satellite.propagate(
    2000, 6, 29, 12, 50, 19)

请尝试一下，让我知道它是否适用于您的机器!