2017-09-24 99 views
-1

我問這對astronomy.se爲好。與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的建議。任何人都可以幫我修復標籤嗎?

+0

你必須,如果你發佈你的代碼的一小部分工作示例得到答案的一個更大的概率。 – Gabriel

回答

3

在錯誤消息中,如果PCi_j格式是用來描述旋轉矩陣,該值應被歸一化到1,和任何比例因子在CDELT值分開編碼。如果想要在旋轉矩陣中包含縮放因子,則應使用CDi_j格式,而不使用任何CDELT值。在v4.0 FITS規範草案中,FITS標準中的相關部分大約在第29頁上。存在https://fits.gsfc.nasa.gov/standard40/fits_standard40draft1.pdf

兩個可能的解決方案:9-11。首先,縮放可以被轉移到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矩陣[已定義]」,因爲非對角線值默認爲零(在上面鏈接的文檔的第32頁定義)這看起來是期望的效果。