Wechat from skyfield import eclipselib from skyfield.api import load, wgs84 from numpy import pi from zoneinfo import ZoneInfo bjt = ZoneInfo("Asia/Shanghai") ts = load.timescale() eph = load("de422.bsp") moon = eph["moon"] earth = eph["earth"] r2d = 180./pi beijing = earth + wgs84.latlon(39.9042, 116.4074, 44) lunar = ['半影月食','月偏食','月全食'] t0 = ts.utc(2026, 1, 1) t1 = ts.utc(2027, 1, 1) t, y, details = eclipselib.lunar_eclipses(t0, t1, eph) idx = 0 for ti, yi in zip(t, y): print( ti.astimezone(bjt).strftime("%Y-%m-%d %H:%M:%S"), f"{lunar[yi]}" ) print(f"\t日月最近夹角：\t{details['closest_approach_radians'][idx]*r2d:.2f}°") print(f"\t月球半径：\t{details['moon_radius_radians'][idx]*r2d:.2f}°") print(f"\t半影半径：\t{details['penumbra_radius_radians'][idx]*r2d:.2f}°") print(f"\t本影半径：\t{de...