Astropy.io.fits库简介：在Python中处理FITS文件的基本方法

Astropy.io.fits是Astropy库中的一个子模块，用于在Python中处理FITS（Flexible Image Transport System）文件。FITS是一种广泛用于天文学的数据格式，包含了天文图像和表格数据。

Astropy.io.fits提供了一组用于读取、写入和操作FITS文件的函数和类。下面是使用Astropy.io.fits处理FITS文件的基本方法的简介。

1. 导入Astropy.io.fits库：

from astropy.io import fits

2. 打开FITS文件：

可以使用open()函数打开FITS文件，并将其保存在一个变量中，以供后续的读取和操作。

hdul = fits.open('filename.fits')

3. 查看FITS文件的信息：

可以使用info()函数来显示FITS文件的基本信息，包括扩展名、数据单元的数量、每个数据单元的类型等。

hdul.info()

4. 读取图像数据：

可以使用getdata()函数读取FITS文件中的图像数据，并保存在一个变量中。

image_data = hdul[0].data

这里的hdul[0]表示 个扩展，即图像数据所在的位置。

5. 读取表格数据：

可以使用getdata()函数读取FITS文件中的表格数据，并保存在一个变量中。

table_data = hdul[1].data

这里的hdul[1]表示第二个扩展，即表格数据所在的位置。

6. 修改图像数据：

可以直接对图像数据进行修改，例如修改某个像素的数值。

image_data[10, 10] = 100

这里将图像数据中第(10, 10)个像素的数值修改为100。

7. 创建新的FITS文件：

可以使用HDUList()和PrimaryHDU()类创建一个新的FITS文件，并将图像或表格数据添加到其中。

new_hdul = fits.HDUList()
new_hdul.append(fits.PrimaryHDU(image_data))
new_hdul.writeto('new_filename.fits')

这里将图像数据添加到一个新的HDUList对象中，并将该对象保存为一个新的FITS文件。

以上是使用Astropy.io.fits处理FITS文件的基本方法的简介。Astropy.io.fits还提供了其他更高级的功能，例如处理多个扩展和更新FITS文件的头文件信息。可以参考Astropy官方文档以获取更详细的信息和示例代码。