##^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^Lance Simms, Stanford University 2009
#! /usr/bin/env python

##This script will run through all of the files specified on the command line
##
##Run either as 
##	imstat_example.py [filename]
#	imstat_example.py [@list_of_files]
import sys,os
from pyraf import iraf
import pyfits
from numarray import *
from pylab import *

##Tell a little about where this file is running from and indicated the 
##operating system class
print 'sys.argv[0] =', sys.argv[0]
print 'sys.argv[1] =', sys.argv[1]
pathname=os.path.dirname(sys.argv[0])
print 'path =', pathname 
print 'full path =', os.path.abspath(pathname)

def run_imstat(input):			#Function Definition
	iraf.images(_doprint=0)		#load the image package
	l=len(input)			#input is a list containing filename
	
	for image in input:
		
		s=iraf.imstat(image,Stdout=1)      #S is huge string containing
		lens=len(s)		           #all output from imstat
		print "Total Number Files:",lens   #lens=length of string
		datamax=zeros(lens,Float32)
		i=3

		for i in range(lens-1): 
			#Now we have array containing five strings
			#[Filename,NPIX,MEAN,STDEV,MIN,MAX]
			Values=s[i+1].split()
			max_count=float(Values[5])
			datamax[i]=max_count
		
	return datamax

if __name__== "__main__":
	datamax=run_imstat(sys.argv[1:]) 

plot(datamax)
xlabel('File Number')
ylabel('Maximum Pixel Value (Counts)')
title('Values For 5-12-05')
#savefig('Datamax.ps')
show()