Problemas con los encabezados FITS

Estoy trabajando con algunas imágenes de NEOWISE que parecen tener etiquetas de encabezado incompletas. SAO me ha informado (porque estoy usando DS9 para verlos) que tengo

"se definió una matriz PC parcial, pero no se ajusta al estándar FITS WCS. Con una combinación de palabras clave CDELT y PC, las palabras clave CDELT contienen el factor de escala, las palabras clave PC, la matriz de rotación, que debe normalizarse. Si desea que la matriz de rotación también contenga el factor de escala, use palabras clave de CD".

Aquí está uno de los encabezados que necesito arreglar:

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

Ya logré arreglar un par de etiquetas inválidas, pero aparte de estas ediciones menores, esto no se ha alterado en el archivo que me dio mi supervisor de investigación. El estándar FITS aquí brinda descripciones de cada una de las etiquetas, pero no estoy seguro de cómo implementar los consejos de SAO. ¿Alguien puede ayudarme a arreglar las etiquetas?

¿Tiene Fit View (FV)? heasarc.gsfc.nasa.gov/ftools/fv
¡Ay, no sabía nada de eso! Ya lo descargué, gracias. Mi problema principal es que no estoy seguro de cómo implementar las sugerencias de SAO. He estado usando astropy para editar los encabezados, así que ese no es el problema, aunque parece que FV lo hará más conveniente.

Respuestas (2)

Le faltan los elementos PC1_2 y PC2_1 de su matriz de transformación de píxeles. Si su transformación de coordenadas es una rotación, usaría

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

para la matriz de PC. La escala de píxeles (en grado/píxel, según CUNITi) entra en CDELTi.

Como se señaló, puede combinar los pasos de PC y CDELT en una sola matriz de CD.

Ver aquí y aquí para más detalles.

Usando Fits Fits View (FV) View (FV), puede modificar el encabezado mediante comandos de terminal o usar la GUI que le permite visualizar todas las palabras clave y valores.