schain Files · schainpy/scripts/sousy.py

Add SkyMapPlotData, operation can access parent kwargs, fix server plot for multiple ReceiverData

Miguel Valdez - - Load All Authors

File last commit:

r635:f36ea42cad01


                r937:6cdcc42f2694

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