##// END OF EJS Templates
mst-isr-eej processing script
Daniel Valdez -
r518:def3e48efb64
parent child
Show More
@@ -1,5 +1,5
1 import os, sys
1 import os, sys
2 import timeit
2 #import timeit
3
3
4 path = os.path.split(os.getcwd())[0]
4 path = os.path.split(os.getcwd())[0]
5 sys.path.append(path)
5 sys.path.append(path)
@@ -13,16 +13,14 controllerObj = Project()
13
13
14 controllerObj.setup(id = '191', name='test01', description=desc)
14 controllerObj.setup(id = '191', name='test01', description=desc)
15
15
16 path='/remote/ewdrifts/RAW_EXP/EW_DRIFT_FARADAY/EW_Drift'
16 path = '/Volumes/HD-PXU2/mst_isr_eej'
17
17
18 path = '/media/administrator/New Volume/DATA/MST_ISR'
18 figpath = '/Users/dsuarez/Pictures/mst_isr_eej/'
19
20 path = '/home/administrator/Documents/mst_isr_eej'
21
19
22 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
20 readUnitConfObj = controllerObj.addReadUnit(datatype='VoltageReader',
23 path=path,
21 path=path,
24 startDate='2014/01/09',
22 startDate='2014/05/01',
25 endDate='2014/01/09',
23 endDate='2014/05/30',
26 startTime='00:00:00',
24 startTime='00:00:00',
27 endTime='23:59:59',
25 endTime='23:59:59',
28 online=0,
26 online=0,
@@ -61,124 +59,124 opObj11.addParameter(name='id', value='100', format='int')
61 opObj11.addParameter(name='wintitle', value='EEJ', format='str')
59 opObj11.addParameter(name='wintitle', value='EEJ', format='str')
62 # opObj11.addParameter(name='zmin', value='20', format='int')
60 # opObj11.addParameter(name='zmin', value='20', format='int')
63 # opObj11.addParameter(name='zmax', value='60', format='int')# opObj11.addParameter(name='ftp', value='1', format='int')
61 # opObj11.addParameter(name='zmax', value='60', format='int')# opObj11.addParameter(name='ftp', value='1', format='int')
64 # opObj11.addParameter(name='save', value='1', format='int')
62 opObj11.addParameter(name='save', value='1', format='int')
65 # opObj11.addParameter(name='figpath', value='/home/operaciones/Pictures/MST-ISR', format='str')
63 opObj11.addParameter(name='figpath', value=figpath, format='str')
66 # opObj11.addParameter(name='wr_period', value='15', format='int')
64 opObj11.addParameter(name='wr_period', value='5', format='int')
67 # opObj11.addParameter(name='ftp', value='1', format='int')
68 # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
69 # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
70 # opObj11.addParameter(name='username', value='wmaster', format='str')
71 # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
72 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
65 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
73 # opObj11.addParameter(name='exp_code', value='22', format='int')
66 opObj11.addParameter(name='exp_code', value='22', format='int')
74 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
67 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
75 # opObj11.addParameter(name='plot_pos', value='0', format='int')
68 # opObj11.addParameter(name='plot_pos', value='0', format='int')
76
69
77
70
78 # opObj11 = procUnitConfObjEEJSpecta.addOperation(name='RTIPlot', optype='other')
71 opObj11 = procUnitConfObjEEJSpecta.addOperation(name='RTIPlot', optype='other')
79 # opObj11.addParameter(name='id', value='101', format='int')
72 opObj11.addParameter(name='id', value='101', format='int')
80 # opObj11.addParameter(name='wintitle', value='EEJ', format='str')
73 opObj11.addParameter(name='wintitle', value='EEJ', format='str')
81 # opObj11.addParameter(name='showprofile', value='0', format='int')
74 opObj11.addParameter(name='showprofile', value='0', format='int')
82 # opObj11.addParameter(name='xmin', value='0', format='int')
75 # opObj11.addParameter(name='xmin', value='0', format='int')
83 # opObj11.addParameter(name='xmax', value='24', format='int')
76 # opObj11.addParameter(name='xmax', value='24', format='int')
84 # opObj11.addParameter(name='save', value='1', format='int')
77 # opObj11.addParameter(name='save', value='1', format='int')
85 # opObj11.addParameter(name='figpath', value='/home/operaciones/Pictures/MST-ISR', format='str')
78 opObj11.addParameter(name='figpath', value=figpath, format='str')
86 # opObj11.addParameter(name='wr_period', value='15', format='int')
79 opObj11.addParameter(name='wr_period', value='5', format='int')
87 # opObj11.addParameter(name='ftp', value='1', format='int')
88 # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
89 # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
90 # opObj11.addParameter(name='username', value='wmaster', format='str')
91 # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
92 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
80 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
93 # opObj11.addParameter(name='exp_code', value='22', format='int')
81 opObj11.addParameter(name='exp_code', value='22', format='int')
94 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
82 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
95 # opObj11.addParameter(name='plot_pos', value='0', format='int')
83 # opObj11.addParameter(name='plot_pos', value='0', format='int')
96
84
85 # opObj11 = procUnitConfObjEEJSpecta.addOperation(name='SendByFTP', optype='other')
86 # opObj11.addParameter(name='ext', value='*.png', format='str')
87 # opObj11.addParameter(name='localfolder', value=figpath, format='str')
88 # opObj11.addParameter(name='remotefolder', value='/home/wmaster/graficos', format='str')
89 # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
90 # opObj11.addParameter(name='username', value='wmaster', format='str')
91 # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
92 # opObj11.addParameter(name='period', value='5', format='int')
93
94
97 ################ MST ####################################
95 ################ MST ####################################
98 procUnitConfObjMST = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
96 procUnitConfObjMST = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
99
97
100 opObj11 = procUnitConfObjMST.addOperation(name='ProfileSelector', optype='other')
98 opObj11 = procUnitConfObjMST.addOperation(name='ProfileSelector', optype='other')
101 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'
99 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'
102 #profileIndex = '0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19'
100 #profileIndex = '0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19'
103 opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
101 opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
104 opObj11.addParameter(name='byblock', value='1', format='bool')
102 opObj11.addParameter(name='byblock', value='1', format='bool')
105
103
106 opObj11 = procUnitConfObjMST.addOperation(name='Decoder', optype='other')
104 opObj11 = procUnitConfObjMST.addOperation(name='Decoder', optype='other')
107 opObj11.addParameter(name='mode',value='3',format='int')
105 opObj11.addParameter(name='mode',value='3',format='int')
108 opObj11.addParameter(name='times',value='10',format='int')
106 opObj11.addParameter(name='times',value='10',format='int')
109
107
110 opObj11 = procUnitConfObjMST.addOperation(name='CohInt', optype='other')
108 opObj11 = procUnitConfObjMST.addOperation(name='CohInt', optype='other')
111 opObj11.addParameter(name='n', value='20', format='int')
109 opObj11.addParameter(name='n', value='20', format='int')
112 opObj11.addParameter(name='byblock', value='1', format='bool')
110 opObj11.addParameter(name='byblock', value='1', format='bool')
113
111
114 procUnitConfObjMSTSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjMST.getId())
112 procUnitConfObjMSTSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjMST.getId())
115 procUnitConfObjMSTSpectra.addParameter(name='nFFTPoints', value='64', format='int')
113 procUnitConfObjMSTSpectra.addParameter(name='nFFTPoints', value='64', format='int')
116 procUnitConfObjMSTSpectra.addParameter(name='nProfiles', value='64', format='int')
114 procUnitConfObjMSTSpectra.addParameter(name='nProfiles', value='64', format='int')
117
115
118 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='IncohInt', optype='other')
116 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='IncohInt', optype='other')
119 opObj11.addParameter(name='n', value='2', format='float')
117 opObj11.addParameter(name='n', value='2', format='float')
120
118
121 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SpectraPlot', optype='other')
119 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='SpectraPlot', optype='other')
122 opObj11.addParameter(name='id', value='200', format='int')
120 opObj11.addParameter(name='id', value='200', format='int')
123 opObj11.addParameter(name='wintitle', value='MST', format='str')
121 opObj11.addParameter(name='wintitle', value='MST', format='str')
124 # opObj11.addParameter(name='zmin', value='35', format='int')
122 # # opObj11.addParameter(name='zmin', value='35', format='int')
125 # opObj11.addParameter(name='zmax', value='60', format='int')
123 # # opObj11.addParameter(name='zmax', value='60', format='int')
126 # opObj11.addParameter(name='save', value='1', format='int')
124 # # opObj11.addParameter(name='save', value='1', format='int')
127 # opObj11.addParameter(name='figpath', value='/home/operaciones/Pictures/MST-ISR', format='str')
125 opObj11.addParameter(name='figpath', value=figpath, format='str')
128 # opObj11.addParameter(name='wr_period', value='5', format='int')
126 opObj11.addParameter(name='wr_period', value='5', format='int')
129 # opObj11.addParameter(name='ftp', value='1', format='int')
127 # # opObj11.addParameter(name='ftp', value='1', format='int')
130 # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
128 # # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
131 # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
129 # # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
132 # opObj11.addParameter(name='username', value='wmaster', format='str')
130 # # opObj11.addParameter(name='username', value='wmaster', format='str')
133 # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
131 # # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
134 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
132 # # opObj11.addParameter(name='ftp_wei', value='0', format='int')
135 # opObj11.addParameter(name='exp_code', value='19', format='int')
133 opObj11.addParameter(name='exp_code', value='19', format='int')
136 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
134 # # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
137 # opObj11.addParameter(name='plot_pos', value='0', format='int')
135 # # opObj11.addParameter(name='plot_pos', value='0', format='int')
138 #
136 # #
139 # opObj11 = procUnitConfObjMSTSpectra.addOperation(name='RTIPlot', optype='other')
137 opObj11 = procUnitConfObjMSTSpectra.addOperation(name='RTIPlot', optype='other')
140 # opObj11.addParameter(name='id', value='201', format='int')
138 opObj11.addParameter(name='id', value='201', format='int')
141 # opObj11.addParameter(name='wintitle', value='MST', format='str')
139 opObj11.addParameter(name='wintitle', value='MST', format='str')
142 # opObj11.addParameter(name='showprofile', value='0', format='int')
140 opObj11.addParameter(name='showprofile', value='0', format='int')
143 # opObj11.addParameter(name='xmin', value='0', format='int')
141 # # opObj11.addParameter(name='xmin', value='0', format='int')
144 # opObj11.addParameter(name='xmax', value='24', format='int')
142 # # opObj11.addParameter(name='xmax', value='24', format='int')
145 # opObj11.addParameter(name='zmin', value='35', format='int')
143 # # opObj11.addParameter(name='zmin', value='35', format='int')
146 # opObj11.addParameter(name='zmax', value='60', format='int')
144 # # opObj11.addParameter(name='zmax', value='60', format='int')
147 # opObj11.addParameter(name='save', value='1', format='int')
145 # # opObj11.addParameter(name='save', value='1', format='int')
148 # opObj11.addParameter(name='figpath', value='/home/operaciones/Pictures/MST-ISR', format='str')
146 opObj11.addParameter(name='figpath', value=figpath, format='str')
149 # opObj11.addParameter(name='wr_period', value='5', format='int')
147 opObj11.addParameter(name='wr_period', value='5', format='int')
150 # opObj11.addParameter(name='ftp', value='1', format='int')
148 # # opObj11.addParameter(name='ftp', value='1', format='int')
151 # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
149 # # opObj11.addParameter(name='server', value='jro-app.igp.gob.pe', format='str')
152 # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
150 # # opObj11.addParameter(name='folder', value='/home/wmaster/graficos', format='str')
153 # opObj11.addParameter(name='username', value='wmaster', format='str')
151 # # opObj11.addParameter(name='username', value='wmaster', format='str')
154 # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
152 # # opObj11.addParameter(name='password', value='mst2010vhf', format='str')
155 # opObj11.addParameter(name='ftp_wei', value='0', format='int')
153 # # opObj11.addParameter(name='ftp_wei', value='0', format='int')
156 # opObj11.addParameter(name='exp_code', value='19', format='int')
154 opObj11.addParameter(name='exp_code', value='19', format='int')
157 # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
155 # # opObj11.addParameter(name='sub_exp_code', value='0', format='int')
158 # opObj11.addParameter(name='plot_pos', value='0', format='int')
156 # # opObj11.addParameter(name='plot_pos', value='0', format='int')
159
157
160 # ################ ISR ####################################
158 # ################ ISR ####################################
161 procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
159 procUnitConfObjISR = controllerObj.addProcUnit(datatype='VoltageProc', inputId=readUnitConfObj.getId())
162
160
163 opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')
161 opObj11 = procUnitConfObjISR.addOperation(name='ProfileSelector', optype='other')
164 # 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'
162 # 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'
165 # opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
163 # opObj11.addParameter(name='profileList', value=profileIndex, format='intlist')
166 opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')
164 opObj11.addParameter(name='profileRangeList', value='20,99', format='intlist')
167 opObj11.addParameter(name='byblock', value='1', format='bool')
165 opObj11.addParameter(name='byblock', value='1', format='bool')
168
166
169 # opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')
167 # opObj11 = procUnitConfObjISR.addOperation(name='ProfileConcat', optype='other')
170 # opObj11.addParameter(name='m', value='5', format='int')
168 # opObj11.addParameter(name='m', value='5', format='int')
171
169
172 opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles
170 opObj11 = procUnitConfObjISR.addOperation(name='Reshaper', optype='other') #Esta Operacion opera sobre bloques y reemplaza el ProfileConcat que opera sobre perfiles
173 opObj11.addParameter(name='shape', value='4,16,6750', format='intlist')
171 opObj11.addParameter(name='shape', value='4,16,6750', format='intlist') # shape = (nchannels, nprofiles, nhieghts)
174
172
175 opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')
173 opObj11 = procUnitConfObjISR.addOperation(name='filterByHeights')
176 opObj11.addParameter(name='window', value='20', format='int')
174 opObj11.addParameter(name='window', value='20', format='int')
177 opObj11.addParameter(name='axis', value='2', format='int')
175 opObj11.addParameter(name='axis', value='2', format='int')
178
176
179 barker3x1 = '1,1,-1,-1,-1,1'
177 barker3x1 = '1,1,-1,-1,-1,1'
180 #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'
178 #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'
181
179
182 opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')
180 opObj11 = procUnitConfObjISR.addOperation(name='Decoder', optype='other')
183 opObj11.addParameter(name='code', value=barker3x1, format='floatlist')
181 opObj11.addParameter(name='code', value=barker3x1, format='floatlist')
184 opObj11.addParameter(name='nCode', value='2', format='int')
182 opObj11.addParameter(name='nCode', value='2', format='int')
@@ -187,20 +185,30 opObj11.addParameter(name='nBaud', value='3', format='int')
187 opObj11.addParameter(name='mode', value='3', format='int')
185 opObj11.addParameter(name='mode', value='3', format='int')
188 opObj11.addParameter(name='times', value='8', format='int')
186 opObj11.addParameter(name='times', value='8', format='int')
189 opObj11.addParameter(name='osamp', value='5', format='int')
187 opObj11.addParameter(name='osamp', value='5', format='int')
190
188
191
189
192 procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())
190 procUnitConfObjISRSpectra = controllerObj.addProcUnit(datatype='SpectraProc', inputId=procUnitConfObjISR.getId())
193 procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')
191 procUnitConfObjISRSpectra.addParameter(name='nFFTPoints', value='16', format='int')
194 procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')
192 procUnitConfObjISRSpectra.addParameter(name='nProfiles', value='16', format='int')
195
193
196 opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')
194 opObj11 = procUnitConfObjISRSpectra.addOperation(name='IncohInt', optype='other')
197 opObj11.addParameter(name='n', value='36', format='float')
195 opObj11.addParameter(name='n', value='36', format='float')
198
196
199 opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')
197 opObj11 = procUnitConfObjISRSpectra.addOperation(name='SpectraPlot', optype='other')
200 opObj11.addParameter(name='id', value='300', format='int')
198 opObj11.addParameter(name='id', value='300', format='int')
201 opObj11.addParameter(name='wintitle', value='ISR', format='str')
199 opObj11.addParameter(name='wintitle', value='ISR', format='str')
202 # opObj11.addParameter(name='save', value='1', format='bool')
200 opObj11.addParameter(name='figpath', value=figpath, format='str')
203 # opObj11.addParameter(name='figpath', value='/home/administrator/Pictures/mst_isr_eej', format='str')
201 opObj11.addParameter(name='wr_period', value='5', format='int')
202 opObj11.addParameter(name='exp_code', value='20', format='int')
203
204 opObj11 = procUnitConfObjISRSpectra.addOperation(name='RTIPlot', optype='other')
205 opObj11.addParameter(name='id', value='301', format='int')
206 opObj11.addParameter(name='wintitle', value='ISR', format='str')
207 opObj11.addParameter(name='showprofile', value='0', format='int')
208 opObj11.addParameter(name='figpath', value=figpath, format='str')
209 opObj11.addParameter(name='wr_period', value='5', format='int')
210 opObj11.addParameter(name='exp_code', value='20', format='int')
211
204
212
205
213
206 print "Escribiendo el archivo XML"
214 print "Escribiendo el archivo XML"
General Comments 0
You need to be logged in to leave comments. Login now