python怎样使用各个日期赤纬_Python PyEphem计算方位角和海拔高度

本文介绍了作者作为PyEphem新手尝试理解和重现一个使用PyEphem计算1998年8月10日特定时刻物体方位角和高度的例子。文章详细展示了代码实现和预期结果之间的差异,探讨了可能的原因如相对论偏转、章动和光的像差调整,并寻求帮助以理解如何获得准确结果。
摘要由CSDN通过智能技术生成

我是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的调整,但它确实很好。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值