def reduce_beehive():
util.mkdir(out_dir)
## Loop through all the different data sets and reduce them.
#for key in dict_suffix.keys():
#for key in ['LS_3wfs_r2', 'LS_5wfs_r2', 'open_r2']:
for key in ['open_r2']:
img = dict_images[key]
suf = dict_suffix[key]
sky = 'beehive_sky2.fits'
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
print(' Sky: ', sky)
img_files = [
data_dir + 'sta{img:03d}{suf:s}.fits'.format(img=ii, suf=suf)
for ii in img
]
scn_files = [
data_dir + 'sta{img:03d}{suf:s}_scan.fits'.format(img=ii, suf=suf)
for ii in img
]
reduce_STA.treat_overscan(img_files)
#reduce_STA.treat_overscan_working(img_files) #BUG
redu.clean_images(scn_files,
out_dir,
rebin=2,
sky_frame=sky_dir + sky,
flat_frame=calib_dir + "domeflat.fits")
return
def reduce_FLD2():
util.mkdir(out_dir)
# Open Loop
img_files = [data_dir + 'obj{0:03d}_o.fits'.format(ii) for ii in fnum_o]
reduce_STA.treat_overscan(img_files)
scan_img_files = [
data_dir + 'obj{0:03d}_o_scan.fits'.format(ii) for ii in fnum_o
]
reduce_fli.clean_images(scan_img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - threeWFS_LS
img_files = [
data_dir + 'obj{0:03d}threeWFS_LS_c.fits'.format(ii) for ii in fnum_c
]
reduce_STA.treat_overscan(img_files)
scan_img_files = [
data_dir + 'obj{0:03d}threeWFS_LS_c_scan.fits'.format(ii)
for ii in fnum_c
]
reduce_fli.clean_images(scan_img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + "flat.fits")
return
def reduce_pleiades():
##########
# Open Loop
##########
data_dir = root_dir + 'Pleiades/'
sky_file = root_dir + 'reduce/sky/flat.fits'
flat_file = root_dir + 'reduce/calib/sky.fits'
out_dir = root_dir + 'reduce/pleiades/'
fnum = np.arange(57, 67)
img_files = ['{0:s}/obj{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1,
sky_frame=sky_file,
flat_frame=flat_file)
##########
# Closed Loop
##########
data_dir = root_dir + 'Pleiades/'
sky_file = root_dir + 'reduce/sky/flat.fits'
flat_file = root_dir + 'reduce/calib/sky.fits'
out_dir = root_dir + 'reduce/pleiades/'
fnum = np.arange(47, 57)
img_files = ['{0:s}/obj{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1,
sky_frame=sky_file,
flat_frame=flat_file)
return
def reduce_beehive():
util.mkdir(out_dir)
# Loop through all the different data sets and reduce them.
for key in dict_suffix.keys():
img = dict_images[key]
suf = dict_suffix[key]
sky = dict_skies[key]
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
print(' Sky: ', sky)
img_files = [
data_dir + 'sta{img:03d}{suf:s}.fits'.format(img=ii, suf=suf)
for ii in img
]
scn_files = [
data_dir + 'sta{img:03d}{suf:s}_scan.fits'.format(img=ii, suf=suf)
for ii in img
]
reduce_STA.treat_overscan(img_files)
reduce_fli.clean_images(scn_files,
out_dir,
rebin=1,
sky_frame=sky_dir + sky,
flat_frame=flat_dir + "flat.fits")
return
def reduce_fld2():
util.mkdir(out_dir)
## Loop through all the different data sets
for key in dict_suffix.keys():
#for key in ['doczskycl']:
img = dict_images[key]
suf = dict_suffix[key]
sky = dict_skies[key]
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
print(' Sky: ', sky)
img_files = [
data_dir + 'sta{img:03d}{suf:s}.fits'.format(img=ii, suf=suf)
for ii in img
]
scn_files = [
data_dir + 'sta{img:03d}{suf:s}_scan.fits'.format(img=ii, suf=suf)
for ii in img
]
reduce_STA.treat_overscan(img_files)
redu.clean_images(scn_files,
out_dir,
rebin=1,
sky_frame=sky_dir + sky,
flat_frame=calib_dir + "flat.fits") #,
# fix_bad_pixels=True, worry_about_edges=True)
return
def reduce_fld2():
util.mkdir(out_dir)
## Loop through all the different data sets and reduce them.
#for key in ['open_RIVB', 'LS_RIVB', 'docz_RIVB']: ## Single key setup
for key in dict_suffix.keys():
img = dict_images[key]
suf = dict_suffix[key]
filt = dict_filt[key]
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
print(' Filter: ', filt)
img_files = [data_dir + 'sta{img:03d}{suf:s}.fits'.format(img=ii, suf=suf) for ii in img]
scn_files = [data_dir + 'sta{img:03d}{suf:s}_scan.fits'.format(img=ii, suf=suf) for ii in img]
reduce_STA.treat_overscan(img_files)
redu.clean_images(scn_files, out_dir, rebin=1,
sky_frame=sky_dir + f"fld2_sky_{filt}.fits",
flat_frame=calib_dir + f"flat_{filt}.fits")#,
# fix_bad_pixels=True, worry_about_edges=True)
return
def reduce_pleiades_binned_ttf():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/flat/'
out_dir = root_dir + 'reduce/pleiades/'
os.chdir(data_dir)
fnum1 = [
143, 144, 145, 146, 147, 148, 151, 152, 153, 154, 157, 157, 158, 161,
162, 165, 166
]
fnum2 = [
169, 170, 173, 174, 175, 176, 181, 182, 187, 188, 189, 190, 202, 203,
208, 209, 214
]
fnum3 = [215, 220, 221, 224, 225, 228, 229]
fnum = fnum1 + fnum2 + fnum3
img_files = ['obj_ttf{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_195.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_pleiades_unbinned_open():
data_dir = '/Users/jlu/data/imaka/2016_11_18/Pleiades_E/open_loop/'
os.chdir(data_dir)
# img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(55, 57)]
img_files = ['obj_056.fits']
reduce_fli.clean_images(img_files, rebin=10)
return
def reduce_pleiades_binned_closed():
data_dir = '/Users/jlu/data/imaka/2016_11_18/Pleiades_E/closed_loop/'
os.chdir(data_dir)
# img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(92, 101)]
# img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(153, 163)]
img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(243, 247)]
img_files += ['obj_{0:03d}.fits'.format(ii) for ii in range(251, 254)]
reduce_fli.clean_images(img_files, rebin=2)
return
def reduce_pleiades_binned_open():
data_dir = '/Users/jlu/data/imaka/2016_11_18/Pleiades_E/open_loop/'
os.chdir(data_dir)
# img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(82, 91)]
# img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(163, 173)]
img_files = ['obj_{0:03d}.fits'.format(ii) for ii in range(247, 251)]
img_files += ['obj_{0:03d}.fits'.format(ii) for ii in range(255, 258)]
reduce_fli.clean_images(img_files, rebin=2)
return
def reduce_pleiades_binned_closed():
sky_dir = '/Users/jlu/data/imaka/2016_11_19/20161118/skies/binned/'
data_dir = '/Users/jlu/data/imaka/2016_11_19/20161118/Pleiades_E/closed_loop/'
os.chdir(data_dir)
fnum = [138, 141, 143, 147, 149, 151]
img_files = ['obj_{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits')
return
def reduce_FLD2():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'FLD2_2/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/FLD2_2/'
util.mkdir(out_dir)
# Open Loop
fnum = [4, 5, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45]
fnum += [51, 54, 57, 60, 63, 66, 69, 72, 75, 78, 81, 85, 88, 89]
fnum += [92, 95, 99, 103, 107, 110, 113, 116, 119, 122, 125, 128]
fnum += [
131, 134, 137, 140, 143, 146, 149, 152, 155, 158, 161, 164, 167, 170
]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop
fnum = [96, 100, 104]
img_files = [data_dir + 'obj{0:04d}_c.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed A
fnum = [6, 8, 11, 14, 17, 20, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53]
fnum += [56, 59, 62, 65, 68, 71, 74, 77, 80, 83, 87, 91, 94, 98, 102, 106]
fnum += [109, 112, 115, 118, 121, 124, 127, 130, 133, 136, 139, 142, 145]
fnum += [148, 151, 154, 157, 160, 163, 166, 169]
img_files = [data_dir + 'obj{0:04d}_cA.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed B
fnum = [10, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58]
fnum += [
61, 64, 67, 70, 73, 76, 79, 82, 86, 90, 93, 97, 101, 105, 108, 111
]
fnum += [
114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 144, 147, 150, 153
]
fnum += [156, 159, 162, 165, 168, 171]
img_files = [data_dir + 'obj{0:04d}_cB.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_orion():
util.mkdir(out_dir)
# Open Loop
img_files = [data_dir + 'obj{0:03d}_o.fits'.format(ii) for ii in fnum_o]
reduce_STA.treat_overscan(img_files)
scan_img_files = [data_dir + 'obj{0:03d}_o_scan.fits'.format(ii) for ii in fnum_o]
reduce_fli.clean_images(scan_img_files, out_dir, rebin=1, sky_frame=sky_dir + 'orion_sky.fits', flat_frame=flat_dir+"flat.fits")
# Closed Loop - 4W
img_files = [data_dir + 'obj{0:03d}LS4WFS_c.fits'.format(ii) for ii in fnum_c_4W]
reduce_STA.treat_overscan(img_files)
scan_img_files = [data_dir + 'obj{0:03d}LS4WFS_c_scan.fits'.format(ii) for ii in fnum_c_4W]
reduce_fli.clean_images(scan_img_files, out_dir, rebin=1, sky_frame=sky_dir + 'orion_sky.fits', flat_frame =flat_dir+"flat.fits")
# Closed Loop - B2
img_files = [data_dir + 'obj{0:03d}LS4WFS_B2_c.fits'.format(ii) for ii in fnum_c_B2]
reduce_STA.treat_overscan(img_files)
scan_img_files = [data_dir + 'obj{0:03d}LS4WFS_B2_c_scan.fits'.format(ii) for ii in fnum_c_B2]
reduce_fli.clean_images(scan_img_files, out_dir, rebin=1, sky_frame=sky_dir + 'orion_sky.fits', flat_frame =flat_dir+"flat.fits")
# Closed Loop - zc
img_files = [data_dir + 'obj{0:03d}LS4WFS_zc21_c.fits'.format(ii) for ii in fnum_c_zc]
reduce_STA.treat_overscan(img_files)
scan_img_files = [data_dir + 'obj{0:03d}LS4WFS_zc21_c_scan.fits'.format(ii) for ii in fnum_c_zc]
reduce_fli.clean_images(scan_img_files, out_dir, rebin=1, sky_frame=sky_dir + 'orion_sky.fits', flat_frame =flat_dir+"flat.fits")
return
def reduce_pleiades():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
#util.mkdir(out_dir)
# # Open Loop
# fnum = [59, 62, 65, 68, 70, 76, 82, 88]#, 94]
# img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame=flat_dir + 'flat.fits')
# # Closed Loop - a
# fnum = [58, 61, 64, 75, 81, 87]#, 93]
# img_files = [data_dir + 'obj{0:04d}_ca.fits'.format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame = flat_dir + 'flat.fits')
# # Closed Loop - b
# fnum = [71, 77, 83, 89]
# img_files = [data_dir + "obj{0:04d}_cb.fits".format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame = flat_dir + 'flat.fits')
# # Closed Loop - c
# fnum = [72, 78, 84, 90]
# img_files = [data_dir + "obj{0:04d}_cc.fits".format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame = flat_dir + 'flat.fits')
# # Closed Loop - d
# fnum = [73, 79, 85]#, 91]
# img_files = [data_dir + "obj{0:04d}_cd.fits".format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame = flat_dir + 'flat.fits')
# # Closed Loop - e
# fnum = [74, 80, 86]#, 92]
# img_files = [data_dir + "obj{0:04d}_ce.fits".format(ii) for ii in fnum]
# reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'pleiades_sky.fits', flat_frame = flat_dir + 'flat.fits')
# TT
fnum = [60, 63, 69]
img_files = [data_dir + "obj{0:04d}_tt.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_pleiadies_east():
data_dir = '/Users/jlu/Google Drive/Instruments/imaka/imaka (1)/Commissioning/2016-11 Observing/20161115/'
data_dir += 'Focal_plane_images/obs_11152016/'
dark_dir = data_dir + 'darks/'
plei_dir = data_dir + 'Pleides_E/'
dark_file = 'off_0S_0E.fits'
# img_files = ['off_0S_0E_1.fits', 'off_0S_0E_2.fits', 'off_10S_0E_1.fits', 'off_10S_0E_2.fits',
# 'off_10S_30E_1.fits', 'off_10S_30E_2.fits']
img_files = ['off_0S_0E_1.fits']
reduce_fli.clean_images(img_files)
return
def reduce_pleiades_binned_tt():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/flat/'
out_dir = root_dir + 'reduce/pleiades/'
os.chdir(data_dir)
fnum = [134, 135, 136, 137, 138, 139, 140]
img_files = ['obj_tt{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_078.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_FLD2():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'FLD2_2/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/FLD2_2/'
util.mkdir(out_dir)
# Open Loop
fnum = [20, 26, 32, 34, 38, 39, 40, 42, 45, 48, 51, 54, 57]
fnum += [60, 63, 66, 69, 72, 75, 78, 81, 84, 87, 90, 93, 95]
fnum += [99, 102, 105, 108, 111, 113, 116, 119, 122, 125, 128]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop
fnum = [19, 21, 22, 23, 24, 25, 27, 28, 29, 30, 31]
img_files = [data_dir + 'obj{0:04d}_c.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed A
fnum = [33, 41, 44, 47, 50, 53, 56, 59, 62, 65, 68]
fnum += [71, 74, 77, 80, 83, 86, 89, 92, 94]
fnum += [98, 101, 104, 107, 110, 112, 115, 118, 121, 124, 127]
img_files = [data_dir + 'obj{0:04d}_cA.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed B
fnum = [35, 36, 37, 43, 46, 49, 52, 55, 58, 61, 64]
fnum += [67, 70, 73, 76, 79, 82, 85, 88, 91, 96]
fnum += [100, 103, 106, 109, 114, 117, 120, 123, 126, 129]
img_files = [data_dir + 'obj{0:04d}_cB.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_FLD2():
util.mkdir(out_dir)
# Open Loop
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum_o]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=None,
flat_frame=None)
# Closed
img_files = [data_dir + 'obj{0:04d}_c.fits'.format(ii) for ii in fnum_c]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=None,
flat_frame=None)
return
def reduce_pleiades_binned_open():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
util.mkdir(out_dir)
os.chdir(data_dir)
fnum = [
9, 10, 13, 14, 17, 18, 21, 22, 28, 29, 32, 33, 36, 37, 40, 41, 46, 47,
50, 51, 54, 55, 58, 59
]
img_files = ['obj_o{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_023.fits',
flat_frame=flat_dir + 'flat.fits')
fnum = [64, 65, 68, 69, 72, 73, 76, 77, 82, 83, 86, 87, 90, 91, 94, 95, 100, 101, 104, 105, 108, 109, 112, \
113, 118, 119, 122, 123, 126, 127, 130, 131]
img_files = ['obj_o{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_025.fits',
flat_frame=flat_dir + 'flat.fits')
fnum = [141, 142, 149, 150]
img_files = ['obj_o{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_078.fits',
flat_frame=flat_dir + 'flat.fits')
fnum = [
179, 180, 185, 186, 193, 194, 200, 201, 206, 207, 212, 213, 218, 219
]
img_files = ['obj_o{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_195.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_fld2():
util.mkdir(out_dir)
## Loop through all the different data sets and reduce them.
#for key in ['set_name']: ## Single key setup
for key in dict_suffix.keys():
img = dict_images[key]
suf = dict_suffix[key]
sky = dict_skies[key]
bin_num = 'bin1' if 'bin1' in key else 'bin2' # key should contain bin info
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
print(' Sky: ', sky)
img_files = [
data_dir + 'sta{img:03d}{suf:s}.fits'.format(img=ii, suf=suf)
for ii in img
]
scn_files = [
data_dir + 'sta{img:03d}{suf:s}_scan.fits'.format(img=ii, suf=suf)
for ii in img
]
reduce_STA.treat_overscan(img_files)
redu.clean_images(scn_files,
out_dir,
rebin=1,
sky_frame=sky_dir + sky,
flat_frame=flat_dir + "flat_" + bin_num +
".fits") #,
# fix_bad_pixels=True, worry_about_edges=True)
return
def reduce_FLD2():
util.mkdir(out_dir)
# Open Loop
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum_o]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - threeWFS_LS
img_files = [
data_dir + 'obj{0:04d}_threeWFS_LS.fits'.format(ii)
for ii in fnum_threeWFS_LS
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - threeWFSLS_B2_c
img_files = [
data_dir + 'obj{0:04d}_threeWFSLS_B2_c.fits'.format(ii)
for ii in fnum_threeWFSLS_B2_c
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - threeWFSMean_B2_c
img_files = [
data_dir + 'obj{0:04d}_threeWFSMean_B2_c.fits'.format(ii)
for ii in fnum_threeWFSMean_B2_c
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
return
def reduce_FLD2():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'FLD2/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/FLD2/'
util.mkdir(out_dir)
# Open Loop
fnum = [63, 67, 71, 75, 80, 84, 88, 92, 96, 100, 104, 108, 112, 124]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop
fnum = [64, 68, 72, 76, 77, 81, 85, 89, 93, 97, 101, 105, 109, 113]
img_files = [data_dir + 'obj{0:04d}_c.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed A
fnum = [65, 69, 73, 78, 82, 86, 90, 94, 98, 102, 106, 110, 114]
img_files = [data_dir + 'obj{0:04d}_cA.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed B
fnum = [66, 70, 74, 79, 83, 87, 91, 95, 99, 103, 107, 111, 115]
img_files = [data_dir + 'obj{0:04d}_cB.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_pleiades_binned_closed():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/flat/'
out_dir = root_dir + 'reduce/pleiades/'
os.chdir(data_dir)
fnum1 = [
7, 8, 12, 15, 16, 19, 20, 26, 27, 30, 31, 34, 35, 38, 39, 44, 45, 48,
49, 52
]
fnum2 = [53, 56, 57, 62, 63]
fnum = fnum1 + fnum2
img_files = ['obj_c{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_025.fits',
flat_frame=flat_dir + 'flat.fits')
fnum1 = [
66, 67, 70, 71, 74, 75, 80, 81, 84, 85, 88, 89, 92, 93, 98, 99, 102,
103, 106
]
fnum2 = [107, 110, 111, 116, 117, 120, 121, 124, 125, 128, 129]
fnum = fnum1 + fnum2
img_files = ['obj_c{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_78.fits',
flat_frame=flat_dir + 'flat.fits')
fnum1 = [
155, 156, 163, 164, 167, 168, 171, 172, 177, 178, 183, 184, 191, 192,
198, 199
]
fnum2 = [204, 205, 210, 211, 216, 217, 222, 223, 226, 227, 230, 231]
fnum = fnum1 + fnum2
img_files = ['obj_c{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_195.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_FLD2():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'FLD2_2/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/FLD2_2/'
util.mkdir(out_dir)
# Open Loop
fnum = [14, 17, 29, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53, 56, 59]
fnum += [62, 65, 68, 71, 74, 77, 80, 83, 86, 89, 92]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed A
fnum = [12, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55]
fnum += [58, 61, 64, 67, 70, 73, 76, 79, 82, 85, 88, 91]
img_files = [data_dir + 'obj{0:04d}_cA.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed D
fnum = [18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63]
fnum += [66, 69, 72, 75, 78, 81, 84, 87, 90, 93]
img_files = [data_dir + 'obj{0:04d}_cB.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_2_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_pleiades():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
util.mkdir(out_dir)
# Open Loop
fnum1 = [
56, 57, 58, 65, 66, 67, 77, 78, 79, 92, 93, 94, 98, 99, 100, 104, 105,
106
]
fnum2 = [
113, 114, 115, 126, 127, 128, 135, 136, 137, 144, 145, 146, 156, 157,
158
]
fnum3 = [
165, 166, 167, 174, 175, 176, 186, 187, 188, 195, 196, 197, 204, 205,
206
]
fnum4 = [
216, 217, 218, 231, 232, 233, 240, 241, 242, 252, 253, 254, 261, 262,
263
]
fnum5 = [
270, 271, 272, 280, 281, 286, 287, 292, 293, 299, 300, 301, 302, 312,
313
]
fnum6 = [314, 321, 322, 323]
fnum = fnum1 + fnum2 + fnum3 + fnum4 + fnum5 + fnum6
img_files = [data_dir + 'obj_o{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# TTF Closed Loop
fnum1 = [
32, 33, 34, 38, 39, 40, 50, 51, 52, 59, 60, 61, 68, 69, 70, 80, 81, 82,
86
]
fnum2 = [
87, 88, 95, 96, 97, 107, 108, 109, 120, 121, 122, 129, 130, 131, 138,
139
]
fnum3 = [
140, 150, 151, 152, 159, 160, 161, 168, 169, 170, 180, 181, 182, 189,
190
]
fnum4 = [
191, 198, 199, 200, 210, 211, 212, 219, 220, 221, 225, 226, 227, 234,
235
]
fnum5 = [
236, 246, 247, 248, 255, 256, 257, 264, 265, 266, 276, 277, 282, 283,
288
]
fnum6 = [289, 294, 295, 306, 307, 308, 315, 316, 317]
fnum = fnum1 + fnum2 + fnum3 + fnum4 + fnum5 + fnum6
img_files = [data_dir + 'obj_ttf{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop
fnum1 = [
27, 28, 29, 30, 31, 35, 36, 37, 44, 45, 46, 47, 48, 49, 53, 54, 55, 62,
63, 64
]
fnum2 = [
74, 75, 76, 89, 90, 91, 101, 102, 103, 110, 111, 112, 123, 124, 125,
132, 133
]
fnum3 = [
134, 141, 142, 143, 153, 154, 155, 162, 163, 164, 171, 172, 173, 183,
184, 185
]
fnum4 = [
192, 193, 194, 201, 202, 203, 213, 214, 215, 222, 223, 224, 228, 229,
230, 237
]
fnum5 = [
238, 239, 249, 250, 251, 258, 259, 260, 267, 268, 269, 278, 279, 284,
285, 290
]
fnum6 = [291, 296, 297, 309, 310, 311, 318, 319, 320]
fnum = fnum1 + fnum2 + fnum3 + fnum4 + fnum5 + fnum6
img_files = [data_dir + 'obj_c{0:03d}.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_FLD2():
util.mkdir(out_dir)
# Open Loop - 30s
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum_o_30]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 3 WFS Small - 30s
img_files = [
data_dir + 'obj{0:04d}_threewfs_small_c.fits'.format(ii)
for ii in fnum_c_3S_30
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 3 WFS Large - 30s
img_files = [
data_dir + 'obj{0:04d}_threeWFS_big_c.fits'.format(ii)
for ii in fnum_c_3L_30
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 4 WFS - 30s
img_files = [
data_dir + 'obj{0:04d}_fourWFS_c.fits'.format(ii) for ii in fnum_c_4_30
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_30_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Open Loop - 60s
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum_o_60]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_60_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 3 WFS Small - 60s
img_files = [
data_dir + 'obj{0:04d}_threewfs_small_c.fits'.format(ii)
for ii in fnum_c_3S_60
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_60_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 3 WFS Large - 60s
img_files = [
data_dir + 'obj{0:04d}_threeWFS_big_c.fits'.format(ii)
for ii in fnum_c_3L_60
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_60_sky.fits',
flat_frame=flat_dir + "flat.fits")
# Closed - 4 WFS - 60s
img_files = [
data_dir + 'obj{0:04d}_fourWFS_c.fits'.format(ii) for ii in fnum_c_4_60
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'FLD2_60_sky.fits',
flat_frame=flat_dir + "flat.fits")
return
def reduce_pleiades():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
util.mkdir(out_dir)
##########
# R-band
##########
# Open Loop 15 s
fnum = [31, 38, 45, 52, 59, 66]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - a, 15s
fnum = [30, 37, 44, 51, 58, 65, 72]
img_files = [data_dir + 'obj{0:04d}_ca.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - b, 15s
fnum = [28, 35, 42, 49, 56, 63, 70]
img_files = [data_dir + "obj{0:04d}_cb.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - c, 15s
fnum = [27, 34, 41, 48, 55, 62, 69]
img_files = [data_dir + "obj{0:04d}_cc.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - d, 15s
fnum = [26, 33, 40, 47, 54, 61, 68]
img_files = [data_dir + "obj{0:04d}_cd.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - e, 15s
fnum = [25, 32, 39, 46, 53, 60, 67]
img_files = [data_dir + "obj{0:04d}_ce.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# TT, 15s
fnum = [29, 36, 43, 50, 57, 64, 71]
img_files = [data_dir + "obj{0:04d}_tt.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
##########
# I-band
##########
# Open Loop 15 s
fnum = [94, 101, 108, 115, 129]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - a, 15s
fnum = [93, 100, 107, 114, 128]
img_files = [data_dir + 'obj{0:04d}_ca.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - b, 15s
fnum = [91, 98, 105, 112, 119]
img_files = [data_dir + "obj{0:04d}_cb.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - c, 15s
fnum = [90, 97, 104, 111, 118, 125]
img_files = [data_dir + "obj{0:04d}_cc.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - d, 15s
fnum = [89, 96, 103, 110, 117, 124]
img_files = [data_dir + "obj{0:04d}_cd.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - e, 15s
fnum = [88, 95, 102, 109, 116, 123]
img_files = [data_dir + "obj{0:04d}_ce.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
# TT, 15s
fnum = [92, 99, 106, 113, 120, 127]
img_files = [data_dir + "obj{0:04d}_tt.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_r_15s.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_FLD2():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'FLD2_2/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/FLD2_2/'
util.mkdir(out_dir)
# Open Loop
fnum = [124, 129, 135, 141, 147, 153, 159, 165, 171, 177, 183, 189, 195]
fnum += [201, 207, 215, 221, 227, 233, 239, 245, 251, 257, 263, 269]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame=flat_dir + 'flat.fits')
# Closed
fnum = [126, 128, 134, 140, 146, 152, 158, 164, 170, 176, 182, 188, 194]
fnum += [200, 206, 212, 214, 220, 226, 232, 238, 244, 250, 256, 262, 268]
img_files = [data_dir + 'obj{0:04d}_c.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame = flat_dir + 'flat.fits')
# Closed A
fnum = [125, 127, 130, 136, 142, 148, 154, 160, 166, 172, 178, 184, 190]
fnum += [196, 202, 208, 216, 222, 228, 234, 240, 246, 252, 258, 264, 270]
img_files = [data_dir + 'obj{0:04d}_cA.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame = flat_dir + 'flat.fits')
# Closed B
fnum = [131, 137, 143, 149, 155, 161, 167, 173, 179, 185, 191, 197, 203]
fnum += [209, 217, 223, 229, 235, 241, 247, 253, 259, 265, 271]
img_files = [data_dir + 'obj{0:04d}_cB.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame = flat_dir + 'flat.fits')
# Closed C
fnum = [132, 138, 144, 150, 156, 162, 168, 174, 180, 186, 192, 198, 204]
fnum += [210, 218, 224, 230, 236, 242, 248, 254, 260, 266, 272]
img_files = [data_dir + 'obj{0:04d}_cC.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame = flat_dir + 'flat.fits')
# Closed D
fnum = [133, 139, 145, 151, 157, 163, 169, 175, 181, 187, 193, 199, 205]
fnum += [211, 219, 225, 231, 237, 243, 249, 255, 261, 267, 273]
img_files = [data_dir + 'obj{0:04d}_cD.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files, out_dir, rebin=1, sky_frame=sky_dir + 'FLD2_2_sky.fits', flat_frame = flat_dir + 'flat.fits')
return
def reduce_pleiades():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
#util.mkdir(out_dir)
# Open Loop
fnum = [22, 28, 34, 40, 46]
fnum += [52, 58, 65, 72, 79, 86]
fnum += [98, 105, 112, 119, 126, 133, 140, 147]
fnum += [159, 166, 173, 180, 187, 195, 202, 209]
fnum += [216, 223, 230, 237]
img_files = [data_dir + 'obj{0:04d}_o.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - a
fnum = [27, 33, 39, 45, 51]
fnum += [57, 64, 71, 78, 85]
fnum += [104, 111, 118, 125, 132, 139, 146]
fnum += [153, 165, 172, 179, 194, 201, 208]
fnum += [215, 222, 229, 236]
img_files = [data_dir + 'obj{0:04d}_ca.fits'.format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - b
fnum = [25, 31, 37, 43, 49]
fnum += [55, 62, 69, 76, 83, 90]
fnum += [102, 109, 116, 123, 130, 137, 144, 151]
fnum += [163, 170, 177, 184, 191, 199, 206, 213]
fnum += [220, 227, 234, 241]
img_files = [data_dir + "obj{0:04d}_cb.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - c
fnum = [24, 30, 36, 42, 48]
fnum += [54, 61, 68, 75, 82, 89]
fnum += [101, 108, 115, 122, 129, 136, 143, 150]
fnum += [162, 169, 176, 183, 190, 198, 205, 212]
fnum += [219, 226, 233, 240]
img_files = [data_dir + "obj{0:04d}_cc.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - d
fnum = [23, 29, 35, 41, 47]
fnum += [53, 60, 67, 74, 81, 88]
fnum += [100, 107, 114, 121, 128, 135, 142, 149]
fnum += [161, 168, 175, 182, 189, 197, 204, 211]
fnum += [218, 225, 232, 239]
img_files = [data_dir + "obj{0:04d}_cd.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
# Closed Loop - f
fnum = [26, 32, 38, 44, 50]
fnum += [56, 63, 70, 77, 84, 91, 92]
fnum += [103, 110, 117, 124, 131, 138, 145, 152]
fnum += [164, 171, 178, 185, 192, 200, 207, 214]
fnum += [221, 228, 235]
img_files = [data_dir + "obj{0:04d}_cf.fits".format(ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky.fits',
flat_frame=flat_dir + 'flat.fits')
return
def reduce_pleiades():
sky_dir = root_dir + 'reduce/sky/'
data_dir = root_dir + 'Pleiades/'
flat_dir = root_dir + 'reduce/calib/'
out_dir = root_dir + 'reduce/pleiades/'
util.mkdir(out_dir)
#####
# Open Loop First Half
#####
###What to do with obj_o594-597??? look at images to see orientation
###obj_o640 has a 0.5s integration time...whut
###obj_o613.fits in log but not in data.... also 624
fnum1 = [
34, 35, 36, 43, 44, 45, 56, 57, 58, 59, 68, 69, 70, 71, 83, 84, 85, 86
]
fnum2 = [95, 96, 97, 98, 480, 582, 583, 584, 585]
fnum = fnum1 + fnum2
img_files = ['{0:s}/obj_o{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_46.fits',
flat_frame=flat_dir + 'flat.fits')
#####
# Open Loop Second Half
#####
fnum1 = [609, 610, 611, 612, 622, 623, 624, 636, 637, 638, 639, 1010, 1011]
fnum2 = [
1012, 1013, 1022, 1023, 1024, 1025, 1034, 1035, 1036, 1037, 1049, 1052,
1061
]
fnum3 = [
1062, 1063, 1064, 1073, 1074, 1075, 1076, 1088, 1089, 1090, 1091, 1104,
1105
]
fnum4 = [1106, 1107]
fnum = fnum1 + fnum2 + fnum3 + fnum4
img_files = ['{0:s}/obj_o{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_598.fits',
flat_frame=flat_dir + 'flat.fits')
#####
# TTF Closed Loop First Half
#####
###obj_ttf541.fits in log but not data
fnum1 = [
37, 38, 39, 49, 50, 51, 60, 61, 62, 63, 75, 76, 77, 78, 87, 88, 89, 90
]
fnum2 = [586, 587, 588, 589]
fnum = fnum1 + fnum2
img_files = [
'{0:s}/obj_ttf{1:03d}.fits'.format(data_dir, ii) for ii in fnum
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_46.fits',
flat_frame=flat_dir + 'flat.fits')
#####
# TTF Closed Loop Second Half
#####
###obj_ttf617 not in file but in log
fnum1 = [
601, 602, 603, 604, 614, 615, 616, 629, 630, 631, 1014, 1015, 1016,
1017
]
fnum2 = [
1026, 1027, 1028, 1029, 1041, 1042, 1043, 1044, 1053, 1054, 1055, 1056,
1065
]
fnum3 = [1066, 1067, 1068, 1080, 1081, 1082, 1083, 1092, 1093, 1094, 1095]
fnum = fnum1 + fnum2 + fnum3
img_files = [
'{0:s}/obj_ttf{1:03d}.fits'.format(data_dir, ii) for ii in fnum
]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_598.fits',
flat_frame=flat_dir + 'flat.fits')
#####
# Closed Loop First Half
#####
#obj_c102-103, 407-416, 573-574 have 1s or 0.1s integration times... not reduced
fnum1 = [
31, 32, 33, 40, 41, 42, 52, 53, 54, 55, 64, 65, 66, 67, 79, 80, 81, 82,
91, 92
]
fnum2 = [
93, 94, 417, 448, 479, 577, 578, 579, 580, 581, 590, 591, 592, 593
]
fnum = fnum1 + fnum2
img_files = ['{0:s}/obj_c{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_46.fits',
flat_frame=flat_dir + 'flat.fits')
#####
# Closed Loop Second Half
#####
###obj_c641 has 0.5s integration time
###obj_c621.fits doesn't exist in file, only in log
fnum1 = [
605, 606, 607, 608, 618, 619, 620, 632, 633, 634, 635, 1006, 1007, 1008
]
fnum2 = [
1009, 1018, 1019, 1020, 1021, 1030, 1031, 1032, 1033, 1045, 1046, 1047,
1048
]
fnum3 = [
1057, 1058, 1059, 1060, 1069, 1070, 1071, 1072, 1084, 1085, 1086, 1087,
1096
]
fnum4 = [1097, 1098, 1099]
fnum = fnum1 + fnum2 + fnum3 + fnum4
img_files = ['{0:s}/obj_c{1:03d}.fits'.format(data_dir, ii) for ii in fnum]
reduce_fli.clean_images(img_files,
out_dir,
rebin=1,
sky_frame=sky_dir + 'pleiades_sky_598.fits',
flat_frame=flat_dir + 'flat.fits')
return