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