如何使用Python3计算三维空间两点间的方位角

阿华AIGC实验室

2026-4-30

计算三维空间两点间方位角的Python3实现

嘿，这个问题在地理空间计算或者3D开发场景里挺常见的～首先得明确：三维空间里的方位角，一般指的是两点在水平面上（XY平面）的投影连线，与正北方向（通常默认Y轴正方向为正北）顺时针旋转的夹角，范围是0°到360°（0°对应正北，90°正东，180°正南，270°正西）。下面一步步给你讲实现方法：

核心步骤拆解

第一步：提取水平坐标差
方位角只和水平位置有关，所以我们可以直接忽略Z轴高度值，计算两点在XY平面的坐标差值：
```
dx = x2 - x1  # 东向差值
dy = y2 - y1  # 北向差值
```
第二步：用math.atan2()计算基础弧度
Python的math.atan2()可以直接根据坐标差计算角度弧度，但要注意参数顺序：我们需要的是从正北（Y轴正方向）出发的角度，所以传入(dx, dy)会得到从Y轴正方向逆时针的弧度值：
```
import math
angle_rad = math.atan2(dx, dy)
```
第三步：转换为角度并调整到0-360°范围
atan2()返回的弧度范围是-π到π（对应-180°到180°），我们需要把负角度转换为正角度，确保结果落在标准的0-360°区间：
```
angle_deg = math.degrees(angle_rad)
azimuth = angle_deg if angle_deg >= 0 else angle_deg + 360
```

完整可复用函数

把上面的步骤整合为一个函数，还可以加入边界判断（比如两点水平重合的情况）：

import math

def calculate_azimuth(p1, p2):
    """
    计算三维空间中从p1到p2的方位角（0-360°，正北为0°）
    
    参数:
        p1 (tuple): 起点坐标 (x1, y1, z1)
        p2 (tuple): 终点坐标 (x2, y2, z2)
    
    返回:
        float/None: 方位角（单位：度），若两点水平重合则返回None
    """
    x1, y1, _ = p1
    x2, y2, _ = p2
    
    dx = x2 - x1
    dy = y2 - y1
    
    # 处理两点水平重合的特殊情况
    if dx == 0 and dy == 0:
        print("两点在水平面上重合，无法计算方位角")
        return None
    
    # 计算弧度并转换为角度
    angle_rad = math.atan2(dx, dy)
    angle_deg = math.degrees(angle_rad)
    
    # 调整到0-360°范围
    azimuth = angle_deg if angle_deg >= 0 else angle_deg + 360
    return azimuth

# 测试示例
if __name__ == "__main__":
    p1 = (0, 0, 5)
    p2 = (10, 10, 5)
    print(calculate_azimuth(p1, p2))  # 输出45.0（东北方向）
    
    p3 = (5, 5, 2)
    p4 = (5, 0, 8)
    print(calculate_azimuth(p3, p4))  # 输出180.0（正南方向）