Question

我有这个脚本生成一个带有多个扩展名的文件，但我想将旧文件中的标题添加到扩展名中。

new_hdul = fits.HDUList()
new_hdul.append(fits.PrimaryHDU(header=headermain))
new_hdul.append(fits.ImageHDU(nod1, header=header1, name='Chop1')) 
new_hdul.append(fits.ImageHDU(nod2, header=header2, name='Chop2'))
new_hdul.append(fits.ImageHDU(diff1, name='Dif'))

现在我试过了：

headermain = fits.getheader(file,0)

和

headermain = fits.open(file).header.copy()

但两个都给我错误说

ValueError：header必须是Header对象

我该如何解决这个问题？

headermain = fits.getheader(file,0)
print(headermain)

请参阅http://pastebin.com/JXki7EPV

Answer 1

通常从文件中将标题作为Header对象获取并不复杂。您astropy.io.fits.open()该文件并从PrimaryHDU中提取标题：

from astropy.io import fits

filename = 'test.fits'

with fits.open(filename) as hdus:
    headermain = hdus[0].header

或getheader：

headermain = fits.getheader(filename) # Defaults to primary header!

，结果将是fits.Header - 您可以在写作期间使用的对象。

但是如果您的文件不是有效的FITS文件，则可能存在问题。如果这不起作用，您可以编辑您的问题并显示这两个功能中的任何一个的输出吗？

print(headermain)

将标题复制到新文件astropy中

1 个答案: