wr_integradorv2.py
214 lines
| 7.7 KiB
| text/x-python
|
PythonLexer
| r1350 | import os,numpy,h5py | |||
| from shutil import copyfile | ||||
| import sys,time | ||||
| def isNumber(str): | ||||
| try: | ||||
| float(str) | ||||
| return True | ||||
| except: | ||||
| return False | ||||
| def getfirstFilefromPath(path,meta,ext): | ||||
| validFilelist = [] | ||||
| fileList = os.listdir(path) | ||||
| if len(fileList)<1: | ||||
| return None | ||||
| # meta 1234 567 8-18 BCDE | ||||
| # H,D,PE YYYY DDD EPOC .ext | ||||
| for thisFile in fileList: | ||||
| if meta =="PE": | ||||
| try: | ||||
| number= int(thisFile[len(meta)+7:len(meta)+17]) | ||||
| except: | ||||
| print("There is a file or folder with different format") | ||||
| if meta == "D": | ||||
| try: | ||||
| number= int(thisFile[8:11]) | ||||
| except: | ||||
| print("There is a file or folder with different format") | ||||
| if not isNumber(str=number): | ||||
| continue | ||||
| if (os.path.splitext(thisFile)[-1].lower() != ext.lower()): | ||||
| continue | ||||
| validFilelist.sort() | ||||
| validFilelist.append(thisFile) | ||||
| if len(validFilelist)>0: | ||||
| validFilelist = sorted(validFilelist,key=str.lower) | ||||
| return validFilelist | ||||
| return None | ||||
| def gettimeutcfromDirFilename(path,file): | ||||
| dir_file= path+"/"+file | ||||
| fp = h5py.File(dir_file,'r') | ||||
| epoc = fp['Metadata'].get('utctimeInit')[()] | ||||
| fp.close() | ||||
| return epoc | ||||
| def getDatavaluefromDirFilename(path,file,value): | ||||
| dir_file= path+"/"+file | ||||
| fp = h5py.File(dir_file,'r') | ||||
| array = fp['Data'].get(value)[()] | ||||
| fp.close() | ||||
| return array | ||||
| r1355 | opt = input ("Ingresa Modo de integracion: ") | |||
| if opt==1: | ||||
| #·········· Velocidad de Pedestal················· | ||||
| w = input ("Ingresa velocidad de Pedestal: ") | ||||
| w = 4 | ||||
| w = float(w) | ||||
| #·········· Resolucion minimo en grados··········· | ||||
| alfa = input ("Ingresa resolucion minima en grados: ") | ||||
| alfa = 1 | ||||
| alfa = float(alfa) | ||||
| #·········· IPP del Experimento ·················· | ||||
| IPP = input ("Ingresa el IPP del experimento: ") | ||||
| IPP = 0.0004 | ||||
| IPP = float(IPP) | ||||
| #·········· MODE ·················· | ||||
| mode = input ("Ingresa el MODO del experimento Time(1) or F(0): ") | ||||
| mode = 1 | ||||
| mode = int(mode) | ||||
| time_Interval_p=0.01 | ||||
| n_perfiles_p = 100 | ||||
| blocksPerFile = 100 | ||||
| tiempo_file_1_adq=25 | ||||
| tiempo_file_1_ped=1 | ||||
| r1350 | ||||
| r1355 | else: | |||
| from readFileconfig import ReadfileWR | ||||
| filename= "/home/developer/Downloads/config_WR.txt" | ||||
| dict= ReadfileWR(filename).getDict() | ||||
| mode = dict['mode'] | ||||
| w = dict['vel_ped_azi'] | ||||
| alfa= dict['resolution'] | ||||
| IPP = dict['ipp']*1e-6 | ||||
| time_Interval_p = dict['t_s_ped'] | ||||
| n_perfiles_p = dict['len_ped'] | ||||
| blocksPerFile = int(dict['b_f_adq']) | ||||
| tiempo_file_1_adq= dict['t_f_adq'] | ||||
| tiempo_file_1_ped= dict['t_f_ped'] | ||||
| r1350 | ||||
| #·········· Tiempo en generar la resolucion min··· | ||||
| #············ MCU ·· var_ang = w * (var_tiempo)··· | ||||
| var_tiempo = alfa/w | ||||
| #·········· Tiempo Equivalente en perfiles········ | ||||
| #·········· var_tiempo = IPP * ( num_perfiles )· | ||||
| num_perfiles = int(var_tiempo/IPP) | ||||
| #··········DATA PEDESTAL·························· | ||||
| r1352 | dir_pedestal = "/home/developer/Downloads/Pedestal/P2021093" | |||
| r1350 | #·········· DATA ADQ······························ | |||
| r1355 | if mode==1: | |||
| r1352 | dir_adq = "/home/developer/Downloads/HDF5_TESTPP2V3/d2021093" # Time domain | |||
| r1350 | else: | |||
| dir_adq = "/home/developer/Downloads/hdf5_test/d2021053" # Frequency domain | ||||
| print( "Velocidad angular :", w) | ||||
| print( "Resolucion minima en grados :", alfa) | ||||
| print( "Numero de perfiles equivalente:", num_perfiles) | ||||
| print( "Mode :", mode) | ||||
| #············ First File············· | ||||
| list_pedestal = getfirstFilefromPath(path=dir_pedestal,meta="PE",ext=".hdf5") | ||||
| list_adq = getfirstFilefromPath(path=dir_adq ,meta="D",ext=".hdf5") | ||||
| print("list_pedestal") | ||||
| #print(list_pedestal) | ||||
| print("list_adq") | ||||
| #print(list_adq) | ||||
| #sys.exit(0) | ||||
| #············ utc time ·············· | ||||
| utc_pedestal= gettimeutcfromDirFilename(path=dir_pedestal,file=list_pedestal[0]) | ||||
| utc_adq = gettimeutcfromDirFilename(path=dir_adq ,file=list_adq[0]) | ||||
| print("utc_pedestal :",utc_pedestal) | ||||
| print("utc_adq :",utc_adq) | ||||
| #sys.exit(0) | ||||
| #·············Relacion: utc_adq (+/-) var_tiempo*nro_file= utc_pedestal | ||||
| r1355 | time_Interval_p = time_Interval_p#0.01 | |||
| n_perfiles_p = n_perfiles_p #100 muestras por achivo del pedestal | ||||
| r1350 | ||||
| if utc_adq>utc_pedestal: | ||||
| nro_file = int((utc_adq - utc_pedestal)/(time_Interval_p*n_perfiles_p)) | ||||
| print("nro_file",nro_file) | ||||
| ff_pedestal = list_pedestal[nro_file] | ||||
| print(ff_pedestal) | ||||
| utc_pedestal = gettimeutcfromDirFilename(path=dir_pedestal,file=ff_pedestal) | ||||
| print(utc_pedestal) | ||||
| nro_key_p = int((utc_adq-utc_pedestal)/time_Interval_p) | ||||
| print(nro_key_p) | ||||
| #sys.exit(0) | ||||
| if utc_adq >utc_pedestal: | ||||
| ff_pedestal = ff_pedestal | ||||
| else: | ||||
| nro_file = nro_file-1 | ||||
| ff_pedestal = list_pedestal[nro_file] | ||||
| utc_pedestal = gettimeutcfromDirFilename(path=dir_pedestal,file=ff_pedestal) | ||||
| angulo = getDatavaluefromDirFilename(path=dir_pedestal,file=ff_pedestal,value="azimuth") | ||||
| nro_key_p = int((utc_adq-utc_pedestal)/time_Interval_p) | ||||
| print("nro_file :",nro_file) | ||||
| print("name_file :",ff_pedestal) | ||||
| print("utc_pedestal_file :",utc_pedestal) | ||||
| print("nro_key_p :",nro_key_p) | ||||
| print("utc_pedestal_init :",utc_pedestal+nro_key_p*time_Interval_p) | ||||
| print("angulo_array :",angulo[nro_key_p]) | ||||
| #4+25+25+25+21 | ||||
| #while True: | ||||
| list_pedestal = getfirstFilefromPath(path=dir_pedestal,meta="PE",ext=".hdf5") | ||||
| list_adq = getfirstFilefromPath(path=dir_adq ,meta="D",ext=".hdf5") | ||||
| nro_file = nro_file #10 | ||||
| nro_key_perfil = nro_key_p | ||||
| r1355 | blocksPerFile = blocksPerFile#100##### aqui se cambia dependiendo de los blqoues por achivo en adq | |||
| r1350 | wr_path = "/home/developer/Downloads/HDF5_WR/" | |||
| # Lectura de archivos de adquisicion para adicion de azimuth | ||||
| # factor de archivos | ||||
| #f_a_p = tiempo_file_1_adq/tiempo_file_1_ped=25/1 = 25 | ||||
| r1355 | tiempo_file_1_adq = tiempo_file_1_adq | |||
| tiempo_file_1_ped = tiempo_file_1_ped | ||||
| r1350 | f_a_p= int(tiempo_file_1_adq/tiempo_file_1_ped) | |||
| r1355 | ||||
| r1350 | for thisFile in range(len(list_adq)): | |||
| print("thisFileAdq",thisFile) | ||||
| angulo_adq = numpy.zeros(blocksPerFile) | ||||
| tmp = 0 | ||||
| for j in range(blocksPerFile): | ||||
| iterador = nro_key_perfil + f_a_p*(j-tmp) | ||||
| #print("iterador",iterador) | ||||
| if iterador < n_perfiles_p: | ||||
| nro_file = nro_file | ||||
| else: | ||||
| nro_file = nro_file+1 | ||||
| dif = blocksPerFile-(nro_key_perfil+f_a_p*(j-tmp-1)) | ||||
| tmp = j | ||||
| nro_key_perfil= f_a_p-dif | ||||
| iterador = nro_key_perfil | ||||
| #print(iterador) | ||||
| print("nro_file",nro_file) | ||||
| ff_pedestal = list_pedestal[nro_file] | ||||
| angulo = getDatavaluefromDirFilename(path=dir_pedestal,file=ff_pedestal,value="azimuth") | ||||
| angulo_adq[j]= angulo[iterador] | ||||
| copyfile(dir_adq+"/"+list_adq[thisFile],wr_path+list_adq[thisFile]) | ||||
| fp = h5py.File(wr_path+list_adq[thisFile],'a') | ||||
| #grp = fp.create_group("Pedestal") | ||||
| grp = fp['Metadata'] | ||||
| #sgrp = grp.create_group('Pedestal') | ||||
| dset = grp.create_dataset("azimuth" , data=angulo_adq) | ||||
| fp.close() | ||||
| print("Angulo",angulo_adq) | ||||
| print("Angulo",len(angulo_adq)) | ||||
| nro_key_perfil=iterador + f_a_p | ||||
| if nro_key_perfil< n_perfiles_p: | ||||
| nro_file = nro_file | ||||
| else: | ||||
| nro_file = nro_file+1 | ||||
| nro_key_perfil= nro_key_p | ||||
