我正在研究一个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秒。有什么方法可以加快速度吗?