def stack():
## EDITED FOR 4F DATA
util.mkdir(stacks_dir)
## Loop through all the different data sets
#for key in ['set_name']: ## 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)
img_files = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean.fits'.format(img=ii, suf=suf)
for ii in img
]
starlists = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean_stars.txt'.format(img=ii, suf=suf)
for ii in img
]
output_root = stacks_dir + 'fld2_stack_' + suf + '_' + filt ## EDITED LINE
redu.shift_and_add(img_files, starlists, output_root, method='mean')
return
def stack_beehive_I():
# I Band - open
images = [
out_dir + 'obj{0:03d}_o_scan_clean.fits'.format(ii) for ii in fnum_o_I
]
starlists = [
out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii)
for ii in fnum_o_I
]
output_root = stacks_dir + 'beehive_stack_open_I'
reduce_fli.shift_and_add(images, starlists, output_root, method='mean')
# I Band - closed
images = [
out_dir + 'obj{0:03d}LS4WFS_c_scan_clean.fits'.format(ii)
for ii in fnum_c_I
]
starlists = [
out_dir + 'obj{0:03d}LS4WFS_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_c_I
]
output_root = stacks_dir + 'beehive_stack_closed_I'
reduce_fli.shift_and_add(images, starlists, output_root, method='mean')
return
def stack():
util.mkdir(stacks_dir)
## Loop through all datasets
for key in dict_suffix.keys():
img = dict_images[key]
suf = dict_suffix[key]
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
img_files = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean.fits'.format(img=ii, suf=suf)
for ii in img
]
starlists = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean_stars.txt'.format(img=ii, suf=suf)
for ii in img
]
output_root = stacks_dir + 'fld2_stack_' + suf
redu.shift_and_add(img_files, starlists, output_root, method='mean')
return
def stack_FLD2():
util.mkdir(stacks_dir)
# Open Loop
open_images = [
out_dir + 'obj{0:03d}_o_scan_clean.fits'.format(ii) for ii in fnum_o
]
open_starlists = [
out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii)
for ii in fnum_o
]
open_output_root = stacks_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop - threeWFS_LS
closed_images = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean.fits'.format(ii)
for ii in fnum_c
]
closed_starlists = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_c
]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFS_LS_c'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_beehive():
util.mkdir(stacks_dir)
# Loop through all the different data sets and reduce them.
#for key in dict_suffix.keys():
for key in ['open_2', 'LS_3wfs_s_2', 'LS_3wfs_w_2', 'LS_5wfs_2']:
img = dict_images[key]
suf = dict_suffix[key]
print('Working on: {1:s} {0:s}'.format(key, suf))
print(' Images: ', img)
img_files = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean.fits'.format(img=ii, suf=suf)
for ii in img
]
starlists = [
out_dir +
'sta{img:03d}{suf:s}_scan_clean_stars.txt'.format(img=ii, suf=suf)
for ii in img
]
output_root = stacks_dir + 'beehive_stack_' + suf
redu.shift_and_add(img_files, starlists, output_root, method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2_2/'
stats_dir = root_dir + 'reduce/stats/'
stacks_dir = root_dir + 'reduce/stacks/'
util.mkdir(stacks_dir)
# Open Loop
open_img_nums = [
14, 17, 29, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53, 56, 59
]
open_img_nums += [62, 65, 68, 71, 74, 77, 80, 83, 86, 89, 92]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in open_img_nums
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums
]
open_output_root = stacks_dir + 'FLD2_2_stack_open'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed A
closed_img_nums = [
12, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55
]
closed_img_nums += [58, 61, 64, 67, 70, 73, 76, 79, 82, 85, 88, 91]
closed_images = [
data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedA'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed D
closed_img_nums = [
18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63
]
closed_img_nums += [66, 69, 72, 75, 78, 81, 84, 87, 90, 93]
closed_images = [
data_dir + 'obj{0:04d}_cD_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cD_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedD'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_FLD2():
util.mkdir(stacks_dir)
# Open Loop
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum_o
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum_o
]
open_output_root = stacks_dir + 'Beehive-W_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop
closed_images = [
data_dir + 'obj{0:04d}_c_clean.fits'.format(ii) for ii in fnum_c
]
closed_starlists = [
data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii) for ii in fnum_c
]
closed_output_root = stacks_dir + 'Beehive-w_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
#def analyze_stacks():
# img_files = [data_dir + 'FLD2_stack_open.fits',
# data_dir + 'FLD2_stack_closed.fits',
# data_dir + 'FLD2_stack_closedA.fits',
# data_dir + 'FLD2_stack_closedB.fits',
# data_dir + 'FLD2_stack_closedC.fits',
# data_dir + 'FLD2_stack_closedD.fits']
#Find stars in image
# reduce_fli.find_stars(img_files, fwhm=5, threshold=10, N_passes=2, plot_psf_compare=False)
# Calc stats on all the stacked images
# reduce_fli.calc_star_stats(img_files, output_stats= stats_dir + 'stats_stacks.fits')
return
def stack_orion():
util.mkdir(stacks_dir)
# Open Loop
open_images = [
out_dir + 'obj{0:03d}_o_scan_clean.fits'.format(ii) for ii in fnum_o
]
open_starlists = [
out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii)
for ii in fnum_o
]
open_output_root = stacks_dir + 'orion_stack_open'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - LS
closed_images = [
out_dir + 'obj{0:03d}LS_c_scan_clean.fits'.format(ii)
for ii in fnum_c_LS
]
closed_starlists = [
out_dir + 'obj{0:03d}LS_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_c_LS
]
closed_output_root = stacks_dir + 'orion_stack_closed_LS'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed Loop - B2
closed_images = [
out_dir + 'obj{0:03d}LS_Bin2_c_scan_clean.fits'.format(ii)
for ii in fnum_c_B2
]
closed_starlists = [
out_dir + 'obj{0:03d}LS_Bin2_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_c_B2
]
closed_output_root = stacks_dir + 'orion_stack_closed_B2'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_orion():
util.mkdir(stacks_dir)
# Open Loop
open_images = [out_dir + 'obj{0:03d}_o_scan_clean.fits'.format(ii) for ii in fnum_o]
open_starlists = [out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii) for ii in fnum_o]
open_output_root = stacks_dir + 'orion_stack_open'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - 4W
closed_images = [out_dir + 'obj{0:03d}LS4WFS_c_scan_clean.fits'.format(ii) for ii in fnum_c_4W]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_4W]
closed_output_root = stacks_dir + 'orion_stack_closed_4W'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - B2
closed_images = [out_dir + 'obj{0:03d}LS4WFS_B2_c_scan_clean.fits'.format(ii) for ii in fnum_c_B2]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_B2_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_B2]
closed_output_root = stacks_dir + 'orion_stack_closed_B2'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
#Closed Loop - z1
closed_images = [out_dir + 'obj{0:03d}LS4WFS_zc11_c_scan_clean.fits'.format(ii) for ii in fnum_c_z1]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_zc11_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_z1]
closed_output_root = stacks_dir + 'orion_stack_closed_z1'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
#Closed Loop - z2
closed_images = [out_dir + 'obj{0:03d}LS4WFS_zc21_c_scan_clean.fits'.format(ii) for ii in fnum_c_z2]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_zc21_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_z2]
closed_output_root = stacks_dir + 'orion_stack_closed_z2'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
#tt
closed_images = [out_dir + 'obj{0:03d}tip_tilt_scan_clean.fits'.format(ii) for ii in fnum_tt]
closed_starlists = [out_dir + 'obj{0:03d}tip_tilt_scan_clean_stars.txt'.format(ii) for ii in fnum_tt]
closed_output_root = stacks_dir + 'orion_stack_tiptilt'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
return
def stack_pleiades_band():
"""
Stack all the images for the Open loop, TTF closed loop, and GLAO closed loop (a)
for both R and I band data.
"""
data_dir = root_dir + 'reduce/pleiades/'
stats_dir = root_dir + 'reduce/stats/'
# Open loop R-band
open_img_nums = [31, 38, 45, 52, 59, 66]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in open_img_nums
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums
]
open_output_root = data_dir + 'pleiades_stack_open_r'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop - a
closed_img_nums = [30, 37, 44, 51, 58, 65, 72]
closed_images = [
data_dir + 'obj{0:04d}_ca_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_ca_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = data_dir + 'pleiades_stack_closed_r'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# TT
tt_img_nums = [29, 36, 43, 50, 57, 64, 71]
tt_images = [
data_dir + "obj{0:04d}_tt_clean.fits".format(ii) for ii in tt_img_nums
]
tt_starlists = [
data_dir + "obj{0:04d}_tt_clean_stars.txt".format(ii)
for ii in tt_img_nums
]
tt_output_root = data_dir + 'pleiades_stack_tt_r'
reduce_fli.shift_and_add(tt_images,
tt_starlists,
tt_output_root,
method='mean')
##########
# I-band
##########
# Open Loop
open_img_nums_i = [94, 101, 108, 115, 129]
open_images_i = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii)
for ii in open_img_nums_i
]
open_starlists_i = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums_i
]
open_output_root_i = data_dir + 'pleiades_stack_open_i'
reduce_fli.shift_and_add(open_images_i,
open_starlists_i,
open_output_root_i,
method='mean')
# Closed Loop - a
closed_img_nums_i = [93, 100, 107, 114, 128]
closed_images_i = [
data_dir + 'obj{0:04d}_ca_clean.fits'.format(ii)
for ii in closed_img_nums_i
]
closed_starlists_i = [
data_dir + 'obj{0:04d}_ca_clean_stars.txt'.format(ii)
for ii in closed_img_nums_i
]
closed_output_root_i = data_dir + 'pleiades_stack_closed_i'
reduce_fli.shift_and_add(closed_images_i,
closed_starlists_i,
closed_output_root_i,
method='mean')
# TT
tt_img_nums_i = [92, 99, 106, 113, 120, 127]
tt_images_i = [
data_dir + "obj{0:04d}_tt_clean.fits".format(ii)
for ii in tt_img_nums_i
]
tt_starlists_i = [
data_dir + "obj{0:04d}_tt_clean_stars.txt".format(ii)
for ii in tt_img_nums_i
]
tt_output_root_i = data_dir + 'pleiades_stack_tt_i'
reduce_fli.shift_and_add(tt_images_i,
tt_starlists_i,
tt_output_root_i,
method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2_2/'
stats_dir = root_dir + 'reduce/stats/'
stacks_dir = root_dir + 'reduce/stacks/'
util.mkdir(stacks_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]
open_images = [data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum]
open_starlists = [data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum]
open_output_root = stacks_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# 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]
closed_images = [data_dir + 'obj{0:04d}_c_clean.fits'.format(ii) for ii in fnum]
closed_starlists = [data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii) for ii in fnum]
closed_output_root = stacks_dir + 'FLD2_stack_closed'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# 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]
closed_images = [data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii) for ii in fnum]
closed_starlists = [data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii) for ii in fnum]
closed_output_root = stacks_dir + 'FLD2_stack_closedA'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# 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]
closed_images = [data_dir + 'obj{0:04d}_cB_clean.fits'.format(ii) for ii in fnum]
closed_starlists = [data_dir + 'obj{0:04d}_cB_clean_stars.txt'.format(ii) for ii in fnum]
closed_output_root = stacks_dir + 'FLD2_stack_closedB'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# 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]
closed_images = [data_dir + 'obj{0:04d}_cC_clean.fits'.format(ii) for ii in fnum]
closed_starlists = [data_dir + 'obj{0:04d}_cC_clean_stars.txt'.format(ii) for ii in fnum]
closed_output_root = stacks_dir + 'FLD2_stack_closedC'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# 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]
closed_images = [data_dir + 'obj{0:04d}_cD_clean.fits'.format(ii) for ii in fnum]
closed_starlists = [data_dir + 'obj{0:04d}_cD_clean_stars.txt'.format(ii) for ii in fnum]
closed_output_root = stacks_dir + 'FLD2_stack_closedD'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2/'
stats_dir = root_dir + 'reduce/stats/'
# Open Loop
open_img_nums = [
63, 67, 71, 75, 80, 84, 88, 92, 96, 100, 104, 108, 112, 124
]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in open_img_nums
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums
]
open_output_root = data_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop
closed_img_nums = [
64, 68, 72, 76, 77, 81, 85, 89, 93, 97, 101, 105, 109, 113
]
closed_images = [
data_dir + 'obj{0:04d}_c_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = data_dir + 'FLD2_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed A
closed_img_nums = [65, 69, 73, 78, 82, 86, 90, 94, 98, 102, 106, 110, 114]
closed_images = [
data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = data_dir + 'FLD2_stack_closedA'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed B
closed_img_nums = [66, 70, 74, 79, 83, 87, 91, 95, 99, 103, 107, 111, 115]
closed_images = [
data_dir + 'obj{0:04d}_cB_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cB_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = data_dir + 'FLD2_stack_closedB'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_FLD2():
util.mkdir(stacks_dir)
# Open Loop
open_images = [out_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum_o]
open_starlists = [out_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum_o]
open_output_root = stacks_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - threeWFS_LS
closed_images = [out_dir + 'obj{0:04d}_threeWFS_LS_clean.fits'.format(ii) for ii in fnum_threeWFS_LS]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFS_LS_clean_stars.txt'.format(ii) for ii in fnum_threeWFS_LS]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFS_LS'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - threeWFSLS_B2_c
closed_images = [out_dir + 'obj{0:04d}_threeWFSLS_B2_c_clean.fits'.format(ii) for ii in fnum_threeWFSLS_B2_c]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFSLS_B2_c_clean_stars.txt'.format(ii) for ii in fnum_threeWFSLS_B2_c]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFSLS_B2_c'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - threeWFSMean_B2_c
closed_images = [out_dir + 'obj{0:04d}_threeWFSMean_B2_c_clean.fits'.format(ii) for ii in fnum_threeWFSMean_B2_c]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFSMean_B2_c_clean_stars.txt'.format(ii) for ii in fnum_threeWFSMean_B2_c]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFSMean_B2_c'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Open Loop - 2 filter
open_images = [out_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum_o_2filt]
open_starlists = [out_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum_o_2filt]
open_output_root = stacks_dir + 'FLD2_stack_open_2filt'
reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - threeWFS_LS - 2 filter
closed_images = [out_dir + 'obj{0:04d}_threeWFS_LS_clean.fits'.format(ii) for ii in fnum_threeWFS_LS_2filt]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFS_LS_clean_stars.txt'.format(ii) for ii in fnum_threeWFS_LS_2filt]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFS_LS_2filt'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - threeWFSLS_B2_c - 2 filter
closed_images = [out_dir + 'obj{0:04d}_threeWFSLS_B2_c_clean.fits'.format(ii) for ii in fnum_threeWFSLS_B2_c_2filt]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFSLS_B2_c_clean_stars.txt'.format(ii) for ii in fnum_threeWFSLS_B2_c_2filt]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFSLS_B2_c_2filt'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - threeWFSMean_B2_c - 2 filter
closed_images = [out_dir + 'obj{0:04d}_threeWFSMean_B2_c_clean.fits'.format(ii) for ii in fnum_threeWFSMean_B2_c_2filt]
closed_starlists = [out_dir + 'obj{0:04d}_threeWFSMean_B2_c_clean_stars.txt'.format(ii) for ii in fnum_threeWFSMean_B2_c_2filt]
closed_output_root = stacks_dir + 'FLD2_stack_threeWFSMean_B2_c_2filt'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2_2/'
stats_dir = root_dir + 'reduce/stats/'
stacks_dir = root_dir + 'reduce/stacks/'
util.mkdir(stacks_dir)
# Open Loop
open_img_nums = [20, 26, 32, 34, 38, 39, 40, 42, 45, 48, 51, 54, 57]
open_img_nums += [60, 63, 66, 69, 72, 75, 78, 81, 84, 87, 90, 93, 95]
open_img_nums += [99, 102, 105, 108, 111, 113, 116, 119, 122, 125, 128]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in open_img_nums
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums
]
open_output_root = stacks_dir + 'FLD2_2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop
closed_img_nums = [19, 21, 22, 23, 24, 25, 27, 28, 29, 30, 31]
closed_images = [
data_dir + 'obj{0:04d}_c_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed A
closed_img_nums = [33, 41, 44, 47, 50, 53, 56, 59, 62, 65, 68]
closed_img_nums += [71, 74, 77, 80, 83, 86, 89, 92, 94]
closed_img_nums += [98, 101, 104, 107, 110, 112, 115, 118, 121, 124, 127]
closed_images = [
data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedA'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed B
closed_img_nums = [35, 36, 37, 43, 46, 49, 52, 55, 58, 61, 64]
closed_img_nums += [67, 70, 73, 76, 79, 82, 85, 88, 91, 96]
fnumclosed_img_nums += [100, 103, 106, 109, 114, 117, 120, 123, 126, 129]
closed_images = [
data_dir + 'obj{0:04d}_cB_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cB_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedB'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_FLD2():
util.mkdir(stacks_dir)
# Open Loop
open_images = [
out_dir + 'obj{0:03d}_o_scan_clean.fits'.format(ii) for ii in fnum_o
]
open_starlists = [
out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii)
for ii in fnum_o
]
open_output_root = stacks_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# EXPERIMENT 1
# Closed - LS_c
img_files = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean.fits'.format(ii)
for ii in fnum_LS_c
]
starlists = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_LS_c
]
output_root = stacks_dir + 'FLD2_stack_closed_LS'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
# Closed - LS_B2_c
img_files = [
out_dir + 'obj{0:03d}threeWFSLS_B2_c_scan_clean.fits'.format(ii)
for ii in fnum_LS_B2_c
]
starlists = [
out_dir + 'obj{0:03d}threeWFSLS_B2_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_LS_B2_c
]
output_root = stacks_dir + 'FLD2_stack_closed_LS_B2'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
# Closed - Mean_B2_c
img_files = [
out_dir + 'obj{0:03d}threeWFSMean_B2_c_scan_clean.fits'.format(ii)
for ii in fnum_Mean_B2_c
]
starlists = [
out_dir + 'obj{0:03d}threeWFSMean_B2_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_Mean_B2_c
]
output_root = stacks_dir + 'FLD2_stack_closed_Mean_B2'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
# EXPERIMENT 2
# Closed - nrej1
img_files = [
out_dir + 'obj{0:03d}threeWFSLS_nrej1_c_scan_clean.fits'.format(ii)
for ii in fnum_nrej1
]
starlists = [
out_dir +
'obj{0:03d}threeWFSLS_nrej1_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_nrej1
]
output_root = stacks_dir + 'FLD2_stack_closed_nrej1'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
# Closed - nrej4
img_files = [
out_dir + 'obj{0:03d}threeWFSLS_nrej4_c_scan_clean.fits'.format(ii)
for ii in fnum_nrej4
]
starlists = [
out_dir +
'obj{0:03d}threeWFSLS_nrej4_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_nrej4
]
output_root = stacks_dir + 'FLD2_stack_closed_nrej4'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
# Closed - nrej7
img_files = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean.fits'.format(ii)
for ii in fnum_nrej7
]
starlists = [
out_dir + 'obj{0:03d}threeWFS_LS_c_scan_clean_stars.txt'.format(ii)
for ii in fnum_nrej7
]
output_root = stacks_dir + 'FLD2_stack_closed_nrej7'
reduce_fli.shift_and_add(img_files, starlists, output_root, method='mean')
open_starlists = [out_dir + 'obj{0:03d}_o_scan_clean_stars.txt'.format(ii) for ii in fnum_o]
open_output_root = stacks_dir + 'orion_stack_open'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - 4W
closed_images = [out_dir + 'obj{0:03d}LS4WFS_c_scan_clean.fits'.format(ii) for ii in fnum_c_4W]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_4W]
closed_output_root = stacks_dir + 'orion_stack_closed_4W'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - B2
closed_images = [out_dir + 'obj{0:03d}LS4WFS_B2_c_scan_clean.fits'.format(ii) for ii in fnum_c_B2]
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_B2_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_B2]
closed_output_root = stacks_dir + 'orion_stack_closed_B2'
reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
<<<<<<< HEAD
# Closed Loop - zc
closed_images = [out_dir + 'obj{0:03dLS4WFS_zc21_c_scan_clean.fits'.format(ii) for ii in fnum_c_zc]
=======
closed_images = [out_dir + 'obj{0:03d}LS4WFS_zc21_c_scan_clean.fits'.format(ii) for ii in fnum_c_zc]
>>>>>>> 9d0fe2dfc407d195227cb24b248c5a5cff886914
closed_starlists = [out_dir + 'obj{0:03d}LS4WFS_zc21_c_scan_clean_stars.txt'.format(ii) for ii in fnum_c_zc]
closed_output_root = stacks_dir + 'orion_stack_closed_zc'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
return
def analyze_stacks():
def stack_FLD2():
util.mkdir(stacks_dir)
# Open Loop - 30 s
open_images = [
out_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum_o_30
]
open_starlists = [
out_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum_o_30
]
open_output_root = stacks_dir + 'FLD2_stack_open_30'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - 3S WFS - 30 S
closed_images = [
out_dir + 'obj{0:04d}_threewfs_small_c_clean.fits'.format(ii)
for ii in fnum_c_3S_30
]
closed_starlists = [
out_dir + 'obj{0:04d}_threewfs_small_c_clean_stars.txt'.format(ii)
for ii in fnum_c_3S_30
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_3S_30'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - 3L WFS - 30 S
closed_images = [
out_dir + 'obj{0:04d}_threeWFS_big_c_clean.fits'.format(ii)
for ii in fnum_c_3L_30
]
closed_starlists = [
out_dir + 'obj{0:04d}_threeWFS_big_c_clean_stars.txt'.format(ii)
for ii in fnum_c_3L_30
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_3L_30'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Closed Loop - 3L WFS - 30 S
closed_images = [
out_dir + 'obj{0:04d}_fourWFS_c_clean.fits'.format(ii)
for ii in fnum_c_4_30
]
closed_starlists = [
out_dir + 'obj{0:04d}_fourWFS_c_clean_stars.txt'.format(ii)
for ii in fnum_c_4_30
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_4_30'
#reduce_fli.shift_and_add(closed_images, closed_starlists, closed_output_root, method='mean')
# Open Loop - 60 s
open_images = [
out_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum_o_60
]
open_starlists = [
out_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum_o_60
]
open_output_root = stacks_dir + 'FLD2_stack_open_60'
#reduce_fli.shift_and_add(open_images, open_starlists, open_output_root, method='mean')
# Closed Loop - 3S WFS - 60 S
closed_images = [
out_dir + 'obj{0:04d}_threewfs_small_c_clean.fits'.format(ii)
for ii in fnum_c_3S_60
]
closed_starlists = [
out_dir + 'obj{0:04d}_threewfs_small_c_clean_stars.txt'.format(ii)
for ii in fnum_c_3S_60
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_3S_60'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed Loop - 3L WFS - 60 S
closed_images = [
out_dir + 'obj{0:04d}_threeWFS_big_c_clean.fits'.format(ii)
for ii in fnum_c_3L_60
]
closed_starlists = [
out_dir + 'obj{0:04d}_threeWFS_big_c_clean_stars.txt'.format(ii)
for ii in fnum_c_3L_60
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_3L_60'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed Loop - 3L WFS - 60 S
closed_images = [
out_dir + 'obj{0:04d}_fourWFS_c_clean.fits'.format(ii)
for ii in fnum_c_4_60
]
closed_starlists = [
out_dir + 'obj{0:04d}_fourWFS_c_clean_stars.txt'.format(ii)
for ii in fnum_c_4_60
]
closed_output_root = stacks_dir + 'FLD2_stack_closed_4_60'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_pleiades_west_ttf():
"""
Stack all the images from the night for the Open loop, TTF closed loop, and GLAO closed loop
images. Calculate statstics on all three.
"""
red_dir = root_dir + 'reduce/pleiades/'
# open second half of night (on west field
open_img_nums = [
609, 610, 611, 612, 622, 623, 624, 636, 637, 638, 639, 1010, 1011
]
open_img_nums += [
1012, 1013, 1022, 1023, 1024, 1025, 1034, 1035, 1036, 1037, 1049, 1052,
1061
]
open_img_nums += [
1062, 1063, 1064, 1073, 1074, 1075, 1076, 1088, 1089, 1090, 1091, 1104,
1105
]
open_img_nums += [1106, 1107]
open_images = [
'{0:s}/obj_o{1:03d}_clean.fits'.format(red_dir, ii)
for ii in open_img_nums
]
open_starlists = [
'{0:s}/obj_o{1:03d}_clean.txt'.format(red_dir, ii)
for ii in open_img_nums
]
open_output_root = 'west_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# TTF closed second half of night (on west field)
ttf_img_nums = [
601, 602, 603, 604, 614, 615, 616, 629, 630, 631, 1014, 1015, 1016,
1017
]
ttf_img_nums += [
1026, 1027, 1028, 1029, 1041, 1042, 1043, 1044, 1053, 1054, 1055, 1056,
1065
]
ttf_img_nums += [
1066, 1067, 1068, 1080, 1081, 1082, 1083, 1092, 1093, 1094, 1095
]
ttf_images = [
'{0:s}/obj_ttf{1:03d}_clean.fits'.format(red_dir, ii)
for ii in ttf_img_nums
]
ttf_starlists = [
'{0:s}/obj_ttf{1:03d}_clean_stars.txt'.format(red_dir, ii)
for ii in ttf_img_nums
]
ttf_output_root = 'west_stack_ttf'
reduce_fli.shift_and_add(ttf_images,
ttf_starlists,
ttf_output_root,
method='mean')
# GLAO closed second half of night (on west field
closed_img_nums = [
605, 606, 607, 608, 618, 619, 620, 632, 633, 634, 635, 1006, 1007, 1008
]
closed_img_nums += [
1009, 1018, 1019, 1020, 1021, 1030, 1032, 1033, 1045, 1046, 1047, 1048
]
closed_img_nums += [
1057, 1058, 1059, 1060, 1069, 1070, 1071, 1072, 1084, 1085, 1086, 1087,
1096
]
closed_img_nums += [1097, 1098, 1099]
closed_images = [
'{0:s}/obj_c{1:03d}_clean.fits'.format(red_dir, ii)
for ii in closed_img_nums
]
closed_starlists = [
'{0:s}/obj_c{1:03d}_clean_stars.txt'.format(red_dir, ii)
for ii in closed_img_nums
]
closed_output_root = 'west_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2_2/'
stats_dir = root_dir + 'reduce/stats/'
stacks_dir = root_dir + 'reduce/stacks/'
util.mkdir(stacks_dir)
# Open Loop
open_img_nums = [4, 5, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45]
open_img_nums += [51, 54, 57, 60, 63, 66, 69, 72, 75, 78, 81, 85, 88, 89]
open_img_nums += [92, 95, 99, 103, 107, 110, 113, 116, 119, 122, 125, 128]
open_img_nums += [
131, 134, 137, 140, 143, 146, 149, 152, 155, 158, 161, 164, 167, 170
]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in open_img_nums
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii)
for ii in open_img_nums
]
open_output_root = stacks_dir + 'FLD2_2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed Loop
closed_img_nums = [96, 100, 104]
closed_images = [
data_dir + 'obj{0:04d}_c_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed A
closed_img_nums = [
6, 8, 11, 14, 17, 20, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53
]
closed_img_nums += [
56, 59, 62, 65, 68, 71, 74, 77, 80, 83, 87, 91, 94, 98, 102, 106
]
closed_img_nums += [
109, 112, 115, 118, 121, 124, 127, 130, 133, 136, 139, 142, 145
]
closed_img_nums += [148, 151, 154, 157, 160, 163, 166, 169]
closed_images = [
data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedA'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed B
closed_img_nums = [
10, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58
]
closed_img_nums += [
61, 64, 67, 70, 73, 76, 79, 82, 86, 90, 93, 97, 101, 105, 108, 111
]
closed_img_nums += [
114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 144, 147, 150, 153
]
closed_img_nums += [156, 159, 162, 165, 168, 171]
closed_images = [
data_dir + 'obj{0:04d}_cB_clean.fits'.format(ii)
for ii in closed_img_nums
]
closed_starlists = [
data_dir + 'obj{0:04d}_cB_clean_stars.txt'.format(ii)
for ii in closed_img_nums
]
closed_output_root = stacks_dir + 'FLD2_2_stack_closedB'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
def stack_FLD2():
data_dir = root_dir + 'reduce/FLD2_2/'
stats_dir = root_dir + 'reduce/stats/'
stacks_dir = root_dir + 'reduce/stacks/'
util.mkdir(stacks_dir)
# Open Loop
fnum = [
51, 53, 56, 59, 62, 65, 68, 71, 74, 77, 80, 83, 86, 89, 92, 96, 99, 102
]
fnum += [
105, 108, 111, 114, 117, 120, 123, 126, 129, 132, 135, 138, 141, 144
]
fnum += [
147, 150, 153, 156, 159, 162, 165, 168, 171, 180, 183, 186, 189, 192
]
fnum += [195, 198, 201, 204, 207, 210, 213]
open_images = [
data_dir + 'obj{0:04d}_o_clean.fits'.format(ii) for ii in fnum
]
open_starlists = [
data_dir + 'obj{0:04d}_o_clean_stars.txt'.format(ii) for ii in fnum
]
open_output_root = stacks_dir + 'FLD2_stack_open'
reduce_fli.shift_and_add(open_images,
open_starlists,
open_output_root,
method='mean')
# Closed
fnum = [
49, 52, 55, 58, 61, 64, 67, 70, 73, 76, 79, 82, 85, 88, 91, 95, 98, 101
]
fnum += [
104, 107, 110, 113, 116, 119, 122, 125, 128, 131, 134, 137, 140, 143,
146
]
fnum += [
149, 152, 155, 158, 161, 164, 167, 170, 173, 176, 179, 182, 185, 188,
191
]
fnum += [194, 197, 200, 203, 206, 209, 212]
closed_images = [
data_dir + 'obj{0:04d}_c_clean.fits'.format(ii) for ii in fnum
]
closed_starlists = [
data_dir + 'obj{0:04d}_c_clean_stars.txt'.format(ii) for ii in fnum
]
closed_output_root = stacks_dir + 'FLD2_stack_closed'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
# Closed A
fnum = [
50, 54, 57, 60, 63, 66, 69, 72, 75, 78, 81, 84, 87, 90, 93, 97, 100,
103, 106
]
fnum += [
109, 112, 115, 118, 121, 124, 127, 130, 133, 136, 139, 142, 145, 148,
151, 154
]
fnum += [
157, 160, 163, 166, 169, 172, 175, 178, 181, 184, 187, 190, 193, 196,
199, 202
]
fnum += [205, 208, 211, 214]
closed_images = [
data_dir + 'obj{0:04d}_cA_clean.fits'.format(ii) for ii in fnum
]
closed_starlists = [
data_dir + 'obj{0:04d}_cA_clean_stars.txt'.format(ii) for ii in fnum
]
closed_output_root = stacks_dir + 'FLD2_stack_closedA'
reduce_fli.shift_and_add(closed_images,
closed_starlists,
closed_output_root,
method='mean')
return
