opencv之离散傅里叶变换（DFT）

 

opencv之离散傅里叶变换

函数解读：

C++: int getOptimalDFTSize(int vecsize)

源码解读；

[cpp]  view plain  copy

 print ?

1. <span style="font-size:18px;">int cv::getOptimalDFTSize( int size0 )  
2. {  
3.    int a = 0, b = sizeof(optimalDFTSizeTab)/sizeof(optimalDFTSizeTab[0]) -1;  
4.    if( (unsigned)size0 >= (unsigned)optimalDFTSizeTab[b] )  
5.        return -1;  
6.    
7.    while( a < b )//二分查找合适的size  
8.     {  
9.        int c = (a + b) >> 1;  
10.        if( size0 <= optimalDFTSizeTab[c] )  
11.            b = c;  
12.        else  
13.            a = c+1;  
14.     }  
15.    
16.     returnoptimalDFTSizeTab[b];  
17. }</span>  

optimalDFTSizeTab定义在namespace cv中，里边的数值为2^x*3^y*5^z

static const int optimalDFTSizeTab[] = {1,2, 3, 4, 5, 6, 8, 9, 10, 12, 15, 16,…, 2123366400, 2125764000};

 copyMakeBorder

C++: void copyMakeBorder(InputArraysrc, OutputArray dst, int top, int bottom, int left,int right, int borderType, const Scalar& value=Scalar())

src: 源图像

dst: 目标图像，和源图像有相同的类型，dst.cols=src.cols+left+right; dst.rows=src.rows+dst.top+dst.bottom

top:

bottom:

left:

right: 以上四个参数指定了在src图像周围附加的像素个数。

borderType: 边框类型

value: 当borderType==BORDER_CONSTANT时需要指定该值。

示例代码：

[cpp]  view plain  copy

 print ?

1. <span style="font-size:18px;">#include <opencv2/core/core.hpp>  
2. #include<opencv2/highgui/highgui.hpp>  
3. #include<opencv2/imgproc/imgproc.hpp>  
4. #include <iostream>  
5.    
6. using namespace cv;  
7. using namespace std;  
8.    
9. int main(){  
10.        Mat src = imread("fruits.jpg");  
11.        if(!src.data)  
12.        {  
13.               return -1;  
14.        }  
15.    
16.        Mat src_gray;  
17.        cvtColor(src,src_gray,CV_RGB2GRAY);//灰度图像做傅里叶变换  
18.    
19.        int m = getOptimalDFTSize(src_gray.rows);//2,3,5的倍数有更高效率的傅里叶转换  
20.        int n = getOptimalDFTSize(src_gray.cols);  
21.    
22.        Mat dst;  
23.        ///把灰度图像放在左上角，在右边和下边扩展图像，扩展部分填充为0；  
24.        copyMakeBorder(src_gray,dst,0,m-src_gray.rows,0,n-src_gray.cols,BORDER_CONSTANT,Scalar::all(0));  
25.        cout<<dst.size()<<endl;  
26.    
27.        //新建一个两页的array，其中第一页用扩展后的图像初始化，第二页初始化为0  
28.        Mat planes[] = {Mat_<float>(dst), Mat::zeros(dst.size(), CV_32F)};  
29.        Mat  completeI;  
30.        merge(planes,2,completeI);//把两页合成一个2通道的mat  
31.    
32.        //对上边合成的mat进行傅里叶变换，支持原地操作，傅里叶变换结果为复数。通道1存的是实部，通道2存的是虚部。  
33.        dft(completeI,completeI);  
34.    
35.        split(completeI,planes);//把变换后的结果分割到各个数组的两页中，方便后续操作  
36.        magnitude(planes[0],planes[1],planes[0]);//求傅里叶变换各频率的幅值，幅值放在第一页中。  
37.    
38.        Mat magI = planes[0];  
39.        //傅立叶变换的幅度值范围大到不适合在屏幕上显示。高值在屏幕上显示为白点，  
40.        //而低值为黑点，高低值的变化无法有效分辨。为了在屏幕上凸显出高低变化的连续性，我们可以用对数尺度来替换线性尺度:  
41.        magI += 1;  
42.        log(magI,magI);//取对数  
43.        magI= magI(Rect(0,0,src_gray.cols,src_gray.rows));//前边对原始图像进行了扩展，这里把对原始图像傅里叶变换取出，剔除扩展部分。  
44.    
45.    
46.        //这一步的目的仍然是为了显示。 现在我们有了重分布后的幅度图，  
47.        //但是幅度值仍然超过可显示范围[0,1] 。我们使用 normalize() 函数将幅度归一化到可显示范围。  
48.        normalize(magI,magI,0,1,CV_MINMAX);//傅里叶图像进行归一化。  
49.    
50.    
51.        //重新分配象限，使（0,0）移动到图像中心，  
52.        //在《数字图像处理》中，傅里叶变换之前要对源图像乘以（-1）^(x+y)进行中心化。  
53.        //这是是对傅里叶变换结果进行中心化  
54.        int cx = magI.cols/2;  
55.        int cy = magI.rows/2;  
56.    
57.        Mat tmp;  
58.        Mat q0(magI,Rect(0,0,cx,cy));  
59.        Mat q1(magI,Rect(cx,0,cx,cy));  
60.        Mat q2(magI,Rect(0,cy,cx,cy));  
61.        Mat q3(magI,Rect(cx,cy,cx,cy));  
62.    
63.         
64.        q0.copyTo(tmp);  
65.        q3.copyTo(q0);  
66.        tmp.copyTo(q3);  
67.    
68.        q1.copyTo(tmp);  
69.        q2.copyTo(q1);  
70.        tmp.copyTo(q2);  
71.    
72.         
73.    
74.        namedWindow("InputImage");  
75.        imshow("InputImage",src);  
76.    
77.        namedWindow("SpectrumImage");  
78.        imshow("SpectrumImage",magI);  
79.    
80.        waitKey();  
81.        return0;  
82. }</span>  

实验结果：

代码已上传至gihub:opencv之离散傅里叶变换

转载：http://blog.csdn.net/ubunfans/article/details/24787569