Example #1
0
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
Example #2
0
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