以3D FFT为中心

Question

我再次提出了一个关于FFT结果居中的问题。在这种情况下，在具有FFTW库的3D空间中。目前我实施了下图所示的概念。它显示了细分的频率空间。每个子立方体（1-8）代表整个空间的八分之一。目的是在所有8个子立方体相遇的2x2x2立方体的中心重新定位低频。在此重定位步骤之后，我只需将原点移动到立方体的中心。

目前它是在前向fft后直接在ifft之前完成的，但我不确定这种方法是否正确，因为当我应用过滤器时我得到了奇怪的结果。如果没有过滤，我会按预期获得OUTPUT = INPUT。

经过几个小时的搜索后，我找不到任何以fft结果为中心的例子。相应的2D方法似乎有效。是否有一种常见的方法来处理它，因为对于低通滤波，我需要中心的低频和边缘的高频（默认情况下不是fftw给出的）。

问候

狼

Answer 1

这在FFTW常见问题解答中有介绍，请参阅http://www.fftw.org/faq/section3.html#centerorigin。你需要做的是将每个元素乘以：（ - 1）^（i + j + k）其中i，j和k是元素索引。