rc('text', usetex=True) #-------------------------------------- rc('xtick',labelsize=22) rc('ytick',labelsize=22) formatter = DateFormatter('%H:%M') """Analyze a series of FIRR sequences""" filters_positions = ['open','F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014','blank'] filters = ['F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014'] list_filters = [filters_positions.index(f) for f in filters] colors=["gold","red","blue","black","orange","yellow","cyan","green","Chartreuse","LightBlue"] series = "Data/Series_2015-07-02" data = FirrSeries(series) date_series = data.date print "Date of measurements:",date_series npos=10 data.get_radiance(filters,npos,method="next",save=False) time_series = data.time_series bt_series = data.bt_series bt_series = ma.masked_equal(bt_series,0) radiance_series = data.radiance_series radiance_series = ma.masked_equal(radiance_series,0) """Figure"""
rc('ytick',labelsize=22) formatter = DateFormatter('%H:%M') """Analyze a series of FIRR sequences""" filters_positions = ['open','F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014','blank'] filters = ['F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014'] list_filters = [filters_positions.index(f) for f in filters] colors=["grey","DarkOrange","Red","Chartreuse","cyan","Gold","LightBlue","green","black","blue"] series = "/media/quentin/LACIE SHARE/FIRR_measurements/2015-06-UQAM/2015-07-02" series = "/media/quentin/LACIE SHARE/FIRR_measurements/2015-02-INO/2015-02-20/Sky" start_meas=datetime(2015,2,20,19,25) end_meas=datetime(2015,2,20,19,55) data = FirrSeries(series,start_meas,end_meas,detector=4,illuminated_pixels=109) date_series = data.date print "Date of measurements:",date_series npos=11 data.get_radiance(filters,npos,method="next",spav="None",non_ill=1,save=True) time_series = data.time_series bt_series = data.bt_series bt_series = ma.masked_equal(bt_series,0) radiance_series = data.radiance_series radiance_series = ma.masked_equal(radiance_series,0) """Figure"""
#-------------------------------------- rc('xtick',labelsize=22) rc('ytick',labelsize=22) formatter = DateFormatter('%H:%M') """Analyze a series of FIRR sequences""" filters_positions = ['open','F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014','blank'] filters = ['F0008','F0009','F0034','F0011','F0007','F0036','F0035','F0010','F0014'] list_filters = [filters_positions.index(f) for f in filters] colors=["gold","red","blue","black","orange","yellow","cyan","green","Chartreuse","LightBlue"] series = "/media/quentin/LACIE SHARE/FIRR_measurements/2015-02-INO/2015-02-22/ZBB_Retrieval/" series = "/media/quentin/LACIE SHARE/FIRR_measurements/2015-02-INO/2015-02-22/Sky" data = FirrSeries(series) date_series = data.date data.get_temperature() plot(data.temp_time,data.temperature) show() npos = 11 print "Date of measurements:",date_series abb = [] hbb = [] abb_non_ill = [] hbb_non_ill = []
new_select = [] count = 0 for k in range(l): count+=1 if count%120==0: if selection[k] == 1: dates_select+=[dates[k]] rad_nadir+=[Toolbox.radiance(t_nadir[k]+273.15,all_filters[j],emiss_wls,emiss,24+273.15) for j in list_filters] n_select = len(dates_select) rad_nadir=reshape(array(rad_nadir),(n_select,11)) dnum = array([date2num(dates[k]) for k in range(l)]) if meas == "SET1": series = FirrSeries("/media/quentin/LACIE SHARE/FIRR_measurements/LR-TECH/SET1") elif meas == "SET2": series = FirrSeries("/media/quentin/LACIE SHARE/FIRR_measurements/LR-TECH/SET2") sequences = series.sequences series.get_temperature() temp_firr = series.temperature temp_firr_time = series.temp_time firr_hbb = temp_firr[:,0] firr_abb = temp_firr[:,1] firr_pm = temp_firr[:,7] rad = [] tms = [] nf = len(sequences)