##// END OF EJS Templates
schain
RSS
update WeatherRadar
update WeatherRadar
avaldezp - - Load All Authors

File last commit:

r1433:24097e0befe4
r1508:6520a1ac8270
Show More
check_pedestal2.py
164 lines | 5.2 KiB | text/x-python | PythonLexer
/ schainpy / scripts / check_pedestal2.py
History | Source | Raw |Copy content |Copy permalink
avaldezp
check rhi
 r1433 # 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/P20220217-172216"
#path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-163824"
#path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173409'
#--------------------------------
#path_ped= "/DATA_RM/TEST_PEDESTAL/P20220401-172744"
path_ped="/DATA_RM/TEST_PEDESTAL/P20220412-140748"
path_ped='/DATA_RM/DATA/TEST@2022-04-12T14:05:11/position/2022-04-12T19-00-00'
# 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
u= 0
# Marca de tiempo temporal
time_ = numpy.zeros([m])
# creacion de
flag=0
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_speed")#azi_speed azi_vel
tmp_ele_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="ele_speed")# nuevo :D ele_speed,ele_vel
try:
time_[i] = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="utc")
except:
time___ = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="utc")
flag=1
k=k +tmp_azi_pos.shape[0]
l=l +tmp_ele_pos.shape[0]
t=t +tmp_azi_vel.shape[0]
if flag==1:
u=u +time___.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])
time__ = numpy.zeros([u])
# Contadores temporales
p=0
r=0
v=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")
#tmp_azi_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_speed")
#tmp_ele_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_speed")
tmp_time = getDatavaluefromDirFilename(path=path_ped,file=filename,value="utc")
# CONDICION POR LEER EN TIEMPO REAL NO OFFLINE
if not flag==1:
if z==(m-1):
tmp_azi_time=numpy.arange(time_[z],time_[z]+1,1.0/(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]
if not flag==1:
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]
i=0
for i in range(tmp_time.shape[0]):
index=v+i
time__[index]=tmp_time[i]
v=v+tmp_time.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)
#print("time__",time__)
##import matplotlib.dates as mdate
##import matplotlib.ticker as mtick
##secs = mdate.epoch2num(time__)
##print(secs)
##plt.plot(secs,azi_pos,label='AZIMUTH_POS',color='blue')
diff= time__[2]-time__[1]
time__=time__ - time__[0]
time__= time__/60.0
#print("diff",diff)
plt.plot(time__,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(time__,ele_pos,label='ELEVATION_POS',color='red')#*10
ax.set_xlim(0, 10)
#ax.set_ylim(-5, 365)
plt.ylabel("Azimuth Position")
plt.xlabel("Muestra")
plt.title('Azimuth Position vs timepo en min ', fontsize=20)
axes = plt.gca()
axes.yaxis.grid()
plt.xticks(fontsize=16)
plt.yticks(fontsize=16)
plt.show()