Ускорить расчет IFFT в Python
Я работаю над программой 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 секунд. Есть ли способ ускорить это?