Beispiel #1
0
    t3refroot = 'ff260'
    t3refcat = '%s.cat' % t3refroot
    t3ccdbfile = 'ccmap_%s_as_ref.txt' % t3refroot
""" 
Do the nonlinearity correction and Create the sky file as the first pass 
"""
if redpass == 'make_sky':
    niri.calib_1(sci_frames, rawroot, outroot, bpmfile, rawdir)
""" Do the final calibration and make the SExtractor catalogs """
if redpass == 'final_calib':
    print ""
    skyfile = '%s_sky.fits' % outroot
    niri.reduce_sci(sci_frames, outroot, skyfile, flatfile)
    ##niri.niri_coadd(outroot) # Don't use this any more
    if skybpm is not None:
        niri.niri_bpm_from_sky(sci_frames, skybpm)
""" Make the SExtractor catalogs """
if redpass == 'make_cats':
    ff_files = []
    for i in sci_frames:
        ff_files.append('ff%d.fits' % i)
    niri.niri_sextractor(ff_files, tilebpm, fixpix=True, catformat='ascii')
"""
Do the tile-based pass of the astrometry.
"""

if redpass == 'astrom_tile':
    if tile == "Tiles2and3":
        swarplist = '%s_swarp.in' % t2root
        #niri.ccmap_tile(t2_frames,t2refcat,t2ccdbfile,tilebpm,swarplist)
        swarpfile = '@%s' % swarplist
Beispiel #2
0
""" 
Do the nonlinearity correction and create the sky file as the first pass 
"""
if redpass=='calib_1':
   niri.calib_1(sci_frames,rawroot,outroot,bpmfile,rawdir)

""" Calibrate the science frames """
if redpass=='calib_2':
   print ""
   skyfile  = '%s_sky.fits' % outroot
   niri.reduce_sci(sci_frames,outroot,skyfile,flatfile)

""" Make the bad pixel mask from the calibrated science frames """
if redpass=='calib_3':
   if skybpm is not None:
      niri.niri_bpm_from_sky(sci_frames,skybpm,outsky='sky_from_ff.fits')
   else:
      print ""
      print "*** No bad pixel mask requested for this set of input files ***"
      print ""

""" 
Prepare the files for running SExtractor
 The input files will be called ff*fits
 The output files will be called fp*fits
"""
if redpass=='calib_4':
   niri.split_and_fix_ff(sci_frames,badpixfile=tilebpm,fixpix=True)

"""
For these Ks-band images, there needs to be a second sky flat correction
Beispiel #3
0
   for i in sci_frames:
      fullnames.append('%s_%d' % (obsdate,i))

""" 
Do the nonlinearity correction and Create the sky file as the first pass 
"""
if redpass=='make_sky':
   niri.calib_1(sci_frames,rawroot,outroot,bpmfile,rawdir,obsdate)

""" Do the final calibration and make the SExtractor catalogs """
if redpass=='final_calib':
   print ""
   skyfile  = '%s_sky.fits' % outroot
   niri.reduce_sci(fullnames,outroot,skyfile,flatfile)
   if skybpm is not None:
      niri.niri_bpm_from_sky(fullnames,skybpm)

""" Make the SExtractor catalogs """
if redpass=='make_cats':
   ff_files = []
   for i in fullnames:
      ff_files.append('ff%s.fits' % i)
   niri.niri_sextractor(ff_files,tilebpm,fixpix=True,catformat='ascii')

"""
Redo the astrometry of all the individual input files onto the full grid
now that we have created it, and then do the final coadd.
"""

if redpass == 'final_coadd':
   """ Set up """