如何使用Skyfield查询行星的地心黄道位置及相关运动参数

环境配置 #

• 安装依赖库：pip install skyfield
• 示例使用de421星历文件，首次运行会自动下载到本地，无需手动处理

完整实现代码 #

from skyfield.api import load
from skyfield.framelib import ecliptic_J2000

# 加载时间工具与星历文件
ts = load.timescale()
eph = load('de421.bsp')

# 定义观测时间（示例为2024年1月1日0时UTC，可按需修改）
t = ts.utc(2024, 1, 1, 0, 0, 0)

# 定义地心观测站与目标行星（示例目标为火星，可替换为其他行星对应key：'mercury','venus','mars','jupiter_barycenter','saturn_barycenter','uranus_barycenter','neptune_barycenter'）
earth = eph['earth']
target_planet = eph['mars']

# 计算目标行星相对地心的位置
astrometric = earth.at(t).observe(target_planet)

# 转换为J2000黄道坐标系下的坐标
lat, lon, distance = astrometric.frame_latlon(ecliptic_J2000)
# 获取黄道坐标系下的速度矢量
velocity = astrometric.velocity.frame_xyz(ecliptic_J2000)

# 输出四类参数
print("黄道经度：", lon.degrees, "度")
print("黄道纬度：", lat.degrees, "度")
print("地心距离：", distance.au, "天文单位 /", distance.km, "公里")
print("黄道坐标系速度分量（AU/天）：", velocity.AU_per_d)

参数说明 #

• 黄道经度/纬度：基于J2000平黄道坐标系，返回结果为带单位的角度对象，调用.degrees可直接获取十进制角度数值
• 地心距离：返回观测时刻行星与地心的直线距离，支持直接调用.au取天文单位值、.km取公里值
• 运行速度：返回黄道坐标系下的三维速度矢量，默认单位为AU/天，可按需求自行转换为km/s等单位（1AU/天 ≈ 1731.46 km/s）

内容的提问来源于stack exchange，提问作者Butterfly Fiziks