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