复制FITS标头

Question

说我读了一个FITS文件

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

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

我想复制其标题，以便可以在不更改原始标题的情况下修改副本

header_copy = cube.header
header_copy.remove('OBJCTY')

但是，这也会修改cube.header。

我该如何制作标题的实际副本，而不是创建指向标题的新指针？

header_copy = ?

Answer 1

Header中的astropy.io.fits类具有copy function。

header_copy = cube.header.copy()
header_copy.remove('OBJCTY')

然后header_copy是对象的实际副本，而不是指向cube.header的指针

if cube.header['OBJECTY'] != header_copy['OBJECTY']:
    print('All ok')

>>> All ok