Public Function 求方位角(ByVal 经度1 As Double, ByVal 纬度1 As Double, ByVal 经度2 As Double, ByVal 纬度2 As Double) As Integer
Dim 平均纬度 As Double
Dim 求方位角s As Double
平均纬度 = (纬度1 + 纬度2) / 2
If 纬度1 - 纬度2 = 0 Then
求方位角s = 90
Else
求方位角s = Atn((经度1 - 经度2) * Cos(弧度(平均纬度)) / (纬度1 - 纬度2)) * 180 / pi
End If
If 纬度1 > 纬度2 Then
求方位角s = 求方位角s + 180
End If
If 求方位角s < 0 Then
求方位角s = 360 + 求方位角s
End If
求方位角 = Round(求方位角s, 0)
End Function