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