（某人终于想起来要更新网站了）

Astropy 作为天文研究领域常用的软件包，其中包括了坐标转换、图像分析、FITS读取等多项工具。

IRSA（NASA/IPAC Infrared Science Archive）是天文工作中经常使用到的红外数据库，可以检索单源的数据也可以导入tbl文件进行批量检索。

鉴于本人科研经历，初入IRSA使用时会遇到tbl不会写导致Multi-object Source经常出错，本文基于本人本科毕业设计《活动星系核的中红外光变研究》总结一部分经验以供参考。😄

ran

放个蓝蓝镇楼（）

IRSA的tbl format

IPAC Table Format给出了tbl格式文件的编写规范。

1
2
3
4
5
6
7
\keyword = value                                                  Keywords (optional)
\ Comment                                                         Comments (optional)
|  column1  |  column2  | column3 |  column4 |    column5       | Column Names (required)
|   double  |   double  |   int   |   real   |     char         | Data Types (standard)
|   unit    |   unit    |   unit  |   unit   |     unit         | Data Units (optional)
|   null    |   null    |   null  |   null   |     null         | Null Values (optional)
  165.466279  -34.704730      5      11.27         K6Ve           Data Rows (1 required)

由于Multi-object Source需要用到天体坐标，在文章中的数据获取中我采取了如下表格样式，包含了源名称及其赤经赤纬。（需要注意的是由于不同观测设备的观测能力和适用波段，给出的坐标通常与真实值存在一定偏差，另外以WISE和ZTF为例，其分别是工作于红外波段和光学波段的设备，因此如果采用Fermi源则需要光学对应体坐标）

1
2
3
4
|source_id|    ra|   dec|
|     char|double|double|
|         |   deg|   deg|
|     null|  null|  null|

IRSA需要的编码样式是ASCII，Astropy已经集成了包含编码在内的多项功能，只需要导入ASCII，里面包括了ipac样式，直接使用即可，见样例。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
import pandas as pd
from astropy.table import Table
from astropy.coordinates import SkyCoord
import astropy.units as u
from astropy.io import ascii

sources = []
ra = []
dec = []     #此处可以自行定义批量输入你需要批量检索的源的名称及其坐标

tbl = Table(
    [sources, ra, dec],
    names=["source_id", "ra", "dec"]
)

tbl["ra"].unit = "deg"
tbl["dec"].unit = "deg"

ascii.write(
    tbl,
    "name.tbl",
    format="ipac",   #Astropy已经含有ipac样式直接引用即可
    overwrite=True
)