从表数据列表创建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中的数据重新组合成新的拟合文件？我尝试使用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')