StarAberrationDiagnostics

PixInsight-Update-Repository für das Script StarAberrationDiagnostics.

Diese Seite dient nur zur Information. In PixInsight wird die unten angegebene Repository-URL verwendet.

Repository-URL für PixInsight

https://tricx.de/pi-scripts/StarAberrationDiagnostics/

Installation in PixInsight

  1. Resources > Updates > Manage Repositories
  2. Repository hinzufügen: https://tricx.de/pi-scripts/StarAberrationDiagnostics/
  3. Check for Updates ausführen
  4. Starten über Script > Tricx > StarAberrationDiagnostics

Was das Script macht

Es diagnostiziert das optische System anhand der Sternformen über das ganze Bildfeld:

  • Vektorkarte der Elongation jedes Sterns, mit Stromlinien und einer Orientierungs-Heatmap
  • Sensorverkippung aus der FWHM über das Bildfeld (Quadranten wie in Siril, 11×11-Raster, Kippachse, Δz und Kippwinkel)
  • Korrektorabstand und Bildfeldwölbung aus dem radialen/tangentialen Muster am Bildrand
  • Komafeld und komafreier Punkt aus der Asymmetrie jedes Sterns, für die Kollimation
  • Der allen Sternen gemeinsame Nachführfehler, gemessen und optional abgezogen
  • Eine Bewertung der Ergebnisse mit Korrekturvorschlägen, in der Reihenfolge, in der man sie abarbeiten sollte
  • Eine Serienanalyse von Aufnahmen vor und nach dem Meridian-Flip, einschließlich der Plate-Solve-Verzeichnung

Screenshot

Dialog von StarAberrationDiagnostics mit der Vorschau der Vektorkarte

Handbuch: So funktioniert die Diagnose

Die folgenden Kapitel erklären jede Berechnung des Scripts: was jeder Fehler mit den Sternen macht und wie das Script ihn mathematisch erfasst. Sie bauen aufeinander auf, lassen sich aber auch einzeln lesen.

  1. Projektstatus und Haftungsausschluss
  2. Projektidee
  3. Sterne als Sonde der Optik
  4. Sterne finden und vermessen
  5. Die Form eines Sterns beschreiben
  6. Nachführ- und Guidingfehler
  7. Vektorkarte, Stromlinien und Heatmap
  8. Fokus, Sensorverkippung und Bildfeldwölbung
  9. Radiale und tangentiale Muster: der Korrektorabstand
  10. Koma und Kollimation
  11. Die Bewertung und die Reihenfolge der Korrekturen
  12. Serien: vor und nach dem Meridian-Flip
  13. Plate-Solve-Verzeichnung (SIP)
  14. Grenzen und gute Praxis

Vertiefung: Koma - wie sie entsteht, mit berechneten Strahlengängen und einem Simulator, und Typische Fehler eines Newton-Teleskops, mit berechneten Sternbildern.

0.Projektstatus und Haftungsausschluss

  • Frühe Phase. Das Script befindet sich in einer frühen Entwicklungsphase. Funktionen, Schwellen und Ausgaben können sich noch ändern.
  • Machbarkeitsstudie. Das Script ist ein Proof of Concept: Es erprobt, ob sich das optische System allein aus den Formen der Sterne zuverlässig diagnostizieren lässt. Es ist kein fertiges oder validiertes Messinstrument.
  • Mit KI entworfen. Das Script und dieses Handbuch wurden mithilfe von KI entworfen.
  • Kann Fehler enthalten. Berechnungen, Urteile und Erklärungen können falsch sein. Prüfen Sie die Ergebnisse mit Ihrer eigenen Erfahrung, bevor Sie Ihre Ausrüstung verstellen.
  • Rückmeldungen sind willkommen. Einen Fehler gefunden, eine Frage oder einen Vorschlag? Melden Sie sich gerne über die Kontaktseite.
  • Keine Gewähr. Das Script und diese Seite werden ohne jede Gewährleistung oder Garantie bereitgestellt. Die Nutzung erfolgt auf eigenes Risiko; für die Ergebnisse und alle Folgen ihrer Verwendung wird keine Haftung übernommen.

1.Projektidee

Das Script untersucht das optische System anhand von Light Frames - den Aufnahmen, die in einer Astrofoto-Nacht ohnehin entstehen. Spezielle Testaufnahmen sollen dafür nicht nötig sein: kein künstlicher Stern, keine bewusst defokussierten Sterne, keine Fokusreihe und kein Aufbau am Tag. Damit gehen unter Umständen Analysemöglichkeiten verloren; das ist eine bewusste Entscheidung.

Was dafür spricht

  • Kein zusätzlicher Aufwand. Jede Aufnahme einer Nacht ist zugleich eine Messung. Die Diagnose kostet keine Beobachtungszeit, und man braucht kein zusätzliches Zubehör.
  • Die Optik im Einsatz. Gemessen wird das System so, wie es tatsächlich fotografiert: mit Kamera, Korrektor und Filter, bei der Temperatur dieser Nacht und in der Lage, in der das Teleskop gerade zum Himmel zeigt.
  • Auch rückblickend. Ältere Aufnahmen lassen sich nachträglich auswerten, etwa vor und nach einer Justage oder vor und nach dem Meridian-Flip.

Was dabei verloren geht

  • Fehler unterhalb des Seeings. Seeing und Nachführung legen eine Unschärfe über jeden Stern. Kleine optische Fehler, die ein künstlicher Stern oder ein stark defokussiertes Sternscheibchen zeigen würde, bleiben darin verborgen.
  • Die Richtung eines Defokus. Eine einzelne fokussierte Aufnahme zeigt, wie stark ein Stern unscharf ist, aber nicht, ob vor oder hinter dem Fokus. Eine Fokusreihe würde das klären; hier hilft nur, in kleinen Schritten zu korrigieren und erneut zu messen.
  • Abhängigkeit vom Sternfeld. Das Script braucht genug gut verteilte Sterne bis in die Ecken. In sternarmen Feldern oder bei großen Nebeln werden die Aussagen unsicher.

Hinweise zu guten Aufnahmen

Das Script gibt Hinweise, welche Aufnahmen sich gut für die Analyse eignen; sie stehen gesammelt im Kapitel Grenzen und gute Praxis. Die meisten davon gelten ebenso für gute Astrofotos: gutes Seeing, eine saubere Nachführung, kurze Belichtungen und ein reiches Sternfeld. Die einzige Besonderheit ist, was ausgewertet wird: das einzelne kalibrierte Bild vor Registrierung, Stacking und Bearbeitung - nicht das fertige Bild.

2.Sterne als Sonde der Optik

Ein Stern ist so weit entfernt, dass er praktisch ein Punkt ist. Alles, was diesen Punkt auf dem Sensor zu einem kleinen Scheibchen, einem Oval oder einem kleinen Kometen macht, passiert unterwegs: in der Atmosphäre, in der Montierung, im Teleskop, im Korrektor, in der Kamera. Das Bild, das eine Punktquelle erzeugt, heißt Punktspreizfunktion (PSF). Weil jeder Stern im Bild durch einen anderen Teil des optischen Systems läuft, bilden die Sterne zusammen eine Karte davon, wie gut die Optik an jeder Stelle des Bildfelds arbeitet.

Das Script liest diese Karte. Für jeden Stern misst es vier Dinge:

  • Größe - wie weit das Licht verteilt ist, als Halbwertsbreite (FWHM, full width at half maximum).
  • Elongation - wie stark der Stern von einem Kreis abweicht.
  • Orientierung - in welche Richtung ein länglicher Stern zeigt.
  • Asymmetrie - ob das Licht zu einer Seite hin ausläuft, wie der Schweif eines Kometen.

Ein einzelner Stern sagt wenig aus: Seeing, Rauschen und Nachbarsterne stören jede Messung. Die Stärke der Methode liegt im Muster über Hunderte von Sternen. Ein Fokusfehler, ein verkippter Sensor, ein falscher Korrektorabstand, ein dejustierter Spiegel und ein Nachführfehler hinterlassen jeweils ein eigenes Muster, und die folgenden Kapitel erklären, welches - und wie das Script sie auseinanderhält.

Sternprofil
Abbildung 1 - Ein Schnitt durch einen Stern: die Helligkeit pro Pixel (Balken) und das angepasste Modell (Kurve). Die FWHM ist die Breite auf halber Höhe über dem Hintergrund.

Im Dialog: Der Dialog

3.Sterne finden und vermessen

Debayern von Farb-Rohbildern

Eine Farbkamera nimmt ein Bayer-Mosaik auf: Jedes Pixel sieht nur Rot, Grün oder Blau. Ein unbearbeitetes Mosaik lässt sich nicht direkt vermessen - das Farbmuster sähe in jedem Stern wie Struktur aus. Die üblichen Debayer-Verfahren ergänzen die fehlenden Farben durch Interpolation, die jeden Stern ein wenig und richtungsabhängig weichzeichnet. Das würde genau das verfälschen, was das Script misst.

Das Script verwendet deshalb das SuperPixel-Verfahren: Jeder 2×2-Block (ein rotes, zwei grüne, ein blaues Pixel) wird zu einem Farbpixel, ganz ohne Interpolation. Das Bild hat die halbe Auflösung, und jede in Mikrometer umgerechnete Länge verwendet den effektiven Pixelabstand, also das Doppelte des physischen.

SuperPixel-Debayering
Abbildung 2 - SuperPixel-Debayering: Jeder 2×2-Block des Mosaiks wird zu einem Pixel.

Detektion

Der StarDetector von PixInsight findet die Kandidaten. Ist die Option Max. stars gesetzt, werden nur die hellsten behalten (nach Fluss sortiert) - sie liefern die zuverlässigsten Fits, und der Fit ist der langsame Teil.

Der PSF-Fit

Für jeden Kandidaten passt DynamicPSF von PixInsight ein elliptisches Modell an die Pixel innerhalb eines Suchfensters an (der search radius). Zwei Modelle stehen zur Wahl. In den Koordinaten u, v entlang der Achsen der Ellipse (gedreht um den Winkel θ) lauten sie:

Gauß:      I(u, v) = B + A · exp( −u²/(2σx²) − v²/(2σy²) )
Moffat:     I(u, v) = B + A / ( 1 + u²/αx² + v²/αy² )β
B Hintergrund, A Amplitude, σ bzw. α die Breite entlang jeder Achse, β die Steilheit der Moffat-Flanken

Das Moffat-Profil hat breitere Flanken als die Gaußkurve und beschreibt reale Sterne (Seeing, Beugung) oft besser; die Gaußkurve ist schneller. Aus den Breiten berechnet das Script die FWHM entlang jeder Achse:

Gauß:      FWHM = 2 √(2 ln 2) · σ ≈ 2,3548 · σ
Moffat:     FWHM = 2 α · √(21/β − 1)

Den Gauß-Faktor für einen Moffat-Fit zu verwenden wäre ein grober Fehler: Je nach β würde er die FWHM um den Faktor 1,5 bis 2,7 und mehr aufblähen.

Qualitätskontrolle

  • Detektionsschwelle. Ein Fit zählt nur oberhalb dieser Signifikanz. Kommt kein einziger Stern durch, senkt das Script die Schwelle schrittweise (×0,6, bis hinunter auf 0,2), statt aufzugeben, und meldet das.
  • MAD-Ausreißer. DynamicPSF liefert für jeden Fit die mittlere absolute Abweichung (MAD) zwischen Modell und Pixeln. Hotpixel, verschmolzene Paare oder sehr schwache Sterne ergeben schlecht passende Modelle mit großer MAD. Ein Fit wird verworfen, wenn
    MAD > Faktor · Median(MAD aller Fits)
    Die Grenze passt sich so dem Rauschen jeder Aufnahme an, statt eine feste Zahl zu sein.
  • Sättigung wird für die Asymmetriemessung geprüft (Kapitel 9): Ein Stern, dessen Maximum 95 % des Wertebereichs erreicht, hat eine abgeflachte Spitze und wird dort übersprungen.

Im Dialog: Setup › General, Setup › Star detection

4.Die Form eines Sterns beschreiben

Die angepasste Ellipse hat eine große Halbachse a (die größere FWHM) und eine kleine Halbachse b, und ihre große Achse zeigt in eine Richtung. Drei Zahlen beschreiben, wie länglich sie ist - sie tragen dieselbe Information, verhalten sich in Rechnungen aber unterschiedlich:

Exzentrizität   e = √(1 − b²/a²)
Elliptizität     ε = 1 − b/a = 1 − √(1 − e²)
Verzerrung       χ = (a² − b²)/(a² + b²) = e²/(2 − e²)
und zurück: e = √(2χ/(1 + χ))
Formmaße
Abbildung 3 - Die drei Maße für dieselben Sterne. Die Exzentrizität steigt schon bei leicht ovalen Sternen schnell an und eignet sich deshalb gut für die Farben der Karte; die Elliptizität liest sich wie ein Prozentwert (0,10 = die kurze Achse ist 10 % kürzer); die Verzerrung lässt sich korrekt addieren.

Die Exzentrizität färbt die Vektorkarte. Die Elliptizität wird überall dort verwendet, wo ein Urteil gefällt wird, weil sie leicht zu lesen ist. Die Verzerrung dient zum Mitteln und Abziehen: Nur ihre Komponenten addieren sich (näherungsweise) linear, wenn zwei Elongationen zusammenkommen, z. B. eine der Optik und eine der Nachführung.

Richtung und das 180°-Problem

PSF-Ellipse
Abbildung 4 - Die angepasste Ellipse in Bildkoordinaten: x nach rechts, y nach unten. Der Winkel ψ der großen Achse wird von der x-Achse aus gezählt; 90° zeigt nach unten.

Alle Richtungen im Script werden in Bildkoordinaten angegeben: x nach rechts, y nach unten, 0° = rechts, 90° = unten. DynamicPSF misst seinen Winkel θ auf dem Bildschirm gegen den Uhrzeigersinn, daher ist die Richtung der großen Achse in Bildkoordinaten ψ = −θ.

Eine Ellipse sieht nach einer Drehung um 180° gleich aus: Ihre Orientierung ist nur modulo 180° definiert. Die Winkel direkt zu mitteln geht deshalb schief - zwei fast waagerechte Sterne bei 10° und 170° ergäben im Mittel senkrechte 90°. Die übliche Abhilfe, die das Script überall verwendet, ist der doppelte Winkel: Jede Orientierung wird zu einem Vektor mit dem Winkel 2ψ, die Vektoren werden gemittelt, und der Winkel des Mittelwerts wird wieder halbiert.

Doppelter Winkel
Abbildung 5 - Orientierungen mitteln: naiv (links) und über den doppelten Winkel (rechts).
χ1 = χ · cos 2ψ,    χ2 = χ · sin 2ψ
mittlere Orientierung = ½ · atan2( Σ w·χ2, Σ w·χ1 )
(χ1, χ2) ist der „Formvektor“ eines Sterns; w ein optionales Gewicht

Im Dialog: Setup › Star detection, Maps › Shape

5.Nachführ- und Guidingfehler

Der Fehler

Folgt die Montierung dem Himmel nicht exakt - periodischer Fehler, schlechtes Guiding, Wind, Kabelzug, ein sich verbiegendes Leitrohr -, bewegt sich das ganze Bild während der Belichtung ein wenig. Jeder Stern wird um den gleichen Betrag in die gleiche Richtung verschmiert, egal wo er im Bildfeld steht. Dieser Fingerabdruck unterscheidet die Nachführung von allen optischen Fehlern, die sich über das Bildfeld ändern.

Die Messung

Das Script bildet den Median der Formvektoren aller Sterne. Der Median ignoriert die Minderheit der Sterne, die von der Optik am Rand verzogen sind, und findet das, was fast alle Sterne gemeinsam haben:

m1 = Median(χ1),   m2 = Median(χ2)
χT = √(m1² + m2²),   Richtung = ½ · atan2(m2, m1)
χT wird für den Bericht in eine Elliptizität zurückgerechnet

Mit der Option Subtract it from all stars wird dieser gemeinsame Anteil vor jeder weiteren Analyse von jedem Stern abgezogen - durch Subtraktion der Vektoren, und genau deshalb wird die Verzerrung χ verwendet:

χ′1 = χ1 − m1,   χ′2 = χ2 − m2
Nachführanteil
Abbildung 6 - Links: Ein Nachführfehler zieht jeden Stern gleich in die Länge. Rechts: Ihn zu entfernen ist eine Vektorsubtraktion in der Ebene der Formvektoren.

Korrigiert werden nur die Formen; die Sterngrößen (FWHM) bleiben wie gemessen. Der Nachführfehler würde sonst die optischen Muster überdecken, vor allem das radiale/tangentiale Muster aus Kapitel 8.

Den Wert lesen

Gemeinsame ElliptizitätUrteil
unter 0,03OK - die Nachführung ist sauber
0,03 - 0,07Hinweis
über 0,07Maßnahme: Guiding (RMS, Kalibrierung), Balance, Kabelzug, Wind und Verbiegung prüfen

Eine Richtung, die über mehrere Aufnahmen gleich bleibt, deutet auf die Montierung; eine, die wechselt, auf Wind oder zufällige Guidingfehler. Die Serienanalyse vergleicht beide Seiten eines Meridian-Flips, bei dem sich die Balance ändert.

Im Dialog: Setup › General, Setup › Star detection, Maps › Shape

6.Vektorkarte, Stromlinien und Heatmap

Die Vektorkarte zeichnet jeden Stern als Ellipse an seiner Position, Länge und Farbe nach der Exzentrizität (blau rund bis rot länglich), die Richtung nach ψ. Hunderte kleiner Ellipsen sind schwer zu lesen, deshalb zeigen zwei Ebenen den großräumigen Trend.

Stromlinien

An jedem Punkt des Bildes berechnet das Script eine geglättete Orientierung aus den umliegenden Sternen: ein gaußgewichtetes Mittel über den doppelten Winkel innerhalb eines Radius R (ein Prozentsatz der Bilddiagonale). Fast runde Sterne haben eine schlecht definierte Richtung, deshalb wird jeder Stern zusätzlich mit dem Quadrat seiner Exzentrizität gewichtet:

w = exp( −d² / (2 (R/2)²) ) · e²,    für alle Sterne im Abstand d < R
Orientierung = ½ · atan2( Σ w · sin 2ψ, Σ w · cos 2ψ )

Eine Stromlinie folgt diesem Feld in kleinen Schritten. Da eine Orientierung keine Pfeilspitze hat, wählt das Script bei jedem Schritt diejenige der beiden Richtungen, die den vorigen Schritt fortsetzt - sonst könnte die Linie zufällig umkehren. Stromlinien zeigen auf einen Blick, ob die Sterne radial, konzentrisch oder alle in eine Richtung ausgerichtet sind.

Orientierungs-Heatmap

Die Heatmap (nach der Seti Astro Suite) passt zwei glatte Polynomflächen eines wählbaren Grades über das ganze Bild an, eine an sin 2θ und eine an cos 2θ aller Sterne, mit auf ±1 normierten Koordinaten. Drei Runden 3σ-Clipping entfernen Ausreißer. An jedem Punkt wird der angepasste Winkel θ = ½ · atan2(s, c) als Farbton eines Farbkreises dargestellt. Eine Ebene (Grad 1) kann nur einen gleichmäßigen Verlauf zeigen; ab Grad 2 lassen sich die gekrümmten Muster der Koma abbilden, das Ergebnis wird aber empfindlicher für Rauschen an den Rändern, wo wenige Sterne es stützen.

Im Dialog: Maps › Shape

7.Fokus, Sensorverkippung und Bildfeldwölbung

Die Physik: Unschärfe durch Defokus

Das Teleskop formt einen Lichtkegel, der in der Brennebene zu einem Punkt zusammenläuft. Sitzt der Sensor um Δz vor oder hinter dieser Ebene, schneidet er den Kegel und nimmt statt eines Punktes ein kleines Scheibchen auf. Aus der Geometrie des Kegels (seine Öffnung ist die Apertur D geteilt durch die Brennweite f, also 1/N):

Unschärfe durch Defokus
Abbildung 7 - Ein um Δz defokussierter Sensor nimmt ein Scheibchen mit dem Durchmesser d = Δz/N auf.
d = Δz / N    ⇔    Δz = N · d,    N = f / D
ein schnelles f/4-System reagiert auf dasselbe Δz doppelt so stark wie ein f/8-System

Auch der fokussierte Stern ist kein Punkt (Seeing, Beugung, Optik): Er hat eine Größe FWHM0. Unschärfen unabhängiger Herkunft addieren sich näherungsweise quadratisch, daher folgt die zusätzliche Unschärfe durch Defokus aus der gemessenen FWHM als

d = √( FWHM² − FWHM0² ) · p,    Δz = N · d
p = (effektiver) Pixelabstand, rechnet Pixel in Mikrometer um

Das ist die Brücke von einer Sterngröße in Pixeln zu einem Abstand in Mikrometern - und sie braucht Pixelabstand, Brennweite und Öffnung, die das Script aus dem Dialog oder dem FITS-Header übernimmt (XPIXSZ, FOCALLEN, APTDIA oder FOCRATIO).

Das Vorzeichen geht verloren. Ein Scheibchen sieht gleich aus, egal ob der Sensor vor oder hinter dem Fokus sitzt. Aus einer Aufnahme ist Δz nur ein Betrag - in welche Richtung eine Kippschraube zu drehen ist, muss man ausprobieren (siehe Abbildung 11).

Wie die Fehler aussehen

Verkippung und Bildfeldwölbung
Abbildung 8 - Links: Ein verkippter Sensor ist entlang einer Linie im Fokus und zu beiden Rändern hin defokussiert, in entgegengesetzte Richtungen. Rechts: Eine gewölbte Bildschale trifft einen ebenen Sensor nur in der Mitte; die Sterne wachsen zu allen Rändern hin gleichmäßig.
  • Defokus macht jeden Stern größer - aber gleichmäßig, er zeigt sich also eher in der absoluten FWHM als in ihrem Muster. Ein verräterisches Zeichen ist eine Mitte, die weicher ist als der Rand: Es wurde auf den Rand fokussiert, oder ein gewölbtes Feld wurde außerhalb der Mitte fokussiert.
  • Sensorverkippung macht eine Seite des Bildes weich und die gegenüberliegende scharf: ein einseitiger Verlauf.
  • Bildfeldwölbung macht den Rand auf allen Seiten gleich weich: eine symmetrische Schüssel.

Methode 1: Quadranten wie in Siril

Wie bei Sirils „Show tilt“ wird das Bild in vier Quadranten sowie einen inneren Kreis und einen äußeren Ring um die Mitte eingeteilt (Radien bezogen auf die halbe Diagonale R). Jeder Bereich erhält das 25%-getrimmte Mittel der Sterngrößen (FWHMx+FWHMy)/2 - das unterste und oberste Viertel wird vor der Mittelung verworfen, was es robust gegen einzelne schlechte Fits macht.

Siril-Quadranten
Abbildung 9 - Die Bereiche der Quadrantenauswertung (O/U = oben/unten, L/R = links/rechts).
Verkippung = max(m1…m4) − min(m1…m4),   in % ihres Mittels
Off-Axis-Aberration = maußen − minnen

Das 11×11-FWHM-Raster ist dieselbe Idee in feinerem Maßstab: ein getrimmtes Mittel pro Zelle, stärker verrauscht, aber detailreicher. Die Kippachse wird aus den Quadranten abgeleitet. Der Überschuss jedes Quadranten gegenüber der Mitte bekommt ein Vorzeichen, der Gradient einer Ebene durch die vier Werte liefert die Richtung des steilsten Anstiegs und, bei bekannter Optik, den Kippwinkel:

eq = sign(Fq² − F0²) · √|Fq² − F0²|,   F0 = minnen
gx = [(eOR + eUR) − (eOL + eUL)] / W,   gy = [(eUL + eUR) − (eOL + eOR)] / H
Achsrichtung = atan2(gy, gx) (mod 180°),   Kippwinkel = atan( N · √(gx² + gy²) )

Pro Quadrant listet die Konsole außerdem den Defokus Δz = N · √(Fq² − F0²) · p und den Winkel atan(Δz / Abstand) zur Quadrantenmitte bei einem Viertel der Breite und Höhe.

Methode 2: die FWHM²-Fläche

Quadranten vermischen Verkippung und Wölbung: Ein gewölbtes Feld macht alle vier Quadranten weich, eine Verkippung nur zwei. Das Script passt deshalb eine glatte Fläche an die quadrierte FWHM jedes Sterns an, mit X und Y als Position relativ zur Mitte, geteilt durch die halbe Diagonale:

FWHM² = c0 + c1·X + c2·Y + c3·(X² + Y²)
FWHM²-Fläche
Abbildung 10 - Ein Schnitt durch die Fläche entlang X. Der symmetrische Term c₃ ist die Wölbung (gestrichelt); die linearen Terme c₁, c₂ kippen die Schüssel zu einer Seite (durchgezogen) - das ist die Sensorverkippung.

Angepasst werden die Quadrate, weil sich Unschärfen quadratisch addieren (siehe oben): Der symmetrische Term misst dann die Wölbung, die linearen Terme die Verkippung, und beide verfälschen sich nicht mehr gegenseitig. Der Fit ist ein robuster Kleinste-Quadrate-Fit: vier Runden, in denen Sterne mit einem Residuum über 3 σ (σ aus dem Median der absoluten Residuen) ausgelassen werden. Aus den Koeffizienten:

FWHMMitte = √c0,   FWHMRand = √(c0 + c3)
G = √(c1² + c2²),   weichere Seite in Richtung atan2(c2, c1)
FWHMweiche Seite = √(c0 + G),   FWHMscharfe Seite = √(c0 − G)
ΔzRand = N · √G · p,   Kippwinkel = atan( ΔzRand / (R · p) )
Wölbungstiefe = N · √|c3| · p
schärfster Punkt = Mitte − (c1, c2) / (2 c3) · R   (falls c3 > 0)

Ein negatives c3 bedeutet, dass der Rand schärfer ist als die Mitte - die Aufnahme ist defokussiert oder wurde auf den Rand fokussiert. Das Script meldet außerdem einen Defokus-Verdacht, wenn die Median-FWHM im inneren Drittel (r < ⅓ R) um mehr als 5 % größer ist als im äußeren Drittel (r > ⅔ R). Solange das der Fall ist, sind die Urteile zu Verkippung und Abstand unzuverlässig - deshalb steht der Fokus in der vorgeschlagenen Reihenfolge an erster Stelle.

Der 3D-Kippplot

Vorzeichen-Mehrdeutigkeit
Abbildung 11 - Zwei Verkippungen in entgegengesetzte Richtungen erzeugen dieselbe Unschärfe. Der 3D-Plot zeichnet deshalb beide Ebenen.

Der Pseudo-3D-Plot setzt die vier Quadrantenwerte an ihre physischen Positionen auf dem Sensor (in mm) als vorzeichenbehaftetes Δz und legt eine Ebene durch die Mitte. Weil das Vorzeichen von Δz unbekannt ist (siehe den Hinweis oben), zeichnet er zwei spiegelbildliche Varianten. Beide sind gleich plausibel; die Kippachse ist für beide dieselbe.

Urteile

WertOKHinweisMaßnahme
Verkippung: (FWHMweich − FWHMscharf) / ihr Mittel< 5 %5 - 10 %> 10 %
Bildfeldwölbung: (FWHMRand − FWHMMitte) / FWHMMitte≤ 15 %15 - 30 %> 30 %
Fokus: Defokus-Verdacht oder c3 < 0zuerst neu fokussieren

Im Dialog: Setup › General, Maps › Size, Maps › Tilt

8.Radiale und tangentiale Muster: der Korrektorabstand

Der Fehler

Ein Komakorrektor, Bildfeldebner oder Reducer ist für genau einen Abstand zum Sensor ausgelegt (den Backfokus, z. B. 55 mm). Im falschen Abstand korrigiert er zu wenig oder zu viel, und die Sterne am Bildrand werden länglich - entweder radial, wie Speichen zur Mitte zeigend, oder tangential, entlang von Kreisen um die Mitte liegend. In einem unkorrigierten System (einem Newton ohne Korrektor) ist radiale Elongation schlicht seine natürliche Koma.

Radial und tangential
Abbildung 12 - Radiale (links) und tangentiale (rechts) Elongation am Bildrand.

Die Messung

Für jeden Stern beim Positionswinkel φ (von der Bildmitte aus gesehen) misst das Script, wie viel seiner Elongation entlang des Radius liegt:

Radiale Elliptizität
Abbildung 13 - Der radiale Anteil der Elliptizität eines Sterns.
eps_rad = ε · cos 2(ψ − φ)
+ε für einen Stern entlang des Radius, −ε quer dazu, 0 bei 45° - Nachführfehler heben sich über einen ganzen Ring auf

Die Sterne werden in sechs Ringe eingeteilt (r/R = 0-⅙, ⅙-⅓, … 5/6-1). Pro Ring listet die Konsole die Anzahl der Sterne, die Median-FWHM, die Elliptizität, eps_rad, die radiale Ausläuferrichtung asym_rad (Kapitel 9) und einen Übereinstimmungswert zwischen Ausläufer und Elongation (+1: der Ausläufer liegt entlang der Elongation, d. h. die Koma formt den Stern). Das Urteil verwendet den Median von eps_rad im äußeren Drittel (r > ⅔ R), mit mindestens 10 Sternen:

eps_rad (äußeres Drittel)MusterFaustregel des Scripts
|eps_rad| ≤ 0,04keinesOK
> +0,04 (Maßnahme ab 0,08)radialmit Korrektor: unterkorrigiert - den Abstand Korrektor-Sensor in kleinen Schritten vergrößern (z. B. 0,5 mm); unkorrigiertes System: normale Koma
< −0,04 (Maßnahme unter −0,08)tangentialden Abstand verkleinern - Defokus zusammen mit Bildfeldwölbung sieht aber genauso aus, daher zuerst den Fokus prüfen

Die Richtung der Korrektur ist eine Faustregel, die für viele Korrektoren gilt; das letzte Wort hat die Anleitung Ihres Korrektors. Ändern Sie den Abstand in kleinen Schritten und messen Sie erneut.

Wie unter- und überkorrigierte Koma im Bildfeld aussieht, zeigt der Simulator auf der Seite Koma - wie sie entsteht.

Im Dialog: Maps › Size, Maps › Shape

9.Koma und Kollimation

Der Fehler

Koma macht Sterne zu kleinen Kometen: ein heller Kern mit einem Lichtfächer auf einer Seite. Bei einem Newton wächst sie linear mit dem Abstand von der optischen Achse und zeigt von ihr weg; ein Komakorrektor beseitigt sie - oder überkorrigiert sie, wenn er im falschen Abstand sitzt. Der Punkt, an dem die Koma verschwindet, markiert die optische Achse. Sind die Spiegel dejustiert (oder ist ein Korrektor dezentriert oder verkippt), wandert dieser komafreie Punkt von der Bildmitte weg.

Warum Koma entsteht, wie stark sie bei welchem Öffnungsverhältnis ist und wie Korrektor und Dejustage das Bild verändern, erklärt die Seite Koma - wie sie entsteht, mit berechneten Strahlengängen und einem interaktiven Simulator.

Warum der Fit sie nicht sieht

Gauß und Moffat sind punktsymmetrische Modelle: Sie legen eine Ellipse um den Kern und können nicht erkennen, auf welcher Seite der Fächer liegt. Das Script misst die Asymmetrie stattdessen direkt an den Pixeln, relativ zum angepassten Zentrum:

Koma-Stern
Abbildung 14 - Ein Stern mit Koma. Der Flussschwerpunkt liegt gegenüber dem Fit-Zentrum zum Ausläufer hin versetzt; das dritte Moment gewichtet das äußere Licht stärker und zeigt ihn deutlicher.
  • Aperturradius rap = 4σ + 1 px (begrenzt auf 3 … 25 px), σ = FWHM/2,3548; Hintergrund = Median eines Rings von rap+2 bis rap+6 px.
  • Gewichte w = Pixelwert − Hintergrund (nur positive).
  • Sterne am Bildrand, ohne Hintergrundring oder mit einem Maximum bei 95 % des Wertebereichs (gesättigt) werden übersprungen.
Schwerpunktversatz = Σ w · d / Σ w
drittes Moment   m3 = Σ w · ρ² · ρ / Σ w,   ρ = d / σ
d = Pixelposition relativ zum Fit-Zentrum; ρ macht m3 unabhängig von der Sterngröße

Das Komafeld

Koma dritter Ordnung ist linear in der Feldposition. Um den komafreien Punkt P0 wird der Asymmetrievektor eines Sterns an der Position P deshalb modelliert als

m3(P) = k · (P − P0) / R
linear in (k, ox, oy):  m3,x = k·X + ox,  m3,y = k·Y + oy
P0 = Mitte − (ox, oy) / k · R
Komafeld
Abbildung 15 - Das Komafeld. Bei k > 0 zeigen die Ausläufer von P₀ weg, bei k < 0 zu P₀ hin.
  • k > 0: Die Ausläufer zeigen von P0 weg - unterkorrigiert (normal für einen Newton ohne Korrektor).
  • k < 0: Die Ausläufer zeigen zu P0 hin - überkorrigiert.

Nötig sind mindestens 30 Sterne mit gültiger Asymmetrie. Der Fit ist robust (3σ-Clipping, drei Runden). Seine Unsicherheit stammt aus einem Bootstrap: Der Fit wird 100-mal an zufälligen Stichproben der Sterne wiederholt; die Streuung von k ist ihre Standardabweichung, die von P0 eine robuste Streuung (MAD), weil einige Stichproben mit k ≈ 0 P0 beliebig weit wegwerfen können. Die Koma gilt als signifikant, wenn |k| > 3 · Streuung(k). Zur Gegenprobe wird dasselbe Modell an die Schwerpunktversätze angepasst.

Kollimation

Das Script empfiehlt eine Kollimation, wenn P0 mehr als 15 % der halben Diagonale von der Mitte entfernt liegt und mehr als das Doppelte seiner eigenen Unsicherheit. Eine Vorsichtsmaßnahme ist eingebaut: Eine ungleichmäßige Nachführdrift fügt eine konstante Asymmetrie hinzu und verschiebt P0 entlang der Drift. Liegt der Versatz innerhalb von 20° zur Nachführ-Elongation, weist der Vorschlag darauf hin - vergleichen Sie mehrere Aufnahmen, bevor Sie die Kollimationsschrauben anfassen.

Koma-Stromlinien

Anders als eine Orientierung ist m3 ein echter Vektor mit Richtung. Die Koma-Stromlinien folgen seinem gaußgeglätteten Feld (mindestens 8 Sterne im Radius) mit Pfeilspitzen und enden dort, wo das Feld schwächer wird als ein eingestellter Prozentsatz seines starken Teils - so täuscht die komafreie Zone keine Richtung vor. Unterkorrigierte Koma sieht aus wie eine Quelle, überkorrigierte wie eine Senke.

Im Dialog: Maps › Coma

10.Die Bewertung und die Reihenfolge der Korrekturen

Der Reiter Assessment macht aus den Werten der vorigen Kapitel Befunde, jeweils bewertet als OK, Note (Hinweis) oder Action (Maßnahme), mit dem, was gefunden wurde, und dem, was zu tun ist. Wo die Richtung einer Korrektur vom optischen System abhängt, entscheidet der Optics type (unbekannt, unkorrigiert, mit Korrektor) über die Formulierung.

Die Maßnahmen werden in die Reihenfolge sortiert, in der man sie am besten abarbeitet, weil sich die Fehler gegenseitig überdecken:

  1. Fokus - eine defokussierte Aufnahme verfälscht die Urteile zu Verkippung und Abstand.
  2. Nachführung - eine gemeinsame Elongation verdeckt das radiale/tangentiale Muster und verschiebt den komafreien Punkt.
  3. Verkippung - eine einseitige Unschärfe, behoben am Kameraanschluss.
  4. Kollimation - verschiebt die optische Achse, um die herum alles andere gemessen wird.
  5. Abstand - der Korrektorabstand, beurteilt aus dem Randmuster.
  6. Koma - Stärke und Vorzeichen des Komafelds.
  7. Bildfeldwölbung - was am Rand übrig bleibt.
  8. Daten - weniger als 150 Sterne machen alle Urteile unsicher.

Ausnahme bei deutlicher Dejustage: Eine verkippte optische Achse kippt auch die Bildebene gegen den Sensor und erscheint dann als Verkippung - das Script kann das nicht von einem schief sitzenden Sensor unterscheiden. Liegt der komafreie Punkt deutlich neben der Bildmitte und wird zugleich eine Verkippung gemeldet, zuerst kollimieren und neu messen, erst dann die Verkippung korrigieren.

Alle Schwellen sind Ausgangspunkte, keine kalibrierten Grenzwerte. Seeing, Belichtungszeit und Sternanzahl verändern, was eine einzelne Aufnahme aussagen kann.

Woher diese Fehler bei einem Newton typischerweise kommen und was jeweils hilft, fasst die Seite Typische Fehler eines Newton-Teleskops zusammen.

Im Dialog: Assessment

11.Serien: vor und nach dem Meridian-Flip

Warum ein Flip hilft

Eine deutsche Montierung muss das Teleskop umschlagen, wenn das Objekt den Meridian überquert. Nach diesem Meridian-Flip ist der Tubus relativ zur Schwerkraft um 180° um seine Achse gedreht. Alles, was unter seinem eigenen Gewicht durchhängt oder sich verschiebt - ein Okularauszug, eine Kamera an einem langen Adapter, ein Hauptspiegel in seiner Fassung („Mirror Flop“), ein lockerer Korrektor -, wandert auf die andere Seite. Alles, was im optischen Zug fest sitzt, bleibt, wie es ist.

Meridian-Flip
Abbildung 16 - Vor und nach dem Flip wirkt die Schwerkraft von entgegengesetzten Seiten auf den optischen Zug.

Ein Vergleich von Aufnahmen beider Seiten beantwortet deshalb eine Frage, die eine einzelne Aufnahme nicht beantworten kann: Ist ein Fehler fest oder mechanisch? Eine Verkippung, die mit der Seite wechselt, verlangt das Festziehen von Okularauszug und Adaptern - ein Justieren der Verkippung würde nur für eine Seite passen. Eine Verkippung, die bleibt, verlangt die Kippjustage.

Auf welcher Seite liegt eine Aufnahme?

  1. Manuell - in der Liste gesetzt.
  2. PIERSIDE - geschrieben von N.I.N.A. und den meisten ASCOM/INDI-Aufnahmeprogrammen (West vor dem Meridian, Teleskop zeigt nach Osten; Ost nach dem Flip).
  3. Ein Kamerawinkel, der sich um 180° dreht - z. B. ROTATOR, wie ihn der ASIAIR aus seinem Plate Solve schreibt (auch POSANGLE, ANGLE, CROTA2). Aufnahmen innerhalb von 45° um den ersten Winkel bilden eine Gruppe, die innerhalb von 45° um den Gegenwinkel die andere; die Stundenwinkel ihrer Mitglieder benennen die Gruppen.
  4. Der Stundenwinkel - aus der Uhrzeit, der Rektaszension α des Objekts und der geografischen Länge λ des Standorts (Ost positiv):
GMST = 280,46061837° + 360,98564736629° · d   (d = Tage seit 2000-01-01 12:00 UT)
Stundenwinkel HA = GMST + λ − α
HA ≤ −0,05 h: West (der Meridian ist noch nicht erreicht) · HA ≥ +0,5 h: Ost · dazwischen: unbekannt, der Flip kann später kommen

Der Vergleich

Jede Aufnahme wird mit denselben Einstellungen vermessen. Pro Seite und Wert bildet das Script den Median und eine robuste Streuung σ = 1,4826 · MAD, die die Standardabweichung schätzt, ohne sich von einzelnen Ausreißern täuschen zu lassen. Der Standardfehler eines Medians beträgt etwa 1,253 σ/√n. Ein Wert ändert sich beim Flip, wenn beides gilt:

|MedianOst − MedianWest| > 3 · 1,253 · √( σW²/nW + σO²/nO )
und die Differenz übersteigt ein praktisches Minimum
Serienvergleich
Abbildung 17 - Ein Wert, der beim Flip springt: Die Differenz der Mediane ist groß im Vergleich zur Streuung innerhalb jeder Seite.
WertPraktisches MinimumEine Änderung beim Flip bedeutet
FWHM Mitte / Rand5 %Nachfokussieren, Temperatur oder rutschender Okularauszug
Verkippung (als Vektor: Größe und Richtung)3 PunkteDurchhang oder Spiel in Kamera, Okularauszug, Adapter
Rand-Elongation eps_rad0,02Spiel in Korrektor oder Auszugsrohr, oder ein anderer Fokus
Nachführ-Elongation0,02Balance (ost-/westlastig), DEC-Spiel, Kabelzug
Koma k-Korrektor oder Abstand mit Spiel
Komafreier Punkt (Vektor)10 % von RMirror Flop, lockerer Fangspiegel oder Korrektor
Plate-Solve-Verzeichnung (Kapitel 12)1 pxein Teil der Optik verschiebt sich - oder eine Konvention des Solvers

Mit weniger als 3 Aufnahmen auf einer Seite gibt es kein Urteil. Vektoren (Verkippung, komafreier Punkt, Verzeichnung) werden als Vektoren verglichen: Ihre Richtung kann sich drehen, während ihre Länge gleich bleibt.

Drift im Laufe der Nacht

Temperatur und Höhe ändern sich langsam, und die Aufnahmen einer Seite folgen zeitlich aufeinander - eine Drift kann deshalb wie ein Sprung beim Flip aussehen. Das Script zieht den Median jeder Seite ab und korreliert den Rest mit der Zeit. Bei mindestens 6 Aufnahmen, einer Korrelation |r| ≥ 0,6 und einer Gesamtänderung über dem praktischen Minimum wird der Wert als driftend gemeldet, mit einer Warnung, wo er eine Änderung beim Flip vorgetäuscht haben könnte. Unterscheidet sich die Fokuserposition (FOCPOS) zwischen den Seiten, vermerkt der Bericht, dass nachfokussiert wurde.

Im Dialog: Series

12.Plate-Solve-Verzeichnung (SIP)

Ein Plate Solve ermittelt, wohin am Himmel jedes Pixel zeigt. Eine reale Optik bildet den Himmel nicht perfekt auf ein ebenes Raster ab; die SIP-Konvention (Simple Imaging Polynomial) speichert die Abweichung als Polynome in den Pixelkoordinaten u, v relativ zu einem Referenzpixel:

u′ = u + Σ Ap,q · up vq,   v′ = v + Σ Bp,q · up vq,   2 ≤ p + q ≤ Ordnung
SIP-Verzeichnung
Abbildung 18 - Quadratische Verzeichnungsterme verbiegen das Pixelraster.

Das misst die Positionen der Sterne - unabhängig von ihren Formen, und ein Nebenprodukt jedes Solves. Die Serienanalyse liest die sechs quadratischen Terme (p + q = 2) aus dem Header. Sie hängen nicht vom Referenzpixel ab und lassen sich deshalb auch dann vergleichen, wenn der Rest der Lösung fehlt. Das Script drückt sie als Verschiebung der vier Bildecken in Pixeln aus und vergleicht diese zwischen den Seiten.

Was sie aussagen kann und was nicht

  • Kein Maß für die Sensorverkippung. Eine Verkippung θ wirkt wie eine leichte Perspektive und verschiebt einen Punkt im Abstand r nur um etwa tan θ · r² / fpx (fpx = Brennweite in Pixeln). Bei 750 mm und 4,29-µm-Pixeln verschiebt 0,1° eine Ecke um 0,07 px - die FWHM-Methode (Kapitel 7) ist weit empfindlicher. Der Bericht zeigt die Verkippung, die eine reine Perspektive für die gemessenen Terme bräuchte; ein Wert von mehreren Grad bedeutet, dass sie andere Ursachen haben.
  • Kein Maß für den Backfokus. Die symmetrische (kissen-/tonnenförmige) Verzeichnung eines Korrektors ist ein Term dritter Ordnung; viele Solver passen nur die zweite Ordnung an.
  • Ein guter Indikator für Bewegung. Ändern sich die quadratischen Terme beim Flip, verschiebt sich etwas im optischen Zug unter der Schwerkraft - typischerweise der Hauptspiegel eines Newtons oder ein Korrektor, der sich im Auszugsrohr bewegt. Wandert auch der komafreie Punkt, stützt das diesen Befund.

Eine Vorzeichenumkehr beim Flip kann auch daher kommen, dass der Solver in einem System rechnet, das sich mit dem Himmel dreht. Um sicherzugehen, lösen Sie pro Seite eine Aufnahme mit dem ImageSolver von PixInsight und vergleichen. Die Terme werden nur verglichen, wenn sich alle Lösungen auf gleich große Bilder beziehen.

Im Dialog: Series

13.Grenzen und gute Praxis

  • Nicht registrierte, unbeschnittene Aufnahmen verwenden. Jede Berechnung setzt voraus, dass die Bildmitte die optische Achse ist. Registrierung, Beschneiden oder Drehen zerstören das.
  • Kalibriert, aber unbearbeitet. Strecken, Dekonvolution, Rauschreduzierung und Sternreduzierung verändern die Sternprofile.
  • Kurze Belichtungen, gutes Seeing. Seeing und Nachführung fügen eine Unschärfe hinzu, die kleine optische Fehler verdeckt. Mehrere Aufnahmen geben mehr Sicherheit als eine.
  • Ein reiches Sternfeld - mindestens 150 vermessene Sterne, gut verteilt bis in die Ecken, wo die Fehler am größten sind.
  • Eine Aufnahme verrät nicht das Vorzeichen eines Defokus. Wo eine Korrektur eine Richtung hat (Verkippung, Abstand), in kleinen Schritten ändern und erneut messen.
  • Geometrische Optik. Die Umrechnungen von Unschärfe in Δz ignorieren Beugung und Seeing-Profile; sie sind Schätzungen als Hilfe für eine Justage, keine Messungen mit Toleranzangabe.

Im Dialog: Der Dialog, Setup › General