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