FFTW3库的应用。

与世无争的帅哥 提交于 2020-01-16 16:11:15

号称最快的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张图像的复数傅里叶变换了。

 

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!