jjzjj

关于 python:FFTW3 上的复杂 numpy 数组直接在 scipy.weave.inline

FFTW3oncomplexnumpyarraydirectlyinscipy.weave.inline我正在尝试在Python中实现基于FFT的亚像素移位(平移)算法。傅里叶位移定理允许将阵列平移一个子像素量:1.正向FFT阵列2.在傅里叶空间中将阵列乘以线性相位斜坡3.逆FFT数组这个算法很容易在python中使用numpy/scipy实现,但是对于256**2数组来说,每班的速度非常慢(~10毫秒)。我试图通过使用scipy.weave.inline直接从python调用c代码来加快速度。我在将复杂的numpy数组传递给FFTW时遇到了麻烦。c代码如下所示:123456789101112
12