c - 使用 fftw 的卷积尺度

Question

我编写了一个简单的程序（在 ANSI C 中）将两个真正的二维数组与 FFTW 进行卷积（使用 r2c 和 c2r 计划）。一切（除了我对比例的理解）都运行良好，输出看起来也不错。

问题是我对如何缩放两个数组有点困惑，所以在处理之后IDFT[DFT(A)DFT(B)]（DFT：离散傅里叶变换。I：逆。）最终结果的比例是正确的。假设图像 A 是宽度a像素的正方形，卷积核是b像素。然后我创建两个宽度的填充方形数组，p=a+b并将我的A和B图像都放入其中并进行卷积。

将最终输出中的所有像素乘以1/p^2正确吗？提前致谢。

Answer 1

既然你想要它作为答案：

http://fftw.org/doc/What-FFTW-Really-Computes.html （我认为网址足够明确）。

一般来说，我发现 FFTW 的文档非常有用，但其中包含您需要了解的一些小细节。例如，在关闭程序之前，应该调用fftw_cleanup()以防止某些内存泄漏。