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()