def conti_test(year, month, day, st_hour, st_minute, ed_hour, ed_minute):
d_h = ed_hour - st_hour
d_m = ed_minute - st_minute
time_range = abs(d_h) * 60 + abs(d_m)
open_folder(month, day, year)
files = select_file(month, day, year, st_hour, st_minute, ed_hour,
ed_minute)
minute = []
for i in range(len(files)):
minute.append(int(files[i][0:2]) * 60 + int(files[i][3:5]))
print('# of files:', len(minute))
print('Time range:', time_range)
if len(minute) != time_range:
print "Files are discontinuous"
return minute
else:
print 'Files are continuous'
return minute
def plot_h5(var, year, month, day, st_hour, st_minute, ed_hour, ed_minute):
open_folder(month, day, year)
files = select_file(st_hour, st_minute, ed_hour, ed_minute)
#m=conti_test(year,month,day,st_hour,st_minute,ed_hour,ed_minute)
#size=0
for fname in files:
with h5.File(fname, 'r') as f:
var1 = f['data']['%s' % var]
#el=f['data']['el']
#az=f['data']['az']
#rev=f['data']['rev']
var2 = f['data']['gpstime']
var1 = np.array(var1)
t = np.array(var2)
var1 = var1[np.logical_not(var1 == 0)]
t = t[np.logical_not(t == 0)]
#size=len(var1)
#t=np.linspace(int(m[i-1]),1+int(m[i-1]),size)
#el,=plt.plot(t,el,'b',label='el')
#az,=plt.plot(t,az,'k',label='az')
#rev,=plt.plot(t,rev,'r--',label='rev')
plt.plot(t, var1, 'b-', linewidth=2)
plt.ylabel('%s (deg)' % var)
plt.xlabel('Time(minute)')
#plt.legend(handles=[el,az,rev])
plt.show()
def plot_h5(var, year, month, day, st_hour, st_minute, ed_hour, ed_minute):
open_folder(month, day, year)
files = select_file(st_hour, st_minute, ed_hour, ed_minute)
numb = len(files)
el = []
az = []
rev = []
for fname in files:
with h5.File(fname, 'r') as f:
var1 = f['dataset'][
0, :] #change dataset link 'dataset' according to your hdf5 configuration
var2 = f['dataset'][1, :]
var3 = f['dataset'][2, :]
el.append(var1)
az.append(var2)
rev.append(var3)
size = len(el)
vars = dict(el=np.array(el).reshape(size, 1),
az=np.array(az).reshape(size, 1),
rev=np.array(rev).reshape(size, 1))
t = np.linspace(0.0, numb, size)
#change variable name (y1,y2,y3) to anything you want
#y1,=plt.plot(t,el,'b',label='y1')
#y2,=plt.plot(t,az,'k',label='y2')
#y3,=plt.plot(t,rev,'r--',label='y3')
plt.plot(t, vars[var], 'b-', linewidth=2)
plt.xlabel('Time(minute)')
plt.ylabel('%s' % var)
#plt.legend(handles=[y1,y2,y3])
plt.show()
def plot_h5(var, year, month, day, st_hour, st_minute, ed_hour, ed_minute):
open_folder(month, day, year)
files = select_file(month, day, year, st_hour, st_minute, ed_hour,
ed_minute)
m = conti_test(year, month, day, st_hour, st_minute, ed_hour, ed_minute)
numb = len(files)
size = 0
i = 0
d_h = st_hour - ed_hour
d_m = st_hour - ed_minute
for fname in files:
with h5.File(fname, 'r') as f:
var1 = f['data']['%s' % var]
## var1 = f['data']['%s' % var]
#el=f['data']['el']
#az=f['data']['az']
t = f['data']['gpstime']
size = len(var1)
# print min(t), max(t)
#t=np.linspace(int(m[i-1]),1+int(m[i-1]),size)
#el,=plt.plot(t,el,'b',label='el')
#az,=plt.plot(t,az,'k',label='az')
#rev,=plt.plot(t,rev,'r--',label='rev')
izero = np.where(t != 0)[0]
t = t[izero]
var1 = var1[izero]
plt.plot(t, var1, 'b-', linewidth=2)
plt.ylabel('%s (deg)' % var)
plt.xlabel('gps time')
#plt.legend(handles=[el,az,rev])
plt.show()