如何使用pyephem查找给定altaz坐标的北角与地平线之间的夹角？

Question

我在三脚架上用相机拍摄了许多天文照片。使用气泡水平仪确保画框的长边与水平线平行，并且我知道每张照片中心的alt / az（和赤道）坐标。现在，我正在编写一些python代码，在每个图像上覆盖一个指示器以标记北向。在给定的alt / az坐标下，我可以使用pyephem获取北极点方向与水平方向之间的角度吗？有什么线索吗？

Answer 1

仅仅为了获得正确的方向，所有这些都可以简化为：

    positionAngle = np.degrees(np.arctan2(CD12, CD11))

因为NumPy（np）可以在所有情况下（象限）正确处理反正切函数。

编辑：如果您想检测图像是否被镜像：计算矩阵的行列式

mirrored = (CD11 * CD22 - CD12 * CD21) < 0

并检查其是否小于0。