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])