schain Files · schainpy/scripts/check_pedestal.py

Script test de pruebas actuales con el update de heading en el block 360, el parametro adicional es heading que se lee del experiment.conf

avaldezp - - Load All Authors

File last commit:

r1429:96f17ea26f9f


                r1759:c8ee9da660f1

Download file

             check_pedestal.py
        
                    133 lines
            
             | 4.2 KiB
            
                | text/x-python
            
             |
                PythonLexer
            
             / schainpy / scripts / check_pedestal.py
          
                    History
                
                 |
                  Source
                 | Raw
                 |Copy content
                 |Copy permalink

        avaldez
    
LAST UPDATE INTEGRACION

              r1388
            
      # Ing-AlexanderValdez

      # Monitoreo de Pedestal

      ############## IMPORTA LIBRERIAS ###################

      import os,numpy,h5py

      import sys,time

      import matplotlib.pyplot as plt

      ####################################################

        avaldez
    
LAST CHANGE WR

              r1393
            
      #################################################################

      # 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'

        avaldezp
    
last update

              r1427
            
      #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"

        avaldezp
    
last script test

              r1428
            
      #path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-163824"

        avaldezp
    
last 30-16:19

              r1405
            
      #path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173409'

        avaldezp
    
last script test

              r1428
            
      #--------------------------------

      path_ped= "/DATA_RM/TEST_PEDESTAL/P20220401-172744"

        avaldezp
    
update

              r1429
            
        avaldez
    
LAST UPDATE INTEGRACION

              r1388
            
      # 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):

        avaldez
    
LAST CHANGE WR

              r1393
            
          print("order:",i)

        avaldez
    
LAST UPDATE INTEGRACION

              r1388
            
          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")

        avaldezp
    
last_update primeras correciones rhi

              r1418
            
          tmp_ele_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="ele_vel")# nuevo :D

        avaldez
    
LAST UPDATE INTEGRACION

              r1388
            
          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:

        avaldezp
    
update

              r1429
            
          #tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_pos")

        avaldezp
    
last script test

              r1428
            
          tmp_ele_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_pos")

        avaldezp
    
update

              r1429
            
          tmp_azi_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_vel")

        avaldezp
    
last script test

              r1428
            
          #tmp_ele_pos   = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_vel")

        avaldez
    
LAST UPDATE INTEGRACION

              r1388
            
          # 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

        avaldezp
    
last update

              r1412
            
      ax.set_xlim(0, 4)

        avaldezp
    
last update

              r1427
            
      ax.set_ylim(-5, 360)

        avaldez
    
LAST UPDATE INTEGRACION

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

	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

avaldez LAST UPDATE INTEGRACION	r1388	# Ing-AlexanderValdez
		# Monitoreo de Pedestal

		############## IMPORTA LIBRERIAS ###################
		import os,numpy,h5py
		import sys,time
		import matplotlib.pyplot as plt
		####################################################
avaldez LAST CHANGE WR	r1393	#################################################################
		# 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'
avaldezp last update	r1427	#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"
avaldezp last script test	r1428	#path_ped = "/DATA_RM/TEST_PEDESTAL/P20220322-163824"
avaldezp last 30-16:19	r1405	#path_ped = '/DATA_RM/TEST_PEDESTAL/P20211111-173409'
avaldezp last script test	r1428

		#--------------------------------

		path_ped= "/DATA_RM/TEST_PEDESTAL/P20220401-172744"
avaldezp update	r1429
avaldez LAST UPDATE INTEGRACION	r1388	# 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):
avaldez LAST CHANGE WR	r1393	print("order:",i)
avaldez LAST UPDATE INTEGRACION	r1388	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")
avaldezp last_update primeras correciones rhi	r1418	tmp_ele_vel = getDatavaluefromDirFilename(path=path_ped,file=LIST[i],value="ele_vel")# nuevo :D
avaldez LAST UPDATE INTEGRACION	r1388
		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:
avaldezp update	r1429	#tmp_azi_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_pos")
avaldezp last script test	r1428	tmp_ele_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_pos")
avaldezp update	r1429	tmp_azi_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="ele_vel")
avaldezp last script test	r1428	#tmp_ele_pos = getDatavaluefromDirFilename(path=path_ped,file=filename,value="azi_vel")
avaldez LAST UPDATE INTEGRACION	r1388	# 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

avaldezp last update	r1412	ax.set_xlim(0, 4)
avaldezp last update	r1427	ax.set_ylim(-5, 360)
avaldez LAST UPDATE INTEGRACION	r1388	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()