mk_thumbs

gistfile1.py

from astLib import astImages
from astLib import astWCS

def mk_thumb_wcs(header,data,ra,dec,window,output):

    try:
        wcs = astWCS.WCS(header,mode='pyfits')
    except:
        wcs = astWCS.WCS(header,mode='pyfits',extensionName='SCI')

    imdict = astImages.clipImageSectionWCS(data, wcs, ra, dec,
        window, returnWCS=True)
    thumb = imdict['data'] #data array
    newwcs = imdict['wcs'] #wcs info
    astImages.saveFITS(output, thumb, imageWCS=newwcs)
    return thumb