wr_marzo_v1_001.py
259 lines
| 11.8 KiB
| text/x-python
|
PythonLexer
avaldezp
|
r1426 | # 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 | ||||
|
avaldezp
|
r1427 | V = 6 | ||
| samp_rate = 10# VERIFICAR | ||||
|
avaldezp
|
r1426 | 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------------------------------- | ||||
|
avaldezp
|
r1427 | #ath = "/DATA_RM/TEST_MARTES_22_1M_1us" | ||
| #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-171722" | ||||
|
avaldezp
|
r1428 | #path = "/DATA_RM/DRONE01ABRIL" | ||
|
avaldezp
|
r1427 | #path = "/DATA_RM/DRONE01ABRIL1429" | ||
| #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-171722" | ||||
|
avaldezp
|
r1428 | #path = "/DATA_RM/DRONE01ABRIL1701" | ||
|
avaldezp
|
r1432 | ##path = "/DATA_RM/DATA/Torre_con_bola_1649092242/rawdata" | ||
| ##path="/DATA_RM/DRONE01ABRIL1727" | ||||
|
avaldezp
|
r1433 | #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" | ||||
|
avaldezp
|
r1432 | #path="/DATA_RM/TEST172956_0411" | ||
| #path = "/DATA_RM/DATA/TEST@2022-04-11T17:29:56/rawdata" | ||||
|
avaldezp
|
r1429 | #path_ped = "/DATA_RM/DRONE01ABRIL1450" | ||
|
avaldezp
|
r1432 | #path_ped="/DATA_RM/TEST_PEDESTAL/P20220401-172744" | ||
|
avaldezp
|
r1433 | #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" | ||||
|
avaldezp
|
r1432 | |||
|
avaldezp
|
r1429 | #path_ped = "/DATA_RM/DATA/Torre_con_bola_1649092242/position/2022-04-04T17-00-00" | ||
|
avaldezp
|
r1426 | #------------------------------------------------------------------------------- | ||
|
avaldezp
|
r1433 | figpath_pp = "/home/soporte/Pictures/PEDMARTES_NEW" | ||
|
avaldezp
|
r1427 | #figpath_pp = "/home/soporte/Pictures/MARTES_22_PP_1M_1us" | ||
| figpath_spec = "/home/soporte/Pictures/MARTES_22_1M_1us" | ||||
|
avaldezp
|
r1432 | figpath_pp_ppi = "/home/soporte/Pictures/PPILUNES11042022" | ||
|
avaldezp
|
r1429 | |||
|
avaldezp
|
r1433 | figpath_pp_rhi = "/DATA_RM/MIERCOLES13" | ||
|
avaldezp
|
r1426 | #--------------------------OPCIONES--------------------------------------------- | ||
|
avaldezp
|
r1433 | plot_ppi = 0 | ||
|
avaldezp
|
r1432 | plot = 0#0 | ||
|
avaldezp
|
r1433 | plot_rhi = 1#1 | ||
|
avaldezp
|
r1432 | integration = 1#1 | ||
|
avaldezp
|
r1426 | 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 ) | ||||
|
avaldezp
|
r1427 | time.sleep(4) | ||
|
avaldezp
|
r1426 | #remotefolder = "/home/wmaster/graficos" | ||
| ################# RANGO DE PLOTEO###################################### | ||||
| dBmin = '20' | ||||
|
avaldezp
|
r1427 | dBmax = '60' | ||
|
avaldezp
|
r1433 | xmin = '14.0' #17.1,17.5 | ||
| xmax = '14.5' #17.2,17.8 | ||||
|
avaldezp
|
r1426 | ymin = '0' #### PONER A 0 | ||
|
avaldezp
|
r1432 | ymax = '1.5' #### PONER A 8 | ||
|
avaldezp
|
r1426 | ########################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, | ||||
|
avaldezp
|
r1433 | startDate="2022/04/13",#today, | ||
| endDate="2022/04/13",#today, | ||||
| startTime='17:37:00',#'17:39:25', | ||||
|
avaldezp
|
r1426 | 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()) | ||||
|
avaldezp
|
r1427 | ''' | ||
| 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') | ||||
|
avaldezp
|
r1432 | opObj11.addParameter(name='h0', value='-1.0', format='float') | ||
|
avaldezp
|
r1426 | |||
| 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 | ||||
|
avaldezp
|
r1428 | opObj11.addParameter(name='maxIndex', value=str(int(num_alturas/20.0)), format='int') | ||
|
avaldezp
|
r1426 | # CUARTA PARTE de 60 Km POR ESO ENTRE 10 - 6 Km | ||
|
avaldezp
|
r1428 | # CUARTA PARTE de 60 Km POR ESO ENTRE 20 - 3 Km | ||
|
avaldezp
|
r1429 | ''' | ||
|
avaldezp
|
r1427 | 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 | ||||
|
avaldezp
|
r1426 | |||
|
avaldezp
|
r1427 | 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') | ||||
|
avaldezp
|
r1429 | ''' | ||
|
avaldezp
|
r1428 | |||
|
avaldezp
|
r1426 | 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) | ||||
|
avaldezp
|
r1427 | opObj11.addParameter(name='ymin', value=ymin) | ||
| opObj11.addParameter(name='ymax', value=ymax) | ||||
|
avaldezp
|
r1426 | 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) | ||||
|
avaldezp
|
r1433 | opObj11.addParameter(name='samp_rate_ped', value='1500',format='int') | ||
| opObj11.addParameter(name='t_Interval_p', value='0.04', format='float') | ||||
|
avaldezp
|
r1429 | #opObj11.addParameter(name='wr_exp', value='PPI') | ||
| opObj11.addParameter(name='wr_exp', value='RHI') | ||||
|
avaldezp
|
r1426 | 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) | ||||
|
avaldezp
|
r1429 | 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') | ||||
|
avaldezp
|
r1434 | opObj11.addParameter(name='ang_max', value='90',format='int') | ||
| opObj11.addParameter(name='ang_min', value='0',format='int') | ||||
|
avaldezp
|
r1429 | opObj11.addParameter(name='save', value=figpath_pp_rhi) | ||
| opObj11.addParameter(name='save_period', value=1) | ||||
|
avaldezp
|
r1426 | |||
| controllerObj.start() | ||||