如何在新的拟合文件中将两个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合并到新的文件（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')

它似乎适用于美国宇航局的样本。