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)
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 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)
# 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 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)
{'ra': (0., 30., 0.), 'dec': (36., 10., 0), 'n': 'S13'},
{'ra': (0., 18., 0.), 'dec': (41., 20., 0), 'n': 'S14'},
{'ra': (0., 28., 0.), 'dec': (43., 40., 0), 'n': 'S15'},
]
for field in sky_fields:
convert_field_coord(field)
for field in square_fields:
convert_field_coord(field)
for field in perimeter_fields:
convert_field_coord(field)
for field in new_sky_fields:
convert_field_coord(field)
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)
radii = np.arange(10., 70., 10.)
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:
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)
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)
