遮蔽角计算公式？

考虑标准大气折射和等效地球半径（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为两点距离，三者单位均为米。