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

Question

我有两个fits文件数据（file1.fits和file2.fits）。第一个（file1.fits）包含80700个重要的数据行，另一个是140000行。它们都有相同的标题。在

$ 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合并到新的fits文件中(新文件.fits)具有与旧标题相同的标题以及新文件.fits是80700+140000=220700？在

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')

它似乎与美国宇航局的样本有关。在