##// END OF EJS Templates
Encontre un bug durante la escritura de datos en este archivo....
Encontre un bug durante la escritura de datos en este archivo. En la linea 358 jroheaderIO.py Si codetype es igual a 0 no deberia pedirme nCode, nBaud

File last commit:

r472:bb102d4de45c
r556:30d732fca7ac
Show More
mst_isr_eej_processing.py
271 lines | 13.7 KiB | text/x-python | PythonLexer
/ schainpy / test / mst_isr_eej_processing.py
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 import os, sys
path = os.path.split(os.getcwd())[0]
sys.path.append(path)
from controller import *
desc = "EWDrifts Experiment Test"
filename = "mst_isr_eej.xml"
controllerObj = Project()
controllerObj.setup(id = '191', name='test01', description=desc)
path='/remote/ewdrifts/RAW_EXP/EW_DRIFT_FARADAY/EW_Drift'
path1 = '/media/New Volume/DATA/MST_ISR'
path2 = '/media/DATA/DATA/RAW_EXP/MST-EEJ'
path = path1 + ',' + path2
path = '/home/dsuarez/.gvfs/data on 10.10.20.13/MST_ISR_EEJ'
gpath = '/media/datos/pictures/mstisr_mayo2014'
Daniel Valdez
r472 path = '/Volumes/New Volume/mst_isr_eej'
gpath = '/Users/dsuarez/Pictures/poster'
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 readUnitConfObj = controllerObj.addReadUnit(datatype='Voltage',
path=path,
Daniel Valdez
r472 startDate='2014/05/15',
endDate='2014/05/15',
startTime='08:00:00',
endTime='16:00:00',
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 delay=3,
Daniel Valdez
r468 set=0,
Daniel Valdez
r472 online=0,
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 walk=1)
opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
#########################################################
################ MST ####################################
#########################################################
Daniel Valdez
r472 # procUnitConfObjMST = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId())
#
# opObj11 = procUnitConfObjMST.addOperation(name='ProfileSelector', optype='other')
# profileIndex = '0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119'
#
# opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
#
# opObj11 = procUnitConfObjMST.addOperation(name='Decoder', optype='other')
#
# opObj11 = procUnitConfObjMST.addOperation(name='CohInt', optype='other')
# opObj11.addParameter(name='n', value='20', format='int')
#
# procUnitConfObjMSTSpectra = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjMST.getId())
# procUnitConfObjMSTSpectra.addParameter(name='nFFTPoints', value='64', format='int')
# procUnitConfObjMSTSpectra.addParameter(name='nProfiles', value='64', format='int')
#
# opObj11 = procUnitConfObjMSTSpectra.addOperation(name='IncohInt', optype='other')
# opObj11.addParameter(name='n', value='2', format='float')
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467
Daniel Valdez
r468 # opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SpectraWriter', optype='other')
# opObj11.addParameter(name='path', value='/media/datos/mst2014')
# opObj11.addParameter(name='blocksPerFile', value='10', format='int')
Daniel Valdez
r472 # opObj11 = procUnitConfObjMSTSpectra.addOperation(name='RTIPlot', optype='other')
# opObj11.addParameter(name='id', value='1000', format='int')
# opObj11.addParameter(name='wintitle', value='MST', format='str')
# opObj11.addParameter(name='showprofile', value='0', format='int')
# opObj11.addParameter(name='xmin', value='0', format='int')
# opObj11.addParameter(name='xmax', value='24', format='int')
# opObj11.addParameter(name='ymin', value='120', format='int')
# opObj11.addParameter(name='ymax', value='190', format='int')
# opObj11.addParameter(name='zmin', value='20', format='int')
# opObj11.addParameter(name='zmax', value='50', format='int')
# opObj11 = procUnitConfObjMSTSpectra.addOperation(name='RTIPlot', optype='other')
# opObj11.addParameter(name='id', value='101', format='int')
# opObj11.addParameter(name='wintitle', value='MST', format='str')
# opObj11.addParameter(name='showprofile', value='0', format='int')
# opObj11.addParameter(name='channelList', value='3', format='intlist')
# #opObj11.addParameter(name='timerange', value='300', format='float')
# opObj11.addParameter(name='xmin', value='8', format='float')
# opObj11.addParameter(name='xmax', value='16', format='float')
# opObj11.addParameter(name='ymin', value='120', format='float')
# opObj11.addParameter(name='ymax', value='190', format='float')
# opObj11.addParameter(name='zmin', value='20', format='float')
# opObj11.addParameter(name='zmax', value='45', format='float')
# opObj11.addParameter(name='save', value='1', format='int')
# opObj11.addParameter(name='figfile', value='rti_mst.pdf', format='str')
# opObj11.addParameter(name='figpath', value='/Users/dsuarez/Pictures/poster', format='str')
Daniel Valdez
r468
Daniel Valdez
r472 # opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SpectraPlot', optype='other')
# opObj11.addParameter(name='id', value='1001', format='int')
# opObj11.addParameter(name='wintitle', value='MST', format='str')
# opObj11.addParameter(name='ymin', value='120', format='int')
# opObj11.addParameter(name='ymax', value='190', format='int')
# opObj11.addParameter(name='zmin', value='20', format='int')
# opObj11.addParameter(name='zmax', value='45', format='int')
Daniel Valdez
r468
Daniel Valdez
r472 # opObj11.addParameter(name='save', value='1', format='int')
# opObj11.addParameter(name='figpath', value=gpath, format='str')
# opObj11.addParameter(name='wr_period', value='5', format='int')
# opObj11.addParameter(name='ftp', value='1', format='int')
# opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
# opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
# opObj11.addParameter(name='username', value='wmaster', format='str')
# opObj11.addParameter(name='password', value='mst2010vhf', format='str')
# opObj11.addParameter(name='ftp_wei', value='0', format='int')
# opObj11.addParameter(name='exp_code', value='19', format='int')
# opObj11.addParameter(name='sub_exp_code', value='0', format='int')
# opObj11.addParameter(name='plot_pos', value='0', format='int')
Daniel Valdez
r468
########################################################
############### ISR ####################################
########################################################
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467
procUnitConfObjISR = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId())
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 profileIndex = '20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99'
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')
opObj11.addParameter(name='m', value='5', format='int')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')
opObj11.addParameter(name='window', value='20', format='int')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 barker3x1 = '1,1,-1,-1,-1,1'
barker3x5 = '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,1,1'
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')
opObj11.addParameter(name='code', value=barker3x5, format='floatlist')
opObj11.addParameter(name='nCode', value='2', format='int')
opObj11.addParameter(name='nBaud', value='15', format='int')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjISR.getId())
procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')
procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')
Daniel Valdez
r472
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467 opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')
opObj11.addParameter(name='n', value='36', format='float')
Daniel Valdez
r472
Daniel Valdez
r468 # opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraWriter', optype='other')
# opObj11.addParameter(name='path', value='/media/datos/isr2014')
# opObj11.addParameter(name='blocksPerFile', value='120', format='int')
Daniel Valdez
r472
Daniel Valdez
r468
opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')
Daniel Valdez
r472 opObj11.addParameter(name='id', value='101', format='int')
Daniel Valdez
r468 opObj11.addParameter(name='wintitle', value='ISR', format='str')
opObj11.addParameter(name='showprofile', value='0', format='int')
Daniel Valdez
r472 opObj11.addParameter(name='channelList', value='3', format='intlist')
#opObj11.addParameter(name='timerange', value='300', format='float')
opObj11.addParameter(name='xmin', value='8', format='float')
opObj11.addParameter(name='xmax', value='16', format='float')
# opObj11.addParameter(name='ymin', value='120', format='float')
# opObj11.addParameter(name='ymax', value='190', format='float')
opObj11.addParameter(name='zmin', value='20', format='float')
opObj11.addParameter(name='zmax', value='55', format='float')
Daniel Valdez
r468 opObj11.addParameter(name='save', value='1', format='int')
Daniel Valdez
r472 opObj11.addParameter(name='figfile', value='rti_mst_isr.pdf', format='str')
opObj11.addParameter(name='figpath', value='/Users/dsuarez/Pictures/poster', format='str')
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467
Daniel Valdez
r468
Daniel Valdez
r472 # opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')
# opObj11.addParameter(name='id', value='2000', format='int')
# opObj11.addParameter(name='wintitle', value='ISR', format='str')
# opObj11.addParameter(name='showprofile', value='0', format='int')
# opObj11.addParameter(name='xmin', value='0', format='int')
# opObj11.addParameter(name='xmax', value='24', format='int')
# opObj11.addParameter(name='zmin', value='30', format='int')
# opObj11.addParameter(name='zmax', value='70', format='int')
#
# opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')
# opObj11.addParameter(name='id', value='2001', format='int')
# opObj11.addParameter(name='wintitle', value='ISR', format='str')
# opObj11.addParameter(name='zmin', value='20', format='int')
# opObj11.addParameter(name='zmax', value='60', format='int')
#
# opObj11.addParameter(name='save', value='1', format='int')
# opObj11.addParameter(name='figpath', value=gpath, format='str')
# opObj11.addParameter(name='wr_period', value='5', format='int')
# opObj11.addParameter(name='ftp', value='1', format='int')
# opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
# opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
# opObj11.addParameter(name='username', value='wmaster', format='str')
# opObj11.addParameter(name='password', value='mst2010vhf', format='str')
# opObj11.addParameter(name='ftp_wei', value='0', format='int')
# opObj11.addParameter(name='exp_code', value='20', format='int')
# opObj11.addParameter(name='sub_exp_code', value='0', format='int')
# opObj11.addParameter(name='plot_pos', value='0', format='int')
#
# #########################################################
# ################ EEJ ####################################
# #########################################################
# #########################################################
#
# procUnitConfObjEEJ = controllerObj.addProcUnit(datatype='Voltage', inputId=readUnitConfObj.getId())
#
# opObj11 = procUnitConfObjEEJ.addOperation(name='ProfileSelector', optype='other')
# opObj11.addParameter(name='profileRangeList', value='120,183', format='intlist')
#
# opObj11 = procUnitConfObjEEJ.addOperation(name='Decoder', optype='other')
# opObj11.addParameter(name='code', value='1,-1', format='floatlist')
# opObj11.addParameter(name='nCode', value='2', format='int')
# opObj11.addParameter(name='nBaud', value='1', format='int')
#
# procUnitConfObjEEJSpecta = controllerObj.addProcUnit(datatype='Spectra', inputId=procUnitConfObjEEJ.getId())
# procUnitConfObjEEJSpecta.addParameter(name='nFFTPoints', value='64', format='int')
# procUnitConfObjEEJSpecta.addParameter(name='nProfiles', value='64', format='int')
#
# opObj11 = procUnitConfObjEEJSpecta.addOperation(name='IncohInt', optype='other')
# opObj11.addParameter(name='n', value='36', format='float')
#
# # opObj11 = procUnitConfObjEEJSpecta.addOperation(name='SpectraWriter', optype='other')
# # opObj11.addParameter(name='path', value='/media/datos/eej2014')
# # opObj11.addParameter(name='blocksPerFile', value='10', format='int')
#
#
# opObj11 = procUnitConfObjEEJSpecta.addOperation(name='RTIPlot', optype='other')
# opObj11.addParameter(name='id', value='3000', format='int')
# opObj11.addParameter(name='wintitle', value='EEJ', format='str')
# opObj11.addParameter(name='showprofile', value='0', format='int')
# opObj11.addParameter(name='xmin', value='0', format='int')
# opObj11.addParameter(name='xmax', value='24', format='int')
# opObj11.addParameter(name='zmin', value='20', format='int')
# opObj11.addParameter(name='zmax', value='50', format='int')
#
# opObj11 = procUnitConfObjEEJSpecta.addOperation(name='SpectraPlot', optype='other')
# opObj11.addParameter(name='id', value='3001', format='int')
# opObj11.addParameter(name='wintitle', value='EEJ', format='str')
# opObj11.addParameter(name='zmin', value='20', format='int')
# opObj11.addParameter(name='zmax', value='50', format='int')
# opObj11.addParameter(name='save', value='1', format='int')
# opObj11.addParameter(name='figpath', value=gpath, format='str')
# opObj11.addParameter(name='wr_period', value='5', format='int')
# opObj11.addParameter(name='ftp', value='1', format='int')
# opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
# opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
# opObj11.addParameter(name='username', value='wmaster', format='str')
# opObj11.addParameter(name='password', value='mst2010vhf', format='str')
# opObj11.addParameter(name='ftp_wei', value='0', format='int')
# opObj11.addParameter(name='exp_code', value='22', format='int')
# opObj11.addParameter(name='sub_exp_code', value='0', format='int')
# opObj11.addParameter(name='plot_pos', value='0', format='int')
Daniel Valdez
Scripts para procesar y graficar datos de MST-ISR-EEJ
r467
print "Escribiendo el archivo XML"
controllerObj.writeXml(filename)
print "Leyendo el archivo XML"
controllerObj.readXml(filename)
controllerObj.createObjects()
controllerObj.connectObjects()
controllerObj.run()