Transforme la rotación de asteroides en coordenadas eclípticas heliocéntricas

Estoy trabajando con datos de la base de datos DAMIT de modelos de formas de asteroides. Los estoy agregando a una visualización en la que el sol está en [0, 0, 0] y los ejes X, Y constituyen el plano de la eclíptica del sistema solar.

Cada modelo de asteroide viene con algunos atributos que definen su orientación y giro:

  • λ (longitud eclíptica),
  • β (latitud eclíptica)
  • P (período de rotación sideral)
  • φ0 (ángulo de rotación inicial)
  • JD0 (fecha inicial).

He aplicado las fórmulas matriciales sugeridas por la gente de DAMIT:

rotación eclíptica

en el cual

definiciones de matrices[2]

Aquí es donde me confundo. Tengo entendido que mi visualización utiliza un sistema de coordenadas de la eclíptica heliocéntrica.

Configuré r_astla ubicación XYZ de un vértice en mi visualización. No estoy seguro de que esto sea correcto.

Luego calculo r_ecla partir de la ecuación anterior. Pero creo que las variables λ y β son ángulos en un sistema de coordenadas eclípticas geocéntricas y necesito hacer una transformación en coordenadas eclípticas heliocéntricas.

¿Mi enfoque es correcto y, de ser así, cómo puedo convertir la rotación del asteroide en el sistema de coordenadas de la eclíptica heliocéntrica de mi visualización?

λ y β son las coordenadas del polo del asteroide en un marco eclíptico heliocéntrico, así que creo que está bien si esto también es lo que usa su visualización
Gracias @astrosnapper. ¿Cómo puedes saber que λ y β son heliocéntricos en lugar de geocéntricos en este caso? De esta tabla de notación de coordenadas eclípticas en Wikipedia, supuse que eran geocéntricas: en.wikipedia.org/wiki/…
Técnicamente, no especifica el origen (y debería), pero dado que el asteroide está en órbita alrededor del Sol y no de la Tierra, lo más probable es que sea heliocéntrico o geocéntrico.
Gracias @astrosnapper. Un punto más de aclaración. ¿Es correcto establecer r_astlas coordenadas rectangulares de la eclíptica heliocéntrica del asteroide? El sitio web lo describe r_astcomo un vector en el "marco de coordenadas co-rotante del asteroide", pero no estoy seguro de lo que eso significa aquí.

Respuestas (1)

r_astes un vector que va desde el centro del asteroide a cada vértice del modelo de forma como si hubieras clavado un juego de varillas en el centro de una papa y apuntando en el eje largo, al polo de arriba y 90 grados a estos para dar direcciones X,Y,Z.

Las matrices anteriores que das transforman estos r_astvectores para las coordenadas de la superficie del asteroide en el plano de la eclíptica en función del tiempo a medida que el asteroide gira, dándote r_ecl. Dice en el artículo de Kaasalainen (creo) que r_ecltiene un origen desplazado hacia el centro del asteroide. Si desea una visualización de cómo se ve el asteroide en relación con la Tierra o el Sol, etc., entonces deberá calcular el vector Sol->asteroide y el vector Sol->Tierra usando, por ejemplo, JPL HORIZONTES y usando el tipo de efemérides de tabla vectorial o un contenedor programático para esto (por ejemplo, astroquery ) o el kit de herramientas NAIF SPICE y obtener o generar un archivo de kernel SPK para el asteroide (el normal, por ejemplo,de430.bspLas efemérides planetarias del JPL manejarán la parte Sol->Tierra)

Gracias por su respuesta, fue muy útil y he progresado mucho en mi proyecto. He hecho una pregunta de seguimiento en caso de que esté interesado en ver a dónde me llevó su respuesta :) astronomy.stackexchange.com/questions/29575/…