我有两个具有不同分辨率的拟合格式的天文图像。我想将较高分辨率的图像卷入较低分辨率的图像。我认为这些是要执行的以下步骤。
对较高分辨率的图像进行去卷积以获得原始信号。
使用新的卷积滤镜（点扩散函数）卷积原始信号以获得最终图像。
我试图避免使用python进行反卷积，因为我不确定这些方法对于射电天文图像是否可靠。
我想过做傅立叶和逆傅里叶变换以解决这个问题。这是我的代码。
import numpy as np
image=fits.open(directory+file1)
test=image[0].data #test is the image convolved with the old beam
#Using Convolution Theorem
test1=np.fft.ifft2(np.fft.fft2(test)*np.fft.fft2(new_beam)/np.fft.fft2(old_beam))
当我尝试绘制时，使用：
plt.imshow(test1, interpolation='none', origin='lower')
plt.xlabel('x [pixels]')
plt.ylabel('y [pixels]')
plt.colorbar()
plt.show()
我的图片看起来不像卷积图像。有人能指出我做错了什么吗？其次，我应该如何规范我的强度值。我知道，对于一个傅立叶变换，你应该做absF（x）/ len（x）。如果我必须执行上述操作，我将如何处理呢？