def run_cosmic_rays(metadata, prev_suffix, curr_suffix, ns=False, multithread=False): from lacosmic import lacos_wifes # now run ONLY ON SCIENCE TARGETS AND STANDARDS sci_obs_list = get_sci_obs_list(metadata) sky_obs_list = get_sky_obs_list(metadata) std_obs_list = get_std_obs_list(metadata) for fn in sci_obs_list+sky_obs_list: in_fn = '%s%s.p%s.fits' % (out_dir, fn, prev_suffix) out_fn = '%s%s.p%s.fits' % (out_dir, fn, curr_suffix) print 'Cleaning cosmics in %s' % in_fn.split('/')[-1] # skip files which are already done #if os.path.isfile(out_fn): # continue if skip_done and os.path.isfile(out_fn): continue lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) if ns: in_fn = '%s%s.s%s.fits' % (out_dir, fn, prev_suffix) out_fn = '%s%s.s%s.fits' % (out_dir, fn, curr_suffix) print 'Cleaning cosmics in %s' % in_fn.split('/')[-1] lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) gc.collect() for fn in std_obs_list: in_fn = '%s%s.p%s.fits' % (out_dir, fn, prev_suffix) out_fn = '%s%s.p%s.fits' % (out_dir, fn, curr_suffix) if skip_done and os.path.isfile(out_fn): continue print 'Cleaning cosmics in %s' % in_fn.split('/')[-1] #lacos_wifes(in_fn, out_fn, niter=1, sig_frac=2.0) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) if ns: in_fn = '%s%s.s%s.fits' % (out_dir, fn, prev_suffix) out_fn = '%s%s.s%s.fits' % (out_dir, fn, curr_suffix) print 'Cleaning cosmics in %s' % in_fn.split('/')[-1] #lacos_wifes(in_fn, out_fn, niter=1, sig_frac=2.0) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) gc.collect() return
def run_cosmic_rays(metadata, prev_suffix, curr_suffix, ns=False, multithread=False): from lacosmic import lacos_wifes # now run ONLY ON SCIENCE TARGETS AND STANDARDS sci_obs_list = get_sci_obs_list(metadata) sky_obs_list = get_sky_obs_list(metadata) std_obs_list = get_std_obs_list(metadata) for fn in sci_obs_list+sky_obs_list: in_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, curr_suffix)) print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) # skip files which are already done #if os.path.isfile(out_fn): # continue if skip_done and os.path.isfile(out_fn): continue lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) if ns: in_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, curr_suffix)) print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) gc.collect() for fn in std_obs_list: in_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, curr_suffix)) if skip_done and os.path.isfile(out_fn): continue print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) #lacos_wifes(in_fn, out_fn, niter=1, sig_frac=2.0) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) if ns: in_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, curr_suffix)) print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) #lacos_wifes(in_fn, out_fn, niter=1, sig_frac=2.0) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) gc.collect() return
from lacosmic import lacos_wifes # now run ONLY ON SCIENCE TARGETS AND STANDARDS sci_obs_list = get_sci_obs_list(metadata) sky_obs_list = get_sky_obs_list(metadata) std_obs_list = get_std_obs_list(metadata) for fn in sci_obs_list+sky_obs_list: in_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, curr_suffix)) print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) # skip files which are already done #if os.path.isfile(out_fn): # continue if skip_done and os.path.isfile(out_fn): continue lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) if ns: in_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, curr_suffix)) print('Cleaning cosmics in %s' % in_fn.split('/')[-1]) lacos_wifes(in_fn, out_fn, wsol_fn=wsol_out_fn, niter=3, sig_clip=10.0, obj_lim=10.0, sig_frac=0.2, multithread=multithread) gc.collect() for fn in std_obs_list: in_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, prev_suffix)) out_fn = os.path.join(out_dir, '%s.p%s.fits' % (fn, curr_suffix)) if skip_done and os.path.isfile(out_fn): continue print('Cleaning cosmics in %s' % in_fn.split('/')[-1])