python中的方位角计算(不带地理定位库)

import math

def calculate_bearing(lat1, lon1, lat2, lon2):
    # 将经纬度转换为弧度
    lat1 = math.radians(lat1)
    lon1 = math.radians(lon1)
    lat2 = math.radians(lat2)
    lon2 = math.radians(lon2)

    # 计算经度差值
    delta_lon = lon2 - lon1

    # 计算方位角
    y = math.sin(delta_lon) * math.cos(lat2)
    x = math.cos(lat1) * math.sin(lat2) - math.sin(lat1) * math.cos(lat2) * math.cos(delta_lon)
    bearing = math.atan2(y, x)

    # 将弧度转换为角度
    bearing = math.degrees(bearing)

    # 将方位角转换为0-360度范围
    bearing = (bearing + 360) % 360

    return bearing