def test_distance_1(self): tgt1 = entity.StaticTarget("vega", vega[0], vega[1]) tgt2 = entity.StaticTarget("altair", altair[0], altair[1]) time1 = self.obs.get_date("2010-10-18 22:00") d_alt, d_az = self.obs.distance(tgt1, tgt2, time1) self.assertEqual(int(math.fabs(d_alt)), 11) self.assertEqual(int(math.fabs(d_az)), 38)
def test_airmass2(self): # now calculate via misc body = entity.StaticTarget('ACTJ0022-0036', '00:22:13.44', '-00:36:25.20') time1 = self.obs.get_date("2010-10-18 21:00") c1 = self.obs.calc(body, time1) self.assertTrue(c1.airmass > 1.2)
def test_observable_1(self): # vega should be visible during this period tgt = entity.StaticTarget("vega", vega[0], vega[1]) time1 = self.obs.get_date("2014-04-29 04:00") time2 = self.obs.get_date("2014-04-29 05:00") is_obs, time_rise, time_set = self.obs.observable( tgt, time1, time2, 15.0, 85.0, 59.9 * 60) print((1, is_obs, time_rise, time_set)) self.assertTrue(is_obs == True)
def test_observable_6(self): # vega should be visible near the beginning but not near the end # during this period (setting) tgt = entity.StaticTarget("vega", vega[0], vega[1]) time1 = self.obs.get_date("2014-04-29 11:00") time2 = self.obs.get_date("2014-04-29 12:00") is_obs, time_rise, time_set = self.obs.observable( tgt, time1, time2, 15.0, 85.0, 60 * 1) # 1 min NOT ok print((6, is_obs, time_rise, time_set)) self.assertTrue(is_obs == False)
def test_observable_2(self): # vega should be visible near the end but not in the beginning # during this period (rising) tgt = entity.StaticTarget("vega", vega[0], vega[1]) time1 = self.obs.get_date("2014-04-28 22:00") time2 = self.obs.get_date("2014-04-28 23:00") is_obs, time_rise, time_set = self.obs.observable( tgt, time1, time2, 15.0, 85.0, 60 * 15) # 15 min ok print((2, is_obs, time_rise, time_set)) self.assertTrue(is_obs == True)
def ftest_airmass(self): # calculate airmass via "observer" module import observer obs = observer.Observer('subaru') obs.almanac('2010/10/18') tgt = observer.tools.Target('ACTJ0022-0036', '00:22:13.44', '-00:36:25.20') am = observer.tools.Airmass(obs.almanac_data, tgt) time1 = self.obs.get_date("2010-10-18 22:00") time1_ut = time1.astimezone(tz.UTC) tup = am.compute_one(tgt.target, time1_ut) amass = tup[4] # now calculate via misc body = entity.StaticTarget('ACTJ0022-0036', '00:22:13.44', '-00:36:25.20') time1 = self.obs.get_date("2010-10-18 22:00") c1 = self.obs.calc(body, time1) self.assertTrue(math.fabs(amass - c1.airmass) < 0.01)
def test_get_body(self): tgt = entity.StaticTarget("vega", vega[0], vega[1]) self.assertTrue(isinstance(tgt.body, calcpos.Body))
start_time = datetime.strptime("2019-06-28 17:00:00", "%Y-%m-%d %H:%M:%S") start_time = start_time.replace(tzinfo=timezone) 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") print('what is t={}'.format(t)) targets = [] site.set_date(start_time) tgt = entity.StaticTarget(name='aaaa', ra='12:34:56.789', dec='12:34:56.78') targets.append(tgt) tgt = entity.StaticTarget(name='bbbb', ra='23:23:23.123', dec='45:45:45.34') targets.append(tgt) # 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)
plot.setup() site = get_site('subaru') tz = site.tz_local 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) tgt = entity.StaticTarget(name='Sn!4$#$@$', ra='20:10:00.00', dec='34:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='Sssn', ra='15:10:00.00', dec='24:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='nnnn', ra='15:10:00.00', dec='14:00:00.00') targets.append(tgt) tgt = entity.StaticTarget(name='mmmmm', ra='05:10:00.00', dec='14:00:00.00') targets.append(tgt)
plot.setup() site = get_site('subaru') tz = site.tz_local 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 num_tgts = len(targets) target_data = [] lengths = [] if num_tgts > 0: for tgt in targets: info_list = site.get_target_info(tgt)
from ginga import toolkit toolkit.use('qt') from ginga.gw import Widgets, Plot plot = AZELPlot(1000, 1000) plot.setup() app = Widgets.Application() topw = app.make_window() plotw = Plot.PlotWidget(plot) topw.set_widget(plotw) topw.add_callback('close', lambda w: w.delete()) topw.show() plot.plot_azel([(-210.0, 60.43, "telescope"),]) tgt3 = entity.StaticTarget(name="Bootes", ra="14:31:45.40", dec="+32:28:38.50") site = get_site('subaru') tz = site.tz_local 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, [common.moon, common.sun, tgt3], start_time, ['white', 'yellow', 'green']) app.mainloop() #END