我是PyEphem的新手,我试图弄清楚它可以做什么以及它是如何工作的。由于我不想将它用作黑盒,并且盲目信任我得到的任何数字,所以我想重新创建一个here解释的示例。
该例子在1998年8月10日23:10UT处计算给定观察者的物体的方位角和高度。给出以下参数:
RA = 16小时41.7分钟,DEC = 36分钟28分钟
观测者的纬度是北纬52°30分钟,西经1°55分钟。
根据示例(我可以在Excel中重新创建)的正确答案是AZ = 269.14634度和ALT = 49.169122度。
我使用pyephem编写了以下代码,试图达到相同的结果:
day = '1998/8/10 23:10:00'
longitude = ephem.degrees('-1.91667')
latitude = ephem.degrees('52.5')
star = ephem.FixedBody()
star._ra = '16:41:42.0'
star._dec = '36:28:00.0'
observer = ephem.Observer()
observer.date = day
observer.lon = longitude
observer.lat = latitude
star.compute(observer)
print 'Observer', observer
print 'RA', star.ra, 'DEC', star.dec
print 'AZ', star.az, 'ALT', star.alt运行这个程序给了我这个输出:
>>>
Observer
RA 16:41:39.23 DEC 36:28:33.5
AZ 269:09:54.9 ALT 49:10:57.7AZ + ALT的结果显然与榜样相差甚远,但却不尽相同。我也对这样的事实感到困惑,即与我输入的内容相比,RA和DEC在打印输出时略有修改。
如果有人能够帮助我阐明为什么结果会有所不同,以及我可以或应该做些什么来复制结果,我将不胜感激。谢谢。
编辑:更正了以下答案中指出的错字。这个问题仍然有效。
编辑2:好的,我已阅读(并理解)为什么赤字和赤纬是由PyEphem从这个link中调整的。我不明白的是,是否有任何方法让PyEphem忽略相对论偏转,章动和光的像差调整,就像忽略大气折射一样?我假设方位角的差异是由于RA和DEC的调整,但它确实很好。