예제 #1
0
import astromatic as ast
import sys

if len(sys.argv)<3:
    print ''
    print 'ERROR: run_sext_scam_2.py requires 2 parameters:'
    print '    1. Input science fits filename'
    print '    2. Zero point'
    print ''
    print 'Example:  python run_sext_scam_2.py 0435_scam_mar14_r.fits 27.685'
    print ''
    exit()

infile = sys.argv[1]
zeropoint = float(sys.argv[2])
whtfile = infile.replace('.fits','_wht.fits')
regfile = infile.replace('.fits','.reg')
outcat  = infile.replace('.fits','.cat')
configfile = 'sext_scam_final.config'
print 'Running SExtractor on %s' % infile

ast.make_cat_suprimecam(infile,outcat=outcat,regfile=regfile,
                        configfile=configfile,whtfile=whtfile,
                        zeropt=zeropoint)

예제 #2
0
import astromatic as ast
import sys

infile = sys.argv[1]
whtfile = infile.replace(".fits", "_wht_init.fits")
regfile = infile.replace(".fits", ".reg")
outcat = infile.replace(".fits", ".cat")
print "Running SExtractor on %s" % infile

ast.make_cat_suprimecam(
    infile, outcat=outcat, regfile=regfile, whtfile=whtfile, configfile="sext_scam_astrom.config", satur=60000.0
)
예제 #3
0
    hdr = pf.getheader(scifile)
except:
    print ''
    print 'ERROR: Could not open header for fits file %s' % scifile
    print ''
    exit()
try:
    zeropoint = hdr['magzpt']
except:
    print ''
    print 'ERROR: Could not get zeropoint in header for fits file %s' % scifile
    print ''
    del hdr
    exit()
del hdr

""" Summarize input information """
print ''
print 'SExtractor input parameters'
print '-------------------------------------------------------'
print 'Science file:   %s' % scifile
print 'Weight file:    %s' % whtfile
print 'Output catalog: %s' % catfile
print 'Regions file:   %s' % regfile
print 'Config file:    %s' % configfile
print 'Zeropoint       %6.3f' % zeropoint
print ''

""" Actually run SExtractor """
astrom.make_cat_suprimecam(scifile,catfile,regfile,configfile=configfile,whtfile=whtfile,zeropt=zeropoint)