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