#import matplotlib
#matplotlib.use('Agg')

import pyfits
import numpy as np
from glob import glob
from astrometry.util.pyfits_utils import *
from astrometry.util.starutil_numpy import *

if __name__ == '__main__':
	for i in range(19, 7, -1):
		fn = 'index-40%02i.fits' % i
		T = fits_table(fn, hdunum=13)
		print 'Index', fn, 'has % 8i' % len(T), 'stars, J <', np.max(T.j_mag),
		hdr = pyfits.open(fn)[0].header
		lo = rad2deg(hdr['SCALE_L'])
		hi = rad2deg(hdr['SCALE_U'])
		print ', and features %.3g to %.3g deg' % (lo, hi),
		print '(%.3g to %.3g arcmin)' % (lo*60., hi*60.)
		
	for i in range(7, -1, -1):
		fns = glob('index-40%02i-*.fits' % i)
		fns.sort()
		N = 0
		for fn in fns:
			T = fits_table(fn, hdunum=13)
			print 'Index', fn, 'has % 8i' % len(T), 'stars, J < %.3f' % np.max(T.j_mag),
			hdr = pyfits.open(fn)[0].header
			lo = rad2deg(hdr['SCALE_L'])
			hi = rad2deg(hdr['SCALE_U'])
			print ' features %.3g to %.3g arcmin' % (lo*60., hi*60.)
			N += len(T)
		print 'Total', N, 'stars'