pacs.py

#!/usr/bin/python
"""
Quick view of PACS data
"""

from herschel import pacsfits
from astropy.io import fits
import argparse
import aplpy
import numpy as np

parser = argparse.ArgumentParser()
parser.add_argument('obsid', default="")
parser.add_argument('-b', '--band', default="blue", choices=("blue", "red"),
                help="PACS band")
parser.add_argument('--datadir', default='./')
args = parser.parse_args()

fitsfile = pacsfits(args.datadir, args.obsid, args.band[0])
hdus = fits.open(fitsfile)

gc = aplpy.FITSFigure(fitsfile, hdu=1)
ra = hdus[0].header['RA_NOM']
dec = hdus[0].header['DEC_NOM']
# Center at the comet position k with radius in degrees
gc.recenter(ra, dec, radius=0.01)
# gc.show_colorscale(pmin=20, vmid=-5e-5, stretch="log")
# gc.show_colorscale(stretch="log", vmid=-0.0001)
# gc.show_colorscale(stretch="arcsinh")
gc.show_contour(colors="white")
gc.show_colorscale(vmin=1e-4)
gc.show_markers(ra, dec, marker="+", facecolor="white", edgecolor="white")
gc.add_colorbar()
gc.save('{}_{}.png'.format(args.obsid, args.band))