FITS标题的问题

时间:2017-09-24 22:32:46

标签: python astropy fits

我也曾在天文学上问过这个问题。

我正在处理一些似乎有不完整标题标签的NEOWISE图像。我已经被SAO(因为我使用DS9来查看它)的建议我已经

“定义了部分PC矩阵,但它不符合FITS WCS标准。使用CDELT和PC关键字的组合,CDELT关键字包含缩放因子,PC关键字,旋转矩阵,应该归一化。如果您希望旋转矩阵也包含比例因子,请使用CD关键字“。

这是我需要解决的标题之一:

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

我已经设法修复了一些无效的标签,但除了这些小编辑之外,这与我的研究主管给我的文件没有任何关系。 FITS标准here给出了每个标签的描述,但我不确定如何实现SAO的建议。任何人都可以帮我修复标签吗?

1 个答案:

答案 0 :(得分:3)

在错误消息之后,如果使用PCi_j格式来描述旋转矩阵,则应将值标准化为1,并且在CDELT值中单独编码任何缩放因子。如果想要在旋转矩阵中包含缩放因子,则应使用CDi_j格式,而不使用任何CDELT值。在v4.0 FITS规范草案中,FITS标准中的相关部分位于第29页的eqs附近。 9-11:https://fits.gsfc.nasa.gov/standard40/fits_standard40draft1.pdf

存在两种可能的解决方案。首先,缩放可以转移到CDELT值,如下所示:

PC1_1 = -1.0
PC2_2 = 1.0
CDELT1 = 0.0333333333333
CDELT2 = 0.0333333333333

或者,可以删除CDELT值并将PC1_1和PC2_2重命名为CD1_1和CD2_2。

请注意,排除PC1_2和PC2_1是有效的,因此消息中注明“部分PC矩阵[is]已定义”,因为非对角线值默认为零(在上面链接的文档的第32页上定义)这似乎是预期的效果。