# 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()