python中如何把两个fits文件中的二维数据分别存入两个points

可以使用astropy库中的fits模块读取fits文件中的数据，并将数据转换为numpy数组，然后将数组中的二维数据分别存入两个points中。

示例代码：

from astropy.io import fits
import numpy as np
from astropy import coordinates as coord
from astropy import units as u

# 读取fits文件中的数据
data1, header1 = fits.getdata('file1.fits', header=True)
data2, header2 = fits.getdata('file2.fits', header=True)

# 将数据转换为numpy数组
data1 = np.array(data1)
data2 = np.array(data2)

# 获取数据的形状
shape1 = data1.shape
shape2 = data2.shape

# 将二维数据分别存入两个points中
points1 = coord.SkyCoord.from_pixel(np.arange(shape1[1]), np.arange(shape1[0]), header1, origin=0)
points2 = coord.SkyCoord.from_pixel(np.arange(shape2[1]), np.arange(shape2[0]), header2, origin=0)

在上述示例代码中，我们首先使用fits模块读取fits文件中的数据，并将数据转换为numpy数组。然后，我们获取了数据的形状，以便将二维数据分别存入两个points中。最后，我们使用astropy库中的coordinates模块创建了两个points对象points1和points2，分别存储了读取的fits文件的二维数据