mapping - 使用测地坐标计算方位角

Question

我正在尝试计算具有由纬度/经度给出的测地线坐标的物体的方位角和指向不同纬度/经度坐标的另一个物体的航向。

问题的粗略草图是

其中三角形代表对象一，箭头代表对象一的方向，正方形代表对象二。

我知道这可以使用 hasrsine 公式来完成，如，但是如何在考虑球体地球的同时做到这一点。