python - 如何在新的拟合文件中将两个 FITS 表合并为一个表？

Question

我有两个拟合文件数据(file1.fits 和 file2.fits)。第一个 (file1.fits) 包含 80,700 行重要数据，另一个包含 140,000 行。它们都有相同的 header 。

$ python
>>> import pyfits 
>>> f1 = pyfits.open('file1.fits')
>>> f2 = pyfits.open('file2.fits')
>>> event1 = f1[1].data
>>> event2 = f2[1].data
>>> len(event1)
80700
>>> len(event2)
140000

如何将 file1.fits 和 file2.fits 合并到新的拟合文件 (newfile.fits) 中，其标题与旧文件相同，并且 newfile.fits 的总行数为 80,700+ 140,000 = 220,700 ？

Answer 1

我尝试过 astropy :

from astropy.table import Table, hstack

t1 = Table.read('file1.fits', format='fits')
t2 = Table.read('file2.fits', format='fits')
new = hstack([t1, t2])
new.write('combined.fits')

它似乎适用于 NASA 的样本。