# Example using astLibs   (see http:// astlib.sourceforge.net/ )
#
# Takes a FITS image and clips a region and returns a thumbnail  
# orientated north up, east left   (from Jim Geach)

from astro_coords import *
from numpy.numarray import *
from astLib import astWCS
from astLib import astImages

image = argv[1]
r,d = argv[2:4]
win = float(argv[4])/3600.
out = argv[5]

def thumb(image, ra, dec,  window, output):
     try:
     	try: wcs = astWCS.WCS(image,mode='image')
     	except: wcs = astWCS.WCS 
(image,mode='image',extensionName='SCI') #usual failure is that there  
is a SCI extension
     	imdict = astImages.clipRotatedImageSectionWCS(getdata(image),  
wcs, ra, dec, window, returnWCS=True) #rotate N-up, E-left
     	thumb = imdict['clippedData'] #data array
     	newwcs = imdict['clippedWCS'] #wcs info
     	astImages.saveFITS(output, thumb, imageWCS=newwcs)
     	return thumb #return the 2d data array
     except:
         print 'error'
         return None #return zilch


if argv[6]=='sex': thumb(image,ra_to_deg(r),dec_to_deg(d),win,out)
if argv[6]=='deg': thumb(image,float(r),float(d),win,out)