# Ing. AVP # 06/10/2021 # ARCHIVO DE LECTURA import os, sys import datetime import time from schainpy.controller import Project #### NOTA########################################### # INPUT : # VELOCIDAD PARAMETRO : V = 2°/seg # MODO PULSE PAIR O MOMENTOS: 0 : Pulse Pair ,1 : Momentos ###################################################### ##### PROCESAMIENTO ################################## ##### OJO TENER EN CUENTA EL n= para el Pulse Pair ## ##### O EL n= nFFTPoints ### ###################################################### ######## BUSCAMOS EL numero de IPP equivalente 1°##### ######## Sea V la velocidad del Pedestal en °/seg##### ######## 1° sera Recorrido en un tiempo de 1/V ###### ######## IPP del Radar 400 useg --> 60 Km ############ ######## n = 1/(V(°/seg)*IPP(Km)) , NUMERO DE IPP ## ######## n = 1/(V*IPP) ############################# #-------------------------VELOCIDAD DEL PEDESTAL Y MODO ------------------------ print("SETUP- RADAR METEOROLOGICO") IPP = 400*1e-6 V = 6 samp_rate = 10# VERIFICAR MODE_TABLE = 1 # PUEDE SER 1 O 0 AXIS = [1,1,1,1] # AZIMUTH 1 ELEVACION 0 SPEED_AXIS = [10,10,10,10] # VELOCIDAD ANGLE_AXIS = [20,25,30,15] # ANGULOS mode_proc = 0 #-----------------------------PATH ADQ Y PEDESTAL------------------------------- #ath = "/DATA_RM/TEST_MARTES_22_1M_1us" #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-171722" #path = "/DATA_RM/DRONE01ABRIL" #path = "/DATA_RM/DRONE01ABRIL1429" #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-171722" #path = "/DATA_RM/DRONE01ABRIL1701" ##path = "/DATA_RM/DATA/Torre_con_bola_1649092242/rawdata" ##path="/DATA_RM/DRONE01ABRIL1727" #path ="/DATA_RM/DATA/TEST@2022-04-11T17:29:56/rawdata" #path ="/DATA_RM/DATA/TEST@2022-04-13T15:10:42/rawdata" path ="/DATA_RM/DATA/TEST@2022-04-13T17:35:06/rawdata" #path="/DATA_RM/TEST172956_0411" #path = "/DATA_RM/DATA/TEST@2022-04-11T17:29:56/rawdata" #path_ped = "/DATA_RM/DRONE01ABRIL1450" #path_ped="/DATA_RM/TEST_PEDESTAL/P20220401-172744" #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220412-140748" path_ped="/DATA_RM/DATA/TEST@2022-04-13T17:35:06/position/2022-04-13T22-00-00" #path_ped = "/DATA_RM/DATA/TEST@2022-04-13T15:10:42/position/2022-04-13T20-00-00" #path_ped = "/DATA_RM/DATA/TEST@2022-04-12T14:05:11/position/2022-04-12T19-00-00" #path_ped="/DATA_RM/TEST_PEDESTAL/P20220411-173017" #path_ped = "/DATA_RM/DATA/Torre_con_bola_1649092242/position/2022-04-04T17-00-00" #------------------------------------------------------------------------------- figpath_pp = "/home/soporte/Pictures/PEDMARTES_NEW" #figpath_pp = "/home/soporte/Pictures/MARTES_22_PP_1M_1us" figpath_spec = "/home/soporte/Pictures/MARTES_22_1M_1us" figpath_pp_ppi = "/home/soporte/Pictures/PPILUNES11042022" figpath_pp_rhi = "/DATA_RM/MIERCOLES13" #--------------------------OPCIONES--------------------------------------------- plot_ppi = 0 plot = 0#0 plot_rhi = 1#1 integration = 1#1 save = 0 plot_spec = 0 #---------------------------SAVE HDF5 PROCESADO/-------------------------------- if save == 1: if mode_proc==0: path_save = '/DATA_RM/TEST_HDF5_PP_MAR22/6v' else: path_save = '/DATA_RM/TEST_HDF5_SPEC_MAR22/6v' print("[SETUP]-RADAR METEOROLOGICO-") print("* PATH data ADQ :", path) print("* PATH data PED :", path_ped) print("* SAMPLE RATE ADQ Mhz :", samp_rate) print("* Velocidad Pedestal :",V,"°/seg") print("* Configuracion del Pedestal *") print("*** AXIS :",AXIS) print("*** SPEED_AXIS:",SPEED_AXIS) print("*** ANGLE_AXIS:",ANGLE_AXIS) num_alturas = int(samp_rate*IPP*1e6) print("* Nro de Altura :",num_alturas) ############################ NRO Perfiles PROCESAMIENTO ################### V=V n= int(1/(V*IPP)) print("* n - NRO Perfiles Proc:", n ) ################################## MODE ################################### print("* Modo de Operacion :",mode_proc) if mode_proc ==0: print("* Met. Seleccionado : Pulse Pair") else: print("* Met. Momentos : Momentos") ################################## MODE ################################### print("* Grabado de datos :",save) if save ==1: if mode_proc==0: print("[ ON ] MODE PULSEPAIR") else: print("[ ON ] MODE FREQUENCY") print("* Integracion de datos :",integration) print("* Ploteo de datos Parameters:", plot) if plot==1: print("* Path PP plot :", figpath_pp ) if plot_ppi==1: print("* Path PPI plot :", figpath_pp_ppi ) time.sleep(4) #remotefolder = "/home/wmaster/graficos" ################# RANGO DE PLOTEO###################################### dBmin = '20' dBmax = '60' xmin = '14.0' #17.1,17.5 xmax = '14.5' #17.2,17.8 ymin = '0' #### PONER A 0 ymax = '1.5' #### PONER A 8 ########################FECHA########################################## str1 = datetime.date.today() today = str1.strftime("%Y/%m/%d") str2 = str1 - datetime.timedelta(days=1) yesterday = str2.strftime("%Y/%m/%d") #------------------------SIGNAL CHAIN------------------------------------------- desc = "USRP_test" filename = "USRP_processing.xml" controllerObj = Project() controllerObj.setup(id = '191', name='Test_USRP', description=desc) #------------------------ UNIDAD DE LECTURA------------------------------------- readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader', path=path, startDate="2022/04/13",#today, endDate="2022/04/13",#today, startTime='17:37:00',#'17:39:25', endTime='23:59:59',#23:59:59', delay=0, #set=0, online=0, walk=1, ippKm = 60) opObj11 = readUnitConfObj.addOperation(name='printInfo') procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId()) ''' opObj10 = procUnitConfObjA.addOperation(name='ScopePlot', optype='external') opObj10.addParameter(name='id', value='10', format='int') opObj10.addParameter(name='zmin', value='0', format='int') opObj10.addParameter(name='zmax', value='3', format='int') opObj10.addParameter(name='type', value='iq') opObj10.addParameter(name='ymin', value='-1200', format='int') opObj10.addParameter(name='ymax', value='1200', format='int') #opObj10.addParameter(name='save', value=figpath, format='str') opObj10.addParameter(name='save_period', value=10, format='int') ''' opObj11 = procUnitConfObjA.addOperation(name='setH0') opObj11.addParameter(name='h0', value='-1.0', format='float') opObj11 = procUnitConfObjA.addOperation(name='selectHeights') opObj11.addParameter(name='minIndex', value='1', format='int') #opObj11.addParameter(name='maxIndex', value='1000', format='int') #opObj11.addParameter(name='maxIndex', value=str(int(num_alturas/4.0)), format='int') # CUARTA PARTE de 60 Km POR ESO ENTRE 4 - 15 Km opObj11.addParameter(name='maxIndex', value=str(int(num_alturas/20.0)), format='int') # CUARTA PARTE de 60 Km POR ESO ENTRE 10 - 6 Km # CUARTA PARTE de 60 Km POR ESO ENTRE 20 - 3 Km ''' procUnitConfObjB = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjA.getId()) procUnitConfObjB.addParameter(name='nFFTPoints', value='32', format='int') procUnitConfObjB.addParameter(name='nProfiles', value='32', format='int') #SpectraPlot opObj11 = procUnitConfObjB.addOperation(name='SpectraPlot', optype='external') opObj11.addParameter(name='id', value='1', format='int') opObj11.addParameter(name='wintitle', value='Spectra', format='str') #opObj11.addParameter(name='xmin', value=-0.01, format='float') #opObj11.addParameter(name='xmax', value=0.01, format='float') opObj11.addParameter(name='zmin', value=dBmin, format='int') opObj11.addParameter(name='zmax', value=dBmax, format='int') opObj11.addParameter(name='ymin', value=ymin, format='int') opObj11.addParameter(name='ymax', value=ymax, format='int') opObj11.addParameter(name='showprofile', value='1', format='int') #opObj11.addParameter(name='save', value=figpath, format='str') opObj11.addParameter(name='save_period', value=10, format='int') ''' if mode_proc ==0: ####################### METODO PULSE PAIR ###################################################################### opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other') opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO QUE LA VELOCIDAD ES 10 GRADOS #opObj11.addParameter(name='removeDC', value=1, format='int') ####################### METODO Parametros ###################################################################### procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId()) if plot==1: opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external') opObj11.addParameter(name='attr_data', value='dataPP_POWER') opObj11.addParameter(name='colormap', value='jet') opObj11.addParameter(name='xmin', value=xmin) opObj11.addParameter(name='xmax', value=xmax) opObj11.addParameter(name='ymin', value=ymin) opObj11.addParameter(name='ymax', value=ymax) opObj11.addParameter(name='zmin', value=dBmin) opObj11.addParameter(name='zmax', value=dBmax) opObj11.addParameter(name='save', value=figpath_pp) opObj11.addParameter(name='showprofile', value=0) opObj11.addParameter(name='save_period', value=10) ####################### METODO ESCRITURA ####################################################################### if save==1: opObj10 = procUnitConfObjB.addOperation(name='HDFWriter') opObj10.addParameter(name='path',value=path_save) #opObj10.addParameter(name='mode',value=0) opObj10.addParameter(name='blocksPerFile',value='100',format='int') opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list') opObj10.addParhirameter(name='dataList',value='dataPP_POWER,dataPP_DOP,utctime',format='list')#,format='list' if integration==1: opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation') opObj11.addParameter(name='path_ped', value=path_ped) opObj11.addParameter(name='samp_rate_ped', value='1500',format='int') opObj11.addParameter(name='t_Interval_p', value='0.04', format='float') #opObj11.addParameter(name='wr_exp', value='PPI') opObj11.addParameter(name='wr_exp', value='RHI') if plot_ppi==1: opObj11 = procUnitConfObjB.addOperation(name='Block360') opObj11.addParameter(name='n', value='10', format='int') opObj11.addParameter(name='mode', value=mode_proc, format='int') # este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180 opObj11= procUnitConfObjB.addOperation(name='WeatherPlot',optype='other') opObj11.addParameter(name='save', value=figpath_pp_ppi) opObj11.addParameter(name='save_period', value=1) if plot_rhi==1: opObj11 = procUnitConfObjB.addOperation(name='Block360') opObj11.addParameter(name='n', value='10', format='int') opObj11.addParameter(name='mode', value=mode_proc, format='int') # este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180 opObj11= procUnitConfObjB.addOperation(name='WeatherRHIPlot',optype='other') opObj11.addParameter(name='ang_max', value='90',format='int') opObj11.addParameter(name='ang_min', value='0',format='int') opObj11.addParameter(name='save', value=figpath_pp_rhi) opObj11.addParameter(name='save_period', value=1) controllerObj.start()