def plot_archive_fields(): client = MongoClient(host='localhost', port=27017) c = client.m31.images 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=(7.5, 7.5)) m.fig.show_rgb(pngpath) m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}, edgecolor=BLUE, lw=3.) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci", "RUNID": {'$in': ['07BC20', '09BC29', '11BC12', '12BH04']}}, edgecolor=RED, lw=3.) # everyone *but* ANDROIDS observations sel = {"INSTRUME": "MegaPrime", "RUNID": {"$nin": androids_runids}, "footprint": {"$exists": 1}, "FILTER": {"$in": ['u', 'g', 'r', 'i']}} fieldnames = c.find(sel).distinct("OBJECT") for fieldname in fieldnames: s = dict(sel) s.update({"OBJECT": fieldname}) doc = c.find_one(s) ra = doc['RA_DEG'] dec = doc['DEC_DEG'] # m.plot_combined_fields(s) m.add_label(ra, dec, fieldname.replace(r"_", r"\_"), size=7, zorder=1000) 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=5) m.save("megacam_archive_fields.png", format='png', dpi=300, transparent=True, adjust_bbox=True)
prof = read_release() m.plot_xvista_profile_ellipse_grid(prof, radii, zorder=100, edgecolor='r') # m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}) # m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}) # m.plot_hst_halo() for field in sky_fields: m.plot_box(field['ra'], field['dec'], width=1., height=1., facecolor='None', edgecolor='orange', lw=1.) # m.add_label(field['ra'], field['dec'], field['n']) for field in perimeter_fields: m.plot_box(field['ra'], field['dec'], width=1., height=1., facecolor='dodgerblue', alpha=0.7, edgecolor='dodgerblue', lw=1.) m.add_label(field['ra'], field['dec'], field['n']) for field in new_sky_fields: m.plot_box(field['ra'], field['dec'], width=1., height=1., facecolor='None', alpha=1., edgecolor='dodgerblue', lw=1.) m.add_label(field['ra'], field['dec'], field['n']) m.fig.recenter(M31RA0, M31DEC0, radius=6) # , width=9, height=9) m.save("cfht_lsb_14b.png", dpi=300, transparent=True, adjust_bbox=True, format='png') def convert_field_coord(field): field['ra'] = 15. * (float(field['ra'][0]) + field['ra'][1] / 60. + field['ra'][1] / 3600.) field['dec'] = float(field['dec'][0]) + field['dec'][1] / 60. \
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)