0

我想使用 Python 从 celestrak.com 的 TLE 文件中提取和解析单个卫星的数据。我使用“readtle”找到了使用 pyephem 执行此操作的代码，但这在 skyfield 中不起作用。如何从 skyfield 的 tle 文件中提取数据？

pyephem 基本上已被弃用。因此，对于 pyephem/skyfield 的新手，我选择使用 skyfield。

我正在尝试从 celestrak.com 的文件中自动提取 TLE 信息。我发现几年前的代码是在 python 2.xx 下使用 pyephem 完成的所以，我将代码转换为 Python 3.7.x。

但是，用于从 TLE 文件中提取单个条目的代码的 readtle 方面在 skyfield 上不起作用。我正在从 celestrak.com 服务器访问文件。

我也尝试安装 pyephem，但出现错误。而且，无论如何，我都想留在skyfield。

我找到的代码如下：

class Satellite:
    def __init__(self, name, l1, l2):
        self.e        = ephem.readtle(name, l1, l2)
    def compute(self):
        self.e.compute(datetime.datetime.utcnow())
        self.long   = math.degrees(float(self.e.sublong))
        self.lat    = math.degrees(float(self.e.sublat))
        self.height = abs(int(self.e.elevation))
        r = 6378150 + self.height
        self.r = r / 1000000.0
        self.x      = -cos(radians(self.lat)) * cos(radians(self.long)) * self.r
        self.y      = sin(radians(self.lat)) * self.r
        self.z      = cos(radians(self.lat)) * sin(radians(self.long)) * self.r
        self.vlists  = sat_shape
        self.label = pyglet.text.Label(self.e.name, y=15, anchor_x="center", color=(255,255,255,200))

    def draw(self):
        glLoadIdentity()
        glTranslatef(0,0,-zoom+6.37815)
        glRotatef(ro, 1, 0, 0)
        glRotatef(rocc, 0, 0, 1)
        glTranslatef(0,0,-6.37815)
        glRotatef(-angle_y, 1, 0, 0)
        glRotatef(angle_x, 0, 1, 0)
        glTranslatef(self.x,self.y,self.z)
        glColor3f(1,0,0)
        glScalef(zoom/100.0, zoom/100.0, zoom/100.0)
        for v in self.vlists:
            v.draw(GL_TRIANGLE_STRIP)
        glScalef(0.02, 0.02, 0.02)
        glRotatef(-angle_x, 0, 1, 0)
        glRotatef(angle_y, 1, 0, 0)
        glRotatef(-rocc, 0, 0, 1)
        glRotatef(-ro, 1, 0, 0)
        self.label.draw()
        self.draw_line()

有没有人有任何想法？

4

1 回答 1

2

是的，您可以使用 loaderSkyfield直接从 Celestrak 加载 TLE 文件：

from skyfield.api import Topos, load
    
stations_url = 'http://celestrak.com/NORAD/elements/stations.txt'
satellites = load.tle(stations_url)
satellite = satellites['ISS (ZARYA)']
print(satellite)

EarthSatellite 'ISS (ZARYA)' number=25544 epoch=2014-01-20T22:23:04Z

文档中提供了更多信息。

于 2019-07-24T18:27:38.177 回答