def plot_known_tnos_singly(ax, extent, date): rate_cut = 'a > 15' print "PLOTTING LOCATIONS OF KNOWN KBOs (using {})".format( parameters.MPCORB_FILE) kbos = mpcread.getKBOs(parameters.MPCORB_FILE) print('Known KBOs: {}'.format(len(kbos))) retkbos = [] for kbo in kbos: kbo.compute(date) # # keep only the ones that'd make it onto this plot if not ((extent[0] <= math.degrees(kbo.ra) <= extent[1]) and (extent[2] <= math.degrees(kbo.dec) <= extent[3])): continue pos = (math.degrees(kbo.ra), math.degrees(kbo.dec)) ax.scatter(pos[0], pos[1], marker='x', facecolor='r') if len(kbo.name) >= 10: # sorted out for Sofia Pro Light ra_shift = -0.45 else: ra_shift = -0.22 ax.annotate(kbo.name, (pos[0] + ra_shift, pos[1] + 0.06), size=7, color='r') retkbos.append(kbo) print('Retained KBOs: {}'.format(len(retkbos))) return ax, retkbos
def plot_known_tnos(ax, ra_cen, dec_cen, width, height, date=newMoons['Oct13']): MPCORB = 'MPCORB.DAT' # TNOs file only: http://www.minorplanetcenter.net/iau/MPCORB/Distant.txt rate_cut = 'a > 15' if os.access(MPCORB,os.F_OK): kbos = mpcread.getKBOs(MPCORB, cond=rate_cut) for kbo in kbos: kbo.compute(ephem.date(date)) ax.scatter(math.degrees(kbo.ra), math.degrees(kbo.dec), marker='+', facecolor='none', edgecolor='g') return ax
def plot_known_tnos_singly(ax, date): MPCORB = 'MPCORB.DAT' # TNOs file only: http://www.minorplanetcenter.net/iau/MPCORB/Distant.txt rate_cut = 'a > 15' if os.access(MPCORB,os.F_OK): kbos = mpcread.getKBOs(MPCORB, cond=rate_cut) for kbo in kbos: kbo.compute(ephem.date(date)) pos = (math.degrees(kbo.ra), math.degrees(kbo.dec)) ax.scatter(pos[0], pos[1], marker='x', facecolor='r') ax.annotate(kbo.name, (pos[0]+.4, pos[1]+0.1), size='xx-small', color='r') return ax
def plot_known_tnos_batch(handles, labels, date): rate_cut = 'a > 15' if os.access(MPCORB, os.F_OK): kbos = mpcread.getKBOs(MPCORB, cond=rate_cut) kbo_ra = [] kbo_dec = [] for kbo in kbos: kbo.compute(ephem.date(date)) kbo_ra.append(math.degrees(kbo.ra)) kbo_dec.append(math.degrees(kbo.dec)) ax.scatter(kbo_ra, kbo_dec, marker='+', facecolor='none', edgecolor='g') return ax
def plot_known_tnos_batch(handles, labels, date): MPCORB = 'MPCORB.DAT' # TNOs file only: http://www.minorplanetcenter.net/iau/MPCORB/Distant.txt rate_cut = 'a > 15' if os.access(MPCORB,os.F_OK): kbos = mpcread.getKBOs(MPCORB, cond=rate_cut) kbo_ra = [] kbo_dec = [] for kbo in kbos: kbo.compute(ephem.date(date)) kbo_ra.append(math.degrees(kbo.ra)) kbo_dec.append(math.degrees(kbo.dec)) handles.append(plt.scatter(kbo_ra, kbo_dec, marker='+', facecolor='none', edgecolor='g')) labels.append('known TNOs') return handles, labels
dec = t.array['DEJ2000'] rects = [Rectangle(xy=(ra[idx] - camera_width_36 / 2.0, dec[idx] - camera_width_36 / 2.0), height=camera_height, width=camera_width_36, edgecolor='m', alpha=0.1, lw=0.1, zorder=-100, fill=False) for idx in xrange(ra.size)] for r in rects: ax.add_artist(r) if PLOT_MPCORB: print "PLOTTING LOCATIONS OF KNOWN KBOs (using {})".format(MPCORB_FILE) kbos = mpcread.getKBOs(MPCORB_FILE) for kbo in kbos: kbo.compute(plot_date) if field_polygon.isInside(math.degrees(float(kbo.ra)), math.degrees(float(kbo.dec))): ax.scatter(math.degrees(kbo.ra), math.degrees(kbo.dec), marker='x', s=4, facecolor='r', edgecolor='r', alpha=0.3) if PLOT_REAL_KBOS: reader = mpc.MPCReader() print "PLOTTING LOCATIONS OF OSSOS KBOs (using directory {})".format(REAL_KBO_AST_DIR) for ast in os.listdir(REAL_KBO_AST_DIR): if not ast.startswith('u'):
width=camera_dimen, edgecolor='k', alpha=0.1, lw=0.1, zorder=-100, fill=False) for idx in xrange(ra.size)] for r in rects: ax.add_artist(r) ra_cen = math.radians(ra_cen) dec_cen = math.radians(dec_cen) width = math.radians(width) height = math.radians(height) MPCORB = os.path.join(os.getenv('HOME','/Users/jjk'), 'MPCORB.DAT') if os.access(MPCORB,os.F_OK): kbos = mpcread.getKBOs(MPCORB) for kbo in kbos: kbo.compute(ephem.date(newMoons['Oct13'])) if not ( kbo.ra < ra_cen + width/2.0 and kbo.ra > ra_cen - width/2.0 and kbo.dec < dec_cen + height/2.0 and kbo.dec > dec_cen - height/2.0 ) : continue ax.scatter(math.degrees(kbo.ra), math.degrees(kbo.dec), marker='h', s=10, facecolor='none', edgecolor='g') savefig('layout.pdf')
rects = [ Rectangle(xy=(ra[idx] - dimen / 2.0, dec[idx] - dimen / 2.0), height=camera_dimen, width=camera_dimen, edgecolor='m', alpha=0.1, lw=0.1, zorder=-100, fill=False) for idx in xrange(ra.size) ] for r in rects: ax.add_artist(r) if PLOT_MPCORB: print "PLOTTING LOCATIONS OF KNOWN KBOs (using {})".format(MPCORB_FILE) kbos = mpcread.getKBOs(MPCORB_FILE) for kbo in kbos: kbo.compute(plot_date) if not ((ra_cen + width / 2.0 > kbo.ra > ra_cen - width / 2.0) and (dec_cen + height / 2.0 > kbo.dec > dec_cen - height / 2.0)): continue ax.scatter(math.degrees(kbo.ra), math.degrees(kbo.dec), marker='h', s=1, facecolor='g', edgecolor='g', alpha=0.3) if PLOT_REAL_KBOS: print "PLOTTING LOCATIONS OF OSSOS KBOs (using directory {})".format(