python - 从表数据列表创建 FITS 文件

Question

我想我在从 fit 文件中选择数据时搞砸了……基本上我所做的是：

File='/SomePath/xxx.fits'
hdulist=fits.open(File)
tbdata=hdulist[1].data

并应用了选择规则：

for i in range (0, len(tbdata)):
     if tbdata[i]['z']<0.2:
         A.append(tbdata[i])

有没有办法将 A 中的数据重新组合成一个新的 fit 文件？我尝试使用 Pandas 进行操作，但找不到正确的方法...

Answer 1

大概是这样的：

from astropy.table import Table
tbdata = Table.read('file.fits')
ok = tbdata['z'] < 0.2  # boolean selection mask
new_tbdata = tbdata[ok]
new_tbdata.write('new_file.fits')