使用Astropy解析FITS文件WCS报错，求助排查问题根源

错误根源分析 #

从报错信息来看，核心问题有两点：

1. FITS头中存在DSS类型的畸变配置，但轴1的系数尺度被设为0，导致WCSLIB无法初始化畸变函数；
2. 你提供的WCS头信息缺少必需的CTYPE1和CTYPE2关键字——这两个参数定义了坐标系统的投影类型（比如天球坐标的切向投影RA---TAN/DEC--TAN），是WCS解析的核心依赖，仅靠CRVAL和CDELT无法构建完整的WCS。

解决方法 #

1. 检查完整的FITS头信息 #

先确认FITS头中所有与畸变相关的关键字（比如DISTORTYPE、DSVAL*、DS*开头的项），你提供的片段里没有显示这些，但报错明确指向DSS畸变：

from astropy.io import fits
hdulist = fits.open('abell.fits')
# 打印全部头信息，查找畸变相关关键字
print(hdulist[0].header[:])

如果找到错误的DSS畸变配置，直接删除这些关键字即可。

2. 让Astropy自动修复WCS问题 #

初始化WCS时添加relax=True（放宽标准检查）和fix=True（自动修复参数错误），可以绕过无效的畸变配置并补全必要参数：

from astropy import wcs
from astropy.io import fits

hdulist = fits.open('abell.fits')
# 启用自动修复和放宽检查
w = wcs.WCS(hdulist[0].header, relax=True, fix=True)

print(w.wcs.name)
w.wcs.print_contents()

3. 手动构建基础线性WCS #

如果你的数据不需要畸变校正，直接手动补全必需的WCS关键字并移除畸变配置：

from astropy import wcs
from astropy.io import fits

hdulist = fits.open('abell.fits')
header = hdulist[0].header.copy()

# 添加天球坐标投影的CTYPE关键字（根据实际投影类型调整）
header['CTYPE1'] = 'RA---TAN'
header['CTYPE2'] = 'DEC--TAN'

# 删除所有DSS畸变相关关键字
for key in list(header.keys()):
    if key.startswith('DS') or key == 'DISTORTYPE':
        del header[key]

# 初始化WCS
w = wcs.WCS(header)

print(w.wcs.name)
w.wcs.print_contents()

内容的提问来源于stack exchange，提问作者Adrian