在Python中，我们可以利用数学库如`math`来编写函数以计算两点间的**方位角（Azimuth Angle）** 和 **俯仰角（Elevation/Pitch Angle）**。以下是具体的解释： ### 背景知识 1. **方位角 (Azimuth)**：从某一点到另一点的方向角度，通常相对于正北方向顺时针测量。 2. **俯仰角 (Pitch/Elevation)**：表示视线相对水平面的倾斜程度。 对于给定的纬度和经度坐标A(lat_a, lon_a)、B(lat_b, lon_b)，以及两者的高度差height_diff，可以按照下面步骤进行计算。 --- ### Python 实现 ```python import math def calculate_angles(lat_a, lon_a, height_a, lat_b, lon_b, height_b): # 地球半径 (单位: 公里) R = 6371 # 将经纬度转换为弧度制 phi_A = math.radians(lat_a) lambda_A = math.radians(lon_a) phi_B = math.radians(lat_b) lambda_B = math.radians(lon_b) # 计算两地之间的Δλ（经度差异） delta_lambda = lambda_B - lambda_A # 使用Haversine公式的一部分获得 Δψ（大地线距离对应的圆心角） y = math.sin(delta_lambda) * math.cos(phi_B) x = math.cos(phi_A)*math.sin(phi_B) - \ math.sin(phi_A)*math.cos(phi_B)*math.cos(delta_lambda) azimuth_rad = math.atan2(y, x) # 方位角(弧度值) azimuth_deg = math.degrees(azimuth_rad) % 360 # 规范化到 [0, 360] # 高度差 height_diff = height_b - height_a # 水平距离 (基于 Haversine formula 近似求解平面距离) a = pow(math.sin((phi_B - phi_A)/2), 2) + math.cos(phi_A) * math.cos(phi_B) * pow(math.sin(delta_lambda / 2), 2) c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a)) horizontal_distance = R * c # 单位: km pitch_angle = math.atan2(height_diff, horizontal_distance*1000) # 弧度制俯仰角 pitch_deg = math.degrees(pitch_angle) # 度数形式 return {"方位角": round(azimuth_deg, 4), "俯仰角": round(pitch_deg, 4)} # 示例调用： result = calculate_angles( lat_a=52.2296756, lon_a=21.0122287, height_a=0, lat_b=41.8919300, lon_b=12.5113300, height_b=1000 ) print(result) ``` 此程序将返回 A 点对 B 点的位置关系，包括**方位角**和**俯仰角**的结果字典。其中方位角取值范围[0°~360°]；而俯仰角则是根据垂直高度差与实际地表距离得出的角度结果。 --- ### 输出示例 假设输入数据分别为 Warsaw(Poland): `(lat_a=52.2296756, lon_a=21.0122287)` 及 Rome(Italy):`(lat_b=41.8919300,lon_b=12.5113300,height_b=1000meters)` 运行上述代码会得到类似如下结果: ``` {'方位角': 187.9917, '俯仰角': 0.5949} ``` 这表明从华沙看罗马位于其南偏西约187.99度，并且有接近+0.59度的小量上扬视角。 ---