如何用Skyfield计算新月首次可见时间以重构古日耳曼历法？

1. Skyfield能否按「首次可见新月」定义计算？ #

Skyfield的almanac模块默认计算的是天文新月（日月黄经差为0的时刻），没有内置的「首次可见新月（蛾眉月）」计算功能，但可以基于Skyfield的天体位置计算能力自行扩展实现。

要计算首次可见新月，需要结合以下关键参数的判断（可通过Skyfield获取数据后推导）：

• 日月角距离：需要达到足够阈值（通常8°-12°，古日耳曼观测可能取10°左右），确保月牙有足够亮度
• 月球地平高度：观测时月球需在地平线以上，且高度不能过低（建议≥3°，避免大气消光导致无法观测）
• 太阳地平高度：太阳需处于地平线以下足够深度（比如≤-6°，对应民用暮光结束），避免日光掩盖月牙
• 观测位置：古日耳曼地区大致纬度范围（比如北纬50°-60°），固定经纬度可标准化计算结果

2. 「首次可见新月」的定义是否足够明确可计算？ #

这个定义并非绝对精确，确实会受观测者视力、天气等因素影响，但古日耳曼历法的观测是基于特定区域（中欧/北欧）的传统实践，因此可以通过设定标准化观测参数来得到可重复、与参考书籍对齐的计算结果：

• 固定观测点的经纬度（比如取古日耳曼核心区域的坐标）
• 统一设定角距离、高度、太阳位置的阈值
• 忽略极端天气等偶然因素，取理想观测条件下的结果

参考书籍《The Lunisolar Calendar of the Germanic Peoples》中的数据，正是基于这类标准化的传统观测规则，你可以调整上述阈值，让计算出的首次可见新月时间晚于冬至后12天，从而修正2034、2060、2061年的闰年判断结果。

简单实现思路示例 #

用Skyfield实现的核心步骤伪代码：

from skyfield.api import load, wgs84

# 加载天体数据
eph = load('de421.bsp')
sun, moon, earth = eph['sun'], eph['moon'], eph['earth']

# 设定古日耳曼观测点（示例：北纬55°，东经10°）
observer = wgs84.latlon(55.0, 10.0)

# 遍历冬至后12天内的时间点，计算参数
for t in time_range:
    astrometric = observer.at(t).observe(moon).apparent()
    sun_astrometric = observer.at(t).observe(sun).apparent()
    
    # 计算日月角距离
    angular_separation = astrometric.separation_from(sun_astrometric).degrees
    
    # 计算月球高度（含大气折射）
    moon_alt, _, _ = astrometric.altaz()
    moon_alt_deg = moon_alt.degrees
    
    # 计算太阳高度
    sun_alt, _, _ = sun_astrometric.altaz()
    sun_alt_deg = sun_alt.degrees
    
    # 判断是否满足首次可见条件
    if angular_separation >= 10 and moon_alt_deg >= 3 and sun_alt_deg <= -6:
        first_visible_new_moon = t
        break

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