当然有合适的简便方案啦！你不用自己啃复杂的公式，Python里不管是用标准库还是专业地理库，都能轻松算出UTM坐标间的方位角～

先说明下：UTM是平面直角坐标系，方位角是水平面上的角度，所以Z坐标（高程）不影响计算结果，我们只需要用到X（东向）和Y（北向）坐标就行。

方案一：用Python标准库math直接计算（零依赖，最简便） #

因为UTM是平面投影，直接用平面几何的方法就能算出方位角，步骤很清晰：

1. 计算两点的东向差值（ΔE = X2 - X1）和北向差值（ΔN = Y2 - Y1）
2. 用atan2计算弧度值，再转成角度
3. 把角度调整到0-360度的范围（避免出现负数角度）

代码示例（代入你的坐标）：

import math

# 你的示例坐标
x1, y1, z1 = 10.3, 20.4, 4.8
x2, y2, z2 = 14.2, 35.1, 5.0

delta_e = x2 - x1
delta_n = y2 - y1

# 处理两点重合的特殊情况
if delta_e == 0 and delta_n == 0:
    raise ValueError("两点坐标完全相同，无法计算方位角")

# 计算弧度并转成角度
azimuth_rad = math.atan2(delta_e, delta_n)
azimuth_deg = math.degrees(azimuth_rad)

# 调整角度到0-360度区间
azimuth_deg = azimuth_deg if azimuth_deg >= 0 else azimuth_deg + 360

print(f"从点1到点2的方位角（正北顺时针）：{azimuth_deg:.2f} 度")

运行这段代码后，你会得到一个0-360度之间的方位角数值，完全符合需求。

方案二：用专业地理库pyproj（适合复杂地理坐标场景） #

如果之后你还需要处理更多地理坐标转换、大地测量相关的操作，推荐用pyproj这个专业库。它能处理UTM和经纬度的转换，再通过大地测量模型计算方位角，结果和平面计算几乎一致（UTM投影误差极小）。

步骤： #

1. 先安装库：pip install pyproj
2. 确定你的UTM带号和半球（比如UTM 32N对应的EPSG代码是32632，需要根据你的实际坐标修改）
3. 把UTM坐标转成WGS84经纬度，再用Geod类计算方位角

代码示例：

from pyproj import CRS, Transformer, Geod

# 定义坐标系：这里假设是UTM 32N WGS84（EPSG:32632），请根据你的坐标带号修改
utm_crs = CRS.from_epsg(32632)
wgs84_crs = CRS.from_epsg(4326)  # WGS84经纬度坐标系

# 创建UTM转经纬度的转换器
transformer = Transformer.from_crs(utm_crs, wgs84_crs, always_xy=True)

# 你的示例坐标
x1, y1, z1 = 10.3, 20.4, 4.8
x2, y2, z2 = 14.2, 35.1, 5.0

# 转换得到经纬度
lon1, lat1 = transformer.transform(x1, y1)
lon2, lat2 = transformer.transform(x2, y2)

# 用WGS84椭球计算方位角
geod = Geod(ellps='WGS84')
azimuth12, _, _ = geod.inv(lon1, lat1, lon2, lat2)

print(f"从点1到点2的方位角：{azimuth12:.2f} 度")

注意： #

如果不知道自己的UTM带号，可以通过坐标的X值判断（UTM带号范围1-60，X值一般在100000-900000之间，带号对应EPSG代码规则：北半球是326+带号，南半球是327+带号）。

总结 #

• 只是单纯计算这组坐标的方位角：选方案一，零依赖，代码简单直接
• 有更多地理坐标处理需求：选方案二，专业且扩展性强

内容的提问来源于stack exchange，提问作者Naveh Elya