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
