importfits.py

Go to the documentation of this file.

#
# This file was generated using xslt from its XML file
#
# Copyright 2009, Associated Universities Inc., Washington DC
#
import sys
import os
from  casac import *
import string
from taskinit import casalog
#from taskmanager import tm
import task_importfits
def importfits(fitsimage='', imagename='', whichrep=0, whichhdu=0, zeroblanks=True, overwrite=False):

        """Convert an image FITS file into a CASA image

        Standard image FITS files are processed into a Casa image

        Keyword arguments:
        fitsimage -- Name of input image FITS file
                default: none; example='3C273XC1.fits'
        imagename -- Name of output CASA image
                default: none; example: imagename='3C273XC1.image'
        whichrep -- If fits image has multiple coordinate reps,
                choose one.
                default: 0 means first; example: whichrep=1
        whichhdu -- If fits file contains multiple images,
                choose this one
                default=0 mean first; example: whichhdu=1
        zeroblanks -- Set blanked pixels to zero (not NaN)
                default=True; example: zeroblanks=True
        overwrite -- Overwrite pre-existing imagename
                default=False; example: overwrite=True

 
        """

#
#    The following is work around to avoid a bug with current python translation
#
        mytmp = {}

        mytmp['fitsimage'] = fitsimage
        mytmp['imagename'] = imagename
        mytmp['whichrep'] = whichrep
        mytmp['whichhdu'] = whichhdu
        mytmp['zeroblanks'] = zeroblanks
        mytmp['overwrite'] = overwrite
        pathname='file:///'+os.environ.get('CASAPATH').split()[0]+'/share/xml/'
        trec = casac.utils().torecord(pathname+'importfits.xml')

        casalog.origin('importfits')
        if trec.has_key('importfits') and casac.utils().verify(mytmp, trec['importfits']) :
            result = task_importfits.importfits(fitsimage, imagename, whichrep, whichhdu, zeroblanks, overwrite)

        else :
          result = False
        return result