天文数据处理Python库Skill astropy

---

name: astropy description: 全面的Python库，用于天文学和天体物理学。当处理天文数据时使用此技能，包括天体坐标、物理单位、FITS文件、宇宙学计算、时间系统、表格、世界坐标系统（WCS）和天文数据分析。当任务涉及坐标转换、单位转换、FITS文件操作、宇宙学距离计算、时间尺度转换或天文数据处理时使用。

Astropy

概述

Astropy是用于天文学的核心Python包，提供天文研究和数据分析的基本功能。使用astropy进行坐标转换、单位和数量计算、FITS文件操作、宇宙学计算、精确时间处理、表格数据操作和天文图像处理。

何时使用此技能

当任务涉及以下内容时使用astropy：

• 在天体坐标系统之间转换（ICRS、Galactic、FK5、AltAz等）
• 处理物理单位和数量（将Jy转换为mJy、秒差距转换为千米等）
• 读取、写入或操作FITS文件（图像或表格）
• 宇宙学计算（光度距离、回溯时间、哈勃参数）
• 使用不同时间尺度（UTC、TAI、TT、TDB）和格式（JD、MJD、ISO）进行精确时间处理
• 表格操作（读取目录、交叉匹配、筛选、连接）
• WCS在像素和世界坐标之间的转换
• 天文常数和计算

快速开始

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)

核心功能

1. 单位和数量（astropy.units）

处理带有单位的物理数量，执行单位转换，并确保计算中的维度一致性。

关键操作：

• 通过将值与单位相乘创建数量
• 使用.to()方法在单位之间转换
• 执行自动单位处理的算术运算
• 使用等效性进行领域特定转换（光谱、多普勒、视差）
• 处理对数单位（星等、分贝）

参见： references/units.md 获取全面文档、单位系统、等效性、性能优化和单位算术。

2. 坐标系统（astropy.coordinates）

表示天体位置并在不同坐标框架之间转换。

关键操作：

• 在任何框架中创建坐标（ICRS、Galactic、FK5、AltAz等）
• 在坐标系统之间转换
• 计算角分离和位置角
• 将坐标与目录匹配
• 包含距离以进行3D坐标操作
• 处理自行运动和径向速度
• 从在线数据库查询命名对象

参见： references/coordinates.md 获取详细坐标框架描述、转换、观察者依赖框架（AltAz）、目录匹配和性能提示。

3. 宇宙学计算（astropy.cosmology）

使用标准宇宙学模型执行宇宙学计算。

关键操作：

• 使用内置宇宙学（Planck18、WMAP9等）
• 创建自定义宇宙学模型
• 计算距离（光度、共动、角直径）
• 计算年龄和回溯时间
• 确定任何红移的哈勃参数
• 计算密度参数和体积
• 执行逆计算（为给定距离查找红移）

参见： references/cosmology.md 获取可用模型、距离计算、时间计算、密度参数和中微子效应。

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

读取、写入和操作FITS（灵活图像传输系统）文件。

关键操作：

• 使用上下文管理器打开FITS文件
• 通过索引或名称访问HDU（头数据单元）
• 读取和修改头（关键字、注释、历史）
• 处理图像数据（NumPy数组）
• 处理表格数据（二进制和ASCII表格）
• 创建新的FITS文件（单扩展或多扩展）
• 使用内存映射处理大文件
• 访问远程FITS文件（S3、HTTP）

参见： references/fits.md 获取全面文件操作、头操作、图像和表格处理、多扩展文件和性能考虑。

5. 表格操作（astropy.table）

处理表格数据，支持单位、元数据和各种文件格式。

关键操作：

• 从数组、列表或字典创建表格
• 以多种格式读写表格（FITS、CSV、HDF5、VOTable）
• 访问和修改列和行
• 排序、筛选和索引表格
• 执行数据库风格操作（连接、分组、聚合）
• 堆叠和连接表格
• 处理带有单位的列（QTable）
• 使用掩码处理缺失数据

参见： references/tables.md 获取表格创建、I/O操作、数据操作、排序、筛选、连接、分组和性能提示。

6. 时间处理（astropy.time）

精确时间表示和在时间尺度与格式之间转换。

关键操作：

• 以各种格式创建时间对象（ISO、JD、MJD、Unix等）
• 在时间尺度之间转换（UTC、TAI、TT、TDB等）
• 使用TimeDelta执行时间算术
• 为观察者计算恒星时
• 计算光行时校正（质心、日心）
• 高效处理时间数组
• 处理掩码（缺失）时间

参见： references/time.md 获取时间格式、时间尺度、转换、算术、观察特性和精度处理。

7. 世界坐标系统（astropy.wcs）

在图像中的像素坐标和世界坐标之间转换。

关键操作：

• 从FITS头读取WCS
• 将像素坐标转换为世界坐标（反之亦然）
• 计算图像足迹
• 访问WCS参数（参考像素、投影、尺度）
• 创建自定义WCS对象

参见： references/wcs_and_other_modules.md 获取WCS操作和转换。

额外功能

references/wcs_and_other_modules.md 文件还涵盖：

NDData和CCDData

用于n维数据集的容器，带有元数据、不确定性、掩码和WCS信息。

建模

用于创建和拟合数学模型到天文数据的框架。

可视化

用于天文图像显示的工具，具有适当的拉伸和缩放。

常数

具有适当单位的物理和天文常数（光速、太阳质量、普朗克常数等）。

卷积

用于平滑和滤波的图像处理内核。

统计

稳健的统计函数，包括sigma裁剪和异常值拒绝。

安装

# 安装astropy
uv pip install astropy

# 安装可选依赖以获取完整功能
uv pip install astropy[all]

常见工作流

在系统之间转换坐标

from astropy.coordinates import SkyCoord
import astropy.units as u

# 创建坐标
c = SkyCoord(ra='05h23m34.5s', dec='-69d45m22s', frame='icrs')

# 转换为银河坐标
c_gal = c.galactic
print(f"l={c_gal.l.deg}, b={c_gal.b.deg}")

# 转换为高度-方位角坐标（需要时间和位置）
from astropy.time import Time
from astropy.coordinates import EarthLocation, AltAz

observing_time = Time('2023-06-15 23:00:00')
observing_location = EarthLocation(lat=40*u.deg, lon=-120*u.deg)
aa_frame = AltAz(obstime=observing_time, location=observing_location)
c_altaz = c.transform_to(aa_frame)
print(f"Alt={c_altaz.alt.deg}, Az={c_altaz.az.deg}")

读取和分析FITS文件

from astropy.io import fits
import numpy as np

# 打开FITS文件
with fits.open('observation.fits') as hdul:
    # 显示结构
    hdul.info()

    # 获取图像数据和头
    data = hdul[1].data
    header = hdul[1].header

    # 访问头值
    exptime = header['EXPTIME']
    filter_name = header['FILTER']

    # 分析数据
    mean = np.mean(data)
    median = np.median(data)
    print(f"Mean: {mean}, Median: {median}")

宇宙学距离计算

from astropy.cosmology import Planck18
import astropy.units as u
import numpy as np

# 计算在z=1.5的距离
z = 1.5
d_L = Planck18.luminosity_distance(z)
d_A = Planck18.angular_diameter_distance(z)

print(f"Luminosity distance: {d_L}")
print(f"Angular diameter distance: {d_A}")

# 在该红移下的宇宙年龄
age = Planck18.age(z)
print(f"Age at z={z}: {age.to(u.Gyr)}")

# 回溯时间
t_lookback = Planck18.lookback_time(z)
print(f"Lookback time: {t_lookback.to(u.Gyr)}")

交叉匹配目录

from astropy.table import Table
from astropy.coordinates import SkyCoord, match_coordinates_sky
import astropy.units as u

# 读取目录
cat1 = Table.read('catalog1.fits')
cat2 = Table.read('catalog2.fits')

# 创建坐标对象
coords1 = SkyCoord(ra=cat1['RA']*u.degree, dec=cat1['DEC']*u.degree)
coords2 = SkyCoord(ra=cat2['RA']*u.degree, dec=cat2['DEC']*u.degree)

# 查找匹配
idx, sep, _ = coords1.match_to_catalog_sky(coords2)

# 按分离阈值筛选
max_sep = 1 * u.arcsec
matches = sep < max_sep

# 创建匹配目录
cat1_matched = cat1[matches]
cat2_matched = cat2[idx[matches]]
print(f"Found {len(cat1_matched)} matches")

最佳实践

1. 始终使用单位：为单位附加数量以避免错误并确保维度一致性
2. 使用上下文管理器处理FITS文件：确保文件正确关闭
3. 优先使用数组而非循环：处理多个坐标/时间作为数组以提高性能
4. 检查坐标框架：在转换前验证框架
5. 使用适当的宇宙学：为分析选择正确的宇宙学模型
6. 处理缺失数据：对包含缺失值的表格使用掩码列
7. 指定时间尺度：为精确计时明确时间尺度（UTC、TT、TDB）
8. 使用QTable处理带单位的表格：当表格列带有单位时
9. 检查WCS有效性：在使用转换前验证WCS
10. 缓存常用值：可以缓存昂贵计算（例如，宇宙学距离）

文档和资源

• 官方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、建模、可视化、常数和实用程序