基本的pyephem代码

Question

我有一个非常简单的代码来获取太阳的经度，但是当我比较 Astrolog 和 Astrodienst 的输出时，它不正确，有 13 分钟的差异。我没有添加观察者，因为我认为默认是格林威治标准时间午夜（这是我想要的）。我做错了什么？

import ephem

start = ephem.date('2015/01/01')
end = ephem.date('2015/12/31')
f2 = open("Sun", 'w')
while start <= end:
    sun = ephem.Sun(start)
    ecl = ephem.Ecliptic(sun)
    f2.write(str(ephem.date(start))+' '+ str(ecl.lon) +'\n')
    start+=1
f2.close()

2015/12/30 的结果示例：

代码 - 2015/12/30 00:00:00 277:43:36.6

Astrodienst - 7°56'39 上限

谢谢

Answer 1

之所以相差13分钟是因为epoch的设置，我在添加的时候

太阳 = ephem.Sun（开始，纪元 = 开始）

结果与瑞士星历相同

基本的pyephem代码

Basic pyephem code

python-2.7

pyephem