我正在开发一个 Python 程序(在 Enthought Canopy 中),它将计算 2 维和 3 维图像的两点自相关函数。
计算的核心是以下代码:
fftx = fftn(newimpad)
del newimpad
ret_int = (fftx * np.conjugate(fftx))
del fftx
ret = ifftn(ret_int)
ret = fftshift(ret)
newimpad
填充图像在哪里。
这行得通。然而,反向变换(ifftn
步骤)非常慢——对于我一直在处理的 2D 大图像,它需要 90 多秒。有什么办法可以加快这个速度吗?