description = """Plot the contents of a FITS image file.""" parser = argparse.ArgumentParser(usage=usage, description=description) parser.add_argument('files', nargs='+') parser.add_argument('--hdu', default=0, type=int, help='Set the HDU number to plot.') args = parser.parse_args() hdulist = pyfits.open(args.files[0]) im = FITSImage.createFromFITS(args.files[0], args.hdu) if isinstance(im, SkyCube): # Integrate over 3rd (energy) dimension im = im.marginalize(2) plt.figure() im.plot(cmap='ds9_b') plt.figure() # Smooth by 0.2 deg im.smooth(0.2).plot(cmap='ds9_b') # Draw an arbitrary contour in Galactic Coordinates
usage = "usage: %(prog)s [options] [FT1 file ...]" description = """Plot the contents of a FITS image file.""" parser = argparse.ArgumentParser(usage=usage,description=description) parser.add_argument('files', nargs='+') parser.add_argument('--hdu', default = 0, type=int, help = 'Set the HDU number to plot.') args = parser.parse_args() hdulist = pyfits.open(args.files[0]) im = FITSImage.createFromFITS(args.files[0],args.hdu) if isinstance(im,SkyCube): # Integrate over 3rd (energy) dimension im = im.marginalize(2) plt.figure() im.plot(cmap='ds9_b') plt.figure() # Smooth by 0.2 deg im.smooth(0.2).plot(cmap='ds9_b')
parser.add_argument('--srclabels_thresh', default=10.0,type=float) parser.add_argument('--format', default='png') parser.add_argument('--hdu', default = 0, type=int, help = 'Set the HDU of the input FITS file from which ' 'the image data will be extracted.') args = parser.parse_args() hdulist = pyfits.open(args.files[0]) model_hdu = None irf_version = None im_mdl = None im = FITSImage.createFromHDU(hdulist[args.hdu]) if args.model_file: model_hdu = pyfits.open(args.model_file)[0] if args.irf: irf_version = [args.irf,None] else: irf_version = get_irf_version(model_hdu.header) im_mdl = FITSImage.createFromHDU(model_hdu) m = re.search('(.+)V5_(.+)',irf_version[0]) if m and 'stacked' in args.model_file: irf_version[0] = m.groups()[0]+'V5' else: irf_version[1] = None
parser.add_argument('--model_file', default=None) parser.add_argument('--prefix', default=None) parser.add_argument('--hdu', default = 0, type=int, help = 'Set the HDU to plot.') args = parser.parse_args() hdulist = pyfits.open(args.files[0]) model_hdu = None irf_version = None im_mdl = None im = FITSImage.createFromHDU(hdulist[args.hdu]) if args.model_file: model_hdu = pyfits.open(args.model_file)[0] irf_version = get_irf_version(model_hdu.header) im_mdl = FITSImage.createFromHDU(model_hdu) #fp = FITSPlotter(im,im_mdl,None,args.prefix) #fp.make_plots_skycube(None,smooth=True,resid_type='fractional', # suffix='_data_map_resid_frac') #fp.make_plots_skycube(4,smooth=True,resid_type='fractional', # suffix='_data_map_slice_resid_frac')