"""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"""
fig,(ax1,ax2) = subplots(1,2,figsize=(18,8))
fig.subplots_adjust(bottom=0.15)
for k in list_filters:
ax1.plot(time_series[:],radiance_series[:,k],"o-",color=colors[k])
ax2.plot(time_series[:],bt_series[:,k]-273.15,"o-",color=colors[k])
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"""
fig,(ax1,ax2) = subplots(1,2,figsize=(18,8))
fig.subplots_adjust(bottom=0.15)
for k in list_filters:
ax1.plot(time_series[:],radiance_series[:,k],"o-",color=colors[k])
ax2.plot(time_series[:],bt_series[:,k]-273.15,"o-",color=colors[k])
