Traducido del inglés al español - www.onlinedoctranslator.

com

Soluciones GPS (2020) 24:94

https://doi.org/10.1007/s10291-020-01010-0

CAJA DE HERRAMIENTAS GPS

Herramientas de software Python para reflectometría interferométrica GNSS

(GNSS-IR)

Ángel Martín1· Raquel Luján1· Ana Belén Anquela1

Recibido: 15 de mayo de 2020 / Aceptado: 7 de julio de 2020

© Springer-Verlag GmbH Alemania, parte de Springer Nature 2020

Abstracto
La reflectometría interferométrica del Sistema Global de Navegación por Satélite (GNSS), también conocida como GNSS-IR, utiliza datos de
antenas GNSS de calidad geodésica para extraer información sobre el entorno que rodea a la antena. El monitoreo de la humedad del suelo
es una de las aplicaciones más importantes de la técnica GNSS-IR. Este manuscrito presenta las principales ideas y decisiones de
implementación necesarias para escribir el código Python para herramientas de software que transforman archivos de observación y
navegación en formato RINEX a un formato apropiado para GNSS-IR (que incluye las observaciones SNR y el acimut y elevación de los
satélites) y para determinar la altura de reflexión y los valores ajustados de fase y amplitud de la onda interferométrica para cada pista de
satélite individual. El objetivo principal del manuscrito es compartir el software con la comunidad científica para presentar a nuevos usuarios
la técnica GNSS-IR.

Palabras claveReflectometría GNSS-IR · Software Python · Humedad del suelo · Relación señal-ruido (SNR)

Introducción independiente de las condiciones climáticas (la técnica es válida en

condiciones de lluvia y niebla) y de iluminación (de día o de noche).
Mediante el uso de una antena geodésica del Sistema Global de Las señales de los satélites GNSS se transmiten en la banda
Navegación por Satélite (GNSS) para monitorear la humedad del suelo L (frecuencia de microondas), por lo que la señal reflejada por
en los 5 a 10 cm superiores de la columna del suelo, la reflectometría las superficies cercanas y registrada por la antena contiene
interferométrica GNSS (GNSS-IR) se ha convertido en una técnica de información sobre el entorno que rodea a la antena. Esta
detección remota interesante y complementaria debido a su ventajas información se puede obtener procesando la relación señal-
sobre las imágenes clásicas de satélite o avión. La resolución GNSS-IR es ruido (SNR) registrada en la antena como interferogramas.
mayor (una escala de aproximadamente 1000 m2alrededor de la La técnica GNSS-IR, desarrollada por primera vez por Larson
antena); Se puede utilizar para monitoreo continuo y es et al. (2008a,b), y basándose en el procedimiento detallado en
Larson et al. (2010), Larson y Nievinski (2013), Chew et al. (2014,
2015,2016), Vey et al. (2016), pequeño y cols. (2016), Wan et al. (
GPS Tool Box es una columna dedicada a resaltar los algoritmos y el 2015), Chen et al. (2016), Roussel et al. (2016), y Zhang et al. (
código fuente utilizados por los ingenieros y científicos de GPS. Si tiene
2017), se resume de la siguiente manera:
un programa o paquete de software interesante que le gustaría
compartir con nuestros lectores, páselo; envíenoslo por correo
electrónico a gpstoolbox@ngs.noaa.gov . Para comentar sobre 1. Las observaciones SNR de los satélites GNSS deberían seleccionarse
cualquiera de los códigos fuente discutidos aquí, o para descargar el para satélites con ángulos de elevación comprendidos entre 5 y 30
código fuente, visite nuestro sitio web en http://www.ngs.noaa.gov/gps-
grados. Este es uno de los principales inconvenientes de la técnica.
toolbox. Esta columna está editada por Stephen Hilla, National Geodetic
Survey, NOAA, Silver Spring, Maryland, y Mike Craymer, Geodetic Survey Puede usarse sólo para satélites de baja elevación ya que este
Division, Natural Resources Canada, Ottawa, Ontario, Canadá. rango es el único útil para obtener un interferograma válido para el
análisis de la humedad del suelo. Los datos de elevaciones de
satélite inferiores a 5 grados se descartan para evitar fuertes
* Ángel Martín
efectos de trayectos múltiples en los datos SNR.
aemartin@upvnet.upv.es
2. Las pistas de satélite ascendentes o descendentes se separan (o se
1
Departamento de Ingeniería Cartográfica, Geodesia y etiquetan para el posprocesamiento), ya que el patrón de interferograma
Fotogrametría, Universitat Politècnica de València, Camino
de Vera s/n, 46022 Valencia, España

Vol.:(0123456789)
 94 Página 2 de 7 Soluciones GPS (2020) 24:94

pueden diferir para las trayectorias ascendente y descendente del mismo En Nievinski y Larson (2014) se desarrolla un código fuente abierto
satélite. en Matlab/Octave que puede producir observaciones GPS simuladas de
3. Los datos SNR se convierten de unidades dB-Hz SNR, fase portadora y pseudodistancia que concuerdan con un modelo
observadas a una escala lineal en voltios mediante la de trayectorias múltiples para aplicaciones de posicionamiento y
expresión SNR.lineal= 10SNR/20(Vey et al.2016). reflectometría cercana a la superficie.
4. La señal reflejadaNRreflejadolineal
se aísla ajustando un En Roesler y Larson (2018), se presenta una herramienta de software
polinomio de segundo orden a la SNRlinealpara gratuita para traducir observaciones GPS (o GNSS) a un formato
eliminar la señal directa de satélite (Wan et al.2015; utilizable para la investigación de reflexiones (escrito en Fortran 77),
Chew et al.2016). mapear zonas de reflexión GNSS-IR alrededor de la antena (escrito en
5. Se calcula un periodograma de Lomb-Scargle a partir de la Matlab) y estimar frecuencias dominantes y altura de reflexión a partir
NRreflejado
lineal de cada pista de satélite para comprobar que sólo de datos GNSS (escritos en Matlab).
se observa una onda primaria clara. Se deben descartar las En este manuscrito, explicamos las ideas principales y las
pistas con múltiples picos o potencia media máxima baja. decisiones de implementación para escribir el código Python para
herramientas de software que implementan los primeros seis pasos
6. Las pistas seleccionadas se pueden modelar como: del procedimiento anterior (el último requiere datos in situ para
obtener una relación lineal precisa), y compartimos las software con
dóndeAyϕson la amplitud y la fase de la onda
la comunidad científica. El software desarrollado funciona en
primaria,λes la longitud de onda de la señal GNSS,
Python 2.7 y Python 3.
mies la elevación del satélite, yhes la altura del
reflector: la distancia vertical entre el centro de fase
de la antena GNSS y la superficie reflectante
Desarrollo de software
horizontal, que se supone que es la distancia entre la
antena y el suelo debido a la baja penetración de la
Las herramientas de software propuestas tienen dos módulos
señal en el suelo (Chew et al. 2014; Roussel et al.2016
principales. El primer módulo transforma los archivos de
; Zhang et al.2017).A yϕse estiman mediante un
observación y navegación en formato RINEX (Gurtner y Estey2015) a
algoritmo de mínimos cuadrados con valores
un archivo único que contiene la época de cada observación, la
iniciales de 1 paraAy 0 paraϕ.
identificación del satélite, la observación SNR y el acimut y elevación
( )
4 h calculados del satélite a partir del archivo de navegación RINEX.
NRreflejado
lineal
=Aporque pecadomi+� (1)
Solo se escribe una línea en el archivo de salida si la elevación del
satélite cae entre 5 y 30 grados. La entrada puede ser un único
7. El paso final es derivar la relación entre las variaciones de la humedad
archivo (un día de observación, por ejemplo) o varios (una semana o
del suelo (más específicamente, el contenido volumétrico de agua,
un mes de datos de observación continua separados por archivos
VWC m3/metro3) y variaciones en la fase estimada de la onda
diarios), pero la salida es siempre un único archivo.
primaria reflejada. Para ello, se necesitan observaciones in situ
La entrada del segundo módulo es la salida del primer
como conjunto de datos de referencia. Esas observaciones se
módulo; genera un archivo de salida por satélite que
pueden obtener a partir de sensores reflectómetros de contenido
contiene la altura de reflexión y los valores ajustados de
de agua convencionales (Vey et al.2016; Larson et al.2010) o de
fase y amplitud de la onda interferométrica para cada pista
muestras de datos de suelo (Martín et al.2020). Considerando una
de satélite individual. El software también genera una salida
relación lineal entre las variaciones de fase estimadas del GNSS-IR y
gráfica para cada pista que contiene la señal SNR directa, la
las variaciones de VWC de referencia, se necesitan varias semanas
señal SNR indirecta, el interferograma calculado y la onda
(o meses) de ambos tipos de datos para obtener una buena
ajustada a la señal SNR indirecta.
relación. La pendiente se ajusta utilizando las pistas satelitales para
las cuales las variaciones de fase presentaron una correlación lineal
Primer módulo: extracción de información de archivos de
más fuerte con las variaciones de humedad del suelo in situ. Por
observación y navegación de RINEX
ejemplo, en Zhang et al. (2017) esta correlación se establece en 0,9,
por lo que sólo se utilizan las trayectorias ascendentes de los
El software desarrollado funciona con el formato versión 3 de
satélites GPS 13, 21, 24 y 30 y las trayectorias descendentes de los
los archivos de observación y navegación RINEX. El software de
satélites GPS 05, 09, 10, 15 y 23. Finalmente, la pendiente que se
código abierto GFZRNX puede traducir datos de la versión 2 de
utilizará para todas las trayectorias válidas debe ser el valor de
RINEX a la versión 3 (Nischan2016). Aunque en el manuscrito
pendiente medio obtenido para las trayectorias satelitales
nos referimos a GNSS, la primera versión del software sólo
altamente correlacionadas.
funciona para la constelación de satélites GPS.

13
Soluciones GPS (2020) 24:94 Página 3 de 7 94

El usuario debe decidir la frecuencia que utilizará antes de ejecutar el Las coordenadas de la antena y el satélite, el acimut y la
software; Esto se puede hacer abriendo un archivo de observación elevación se calculan y almacenan en el archivo de salida. El
RINEX e inspeccionando la información del encabezado. Una vez módulo que calcula el azimut y la elevación requiere dos
seleccionada la frecuencia, el usuario deberá determinar las módulos adicionales, uno para calcular el tiempo juliano y
observaciones relacionadas con esta frecuencia y almacenarlas en el GPS y otro para calcular las coordenadas geodésicas
archivo de observaciones RINEX que será utilizado en el proceso. Esas esféricas de las coordenadas geocéntricas y calcular el
observaciones son la SNR y el pseudorango o código de observación acimut y la elevación a partir de las coordenadas
relacionado con la frecuencia seleccionada. geodésicas de la antena y el satélite.
El algoritmo lee los archivos de observación y navegación
almacenados en una carpeta (\data\input\ por defecto). Los
nombres y extensiones de dichos archivos deben seguir el formato
estándar RINEX: El nombre se compone de cuatro caracteres para el
Segundo módulo: cálculo de la altura del reflector GPS-IR y
nombre de la estación, tres para el día del año y un carácter para
ajuste de la onda a la SNR observada
identificar la sesión. La extensión se compone de dos caracteres
correspondientes a los dos últimos dígitos del año y una "o" para
El segundo módulo de software utiliza la salida del primer
los archivos de observación o una "n" para el archivo de navegación
módulo como entrada. El programa identifica todos los tracks
GPS. El algoritmo ordena los archivos de entrada en secuencia
válidos por satélite y resolverá los pasos 2 a 6 descritos en la
cronológica y los abre en orden.
sección de introducción.
El algoritmo lee cada observación de satélite GPS para cada
Primero, el software crea una estructura de carpetas de
época de cada archivo de observación y escribe en el archivo de
salida: una carpeta para cada satélite identificado por su
salida (ubicado en la carpeta data\output1\ de forma
identificador numérico, creando 32 carpetas numeradas del 1 al
predeterminada) los siguientes datos solo de satélites con
32. Si estas carpetas ya existen, el programa fallará, ya que se
elevaciones entre 5 y 30 grados:
crean durante la ejecución; esto evita que los resultados
anteriores se sobrescriban o generen más archivos dentro de
1. Un identificador numérico relacionado con las épocas del archivo
las carpetas con cada ejecución sucesiva.
observado. Este número de identificador entero comienza con
Antes de ejecutar el software, el usuario debe configurar algunos
cero para todas las observaciones de la primera época del
parámetros de entrada internos:
primer archivo de observación, aumenta una por una y termina
con la última época del último archivo de observación. Este
1. El intervalo de tiempo entre observaciones: Se puede obtener del
número numérico se puede utilizar como identificador en el
encabezado del archivo de observaciones.
segundo módulo de software. El software escribe tantas líneas,
2. Las elevaciones mínima y máxima del satélite: Estos parámetros
con el mismo identificador numérico, como observaciones de
se vuelven a configurar, por defecto, en 5 y 30 grados.
diferentes satélites en la misma época.
3. El azimut mínimo y máximo de las trayectorias del satélite a
2. El tiempo relacionado con las líneas anteriores. El año, mes, día y
considerar: Esto permite al usuario enfocar un área
un número flotante que contiene la hora (calculado con la
determinada alrededor de la antena. Estos parámetros están
información de horas enteras, minutos enteros y segundos
configurados por defecto en 0 y 360 grados.
flotantes del archivo de observación RINEX).
4. La altura de la antena: medida desde el suelo hasta la marca de
3. El identificador numérico GPS satelital de dos dígitos.
referencia de la antena en metros.
4. La SNR observada.
5. Longitud de la frecuencia de trabajo: 0,1904 m para L1,
5. Basado en la época de observación, la pseudodistancia o código
0,2443 m para L2.
de observación al satélite, las coordenadas de la estación
6. Un rango mínimo de ángulo de elevación que debe cubrir el satélite
ubicadas en el encabezado del archivo de observación y el
para considerar el track como válido: establecido por defecto en 10
archivo de navegación, el algoritmo calcula el acimut y la
grados.
elevación del satélite desde la antena. Este cálculo debe
7. Un número mínimo de épocas sin observaciones para asumir que las
determinar el tiempo de emisión y las coordenadas del satélite
observaciones pertenecen a la misma ruta satelital o a una nueva: en
en el sistema de referencia Earth-Fixed-Earth-Center (ECEF).
ciertos casos, parte de la información de la época se pierde en el archivo
Para el tiempo de emisión se utiliza un “algoritmo basado en
de observación debido a interrupciones de la señal o un mal
pseudodistancia” iterativo (Sanz et al. 2013). El algoritmo
funcionamiento, por lo que el software identifica las observaciones
descrito en Leick et al. (2015) en la pág. 240 se utiliza para el
satelitales. que están separados temporalmente por menos del número
cálculo de las coordenadas del satélite ECEF. Finalmente, se
de épocas establecidas en este parámetro como pertenecientes a la
tiene en cuenta la rotación de la Tierra durante el viaje de la
misma pista. El valor depende del intervalo de tiempo entre épocas; se
señal para obtener las coordenadas finales del satélite en el
puede establecer en 1 si el intervalo de tiempo
tiempo del receptor. Basado en el ECEF

13
 94 Página 4 de 7 Soluciones GPS (2020) 24:94

entre períodos de observación es de 30 s, 7 si el intervalo es de 5 s o 6. número de identificación del satélite

35 para 1 s. 7. acimut inicial de la trayectoria del satélite

8. acimut final de la trayectoria del satélite

El software selecciona cada pista y la etiqueta como pista de 9. elevación inicial de la trayectoria del satélite

satélite ascendente o poniente; luego realiza los procesos 10. elevación final de la trayectoria del satélite

descritos en los siguientes párrafos. 11. ajustadoAvalor

Los pasos 3, 4, 5 y 6 descritos en la sección de introducción se 12. desviación estándar del ajustadoAvalor
realizan para cada pista; finalmente,AyϕLos parámetros y sus 13. ajustadoϕvalor (en grados)
desviaciones estándar se estiman mediante el algoritmo de 14. desviación estándar del ajustadoϕvalor (en grados)
mínimos cuadrados. 15. altura del reflector calculada
Una de las mayores ventajas de utilizar Python son las 16. Campo de texto de satélite en caída/satélite en ascenso

numerosas bibliotecas científicas gratuitas disponibles para un fácil

acceso. El proceso más complejo, en este caso, es el cálculo del También se crea un segundo archivo de salida, que contiene
periodograma de Lomb-Scargle: paso 5 en la sección de la misma información pero solo para las pistas consideradas
introducción. Para ello se utiliza la función LombScargle de la "buenas o válidas". Para todas las pistas válidas e inválidas, se
biblioteca Astropy (https://docs.astropy.org/en/stable/timeseries / genera un gráfico que contiene cuatro cifras: la señal completa
lombscargle.html). La entrada para calcular el periodograma de (SNR), la señal reflejada (NRlineal
reflejado), el periodograma de Lomb-

Lomb-Scargle para cada pista es el seno del ángulo de elevación del Scargle para la señal reflejada y la señal reflejada con la onda
satélite en el eje X, y elNRreflejadoen el eje Y. Con esta configuración,
lineal
ajustada. El nombre de la figura incluye el identificador
el resultado convierte la frecuencia en la altura del reflector de la numérico, año, mes, día, identificador de satélite, el azimut
antena en metros en el eje X de salida. Sin embargo, debido al uso inicial y el azimut final para permitir al usuario identificar de
del seno del ángulo de elevación del satélite, el espaciado de la forma clara y confiable la figura con la línea correspondiente en
cuadrícula en el eje X es irregular. Para determinar el espaciado de los archivos de salida.
cuadrícula apropiado a usar, la biblioteca introduce una opción a
través de palabras clave pasadas alapagado automático()método.
De forma predeterminada (después de algunas pruebas Ejemplo de GPS-IR
experimentales), la frecuencia más alta se fija en dos veces la
frecuencia promedio de Nyquist:factor_nyquist=2en el código. Los datos que proporcionamos con el software forman parte de un
Además, la altura máxima del reflector permitida en el experimento realizado en las instalaciones del Centro de
periodograma de salida se fija en 2,5 m, pero este parámetro se Experiencias Cajamar, Paiporta, Valencia, España, Fig.1, Martín et al.
puede cambiar si la elevación de la antena es mayor. Se puede
encontrar una discusión teórica sobre la extracción de frecuencia
de GNSS-IR en Roesler y Larson (2018).
El paso más importante de este módulo es establecer las
condiciones para seleccionar pistas “buenas”. En este caso, en
base a diferentes experimentos, un track se considera válido si
el track del satélite contiene más de 30 min de observación y
más que el rango de elevación mínimo previamente establecido
que debe cubrir el satélite, la potencia de la frecuencia
dominante es 6 veces mayor que el ruido de fondo de los
medios, y el ajuste de la onda teórica alNRreflejado lineal
La señal presenta un vector residual con una media inferior a 1,3 V/
Voltios y una desviación estándar inferior a 25 V/Voltios.
Se genera un archivo de salida, en el que cada línea
corresponde a una pista y contiene las siguientes columnas:

1. identificador numérico relacionado con las épocas del archivo

observado
2. campo de texto sí/no relacionado con las condiciones anteriores
(sí indica pistas “buenas o válidas”)
3. época del tiempo de observación
4. mes del tiempo de observación Figura 1Antena GNSS de calidad geodésica ubicada en la zona del
5. día del tiempo de observación experimento

13
Soluciones GPS (2020) 24:94 Página 5 de 7 94

Figura 2Datos SNR observados en dB-Hz Figura 4 Periodograma de Lomb-Scargle para la señal reflejada SNR

Fig. 3Datos SNR reflejados o indirectos en voltios

figura 5SNR señal reflejada con la onda ajustada en Voltios

(2020). Proporcionamos 7 días (de un total de 66 en el

experimento) de observaciones GPS con un receptor geodésico
Con los criterios establecidos para seleccionar las trayectorias
GNSS (Trimble R10, del Departamento de Ingeniería
correctas, la matriz de confusión obtenida para todos los satélites
Cartográfica, Geodesia y Fotogrametría de la Universitat
durante los 7 días de observación presenta los siguientes valores: el 45%
Politècnica de València). Se proporcionan observaciones RINEX
de las trayectorias durante los siete días son válidas y clasificadas como
(con una frecuencia de muestreo de 5 segundos) y archivos de
tales, el 43% no son válidas y se clasifican como Así, el 9% son pistas no
navegación. Para el ejemplo se utiliza la frecuencia L1.
válidas clasificadas como válidas, y sólo el 3% de las pistas son válidas
El experimento se realizó del 3 de diciembre de 2018 al 6 de
pero clasificadas como no válidas.
febrero de 2019. Los archivos de observaciones incluidos con el
software son archivos de datos diarios del 14 al 20 de enero de
2020 y la altura de la antena al suelo fue de 1,8 m. Conclusiones y comentarios finales.
No hubo lluvias durante la observación. Además del
software, se proporciona el archivo de salida del primer Pretendemos compartir el software con la comunidad científica
módulo de software, junto con las figuras y los dos archivos de para presentar a nuevos usuarios la técnica GNSS-IR. Esta técnica
salida del satélite número 15 del segundo módulo, para que el puede usarse no sólo para monitorear la humedad del suelo sino
usuario pueda validar el funcionamiento del software. Cifras2, 3 también, por ejemplo, para monitorear el contenido de agua de la
,4,5mostrar la pista con identificador época 88,518, enero vegetación (Wan et al.2015), mediciones de la profundidad de la
19 de octubre de 2019, azimut inicial de 281 grados y azimut nieve (Larson et al.2009), o mediciones de mareas (Roussel et al.
final de 300,1 grados. En concreto, la Fig.2son los datos SNR en 2015), lo que ilustra que los usuarios pueden adaptar fácilmente el
dB-Hz, Fig.3son los datos SNR indirectos en voltios, Fig.4es el software para otros fines.
periodograma de Lomb-Scargle para la señal reflejada SNR, y la El software también se puede ampliar para que funcione con
Fig.5Representa la señal reflejada SNR con la onda ajustada en observaciones y frecuencias de los satélites Galileo, GLONASS y Beidou
voltios. GNSS, modificando únicamente el primer módulo de software.

13
 94 Página 6 de 7 Soluciones GPS (2020) 24:94

Disponibilidad de datos y material. Sensores como solución para el monitoreo de la humedad del suelo. Ciencia
del sistema hidrol terrestre.https://doi.org/10.5194/hess-24-3573-2020
Nievinski GG, Larson KM (2014) Un GPS multitrayectoria de código abierto
El software está disponible en el sitio web de GPS Toolbox en Simulador en Matlab/Octave. Solución GPS 18:473–481.https://doi.
https://geodesy.noaa.gov/gps-toolbox/. org/10.1007/s10291-014-0370-z
Nischan T (2016) GFZRNX: conversión de datos RINEX GNSS y
AgradecimientosLos autores quieren agradecer al personal del Centro caja de herramientas de manipulación (Versión 1.05). Servicio de datos GFZ.https://

de Experiencias Cajamar, y en especial a Carlos Baixauli, su apoyo y doi. org/10.5880/GFZ.1.1.2016.002

colaboración en el experimento Paiporta. Los autores también quieren Roesler C, Larson KM (2018) Herramientas de software para interferómetro GNSS

agradecer a Alfred Leick y Steve Hilla por sus valiosos comentarios y Reflectometría rica (GNSS-IR). Solución GPS.https://doi.org/10.1007/
sugerencias. s10291-018-0744-8
Roussel N, Ramilien G, Frappart F, Darrozes J, Gay A, Biancale R,
Striebig N, Hanquiez V, Bertin X, Allain A (2015) Monitoreo y
estimación del nivel del mar utilizando un único receptor
geodésico. Entorno de sensores remotos 171:261–277.https://
Referencias doi.org/10.1016/j. rse.2015.10.011
Roussel N, Frappart F, Ramillien G, Darroes J, Baup F, Lestarquit L,
Ha MC (2016) Detección de variaciones de humedad del suelo utilizando
Chen Q, Won D, Akos DM, Small EE (2016) Vegetación mediante GPS
datos GPS y GLONASS SNR para ángulos de elevación que van desde 2º a
Reflectometría interferométrica: resultados experimentales con
70º. IEEE J Selec Top Appl Earth Obs Rem Sens 9(10):4781– 4794.https://
una antena polarizada horizontal. IEEE J Select Top Appl Earth Obs
doi.org/10.1109/JSTARS.2016.2537847
Rem Sens 9(10):4771–4780.https://doi.org/10.1109/JSTAR
Sanz J, Juan JM, Hernández-Pajares M (2013) Procesamiento de datos GNSS.
S.2016.2565687
Tomo I: fundamentos y algoritmos. Comunicaciones de la Agencia
Chew CC, Small EE, Larson KM, Zavorotny VU (2014) Efectos de
Espacial Europea, 223 págs.
Humedad del suelo cerca de la superficie en datos GPS SNR: algoritmo
Small EE, Larson KM, Chew CC, Dong J, Ochsner TE (2016) Validación
de desarrollo y recuperación para la humedad del suelo. IEEE T Geosci
ción de recuperaciones de humedad del suelo GPS-IR: comparación de
Rem Sens 52(1):537–543.https://doi.org/10.1109/TGRS.2013.2242332
diferentes algoritmos para eliminar los efectos de la vegetación. IEEE J
Chew CC, Small EE, Larson KM, Zavorotny UZ (2015) Vegetación
Selec Top Appl Earth Obs Rem Sens 9(10):4759–4770.https://doi.org/
detección mediante reflectometría interferométrica GPS: efectos
10.1109/JSTARS.2015.2504527
teóricos de los parámetros del dosel en los datos de la relación señal-
Vey S, Güntner A, Wickert J, Blume T, Ramatschi M (2016) Long-
ruido. IEEE T Geosci Rem Sens 53(5):2755–2764.https://doi.org/10.1109/
Dinámica de la humedad del suelo a término derivada de la reflectometría
TGRS.2014.2364513
interferométrica GNSS: un estudio de caso para Sutherland, Sudáfrica.
Chew CC, Small EE, Larson KM (2016) Un algoritmo para la humedad del suelo
Solución GPS 20:641–654.https://doi.org/10.1007/s10291-015-0474-0 Wan W,
estimación mediante reflectometría interferométrica GPS para suelos
Larson KM, Small EE, Chew CC, Braun JJ (2015) Usando
desnudos y con vegetación. Solución GPS 20(3):525–537.https://doi.org/
Receptores GPS geodésicos para medir el contenido de agua de la vegetación.
10.1007/s10291-015-0462-4
Solución GPS 19:237–248.https://doi.org/10.1007/s10291-014-0383-7 Zhang S,
Gurtner W, Estey L (2015) RINEX: el intercambio independiente del receptor
Roussel N, Boniface K, Ha MC, Frappart F, Darrozes J, Baup
versión de formato 3.03.ftp://igs.org/pub/data/format/rinex303.pdf
F, Calvet JC (2017) Uso de datos SNR GNSS reflejados para
Larson KM, Small EE, Gutmann ED, Bilich AL, Axelrad A, Braun
recuperar la humedad del suelo o la altura de la vegetación de un
JJ (2008a) Uso de rutas múltiples de GPS para medir las fluctuaciones de
cultivo de trigo. Hydrol Earth Syst Sci 21:4767–4784.https://doi.org/
humedad del suelo: resultados iniciales. Solución GPS 12(3):173–177.https://
10.5194/hess-21-4767-2017
doi. org/10.1007/s10291-007-0076-6
Larson KM, Small EE, Gutmann ED, Bilich AL, Braun JJ, Zavorotny
Nota del editorSpringer Nature se mantiene neutral con respecto a reclamos
VU (2008b) Uso de receptores GPS como red de humedad del suelo para
jurisdiccionales en mapas publicados y afiliaciones institucionales.
estudios del ciclo del agua. Geophys Res Lett 35:L24405.https://doi. org/
10.1029/2008GL036013
Larson KM, Gutmann E, Zavorotny VU, Braun J, Williams M, Niev-
inski FG (2009) ¿Podemos medir la profundidad de la nieve con
Ángel Martínrecibió su doctorado.
receptores GPS? Geophys Res Lett 36:L17502.https://doi.org/
Licenciado en Geodesia y
10.1029/2009G L039430
Cartografía por la Universitat
Larson KM, Braun JJ, Small EE, Zavorotny VU (2010) GPS mul-
Politècnica de València (UPV) en
tipath y su relación con el contenido de humedad del suelo cerca de la
2001. Es científico titular y profesor
superficie. IEEE J Selec Top Appl Earth Obs Rem Sens 3(1):91–99.https://
titular del Departamento de
doi. org/10.1109/JSTARS.2009.2033612
Ingeniería Cartográfica, Geodesia y
Larson KM, Nievinski FG (2013) Detección de nieve por GPS: resultados de la
Fotogrametría de la UPV. Sus
Observatorio del límite de placas EarthScope. Solución GPS 17(1):41–52.
principales áreas de investigación
https://doi.org/10.1007/s10291-012-0259-7
son GNSS, Geodesia Física y Big
Leick A, Rapoport L, Tatarnikov D (2015) Levantamiento por satélite GPS, 4to
Data y Machine Learning aplicados a
ed. Wiley, Hoboken, página 840
las Geociencias.
Martín A, Ibañez S, Baixauli C, Blanc S, Anquela AB (2020) Multi-
Reflectometría interferométrica de constelación con mercado masivo.

13
Soluciones GPS (2020) 24:94 Página 7 de 7 94

Raquel LujánRecibió su doble Ana Belén Anquelarecibió su

máster por la Universitat doctorado. Licenciada en
Politècnica de València (UPV) y la Geodesia y Cartografía por la
Universidad de Ciencias Aplicadas Universitat Politècnica de València
de Karlsruhe en 2019. Es Ph.D. (UPV) en 2001, donde es científica
Estudiante del Departamento de titular y profesora profesora del
Ingeniería Cartográfica, Geodesia Departamento de Ingeniería
y Fotogrametría de la UPV donde Cartográfica, Geodesia y
trabaja en el campo de los GNSS. Fotogrametría. Sus principales
áreas de investigación son GNSS y
monitoreo geodésico. También es
directora de la escuela de
geomática de la UPV.

13
Reproducido con permiso del propietario de los derechos de autor. Reproducción adicional

prohibido sin permiso.