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)
还没有任何评论哟~
