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=(3.5, 3.5)) m.fig.show_rgb(pngpath) # m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}, # edgecolor=RED, lw=1.5) # m.plot_combined_fields({"INSTRUME": "WIRCam", # "TYPE": "sci", # "RUNID": {'$in': ['07BC20', '09BC29']}}, # edgecolor=RED, lw=3.) # m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}, # edgecolor=BLUE, lw=3.) # m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sky"}, # edgecolor=BLUE, lw=0.8) # m.plot_phat(union=True, edgecolor=GREEN, facecolor='None', # lw=2.) 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.plot_box(45. / 60. * 15., 42., 2.6, 1.9, color='k', lw=2.) # m.plot_box(37. / 60. * 15., 39.9, 2.6, 1.9, color='k', lw=2.) c1 = ICRSCoordinates(ra=(47. / 60. * 15.) * u.degree, dec=42.9 * u.degree) c2 = ICRSCoordinates(ra=M31RA0 * u.degree, dec=M31DEC0 * u.degree) c3 = ICRSCoordinates(ra=(38. / 60. * 15.) * u.degree, dec=39.6 * u.degree) m.plot_box(c1.ra.degree, c1.dec.degree, 2.6, 1.9, color='k', lw=2.) m.plot_box(c2.ra.degree, c2.dec.degree, 2.6, 1.9, color='k', lw=2.) m.plot_box(c3.ra.degree, c3.dec.degree, 2.6, 1.9, color='k', lw=2.) # m.save("oxford_map_zoomin.pdf", format='pdf', # dpi=300, transparent=True, adjust_bbox=True) m.fig.recenter(M31RA0, M31DEC0, radius=4) # , width=9, height=9) m.save("dragonfly_map.png", format='png', dpi=300, transparent=True, adjust_bbox=True) print c1.ra.hms, c1.dec.dms print c2.ra.hms, c2.dec.dms print c3.ra.hms, c3.dec.dms print "North: {0:8f} {1:8f}".format(c1.ra.degree, c1.dec.degree) print "Central: {0:8f} {1:8f}".format(c2.ra.degree, c2.dec.degree) print "South: {0:8f} {1:8f}".format(c3.ra.degree, c3.dec.degree)
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) # box(0:51:44.130,+41:11:54.49,3600",3600",0) # text={SE1} # box(0:50:02.796,+40:16:10.59,3600",3600",0) # text={SE2} # box(0:46:18.611,+39:21:42.43,3600",3600",0) # text={SE4} # box(0:51:07.630,+39:19:35.84,3600",3600",0) # text={SE3} # box(0:41:27.670,+43:08:55.63,3600",3600",0) # text={NW1} # box(0:36:23.986,+42:12:54.65,3600",3600",0) # text={NW2} # box(0:32:44.822,+41:15:51.34,3600",3600",0) # text={NW3} # box(0:31:14.572,+40:18:11.44,3600",3600",0) # text={NW4} # box(0:36:17.710,+43:08:30.21,3600",3600",0) # text={NW5} # box(0:31:27.235,+42:11:42.97,3600",3600",0) # text={NW6} # box(0:31:18.318,+43:06:58.70,3600",3600",0) # text={NW7} fields = [ {"ra": (0, 51.0, 44.130), "dec": (41, 11, 54.49), "n": "SE1"}, {"ra": (0, 50.0, 02.796), "dec": (40, 16, 10.59), "n": "SE2"}, {"ra": (0, 46.0, 18.611), "dec": (39, 21, 42.43), "n": "SE4"}, {"ra": (0, 51.0, 07.630), "dec": (39, 19, 35.84), "n": "SE3"}, {"ra": (0, 41.0, 27.670), "dec": (43, 8, 55.63), "n": "NW1"}, {"ra": (0, 36.0, 23.986), "dec": (42, 12, 54.65), "n": "NW2"}, {"ra": (0, 32.0, 44.822), "dec": (41, 15, 51.34), "n": "NW3"}, {"ra": (0, 31.0, 14.572), "dec": (40, 18, 11.44), "n": "NW4"}, {"ra": (0, 36.0, 17.710), "dec": (43, 8, 30.21), "n": "NW5"}, {"ra": (0, 31.0, 27.235), "dec": (42, 11, 42.97), "n": "NW6"}, {"ra": (0, 31.0, 18.318), "dec": (43, 6, 58.70), "n": "NW7"}, ] for field in fields: field["ra"] = 15.0 * (float(field["ra"][0]) + field["ra"][1] / 60.0 + field["ra"][1] / 3600.0) field["dec"] = float(field["dec"][0]) + field["dec"][1] / 60.0 + field["dec"][1] / 3600.0 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=(3.5, 3.5)) m.fig.show_rgb(pngpath) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sky", "RUNID": {"$in": ["07BC20", "07BH47"]}}) m.plot_combined_fields( {"INSTRUME": "WIRCam", "TYPE": "sky", "RUNID": {"$in": ["09BH52", "09BC29", "11BC12", "12BH04"]}} ) m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}) # m.plot_phat(union=True) # m.plot_phat_fields(union=False, edgecolor='None', facecolor='r', # zorder=10, alpha=0.2, band="F110W") m.plot_hst_halo() for field in fields: m.plot_box(field["ra"], field["dec"], width=1.0, height=1.0, facecolor="None", edgecolor="orange", lw=2.0) m.save("cfht_lsb_14b.pdf", dpi=300, transparent=True, adjust_bbox=True, format="pdf")
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=(3.5, 3.5)) m.fig.show_rgb(pngpath) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}, edgecolor=RED, lw=1.5) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci", "RUNID": {'$in': ['07BC20', '09BC29']}}, edgecolor=RED, lw=3.) # m.plot_fields({"INSTRUME": "WIRCam", "TYPE": "sky", # "RUNID": {"$in": ['07BC20', '07BH47']}}, # edgecolor=RED, lw=0.8) # m.plot_fields({"INSTRUME": "WIRCam", "TYPE": "sky", # "RUNID": {"$in": ['09BC29']}}, # edgecolor=RED, lw=0.8) m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}, edgecolor=BLUE, lw=3.) m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sky"}, edgecolor=BLUE, lw=0.8) m.plot_phat(union=True, edgecolor=GREEN, facecolor='None', lw=2.) 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.save("oxford_map_zoomin.pdf", format='pdf', dpi=300, transparent=True, adjust_bbox=True) m.fig.recenter(M31RA0, M31DEC0, radius=4) # , width=9, height=9) m.save("oxford_map.pdf", format='pdf', dpi=300, transparent=True, adjust_bbox=True)
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)
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=(3.5, 3.5)) m.fig.show_rgb(pngpath) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}, edgecolor="r") m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sky", "RUNID": {"$in": ["07BC20", "07BH47"]}}, edgecolor="r") m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}, edgecolor="b") m.plot_phat(union=True) m.plot_hst_halo(label=True) m.save("hst_footprint.pdf", format="pdf", dpi=300, transparent=True, adjust_bbox=True)
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=(3.5, 3.5)) m.fig.show_rgb(pngpath) m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sci"}, edgecolor='r') m.plot_combined_fields({"INSTRUME": "WIRCam", "TYPE": "sky", "RUNID": {"$in": ['07BC20', '07BH47']}}, edgecolor='r') m.plot_combined_fields({"INSTRUME": "MegaPrime", "lsb_mosaic.kind": "sci"}, edgecolor='b') m.plot_phat(union=True) # Plot archival fields m.plot_narrowband_fields(['SW1', 'SW2', 'NGC205'], facecolor='g', alpha=0.7) # Plot observed fields observed = ["AGB_%i" % i for i in range(1, 17)] m.plot_narrowband_fields(observed, facecolor='orange', alpha=0.6) # Plot un-observed fields unobserved = ["AGB_%i" % i for i in range(17, 21)] m.plot_narrowband_fields(unobserved, edgecolor='purple', linestyle='solid') m.save("narrowband_footprint.pdf", format='pdf', dpi=300, transparent=True, adjust_bbox=True)
# 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. \ + field['dec'][1] / 3600. if __name__ == '__main__': main()
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)