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 ( , y ) para un planeta del software Stellarium. El panel de información actual de un planeta muestra su velocidad orbital (que es la velocidad resultante, ) pero no le da al individuo , y 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!
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 UTC
cuá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])
load = Loader('pathname')
para que pueda ejecutar desde cualquier carpeta sin que descargue un archivo nuevo (enorme) en cada carpeta.
usuario21
astrosnapper
astroquery.jplhorizons
módulo .UH oh
UH oh
usuario21
Siddharth Bhatnagar
Siddharth Bhatnagar