##// END OF EJS Templates
schain
RSS
scripts to processing rawdata of MST-ISR-EEJ experiment
Daniel Valdez -
r519:e595cfcb98a5
parent child
Show More
Show file before | Show file after | Hide comments
@@ -0,0 +1,111
1 import os, sys
2 #import timeit
3 import datetime
4
5 path = os.path.split(os.getcwd())[0]
6 sys.path.append(path)
7
8 from controller import *
9 dt1 = datetime.datetime.now()
10 desc = "MST-ISR-EEJ Experiment Test"
11 filename = "eej_blocks.xml"
12
13 controllerObj = Project()
14
15 controllerObj.setup(id = '191', name='test01', description=desc)
16
17 path = '/media/signalchain/HD-PXU2/mst_isr_eej'
18
19 figpath = '/home/signalchain/Pictures/mst_isr_eej/eej'
20
21 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
22 path=path,
23 startDate='2014/05/01',
24 endDate='2014/05/30',
25 startTime='00:00:00',
26 endTime='23:59:59',
27 online=0,
28 delay=10,
29 walk=0,
30 getblock=1)
31
32 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
33 # ################ EEJ ####################################
34 procUnitConfObjEEJ = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
35
36 opObj11 = procUnitConfObjEEJ.addOperation(name='ProfileSelector', optype='other')
37 opObj11.addParameter(name='profileRangeList', value='120,183', format='intlist')
38 opObj11.addParameter(name='byblock', value='1', format='bool')
39
40 opObj11 = procUnitConfObjEEJ.addOperation(name='Decoder', optype='other')
41 opObj11.addParameter(name='code', value='1,-1', format='floatlist')
42 opObj11.addParameter(name='nCode', value='2', format='int')
43 opObj11.addParameter(name='nBaud', value='1', format='int')
44 opObj11.addParameter(name='mode', value='3', format='int')
45 opObj11.addParameter(name='times', value='32', format='int')
46
47 # opObj11 = procUnitConfObjEEJ.addOperation(name='CohInt', optype='other')
48 # opObj11.addParameter(name='n', value='2', format='int')
49
50 procUnitConfObjEEJSpecta = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjEEJ.getId())
51 procUnitConfObjEEJSpecta.addParameter(name='nFFTPoints', value='64', format='int')
52 procUnitConfObjEEJSpecta.addParameter(name='nProfiles', value='64', format='int')
53
54 opObj11 = procUnitConfObjEEJSpecta.addOperation(name='IncohInt', optype='other')
55 #opObj11.addParameter(name='timeInterval', value='10', format='float')
56 opObj11.addParameter(name='n', value='36', format='float')
57
58 opObj11 = procUnitConfObjEEJSpecta.addOperation(name='SpectraPlot', optype='other')
59 opObj11.addParameter(name='id', value='100', format='int')
60 opObj11.addParameter(name='wintitle', value='EEJ', format='str')
61 # opObj11.addParameter(name='zmin', value='20', format='int')
62 # opObj11.addParameter(name='zmax', value='60', format='int')# opObj11.addParameter(name='ftp', value='1', format='int')
63 opObj11.addParameter(name='save', value='1', format='int')
64 opObj11.addParameter(name='figpath', value=figpath, format='str')
65 opObj11.addParameter(name='wr_period', value='5', format='int')
66 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
67 opObj11.addParameter(name='exp_code', value='22', format='int')
68 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
69 # opObj11.addParameter(name='plot_pos', value='0', format='int')
70
71
72 opObj11 = procUnitConfObjEEJSpecta.addOperation(name='RTIPlot', optype='other')
73 opObj11.addParameter(name='id', value='101', format='int')
74 opObj11.addParameter(name='wintitle', value='EEJ', format='str')
75 opObj11.addParameter(name='showprofile', value='0', format='int')
76 # opObj11.addParameter(name='xmin', value='0', format='int')
77 # opObj11.addParameter(name='xmax', value='24', format='int')
78 # opObj11.addParameter(name='save', value='1', format='int')
79 opObj11.addParameter(name='figpath', value=figpath, format='str')
80 opObj11.addParameter(name='wr_period', value='5', format='int')
81 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
82 opObj11.addParameter(name='exp_code', value='22', format='int')
83 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
84 # opObj11.addParameter(name='plot_pos', value='0', format='int')
85
86 opObj11 = procUnitConfObjEEJSpecta.addOperation(name='SendByFTP', optype='other')
87 opObj11.addParameter(name='ext', value='*.png', format='str')
88 opObj11.addParameter(name='localfolder', value=figpath, format='str')
89 opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')
90 opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
91 opObj11.addParameter(name='username', value='wmaster', format='str')
92 opObj11.addParameter(name='password', value='mst2010vhf', format='str')
93 opObj11.addParameter(name='period', value='5', format='int')
94
95 print "Escribiendo el archivo XML"
96 controllerObj.writeXml(filename)
97 print "Leyendo el archivo XML"
98 controllerObj.readXml(filename)
99
100 controllerObj.createObjects()
101 controllerObj.connectObjects()
102
103 #timeit.timeit('controllerObj.run()', number=2)
104
105 controllerObj.run()
106 #print fib(5)
107
108 dt2 = datetime.datetime.now()
109 print "======================="
110 print dt2-dt1
111 print "=======================" No newline at end of file
Show file before | Show file after | Hide comments
@@ -0,0 +1,115
1 import os, sys
2 #import timeit
3 import datetime
4
5 path = os.path.split(os.getcwd())[0]
6 sys.path.append(path)
7
8 from controller import *
9 dt1 = datetime.datetime.now()
10 desc = "MST-ISR-EEJ Experiment Test"
11 filename = "isr_blocks.xml"
12
13 controllerObj = Project()
14
15 controllerObj.setup(id = '191', name='test01', description=desc)
16
17 path = '/media/signalchain/HD-PXU2/mst_isr_eej'
18
19 figpath = '/home/signalchain/Pictures/mst_isr_eej/isr'
20
21 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
22 path=path,
23 startDate='2014/05/01',
24 endDate='2014/05/30',
25 startTime='00:00:00',
26 endTime='23:59:59',
27 online=0,
28 delay=10,
29 walk=0,
30 getblock=1)
31
32 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
33
34 procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
35
36 opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')
37 # 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'
38 # opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
39 opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')
40 opObj11.addParameter(name='byblock', value='1', format='bool')
41
42 # opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')
43 # opObj11.addParameter(name='m', value='5', format='int')
44
45 opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles
46 opObj11.addParameter(name='shape', value='4,16,6750', format='intlist') # shape = (nchannels, nprofiles, nhieghts)
47
48 opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')
49 opObj11.addParameter(name='window', value='20', format='int')
50 opObj11.addParameter(name='axis', value='2', format='int')
51
52 barker3x1 = '1,1,-1,-1,-1,1'
53 #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'
54
55 opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')
56 opObj11.addParameter(name='code', value=barker3x1, format='floatlist')
57 opObj11.addParameter(name='nCode', value='2', format='int')
58 #opObj11.addParameter(name='nBaud', value='15', format='int')
59 opObj11.addParameter(name='nBaud', value='3', format='int')
60 opObj11.addParameter(name='mode', value='3', format='int')
61 opObj11.addParameter(name='times', value='8', format='int')
62 opObj11.addParameter(name='osamp', value='5', format='int')
63
64
65 procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())
66 procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')
67 procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')
68
69 opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')
70 opObj11.addParameter(name='n', value='36', format='float')
71
72 opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')
73 opObj11.addParameter(name='id', value='300', format='int')
74 opObj11.addParameter(name='wintitle', value='ISR', format='str')
75 opObj11.addParameter(name='figpath', value=figpath, format='str')
76 opObj11.addParameter(name='wr_period', value='5', format='int')
77 opObj11.addParameter(name='exp_code', value='20', format='int')
78
79 opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')
80 opObj11.addParameter(name='id', value='301', format='int')
81 opObj11.addParameter(name='wintitle', value='ISR', format='str')
82 opObj11.addParameter(name='showprofile', value='0', format='int')
83 opObj11.addParameter(name='figpath', value=figpath, format='str')
84 opObj11.addParameter(name='wr_period', value='5', format='int')
85 opObj11.addParameter(name='exp_code', value='20', format='int')
86
87
88 opObj11 = procUnitConfObjISRSpectra.addOperation(name='SendByFTP', optype='other')
89 opObj11.addParameter(name='ext', value='*.png', format='str')
90 opObj11.addParameter(name='localfolder', value=figpath, format='str')
91 opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')
92 opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
93 opObj11.addParameter(name='username', value='wmaster', format='str')
94 opObj11.addParameter(name='password', value='mst2010vhf', format='str')
95 opObj11.addParameter(name='period', value='5', format='int')
96
97
98
99 print "Escribiendo el archivo XML"
100 controllerObj.writeXml(filename)
101 print "Leyendo el archivo XML"
102 controllerObj.readXml(filename)
103
104 controllerObj.createObjects()
105 controllerObj.connectObjects()
106
107 #timeit.timeit('controllerObj.run()', number=2)
108
109 controllerObj.run()
110 #print fib(5)
111
112 dt2 = datetime.datetime.now()
113 print "======================="
114 print dt2-dt1
115 print "======================="
Show file before | Show file after | Hide comments
@@ -0,0 +1,121
1 import os, sys
2 #import timeit
3 import datetime
4
5 path = os.path.split(os.getcwd())[0]
6 sys.path.append(path)
7
8 from controller import *
9 dt1 = datetime.datetime.now()
10 desc = "MST-ISR-EEJ Experiment Test"
11 filename = "mst_blocks.xml"
12
13 controllerObj = Project()
14
15 controllerObj.setup(id = '191', name='test01', description=desc)
16
17 path = '/media/signalchain/HD-PXU2/mst_isr_eej'
18
19 figpath = '/home/signalchain/Pictures/mst_isr_eej/mst'
20
21 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
22 path=path,
23 startDate='2014/05/01',
24 endDate='2014/05/30',
25 startTime='00:00:00',
26 endTime='23:59:59',
27 online=0,
28 delay=10,
29 walk=0,
30 getblock=1)
31
32 opObj11 = readUnitConfObj.addOperation(name='printNumberOfBlock')
33
34 procUnitConfObjMST = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
35
36 opObj11 = procUnitConfObjMST.addOperation(name='ProfileSelector', optype='other')
37 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'
38 #profileIndex = '0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19'
39 opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
40 opObj11.addParameter(name='byblock', value='1', format='bool')
41
42 opObj11 = procUnitConfObjMST.addOperation(name='Decoder', optype='other')
43 opObj11.addParameter(name='mode',value='3',format='int')
44 opObj11.addParameter(name='times',value='10',format='int')
45
46 opObj11 = procUnitConfObjMST.addOperation(name='CohInt', optype='other')
47 opObj11.addParameter(name='n', value='20', format='int')
48 opObj11.addParameter(name='byblock', value='1', format='bool')
49
50 procUnitConfObjMSTSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjMST.getId())
51 procUnitConfObjMSTSpectra.addParameter(name='nFFTPoints', value='64', format='int')
52 procUnitConfObjMSTSpectra.addParameter(name='nProfiles', value='64', format='int')
53
54 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='IncohInt', optype='other')
55 opObj11.addParameter(name='n', value='2', format='float')
56
57 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SpectraPlot', optype='other')
58 opObj11.addParameter(name='id', value='200', format='int')
59 opObj11.addParameter(name='wintitle', value='MST', format='str')
60 # # opObj11.addParameter(name='zmin', value='35', format='int')
61 # # opObj11.addParameter(name='zmax', value='60', format='int')
62 # # opObj11.addParameter(name='save', value='1', format='int')
63 opObj11.addParameter(name='figpath', value=figpath, format='str')
64 opObj11.addParameter(name='wr_period', value='5', format='int')
65 # # opObj11.addParameter(name='ftp', value='1', format='int')
66 # # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
67 # # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
68 # # opObj11.addParameter(name='username', value='wmaster', format='str')
69 # # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
70 # # opObj11.addParameter(name='ftp_wei', value='0', format='int')
71 opObj11.addParameter(name='exp_code', value='19', format='int')
72 # # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
73 # # opObj11.addParameter(name='plot_pos', value='0', format='int')
74 # #
75 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='RTIPlot', optype='other')
76 opObj11.addParameter(name='id', value='201', format='int')
77 opObj11.addParameter(name='wintitle', value='MST', format='str')
78 opObj11.addParameter(name='showprofile', value='0', format='int')
79 # # opObj11.addParameter(name='xmin', value='0', format='int')
80 # # opObj11.addParameter(name='xmax', value='24', format='int')
81 # # opObj11.addParameter(name='zmin', value='35', format='int')
82 # # opObj11.addParameter(name='zmax', value='60', format='int')
83 # # opObj11.addParameter(name='save', value='1', format='int')
84 opObj11.addParameter(name='figpath', value=figpath, format='str')
85 opObj11.addParameter(name='wr_period', value='5', format='int')
86 # # opObj11.addParameter(name='ftp', value='1', format='int')
87 # # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
88 # # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
89 # # opObj11.addParameter(name='username', value='wmaster', format='str')
90 # # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
91 # # opObj11.addParameter(name='ftp_wei', value='0', format='int')
92 opObj11.addParameter(name='exp_code', value='19', format='int')
93 # # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
94 # # opObj11.addParameter(name='plot_pos', value='0', format='int')
95
96 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SendByFTP', optype='other')
97 opObj11.addParameter(name='ext', value='*.png', format='str')
98 opObj11.addParameter(name='localfolder', value=figpath, format='str')
99 opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')
100 opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
101 opObj11.addParameter(name='username', value='wmaster', format='str')
102 opObj11.addParameter(name='password', value='mst2010vhf', format='str')
103 opObj11.addParameter(name='period', value='5', format='int')
104
105 print "Escribiendo el archivo XML"
106 controllerObj.writeXml(filename)
107 print "Leyendo el archivo XML"
108 controllerObj.readXml(filename)
109
110 controllerObj.createObjects()
111 controllerObj.connectObjects()
112
113 #timeit.timeit('controllerObj.run()', number=2)
114
115 controllerObj.run()
116 #print fib(5)
117
118 dt2 = datetime.datetime.now()
119 print "======================="
120 print dt2-dt1
121 print "=======================" No newline at end of file
General Comments 0
You need to be logged in to leave comments. Login now