Gabor函数应用

原文：http://blog.csdn.net/yao_zhuang/article/details/2532279

 

下载cvgabor.cpp和cvgabor.h到你的C/C++工程目录下

注：在我的资源中有改进过的cvgabor类

    相关链接为：http://download.csdn.net/source/490114

特别注意：使用该类需要opencv库的支持，如何配置环境参见：http://www.opencv.org.cn/index.php/Template:Install

它有如下的功能：

生成特定方向和尺度的gabor

生成可以显示或者保存的gabor核的实部，虚部

图像的实部，虚部或者主要（Magnitude）响应

响应可以保存在XML文件中

 

[cpp] view plain copy print ?

1. #include "cvgabor.h"   
2. int main(){  
3.   
4.     //创建一个方向是PI/4而尺度是3的gabor   
5.     double Sigma = 2*PI;    
6.     double F = sqrt(2.0);    
7.     CvGabor *gabor1 = new CvGabor;    
8.     gabor1->Init(PI/4, 3, Sigma, F);  
9.       
10.     //获得实部并显示它   
11.     IplImage *kernel = cvCreateImage( cvSize(gabor1->get_mask_width(), gabor1->get_mask_width()), IPL_DEPTH_8U, 1);  
12.     kernel = gabor1->get_image(CV_GABOR_REAL);  
13.     cvNamedWindow("Gabor Kernel", 1);  
14.     cvShowImage("Gabor Kernel", kernel);  
15.     cvWaitKey(0);  
16.   
17.     //载入一个图像并显示   
18.     IplImage *img = cvLoadImage( "D:/Demo.jpg", CV_LOAD_IMAGE_GRAYSCALE );  
19.     cvNamedWindow("Original Image", 1);  
20.     cvShowImage("Original Image", img);  
21.     cvWaitKey(0);  
22.   
23.     //获取载入图像的gabor滤波响应的实部并且显示   
24.     IplImage *reimg = cvCreateImage(cvSize(img->width,img->height), IPL_DEPTH_8U, 1);  
25.     gabor1->conv_img(img, reimg, CV_GABOR_REAL);  
26.     cvNamedWindow("Real Response", 1);  
27.     cvShowImage("Real Response",reimg);  
28.     cvWaitKey(0);  
29.     cvDestroyWindow("Real Response");  
30.   
31.     //获取载入图像的gabor滤波响应的虚部并且显示   
32. //  IplImage *reimg = cvCreateImage(cvSize(img->width,img->height), IPL_DEPTH_8U, 1);   
33.     gabor1->conv_img(img, reimg, CV_GABOR_IMAG);  
34.     cvNamedWindow("Imaginary Response", 1);  
35.     cvShowImage("Imaginary Response",reimg);  
36.     cvWaitKey(0);  
37.     cvDestroyWindow("Imaginary Response");  
38.   
39.     //获取载入图像的gabor滤波响应的模并且显示   
40. //  IplImage *reimg = cvCreateImage(cvSize(img->width,img->height), IPL_DEPTH_8U, 1);   
41.     gabor1->conv_img(img, reimg, CV_GABOR_MAG);  
42.     cvNamedWindow("Magnitude Response", 1);  
43.     cvShowImage("Magnitude Response",reimg);  
44.     cvWaitKey(0);  
45.     cvDestroyWindow("Magnitude Response");  
46.   
47.     /* 
48.     //这个响应可以被取样为8位的灰度图。如果你要原始的浮点类型的数据，你可以这样做 
49.     IplImage *reimg = cvCreateImage(cvSize(img->width,img->height), IPL_DEPTH_32F, 1); 
50.     gabor1->conv_img(img, reimg, CV_GABOR_MAG); 
51.     //然而，这些浮点数据是不能够以上面灰度图的形式简单的显示，但是它可以被保存在一个XML文件中。 
52.     cvSave( "reimg.xml", (IplImage*)reimg, NULL, NULL, cvAttrList(0,0)); 
53.     */  
54.   
55.   
56.   
57. }