def get_wl_identy(filename, coordlistname): iraf.onedspec() print 'run get_wl_identy identify...' print 'the input file is ' + filename print 'implot %s' % filename iraf.implot(image = filename) print 'run identify...' iraf.identify(images = filename , section = 'middle line', database = 'database' , coordlist = coordlistname, units = '', nsum = 10 , match = -3.0, maxfeatures = 50, zwidth = 100.0 , ftype = 'emission', fwidth = 18.0, cradius = 5.0 , threshold = 0.0, minsep = 2.0, function = 'spline3' , order = 1, sample = '*', niterate = 0 , low_reject = 3.0, high_reject = 3.0, grow = 0.0 , autowrite = False, graphics = 'stdgraph', cursor = '' , crval = '', cdelt = '')
def shift_align(filenames): if len(filenames) == 0: return '' fno = 1 for i in filenames: iraf.display(image = i, frame = fno) fno += 1 print 'The fits number is', len(filenames) if len(filenames) > 1: print 'Please check the fits pic and determine whether you need to shift the pic' print 'imexamine...' iraf.imexamine() valget = raw_input('Are you need to shift the fits picture?(y)') if valget == 'y' or valget == 'Y': nfilenames = [filenames[0]] for i2 in filenames: print 'implot ' + i2 iraf.implot(image = i2) fno = 1 iraf.immatch() while fno < len(filenames): valget = raw_input('Please input ' + filenames[fno] + ' xshift, yshift: ') valget = valget.split() fxshift = float(valget[0]) fyshift = float(valget[1]) ims_key = 1 if os.path.isfile('shift_' + filenames[fno]): valget = raw_input('shift_' + filenames[fno] + ' is already exist, are you want to delete the old one?(y):') if valget == 'y' or valget == 'Y' or valget == 'yes' or valget == 'Yes': os.remove('shift_' + filenames[fno]) else: ims_key = 0 if ims_key == 1: print 'what happend????????????????' iraf.imshift(input = filenames[fno] , output = 'shift_' + filenames[fno], xshift = fxshift , yshift = fyshift, shifts_file = '', interp_type = 'linear' , boundary_type = 'nearest', constant = 0.0) nfilenames.append('shift_' + filenames[fno]) fno += 1 for i2 in xrange(len(nfilenames)): iraf.display(image = nfilenames[i2], frame = i2 + 1) iraf.imexamine() return cos_clear(nfilenames) return cos_clear(filenames)