复制FITS头

Question

比如说我读了一份文件

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

来自astropy.io.fits的astropy.io.fits类有一个复制功能。

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