PyEphem：RA和分离问题

Question

我正在尝试使用PyEphem库进行天文学研究，并尝试在那里重新创建/理解一些基本计算。

有一个称为间隔的函数，给定两个行星并指定一个日期/时间，它可以计算出这两个行星之间相对于其在平面上的“ x”投影的间隔（角度）。

因此对于水星和火星行星的2018/1/1，我们有：

import ephem
import math 

mercury = ephem.Mercury('2018/1/1')
mars = ephem.Mars('2018/1/1')

s1 = ephem.separation(mercury, mars)
print(math.degrees(s1))

返回哪个

33.792384499568264

但是，如果我想在不使用分离函数的情况下进行计算，那么计算将很简单，即1的“右提升”减去另一个的“右提升”：

math.degrees(mercury.ra) - math.degrees(mars.ra)

返回哪个

35.114532008671574

为什么角度不同？根据PyEphem的说法，由于我不包括观察者的经纬度，因此所有计算都应该以地心为中心。

是否有人熟悉PyEphem背后的计算，或者是另一个具有内置星历表的库，可以产生一致的分离结果？

Answer 1

好问题！您将要查阅有关球三角学的教科书。正如您无法仅根据其经度来预测两个位置之间的英里数一样，因为它们可能都在赤道上并且相距很远，或者可能在北极附近并且彼此之间只有几英尺的距离，所以您无法预测仅将两个对象之间的分离与右提升（=“天空经度”）分开，但必须添加（a）两个对象的偏角，以及（b）使用特殊公式来计算球体上两个位置之间的弧度沿着大圈路线。