sys.exit(1) # Loop over fits files on command line for file in sys.argv[1:]: print file # If file does not exist, do not go further if not os.path.exists(file): WriteLog(file, 'does not exist') continue # Extract fits header values using pyfits # header = [<instrument>,<tbl_prefix>,primaryHDU,ext1HDU,ext2HDU,..,extNHDU] try: header = extractHeader(str(file)) except HeaderException, e: WriteLog(file, str(e)) continue # Pyfits do not recognize this a a fits file. Make a note and skip file. if header[0] == 'ERROR': WriteLog(file, 'not a valid fits file') continue instrument = header[0] tbl_prefix = header[1] # filename (stripped of .fits extension) fits_nm = str(os.path.basename(file))[:-5]
# file path + file name fits_file = str(file) # filename fits_nm = str(os.path.basename(file)) # Check file extention corresponds to a fits file if not (fits_nm[-5:] == '.fits'): writeLog(logpath,file,'not .fits ext') continue # file name prefix (night identifier) night_nm = fits_nm[:-9] # Extract fits header values using pyfits header = extractHeader(fits_file) header['filename'] = fits_nm[:-5] # Either the ccd shutter is closed or some fits keywords # are missing. We skip this image. if not header['imaging']: writeLog(logpath,file,'Shutter Closed') continue # Skip spectroscopy images if (header['mode'] == 'SPECTROSCOPY'): writeLog(logpath,file,'Spectroscopy') continue # Execute an external fortran program that will return
for file in sys.argv[1:]: # If file does not exist, do not go further if not os.path.exists(file): WriteLog(file,'does not exist') continue # Only proceed if filename matches standard NOT format if not re.search(r'[A-Z]{2}[a-z]{2}[0-9]{6}', file): WriteLog(file,'filename does not match standard format') continue # Extract fits header values using pyfits # header = [<instrument>,<tbl_prefix>,primaryHDU,ext1HDU,ext2HDU,..,extNHDU] try: header = extractHeader(str(file)) except HeaderException, e: WriteLog(file,str(e)) continue # Pyfits do not recognize this a a fits file. Make a note and skip file. if header[0] == 'ERROR': WriteLog(file,'not a valid fits file') continue instrument = header[0] tbl_prefix = header[1] # filename (stripped of .fits extension) fits_nm = str(os.path.basename(file))[:-5]
# file path + file name fits_file = str(file) # filename fits_nm = str(os.path.basename(file)) # Check file extention corresponds to a fits file if not (fits_nm[-5:] == '.fits'): writeLog(logpath, file, 'not .fits ext') continue # file name prefix (night identifier) night_nm = fits_nm[:-9] # Extract fits header values using pyfits header = extractHeader(fits_file) header['filename'] = fits_nm[:-5] # Either the ccd shutter is closed or some fits keywords # are missing. We skip this image. if not header['imaging']: writeLog(logpath, file, 'Shutter Closed') continue # Skip spectroscopy images if (header['mode'] == 'SPECTROSCOPY'): writeLog(logpath, file, 'Spectroscopy') continue # Execute an external fortran program that will return # the string 'day' or 'night' depending whether the