WindProfiler_SA01.py
128 lines
| 5.7 KiB
| text/x-python
|
PythonLexer
|
r502 | # DIAS 19 Y 20 FEB 2014 | |
# Comprobacion de Resultados DBS con SA | |||
import os, sys | |||
path = os.path.split(os.getcwd())[0] | |||
|
r763 | path = os.path.split(path)[0] | |
|
r502 | ||
|
r763 | sys.path.insert(0, path) | |
from schainpy.controller import Project | |||
|
r502 | ||
desc = "SA Experiment Test" | |||
filename = "SA2014050.xml" | |||
controllerObj = Project() | |||
controllerObj.setup(id = '191', name='test01', description=desc) | |||
#Experimentos | |||
#2014050 19 Feb 2014 | |||
|
r763 | path = '/media/joscanoa/84A65E64A65E5730/soporte/Data/MST/SA/d2014050' | |
pathFigure = '/media/joscanoa/84A65E64A65E5730/soporte/workspace/Graficos/SA/prueba1/' | |||
|
r511 | xmin = '15.5' | |
xmax = '24' | |||
startTime = '15:30:00' | |||
filehdf5 = "SA_2014050.hdf5" | |||
|
r502 | ||
#2014051 20 Feb 2014 | |||
|
r511 | # path = '/home/soporte/Data/MST/SA/d2014051' | |
# pathFigure = '/home/soporte/workspace/Graficos/SA/new/' | |||
# xmin = '0.0' | |||
# xmax = '8.0' | |||
# startTime = '00:00:00' | |||
# filehdf5 = "SA_2014051.hdf5" | |||
|
r502 | ||
readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader', | |||
path=path, | |||
startDate='2014/01/01', | |||
endDate='2014/03/31', | |||
startTime=startTime, | |||
endTime='23:59:59', | |||
online=0, | |||
delay=5, | |||
|
r855 | walk=0, | |
getblock=1, | |||
blocksize=32768) | |||
|
r502 | ||
opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock') | |||
#-------------------------------------------------------------------------------------------------- | |||
procUnitConfObj0 = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId()) | |||
opObj11 = procUnitConfObj0.addOperation(name='Decoder', optype='other') | |||
opObj11 = procUnitConfObj0.addOperation(name='CohInt', optype='other') | |||
|
r855 | # opObj11.addParameter(name='n', value='600', format='int') | |
opObj11.addParameter(name='n', value='256', format='int') | |||
|
r502 | ||
opObj11 = procUnitConfObj0.addOperation(name='selectHeightsByIndex') | |||
opObj11.addParameter(name='minIndex', value='10', format='float') | |||
opObj11.addParameter(name='maxIndex', value='60', format='float') | |||
#--------------------------------------------------------------------------------------------------- | |||
procUnitConfObj1 = controllerObj.addProcUnit(datatype='CorrelationProc', inputId=procUnitConfObj0.getId()) | |||
|
r855 | procUnitConfObj1.addParameter(name='pairsList', value='(0,0),(1,1),(2,2),(3,3),(1,0),(2,3)', format='pairsList') | |
# procUnitConfObj1.addParameter(name='removeDC', value='1', format='bool') | |||
# #procUnitConfObj1.addParameter(name='lagT', value='0,1,2,3', format='intlist') | |||
# | |||
# opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other') | |||
# opObj12.addParameter(name='id', value='1', format='int') | |||
# opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str') | |||
# opObj12.addParameter(name='save', value='1', format='bool') | |||
# opObj12.addParameter(name='zmin', value='0', format='int') | |||
# opObj12.addParameter(name='zmax', value='1', format='int') | |||
# opObj12.addParameter(name='figpath', value = pathFigure, format='str') | |||
# | |||
# opObj12 = procUnitConfObj1.addOperation(name='removeNoise') | |||
# opObj12.addParameter(name='mode', value='2', format='int') | |||
# opObj12 = procUnitConfObj1.addOperation(name='calculateNormFactor') | |||
# | |||
# opObj12 = procUnitConfObj1.addOperation(name='CorrelationPlot', optype='other') | |||
# opObj12.addParameter(name='id', value='2', format='int') | |||
# opObj12.addParameter(name='wintitle', value='CrossCorrelation Plot', format='str') | |||
# opObj12.addParameter(name='save', value='1', format='bool') | |||
# opObj12.addParameter(name='zmin', value='0', format='int') | |||
# opObj12.addParameter(name='zmax', value='1', format='int') | |||
# opObj12.addParameter(name='figpath', value = pathFigure, format='str') | |||
# | |||
# #--------------------------------------------------------------------------------------------------- | |||
|
r502 | procUnitConfObj2 = controllerObj.addProcUnit(datatype='ParametersProc', inputId=procUnitConfObj1.getId()) | |
|
r855 | opObj20 = procUnitConfObj2.addOperation(name='SALags', optype='other') | |
# | |||
|
r502 | opObj21 = procUnitConfObj2.addOperation(name='WindProfiler', optype='other') | |
opObj21.addParameter(name='technique', value='SA', format='str') | |||
|
r855 | # # opObj21.addParameter(name='correctFactor', value='-1', format='float') | |
|
r502 | opObj21.addParameter(name='positionX', value='36,0,36,0', format='floatlist') | |
opObj21.addParameter(name='positionY', value='36,0,0,36', format='floatlist') | |||
opObj21.addParameter(name='azimuth', value='51.06', format='float') | |||
|
r855 | ||
# opObj22 = procUnitConfObj2.addOperation(name='WindProfilerPlot', optype='other') | |||
# opObj22.addParameter(name='id', value='4', format='int') | |||
# opObj22.addParameter(name='wintitle', value='Wind Profiler', format='str') | |||
|
r502 | # opObj22.addParameter(name='save', value='1', format='bool') | |
# opObj22.addParameter(name='figpath', value = pathFigure, format='str') | |||
|
r855 | # opObj22.addParameter(name='zmin', value='-15', format='int') | |
# opObj22.addParameter(name='zmax', value='15', format='int') | |||
# opObj22.addParameter(name='zmin_ver', value='-80', format='float') | |||
# opObj22.addParameter(name='zmax_ver', value='80', format='float') | |||
# opObj22.addParameter(name='SNRmin', value='-20', format='int') | |||
# opObj22.addParameter(name='SNRmax', value='40', format='int') | |||
# opObj22.addParameter(name='SNRthresh', value='-3.5', format='float') | |||
# opObj22.addParameter(name='xmin', value=xmin, format='float') | |||
# opObj22.addParameter(name='xmax', value=xmax, format='float') | |||
|
r502 | ||
#----------------------------------------------------------------------------------- | |||
print "Escribiendo el archivo XML" | |||
controllerObj.writeXml(filename) | |||
print "Leyendo el archivo XML" | |||
controllerObj.readXml(filename) | |||
controllerObj.createObjects() | |||
controllerObj.connectObjects() | |||
controllerObj.run() |