##// END OF EJS Templates
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014
This is the new organization by packages and scripts for Signal Chain, this version contains new features and bugs fixed until August 2014

File last commit:

r472:bb102d4de45c
r487:89975db10cce
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()