def main(): log = logging.getLogger('andromap') log.setLevel(logging.DEBUG) ch = logging.StreamHandler() ch.setLevel(logging.DEBUG) formatter = logging.Formatter( '%(asctime)s - %(name)s - %(levelname)s - %(message)s') ch.setFormatter(formatter) log.addHandler(ch) pngpath = os.path.expanduser("~/andromap/Elixir_B3_r.resamp.inverted.png") fitspath = os.path.expanduser("~/andromap/Elixir_B3_r.resamp.fits") m = Andromap(fitspath, figsize=(6.5, 6.5)) m.fig.show_rgb(pngpath) client = MongoClient(host='localhost', port=27017) c = client.m31.images # WIRCam mosaic footprint. m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci", "RUNID": {'$in': ['07BC20', '09BC29']}}, edgecolor=RED, lw=1, alpha=0.8) wircam_sky_sel = {"INSTRUME": "WIRCam", "TYPE": "sky", "RUNID": {"$in": ['07BC20', '07BH47', '09BC29']}} for field in c.find(wircam_sky_sel).distinct('OBJECT'): m.plot_combined_fields( {"INSTRUME": "WIRCam", "TYPE": "sky", "OBJECT": field}, edgecolor=ORANGE, lw=1, alpha=0.8) for field in MEGACAM_FIELD_NAMES: m.plot_combined_fields( {"INSTRUME": "MegaPrime", 'OBJECT': field}, edgecolor=BLUE, lw=2.) coord = m.compute_mean_coordinate({"OBJECT": field}) m.add_label(coord.ra.value, coord.dec.value, field.split('_')[-1], size=14, weight='heavy', ha='center', va='center', color=BLUE) for field in MEGACAM_SKY_NAMES: m.plot_combined_fields( {"INSTRUME": "MegaPrime", 'OBJECT': field}, edgecolor=BLUE, lw=1) coord = m.compute_mean_coordinate({"OBJECT": field}) m.add_label(coord.ra.value, coord.dec.value, field.split('_')[-1], size=14, weight='heavy', ha='center', va='center', color=BLUE) m.plot_phat(union=True, edgecolor=GREEN, facecolor='None', lw=1, alpha=0.8) radii = np.arange(10., 50., 10.) prof = read_release() m.plot_xvista_profile_ellipse_grid(prof, radii, zorder=100, edgecolor='0.2', alpha=0.7) m.fig.recenter(M31RA0, M31DEC0, radius=4) # , width=9, height=9) m.save("megacam_lsb_fields.pdf", format='pdf', dpi=300, transparent=True, adjust_bbox=True)