##// END OF EJS Templates
schain
RSS
Update sophy script
jespinoza -
r1716:d0250cd06eec
parent child
Show More
Side by Side Unified
Context file:
r1716:d0250cd06eec Loading diff...:
Show file before | Show file after | Hide comments
@@ -1,478 +1,489
1 1 # SOPHY PROC script
2 2 import os, sys, json, argparse
3 3 import multiprocessing
4 4 import datetime
5 5 import time
6 6
7 7 PATH = '/data'
8 8 #PATH = "/media/soporte/TOSHIBAEXT/sophy/"
9
10
11
9 12 # SNR ZMIN -40 A ZMAX -20
10 13 PARAM = {
11 14 'S': {'zmin': -45, 'zmax':-25, 'colormap': 'jet' , 'label': 'Power', 'wrname': 'power','cb_label': 'dBm', 'ch':0},
12 15 'SNR':{'zmin': -10, 'zmax': 15, 'colormap': 'jet' , 'label': 'SNR', 'wrname': 'snr','cb_label': 'dB', 'ch':0},
13 16 'V': {'zmin': -12, 'zmax': 12, 'colormap': 'sophy_v', 'label': 'Velocity', 'wrname': 'velocity', 'cb_label': 'm/s', 'ch':0},
14 17 'R': {'zmin': 0.5, 'zmax': 1 , 'colormap': 'sophy_r', 'label': 'RhoHV', 'wrname':'rhoHV', 'cb_label': '', 'ch':0},
15 18 'P': {'zmin': -180,'zmax': 180,'colormap': 'sophy_p', 'label': 'PhiDP', 'wrname':'phiDP' , 'cb_label': 'degrees', 'ch':0},
16 19 'D': {'zmin': -9 , 'zmax': 12, 'colormap': 'sophy_d', 'label': 'ZDR','wrname':'differential_reflectivity' , 'cb_label': 'dB','ch':0},
17 20 'Z': {'zmin': -20, 'zmax': 80, 'colormap': 'sophy_z', 'label': 'Reflectivity ', 'wrname':'reflectivity', 'cb_label': 'dBz','ch':0},
18 21 'W': {'zmin': 0 , 'zmax': 12, 'colormap': 'sophy_w', 'label': 'Spectral Width', 'wrname':'spectral_width', 'cb_label': 'm/s', 'ch':0}
19 22 }
20 23
21 24 META = ['heightList', 'data_azi', 'data_ele', 'mode_op', 'latitude', 'longitude', 'altitude', 'heading', 'radar_name',
22 25 'institution', 'contact', 'h0', 'range_unit', 'prf', 'prf_unit', 'variable', 'variable_unit', 'n_pulses',
23 26 'pulse1_range', 'pulse1_width', 'pulse2_width', 'pulse1_repetitions', 'pulse2_repetitions', 'pulse_width_unit',
24 27 'snr_threshold', 'data_noise']
25 28
26 29
27 30 def max_index(r, sample_rate, ipp, h0):
28 31
29 32 return int(sample_rate*ipp*1e6 * r / 60) + int(sample_rate*ipp*1e6 * -h0 / 60)
30 33
31 34 def main(args):
32 35
33 36 experiment = args.experiment
34 37 fp = open(os.path.join(PATH, experiment, 'experiment.json'))
35 38 conf = json.loads(fp.read())
36 39
37 40 ipp_km = conf['usrp_tx']['ipp']
38 41 ipp = ipp_km * 2 /300000
39 42 sample_rate = conf['usrp_rx']['sample_rate']
40 43 speed_axis = conf['pedestal']['speed']
41 44 if args.angles:
42 45 angles = args.angles
43 46 else:
44 47 angles = conf['pedestal']['table']
45 48 time_offset = args.time_offset
46 49 parameters = args.parameters
47 50 start_date = conf['name'].split('@')[1].split('T')[0].replace('-', '/')
48 51 end_date = start_date
49 52 if args.start_time:
50 53 start_time = args.start_time
51 54 else:
52 55 start_time = conf['name'].split('@')[1].split('T')[1].replace('-', ':')
53 56
54 57 if args.end_time:
55 58 end_time = args.end_time
56 59 else:
57 60 end_time = '23:59:59'
58 61
59 62 N = int(1.0/(abs(speed_axis[0])*ipp)) # 1 GRADO DE RESOLUCION
60 63
61 64 path = os.path.join(PATH, experiment, 'rawdata')
62 65 path_ped = os.path.join(PATH, experiment, 'position')
63 66 if args.label:
64 67 label = '-{}'.format(args.label)
65 68 else:
66 69 label = ''
67 70 path_plots = os.path.join(PATH, experiment, 'plots{}'.format(label))
68 71 path_save = os.path.join(PATH, experiment, 'param{}'.format(label))
69 72 RMIX = 4.8#5.8 #4.8#5.68#4.8#4.8#2.64#10#2.64
70 73 H0 = -1.68 #-1.68# -1.2#-1.68#-1.2#0.5#-1.2
71 74 MASK = 0.6
72 75 #MASK = 0.4#0.35
73 76
74 77 from schainpy.controller import Project
75 78
76 79 project = Project()
77 80 project.setup(id='1', name='Sophy', description='sophy proc')
78 81
79 82 reader = project.addReadUnit(datatype='DigitalRFReader',
80 83 path=path,
81 84 startDate=start_date,
82 85 endDate=end_date,
83 86 startTime=start_time,
84 87 endTime=end_time,
85 88 delay=30,
86 89 online=args.online,
87 90 walk=1,
88 91 ippKm = ipp_km,
89 92 getByBlock = 1,
90 93 nProfileBlocks = N,
91 94 )
92 95
93 96 if not conf['usrp_tx']['enable_2']: # One Pulse
94 97 n_pulses = 1
95 98 pulse_1_width = conf['usrp_tx']['pulse_1']
96 99 pulse_1_repetitions = conf['usrp_tx']['repetitions_1']
97 100 pulse_2_width = 0
98 101 pulse_2_repetitions = 0
99 102
100 103 voltage = project.addProcUnit(datatype='VoltageProc', inputId=reader.getId())
101 104
102 105 if conf['usrp_tx']['code_type_1'] != 'None':
103 106 codes = [ c.strip() for c in conf['usrp_tx']['code_1'].split(',')]
104 107 code = []
105 108 for c in codes:
106 109 code.append([int(x) for x in c])
107 110 op = voltage.addOperation(name='Decoder', optype='other')
108 111 op.addParameter(name='code', value=code)
109 112 op.addParameter(name='nCode', value=len(code), format='int')
110 113 op.addParameter(name='nBaud', value=len(code[0]), format='int')
111 114
112 115 op = voltage.addOperation(name='CohInt', optype='other') #Minimo integrar 2 perfiles por ser codigo complementario
113 116 op.addParameter(name='n', value=len(code), format='int')
114 117 ncode = len(code)
115 118 else:
116 119 ncode = 1
117 120 code = ['0']
118 121
119 122 op = voltage.addOperation(name='setH0')
120 123 op.addParameter(name='h0', value=H0)
121 124
122 125 if args.range > 0:
123 126 op = voltage.addOperation(name='selectHeights')
124 127 op.addParameter(name='minIndex', value='0', format='int')
125 128 op.addParameter(name='maxIndex', value=max_index(args.range, sample_rate, ipp, H0), format='int')
126 129
127 130 op = voltage.addOperation(name='PulsePair_vRF', optype='other')
128 131 op.addParameter(name='n', value=int(N)/ncode, format='int')
129 132 if args.rmDC:
130 133 op.addParameter(name='removeDC', value=1, format='int')
131 134
132 135 proc = project.addProcUnit(datatype='ParametersProc', inputId=voltage.getId())
133 136
134 137 opObj10 = proc.addOperation(name="WeatherRadar")
135 138 opObj10.addParameter(name='tauW',value=(1e-6/sample_rate)*len(code[0]))
136 139 # opObj10.addParameter(name='Pt',value=((1e-6/sample_rate)*len(code[0])/ipp)*200)
137 140 opObj10.addParameter(name='Pt',value=200)
138 141
139 142 op = proc.addOperation(name='PedestalInformation')
140 143 op.addParameter(name='path', value=path_ped, format='str')
141 144 op.addParameter(name='interval', value='0.04')
142 145 op.addParameter(name='time_offset', value=time_offset)
143 146 op.addParameter(name='mode', value=args.mode)
144 147
145 148 for param in parameters:
146 149 op = proc.addOperation(name='Block360')
147 150 op.addParameter(name='runNextOp', value=True)
148 151 op.addParameter(name='attr_data', value='data_param')
149 152 op.addParameter(name='angles', value=angles)
150 153
151 154 op= proc.addOperation(name='WeatherParamsPlot')
152 155 if args.save: op.addParameter(name='save', value=path_plots, format='str')
153 156 op.addParameter(name='save_period', value=-1)
154 157 op.addParameter(name='show', value=args.show)
155 158 op.addParameter(name='channels', value='0,')
156 159 op.addParameter(name='zmin', value=PARAM[param]['zmin'])
157 160 op.addParameter(name='zmax', value=PARAM[param]['zmax'])
161 op.addParameter(name='yrange', value=20, format='int')
162 op.addParameter(name='xrange', value=args.range, format='int')
158 163 op.addParameter(name='attr_data', value=param, format='str')
159 164 op.addParameter(name='labels', value=[PARAM[param]['label'], PARAM[param]['label']])
160 165 op.addParameter(name='save_code', value=param)
161 166 op.addParameter(name='cb_label', value=PARAM[param]['cb_label'])
162 167 op.addParameter(name='colormap', value=PARAM[param]['colormap'])
163 168 op.addParameter(name='bgcolor', value='black')
164 169 op.addParameter(name='localtime', value=False)
165 170 op.addParameter(name='shapes', value='./shapes')
171 op.addParameter(name='latitude', value=conf['latitude'], format='float')
172 op.addParameter(name='longitude', value=conf['longitude'], format='float')
173 op.addParameter(name='map', value=True)
174
166 175 if MASK: op.addParameter(name='mask', value=MASK, format='float')
167 176 if args.server:
168 177 op.addParameter(name='server', value='190.187.237.239:4444')
169 178 op.addParameter(name='exp_code', value='400')
170 179
171 180 desc = {
172 181 'Data': {
173 182 'data_param': {PARAM[param]['wrname']: ['H', 'V']},
174 183 'utctime': 'time'
175 184 },
176 185 'Metadata': {
177 186 'heightList': 'range',
178 187 'data_azi': 'azimuth',
179 188 'data_ele': 'elevation',
180 189 'mode_op': 'scan_type',
181 190 'h0': 'range_correction',
182 191 'dataPP_NOISE': 'noise',
183 192 }
184 193 }
185 194
186 195 if args.save:
187 196 writer = proc.addOperation(name='HDFWriter')
188 197 writer.addParameter(name='path', value=path_save, format='str')
189 198 writer.addParameter(name='Reset', value=True)
190 199 writer.addParameter(name='setType', value='weather')
191 200 writer.addParameter(name='description', value=json.dumps(desc))
192 201 writer.addParameter(name='blocksPerFile', value='1',format='int')
193 202 writer.addParameter(name='metadataList', value=','.join(META))
194 203 writer.addParameter(name='dataList', value='data_param,utctime')
195 204 writer.addParameter(name='weather_var', value=param)
196 205 writer.addParameter(name='mask', value=MASK, format='float')
197 206 writer.addParameter(name='localtime', value=False)
198 207 # meta
199 208 writer.addParameter(name='latitude', value='-12.040436')
200 209 writer.addParameter(name='longitude', value='-75.295893')
201 210 writer.addParameter(name='altitude', value='3379.2147')
202 211 writer.addParameter(name='heading', value='0')
203 212 writer.addParameter(name='radar_name', value='SOPHy')
204 213 writer.addParameter(name='institution', value='IGP')
205 214 writer.addParameter(name='contact', value='dscipion@igp.gob.pe')
206 215 writer.addParameter(name='created_by', value='Signal Chain (https://pypi.org/project/schainpy/)')
207 216 writer.addParameter(name='range_unit', value='km')
208 217 writer.addParameter(name='prf', value=1/ipp)
209 218 writer.addParameter(name='prf_unit', value='hertz')
210 219 writer.addParameter(name='variable', value=PARAM[param]['label'])
211 220 writer.addParameter(name='variable_unit', value=PARAM[param]['cb_label'])
212 221 writer.addParameter(name='n_pulses', value=n_pulses)
213 222 writer.addParameter(name='pulse1_range', value=RMIX)
214 223 writer.addParameter(name='pulse1_width', value=pulse_1_width)
215 224 writer.addParameter(name='pulse2_width', value=pulse_2_width)
216 225 writer.addParameter(name='pulse1_repetitions', value=pulse_1_repetitions)
217 226 writer.addParameter(name='pulse2_repetitions', value=pulse_2_repetitions)
218 227 writer.addParameter(name='pulse_width_unit', value='microseconds')
219 228 writer.addParameter(name='snr_threshold', value=MASK)
220 229
221 230
222 231 else: #Two pulses
223 232 n_pulses = 1
224 233 pulse_1_width = conf['usrp_tx']['pulse_1']
225 234 pulse_1_repetitions = conf['usrp_tx']['repetitions_1']
226 235 pulse_2_width = conf['usrp_tx']['pulse_2']
227 236 pulse_2_repetitions = conf['usrp_tx']['repetitions_2']
228 237
229 238 voltage1 = project.addProcUnit(datatype='VoltageProc', inputId=reader.getId())
230 239
231 240 op = voltage1.addOperation(name='ProfileSelector')
232 241 op.addParameter(name='profileRangeList', value='0,{}'.format(conf['usrp_tx']['repetitions_1']-1))
233 242
234 243 if conf['usrp_tx']['code_type_1'] != 'None':
235 244 codes = [ c.strip() for c in conf['usrp_tx']['code_1'].split(',')]
236 245 code = []
237 246 for c in codes:
238 247 code.append([int(x) for x in c])
239 248 op = voltage1.addOperation(name='Decoder', optype='other')
240 249 op.addParameter(name='code', value=code)
241 250 op.addParameter(name='nCode', value=len(code), format='int')
242 251 op.addParameter(name='nBaud', value=len(code[0]), format='int')
243 252 ncode = len(code)
244 253 else:
245 254 ncode = 1
246 255 code = ['0']
247 256
248 257 op = voltage1.addOperation(name='CohInt', optype='other') #Minimo integrar 2 perfiles por ser codigo complementario
249 258 op.addParameter(name='n', value=2, format='int')
250 259
251 260 op = voltage1.addOperation(name='setH0')
252 261 op.addParameter(name='h0', value=H0, format='float')
253 262
254 263 if args.range > 0:
255 264 op = voltage1.addOperation(name='selectHeights')
256 265 op.addParameter(name='minIndex', value=max_index(0, sample_rate, ipp, H0), format='int')
257 266 op.addParameter(name='maxIndex', value=max_index(RMIX, sample_rate, ipp, H0), format='int')
258 267
259 268 #op = voltage1.addOperation(name='setH0')
260 269 #op.addParameter(name='h0', value=H0, format='float')
261 270
262 271 op = voltage1.addOperation(name='PulsePair_vRF', optype='other')
263 272 op.addParameter(name='n', value=int(conf['usrp_tx']['repetitions_1'])/ncode, format='int')
264 273 if args.rmDC:
265 274 op.addParameter(name='removeDC', value=1, format='int')
266 275
267 276 proc1 = project.addProcUnit(datatype='ParametersProc', inputId=voltage1.getId())
268 277 proc1.addParameter(name='runNextUnit', value=True)
269 278
270 279 opObj10 = proc1.addOperation(name="WeatherRadar")
271 280 opObj10.addParameter(name='CR_Flag',value=True)
272 281 opObj10.addParameter(name='tauW',value=(1e-6/sample_rate)*len(code[0]))
273 282 #opObj10.addParameter(name='Pt',value=((1e-6/sample_rate)*len(code[0])/ipp)*200)
274 283 opObj10.addParameter(name='Pt',value=200)
275 284 #opObj10.addParameter(name='min_index',value=0)
276 285 opObj10.addParameter(name='min_index',value=max_index(0, sample_rate, ipp, H0))
277 286 #opObj10.addParameter(name='sesgoZD',value=7.73)
278 287
279 288
280 289 op = proc1.addOperation(name='PedestalInformation')
281 290 op.addParameter(name='path', value=path_ped, format='str')
282 291 op.addParameter(name='interval', value='0.04')
283 292 op.addParameter(name='time_offset', value=time_offset)
284 293 op.addParameter(name='mode', value=args.mode)
285 294
286 295 op = proc1.addOperation(name='Block360')
287 296 op.addParameter(name='attr_data', value='data_param')
288 297 op.addParameter(name='runNextOp', value=True)
289 298 op.addParameter(name='angles', value=angles)
290 299 #op.addParameter(name='horario',value=False)
291 300
292 301 voltage2 = project.addProcUnit(datatype='VoltageProc', inputId=reader.getId())
293 302
294 303 op = voltage2.addOperation(name='ProfileSelector')
295 304 op.addParameter(name='profileRangeList', value='{},{}'.format(conf['usrp_tx']['repetitions_1'], conf['usrp_tx']['repetitions_1']+conf['usrp_tx']['repetitions_2']-1))
296 305
297 306 if conf['usrp_tx']['code_type_2']:
298 307 codes = [ c.strip() for c in conf['usrp_tx']['code_2'].split(',')]
299 308 code = []
300 309 for c in codes:
301 310 code.append([int(x) for x in c])
302 311 op = voltage2.addOperation(name='Decoder', optype='other')
303 312 op.addParameter(name='code', value=code)
304 313 op.addParameter(name='nCode', value=len(code), format='int')
305 314 op.addParameter(name='nBaud', value=len(code[0]), format='int')
306 315
307 316 op = voltage2.addOperation(name='CohInt', optype='other') #Minimo integrar 2 perfiles por ser codigo complementario
308 317 op.addParameter(name='n', value=len(code), format='int')
309 318 ncode = len(code)
310 319 else:
311 320 ncode = 1
312 321
313 322 op = voltage2.addOperation(name='setH0')
314 323 op.addParameter(name='h0', value=H0, format='float')
315 324
316 325 if args.range > 0:
317 326 print('voltage 2', flush=True)
318 327 op = voltage2.addOperation(name='selectHeights')
319 328 op.addParameter(name='minIndex', value=max_index(RMIX, sample_rate, ipp, H0), format='int')
320 329 op.addParameter(name='maxIndex', value=max_index(args.range, sample_rate, ipp, H0), format='int')
321 330
322 331 #op = voltage2.addOperation(name='setH0')
323 332 #op.addParameter(name='h0', value=H0, format='float')
324 333
325 334 op = voltage2.addOperation(name='PulsePair_vRF', optype='other')
326 335 op.addParameter(name='n', value=int(conf['usrp_tx']['repetitions_2'])/ncode, format='int')
327 336 if args.rmDC:
328 337 op.addParameter(name='removeDC', value=1, format='int')
329 338
330 339 proc2 = project.addProcUnit(datatype='ParametersProc', inputId=voltage2.getId())
331 340 proc2.addParameter(name='runNextUnit', value=True)
332 341
333 342 opObj10 = proc2.addOperation(name="WeatherRadar")
334 343 opObj10.addParameter(name='CR_Flag',value=True,format='bool')
335 344 opObj10.addParameter(name='tauW',value=(1e-6/sample_rate)*len(code[0]))
336 345 #opObj10.addParameter(name='Pt',value=((1e-6/sample_rate)*len(code[0])/ipp)*200)
337 346 opObj10.addParameter(name='Pt',value=200)
338 347 opObj10.addParameter(name='min_index',value=max_index(RMIX, sample_rate, ipp, H0))
339 348 #opObj10.addParameter(name='sesgoZD',value=7.73)
340 349
341 350 op = proc2.addOperation(name='PedestalInformation')
342 351 op.addParameter(name='path', value=path_ped, format='str')
343 352 op.addParameter(name='interval', value='0.04')
344 353 op.addParameter(name='time_offset', value=time_offset)
345 354 op.addParameter(name='mode', value=args.mode)
346 355
347 356 op = proc2.addOperation(name='Block360')
348 357 op.addParameter(name='attr_data', value='data_param')
349 358 op.addParameter(name='runNextOp', value=True)
350 359 op.addParameter(name='angles', value=angles)
351 360 #op.addParameter(name='horario',value=False)
352 361
353 362 merge = project.addProcUnit(datatype='MergeProc', inputId=[proc1.getId(), proc2.getId()])
354 363 merge.addParameter(name='attr_data', value='data_param')
355 364 merge.addParameter(name='mode', value='7') #RM
356 365
357
358 366 for param in parameters:
359 367
360 368 if args.plot:
361 369 op= merge.addOperation(name='WeatherParamsPlot')
362 370 if args.save:
363 371 op.addParameter(name='save', value=path_plots, format='str')
364 372 op.addParameter(name='save_period', value=-1)
365 373 op.addParameter(name='show', value=args.show)
366 374 #op.addParameter(name='channels', value='0,1')
367 375 op.addParameter(name='channels', value='0,')
368 376 op.addParameter(name='zmin', value=PARAM[param]['zmin'], format='int')
369 377 op.addParameter(name='zmax', value=PARAM[param]['zmax'], format='int')
370 op.addParameter(name='ymax', value=20, format='int')
371 op.addParameter(name='xmin', value=-50, format='int')
372 op.addParameter(name='xmax', value=50, format='int')
378 op.addParameter(name='yrange', value=20, format='int')
379 op.addParameter(name='xrange', value=args.range, format='int')
373 380 op.addParameter(name='attr_data', value=param, format='str')
374 381 op.addParameter(name='labels', value=[[PARAM[param]['label']], [PARAM[param]['label']]])
375 382 op.addParameter(name='save_code', value=param)
376 383 op.addParameter(name='cb_label', value=PARAM[param]['cb_label'])
377 384 op.addParameter(name='colormap', value=PARAM[param]['colormap'])
378 385 op.addParameter(name='bgcolor', value='black')
379 386 op.addParameter(name='localtime', value=False)
380 387 op.addParameter(name='shapes', value='./shapes')
388 op.addParameter(name='latitude', value=conf['latitude'], format='float')
389 op.addParameter(name='longitude', value=conf['longitude'], format='float')
390 op.addParameter(name='map', value=False)
391
381 392 if MASK: op.addParameter(name='mask', value=MASK, format='float')
382 393 if args.server:
383 394 op.addParameter(name='server', value='190.187.237.239:4444')
384 395 op.addParameter(name='exp_code', value='400')
385 396
386 397 desc = {
387 398 'Data': {
388 399 'data_param': {PARAM[param]['wrname']: ['H', 'V']},
389 400 'utctime': 'time'
390 401 },
391 402 'Metadata': {
392 403 'heightList': 'range',
393 404 'data_azi': 'azimuth',
394 405 'data_ele': 'elevation',
395 406 'mode_op': 'scan_type',
396 407 'h0': 'range_correction',
397 408 'dataPP_NOISE': 'noise',
398 409 }
399 410 }
400 411
401 412 if args.save:
402 413 writer = merge.addOperation(name='HDFWriter')
403 414 writer.addParameter(name='path', value=path_save, format='str')
404 415 writer.addParameter(name='Reset', value=True)
405 416 writer.addParameter(name='setType', value='weather')
406 417 writer.addParameter(name='setChannel', value='0') #new parameter choose ch 0 H or ch 1 V
407 418 writer.addParameter(name='description', value=json.dumps(desc))
408 419 writer.addParameter(name='blocksPerFile', value='1',format='int')
409 420 writer.addParameter(name='metadataList', value=','.join(META))
410 421 writer.addParameter(name='dataList', value='data_param,utctime')
411 422 writer.addParameter(name='weather_var', value=param)
412 423 writer.addParameter(name='mask', value=MASK, format='float')
413 424 writer.addParameter(name='localtime', value=False)
414 425 # meta
415 writer.addParameter(name='latitude', value='-12.040436')
416 writer.addParameter(name='longitude', value='-75.295893')
417 writer.addParameter(name='altitude', value='3379.2147')
418 writer.addParameter(name='heading', value='0')
426 writer.addParameter(name='latitude', value=conf['latitude'])
427 writer.addParameter(name='longitude', value=conf['longitude'])
428 writer.addParameter(name='altitude', value=conf['altitude'])
429 writer.addParameter(name='heading', value=conf['heading'])
419 430 writer.addParameter(name='radar_name', value='SOPHy')
420 431 writer.addParameter(name='institution', value='IGP')
421 432 writer.addParameter(name='contact', value='dscipion@igp.gob.pe')
422 433 writer.addParameter(name='created_by', value='Signal Chain (https://pypi.org/project/schainpy/)')
423 434 writer.addParameter(name='range_unit', value='km')
424 435 writer.addParameter(name='prf', value=1/ipp)
425 436 writer.addParameter(name='prf_unit', value='hertz')
426 437 writer.addParameter(name='variable', value=PARAM[param]['label'])
427 438 writer.addParameter(name='variable_unit', value=PARAM[param]['cb_label'])
428 439 writer.addParameter(name='n_pulses', value=n_pulses)
429 440 writer.addParameter(name='pulse1_range', value=RMIX)
430 441 writer.addParameter(name='pulse1_width', value=pulse_1_width)
431 442 writer.addParameter(name='pulse2_width', value=pulse_2_width)
432 443 writer.addParameter(name='pulse1_repetitions', value=pulse_1_repetitions)
433 444 writer.addParameter(name='pulse2_repetitions', value=pulse_2_repetitions)
434 445 writer.addParameter(name='pulse_width_unit', value='microseconds')
435 446 writer.addParameter(name='snr_threshold', value=MASK)
436 447
437 448 return project
438 449
439 450 if __name__ == '__main__':
440 451
441 452 parser = argparse.ArgumentParser(description='Script to process SOPHy data.')
442 453 parser.add_argument('experiment',
443 454 help='Experiment name')
444 455 parser.add_argument('--parameters', nargs='*', default=['S'],
445 456 help='Variables to process: P, Z, V')
446 457 parser.add_argument('--angles', nargs='*', default=[], type=int,
447 458 help='Angles to process')
448 459 parser.add_argument('--time_offset', default=0,
449 460 help='Fix time offset')
450 parser.add_argument('--range', default=0, type=float,
461 parser.add_argument('--range', default=60, type=float,
451 462 help='Max range to plot')
452 463 parser.add_argument('--save', action='store_true',
453 464 help='Create output files')
454 465 parser.add_argument('--plot', action='store_true',
455 466 help='Create plot files')
456 467 parser.add_argument('--show', action='store_true',
457 468 help='Show matplotlib plot.')
458 469 parser.add_argument('--online', action='store_true',
459 470 help='Set online mode.')
460 471 parser.add_argument('--server', action='store_true',
461 472 help='Send to realtime')
462 473 parser.add_argument('--start_time', default='',
463 474 help='Set start time.')
464 475 parser.add_argument('--end_time', default='',
465 476 help='Set end time.')
466 477 parser.add_argument('--label', default='',
467 478 help='Label for plot & param folder')
468 479 parser.add_argument('--mode', default=None,
469 480 help='Type of scan')
470 481 parser.add_argument('--rmDC', action='store_true',
471 482 help='Apply remove DC.')
472 483 args = parser.parse_args()
473 484
474 485 project = main(args)
475 486 project.start()
476 487
477 488 #python sophy_A.py HYO_CC4_CC64_COMB@2022-12-27T00-00-32 --parameters Z --plot --save --show --rmDC --label Z_04 --range 60 --start_time "22:00:00"
478 489 # colocar siempre el range que asume 0 y no hace la seleccion de alturas
General Comments 0
You need to be logged in to leave comments. Login now