def load_echelogram(ref_date, band):
from echellogram import Echellogram
echel_name = get_master_calib_abspath("fitted_echellogram_sky_%s_%s.json" % (band, ref_date))
echel = Echellogram.from_json_fitted_echellogram_sky(echel_name)
return echel
r = pickle.load(open("flat_info_%s_%s.pickle" % (igrins_log.date, band)))
r2 = pickle.load(open("thar_%s_%s.pickle" % (igrins_log.date, band)))
bpix_mask = r["flat_bpix_mask"]
trace_sol = r["bottomup_solutions"]
from apertures import Apertures
apertures = Apertures(igrins_orders[band],
trace_sol)
# load echellogram data
from echellogram import Echellogram
echel_name = "fitted_echellogram_sky_%s_%s.json" % (band, thar_init["ref_date"])
echel = Echellogram.from_json_fitted_echellogram_sky(echel_name)
th = np.genfromtxt("ThArlines.dat")
wvl_thar = th[:,0]/1.e4
s_thar = np.clip(th[:,1], a_min=20, a_max=np.inf)
# line_list : dict of (order, (pixel coord list, wavelengths))
wvl_list = {}
pixel_list = {}
for o, s in zip(igrins_orders[band], thar_init["match_list"]):
lineid_list = s[0] # [s1[0] for s1 in s]
wvl = wvl_thar[lineid_list]
wvl_list[o] = wvl
x = [s1[0] for s1 in s[1]]
pixel_list[o] = x
igrins_orders["K"] = range(72, 94)
r = pickle.load(open("flat_info_%s_%s.pickle" % (igrins_log.date, band)))
r2 = pickle.load(open("thar_%s_%s.pickle" % (igrins_log.date, band)))
bpix_mask = r["flat_bpix_mask"]
trace_sol = r["bottomup_solutions"]
from apertures import Apertures
apertures = Apertures(igrins_orders[band], trace_sol)
# load echellogram data
from echellogram import Echellogram
echel_name = "fitted_echellogram_sky_%s_%s.json" % (band,
thar_init["ref_date"])
echel = Echellogram.from_json_fitted_echellogram_sky(echel_name)
th = np.genfromtxt("ThArlines.dat")
wvl_thar = th[:, 0] / 1.e4
s_thar = np.clip(th[:, 1], a_min=20, a_max=np.inf)
# line_list : dict of (order, (pixel coord list, wavelengths))
wvl_list = {}
pixel_list = {}
for o, s in zip(igrins_orders[band], thar_init["match_list"]):
lineid_list = s[0] # [s1[0] for s1 in s]
wvl = wvl_thar[lineid_list]
wvl_list[o] = wvl
x = [s1[0] for s1 in s[1]]
pixel_list[o] = x
