Пример #1
0
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"""
Пример #2
0
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"""
Пример #3
0
#--------------------------------------

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 = []
Пример #4
0
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)