schain Files · schainpy/scripts/USRP_ADQ_PP.py

last 30-16:19

avaldez - - Load All Authors

File last commit:

r1384:2acdbd5998b4


                r1405:49f004b59e80

Download file

             USRP_ADQ_PP.py
        
                    147 lines
            
             | 6.7 KiB
            
                | text/x-python
            
             |
                PythonLexer
            
             / schainpy / scripts / USRP_ADQ_PP.py
          
                    History
                
                 |
                  Source
                 | Raw
                 |Copy content
                 |Copy permalink

        avaldez
    
19 DE AGOSTO 2021 RM

              r1367
            
      #!python

      '''

      '''

      import os, sys

      import datetime

      import time

      #path = os.path.dirname(os.getcwd())

      #path = os.path.dirname(path)

      #sys.path.insert(0, path)

      from schainpy.controller import Project

      desc = "USRP_test"

      filename = "USRP_processing.xml"

      controllerObj = Project()

      controllerObj.setup(id = '191', name='Test_USRP', description=desc)

      ############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############

      #######################################################################

      ######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################

      #######################################################################

      #path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'

      #path = '/DATA_RM/TEST_INTEGRACION'

        david c
    
update de procesamiento y revision de ploteo

              r1378
            
      #path = '/DATA_RM/TEST_ONLINE'

        avaldez
    
update change RM

              r1384
            
      #path = '/DATA_RM/TEST_INTEGRACION/ADQ_OFFLINE/'

      # ULTIMO TEST 22 DE SEPTIEMBRE

      path = '/DATA_RM/USRP_22'

      #path_pp = '/DATA_RM/TEST_HDF5'

      # UTIMO TEST 22 DE SEPTIEMBRE

      path_pp = '/DATA_RM/TEST_HDF5_PP_22'

      ######################################################

      ##### OJO TENER EN CUENTA EL n= para el Pulse Pair ###

      ######################################################

      ######## 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*IPP) , NUMERO DE IPP #############

      ######## n   = 1/(V*IPP) #############################

      V=2

      IPP=400*1e-6

      n= 1/(V*IPP)

      print("n numero de Perfiles a procesar con Pulse Pair: ", n)

        avaldez
    
19 DE AGOSTO 2021 RM

              r1367
            
      figpath = '/home/soporte/Pictures/TEST_INTEGRACION_IMG'

      #remotefolder = "/home/wmaster/graficos"

      #######################################################################

      ################# RANGO DE PLOTEO######################################

      #######################################################################

      dBmin = '-5'

      dBmax = '20'

      xmin = '0'

      xmax ='24'

      ymin = '0'

      ymax = '600'

      #######################################################################

      ########################FECHA##########################################

      #######################################################################

      str = datetime.date.today()

      today = str.strftime("%Y/%m/%d")

      str2 = str - datetime.timedelta(days=1)

      yesterday = str2.strftime("%Y/%m/%d")

      #######################################################################

      ######################## UNIDAD DE LECTURA#############################

      #######################################################################

      readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',

                                                  path=path,

                                                  startDate="2021/01/01",#today,

                                                  endDate="2021/12/30",#today,

                                                  startTime='00:00:00',

                                                  endTime='23:59:59',

                                                  delay=0,

                                                  #set=0,

        david c
    
update de procesamiento y revision de ploteo

              r1378
            
                                                  online=0,

        avaldez
    
19 DE AGOSTO 2021 RM

              r1367
            
                                                  walk=1,

                                                  ippKm = 60)

      opObj11 = readUnitConfObj.addOperation(name='printInfo')

      #opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')

      #######################################################################

      ################ OPERACIONES DOMINIO DEL TIEMPO########################

      #######################################################################

      procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())

      #

      # codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\

      #              '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'

      #opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')

      #opObj11.addParameter(name='frequency', value='70312500')

      opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')

        avaldez
    
update change RM

              r1384
            
      opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO  QUE LA VELOCIDAD ES 10 GRADOS

        avaldez
    
19 DE AGOSTO 2021 RM

              r1367
            
      opObj11.addParameter(name='removeDC', value=1, format='int')

      # Ploteo TEST

      '''

      opObj11 = procUnitConfObjA.addOperation(name='PulsepairPowerPlot', optype='other')

      opObj11 = procUnitConfObjA.addOperation(name='PulsepairSignalPlot', optype='other')

      opObj11 = procUnitConfObjA.addOperation(name='PulsepairVelocityPlot', optype='other')

      #opObj11.addParameter(name='xmax', value=8)

      opObj11 = procUnitConfObjA.addOperation(name='PulsepairSpecwidthPlot', optype='other')

      '''

      # OJO SCOPE

      #opObj10 = procUnitConfObjA.addOperation(name='ScopePlot', optype='external')

      #opObj10.addParameter(name='id', value='10', format='int')

      ##opObj10.addParameter(name='xmin', value='0', format='int')

      ##opObj10.addParameter(name='xmax', value='50', format='int')

      #opObj10.addParameter(name='type', value='iq')

      ##opObj10.addParameter(name='ymin', value='-5000', format='int')

      ##opObj10.addParameter(name='ymax', value='8500', format='int')

      #opObj11.addParameter(name='save', value=figpath, format='str')

      #opObj11.addParameter(name='save_period', value=10, format='int')

      #opObj10 = procUnitConfObjA.addOperation(name='setH0')

      #opObj10.addParameter(name='h0', value='-5000', format='float')

      #opObj11 =  procUnitConfObjA.addOperation(name='filterByHeights')

      #opObj11.addParameter(name='window', value='1', format='int')

      #codigo='1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,1,1,-1,-1,-1'

      #opObj11 = procUnitConfObjSousy.addOperation(name='Decoder', optype='other')

      #opObj11.addParameter(name='code', value=codigo, formatyesterday='floatlist')

      #opObj11.addParameter(name='nCode', value='1', format='int')

      #opObj11.addParameter(name='nBaud', value='28', format='int')

      #opObj11 = procUnitConfObjA.addOperation(name='CohInt', optype='other')

      #opObj11.addParameter(name='n', value='100', format='int')

      #######################################################################

      ########## OPERACIONES ParametersProc########################

      #######################################################################

      procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())

      opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')

      opObj10.addParameter(name='path',value=path_pp)

      #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.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'

      controllerObj.start()

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages

avaldez 19 DE AGOSTO 2021 RM	r1367	#!python
		'''
		'''

		import os, sys
		import datetime
		import time

		#path = os.path.dirname(os.getcwd())
		#path = os.path.dirname(path)
		#sys.path.insert(0, path)

		from schainpy.controller import Project

		desc = "USRP_test"
		filename = "USRP_processing.xml"
		controllerObj = Project()
		controllerObj.setup(id = '191', name='Test_USRP', description=desc)

		############## USED TO PLOT IQ VOLTAGE, POWER AND SPECTRA #############

		#######################################################################
		######PATH DE LECTURA, ESCRITURA, GRAFICOS Y ENVIO WEB#################
		#######################################################################
		#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'
		#path = '/DATA_RM/TEST_INTEGRACION'
david c update de procesamiento y revision de ploteo	r1378	#path = '/DATA_RM/TEST_ONLINE'
avaldez update change RM	r1384	#path = '/DATA_RM/TEST_INTEGRACION/ADQ_OFFLINE/'
		# ULTIMO TEST 22 DE SEPTIEMBRE
		path = '/DATA_RM/USRP_22'
		#path_pp = '/DATA_RM/TEST_HDF5'
		# UTIMO TEST 22 DE SEPTIEMBRE
		path_pp = '/DATA_RM/TEST_HDF5_PP_22'
		######################################################
		##### OJO TENER EN CUENTA EL n= para el Pulse Pair ###
		######################################################
		######## 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*IPP) , NUMERO DE IPP #############
		######## n = 1/(V*IPP) #############################
		V=2
		IPP=400*1e-6
		n= 1/(V*IPP)
		print("n numero de Perfiles a procesar con Pulse Pair: ", n)



avaldez 19 DE AGOSTO 2021 RM	r1367
		figpath = '/home/soporte/Pictures/TEST_INTEGRACION_IMG'
		#remotefolder = "/home/wmaster/graficos"
		#######################################################################
		################# RANGO DE PLOTEO######################################
		#######################################################################
		dBmin = '-5'
		dBmax = '20'
		xmin = '0'
		xmax ='24'
		ymin = '0'
		ymax = '600'
		#######################################################################
		########################FECHA##########################################
		#######################################################################
		str = datetime.date.today()
		today = str.strftime("%Y/%m/%d")
		str2 = str - datetime.timedelta(days=1)
		yesterday = str2.strftime("%Y/%m/%d")
		#######################################################################
		######################## UNIDAD DE LECTURA#############################
		#######################################################################
		readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
		path=path,
		startDate="2021/01/01",#today,
		endDate="2021/12/30",#today,
		startTime='00:00:00',
		endTime='23:59:59',
		delay=0,
		#set=0,
david c update de procesamiento y revision de ploteo	r1378	online=0,
avaldez 19 DE AGOSTO 2021 RM	r1367	walk=1,
		ippKm = 60)

		opObj11 = readUnitConfObj.addOperation(name='printInfo')
		#opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
		#######################################################################
		################ OPERACIONES DOMINIO DEL TIEMPO########################
		#######################################################################

		procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())

		#
		# codigo64='1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,0,0,1,0,1,1,1,0,0,0,1,0,'+\
		# '1,1,1,0,1,1,0,1,1,1,1,0,0,0,1,0,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1,0,0,0,1,0,0,1,0,0,0,0,1,1,1,0,1,1,1,1,0,1,1,0,1,0,0,0,1,1,1,0,1'

		#opObj11 = procUnitConfObjA.addOperation(name='setRadarFrequency')
		#opObj11.addParameter(name='frequency', value='70312500')
		opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
avaldez update change RM	r1384	opObj11.addParameter(name='n', value=int(n), format='int')#10 VOY A USAR 250 DADO QUE LA VELOCIDAD ES 10 GRADOS
avaldez 19 DE AGOSTO 2021 RM	r1367	opObj11.addParameter(name='removeDC', value=1, format='int')
		# Ploteo TEST
		'''
		opObj11 = procUnitConfObjA.addOperation(name='PulsepairPowerPlot', optype='other')
		opObj11 = procUnitConfObjA.addOperation(name='PulsepairSignalPlot', optype='other')
		opObj11 = procUnitConfObjA.addOperation(name='PulsepairVelocityPlot', optype='other')
		#opObj11.addParameter(name='xmax', value=8)
		opObj11 = procUnitConfObjA.addOperation(name='PulsepairSpecwidthPlot', optype='other')
		'''
		# OJO SCOPE
		#opObj10 = procUnitConfObjA.addOperation(name='ScopePlot', optype='external')
		#opObj10.addParameter(name='id', value='10', format='int')
		##opObj10.addParameter(name='xmin', value='0', format='int')
		##opObj10.addParameter(name='xmax', value='50', format='int')
		#opObj10.addParameter(name='type', value='iq')
		##opObj10.addParameter(name='ymin', value='-5000', format='int')
		##opObj10.addParameter(name='ymax', value='8500', format='int')
		#opObj11.addParameter(name='save', value=figpath, format='str')
		#opObj11.addParameter(name='save_period', value=10, format='int')

		#opObj10 = procUnitConfObjA.addOperation(name='setH0')
		#opObj10.addParameter(name='h0', value='-5000', format='float')

		#opObj11 = procUnitConfObjA.addOperation(name='filterByHeights')
		#opObj11.addParameter(name='window', value='1', format='int')

		#codigo='1,1,-1,1,1,-1,1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,-1,-1,-1,1,1,1,1,-1,-1,-1'
		#opObj11 = procUnitConfObjSousy.addOperation(name='Decoder', optype='other')
		#opObj11.addParameter(name='code', value=codigo, formatyesterday='floatlist')
		#opObj11.addParameter(name='nCode', value='1', format='int')
		#opObj11.addParameter(name='nBaud', value='28', format='int')

		#opObj11 = procUnitConfObjA.addOperation(name='CohInt', optype='other')
		#opObj11.addParameter(name='n', value='100', format='int')

		#######################################################################
		########## OPERACIONES ParametersProc########################
		#######################################################################

		procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
		opObj10 = procUnitConfObjB.addOperation(name='HDFWriter')
		opObj10.addParameter(name='path',value=path_pp)
		#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.addParameter(name='dataList',value='dataPP_POW,dataPP_DOP,utctime',format='list')#,format='list'

		controllerObj.start()