提取第三個軸，數據立方體pyfits

Question

我有一個數據立方體fitsfile，包含一個規則（512x512）2d網格的強度。然而，還有第三軸，該軸上的每個值都對應一個頻率。在每個頻率下，拍攝這樣的2d照片。

我一直在試圖找到適合文件中的這些頻率值，但到目前爲止我還沒有運氣。我知道有140個不同的頻率通道，我通過查看適配文件的頭文件來了解單元等。 問題是找到那些140個通道的值。

順便編程在PYTHON中。

Answer 1

它通常以相同的方式給出軸「n」的座標。參考值「CRVALn」，參考位置「CRPIXn」，每個像素/通道「CDELTn」的增量以及在這些關鍵字「CTYPEn」中給出的座標類型。所以從這一點，你可以在你的立方體剛剛制定出每個「通道」的值作爲

(channel_number - CRPIXn) * CRDELTn + CRVALn 

既然你說你解決它，也許你可以發佈你對他人的利益的解決方案？

---

編輯/添加：如何在這個表達

更多信息FITS文件： 「三天打魚兩天曬譜座標三個代表」 雲英巖等 http://fits.gsfc.nasa.gov/fits_wcs.html

尤其是紙。 （2005）（link to NASA ADS）