Herramientas para el procesamiento de datos GPS de la red LISN

2.2.2 Rinex.py
Este módulo engloba funciones y clases para leer archivos en formato Rinex y
guarda la información de los observables (Pseudoranges, Phases, Doppler Frequency y
SNR) por satélite y frecuencia observada, también se guarda la información de la cabecera
para ser utilizada posteriormente.
El módulo soporta archivos en formato RINEX (estándar *.yyo) y CRINEX (con
compresión hatanaka *.yyd), los cuales pueden incluso tener una segunda compresión
(*.gz, *.tar.gz).
Adicionalmente, la clase GPSdata tiene la función “calctec” que permite calcular el
TEC (slant TEC y TEC equivalente).
Este módulo está basado en el módulo “gpsdata” creado por Nick Matteo
kundor@kundor.org
read_rinex
Función
count_rinex
Descripción
Esta función lee un archivo en formato Rinex y crea un objeto del tipo
GPSData con la información de los observables encontrados.
Versión simplificada de read_rinex que solo cuenta el número de records que
contiene el archivo, no crea el objeto GPSData.
Tabla 3 Principales funciones del modulo Rinex.py
Clase
GPSData(list)
SatBias(dict)
LeapSeconds(dict)
Descripción
Objeto del tipo lista donde se almacenan los records, uno por cada tiempo de
observación, con información de los observables encontrados por satélite.
Tiene el atributo “header” con información de la cabecera del archivo.
Esta clase busca los archivos *.dcb que contienen información del bias de los
satélites y crea un diccionario con los valores del bias por satélite, sino
encuentra los archivos intenta bajarlos del sitio ftp que provee la UNIBE
(ftp.unibe.ch/aiub/CODE).
Usa el archivo leapseco.dat para formar un diccionario con las fechas y los
segundos intercalares, sino encuentra el archivo intenta bajar y actualizar el
archivo desde la web.
Tabla 4 Principales clases del modulo Rinex.py
2.2.2.1 Función “calctec”
Esta función pertenece a la clase “GPSData” y permite el cálculo del TEC a partir
de los datos observables, esta función está basada en el programa “tec_fromRT.f” provisto
por el Boston College.
Las principales características de esta función son:
Calcula el TEC relativo (utilizando los pseudoranges) y absoluto (utilizando las
phases) por satélite.
Calcula la elevación, azimut, latitud y longitud de cada satélite observado,
utilizando la clase “gpselem” del módulo “gpsProcess.py”.
Calcula el slant TEC.
Corrige los saltos y valores erróneos encontrados en los datos.
Corrige los cycle-slip encontrados.
Estima el bias del receptor, minimizando la varianza del TEC equivalente entre
las 3 y 6 a.m (hora local). Para la minimización se utiliza el método de Brent.
Corrige el bias del receptor hallado de tal forma que el TEC equivalente sea
siempre mayor a cero.
7