def test_distance_1(self): tgt1 = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) tgt2 = entity.StaticTarget(name="altair", ra=altair[0], dec=altair[1]) time1 = self.obs.get_date("2010-10-18 22:30") d_alt, d_az = self.obs.distance(tgt1, tgt2, time1) self.assertEquals(str(d_alt)[:7], '-9.9657') self.assertEquals(str(d_az)[:7], '36.1910')
def test_airmass1(self): # now calculate via misc body = entity.StaticTarget(name='ACTJ0022-0036', ra='00:22:13.44', dec='-00:36:25.20') time1 = self.obs.get_date("2010-10-18 21:00") c1 = self.obs.calc(body, time1) self.assert_(c1.airmass > 1.13)
def test_observable_6(self): # vega should be visible near the beginning but not near the end # during this period (setting) tgt = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) time1 = self.obs.get_date("2014-04-29 10:00") time2 = self.obs.get_date("2014-04-29 11:00") is_obs, t1, t2 = self.obs.observable(tgt, time1, time2, 15.0, 85.0, 60 * 15) # 15 min NOT ok self.assert_(is_obs == False)
def test_observable_4(self): # vega should be visible near the beginning but not near the end # during this period (setting) tgt = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) time1 = self.obs.get_date("2014-04-29 09:30") time2 = self.obs.get_date("2014-04-29 10:30") is_obs, t1, T2 = self.obs.observable(tgt, time1, time2, 15.0, 85.0, 60 * 30) # 30 min ok self.assert_(is_obs == True)
def test_observable_3(self): # vega should be visible near the end but not in the beginning # during this period (rising) tgt = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) time1 = self.obs.get_date("2014-04-28 22:30") time2 = self.obs.get_date("2014-04-28 23:30") is_obs, t1, t2 = self.obs.observable(tgt, time1, time2, 15.0, 85.0, 60 * 50) # 50 min NOT ok self.assert_(is_obs == False)
def test_observable_1(self): # vega should be visible during this period tgt = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) time1 = self.obs.get_date("2014-04-29 04:30") time2 = self.obs.get_date("2014-04-29 05:30") is_obs, t1, t2 = self.obs.observable(tgt, time1, time2, 15.0, 85.0, 59.9 * 60) ## print((1, is_obs, ## t1.astimezone(self.obs.tz_local).strftime("%H:%M"), ## t2.astimezone(self.obs.tz_local).strftime("%H:%M"))) self.assert_(is_obs == True)
self._plot_target(observer, target, time_start, colors[i]) i = (i+1) % len(colors) self.redraw() if __name__ == '__main__': from obsplan import entity import pytz from matplotlib.backends.backend_qt4agg import FigureCanvasQTAgg as FigureCanvas from PyQt4 import QtGui import ephem app = QtGui.QApplication([]) plot = AZELPlot(10, 10) plot.setup() plot.plot_azel([(-210.0, 60.43, "telescope"),]) tgt3 = entity.StaticTarget(name="Bootes", ra="14:31:45.40", dec="+32:28:38.50") tz = pytz.timezone('US/Hawaii') site = entity.Observer('subaru', longitude='-155:28:48.900', latitude='+19:49:42.600', elevation=4163, pressure=615, temperature=0, timezone='US/Hawaii') start_time = datetime.strptime("2015-03-27 20:05:00", "%Y-%m-%d %H:%M:%S") start_time = tz.localize(start_time) plot.plot_targets(site, [entity.moon, entity.sun, tgt3], start_time, ['white', 'yellow', 'green'])
def test_get_tgt(self): tgt = entity.StaticTarget(name="vega", ra=vega[0], dec=vega[1]) self.assert_(isinstance(tgt.body, ephem.Body))
temperature=0, timezone='US/Hawaii') start_time = datetime.strptime("2015-03-30 18:30:00", "%Y-%m-%d %H:%M:%S") start_time = tz.localize(start_time) t = start_time # if schedule starts after midnight, change start date to the # day before if 0 <= t.hour < 12: t -= timedelta(0, 3600*12) ndate = t.strftime("%Y/%m/%d") targets = [] site.set_date(t) tgt = entity.StaticTarget(name='S5', ra='14:20:00.00', dec='48:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='Sf', ra='09:40:00.00', dec='43:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='Sm', ra='10:30:00.00', dec='36:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='Sn', ra='15:10:00.00', dec='34:00:00.00') targets.append(tgt) # make airmass plot plot.plot_targets(site, targets, tz) if outfile == None: canvas = plot.make_canvas() canvas.show() else: