def get_discoveries(blockID='O13AE'):
retval = []
path = 'vos:OSSOS/measure3/2013A-E/track/submitted/' # rather than discoveries, for longer arcs
# path = 'vos:OSSOS/measure3/2013A-E/track/radec_corrected/'
discoveries = storage.listdir(path)
names = get_names(path, blockID)
print 'Discoveries:', len(names)
for i, kbo in enumerate(names):
try:
arclen, orbit = parse(kbo, discoveries, path=path)
retval.append(orbit)
except:
continue
return retval
def plot_ossos_discoveries(ax, blockID='O13AE', date="2013/04/09 08:50:00"):
path = 'vos:OSSOS/measure3/2013A-E/track/discoveries/'
discoveries = storage.listdir(path)
names = get_names(path, blockID)
for kbo in names:
arclen, orbit = parse(kbo, discoveries, path=path)
if date != "2013/04/09 08:50:00": # date is the new moon for a given month
orbit.predict(date.replace('/', '-'))
ra = orbit.coordinate.ra.degrees
dec = orbit.coordinate.dec.degrees
else: # specific date on which discovery was made: use discovery locations
ra = orbit.observations[0].coordinate.ra.degrees
dec = orbit.observations[0].coordinate.dec.degrees
ax.scatter(ra, dec, marker='+', facecolor='k')
ax.annotate(kbo.replace(blockID,''),
(orbit.coordinate.ra.degrees+.05, orbit.coordinate.dec.degrees-0.2),
size='xx-small',
color='k')
return ax