schain Files · schainpy/scripts/sousy.py

Se comenta el retorno para evitar la interrupción del programa debido al uso del...

Se comenta el retorno para evitar la interrupción del programa debido al uso del header antiguo, en la ultima version de Signal Chain se genera un error

Miguel Valdez - - Load All Authors

File last commit:

r635:f36ea42cad01


                r725:55bd4dcd7f84

Download file

             sousy.py
        
                    104 lines
            
             | 4.4 KiB
            
                | text/x-python
            
             |
                PythonLexer
            
             / schainpy / scripts / sousy.py
          
                    History
                
                 |
                  Source
                 | Raw
                 |Copy content
                 |Copy permalink

        Miguel Valdez
    
Script updated

              r635
            
      #!python

      '''

      Created on Jul 7, 2015

        Daniel Valdez
    
script to processing of sousy radar data

              r483
            
        Miguel Valdez
    
Script updated

              r635
            
      @author: Miguel Urco

      '''

        Daniel Valdez
    
script to processing of sousy radar data

              r483
            
      import os, sys

        Miguel Valdez
    
Script updated

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

      path = os.path.dirname(path)

      sys.path.insert(0, path)

        Daniel Valdez
    
script to processing of sousy radar data

              r483
            
        Miguel Valdez
    
Script updated

              r635
            
      from schainpy.controller import Project

        Daniel Valdez
    
script to processing of sousy radar data

              r483
            
      desc = "Sousy_test"

      filename = "sousy_processing.xml"

      controllerObj = Project()

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

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

      #path2= '/media/'

      #path2='/media/New Volume/LowTroposphere'

      #path1='/media/New Volume/LT_shortpulse'

      #path = path1 + ',' + path2

      path='G:\\LowTroposphere'

      path = '/media/signalchain/FVillanuevaR/LowTroposphere'

      wr_path = '/media/signalchain/datos/sousy'

      figures_path = '/home/signalchain/Pictures/sousy'

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

                                                  path=path,

                                                  startDate='2014/07/08',

                                                  endDate='2014/07/08',

                                                  startTime='10:00:00',

                                                  endTime='17:59:59',

                                                  delay=0,

                                                  set=0,

                                                  online=0,

                                                  walk=1)

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

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

      ################ SOUSY###################################

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

      #

      procUnitConfObjSousy = controllerObj.addProcUnit(datatype='Voltage', 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 = procUnitConfObjSousy.addOperation(name='setRadarFrequency')

      opObj11.addParameter(name='frequency', value='53.5e6', format='float')

      opObj11 =  procUnitConfObjSousy.addOperation(name='filterByHeights')

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

      codigo='1,-1'

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

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

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

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

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

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

      procUnitConfObjSousySpectra = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjSousy.getId())

      procUnitConfObjSousySpectra.addParameter(name='nFFTPoints', value='64', format='int')

      procUnitConfObjSousySpectra.addParameter(name='nProfiles', value='64', format='int')

      opObj13 = procUnitConfObjSousySpectra.addOperation(name='removeDC')

      opObj13.addParameter(name='mode', value='2', format='int')

      opObj11 = procUnitConfObjSousySpectra.addOperation(name='IncohInt', optype='other')

      opObj11.addParameter(name='n', value='1', format='float')

      # 

      # opObj11 = procUnitConfObjSousySpectra.addOperation(name='RTIPlot', optype='other')

      # opObj11.addParameter(name='id', value='101', format='int')

      # opObj11.addParameter(name='wintitle', value='Sousy_RTIPlot', format='str')

      # opObj11.addParameter(name='zmin', value='30', format='int')

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

      # opObj11.addParameter(name='ymin', value='0', format='int')

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

      # opObj11.addParameter(name='xmin', value='10', format='float')

      # opObj11.addParameter(name='xmax', value='18', format='float')

      # opObj11.addParameter(name='showprofile', value='0', format='int')

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

      # #opObj11.addParameter(name='figfile', value='rti0_sousy.png', format='str')

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

      opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraWriter', optype='other')

      opObj11.addParameter(name='path', value=wr_path)

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

      print "Escribiendo el archivo XML"

      controllerObj.writeXml(filename)

      print "Leyendo el archivo XML"

      controllerObj.readXml(filename)

      controllerObj.createObjects()

      controllerObj.connectObjects()

      controllerObj.run()

	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

Miguel Valdez Script updated	r635	#!python
		'''
		Created on Jul 7, 2015
Daniel Valdez script to processing of sousy radar data	r483
Miguel Valdez Script updated	r635	@author: Miguel Urco
		'''
Daniel Valdez script to processing of sousy radar data	r483	import os, sys

Miguel Valdez Script updated	r635	path = os.path.dirname(os.getcwd())
		path = os.path.dirname(path)
		sys.path.insert(0, path)
Daniel Valdez script to processing of sousy radar data	r483
Miguel Valdez Script updated	r635	from schainpy.controller import Project
Daniel Valdez script to processing of sousy radar data	r483
		desc = "Sousy_test"
		filename = "sousy_processing.xml"

		controllerObj = Project()

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

		#path = '/media/data/data/vientos/57.2063km/echoes/NCO_Woodman'
		#path2= '/media/'
		#path2='/media/New Volume/LowTroposphere'
		#path1='/media/New Volume/LT_shortpulse'
		#path = path1 + ',' + path2
		path='G:\\LowTroposphere'

		path = '/media/signalchain/FVillanuevaR/LowTroposphere'
		wr_path = '/media/signalchain/datos/sousy'
		figures_path = '/home/signalchain/Pictures/sousy'

		readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage',
		path=path,
		startDate='2014/07/08',
		endDate='2014/07/08',
		startTime='10:00:00',
		endTime='17:59:59',
		delay=0,
		set=0,
		online=0,
		walk=1)

		opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
		#########################################################
		################ SOUSY###################################
		#########################################################
		#
		procUnitConfObjSousy = controllerObj.addProcUnit(datatype='Voltage', 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 = procUnitConfObjSousy.addOperation(name='setRadarFrequency')
		opObj11.addParameter(name='frequency', value='53.5e6', format='float')



		opObj11 = procUnitConfObjSousy.addOperation(name='filterByHeights')
		opObj11.addParameter(name='window', value='2', format='int')

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

		opObj11 = procUnitConfObjSousy.addOperation(name='CohInt', optype='other')
		opObj11.addParameter(name='n', value='2048', format='int')

		procUnitConfObjSousySpectra = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjSousy.getId())
		procUnitConfObjSousySpectra.addParameter(name='nFFTPoints', value='64', format='int')
		procUnitConfObjSousySpectra.addParameter(name='nProfiles', value='64', format='int')

		opObj13 = procUnitConfObjSousySpectra.addOperation(name='removeDC')
		opObj13.addParameter(name='mode', value='2', format='int')

		opObj11 = procUnitConfObjSousySpectra.addOperation(name='IncohInt', optype='other')
		opObj11.addParameter(name='n', value='1', format='float')
		#
		# opObj11 = procUnitConfObjSousySpectra.addOperation(name='RTIPlot', optype='other')
		# opObj11.addParameter(name='id', value='101', format='int')
		# opObj11.addParameter(name='wintitle', value='Sousy_RTIPlot', format='str')
		# opObj11.addParameter(name='zmin', value='30', format='int')
		# opObj11.addParameter(name='zmax', value='100', format='int')
		# opObj11.addParameter(name='ymin', value='0', format='int')
		# opObj11.addParameter(name='ymax', value='10', format='int')
		# opObj11.addParameter(name='xmin', value='10', format='float')
		# opObj11.addParameter(name='xmax', value='18', format='float')
		# opObj11.addParameter(name='showprofile', value='0', format='int')
		# opObj11.addParameter(name='save', value='1', format='int')
		# #opObj11.addParameter(name='figfile', value='rti0_sousy.png', format='str')
		# opObj11.addParameter(name='figpath', value=figures_path, format='str')

		opObj11 = procUnitConfObjSousySpectra.addOperation(name='SpectraWriter', optype='other')
		opObj11.addParameter(name='path', value=wr_path)
		opObj11.addParameter(name='blocksPerFile', value='100', format='int')

		print "Escribiendo el archivo XML"
		controllerObj.writeXml(filename)
		print "Leyendo el archivo XML"
		controllerObj.readXml(filename)

		controllerObj.createObjects()
		controllerObj.connectObjects()
		controllerObj.run()