在python中，如何将ascii表转换为fits表？

Question

我有一个数据表在一个ascii文件头，我想把我的ascii表转换成一个有标题的fits文件

#ID ra  dec x   y   Umag    Bmag    Vmag    Rmag    Imag    
1.0 53.146  -27.8123    3422.98 3823.58 24.4528 24.7995 23.6266 22.64   20.8437 
2.0 53.1064 -27.801         3953.49 3994.62 23.3284 22.6716 22.1762 21.6189 21.2141 
3.0 53.1322 -27.7829    3608.34 4269.29 21.2676 20.1937 19.6743 19.0707 18.6983 
4.0 53.1017 -27.8022    4017.09 3975.24 23.6987 22.84   21.9946 21.0781 19.8616 
5.0 53.118  -27.8021    3798.98 3978.42 23.3087 22.1932 21.2205 20.1842 18.6448     
6.0 53.1479 -27.8239    3397.92 3648.27 25.0347 24.598  23.7259 22.9945 21.9228     
7.0 53.1334 -27.7758    3592.51 4375.76 21.5159 20.4777 19.6065 18.6609 17.188      
8.0 53.1048 -27.8259    3974.47 3617.5  22.3266 22.3517 22.0677 21.7664 21.6781     
9.0 53.1249 -27.8109    3706.47 3843.89 24.0539 23.3009 22.4001 21.4732 20.1244     
10.0 53.1207 -27.7822   3763.3  4278.76 24.417  23.7635 22.9405 22.1379 21.5564     
11.0 53.0886 -27.7611   4193.25 4596.77 22.012  22.3081 22.125  21.9488 21.9071     
12.0 53.1272 -27.7498   3676.7  4768.82 19.3631 19.7458 19.5979 19.4438 19.4002 

你知道我怎么用python做到这一点吗？ 干杯。在

Answer 1

将ASCII文件转换为FITS文件可以以一种更简单、更省时的方式完成。 我们可以利用这一点：

from astropy.io import ascii
test_input = '../Data/test.txt'
text_file = ascii.read(test_input)
text_file.write('text_to_fits.fits')

如果没有提供绝对路径，则将在当前工作目录中创建FITS文件。在

所需安装：（Astropy库）

^{pr2}$

您可以在下面的链接中查看其他基本要求。例如，纽比。 如果不存在的话，可以下载上面的语句。所以不用担心。在

Astropy网站链接：http://docs.astropy.org/en/stable/install.html

创建配合文件后，可以按如下方式打开它：

import astropy.table as at
fits_file = '../Data/text_to_fits.fits'
table = at.Table.read(fits_file)

为了操作和使用FITS文件，可以尝试使用Astropy。 它最初是为天文学相关数据集编写的，在那里FITS文件非常常用，当涉及到其他应用程序时，它也是一个非常有用的库。在

文档：http://docs.astropy.org/en/stable/index.html

从输入到输出：

输入ASCII文件：

col1 col2 col3
1234 2345 3456

输出适合文件：

Out[10]: 
<Table length=1>
 col1  col2  col3
int64 int64 int64
----- ----- -----
    1     2     3

希望这有帮助！在

Answer 2

尽管问题是如何通过Python实现，但更简单的方法是使用Topcat（目录和表操作工具）。在

此软件的网站可在THIS LINK HERE中找到

Topcat允许您加载几乎所有格式（ASCII、FITS、FITS+、CSV等）。一旦加载了表格，在本例中是ASCII表格，您可以将其保存为FITS文件，就像将.pptx文档保存到.pdf文档一样！保存为FITS文件后，就可以打开FITS文件，并且可以手动输入列的标题！

这是非常快速和有用的，而且可以立即绘制列，而无需编写任何代码。

Answer 3

我找到了解决自己问题的方法： 将文件作为数组读取

import pyfits
from scipy.io import *
M=read_array(filename)

将每个列指定为标题名称

^{pr2}$

将标题和列写入fits文件：

tbhdu = pyfits.new_table(cols)
hdu = pyfits.PrimaryHDU(data=M)
thdulist = pyfits.HDUList([hdu,tbhdu])
thdulist.writeto(outfilename)
thdulist.close()

成功了！！在