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()
