This needs to be run from the directory "reduce_acs.py" is in. Also requires ureka to be activated. """ import os import shutil import reduce_acs ref_scale = 0.03 ref_filter = 'F606W' refs_dir = 'drizzled_refs' asnd, _ = reduce_acs.find_target_associations('/Volumes/ErikDroboData/data/ACS/pisc_dwarfs/', {'HI22': 'PiscA', 'HI23': 'PiscB'}) for i, (objnm, asns) in enumerate(asnd.items()): print('Making ref for', objnm, '- object', i+1, 'of', len(asnd)) refasn = [asn for asn in asns if asn.filter == ref_filter] if len(refasn) == 0: raise ValueError('Found no association with filter {0} for obj {1}'.format(ref_filter, objnm)) elif len(refasn) > 1: print('Found multiple associations with filter', ref_filter, 'using first one') refasn = refasn[:1] outnm = objnm + '_ref' workingdir = os.path.join(refs_dir, objnm) reduce_acs.do_drizzle(refasn, outnm, workingdir, final_scale=ref_scale, final_wcs=True) shutil.move(os.path.join(workingdir, outnm) + '_drc_sci.fits', os.path.join(refs_dir, outnm) + '_drc_sci.fits') shutil.rmtree(workingdir)
refs_dir = 'drizzled_refs' pixfrac = 0.6 build = False dirprefix = 'drizzled_build_' if build else 'drizzled_' final_wht_type = None asnd, _ = reduce_acs.find_target_associations('raw', {'HI22': 'PiscA', 'HI23': 'PiscB'}) for objnm, asns in asnd.items(): reffn = os.path.abspath(os.path.join(refs_dir, objnm + '_ref_drc_sci.fits')) for i, asn in enumerate(asns): print('Starting', i+1, 'of', len(asns), 'associations for ', objnm) if final_wht_type is None: whtending = '' else: whtending = '_' + final_wht_type reduce_acs.do_drizzle([asn], objnm + '_' + asn.filter + whtending, dirprefix, final_refimage=reffn, final_wcs=True, final_pixfrac=pixfrac, final_wht_type=final_wht_type, build=build) print('Doing combination of all associations for', objnm) reduce_acs.do_drizzle(asns, objnm + '_allfilters_ivm', dirprefix, final_refimage=reffn, final_wcs=True, final_pixfrac=pixfrac, final_wht_type='IVM', build=build) reduce_acs.do_drizzle(asns, objnm + '_allfilters_err', dirprefix, final_refimage=reffn, final_wcs=True, final_pixfrac=pixfrac, final_wht_type='ERR', build=build)