def plot(galactic, rect): pf = MonteCarlo.make_allsky_isotropic_pyfits(galactic=galactic, proj='AIT') h,d=pf[0].header,pf[0].data ax = make_allsky_axes_from_header(fig, rect=rect, header=h, lon_center=0.) im = ax.imshow(pf[0].data, origin="lower", cmap=plt.cm.gray_r) ax.grid() overlay(ax)
def plot_allsky(fig, pyfits, label, imshow_kwargs=dict(), **kwargs): data = p[0].data if len(p[0].data.shape) == 3: data = data.sum(axis=0) print data.shape header = p[0].header ax = make_allsky_axes_from_header(fig, header=header, lon_center=0, **kwargs) ax.imshow(data, origin="lower", cmap=colormaps.b, **imshow_kwargs) P.title(label) return ax
proj_list = ["CYP", "CEA", "CAR", "MER", "SFL", "PAR", "MOL", "AIT"] DO_HEALPIX = True if DO_HEALPIX: healpix_data = get_LAB_healpix_data() else: healpix_data = None for proj in proj_list: fig = plt.figure() rect = 111 coord, lon_center = "fk5", 180 header = allsky_header(coord=coord, proj=proj, lon_center=lon_center, cdelt=0.2) ax = make_allsky_axes_from_header(fig, rect, header, lon_center=lon_center) ax.set_title("proj = %s" % proj, position=(0.5, 1.1)) do_allsky(ax, "gal") if healpix_data is not None: d = healpix_data.get_projected_map(header) im = ax.imshow(d**.5, origin="lower", cmap="gist_heat_r") c1, c2 = im.get_clim() im.set_clim(c1, c2*0.8) plt.show()
fits_name = "lambda_mollweide_halpha_fwhm06_0512.fits" f = pyfits.open(fits_name) data = f[1].data header = f[1].header # used fixed header del header["LONPOLE"] del header["LATPOLE"] fig = plt.figure(1, figsize=(8, 5)) fig.clf() from pywcsgrid2.allsky_axes import make_allsky_axes_from_header ax = make_allsky_axes_from_header(fig, rect=111, header=header, lon_center=0.) if data is not None: im = ax.imshow(data, origin="lower", cmap=plt.cm.gray_r) im.set_clip_path(ax.patch) im.set_clim(0, 30) ef = matplotlib.patheffects.withStroke(foreground="w", linewidth=3) axis = ax.axis["lat=0"] axis.major_ticklabels.set_path_effects([ef]) axis.label.set_path_effects([ef]) axis.set_zorder(5) ax.grid() plt.show()
proj_list = ["CYP", "CEA", "CAR", "MER", "SFL", "PAR", "MOL", ] DO_HEALPIX = True if DO_HEALPIX: healpix_data = get_LAB_healpix_data() else: healpix_data = None for proj in proj_list: fig = plt.figure() rect = 111 coord, lon_center = "fk5", 180 header = allsky_header(coord=coord, proj=proj, lon_center=lon_center, cdelt=0.2) ax = make_allsky_axes_from_header(fig, rect, header, lon_center=lon_center) ax.set_title("proj = %s" % proj, position=(0.5, 1.1)) do_allsky(ax, "gal") if healpix_data is not None: d = healpix_data.get_projected_map(header) im = ax.imshow(d**.5, origin="lower", cmap="gist_heat_r") c1, c2 = im.get_clim() im.set_clim(c1, c2*0.8) plt.show()