利用python计算雷达遮挡角

考虑标准大气折射和等效地球半径（8500km），雷达遮挡角计算公式用python计算程序如下

"""定义计算雷达遮蔽角的函数"""
def arfa_rda(h_grd, distance, h_rda):
    arfa = math.atan((math.sqrt(math.pow((h_grd +8500000),2)-\
                math.pow(distance,2))-(h_rda+8500000))/distance)
    return arfa   #返回为弧度值

其中h_grd为地形高度，h_rda为雷达站的高度，distance为两点距离，三者单位均为米。

参考文献：

周嘉健, 徐黄飞, 邹庆彪, 汤晶晶, 刘艳中. 基于SRTM地形数据天气雷达地形遮挡分析系统开发及应用[J]. 热带气象学报, 2021, 37(2): 258-267. doi: 10.16032/j.issn.1004-4965.2021.024