Probleme mit FITS-Headern

Ich arbeite mit einigen NEOWISE-Bildern, die unvollständige Header-Tags zu haben scheinen. Ich wurde vom SAO darauf hingewiesen (weil ich DS9 verwende, um sie anzuzeigen), dass ich habe

"eine partielle PC-Matrix definiert, die jedoch nicht dem FITS-WCS-Standard entspricht. Bei einer Kombination aus CDELT- und PC-Schlüsselwörtern enthalten die CDELT-Schlüsselwörter den Skalierungsfaktor, die PC-Schlüsselwörter die Rotationsmatrix, die normalisiert werden sollten. Wenn Sie möchten, dass die Rotationsmatrix auch den Skalierungsfaktor enthält, verwenden Sie CD-Schlüsselwörter".

Hier ist einer der Header, die ich reparieren muss:

SIMPLE  =                    T / conforms to FITS standard                      
BITPIX  =                  -32 / array data type                                
NAXIS   =                    2 / number of array dimensions                     
NAXIS1  =                 1459                                                  
NAXIS2  =                  903                                                  
WCSAXES =                    2 / Number of coordinate axes                      
CRPIX1  =                729.5 / Reference pixel for axis 1                     
CRPIX2  =                451.5 / Reference pixel for axis 2                     
PC1_1   =     -0.0333333333333 / Coordinate transformation matrix element       
PC2_2   =      0.0333333333333 / Coordinate transformation matrix element       
CDELT1  =                   1. / [deg] Coordinate increment at reference point  
CDELT2  =                   1. / [deg] Coordinate increment at reference point  
CUNIT1  = 'deg     '           / Units of coordinate increment and value        
CUNIT2  = 'deg     '           / Units of coordinate increment and value        
CTYPE1  = 'GLON-AIT'           / Projection for axis1                           
CTYPE2  = 'GLAT-AIT'           / Projection for axis1                           
CRVAL1  =                   0. / [deg] Coordinate value at reference point      
CRVAL2  =                   0. / [deg] Coordinate value at reference point      
LONPOLE =                   0. / [deg] Native longitude of celestial pole       
LATPOLE =                  90. / [deg] Native latitude of celestial pole        
RADESYS = 'ICRS    '           / Equatorial coordinate system                   
BAND    = 'W1      '           / Wise Band                                      
END

Ich habe es bereits geschafft, ein paar der ungültigen Tags zu reparieren, aber abgesehen von diesen geringfügigen Änderungen ist dies unverändert gegenüber der Datei, die mir mein Forschungsleiter gegeben hat. Der FITS-Standard hier enthält Beschreibungen der einzelnen Tags, aber ich bin mir nicht sicher, wie ich den Rat der SAO implementieren soll. Kann mir jemand helfen, die Tags zu reparieren?

Haben Sie Fits View (FV)? heasarc.gsfc.nasa.gov/ftools/fv
Ach, das wusste ich nicht! Ich habe es jetzt heruntergeladen, danke. Mein Hauptproblem ist, dass ich nicht sicher bin, wie ich die Vorschläge des SAO umsetzen soll. Ich habe Astropy verwendet, um die Header zu bearbeiten, das ist also nicht das Problem, obwohl FV so aussieht, als würde es das bequemer machen.

Antworten (2)

Ihnen fehlen die Elemente PC1_2 und PC2_1 Ihrer Pixeltransformationsmatrix. Wenn Ihre Koordinatentransformation eine Rotation ist, würden Sie verwenden

PC1_1 = cos(t)
PC1_2 = sin(t)
PC2_1 = -sin(t)
PC2_2 = cos(t)

für die PC-Matrix. Die Pixelskala (in Grad/Pixel, nach CUNITi) geht in CDELTi ein.

Wie bereits erwähnt, können Sie die PC- und CDELT-Schritte in einer einzigen CD-Matrix kombinieren.

Weitere Einzelheiten finden Sie hier und hier .

Mit Fits Fits View (FV) View (FV) können Sie den Header entweder mit Terminalbefehlen ändern oder die GUI verwenden, mit der Sie alle Schlüsselwörter und Werte visualisieren können.