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 = []
time = []
for seq in data.sequences[40:]:
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)
#nf = 30
print nf
counts_nadir = ma.zeros([11,nf])
