如何用Python获取1月1日UT午夜以来的天数（用于太阳赤纬计算）

获取1月1日UT午夜以来的天数N #

要得到符合要求的N值，核心是计算当年1月1日UTC午夜到当前UTC时间的时间差，并转换为带小数的总天数。用Python的datetime库可以轻松实现：

import datetime
import math

def get_N():
    # 获取当前UTC时间
    now_utc = datetime.datetime.utcnow()
    # 构造当年1月1日UTC午夜的时间对象
    jan1_utc = datetime.datetime(now_utc.year, 1, 1, 0, 0, 0)
    # 计算时间差
    time_diff = now_utc - jan1_utc
    # 将时间差转换为总天数（包含小数，精确到时分秒）
    N = time_diff.total_seconds() / 86400  # 86400秒 = 1天
    return N

结合公式计算太阳赤纬 #

根据太阳赤纬计算公式，将N代入后即可得到赤纬值（单位：度）：

def calculate_solar_declination():
    N = get_N()
    # 将角度转换为弧度进行计算
    deg_to_rad = math.pi / 180
    rad_to_deg = 180 / math.pi
    
    # 公式各部分计算
    term1 = math.sin(-23.44 * deg_to_rad)
    term2 = math.cos( (360/365)*(N + 10) * deg_to_rad + (360/math.pi)*0.0167*math.sin( (360/365)*(N - 2) * deg_to_rad ) )
    declination_rad = math.asin(term1 * term2)
    # 转换回角度
    declination_deg = declination_rad * rad_to_deg
    return declination_deg

调用示例 #

if __name__ == "__main__":
    current_N = get_N()
    current_declination = calculate_solar_declination()
    print(f"当前N值: {current_N:.4f}")
    print(f"当前太阳赤纬: {current_declination:.2f}°")

关键说明 #

• 必须使用UTC时间：因为N的定义是基于UT（协调世界时）的，用本地时间会产生时区偏移误差。
• 小数天数的意义：total_seconds() / 86400会把小时、分钟、秒转换为天的小数部分，比如中午12点对应0.5天，满足公式对N的精度要求。

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