Question

我正在尝试修改FITS图像并相应地修改其标题，然后将其保存到新的FITS文件中

from astropy.io import fits
from astropy.utils.data import get_pkg_data_filename

# Open FITS file
image_file = get_pkg_data_filename('tutorials/FITS-images/HorseHead.fits')
cube = fits.open(image_file)[0]

# Modify the cube
cube.header.remove('NAXIS2')

# ...

cube.header['NAXIS2'] = 5

# Save the new FITS file
hdu = fits.PrimaryHDU(cube.data)
hdu.header = cube.header
hdu.writeto('new.fits', overwrite=True)

但是，这将返回以下错误：

VerifyError: 
Verification reported errors:
HDU 0:
    'NAXIS2' card at the wrong place (card 159).
    'EXTEND' card at the wrong place (card 160).
Note: astropy.io.fits uses zero-based indexing.

如何将卡以正确的顺序放在页眉中？还是有办法直接在正确的位置放置卡？

注意：这显然不是修改标头中字段的最佳方法，但这是再现我所得到的错误的最短示例。我不想知道如何修改字段，而是想知道如何在标题的正确位置获取新字段。

Answer 1

Header类具有insert功能，可以在给定卡之前或之后添加卡。

cube.header.insert('EXTEND', ('NAXIS2', 5))

将在值'NAXIS2'之前添加值5的{{1}}卡。

通过使用'EXTEND'关键字，也可以在给定卡之后立即设置新卡。所以做同样事情的另一种方法是：

after=True

使用此修改，cube.header.insert('NAXIS1', ('NAXIS2', 5), after=True)不再引发错误。

重新排列FITS标头中的卡

1 个答案: