LEO-Metrik zur Bestimmung der Bodenspur eines Satelliten

Ich versuche, ein Verhalten auf LEO-Satelliten mit der Position über der Erde zu korrelieren. Das Verhalten findet über einen bestimmten Zeitraum statt. Ich versuche, eine Metrik zu entwickeln, die etwas über die Position relativ zur Position am Äquator und im Moment nicht über die tatsächliche Länge / Breite identifiziert. Der Grund ist, dass die Bodenspur der Umlaufbahn ihn in eine ganz andere Position gebracht haben könnte. Hier sind zum Beispiel ein paar ISS-Bodentracks von STK. Diese beiden Bodenspuren sind für meine Zwecke erheblich unterschiedlich, und jede Metrik sollte für jede unterschiedliche Werte ergeben, obwohl die Breite/Länge der ISS in jedem Fall nahezu identisch ist.

Geben Sie hier die Bildbeschreibung ein

Geben Sie hier die Bildbeschreibung ein

Jede Metrik, die ich finde, muss leicht berechenbar sein, entweder mit einer Python-Bibliothek oder etwas, das direkt berechnet werden kann. Man kann davon ausgehen, dass ich die TLE-Dateien für den Satelliten habe. Der Grund dafür ist, dass ich zu vielen davon statistische Daten erstelle und wirklich nicht die Zeit habe, jede dieser Instanzen in ein eigenständiges Programm zu stecken. Die Verwendung einer Web-API ist ebenfalls akzeptabel.

Worüber ich nachgedacht habe, ist eine Art Längenmaß des aufsteigenden Knotens, insbesondere wenn es reibungslos verfolgt werden kann, aber ich finde es schwierig herauszufinden, was das berechnen könnte. Irgendwelche anderen Vorschläge da draußen?

Grundsätzlich würde alles, was für die ISS funktionieren würde, für meine Bedürfnisse funktionieren. Meine Satelliten sind in LEO, gelegentlich Schub, meist kreisförmig, aber insgesamt ziemlich stabil.

Ich bin mir nicht sicher, wonach Sie suchen ... wenn Sie die Satellitenposition bezüglich der Erdoberfläche beschreiben möchten, dann scheint Lat/Lon/Alt eine natürliche Wahl zu sein ... könnten Sie etwas mehr erklären?
@Chris: Grundsätzlich möchte ich eine Zahl, die die beiden Bodenspuren voneinander unterscheiden kann, obwohl sie eine nahezu identische Breite / Länge haben
Der Hauptunterschied, den ich sehe, ist, dass einer ein aufsteigender Pass ist, während der andere absteigt - ist das das, was Sie unterscheiden möchten?
@Chris: Grundsätzlich möchte ich, dass eine Metrik die Bodenspur misst, die sie erstellt. Die Metrik sollte nicht viel anders sein, wenn ich die Zeit um 10 Minuten in beide Richtungen geändert habe, mit der entsprechenden Entfernungsänderung.
Hmmm ... ich versuche immer noch, mich darum zu kümmern. Angenommen, Sie sind nicht so besorgt über die Form der Bodenspur wie über ihre aktuelle "Position", suchen Sie nach etwas wie der Ortszeit am aufsteigenden Knoten? Oder sogar nur der Längengrad des aufsteigenden Knotens selbst?
Ich möchte nicht, dass es an die Tageszeit gebunden ist. Der Längengrad des aufsteigenden Knotens wäre jedoch von Interesse.
Würde der Steuerkurs des Raumfahrzeugs diesem Zweck dienen? Dann hätten Sie ein Lat/Lon/Heading-Trio für Ihre Position, das die Bodenspur eindeutig identifiziert und möglicherweise vom Boden aus beobachtbar ist.
Nun, es würde einen Punkt eindeutig identifizieren, aber ich möchte wirklich ein Diagramm der Leistung eines Subsystems im Vergleich zur Bodenspur erstellen. Lat/Long/Heading sind zu viele Datenpunkte...

Antworten (1)

Ich glaube, ich verstehe, wonach du suchst. Sie erwähnen den Längengrad des aufsteigenden Knotens, sind sich aber nicht sicher, wie Sie dies tatsächlich implementieren sollen. Hier ist wie.

Erstens der Längengrad des aufsteigenden Knotens (auch bekannt als Rektaszension des aufsteigenden Knotens) Ω für einen erdumkreisenden Satelliten ist definiert als der Winkel zwischen einem Referenzvektor (typischerweise der X Achse in einem erdzentrierten Trägheitsrahmen, der wiederum typischerweise mit einem Vektor vom Erdmittelpunkt zum ersten Punkt des Widders zusammenfällt ) und dem Punkt, an dem Ihre Umlaufbahn den Äquator kreuzt, während Sie aufsteigen (von Süden nach Norden gehen).

RAAN

Sie haben Glück, denn der zweizeilige Elementsatz-Dateityp enthält diesen Winkel direkt ( Spalten 18-25 in der zweiten Zeile ). Leider hat dies (noch) keinen Zusammenhang mit der Orientierung zur Erdoberfläche.

Für die Umrechnung von einem "Trägheits" Ω etwas relativ zur Erde, müssen Sie wissen, wie spät es ist. Siehe unten

Stundenwinkel

Hier zeigt der grüne Vektor auf das interessierende Objekt (hier der aufsteigende Knoten), der graue Vektor ist Ihre Bezugsachse (erster Punkt des Widders) und der gelbe Vektor ist der Nullmeridian (0 Grad Länge).

Um also den Winkel zu erhalten, an dem Sie interessiert sind (erdbezogener Längengrad des aufsteigenden Knotens oder der Winkel zwischen dem gelben und dem grünen Vektor), müssen Sie GMST davon subtrahieren Ω . Ihre TLE gibt Ihnen die Weltzeit (UT) für die Epoche an, aber Sie möchten vielleicht Ihre Umlaufbahn propagieren und diese für andere Zeiten berechnen (denken Sie nur daran, dass Ihr Knoten in LEO driften wird). Der Wechsel von UT zu GMST ist nicht besonders schwierig, aber da Sie Python in der Frage erwähnt haben, ist hier ein Paket, das diese Konvertierung für Sie durchführt . Um Stunden in einen Längengrad umzuwandeln, ist es einfach λ = GMT 24 360 .

Nun die Menge Ω λ wird der Längengrad (in Bezug auf die Erdoberfläche) Ihres aufsteigenden Knotens sein.

Eine Einschränkung ist, dass die Verwendung von GMST hier nur einige andere Effekte (Präzession und Nutation der Erdpole) ignoriert, aber im Allgemeinen sind diese Effekte im Vergleich zur GMST-Korrektur sehr gering, und wenn man bedenkt, dass Sie sowieso nur TLEs verwenden, ist es wahrscheinlich genau genug für was Sie wollen.

Eine andere Sache, die Sie im Auge behalten sollten - in STK wird diese Bodenspur in jeder Periode "schreiten", aber das liegt nur an der Art und Weise, wie sie angezeigt wird. Die tatsächliche Länge des aufsteigenden Knotens wird sich mit der Zeit kontinuierlich ändern, wenn sich die Erde unter Ihrer Umlaufbahn dreht (und wenn Ihr Flugzeug präzediert), also denken Sie daran.

Also verwende ich nur die RA des aufsteigenden Knotens und mache eine Berechnung basierend auf der Rotation der Erde. Brillant! Danke für die Hilfe!