Konvertieren von FITS (NuSTAR)-Koordinaten in WCS

Ich versuche, einige Pixelkoordinaten, die ich habe, in WCS-Koordinaten umzuwandeln, idealerweise in eine WCS-Region, um sie in weiteren Analysen zu verwenden.

Bisher konnte ich eine NuSTAR FITS-Datei laden und analysieren, einige Analysen durchführen, um meine Auswahl zu treffen und die Pixelkoordinaten für das Bild zu erhalten. Zum Beispiel würde ich nach meiner Analyse am Ende ein Rechteck mit xWerten von Pixel 480 bis 518 und yWerten von 478 bis 516 auswählen.

Danach habe ich versucht, sie mit dem Astropy WCS-Modul zu konvertieren, aber es scheint nicht die erforderlichen Daten im Header zu finden, um die Konvertierung tatsächlich durchzuführen, und sagt nur, dass Pixel 480 die Koordinate 480 ergibt. Das oder ich mache es gerade Stimmt etwas nicht.

Also habe ich den FITS-Header selbst durchgesehen und festgestellt, dass einige Schlüssel ( TCRPX, TCRVL, TCDLT) das Referenzpixel, die Gradkoordinate des Referenzpixels und die Skalierungswerte der Pixelachse angeben.

Aber jetzt stecke ich wirklich fest, da ich eine Regionsdatei wie die von ds9 erzeugte im Format von:

# Region file format: DS9 version 4.1
global color=green dashlist=8 3 width=1 font="helvetica 10 normal roman" select=1 highlite=1 dash=0 fixed=0 edit=1 move=1 delete=1 include=1 source=1
fk5
circle(19:58:21.593,+35:11:57.46,10")

Um diese dann durch das nuproductsTool laufen zu lassen.

TL; DR: Weiß jemand, wie man die physikalischen Pixelkoordinaten von NuSTAR FITS in FK5-Koordinaten umwandelt?

BEARBEITEN: Etwas Code geschrieben, um die Koordinaten zu extrahieren . Es rechnet die Pixel Pixel im Glossar erklärt grad fein um, konvertiert aber die FK5 a angle to sexagesimal funktioniert aus irgendeinem Grund nicht richtig.

EDIT2: Alles funktioniert jetzt, hier ist ein Notizbuch, in dem ich alles durchlaufe, falls jemand auf ein ähnliches Problem stößt. Ich extrahiere die Rektaszension und Deklination in Grad aus den Pixelkoordinaten, die in der FITS-Datei angegeben sind, und konvertiere diese Grade dann in Sexagesimalwinkel (h:m:s). Nach meinen Tests scheint es ziemlich gut zu funktionieren, die Ergebnisse sind ein bisschen abweichend, da ich annehme, dass es sich um einige Gleitkomma-Arithmetikfehler handelt.

In welcher Programmiersprache arbeitest du? Außerdem ist mir nicht ganz klar, was Sie als Ausgabe wollen, und ich bin mit diesem 'nuproducts'-Tool nicht vertraut. Warum sind die Grade, in die du die Pixelkoordinaten umgerechnet hast, wie du es im EDIT sagst, nicht ausreichend?
Ich arbeite jetzt in Julia. Grundsätzlich müssen NuSTAR-Daten kalibriert werden, bevor sie für die wissenschaftliche Analyse geeignet sind, und diese Kalibrierung umfasst einige Schritte. In einem der Schritte müssen Sie eine Quellregion auswählen (als .reg-Datei, die normalerweise durch manuelles Auswählen eines Bereichs mit erstellt wird DS9), bevor Sie die endgültigen, vollständig kalibrierten/gereinigten Daten erhalten. Die FITS-Dateien speichern Daten als XY-Pixelkoordinaten, aber die Regionsdatei muss in WCS sein, also konvertiere ich X, Y in FK5-Winkel. Das Problem ist jetzt, dass die Umwandlung des Winkels, den ich in Grad erhalte (was korrekt ist), in hh: mm: ss anders ist als das, was DS9 ausgibt
Wie groß ist die Diskrepanz zwischen der DS9-Ausgabe und dem, was Sie erhalten?
Aha, ich habe herausgefunden, was ich falsch gemacht habe, ich war dumm und habe Grad in h: m: s umgewandelt, als wären sie Winkel auf einem Kreis, was für Rektaszension nicht funktioniert, da es von -90 bis 90 geht, nicht 0 bis 360. Jetzt, wo ich es repariert habe ( hier ist ein Notebook, wo es funktioniert ), ist es immer noch ein bisschen weg von etwa 10 ^ -5, aber ich denke, das sind nur einige Gleitkommafehler, die sich auf dem Weg summieren.

Antworten (2)

Sie können Ihre FITS-Dateien auf astrometry.net hochladen (oder eine API verwenden) und die Koordinaten abrufen. Sie können optional neue FITS-Dateien mit den in den Metadaten enthaltenen Koordinaten zurückerhalten. Von dort aus können Sie mit Ihrer Bearbeitung fortfahren.

Und anstatt alles hochzuladen, können Sie auch eine Kopie des astrometry.net-Pakets lokal installieren, wenn Sie Linux verwenden.

Ich habe diese Übung mit IRAF mit den Funktionen geomap und geotran gemacht. Es gibt eine PYRAF-Version der IRAF, aber ich habe die Python-Version noch nie ausprobiert. Mit diesen Routinen können Sie die Transformationsgleichung unabhängig vom Header finden. Was Sie wissen müssen, sind einige Standardobjekte mit bekannten WCS- und Pixelkoordinaten.

Konvertieren von FITS (NuSTAR)-Koordinaten in WCS
AntwortenHierDatenschutz-Bestimmungen1y, 0v