号称最快的FFT算法。是个C语言程序集。
下面是具体应用方法和注意事项:
以复数傅里叶变换为例了。我用到的就是复数的计算。
1.定义输入输出变量
fftw_complex * fttInput; //输入需要变换的复数数据
fftw_complex * fttOutput;//输出变换后的复数结果
2.定义变换方案
fftw_plan p;//创建变换方案。例如我做的是对图像的傅里叶变换 就选择p = fftw_plan_dft_2d。对一维数据做变换就选择fftw_plan_dft_1d.
3.实施变换
fftw_execute(p1);//实施变换
4.注意事项:一个fftw_plan只能对固定的in/out进行变换, 但可以在变换后改变in的内容(大小不变)以用同一个方案执行新的变换。这个尤其重要。在我的试验中发现第一次对某个图像做傅里叶变换时,性能并没有提升反而慢于opencv,然而对同样大小图像内容不同做多次变换(不重新创建plan)时,性能提升明显。所以提速关键应该在此吧。还有FFTW_FORWARD是正变换,FFTW_BACKWARD是逆变换。
上完整代码如下:
fftw_complex * fttinput;
fftw_complex * fttoutput;
fftw_plan p;
fttinput = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * w* h);
fttoutput = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * w* h);
p = fftw_plan_dft_2d(h, w,fttinput, fttoutput,FFTW_FORWARD, FFTW_ESTIMATE);
for(int i = 0; i < 100; i ++ )
{
//对某组图像做傅里叶变换 不断改变fttinput的内容 用fftw_execute执行傅里叶变换。不断刷新结果
//结果不断刷新在fttoutput 中。
fftw_execute(p);
}
//这样就完成了对100张图像的复数傅里叶变换了。
来源:CSDN
作者:若曦爹
链接:https://blog.csdn.net/gbz3300255/article/details/104003311