##// END OF EJS Templates
AMISR Reader integration with Signal Chain Blocks, this time just only for Voltages to Profile Selection and Plotting Scope(Power,IQ) and Power Profile(dB). There is thwo python scripts as experiment's test.
AMISR Reader integration with Signal Chain Blocks, this time just only for Voltages to Profile Selection and Plotting Scope(Power,IQ) and Power Profile(dB). There is thwo python scripts as experiment's test.

File last commit:

r472:bb102d4de45c
r474:e3c42f4268a9
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()