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

Question

我有两个fits文件数据(file1.fits和file2.fits)。第一个文件(file1.fit)包含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.fit和file2.fit合并到新的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')

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