schain Commit - r1412:7050d08f9cdf · Jicamarca Repository

last update

avaldezp -

r1412:7050d08f9cdf

parent child

Context file:

r1412:7050d08f9cdf

Collapse all files

The requested changes are too big and content was truncated. Show full diff

schainpy/model/proc/jroproc_parameters.py 0 0

	1	NO CONTENT: modified file			NO CONTENT: modified file
The requested commit or file is too big and content was truncated. Show full diff

schainpy/scripts/basic_proc00004.py +2 -2

             # Ing. AVP
             # 04/01/2022
             # ARCHIVO DE LECTURA
             #---- DATA RHI --- 23 DE NOVIEMBRE DEL 2021 --- 23/11/2021---
             #---- PEDESTAL ----------------------------------------------
             #------- HORA 143826 /DATA_RM/TEST_PEDESTAL/P20211123-143826 14:38-15:10
             #---- RADAR ----------------------------------------------
             #-------  14:26-15:00
             #-------  /DATA_RM/DRONE/2MHZ_5V_ELEVACION/
             #-------  /DATA_RM/DRONE/2MHZ_5V_ELEVACION/ch0/2021-11-23T19-00-00
             import os, sys
             import datetime
             import time
             import numpy
             from  ext_met import  getfirstFilefromPath,getDatavaluefromDirFilename
             from schainpy.controller import Project
             #-----------------------------------------------------------------------------------------
             print("[SETUP]-RADAR METEOROLOGICO-")
             path_ped  = "/DATA_RM/TEST_PEDESTAL/P20211123-143826"
             print("PATH PEDESTAL :",path_ped)
             path_adq  = "/DATA_RM/DRONE/2MHZ_5V_ELEVACION/"
             print("PATH DATA     :",path_adq)
             figpath_pp_rti  = "/home/soporte/Pictures/TEST_PP_RHI"
             print("PATH PP RTI   :",figpath_pp_rti)
             figpath_pp_rhi  = "/home/soporte/Pictures/TEST_PP_RHI"
             print("PATH PP RHI   :",figpath_pp_rhi)
             path_pp_save_int  = "/DATA_RM/TEST_SAVE_PP_INT_RHI"
             print("PATH SAVE PP INT   :",path_pp_save_int)
             print(" ")
             #-------------------------------------------------------------------------------------------
             print("SELECCIONAR MODO: PPI        (0) O RHI       (1)")
             mode_wr     = 1
             if mode_wr==0:
                 print("[ ON ] MODE PPI")
                 list_ped     = getfirstFilefromPath(path=path_ped,meta="PE",ext=".hdf5")
                 ff_pedestal  = list_ped[2]
                 azi_vel      = getDatavaluefromDirFilename(path=path_ped,file=ff_pedestal,value="azi_vel")
                 V            = round(azi_vel[0])
                 print("VELOCIDAD AZI :", int(numpy.mean(azi_vel)),"°/seg")
             else:
                 print("[ ON ] MODE RHI")
                 list_ped     = getfirstFilefromPath(path=path_ped,meta="PE",ext=".hdf5")
                 ff_pedestal  = list_ped[2]
                 ele_vel      = getDatavaluefromDirFilename(path=path_ped,file=ff_pedestal,value="ele_vel")
                 V            = round(ele_vel[0])
                 V            = 10.0
                 print("VELOCIDAD ELE :", int(numpy.mean(ele_vel)),"°/seg")
             print(" ")
             #---------------------------------------------------------------------------------------
             print("SELECCIONAR MODO: PULSE PAIR (0) O FREQUENCY (1)")
             mode_proc   = 0
             if mode_proc==0:
                 print("[ ON ] MODE PULSEPAIR")
             else:
                 print("[ ON ] MODE FREQUENCY")
             ipp  = 60.0
             print("IPP(Km.)      :  %1.2f"%ipp)
             ipp_sec = (ipp*1.0e3/150.0)*1.0e-6
             print("IPP(useg.)    : %1.2f"%(ipp_sec*(1.0e6)))
             VEL=V
             n= int(1/(VEL*ipp_sec))
             print("N° Profiles   : ", n)
             #--------------------------------------------
             plot_rti    = 0
             plot_rhi    = 1
             integration = 1
             save        = 0
             #---------------------------RANGO DE PLOTEO----------------------------------
             dBmin = '1'
             dBmax = '85'
             xmin = '17'
             xmax = '17.25'
             ymin = '0'
             ymax = '600'
             #----------------------------------------------------------------------------
             time.sleep(3)
             #---------------------SIGNAL CHAIN ------------------------------------
             desc = "USRP_WEATHER_RADAR"
             filename = "USRP_processing.xml"
             controllerObj = Project()
             controllerObj.setup(id = '191', name='Test_USRP', description=desc)
             #---------------------UNIDAD DE LECTURA--------------------------------
             readUnitConfObj = controllerObj.addReadUnit(datatype='DigitalRFReader',
                                                         path=path_adq,
-                                                        startDate="2021/11/10",#today,
+                                                        startDate="2021/11/23",#today,
                                                         endDate="2021/12/30",#today,
-                                                        startTime='17:10:25',
+                                                        startTime='14:38:23',
                                                         endTime='23:59:59',
                                                         delay=0,
                                                         #set=0,
                                                         online=0,
                                                         walk=1,
                                                         ippKm=ipp)
             procUnitConfObjA = controllerObj.addProcUnit(datatype='VoltageProc',inputId=readUnitConfObj.getId())
             opObj11 = procUnitConfObjA.addOperation(name='selectHeights')
             opObj11.addParameter(name='minIndex', value='1', format='int')
             #    opObj11.addParameter(name='maxIndex', value='10000', format='int')
             opObj11.addParameter(name='maxIndex', value='400', format='int')
             if mode_proc==0:
                 opObj11 = procUnitConfObjA.addOperation(name='PulsePair', optype='other')
                 opObj11.addParameter(name='n', value=int(n), format='int')
                 procUnitConfObjB= controllerObj.addProcUnit(datatype='ParametersProc',inputId=procUnitConfObjA.getId())
                 if integration==1:
                     opObj11 = procUnitConfObjB.addOperation(name='PedestalInformation')
                     opObj11.addParameter(name='path_ped', value=path_ped)
                     opObj11.addParameter(name='t_Interval_p', value='0.01', format='float')
                 if plot_rhi==1:
                     opObj11 = procUnitConfObjB.addOperation(name='Block360')
                     opObj11.addParameter(name='n', value='10', format='int')
                     opObj11.addParameter(name='mode', value=mode_proc, format='int')
                     # este bloque funciona bien con divisores de 360 no olvidar 0 10 20 30 40 60 90 120 180
                     opObj11= procUnitConfObjB.addOperation(name='WeatherRHIPlot',optype='other')
                     opObj11.addParameter(name='save', value=figpath_pp_rhi)
                     opObj11.addParameter(name='save_period', value=1)
             controllerObj.start()

schainpy/scripts/check_pedestal.py +3 -2

             # Ing-AlexanderValdez
             # Monitoreo de Pedestal
             ############## IMPORTA LIBRERIAS ###################
             import os,numpy,h5py
             import sys,time
             import matplotlib.pyplot as plt
             ####################################################
             #################################################################
             # LA FECHA 21-10-20 CORRESPONDE A LAS PRUEBAS DEL DIA MIERCOLES
             # 1:15:51 pm hasta 3:49:32 pm
             #################################################################
             #path_ped = '/DATA_RM/TEST_PEDESTAL/P20211012-082745'
             path_ped = '/DATA_RM/TEST_PEDESTAL/P20211020-131248'
             path_ped = '/DATA_RM/TEST_PEDESTAL/P20211110-171003'
             path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173856'
+            path_ped = '/DATA_RM/TEST_PEDESTAL/P20211123-143826'
             #path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173409'
             # Metodo para verificar numero
             def isNumber(str):
                 try:
                     float(str)
                     return True
                 except:
                     return False
             # Metodo para extraer el arreglo
             def getDatavaluefromDirFilename(path,file,value):
                 dir_file= path+"/"+file
                 fp      = h5py.File(dir_file,'r')
                 array    = fp['Data'].get(value)[()]
                 fp.close()
                 return array
             # LISTA COMPLETA DE ARCHIVOS HDF5 Pedestal
             LIST= sorted(os.listdir(path_ped))
             m=len(LIST)
             print("TOTAL DE ARCHIVOS DE PEDESTAL:",m)
             # Contadores temporales
             k=  0
             l=  0
             t=  0
             # Marca de tiempo temporal
             time_        = numpy.zeros([m])
             # creacion de
             for i in range(m):
                 print("order:",i)
                 tmp_azi_pos = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="azi_pos")
                 tmp_ele_pos = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="ele_pos")
                 tmp_azi_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="azi_vel")
                 tmp_ele_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="azi_vel")# nuevo :D
                 time_[i]    = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="utc")
                 k=k +tmp_azi_pos.shape[0]
                 l=l +tmp_ele_pos.shape[0]
                 t=t +tmp_azi_vel.shape[0]
             print("TOTAL DE MUESTRAS, ARCHIVOS X100:",k)
             time.sleep(5)
             ######CREACION DE ARREGLOS CANTIDAD DE VALORES POR MUESTRA#################
             azi_pos     = numpy.zeros([k])
             ele_pos     = numpy.zeros([l])
             time_azi_pos= numpy.zeros([k])
             # Contadores temporales
             p=0
             r=0
             z=0
             # VARIABLES TMP para almacenar azimuth, elevacion y tiempo
             #for filename in sorted(os.listdir(path_ped)):
             # CONDICION POR LEER EN TIEMPO REAL NO OFFLINE
             for filename in LIST:
                 tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_pos")
                 tmp_ele_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_pos")
                 # CONDICION POR LEER EN TIEMPO REAL NO OFFLINE
                 if z==(m-1):
                     tmp_azi_time=numpy.arange(time_[z],time_[z]+1,1/(tmp_azi_pos.shape[0]))
                 else:
                     tmp_azi_time=numpy.arange(time_[z],time_[z+1],(time_[z+1]-time_[z])/(tmp_azi_pos.shape[0]))
                 print(filename,time_[z])
                 print(z,tmp_azi_pos.shape[0])
                 i=0
                 for i in range(tmp_azi_pos.shape[0]):
                     index=p+i
                     azi_pos[index]=tmp_azi_pos[i]
                     time_azi_pos[index]=tmp_azi_time[i]
                 p=p+tmp_azi_pos.shape[0]
                 i=0
                 for i in range(tmp_ele_pos.shape[0]):
                     index=r+i
                     ele_pos[index]=tmp_ele_pos[i]
                 r=r+tmp_ele_pos.shape[0]
                 z+=1
             ######## GRAFIQUEMOS Y VEAMOS LOS DATOS DEL Pedestal
             fig, ax = plt.subplots(figsize=(16,8))
             print(time_azi_pos.shape)
             print(azi_pos.shape)
             t=numpy.arange(time_azi_pos.shape[0])*0.01/(60.0)
             plt.plot(t,azi_pos,label='AZIMUTH_POS',color='blue')
             # AQUI ESTOY ADICIONANDO LA POSICION EN elevaciont=numpy.arange(len(ele_pos))*0.01/60.0
             t=numpy.arange(len(ele_pos))*0.01/60.0
             plt.plot(t,ele_pos,label='ELEVATION_POS',color='red')#*10
-            #ax.set_xlim(0, 9)
+            ax.set_xlim(0, 4)
-            ax.set_ylim(-5, 20)
+            ax.set_ylim(-5, 50)
             plt.ylabel("Azimuth Position")
             plt.xlabel("Muestra")
             plt.title('Azimuth Position vs Muestra ', fontsize=20)
             axes = plt.gca()
             axes.yaxis.grid()
             plt.xticks(fontsize=16)
             plt.yticks(fontsize=16)
             plt.show()

schainpy/scripts/test_wradlib_RHI.py +1 0

             import numpy as np
             import matplotlib.pyplot as plt
             import wradlib as wrl
             import warnings
             # libreia nueva
+            #export WRADLIB_DATA="/home/soporte/Downloads/2014-06-09--185000.rhi.mvol"
             from mpl_toolkits.axisartist.grid_finder import FixedLocator, DictFormatter
             warnings.filterwarnings('ignore')
             # lectura de gaMIC hdf5 file
             filename        = wrl.util.get_wradlib_data_file("/home/soporte/Downloads/2014-06-09--185000.rhi.mvol")
             data1, metadata = wrl.io.read_gamic_hdf5(filename)
             print(data1)
             data1 = data1['SCAN0']['ZH']['data']
             print(data1)
             print("SHAPE Data",np.array(data1).shape)
             r = metadata['SCAN0']['r']
             print("r",r)
             print("longitud r",len(r))
             th = metadata['SCAN0']['el']
             print("th",th)
             print("longitud th",len(th))
             az = metadata['SCAN0']['az']
             print("az",az)
             site = (metadata['VOL']['Longitude'], metadata['VOL']['Latitude'],
                     metadata['VOL']['Height'])
             print("Longitud,Latitud,Altura",site)
             ma1 =  np.array(data1)
             '''
             mask_ind = np.where(data1 <= np.nanmin(data1))
             data1[mask_ind] = np.nan
             ma1 = np.ma.array(data1, mask=np.isnan(data1))
             '''
             #cgax, pm = wrl.vis.plot_rhi(ma1,r=r,th=th,rf=1e3)
             fig = plt.figure(figsize=(10,8))
             cgax, pm = wrl.vis.plot_rhi(ma1,r=r,th=th,rf=1e3,fig=fig, ax=111,proj='cg')
             caax = cgax.parasites[0]
             paax = cgax.parasites[1]
             cgax.set_ylim(0, 14)
             #caax = cgax.parasites[0]
             #paax = cgax.parasites[1]
             #cgax, pm = wrl.vis.plot_rhi(ma1, r=r, th=th, rf=1e3, fig=fig, ax=111, proj='cg')
             txt = plt.title('Simple RHI',y=1.05)
             #cbar = plt.gcf().colorbar(pm, pad=0.05, ax=paax)
             cbar = plt.gcf().colorbar(pm, pad=0.05)
             cbar.set_label('reflectivity [dBZ]')
             caax.set_xlabel('x_range [km]')
             caax.set_ylabel('y_range [km]')
             plt.text(1.0, 1.05, 'azimuth', transform=caax.transAxes, va='bottom',ha='right')
             gh = cgax.get_grid_helper()
             # set theta to some nice values
             locs = [0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13., 14.,
 ., 16., 17., 18., 20., 22., 25., 30., 35.,  40., 50., 60., 70., 80., 90.]
             gh.grid_finder.grid_locator1 = FixedLocator(locs)
             gh.grid_finder.tick_formatter1 = DictFormatter(dict([(i, r"${0:.0f}^\circ$".format(i)) for i in locs]))
             plt.show()

General Comments 0

Write
Preview

You need to be logged in to leave comments. Login now

No TODOs yet

	Site-wide shortcuts
/	Use quick search box
g h	Goto home page
g g	Goto my private gists page
g G	Goto my public gists page
g 0-9	Goto bookmarked items from 0-9
n r	New repository page
n g	New gist page

	Repositories
g s	Goto summary page
g c	Goto changelog page
g f	Goto files page
g F	Goto files page with file search activated
g p	Goto pull requests page
g o	Goto repository settings
g O	Goto repository access permissions settings
t s	Toggle sidebar on some pages