切片RGB适配imag时保留WCS坐标

Question

我有一个RGB拟合图像的集合，我想用坐标轴上的WCS绘制蓝色的框架。我用以下方法将图像切成蓝色：

from astropy.wcs import WCS
from astropy.io import fits

image1 = fits.open('fits_file.fits')
data1 = image1[0].data
image1_slice_blue = data1[2,:,:]

当我绘制这个时，轴对应于像素的数量，而不是WCS坐标。 我尝试使用以下方法从文件中获取WCS信息：

wcs1 = WCS(image1[0].header)

但是我有一个问题，因为NAXIS=3（图像宽度、高度和平面），但是WCSAXS=2（RA和Dec）。错误显示在核心WCS关键字中检测到3个维度，并告诉我使用naxis-kwarg选择/减少维度。我又运行了一次，但是：

wcs1 = WCS(image1[0].header, naxis=2)

但这也不起作用，因为找不到CTYPE3。有人能帮我解决这个问题吗？谢谢

这是我试图分割的图像：http://nova.astrometry.net/user_images/2677994#annotated

我无法共享实际文件，但它是“新文件”-图像.fits“在RHS上归档。你知道吗