我已经在 python 中使用 Skyfield 实现了 sgp4 算法。
我需要 1 天时间间隔的位置向量,周期为 1 秒。为此,我必须反复计算 sgp4,每次增加 1 秒。
有没有办法返回一个包含所有位置的向量?
最佳答案
幸运的是,Skyfield 是为 NumPy 数组的时间对象设计的,它将生成与数组中的对象一样多的位置。以下是描述该功能的文档部分:
http://rhodesmill.org/skyfield/time.html#date-arrays
这是一个完成您描述的任务的示例:
from numpy import arange
from skyfield.api import load, EarthSatellite
iss_tle0 = """\
1 25544U 98067A 18184.80969102 .00001614 00000-0 31745-4 0 9993
2 25544 51.6414 295.8524 0003435 262.6267 204.2868 15.54005638121106
"""
my_earthsat = EarthSatellite(*iss_tle0.splitlines())
tz_offset = -4
ts = load.timescale()
t = ts.utc(2018, 6, 21, tz_offset, 0, arange(24 * 60 * 60))
astrometrics = my_earthsat.at(t)
这段代码目前不是特别快,因为地球卫星计算尚未优化,但应该比在循环中手动构建多个时间对象更快。
关于python - Skyfield:在给定的时间间隔内以 1 秒的周期获得 sgp4 结果,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/49494082/