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

Question

我有两个适合文件数据（file1.fits 和 file2.fits）。第一个 (file1.fits) 包含 80,700 行重要数据，另一个是 140,000 行。他们都有相同的标题。

$ 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 合并到新的 fit 文件 (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 的样本。