使用Skyfield计算地金最近/最远距离时触发ValueError问题求助

核心错误点 #

1. 未定义的position变量：distance_at函数里直接调用position.radec()，但position既不是函数参数，也没有在函数内计算金星与地球的相对位置，这是触发错误的主要原因。
2. 重复加载星历文件：eph = load('de421.bsp')放在while循环内，每次循环都重新加载，既浪费资源又没必要。
3. 函数逻辑与需求不匹配：当前只调用了find_maxima找最远点，但需求是同时找最近（最小值）和最远（最大值）位置。
4. 冗余变量：earth_d, planet_d定义后未使用，属于无效代码。
5. 潜在的空结果问题：如果某一年没有找到极值点，zip(times, distance)会因空序列触发错误，需要先判断结果是否存在。

修复后的完整代码 #

from skyfield.api import load
from skyfield.searchlib import find_maxima, find_minima

# 只加载一次星历文件和时间尺度，提升效率
ts = load.timescale()
eph = load('de421.bsp')
sun, earth, venus = eph['sun'], eph['earth'], eph['venus']

def distance_at(t):
    # 计算地球与金星的相对位置，获取距离
    position = venus.at(t) - earth.at(t)
    return position.distance().au  # 返回天文单位的距离

# 设置搜索步长
distance_at.step_days = 15.0

# 遍历2000到2023年
for year in range(2000, 2024):
    t0 = ts.utc(year)
    t1 = ts.utc(year + 1)
    
    print(f"\n=== {year}年 金星与地球距离极值 ===")
    
    # 找最远点（最大值）
    max_times, max_distances = find_maxima(t0, t1, distance_at)
    if max_times:
        print("最远位置：")
        for t, dist in zip(max_times, max_distances):
            print(f"  {t.utc_strftime('%Y-%m-%d %H:%M:%S')}，距离：{dist:.4f} AU")
    else:
        print("  未找到最远位置")
    
    # 找最近点（最小值）
    min_times, min_distances = find_minima(t0, t1, distance_at)
    if min_times:
        print("最近位置：")
        for t, dist in zip(min_times, min_distances):
            print(f"  {t.utc_strftime('%Y-%m-%d %H:%M:%S')}，距离：{dist:.4f} AU")
    else:
        print("  未找到最近位置")

关键修复说明 #

• 正确计算相对位置：在distance_at函数内，通过venus.at(t) - earth.at(t)得到两者的相对位置向量，再调用distance().au获取以天文单位为单位的距离。
• 优化资源加载：将星历文件和时间尺度的加载移到循环外，避免重复IO操作。
• 覆盖需求场景：同时调用find_maxima和find_minima分别获取最远和最近位置。
• 增加空结果判断：在遍历结果前先检查times是否非空，避免zip处理空序列时触发错误。
• 代码结构优化：用for循环替代while循环，更符合Python风格，逻辑更清晰。

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