| @@ -7,15 +7,18 PATH = '/DATA_RM/DATA' | |||
|
|
7 | 7 | #PATH = '/media/jespinoza/Elements' |
|
|
8 | 8 | #PATH = '/media/jespinoza/data/SOPHY' |
|
|
9 | 9 | PATH = '/home/soporte/Documents/HUANCAYO' |
|
|
10 | PATH = '/home/soporte/Documents/EVENTO/' | |
|
|
11 | #HYO_PM@2022-05-31T12-00-17 | |
|
|
10 | 12 | |
|
|
11 | 13 | PARAM = { |
|
|
12 | 14 | 'S': {'name': 'dataPP_POWER', 'zmin': -45, 'zmax': -15, 'colormap': 'jet', 'label': 'Power', 'wrname': 'power','cb_label': 'dBm', 'ch':0}, |
|
|
13 |
'V': {'name': 'dataPP_DOP', 'zmin': -10, 'zmax': 10, 'colormap': ' |
|
|
|
15 | #'V': {'name': 'dataPP_DOP', 'zmin': -10, 'zmax': 10, 'colormap': 'sophy_v', 'label': 'Velocity', 'wrname': 'velocity', 'cb_label': 'm/s', 'ch':0}, | |
|
|
16 | 'V': {'name': 'velRadial_V', 'zmin': -10, 'zmax': 10, 'colormap': 'sophy_v', 'label': 'Velocity', 'wrname': 'velocity', 'cb_label': 'm/s', 'ch':0}, | |
|
|
14 | 17 | 'R': {'name': 'RhoHV_R', 'zmin': 0, 'zmax': 1, 'colormap': 'jet', 'label': 'RhoHV', 'wrname':'rhoHV', 'cb_label': '*', 'ch':0}, |
|
|
15 | 18 | 'P': {'name': 'PhiD_P', 'zmin': -180,'zmax': 180,'colormap': 'RdBu_r', 'label': 'PhiDP', 'wrname':'phiDP' , 'cb_label': 'ΒΊ', 'ch':0}, |
|
|
16 | 19 | 'D': {'name': 'Zdb_D', 'zmin': -20, 'zmax': 80, 'colormap': 'gist_ncar','label': 'ZDR','wrname':'differential_reflectivity' , 'cb_label': 'dBz','ch':0}, |
|
|
17 |
'Z': {'name': 'Zdb', 'zmin': - |
|
|
|
18 |
'W': {'name': 'Sigmav_W', 'zmin': -20, 'zmax': 60, 'colormap': ' |
|
|
|
20 | 'Z': {'name': 'Zdb', 'zmin': -30, 'zmax': 80, 'colormap': 'sophy_r','label': 'Reflectivity', 'wrname':'reflectivity', 'cb_label': 'dBz','ch':1}, | |
|
|
21 | 'W': {'name': 'Sigmav_W', 'zmin': -20, 'zmax': 60, 'colormap': 'sophy_w','label': 'Spectral Width', 'wrname':'spectral_width', 'cb_label': 'hz', 'ch':1} | |
|
|
19 | 22 | } |
|
|
20 | 23 | |
|
|
21 | 24 | def max_index(r, sample_rate, ipp): |
| @@ -42,14 +45,14 def main(args): | |||
|
|
42 | 45 | start_time = args.start_time |
|
|
43 | 46 | else: |
|
|
44 | 47 | start_time = experiment.split('@')[1].split('T')[1].replace('-', ':') |
|
|
45 |
#start_time = '1 |
|
|
|
48 | #start_time = '16:15:00' | |
|
|
46 | 49 | end_time = '23:59:59' |
|
|
47 | 50 | N = int(1/(speed_axis[0]*ipp)) # 1 GRADO DE RESOLUCION |
|
|
48 | 51 | path = os.path.join(PATH, experiment, 'rawdata') |
|
|
49 | 52 | path_ped = os.path.join(PATH, experiment, 'position') |
|
|
50 | path_plots = os.path.join(PATH, experiment, 'plots') | |
|
|
53 | path_plots = os.path.join(PATH, experiment, 'plots_R5km_M1.62km') | |
|
|
51 | 54 | path_save = os.path.join(PATH, experiment, 'param') |
|
|
52 | RMIX = 1 | |
|
|
55 | RMIX = 1.62 | |
|
|
53 | 56 | |
|
|
54 | 57 | from schainpy.controller import Project |
|
|
55 | 58 | |
| @@ -60,7 +63,7 def main(args): | |||
|
|
60 | 63 | path=path, |
|
|
61 | 64 | startDate=start_date, |
|
|
62 | 65 | endDate=end_date, |
|
|
63 |
start |
|
|
|
66 | startTime=start_time, | |
|
|
64 | 67 | endTime=end_time, |
|
|
65 | 68 | delay=30, |
|
|
66 | 69 | online=args.online, |
| @@ -109,7 +112,7 def main(args): | |||
|
|
109 | 112 | op.addParameter(name='mode', value='PPI') |
|
|
110 | 113 | |
|
|
111 | 114 | for param in parameters: |
|
|
112 |
op = proc.addOperation(name='Block360 |
|
|
|
115 | op = proc.addOperation(name='Block360') | |
|
|
113 | 116 | op.addParameter(name='attr_data', value=PARAM[param]['name']) |
|
|
114 | 117 | op.addParameter(name='runNextOp', value=True) |
|
|
115 | 118 | |
| @@ -172,7 +175,7 def main(args): | |||
|
|
172 | 175 | op.addParameter(name='n', value=2, format='int') |
|
|
173 | 176 | |
|
|
174 | 177 | op = voltage1.addOperation(name='setH0') |
|
|
175 | op.addParameter(name='h0', value='-1.6') | |
|
|
178 | op.addParameter(name='h0', value='-1.68') | |
|
|
176 | 179 | |
|
|
177 | 180 | if args.range > 0: |
|
|
178 | 181 | op = voltage1.addOperation(name='selectHeights') |
| @@ -200,11 +203,11 def main(args): | |||
|
|
200 | 203 | op.addParameter(name='path', value=path_ped, format='str') |
|
|
201 | 204 | op.addParameter(name='interval', value='0.04') |
|
|
202 | 205 | op.addParameter(name='time_offset', value=time_offset) |
|
|
203 | op.addParameter(name='az_offset', value=-26.2) | |
|
|
206 | #op.addParameter(name='az_offset', value=-26.2) | |
|
|
204 | 207 | op.addParameter(name='mode', value='PPI') |
|
|
205 | 208 | |
|
|
206 | 209 | for param in parameters: |
|
|
207 |
op = proc1.addOperation(name='Block360 |
|
|
|
210 | op = proc1.addOperation(name='Block360') | |
|
|
208 | 211 | op.addParameter(name='attr_data', value=PARAM[param]['name']) |
|
|
209 | 212 | op.addParameter(name='runNextOp', value=True) |
|
|
210 | 213 | |
| @@ -231,7 +234,7 def main(args): | |||
|
|
231 | 234 | ncode = 1 |
|
|
232 | 235 | |
|
|
233 | 236 | op = voltage2.addOperation(name='setH0') |
|
|
234 | op.addParameter(name='h0', value='-1.6') | |
|
|
237 | op.addParameter(name='h0', value='-1.68') | |
|
|
235 | 238 | |
|
|
236 | 239 | if args.range > 0: |
|
|
237 | 240 | op = voltage2.addOperation(name='selectHeights') |
| @@ -244,7 +247,7 def main(args): | |||
|
|
244 | 247 | proc2 = project.addProcUnit(datatype='ParametersProc', inputId=voltage2.getId()) |
|
|
245 | 248 | |
|
|
246 | 249 | opObj10 = proc2.addOperation(name="WeatherRadar") |
|
|
247 | opObj10.addParameter(name='variableList',value='Reflectividad,AnchoEspectral') | |
|
|
250 | opObj10.addParameter(name='variableList',value='Reflectividad,VelocidadRadial,AnchoEspectral') | |
|
|
248 | 251 | print("tauw2",(1e-6/sample_rate)*len(code[0])) |
|
|
249 | 252 | print("Pt2",((1e-6/sample_rate)*len(code[0])/ipp)*200) |
|
|
250 | 253 | opObj10.addParameter(name='tauW',value=(1e-6/sample_rate)*len(code[0])) |
| @@ -255,11 +258,11 def main(args): | |||
|
|
255 | 258 | op.addParameter(name='path', value=path_ped, format='str') |
|
|
256 | 259 | op.addParameter(name='interval', value='0.04') |
|
|
257 | 260 | op.addParameter(name='time_offset', value=time_offset) |
|
|
258 | op.addParameter(name='az_offset', value=-26.2) | |
|
|
261 | #op.addParameter(name='az_offset', value=-26.2) | |
|
|
259 | 262 | op.addParameter(name='mode', value='PPI') |
|
|
260 | 263 | |
|
|
261 | 264 | for param in parameters: |
|
|
262 |
op = proc2.addOperation(name='Block360 |
|
|
|
265 | op = proc2.addOperation(name='Block360') | |
|
|
263 | 266 | #op.addParameter(name='axis', value=','.join(axis)) |
|
|
264 | 267 | op.addParameter(name='attr_data', value=PARAM[param]['name']) |
|
|
265 | 268 | op.addParameter(name='runNextOp', value=True) |
| @@ -280,6 +283,10 def main(args): | |||
|
|
280 | 283 | op.addParameter(name='save_code', value=param) |
|
|
281 | 284 | op.addParameter(name='cb_label', value=PARAM[param]['cb_label']) |
|
|
282 | 285 | op.addParameter(name='colormap', value=PARAM[param]['colormap']) |
|
|
286 | op.addParameter(name='bgcolor',value='black') | |
|
|
287 | op.addParameter(name='snr_threshold',value=0.4) | |
|
|
288 | ||
|
|
289 | ||
|
|
283 | 290 | |
|
|
284 | 291 | desc = { |
|
|
285 | 292 | 'Data': { |
General Comments 0
You need to be logged in to leave comments.
Login now