目前,我正在制作一个脚本,给定一组天体坐标,它会告诉你在接下来的几天里,当那个点对特定的望远镜可见时。标准很简单,在水平坐标系中,物体的高度必须在 30 到 65 度之间(这里的变量"暴击"表示这一点,但以弧度为单位)。所以我有一组名为"Ant"的望远镜参数,然后使用Pyephem:
#imported ephem as ep
obj= ep.FixedBody()
obj._ra= E.ra
obj._dec= E.dec
obj._epoch = E.epoch
Ant.date = ep.now()
for d in range(days):
for i in range(24):
for j in range (60):
Ant.date += ep.minute
obj.compute(Ant)
crit= float(obj.alt)
if crit>=0.523599 and crit <=1.13446:
print "Visible at %s" %Ant.date
这导致打印大量"可见于2016/7/11 19:41:21",每分钟打印1个。例如,我只想让它打印类似"在 2016/7/11 19:41:21 进入可见性,在 2016/7/11 23:41:00 离开"之类的内容。任何想法将不胜感激。
免责声明:对不起,不是母语为英语的人。
您需要跟踪它是否已经在范围内。因此,例如,在开始时,您将初始化它:
is_visible = False
您的 if 语句可能如下所示:
if crit>=0.523599 and crit <=1.13446:
if not is_visible:
print "Visible at %s" %Ant.date
is_visible = True
else:
if is_visible:
print "No longer visible at %s" % Ant.date
is_visible = False