gdal - PyEphem/GDAL 在特定投影的特定范围内查找角度和仰角

Question

例如，我正在创建一条 5000 米长的线，该线在具有高度的观察点与直接位于观察者和被观察物体之间的路径中的点之间。

就我而言，现在是太阳。

我找不到一种预先存在的方法来找到椭球上两点之间的斜率。

一种从 PyEphem 提供的任何结果或转换中找到全局椭球投影（如 WGS94 墨卡托）中的北角和斜率的方法会让我快速上手。