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:
