astropy

用途

Astropy 是面向天文学与天体物理研究的核心 Python 包，为天文研究与数据分析提供基础功能。当任务涉及天体坐标转换、单位与量值计算、FITS 文件操作、宇宙学计算、精确时间处理、表格数据操作以及天文图像处理时，可使用该技能。

主要能力

单位与量值（astropy.units）

处理带单位的物理量，执行单位换算，并确保计算中的量纲一致性。主要操作包括：将数值与单位相乘创建量值；使用 .to() 方法进行单位转换；算术运算时自动处理单位；使用等价关系进行特定领域换算（光谱、多普勒、视差）；处理对数单位（星等、分贝）。

坐标系统（astropy.coordinates）

表示天体位置并在不同坐标框架之间转换。主要操作包括：使用 SkyCoord 在任意框架（ICRS、Galactic、FK5、AltAz 等）中创建坐标；在坐标系统之间转换；计算角距离与位置角；将坐标与星表匹配；包含距离以进行三维坐标操作；处理自行与径向速度；从在线数据库查询命名天体。

宇宙学计算（astropy.cosmology）

使用标准宇宙学模型执行宇宙学计算。主要操作包括：使用内置宇宙学模型（Planck18、WMAP9 等）；创建自定义宇宙学模型；计算距离（光度距离、共动距离、角直径距离）；计算年龄与回溯时间；确定任意红移处的哈勃参数；计算密度参数与体积；执行反向计算（根据给定距离求红移）。

FITS 文件处理（astropy.io.fits）

读取、写入和操作 FITS（Flexible Image Transport System）文件。主要操作包括：使用上下文管理器打开 FITS 文件；按索引或名称访问 HDU（Header Data Unit）；读取和修改头信息（关键字、注释、历史）；处理图像数据（NumPy 数组）；处理表格数据（二进制表和 ASCII 表）；创建新的 FITS 文件（单扩展或多扩展）；对大文件使用内存映射；访问远程 FITS 文件（S3、HTTP）。

表格操作（astropy.table）

处理支持单位、元数据和多种文件格式的表格数据。主要操作包括：从数组、列表或字典创建表格；以多种格式（FITS、CSV、HDF5、VOTable）读写表格；访问和修改列与行；排序、筛选和索引表格；执行数据库风格操作（连接、分组、聚合）；堆叠与拼接表格；处理带单位的列（QTable）；使用掩码处理缺失数据。

时间处理（astropy.time）

精确表示时间并在时间尺度与格式之间转换。主要操作包括：以多种格式（ISO、JD、MJD、Unix 等）创建 Time 对象；在时间尺度（UTC、TAI、TT、TDB 等）之间转换；使用 TimeDelta 进行时间算术；为观测者计算恒星时；计算光行时修正（质心、日心）；高效处理时间数组；处理掩码（缺失）时间。

世界坐标系统（astropy.wcs）

在图像像素坐标与世界坐标之间转换。主要操作包括：从 FITS 头读取 WCS；将像素坐标转换为世界坐标（反之亦然）；计算图像覆盖范围；访问 WCS 参数（参考像素、投影、比例）；创建自定义 WCS 对象。

其他能力

references/wcs_and_other_modules.md 还涵盖：NDData 与 CCDData（带元数据、不确定性、掩码和 WCS 信息的 N 维数据集容器）；建模（创建并拟合数学模型到天文数据的框架）；可视化（具有适当拉伸与缩放的天文图像显示工具）；常数（带正确单位的物理与天文常数，如光速、太阳质量、普朗克常数等）；卷积（用于平滑和滤波的图像处理核）；统计（包括 sigma 裁剪和异常值剔除的稳健统计函数）。

使用方式

安装命令：

# 安装 astropy
uv pip install astropy

# 安装完整功能所需的可选依赖
uv pip install astropy[all]

快速开始示例：

import astropy.units as u
from astropy.coordinates import SkyCoord
from astropy.time import Time
from astropy.io import fits
from astropy.table import Table
from astropy.cosmology import Planck18

# 单位与量值
distance = 100 * u.pc
distance_km = distance.to(u.km)

# 坐标
coord = SkyCoord(ra=10.5*u.degree, dec=41.2*u.degree, frame='icrs')
coord_galactic = coord.galactic

# 时间
t = Time('2023-01-15 12:30:00')
jd = t.jd  # 儒略日

# FITS 文件
data = fits.getdata('image.fits')
header = fits.getheader('image.fits')

# 表格
table = Table.read('catalog.fits')

# 宇宙学
d_L = Planck18.luminosity_distance(z=1.0)

常见工作流包括：坐标系统间转换（如 ICRS 转 Galactic、转 AltAz 需时间与地点）；读取与分析 FITS 文件（打开文件、查看结构、获取图像数据与头信息、访问头关键字、分析数据）；宇宙学距离计算（计算光度距离、角直径距离、宇宙年龄、回溯时间）；星表交叉匹配（读取星表、创建坐标对象、使用 match_to_catalog_sky 查找匹配、按分离阈值筛选、创建匹配后的星表）。

适用场景

适用于天文研究与数据分析中的坐标转换、单位换算、FITS 文件读写、宇宙学计算、时间处理、星表操作、WCS 变换等任务。

限制

资料未提供关于收费、价格、平台限制或具体性能指标的信息。

最佳实践

始终使用单位：为量值附加单位以避免错误并确保量纲一致性；对 FITS 文件使用上下文管理器：确保文件正确关闭；优先使用数组而非循环：以数组方式处理多个坐标/时间以获得更好性能；检查坐标框架：转换前验证框架；使用适当的宇宙学模型：为分析选择正确的宇宙学模型；处理缺失数据：对含缺失值的表格使用掩码列；指定时间尺度：对精确计时明确时间尺度（UTC、TT、TDB）；对带单位的表格使用 QTable：当表格列具有单位时；检查 WCS 有效性：使用变换前验证 WCS；缓存常用值：可缓存开销较大的计算（如宇宙学距离）。

文档与资源

官方 Astropy 文档：https://docs.astropy.org/en/stable/ 教程：https://learn.astropy.org/ GitHub：https://github.com/astropy/astropy

参考文件：references/units.md（单位、量值、换算与等价关系）；references/coordinates.md（坐标系统、转换与星表匹配）；references/cosmology.md（宇宙学模型与计算）；references/fits.md（FITS 文件操作与处理）；references/tables.md（表格创建、I/O 与操作）；references/time.md（时间格式、尺度与计算）；references/wcs_and_other_modules.md（WCS、NDData、建模、可视化、常数与工具）。