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:
He aplicado las fórmulas matriciales sugeridas por la gente de DAMIT:
en el cual
Aquí es donde me confundo. Tengo entendido que mi visualización utiliza un sistema de coordenadas de la eclíptica heliocéntrica.
Configuré r_ast
la ubicación XYZ de un vértice en mi visualización. No estoy seguro de que esto sea correcto.
Luego calculo r_ecl
a 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?
r_ast
es 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_ast
vectores 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_ecl
tiene 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.bsp
Las efemérides planetarias del JPL manejarán la parte Sol->Tierra)
astrosnapper
ty.
astrosnapper
ty.
r_ast
las coordenadas rectangulares de la eclíptica heliocéntrica del asteroide? El sitio web lo describer_ast
como un vector en el "marco de coordenadas co-rotante del asteroide", pero no estoy seguro de lo que eso significa aquí.