basic_proc00003.py
184 lines
| 8.5 KiB
| text/x-python
|
PythonLexer
avaldezp
|
r1409 | # Ing. AVP | ||
| # 04/01/2022 | ||||
| # ARCHIVO DE LECTURA | ||||
| import os, sys | ||||
| import datetime | ||||
| import time | ||||
| import numpy | ||||
|
avaldezp
|
r1416 | import json | ||
|
avaldezp
|
r1409 | from ext_met import getfirstFilefromPath,getDatavaluefromDirFilename | ||
| from schainpy.controller import Project | ||||
| #----------------------------------------------------------------------------------------- | ||||
|
avaldezp
|
r1416 | # path_ped = "/DATA_RM/TEST_PEDESTAL/P20211110-171003" | ||
| ## print("PATH PEDESTAL :",path_ped) | ||||
|
avaldezp
|
r1409 | print("[SETUP]-RADAR METEOROLOGICO-") | ||
|
avaldezp
|
r1416 | path_ped = "/DATA_RM/TEST_PEDESTAL/P20211111-173856" | ||
|
avaldezp
|
r1409 | print("PATH PEDESTAL :",path_ped) | ||
|
avaldezp
|
r1416 | path_adq = "/DATA_RM/11" | ||
|
avaldezp
|
r1409 | print("PATH DATA :",path_adq) | ||
|
avaldezp
|
r1416 | |||
|
avaldezp
|
r1409 | figpath_pp_rti = "/home/soporte/Pictures/TEST_PP_RTI" | ||
| print("PATH PP RTI :",figpath_pp_rti) | ||||
| figpath_pp_ppi = "/home/soporte/Pictures/TEST_PP_PPI" | ||||
| print("PATH PP PPI :",figpath_pp_ppi) | ||||
|
avaldezp
|
r1416 | path_pp_save_int = "/DATA_RM/TEST_NEW_FORMAT" | ||
|
avaldezp
|
r1409 | print("PATH SAVE PP INT :",path_pp_save_int) | ||
| print(" ") | ||||
| #------------------------------------------------------------------------------------------- | ||||
| print("SELECCIONAR MODO: PPI (0) O RHI (1)") | ||||
| mode_wr = 0 | ||||
| if mode_wr==0: | ||||
| print("[ ON ] MODE PPI") | ||||
| list_ped = getfirstFilefromPath(path=path_ped,meta="PE",ext=".hdf5") | ||||
| ff_pedestal = list_ped[2] | ||||
| azi_vel = getDatavaluefromDirFilename(path=path_ped,file=ff_pedestal,value="azi_vel") | ||||
| V = round(azi_vel[0]) | ||||
| print("VELOCIDAD AZI :", int(numpy.mean(azi_vel)),"°/seg") | ||||
| else: | ||||
| print("[ ON ] MODE RHI") | ||||
| list_ped = getfirstFilefromPath(path=path_ped,meta="PE",ext=".hdf5") | ||||
| ff_pedestal = list_ped[2] | ||||
| V = round(ele_vel[0]) | ||||
| ele_vel = getDatavaluefromDirFilename(path=path_ped,file=ff_pedestal,value="ele_vel") | ||||
| print("VELOCIDAD ELE :", int(numpy.mean(ele_vel)),"°/seg") | ||||
| print(" ") | ||||
| #--------------------------------------------------------------------------------------- | ||||
| print("SELECCIONAR MODO: PULSE PAIR (0) O FREQUENCY (1)") | ||||
| mode_proc = 0 | ||||
| if mode_proc==0: | ||||
| print("[ ON ] MODE PULSEPAIR") | ||||
| else: | ||||
| print("[ ON ] MODE FREQUENCY") | ||||
| ipp = 60.0 | ||||
| print("IPP(Km.) : %1.2f"%ipp) | ||||
| ipp_sec = (ipp*1.0e3/150.0)*1.0e-6 | ||||
| print("IPP(useg.) : %1.2f"%(ipp_sec*(1.0e6))) | ||||
| VEL=V | ||||
| n= int(1/(VEL*ipp_sec)) | ||||
| print("N° Profiles : ", n) | ||||
| #--------------------------------------------------------------------------------------- | ||||
| plot_rti = 0 | ||||
|
avaldezp
|
r1416 | plot_ppi = 0 | ||
|
avaldezp
|
r1409 | integration = 1 | ||
|
avaldezp
|
r1416 | save = 1 | ||
|
avaldezp
|
r1409 | #---------------------------RANGO DE PLOTEO---------------------------------- | ||
| dBmin = '1' | ||||
| dBmax = '85' | ||||
|
avaldezp
|
r1416 | xmin = '14' | ||
| xmax = '16' | ||||
|
avaldezp
|
r1409 | ymin = '0' | ||
| ymax = '600' | ||||
| #---------------------------------------------------------------------------- | ||||
| time.sleep(3) | ||||
| #---------------------SIGNAL CHAIN ------------------------------------ | ||||
|
avaldezp
|
r1416 | desc_wr= { | ||
| 'Data': { | ||||
| 'dataPP_POW': 'Power', | ||||
| 'utctime': 'Time', | ||||
| 'azimuth': 'az', | ||||
| 'elevation':'el' | ||||
| }, | ||||
| 'Metadata': { | ||||
| 'heightList': 'range', | ||||
| 'channelList': 'Channels' | ||||
| } | ||||
| } | ||||
|
avaldezp
|
r1409 | desc = "USRP_WEATHER_RADAR" | ||
| 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_adq, | ||||
|
avaldezp
|
r1416 | startDate="2021/11/11",#today, | ||
|
avaldezp
|
r1409 | endDate="2021/12/30",#today, | ||
|
avaldezp
|
r1416 | startTime='17:39:17', | ||
|
avaldezp
|
r1409 | endTime='23:59:59', | ||
| delay=0, | ||||
| #set=0, | ||||
| online=0, | ||||
| walk=1, | ||||
| ippKm=ipp) | ||||
| procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc',inputId=readUnitConfObj.getId()) | ||||
| opObj11 = procUnitConfObjA.addOperation(name='selectHeights') | ||||
| opObj11.addParameter(name='minIndex', value='1', format='int') | ||||
| # opObj11.addParameter(name='maxIndex', value='10000', format='int') | ||||
| opObj11.addParameter(name='maxIndex', value='400', format='int') | ||||
| if mode_proc==0: | ||||
| opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other') | ||||
| opObj11.addParameter(name='n', value=int(n), format='int') | ||||
| procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId()) | ||||
|
avaldezp
|
r1411 | # REVISAR EL test_sim00013.py | ||
|
avaldezp
|
r1409 | if plot_rti==1: | ||
| opObj11 = procUnitConfObjB.addOperation(name='GenericRTIPlot',optype='external') | ||||
| opObj11.addParameter(name='attr_data', value='dataPP_POW') | ||||
| opObj11.addParameter(name='colormap', value='jet') | ||||
| opObj11.addParameter(name='xmin', value=xmin) | ||||
| opObj11.addParameter(name='xmax', value=xmax) | ||||
| opObj11.addParameter(name='zmin', value=dBmin) | ||||
| opObj11.addParameter(name='zmax', value=dBmax) | ||||
| opObj11.addParameter(name='save', value=figpath_pp_rti) | ||||
| opObj11.addParameter(name='showprofile', value=0) | ||||
| opObj11.addParameter(name='save_period', value=50) | ||||
| if integration==1: | ||||
| opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation') | ||||
| opObj11.addParameter(name='path_ped', value=path_ped) | ||||
| opObj11.addParameter(name='t_Interval_p', value='0.01', format='float') | ||||
|
avaldezp
|
r1416 | opObj11.addParameter(name='wr_exp', value='PPI') | ||
| #------------------------------------------------------------------------------ | ||||
| ''' | ||||
| opObj11.addParameter(name='Datatype', value='RadialSet') | ||||
| opObj11.addParameter(name='Scantype', value='PPI') | ||||
| opObj11.addParameter(name='Latitude', value='-11.96') | ||||
| opObj11.addParameter(name='Longitud', value='-76.54') | ||||
| opObj11.addParameter(name='Heading', value='293') | ||||
| opObj11.addParameter(name='Height', value='293') | ||||
| opObj11.addParameter(name='Waveform', value='OFM') | ||||
| opObj11.addParameter(name='PRF', value='2000') | ||||
| opObj11.addParameter(name='CreatedBy', value='WeatherRadarJROTeam') | ||||
| opObj11.addParameter(name='ContactInformation', value='avaldez@igp.gob.pe') | ||||
| ''' | ||||
|
avaldezp
|
r1409 | 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 save==1: | ||||
| opObj10 = procUnitConfObjB.addOperation(name='HDFWriter') | ||||
| opObj10.addParameter(name='path',value=path_pp_save_int) | ||||
| opObj10.addParameter(name='mode',value="weather") | ||||
|
avaldezp
|
r1416 | opObj10.addParameter(name='type_data',value='F') | ||
|
avaldezp
|
r1409 | opObj10.addParameter(name='blocksPerFile',value='360',format='int') | ||
|
avaldezp
|
r1416 | #opObj10.addParameter(name='metadataList',value='utctimeInit,paramInterval,channelList,heightList,flagDataAsBlock',format='list') | ||
| opObj10.addParameter(name='metadataList',value='heightList,channelList,Typename,Datatype,Scantype,Latitude,Longitud,Heading,Height,Waveform,PRF,CreatedBy,ContactInformation',format='list') | ||||
| #-------------------- | ||||
| opObj10.addParameter(name='Typename', value='Differential_Reflectivity') | ||||
| opObj10.addParameter(name='Datatype', value='RadialSet') | ||||
| opObj10.addParameter(name='Scantype', value='PPI') | ||||
| opObj10.addParameter(name='Latitude', value='-11.96') | ||||
| opObj10.addParameter(name='Longitud', value='-76.54') | ||||
| opObj10.addParameter(name='Heading', value='293') | ||||
| opObj10.addParameter(name='Height', value='293') | ||||
| opObj10.addParameter(name='Waveform', value='OFM') | ||||
| opObj10.addParameter(name='PRF', value='2000') | ||||
| opObj10.addParameter(name='CreatedBy', value='WeatherRadarJROTeam') | ||||
| opObj10.addParameter(name='ContactInformation', value='avaldez@igp.gob.pe') | ||||
| #--------------------------------------------------- | ||||
| #opObj10.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,azimuth,elevation,utctime',format='list')#,format='list' | ||||
| #opObj10.addParameter(name='metadataList',value='utctimeInit,timeZone,paramInterval,profileIndex,channelList,heightList,flagDataAsBlock',format='list') | ||||
|
avaldezp
|
r1409 | |||
|
avaldezp
|
r1416 | opObj10.addParameter(name='dataList',value='dataPP_POW,azimuth,elevation,utctime',format='list')#,format='list' | ||
| opObj10.addParameter(name='description',value=json.dumps(desc_wr)) | ||||
|
avaldezp
|
r1409 | |||
| controllerObj.start() | ||||