Actualmente estoy tratando de determinar la posición del sol y la luna a partir de jpl ephemeris DE200, que se refiere al ecuador dinámico y al equinoccio de 2000 y usa el tiempo tdb. Estoy usando el sistema j2000 ECI. Las cosas que me preocupan son las siguientes:
1.¿Se da la posición de la Luna con respecto a la Tierra o con respecto al baricentro? Además, cuando encontramos la posición del sol usando la posición del baricentro tierra-luna y la coordenada lunar, ¿será la posición del sol con respecto a la Tierra?
2. Si estoy usando la hora UTC, ¿cómo obtener exactamente el TDB o TAi, ya que está oscilando? (Después de eso, quiero obtener el JD para TDB)
3.Si estoy usando el sistema de coordenadas J2000 ECI, ¿debo tener en cuenta el efecto de la nutación y la precesión? ¿Y como hacer eso?
La posición de la Luna se da con respecto al baricentro de la Tierra y la Luna, cuyo NAIF ID es 3 . Tenga en cuenta que el baricentro del sistema Tierra-Luna está muy cerca del centro de la Tierra (de lo contrario, a los humanos no les habría llevado tanto tiempo averiguar de dónde venían las mareas). Debería poder enumerar todos los cuerpos de un archivo BSP determinado mirando al principio del archivo. Si está ejecutando Linux o MacOS, abra una terminal e inspeccione ese archivo con od -vc de438.bsp | less
. Tus ojos tardarán unos segundos en darse cuenta de lo que está pasando, pero luego es algo legible para los humanos.
Los tiempos en estos núcleos se dan en "tiempo de efemérides". Con toda honestidad, olvidé contra qué está vinculado el tiempo de las efemérides. Este tutorial de SPK menciona que el "archivo de segundo bisiesto" .LSK
permite una conversión entre ET y UTC, por lo que sospecho que el tiempo de las efemérides podría ser UTC sin los segundos intercalares. Si solo necesitas una lista de segundos bisiestos (que es lo que obtienes en los archivos naif000xx.lsk), te recomendaría usar el IETF life from here , ya que son los que anuncian los segundos bisiestos con unos meses de antelación. Si por suerte está utilizando un lenguaje de programación llamado Rust, es posible que desee consultar la biblioteca llamada hifitimeque ya incluye la fecha juliana modificada y la gestión de la fecha UTC con segundos bisiestos (divulgación completa, yo lo escribí). (Dado que es Rust, debería ser razonablemente fácil usar esta biblioteca desde C o C++, y un poco menos fácil de usar en Python).
ECI significa marco inercial centrado en la tierra. Como marco inercial (es decir, no tiene aceleración), no se aplican ni la nutación ni la precesión. Sin embargo, el marco ECEF (Earth Centered Earth Fixed) necesita manejar el tiempo solar aparente de Greenwich, la nutación, la precesión y el movimiento polar. De acuerdo con el marco IAU2000, el marco "ECI" pasó a llamarse CRS para "Sistema de referencia celestial" y el marco ECEF ahora es TRS para "sistema de referencia terrestre". El marco ECEF/TRS solo es necesario si necesita armónicos esféricos o algunos parámetros geodésicos específicos (la latitud, la longitud y la altura geodésica son con respecto al elipsoide de referencia y, por lo tanto, solo requieren conocer el coeficiente de aplanamiento del cuerpo). Por ejemplo, si desea calcular la visibilidad de un vehículo en el marco ECI desde el punto de vista de un observador en el marco ECEF, deberá incluir todos los parámetros de dicho marco. Finalmente, tenga en cuenta que el marco IAU2000 permite un cálculo decentemente fácil del tiempo sideral aparente de Greenwich, la precesión y la nutación, pero el movimiento polar esno predecible y como tal es publicado regularmente por BIPM (Bureau International des Poids et Mesures). Dicho esto, el movimiento polar es increíblemente pequeño (creo que el eje cambia 0,002 grados por siglo), por lo que es probable que no necesite tenerlo en cuenta para la mayoría de los cálculos a corto y mediano plazo. Sobre el tema de la transformación del marco de referencia IAU2000, recomiendo leer la Nota técnica n.º 29 del IERS de acceso abierto de McCarthy y Capitaine (cuyo título completo es: Consecuencias prácticas de la Resolución B1.6 "Modelo de nutación de precesión IAU2000", Resolución B1 .7 "Definición de Polo Celeste Intermedio", y Resolución B1.8 "Definición y Uso de Origen de Efemérides Celeste y Terrestre").
Como mencionó uhoh en el comentario anterior, DE200 es unos años más antiguo, y lo alentaría a usar un kernel de efemérides planetaria más reciente. Por ejemplo, DE438 se publicó el lunes pasado (4 de junio de 2018), cf. ftp://ssd.jpl.nasa.gov/pub/eph/planets/ascii/de438/README.txt . Además, si solo necesita fechas de enero de 1900 a diciembre de 2100, probablemente debería usar de438s.bsp (cf. ftp://ssd.jpl.nasa.gov/pub/eph/planets/bsp/ ).
UH oh
Alexandru Lapusneanu
UH oh
usuario7073