Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
AstroDrizzle(files, output=root, clean=True, final_scale=None, final_pixfrac=1, context=False, final_bits=576, preserve=False, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, build=False, final_wht_type='IVM')
#cat = make_drz_catalog(root=root, threshold=2)
cat = make_SEP_catalog(root=root, threshold=2)
if False:
files0 = ['icou09fvq_flt.fits', 'icou09fyq_flt.fits', 'icou09gpq_flt.fits',
'icou09h3q_flt.fits']
files1 = ['icou10emq_flt.fits', 'icou10eqq_flt.fits', 'icou10euq_flt.fits',
'icou10frq_flt.fits']
all_files = list(np.append(files0, files1))
AstroDrizzle(all_files, output='total', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_rot=0, build=False, final_wht_type='IVM')
AstroDrizzle(files0, output='group0', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
AstroDrizzle(files1, output='group1', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
im0 = pyfits.open('group0_drz_sci.fits')
im1 = pyfits.open('group1_drz_sci.fits')
imt = pyfits.open('total_drz_sci.fits')
# Overwrite image
for i in range(N):
images[i].flush()
if drizzle:
files = [flt.filename() for flt in images]
bits = 576
if root.startswith('par'):
pixfrac=1.0
else:
pixfrac=0.8
AstroDrizzle(files, output=root,
clean=True, final_pixfrac=pixfrac, context=False,
resetbits=0, final_bits=bits, driz_sep_bits=bits,
preserve=False, driz_separate=False,
driz_sep_wcs=False, median=False, blot=False,
driz_cr=False, driz_cr_corr=False, build=False,
final_wht_type='IVM')
clean_drizzle(root)
#cat = make_drz_catalog(root=root)
cat = make_SEP_catalog(root=root)
elif ref_catalog == 'WISE':
align_mag_limits = [15,20]
else:
ref_catalog = 'USER'
print('{0}: First Drizzle'.format(direct['product']))
### Clean up
for ext in ['.fits', '.log']:
file = '{0}_wcs.{1}'.format(direct['product'], ext)
if os.path.exists(file):
os.remove(file)
### First drizzle
if len(direct['files']) > 1:
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, context=False, preserve=False,
skysub=True, driz_separate=True, driz_sep_wcs=True,
median=True, blot=True, driz_cr=True,
driz_cr_snr=driz_cr_snr, driz_cr_scale=driz_cr_scale,
driz_cr_corr=False, driz_combine=True,
final_bits=bits, coeffs=True, build=False,
final_wht_type='IVM', **drizzle_params)
else:
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, final_scale=None, final_pixfrac=1,
context=False, final_bits=bits, preserve=False,
driz_separate=False, driz_sep_wcs=False,
median=False, blot=False, driz_cr=False,
driz_cr_corr=False, driz_combine=True,
build=False, final_wht_type='IVM', **drizzle_params)
if 'reference' in group:
AstroDrizzle(group['files'], output=group['product'],
clean=True, context=context, preserve=False,
skysub=skysub, skyuser=skyuser, skymethod=skymethod,
driz_separate=False, driz_sep_wcs=False,
median=False, blot=False, driz_cr=False,
driz_cr_corr=False, driz_combine=True,
final_bits=bits, coeffs=True, build=False,
final_wht_type=final_wht_type,
final_wt_scl=final_wt_scl,
final_pixfrac=pixfrac,
final_wcs=True, final_refimage=group['reference'],
resetbits=0)
else:
AstroDrizzle(group['files'], output=group['product'],
clean=True, context=context, preserve=False,
skysub=skysub, skyuser=skyuser, skymethod=skymethod,
driz_separate=False, driz_sep_wcs=False,
median=False, blot=False, driz_cr=False,
driz_cr_corr=False, driz_combine=True,
final_bits=bits, coeffs=True, build=False,
final_wht_type=final_wht_type,
final_wt_scl=final_wt_scl,
final_pixfrac=pixfrac,
final_wcs=final_wcs, final_rot=final_rot,
final_scale=scale,
final_ra=final_ra, final_dec=final_dec,
final_outnx=final_outnx, final_outny=final_outny,
resetbits=0)
clean_drizzle(group['product'])
### Redrizzle with new background subtraction
if isACS:
skyfile=''
else:
skyfile = '/tmp/{0}.skyfile'.format(grism['product'])
fp = open(skyfile,'w')
fp.writelines(['{0} 0.0\n'.format(f) for f in grism['files']])
fp.close()
if 'par' in grism['product']:
pixfrac=1.0
else:
pixfrac=0.8
AstroDrizzle(grism['files'], output=grism['product'], clean=True,
context=isACS, preserve=False, skysub=True, skyfile=skyfile,
driz_separate=gris_cr_corr, driz_sep_wcs=gris_cr_corr, median=gris_cr_corr,
blot=gris_cr_corr, driz_cr=gris_cr_corr, driz_cr_corr=gris_cr_corr,
driz_cr_snr=driz_cr_snr, driz_cr_scale=driz_cr_scale,
driz_combine=True, driz_sep_bits=bits, final_bits=bits,
coeffs=True, resetbits=4096, final_pixfrac=pixfrac,
build=False, final_wht_type='IVM')
clean_drizzle(grism['product'])
### Add direct filter to grism FLT headers
set_grism_dfilter(direct, grism)
return True
else:
all_groups = exposure_groups
for group in all_groups:
if len(group['files']) == 0:
continue
isACS = '_flc' in group['files'][0]
if bits is None:
if isACS:
bits = 64+32
else:
bits = 576
if 'reference' in group:
AstroDrizzle(group['files'], output=group['product'],
clean=True, context=context, preserve=False,
skysub=skysub, skyuser=skyuser, skymethod=skymethod,
driz_separate=False, driz_sep_wcs=False,
median=False, blot=False, driz_cr=False,
driz_cr_corr=False, driz_combine=True,
final_bits=bits, coeffs=True, build=False,
final_wht_type=final_wht_type,
final_wt_scl=final_wt_scl,
final_pixfrac=pixfrac,
final_wcs=True, final_refimage=group['reference'],
resetbits=0)
else:
AstroDrizzle(group['files'], output=group['product'],
clean=True, context=context, preserve=False,
skysub=skysub, skyuser=skyuser, skymethod=skymethod,
driz_separate=False, driz_sep_wcs=False,
xsh=out_shift[0], ysh=out_shift[1],
rot=out_rot, scale=out_scale,
wcsname=ref_catalog, force=True,
reusename=True, verbose=True,
sciext='SCI')
### Bug in astrodrizzle? Dies if the FLT files don't have MJD-OBS
### keywords
im = pyfits.open(file, mode='update')
im[0].header['MJD-OBS'] = im[0].header['EXPSTART']
im.flush()
### Second drizzle with aligned wcs, refined CR-rejection params
### tuned for WFC3/IR
if len(direct['files']) == 1:
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, final_pixfrac=0.8, context=False,
resetbits=4096, final_bits=bits, driz_sep_bits=bits,
preserve=False, driz_cr_snr=driz_cr_snr,
driz_cr_scale=driz_cr_scale, driz_separate=False,
driz_sep_wcs=False, median=False, blot=False,
driz_cr=False, driz_cr_corr=False,
build=False, final_wht_type='IVM', **drizzle_params)
else:
if 'par' in direct['product']:
pixfrac=1.0
else:
pixfrac=0.8
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, final_pixfrac=pixfrac, context=isACS,
resetbits=4096, final_bits=bits, driz_sep_bits=bits,
AstroDrizzle(files, output=root, clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, build=False, final_wht_type='IVM')
else:
AstroDrizzle(files, output=root, clean=True, final_scale=None, final_pixfrac=1, context=False, final_bits=576, preserve=False, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, build=False, final_wht_type='IVM')
#cat = make_drz_catalog(root=root, threshold=2)
cat = make_SEP_catalog(root=root, threshold=2)
if False:
files0 = ['icou09fvq_flt.fits', 'icou09fyq_flt.fits', 'icou09gpq_flt.fits',
'icou09h3q_flt.fits']
files1 = ['icou10emq_flt.fits', 'icou10eqq_flt.fits', 'icou10euq_flt.fits',
'icou10frq_flt.fits']
all_files = list(np.append(files0, files1))
AstroDrizzle(all_files, output='total', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_rot=0, build=False, final_wht_type='IVM')
AstroDrizzle(files0, output='group0', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
AstroDrizzle(files1, output='group1', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
im0 = pyfits.open('group0_drz_sci.fits')
im1 = pyfits.open('group1_drz_sci.fits')
imt = pyfits.open('total_drz_sci.fits')
#cat = make_drz_catalog(root=root, threshold=2)
cat = make_SEP_catalog(root=root, threshold=2)
if False:
files0 = ['icou09fvq_flt.fits', 'icou09fyq_flt.fits', 'icou09gpq_flt.fits',
'icou09h3q_flt.fits']
files1 = ['icou10emq_flt.fits', 'icou10eqq_flt.fits', 'icou10euq_flt.fits',
'icou10frq_flt.fits']
all_files = list(np.append(files0, files1))
AstroDrizzle(all_files, output='total', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_rot=0, build=False, final_wht_type='IVM')
AstroDrizzle(files0, output='group0', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
AstroDrizzle(files1, output='group1', clean=True, context=False, preserve=False, skysub=True, driz_separate=False, driz_sep_wcs=False, median=False, blot=False, driz_cr=False, driz_cr_corr=False, driz_combine=True, final_bits=576, coeffs=True, resetbits=0, final_refimage='total_drz_sci.fits', build=False, final_wht_type='IVM')
im0 = pyfits.open('group0_drz_sci.fits')
im1 = pyfits.open('group1_drz_sci.fits')
imt = pyfits.open('total_drz_sci.fits')
if len(direct['files']) == 1:
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, final_pixfrac=0.8, context=False,
resetbits=4096, final_bits=bits, driz_sep_bits=bits,
preserve=False, driz_cr_snr=driz_cr_snr,
driz_cr_scale=driz_cr_scale, driz_separate=False,
driz_sep_wcs=False, median=False, blot=False,
driz_cr=False, driz_cr_corr=False,
build=False, final_wht_type='IVM', **drizzle_params)
else:
if 'par' in direct['product']:
pixfrac=1.0
else:
pixfrac=0.8
AstroDrizzle(direct['files'], output=direct['product'],
clean=True, final_pixfrac=pixfrac, context=isACS,
resetbits=4096, final_bits=bits, driz_sep_bits=bits,
preserve=False, driz_cr_snr=driz_cr_snr,
driz_cr_scale=driz_cr_scale, build=False,
final_wht_type='IVM', **drizzle_params)
### Make DRZ catalog again with updated DRZWCS
clean_drizzle(direct['product'])
if isWFPC2:
thresh = 8
else:
thresh = 1.6
#cat = make_drz_catalog(root=direct['product'], threshold=thresh)
cat = make_SEP_catalog(root=direct['product'], threshold=thresh)