python3 numpy二维方法_使用Python numpy 进行二维傅里叶变换 和 图片自相关

本篇主要记录一下

使用Python numpy库 进行求解图像的FFT 以及求解图片的自相关函数

研究背景：

1)某M姓知名数学矩阵分析工具由于灯塔国的政策导致不许使用，平常在M中跑的函数需要转到Python中；

2)笔者的研究课题需要研究一个图片的自相关函数(Auto correlation)

所以使用Python实现了一下

废话结束下面是正文：

#FFT:

#使用 numpy.fft.fft2(正变换) 和numpy.fft.ifft2函数(反变换)

#与Mxxxxx一样，需要np.fft.fftshift 和np.fft.ifftshift调整(相关文章一大堆)

#(画图部分以及完整代码在下面)

f = np.fft.fft2(image)

fshift = np.fft.fftshift(f)

求图片的自相关：网上的资料较少，先给出公式：

中间的圆圈和叉的组合符号 是相关符号 使用傅里叶变换的卷积定理(相关定理) 不同的教科书叫法不同，找一本图像处理的书都有讲到 下面截取自超星阅读中 某一本数字图像处理的内容：

笔者只需要求自相关即可 所以上文中共轭相乘变成了求FT模的平方

下面是公式：

#Auto correlation：

AC = np.fft.ifft2(np.abs(f) ** 2 )

AC = np.fft.fftshift(AC)

完整代码如下&#