Diferencias entre TLE y RINEX gps

Tengo archivo de navegación TLE y RINEX para satélite GPS PRN 20, del 01.04.17r. 19:56:56 y 20:00:00 respectivamente. Consulte también el formato de intercambio independiente del receptor, versión 3.03 .

Parece que el valor de la longitud del nodo ascendente (Ω0) difiere para ambos, mientras que todos los demás parámetros son relativamente similares.

Archivo RINEX:

M0 = -0.41907076212790001
e  =  0.0044262305600569997
ω  =  1.5754052384999999
i  =  0.92697366533970005
Ω0 = -0.4446281969834

Archivo TLE:

M0 =  5.8389169820721856
e  =  0.0042077
ω  =  1.5770882387483363
i  =  0.92682044735729674
Ω0 =  2.7589292150900362

¿Cuál puede ser la razón de esto?

Si bien los valores de M0están cerca del módulo 2π, los valores de Ω0 parecen diferir solo en π, o alrededor de 180°.

Comprobando también otros conjuntos de datos, el resultado es similar, aunque la diferencia entre los valores difiere.

TLE:

GPS BIIR-4  (PRN 20)    
1 26360U 00025A   17091.83120691 -.00000002  00000-0  00000-0 0  9995
2 26360  53.1029 158.0750 0042077  90.3605 334.5453  2.00565359123815

RINEX:

20 17  4  1 20  0  0.0 4.665344022214D-04 1.591615728103D-12 0.000000000000D+00
5.600000000000D+01-1.231250000000D+01 5.538802141913D-09-4.190707621279D-01
-5.606561899185D-07 4.426230560057D-03 6.150454282761D-06 5.153646507263D+03
 5.904000000000D+05-6.705522537231D-08-4.446281969834D-01 3.725290298462D-08
 9.269736653397D-01 2.442500000000D+02 1.575405238500D+00-8.669646839876D-09
 2.714398779890D-11 0.000000000000D+00 1.942000000000D+03 0.000000000000D+00
 2.000000000000D+00 0.000000000000D+00-8.381903171539D-09 5.600000000000D+01
 5.899500000000D+05 4.000000000000D+00
Como sugerencia, ¿has comprobado que PRN20 de TLE y RINEX son el mismo satélite?. el desplazamiento de METRO 0 y Ω 0 que está informando son típicos entre diferentes satélites de una constelación. Además, su atajo de zapatos TLE y RINEX PRN20 mientras comenzaba la pregunta sobre PRN1 (¿un error tipográfico tal vez?)
¿De dónde sacaste el archivo RINEX?
El archivo rinex es de garner.ucsd.edu , TLE de celestrak.com/NORAD/elements/gps-ops.txt . @Julio, ¿cómo puedo verificar si PRN20 es el mismo satélite para ambos archivos? Sé que los números de satélite son diferentes para los archivos glonass rinex y tle, pero pensé que para el GPS es lo mismo. Sí, PRN1 es un error tipográfico.
@MartynaKasyr, si fuera usted, buscaría otro satélite en su TLE que coincida con los valores de los elementos orbitales actuales de su archivo RINEX para confirmar (o no) la discrepancia de los números PRN entre TLE y RINEX. También hay una pregunta abierta sobre este tema aquí space.stackexchange.com/questions/20887/…

Respuestas (2)

Están en diferentes marcos de referencia. El conjunto de elementos TLE se define en relación con el marco de referencia True Equator, Mean Equinox ( TEME ), y los elementos en el mensaje RINEX Nav se definen en relación con Earth Centered, Earth Fixed (ECEF). Dado que ambos marcos usan el ecuador verdadero, otros elementos (como la inclinación) no varían tanto como la ubicación del nodo.

Están en diferentes marcos de referencia, pero hay más que eso. Incluso una vez que haya convertido TEME a ECEF, seguirá habiendo diferencias, porque las cosas que tienen el mismo nombre en realidad tienen significados muy diferentes. RINEX es solo un formato de archivo; la carga recae sobre los usuarios para comprender que los mensajes de navegación GPS y los TLE utilizan los mismos nombres de campo para diferentes propósitos.

La especificación IS-GPS-200 (Revisión M, con fecha de mayo de 2021), sección 20.3.3.4, "Subtramas 2 y 3", de NAVSTAR GPS Space Segment/Navigation User Interface Specification , proporciona 12 páginas de instrucciones sobre cómo se requiere que el segmento de usuario realice cálculos de efemérides, que no se alinean del todo con la forma en que alguien más usa esos términos. En particular, tenga en cuenta estas declaraciones (énfasis mío):

Los parámetros de efemérides describen la órbita durante los intervalos de ajuste de la curva descritos en la sección 20.3.4. La Tabla 20-II da la definición de los parámetros orbitales usando la terminología típica de los parámetros orbitales Keplerianos; sin embargo, cabe señalar que los valores de los parámetros transmitidos son tales que proporcionan el mejor ajuste de trayectoria en coordenadas centradas en la Tierra, fijas en la Tierra (ECEF) para cada intervalo de ajuste específico. El usuario no interpretará los valores de las coordenadas intermedias como pertenecientes a ningún sistema de coordenadas convencional. El usuario deberá calcular las coordenadas ECEF de posición para el centro de fase de las antenas de los SV utilizando una variación de las ecuaciones que se muestran en la Tabla 20-IV. Los parámetros de las subtramas 2 y 3 tienen apariencia kepleriana; los valores de estos parámetros, sin embargo, se producen... a través de un ajuste de curva de mínimos cuadrados de las efemérides propagadas del centro de fase de las antenas de los SV

Los intervalos de ajuste son de tres o cuatro horas. Los números en los mensajes de navegación GPS brindan los valores que mejor se ajustan a lo largo de todo ese intervalo, por lo que no son lo mismo que los elementos osculadores de las órbitas GPS en una época específica. Tienes que leer los documentos para saber cómo manejarlos, y tienes que usar las ecuaciones que te indiquen (incluyendo, por ejemplo, una forma específica de iterar la ecuación de Kepler hacia una solución), o la respuesta que obtengas no tendrá la mismas estadísticas de error que estas tablas existen para permitirle lograr.

Estas advertencias son aún más ciertas acerca de los conjuntos de elementos de dos líneas. Esos números están aún más lejos de oscular. Se llaman elementos medios . Se podría argumentar que la palabra debería aplicarse a la forma en que lo hace el GPS, pero esa no es la forma en que se usa la palabra en la mecánica orbital. "Media" generalmente solo significa "promedio", pero hay muchas maneras de promediar algo, y el uso técnico del término en astrodinámica es más específico y más general de lo que cabría esperar. En palabras de HG Walter (1967), "por elementos medios entendemos los elementos osculantes de los que se han sustraído perturbaciones del potencial terrestre de periodo corto y de periodo largo". Para convertir de elementos medios a elementos osculadores, debe saber exactamente qué perturbaciones se eligieron para la resta, y debe calcular esas correcciones y agregarlas nuevamente, o obtendrá una respuesta incorrecta.

Los elementos de la órbita en los TLE se especifican en "forma Kozai", lo que significa usar la lista exacta de perturbaciones elegidas para su eliminación en Yoshihide Kozai, The motion of a close earth satellite , Astronomical Journal 64 (367-377) 1959. Por ejemplo, el valor colocarse en el campo denominado "semieje mayor" no es lo habitual a (que el periódico vinculado llama a 0 ), sino más bien el "medio" a ¯ , lo que equivale a 0 veces el factor

1 A 2 pags 2 ( 1 3 2 pecado 2 i ) ( 1 mi 2 )

dónde pags = a ( 1 mi 2 ) y A 2 es un coeficiente de expansión de la gravedad que es igual a 3 2 j 2 R 2 , dónde j 2 es lo familiar y R es el radio ecuatorial de la tierra.

Esta es, con mucho, la fórmula más simple de este tipo. Las expresiones equivalentes para el argumento del perigeo ( ω ¯ ) y ascensión recta del nodo ascendente ( Ω ¯ ) ocupan varias páginas. Tenga en cuenta que la conversión de Kozai A 2 al ahora común j 2 Obtuve de Brouwer (1959) , que describe un conjunto diferente y mucho más complicado de elementos medios, y tiene una tabla que muestra cómo traducir entre nueve formas diferentes de coeficientes de expansión de la gravedad. El artículo de Brouwer termina con la declaración,

Ahora me arrepiento de haber presentado k 2 , k 4 en mi trabajo en 1946. La razón principal fue que dan una forma particularmente simple para la expresión del potencial en el plano ecuatorial. Si hubiera podido prever el aumento del interés por el tema y la confusión a la que estaba contribuyendo, habría elegido... la forma alternativa que utilizó Vinti (1959) . Tengo la intención de volver a este formulario y recomendarlo a otros autores.

La súplica de Brouwer parece haber funcionado, ya que la astrodinámica parece haberse asentado firmemente en la forma que Vinti le atribuye a R. H. Merson y DG King-Hele , Uso de satélites artificiales para explorar el campo gravitacional de la Tierra: resultados del Sputnik 2 (1957β), Nature 182 (640) 1958, el familiar

V ( r , θ ) = GRAMO METRO r [ 1 pags = 2 j pags ( R r ) pags PAGS pags ( pecado θ ) ]