Dificultad al tratar de usar el archivo SPK ast343de430.bsp

Estoy buscando un archivo NAIF SPK (por lo general, terminan en .bsp) para un grupo de asteroides en el marco de referencia J2000. encontré ast343de430.bspen ftp://ssd.jpl.nasa.gov/pub/eph/small_bodies/asteroids_de430/ast343de430.bsp. Parecen estar referenciados a J2000 y tienen 10(sol) como su referencia central, pero cuando voy a graficarlos se ven torcidos.

Cuando superpongo sus rastros de430.bsp, obtengo rastros de aspecto epicicloidal para cada cuerpo individual y la nube de puntos parece formar la forma característica del cinturón de asteroides pero compensada por 2-3 AU.

Estoy usando un lector SPK básico que he desarrollado yo mismo. El código está aquí (con licencia MIT, ya que surgió la pregunta a continuación). Este lector está de acuerdo con los archivos satelitales y de cuerpo grande publicados públicamente ( https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets/ , https://naif.jpl.nasa.gov/pub /naif/generic_kernels/spk/satellites/ )

¿Alguien ha usado este archivo y tuvo algún problema con él? ¡Gracias!

@uhoh No sé específicamente sobre los archivos de asteroides, pero puedo cargar todos los archivos de cuerpos grandes publicados públicamente + pluto ( naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/planets , naif.jpl.nasa. gov/pub/naif/generic_kernels/spk/satellites ) sin problemas. ¡Gracias!
Lo sentimos, esto es un poco fuera de tema, pero ¿su lector SPK es de código abierto por casualidad? Actualmente estoy tratando de leer un lector SPK / DAF yo mismo, y me está costando un poco.
@ChrisR MIT con licencia aquí: github.com/kghose/groho . Los bits del lector están aquí: github.com/kghose/groho/tree/master/src/orrery
@uhoh, encontré jplephem de Brandon Rhodes bastante agradable de usar.
Sugeriría usar CSPICE ( naif.jpl.nasa.gov/naif ) para leer archivos del kernel SPICE (SPK). Si no, el skyfield de Brandon Rhodes ( github.com/skyfielders/python-skyfield ), que es más reciente que pyephem. Explico el formato SPICE aquí: github.com/skyfielders/python-skyfield/issues/19 y he escrito un lector de archivos SPICE para Mathematica (pero escrito en Perl) en github.com/barrycarter/bcapps/blob/master/ASTRO/ bc-xsp2math.pl (otros archivos en ese directorio también pueden ser útiles). astronomy.stackexchange.com/questions/13488 también puede ser útil.
@barrycarter, ¿tiene un software goto que usa para verificar/verificar las cosas? Es un proyecto de pasatiempo y una de las cosas que disfruto hacer (más o menos...) es hacerlo todo yo mismo. Entonces, si a mi lector le falta algo que necesito, me gusta agregarlo yo mismo. En este caso, me gustaría estar seguro de que soy yo quien está torcido en lugar del archivo. ¡Gracias!
@KaushikGhose Usaría CSPICE (y HORIZONTES) como mi cheque. Si no están de acuerdo con sus resultados, probablemente sea su código el que está errado, no el de ellos.
Amigos, hubo un error en mi programa y pude continuar una vez que pude usar jplephempara trazar los datos del archivo. ¡Gracias por toda tu ayuda! Puedo escribir esto como una respuesta básicamente simplemente señalando la naturaleza no autónoma del archivo (hace referencia a todo 10pero no se contiene 10a sí mismo) y recomendándolo jplephemcomo una forma de verificación cruzada de los cálculos. ¡Solo necesitaba toda la camaradería de ustedes para recordar jplephem! ¡Gracias!

Respuestas (1)

El ast343de430.bsparchivo parece estar bien. Solo hay una cosa a tener en cuenta: todos los cuerpos en este archivo están referenciados al sol ( 10), pero las efemérides del sol no están incluidas. (La razón se da aquí : el autor de este archivo considera que la posición del Sol está mejor definida que la SSB). Por esta razón .bsp, se debe usar un archivo adicional que contenga la SSB, a menos que esté satisfecho con usar el Sol como origen del sistema de coordenadas.

Sin embargo, creo que definitivamente puedes aceptar este ;-)