def test_gsm_observer():
""" Test GSMObserver() is working
"""
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.view(logged=True)
ov.view_observed_gsm(logged=True)
def test_gsm_observer():
""" Test GSMObserver() is working
"""
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
ov.date = datetime(2000, 1, 1, 23, 0)
ov.generate(50)
ov.view(logged=True)
ov.view_observed_gsm(logged=True)
def test_observed_mollview():
""" Generate animated maps of observing coverage over 24 hours """
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
obs = []
for ii in range(0, 24, 1):
ov.date = datetime(2000, 1, 1, ii, 0)
ov.generate(50)
sky = ov.view_observed_gsm(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/galactic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'E'], min=9, max=20)
plt.savefig('generated_sky/ecliptic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'C'], min=9, max=20)
plt.savefig('generated_sky/equatorial-%02d.png' % ii)
plt.close()
ov.view(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/ortho-%02d.png' % ii)
plt.close()
print ii
os.system('convert -delay 20 generated_sky/ortho-*.png ortho.gif')
os.system('convert -delay 20 generated_sky/galactic-*.png galactic.gif')
os.system('convert -delay 20 generated_sky/ecliptic-*.png ecliptic.gif')
os.system('convert -delay 20 generated_sky/equatorial-*.png equatorial.gif')
def test_observed_mollview():
""" Generate animated maps of observing coverage over 24 hours """
(latitude, longitude, elevation) = ('37.2', '-118.2', 1222)
ov = GSMObserver()
ov.lon = longitude
ov.lat = latitude
ov.elev = elevation
obs = []
for ii in range(0, 24, 1):
ov.date = datetime(2000, 1, 1, ii, 0)
ov.generate(50)
sky = ov.view_observed_gsm(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/galactic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'E'], min=9, max=20)
plt.savefig('generated_sky/ecliptic-%02d.png' % ii)
plt.close()
hp.mollview(sky, coord=['G', 'C'], min=9, max=20)
plt.savefig('generated_sky/equatorial-%02d.png' % ii)
plt.close()
ov.view(logged=True, show=False, min=9, max=20)
plt.savefig('generated_sky/ortho-%02d.png' % ii)
plt.close()
print ii
os.system('convert -delay 20 generated_sky/ortho-*.png ortho.gif')
os.system('convert -delay 20 generated_sky/galactic-*.png galactic.gif')
os.system('convert -delay 20 generated_sky/ecliptic-*.png ecliptic.gif')
os.system(
'convert -delay 20 generated_sky/equatorial-*.png equatorial.gif')