def run_sky_sub(metadata, prev_suffix, curr_suffix, ns=False):
if ns:
run_sky_sub_ns(metadata, prev_suffix, curr_suffix)
else:
# subtract sky frames from science objects
for obs in metadata['sci']:
if len(obs['sky']) > 0:
sky_fn = obs['sky'][0]
sky_proc_fn = os.path.join(out_dir, '%s.p%s.fits' % (sky_fn, prev_suffix))
for fn in obs['sci']:
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('Subtracting sky frame for %s' % in_fn.split('/')[-1])
# subtract scaled sky frame!
pywifes.scaled_imarith_mef(
in_fn, '-', sky_proc_fn, out_fn,
scale='exptime')
else:
for fn in obs['sci']:
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('Copying image %s' % in_fn.split('/')[-1])
# subtract scaled sky frame!
pywifes.imcopy(in_fn, out_fn)
# copy stdstar frames
std_obs_list = get_std_obs_list(metadata)
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('Copying standard star image %s' % in_fn.split('/')[-1])
pywifes.imcopy(in_fn, out_fn)
return
def run_sky_sub(metadata, prev_suffix, curr_suffix, ns=False):
if ns:
run_sky_sub_ns(metadata, prev_suffix, curr_suffix)
else:
# subtract sky frames from science objects
for obs in metadata['sci']:
if len(obs['sky']) > 0:
sky_fn = obs['sky'][0]
sky_proc_fn = '%s%s.p%s.fits' % (out_dir, sky_fn, prev_suffix)
for fn in obs['sci']:
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 'Subtracting sky frame for %s' % in_fn.split('/')[-1]
# subtract scaled sky frame!
pywifes.scaled_imarith_mef(
in_fn, '-', sky_proc_fn, out_fn,
scale='exptime')
else:
for fn in obs['sci']:
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 'Copying image %s' % in_fn.split('/')[-1]
# subtract scaled sky frame!
pywifes.imcopy(in_fn, out_fn)
# copy stdstar frames
std_obs_list = get_std_obs_list(metadata)
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 'Copying standard star image %s' % in_fn.split('/')[-1]
pywifes.imcopy(in_fn, out_fn)
return
def run_sky_sub_ns(metadata, prev_suffix, curr_suffix):
sci_obs_list = get_sci_obs_list(metadata)
std_obs_list = get_std_obs_list(metadata)
ns_proc_list = sci_obs_list + std_obs_list
for fn in ns_proc_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)
sky_fn = '%s%s.s%s.fits' % (out_dir, fn, prev_suffix)
print 'Subtracting N+S sky frame for %s' % in_fn.split('/')[-1]
pywifes.scaled_imarith_mef(in_fn, '-', sky_fn, out_fn, scale='exptime')
return
def run_sky_sub_ns(metadata, prev_suffix, curr_suffix):
sci_obs_list = get_sci_obs_list(metadata)
std_obs_list = get_std_obs_list(metadata)
ns_proc_list = sci_obs_list+std_obs_list
for fn in ns_proc_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)
sky_fn = '%s%s.s%s.fits' % (out_dir, fn, prev_suffix)
print 'Subtracting N+S sky frame for %s' % in_fn.split('/')[-1]
pywifes.scaled_imarith_mef(
in_fn, '-', sky_fn, out_fn,
scale='exptime')
return
gc.collect()
return
#------------------------------------------------------
# Sky subtraction!
def run_sky_sub_ns(metadata, prev_suffix, curr_suffix):
sci_obs_list = get_sci_obs_list(metadata)
std_obs_list = get_std_obs_list(metadata)
ns_proc_list = sci_obs_list+std_obs_list
for fn in ns_proc_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))
sky_fn = os.path.join(out_dir, '%s.s%s.fits' % (fn, prev_suffix))
print('Subtracting N+S sky frame for %s' % in_fn.split('/')[-1])
pywifes.scaled_imarith_mef(
in_fn, '-', sky_fn, out_fn,
scale='exptime')
return
def run_sky_sub(metadata, prev_suffix, curr_suffix, ns=False):
if ns:
run_sky_sub_ns(metadata, prev_suffix, curr_suffix)
else:
# subtract sky frames from science objects
for obs in metadata['sci']:
if len(obs['sky']) > 0:
sky_fn = obs['sky'][0]
sky_proc_fn = os.path.join(out_dir, '%s.p%s.fits' % (sky_fn, prev_suffix))
for fn in obs['sci']:
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))
