¿Podría alguien explicarme cómo puedo extraer los valores de los componentes de velocidad de Stellarium para un planeta determinado?

Para un proyecto en el que estoy trabajando actualmente, me gustaría saber cómo se pueden extraer los valores de los componentes de velocidad ( v X , v y y v z ) para un planeta del software Stellarium. El panel de información actual de un planeta muestra su velocidad orbital (que es la velocidad resultante, v ) pero no le da al individuo v X , v y y v z valores que resultan en esa velocidad.

Lo que necesito esencialmente es establecer un rumbo de velocidad hacia el planeta.

¡Cualquier ayuda (basada en código/matemática) sería muy apreciada! ¡Gracias de antemano!

Stellarium realmente no es la herramienta adecuada aquí. Considere HORIZONTES: ssd.jpl.nasa.gov/horizons.cgi#top o la lista en astronomy.stackexchange.com/questions/13488
Segunda sugerencia de @barrycarter de HORIZONTES. Tiene una opción de salida de 'Tabla de vectores' que permite la generación de una tabla de vectores de estado cartesiano (x,y,z, velocidad en x,y,z) de cualquier objeto con respecto a cualquier cuerpo principal. Se puede consultar a través de Python a través del astroquery.jplhorizonsmódulo .
@astrosnapper ¿está seguro de que esto no es posible con Stellarium? Estas son excelentes sugerencias, pero a veces las personas necesitan usar las herramientas que tienen a mano por el momento.
@barrycarter ídem.
@uhoh, no dije que no se pudiera hacer con Stellarium. Stellarium obviamente tiene que calcular estas velocidades individuales para mostrarlas. Sin embargo, el objetivo principal de Stellarium es mostrar las cosas de forma agradable. Conseguir que haga algo para lo que no fue diseñado parece difícil, como mínimo. Si vas a F2 en la ventana de configuración y eliges información, puedes hacer clic en "todo disponible", que te muestra mucha más información sobre los planetas, pero aún no creo que te dé lo que quieres.
@barrycarter Muchas gracias por indicarme el recurso HORIZONTES. ¡Parece que esto podría resolver mi problema!
@astrosnapper Muchas gracias a ti también.

Respuestas (1)

Para dar seguimiento a los comentarios útiles debajo de su pregunta, aquí hay un programa que hará lo que necesita, pero sin usar Stellarium.

Si realmente necesita usar Stellarium, es de esperar que alguien encuentre la manera de ayudarlo.

Esto usa el paquete Skyfield de Python y puede ver mucha documentación en ese sitio y en otros lugares. Ahora mismo he especificado las efemérides DE421, pero hay muchas para elegir. Este solo funciona +/- 100 o 200 años creo, otros van por miles de años en el pasado o en el futuro.

Si desea utilizar la interfaz web Horizons de NASA JPL en línea, aquí tiene algunas instrucciones . Utiliza algunas de las mismas Efemérides de desarrollo que descarga Skyfield para usted.

Aquí hay algunos resultados para aproximadamente 2019-10-02 07:18:59 UTCcuál es el momento en que publicó su pregunta.

                    relative to solar system barycenter
               x            y           z            vx        vy          vz
Sun        -453731.5    1039601.2    451057.1     -0.015      -0.002       0.0
Earth     1.4764e+08   2.1455e+07   9.3005e+06    -4.9362     26.9295     11.6751
Jupiter  -2.2857e+07  -7.2214e+08  -3.0898e+08    12.9012      0.3382     -0.1690
Mumbai    1.4763e+08   2.1454e+07   9.3025e+06    -4.8369     26.5009     11.6749

                      relative to Earth's geocenter

Sun      -1.4809e+08  -2.0416e+07  -8.8494e+06     4.9215    -26.9312    -11.6754
Earth     0.           0.           0.             0.          0.          0.
Jupiter  -1.7050e+08  -7.4360e+08  -3.1828e+08    17.8374    -26.5912    -11.844
Mumbai   -5873.8303   -1362.2932     2072.3226  9.9338e-02    -4.2861e-01 -1.9206e-04

                           relative to Mumbai

Sun      -1.4809e+08  -2.0414e+07  -8.8515e+06     4.8222    -26.5026    -11.6752
Earth       5873.83     1362.29     -2072.32      -9.9338e-02  4.2861e-01  1.9206e-04
Jupiter  -1.7049e+08  -7.4359e+08  -3.1828e+08    17.7381    -26.1626    -11.8439
Mumbai    0.           0.           0.             0.          0.          0.

Aquí hay un script de Python 3 para esa salida como ejemplo

import numpy as np
import matplotlib.pyplot as plt
from skyfield.api import Topos, Loader, EarthSatellite

load    = Loader('~/Documents/fishing/SkyData')  # avoids multiple copies of large files
data    = load('de421.bsp')
ts      = load.timescale()

sun     = data['sun']
earth   = data['earth']
jupiter = data['jupiter barycenter']
mumbai  = earth + Topos(latitude_degrees  = 18.98,
                        longitude_degrees = 72.83,
                        elevation_m       =  20.0)
time    = ts.utc(2019, 10, 2, 7, 18, 59)

things     = (sun, earth, jupiter, mumbai)
names      = ('Sun', 'Earth', 'Jupiter', 'Mumbai')
positions  = [thing.at(time).position.km for thing in things]
velocities = [thing.at(time).velocity.km_per_s for thing in things]

for name, position, velocity in zip(names, positions, velocities):
    print(name, position, velocity)

print("relative to Earth's Geocenter")
for name, position, velocity in zip(names, positions, velocities):
    print(name, position-positions[1], velocity-velocities[1])

print("relative to Mumbai")
for name, position, velocity in zip(names, positions, velocities):
    print(name, position-positions[3], velocity-velocities[3])
¡Muchas gracias por compartir el código Python y por los enlaces HORIZONS y DE! Además, DE421, la efehermis que has recomendado, cubre los años 1900 a 2200. Sucede que para mi proyecto, el período de tiempo que me interesa es entre 1781 y 1881, para Urano. ¿Qué DE crees que sería más apropiado para esto? DE414?
@SiddharthBhatnagar No soy un experto en esto, hay resúmenes de varios en Wikipedia y en JPL , por ejemplo, DE405 cubre 1600 a 2200 (60 MB) y DE406 cubre 3000 a. C. a 3000 d. C. (290 MB). Se tarda un poco en descargar los grandes, pero no creo que afecte la velocidad mientras corres. Asegúrese de usarlo load = Loader('pathname')para que pueda ejecutar desde cualquier carpeta sin que descargue un archivo nuevo (enorme) en cada carpeta.
Tiene sentido, lo tendré en cuenta. Gracias de nuevo.