Advertisement

openCV 傅立叶逆变换 inverse fourier transform

阅读量:

最近在写一个people-detection的project, 其中需要用到openCV的傅立叶变换,http://docs.opencv.org/modules/core/doc/operations_on_arrays.html#dft 这里面的dft函数。

调用dft(img, img)便可以得到img的傅立叶变换的值,并替换原img的值。

这时候如果需要逆变换的值,那么你便需要调用idft()或者dft().

这步的结果一直有问题,因为我在调用的时候少加了一个DFT_SCALE!!!

正确的方法:

复制代码
    dft(img, img, DFT_SCALE|DFT_INVERSE);

看我上面给的网页里面官方说明, DFT_SCALE在求逆变换时一般是同时存在的,要不然得到的是放大了N(N为元素总数)倍。

实在搞不明白OpenCV要搞这么一出,浪费了我不知道多少个小时,不过这也告诉我,调用函数时认真读API的文档是很有必要的。

下面是我测试结果是否正确的源码, 在ubuntu的终端编译指令为: g++ -o test test.C pkg-config opencv --cflags --libs:

复制代码
 #include "opencv2/imgproc/imgproc.hpp"

    
 #include "opencv2/objdetect/objdetect.hpp"
    
 #include "opencv2/highgui/highgui.hpp"
    
 #include <iostream>
    
  
    
 using namespace std;
    
 using namespace cv;
    
  
    
 int main(){
    
     Mat img(3, 3, CV_64FC2);
    
     int k = 0;
    
     for(unsigned i = 0; i < 3; i++){
    
     for(unsigned j = 0; j < 3; j++){
    
         img.at<std::complex<double> >(i,j) = k++;
    
     }
    
     }
    
     for(unsigned i = 0; i < 3; i++){
    
     for(unsigned j = 0; j < 3; j++){
    
         cout << img.at<std::complex<double> >(i,j) << endl;
    
     }
    
     dft(img, img);
    
     dft(img, img, DFT_SCALE|DFT_INVERSE);
    
     for(unsigned i = 0; i < 3; i++){
    
     for(unsigned j = 0; j < 3; j++){
    
         cout << img.at<std::complex<double> >(i,j) << endl;
    
     }
    
     }
    
 }
复制代码
复制代码
复制代码

全部评论 (0)

还没有任何评论哟~