python坐标方位角的计算

在Python中，计算坐标方位角通常涉及到地理空间分析，特别是当处理地图投影中的方向时。方位角（Azimuth）是指从某一点到另一点的直线角度，通常是以正北（0°）为基准，顺时针测量到该直线的角度。以下是一个简单的例子，假设我们有两点的经纬度：

from math import atan2, degrees

def calculate_azimuth(point1, point2):
    # 假设point1是经度lat1, 纬度lon1，point2是经度lat2, 纬度lon2
    lat1, lon1 = point1
    lat2, lon2 = point2
    
    # 将纬度转换成弧度
    y = lat1 * math.pi / 180
    y2 = lat2 * math.pi / 180
    
    # 经度差转换成弧度
    x = (lon2 - lon1) * math.pi / 180
    
    # 使用反正切函数计算方位角
    azimuth = degrees(atan2(math.sin(x), math.cos(y) * math.tan(y2)))
    
    # 如果结果在西半球（负经度），需要将其调整到0-360范围内
    if lon1 > lon2 or (lon1 == lon2 and lat1 < lat2):  # 逆时针方向
        azimuth += 360
    else:
        azimuth %= 360
    
    return azimuth

# 示例
point1 = [51.5074, -0.1278]  # 伦敦坐标
point2 = [40.7128, -74.0060]  # 纽约坐标
azimuth_degrees = calculate_azimuth(point1, point2)
print(f"方位角：{azimuth_degrees} 度")