Für ein Projekt, an dem ich gerade arbeite, würde ich gerne wissen, wie man Geschwindigkeitskomponentenwerte extrahieren könnte ( , Und ) für einen Planeten aus der Stellarium-Software. Das aktuelle Informationsfeld für einen Planeten zeigt seine Umlaufgeschwindigkeit (das ist die resultierende Geschwindigkeit, ) gibt aber nicht das Individuum an , Und Werte, die zu dieser Geschwindigkeit führen.
Was ich im Wesentlichen brauche, ist, einen Geschwindigkeitskurs für den Planeten festzulegen.
Jede Hilfe (codebasiert/mathematisch) wäre sehr willkommen! Vielen Dank im Voraus!
Um die hilfreichen Kommentare unter Ihrer Frage weiterzuverfolgen, hier ist ein Programm, das das tut, was Sie brauchen, aber ohne Stellarium zu verwenden.
Wenn Sie Stellarium wirklich brauchen, wird hoffentlich jemand einen Weg finden, Ihnen zu helfen.
Dies verwendet das Python-Paket Skyfield und Sie können eine Menge Dokumentation auf dieser Seite und anderswo sehen. Im Moment habe ich die Ephemeride DE421 spezifiziert, aber es gibt viele zur Auswahl. Dieser funktioniert meiner Meinung nach nur +/- 100 oder 200 Jahre, andere gehen Tausende von Jahren in die Vergangenheit oder Zukunft.
Wenn Sie die Horizons- Weboberfläche von NASA JPL online verwenden möchten , finden Sie hier einige Anweisungen . Es verwendet einige der gleichen Entwicklungs-Ephemeriden , die Skyfield für Sie herunterlädt.
Hier ist eine Ausgabe für ungefähr 2019-10-02 07:18:59 UTC
den Moment, in dem Sie Ihre Frage gepostet haben.
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.
Hier ist als Beispiel ein Python 3-Skript für diese Ausgabe
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')
, dass Sie es von jedem Ordner aus ausführen können, ohne dass es eine neue (riesige) Datei in jedem Ordner herunterlädt.
Benutzer21
Astroschnapper
astroquery.jplhorizons
Modul über Python abfragbaräh
äh
Benutzer21
Siddharth Bhatnagar
Siddharth Bhatnagar