我有一个天空区域的2个不同的拟合文件，fitsFile2是适合文件1所涵盖的一个小区域。我想在同一个绘图中绘制它们并将fitsFile2叠加在fitsFile1上的正确位置，但结果是基于它们的像素数而不是坐标。我已经使用WCS和aplpy进行坐标，到目前为止还没有运气。我得到这张图片：2 fits files
蓝色地图应该是另一个地图的一小部分，但因为它有更多的像素数，即使它覆盖了天空的一小部分，它也会在地图中变得更大。
我该如何解决这个坐标问题？感谢。
这是我尝试过的一个例子（使用wcs）：
hdulist2 = fits.open(FitsFile2Path)
hdu2 = hdulist2[0]
wcs2 = WCS(hdu2.header)
hdulist1 = fits.open(FitsFile1Path)
hdu1 = IRhdulist1[0]
wcs1 = WCS(hdu1.header)
fig = plt.figure(num=None, figsize=(9,6), dpi=120, facecolor= 'w', edgecolor ='k')
ax = fig.add_subplot(111, projection=wcs2)
ax.set_xlim(-0.5, hdu1.data.shape[1] - 0.5)
ax.set_ylim(-0.5, hdu1.data.shape[0] - 0.5)
ra = ax.coords[0]
dec = ax.coords[1]
ra.set_axislabel('RA (degree)')
dec.set_axislabel('Dec (degree)')
ra.grid(color='black', alpha=0.5, linestyle='solid')
dec.grid(color='black', alpha=0.5, linestyle='solid')
ra.display_minor_ticks(True)
dec.display_minor_ticks(True)
ra.set_minor_frequency(10)
plt.imshow(hdu2.data, origin='lower', cmap = 'BrBG', interpolation='nearest')
plt.imshow(hdu1.data, origin='lower')
plt.title('Test' , fontsize = 12)
plt.savefig(saveFilePath)