使用python计算太阳高度角

爱吃小石榴16

于 2024-04-26 14:45:00 发布

阅读量1.4k

点赞数 20

文章标签：学习 python 开发语言网络

本文链接：https://blog.csdn.net/2401_84618514/article/details/138211325

版权

太阳高度角：对于地球上的某个地点，太阳高度角是指太阳光的入射方向和地平面之间的夹角，专业上讲太阳高度角是指某地太阳光线与通过该地与地心相连的地表切面的夹角。太阳高度角简称高度角。当太阳高度角为90°时，此时太阳辐射强度最大；当太阳斜射地面时，太阳辐射强度就小。（摘自百度百科）

具体相关定义、资料和图片都比较多，这里就不一一提了，下面直接说计算方法。

根据《制定地方大气污染物排放标准的技术方法 GB-T3840-91》规范中附录B提到的方法来计算太阳高度角h0。

具体计算公式摘录如下：

太阳高度角计算截图

1 计算太阳倾角

需输入的变量只有1个：一年中的日期序数，取值范围从0到364。

我在实际计算中计算的天数比较多，直接使用pandas的df.index.dayofyear函数计算得到的，另外一个计算思路是使用月份天数进行累加。根据上面的式B1即可求出太阳倾角。

在三角函数计算中要将角度转为弧度，直接使用math工具包中的radians即可，具体代码如下：

from math import sin, asin, cos, radians, pi, degrees

f0 = 360*Dn/365
f1 = 0.006918
f2 = 0.399912*cos(radians(f0))
f3 = 0.070257*sin(radians(f0))
f4 = 0.006758*cos(radians(f0*2))
f5 = 0.000907*sin(radians(f0*2))
f6 = 0.002697*cos(radians(f0*3))
f7 = 0.001480*sin(radians(f0*3))
f = (f1-f2+f3-f4+f5-f6+f7)*180/pi

其中有个问题，如果是闰年366天的话，日期序数是否是到365？

2 计算太阳高度角h0

需输入的变量包括：经度、纬度、小时数、太阳倾角。

其中这里的小时数即0~23（北京时间），其中仍需要注意的是弧度的转化。

结合上面的代码，最终计算太阳高度角如下：

from math import sin, asin, cos, radians, pi, degrees

def get_h0(lon, lat, Dn, Hour):
    """计算太阳高度角
    Arguments:
        lon {float} -- [经度(度)]
        lat {float} -- [纬度(度)]
        Dn {int} -- [一年中日期序数，0,1,2,3......，364]
        Hour {int} -- [小时(北京时间)]
    Returns:
        float -- [太阳高度角(度)]
    """
    # 第一步：计算太阳倾角
    f0 = 360*Dn/365
    f1 = 0.006918
    f2 = 0.399912*cos(radians(f0))
    f3 = 0.070257*sin(radians(f0))
    f4 = 0.006758*cos(radians(f0*2))
    f5 = 0.000907*sin(radians(f0*2))
    f6 = 0.002697*cos(radians(f0*3))
    f7 = 0.001480*sin(radians(f0*3))
    f = (f1-f2+f3-f4+f5-f6+f7)*180/pi

    # 第二步：计算太阳高度角
    p3 = radians(15*Hour + lon - 300)
    lat = radians(lat)
    f = radians(f)  # 转为弧度，切记！
    h0 = asin(sin(lat)*sin(f) + cos(lat)*cos(f)*cos(p3))

    return degrees(h0)

选取经纬120oE、30oN的位置，第61天（3月2日）来进行测试。

print(get_h0(120, 30, 60, 2))  # 02时
print(get_h0(120, 30, 60, 8))  # 08时
print(get_h0(120, 30, 60, 14))  # 14时
print(get_h0(120, 30, 60, 20))  # 20时

# -53.98284690230641  # 02时
# 21.349321051146905  # 08时
# 42.71321337396969  # 14时
# -29.641067525563916  # 20时

最终02时、08时、14时、20时太阳高度角依次为-53、21、42、-29，夜间太阳高度角为负值。

太阳高度角计算完毕~~