Könnte mir bitte jemand erklären, wie ich Geschwindigkeitskomponentenwerte aus Stellarium für einen bestimmten Planeten extrahieren kann?

Für ein Projekt, an dem ich gerade arbeite, würde ich gerne wissen, wie man Geschwindigkeitskomponentenwerte extrahieren könnte ( v X , v j Und v z ) für einen Planeten aus der Stellarium-Software. Das aktuelle Informationsfeld für einen Planeten zeigt seine Umlaufgeschwindigkeit (das ist die resultierende Geschwindigkeit, v ) gibt aber nicht das Individuum an v X , v j Und v z 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!

Stellarium ist hier wirklich nicht das richtige Werkzeug. Betrachten Sie HORIZONS: ssd.jpl.nasa.gov/horizons.cgi#top oder die Liste unter astronomy.stackexchange.com/questions/13488
Zweiter Vorschlag von @barrycarter für HORIZONS. Es verfügt über eine 'Vektortabelle'-Ausgabeoption, die die Generierung einer kartesischen Zustandsvektortabelle (x,y,z, Geschwindigkeit in x,y,z) eines beliebigen Objekts in Bezug auf einen beliebigen Hauptkörper ermöglicht. Es ist über das astroquery.jplhorizonsModul über Python abfragbar
@astrosnapper ist es sicher , dass dies mit Stellarium nicht möglich ist? Dies sind ausgezeichnete Vorschläge, aber manchmal müssen die Leute die Werkzeuge verwenden, die sie im Moment zur Hand haben.
@barrycarter dito.
@uhoh Ich habe nicht gesagt, dass es mit Stellarium nicht geht. Stellarium muss diese individuellen Geschwindigkeiten natürlich berechnen, um sie anzuzeigen. Das Hauptziel von Stellarium ist es jedoch, die Dinge schön darzustellen. Es scheint zumindest schwierig zu sein, es dazu zu bringen, etwas zu tun, wofür es nicht entwickelt wurde. Wenn Sie zu F2 im Konfigurationsfenster gehen und Informationen auswählen, können Sie auf "Alle verfügbaren" klicken, was Ihnen viel mehr Informationen zu Planeten zeigt, aber ich glaube immer noch nicht, dass es Ihnen das gibt, was Sie wollen.
@barrycarter Vielen Dank für den Hinweis auf die HORIZONS-Ressource. Sieht so aus, als könnte dies mein Problem lösen!
@astrosnapper Vielen Dank auch an dich.

Antworten (1)

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 UTCden 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])
Vielen Dank für das Teilen des Python-Codes und für die HORIZONS- und DE-Links! Außerdem deckt DE421, die von Ihnen empfohlene Ephermis, die Jahre 1900 bis 2200 ab. Zufällig liegt der Zeitraum, der mich für mein Projekt interessiert, zwischen 1781 und 1881 für Uranus. Welches DE wäre Ihrer Meinung nach dafür am geeignetsten? DE414?
@SiddharthBhatnagar Ich bin kein Experte dafür, es gibt Zusammenfassungen der verschiedenen in Wikipedia und bei JPL , zum Beispiel deckt DE405 1600 bis 2200 (60 MB) und DE406 ab 3000 v. Chr. bis 3000 n. Chr. (290 MB). Es dauert eine Weile, die großen herunterzuladen, aber ich glaube nicht, dass dies die Geschwindigkeit beim Laufen beeinflusst. Stellen Sie sicher, dass Sie es so verwenden 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.
Macht Sinn, werde das im Hinterkopf behalten. Danke noch einmal.