用VLFeat库进行SIFT特征提取(C++ 实现)

说明:

                                                         ~~~~代码来自网络,无法追溯源头,其中部分稍有改动,以适应系统环境~~~
                                                                   环境:win 7 64专业版,VS2013,opencv2.4.9,Vlfeat 0.9.20.
                                                     配置:可参考我之前的博文http://blog.csdn.net/lilai619/article/details/50538303

具体改动:用Mat类型替代原始代码中的IPlImage

#include <iostream>
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp> 
#include <stdio.h>//标准输入输出头文件,编写C语言中,常用到printf()和scanf()函数.
#include <string.h>//是C版本的头文件,包含比如strcpy、strcat之类的字符串处理函数。

#define TRUE 1//define 用于简单的字符替换

extern "C"{
#include <vl/generic.h>
#include <vl/sift.h>
}

using namespace cv;
using namespace std;

int main(int argc, char* argv[])
{

	cout << "SIFT特征提取演示程序" << endl;
	VL_PRINT("Hello world!\n");

	Mat img_1 = imread("1.jpg", 0);
	int noctaves = 4, nlevels = 2, o_min = 0;

	// noctaves=(int)(log(min)/log(2));  
	//在sift.h头文件中有 typedef float vl_sift_pix 声明

	vl_sift_pix *ImageData = new vl_sift_pix[img_1.cols*img_1.rows];
	for (int i = 0; i < img_1.rows; i++)
	{
		uchar *Pixel = img_1.ptr<uchar>(i);
		for (int j = 0; j < img_1.cols; j++)
		{
			ImageData[i*img_1.cols + j] = Pixel[j];
		}
	}
	// 定义VlSiftFilt结构体指针
	VlSiftFilt *SiftFilt = NULL;
	// 创建一个新的sift滤波器
	SiftFilt = vl_sift_new(img_1.cols, img_1.rows, noctaves, nlevels, o_min);

	int KeyPoint = 0;
	int idx = 0;
	if (vl_sift_process_first_octave(SiftFilt, ImageData) != VL_ERR_EOF)
	{
		while (TRUE)
		{
			//计算每组中的关键点  
			vl_sift_detect(SiftFilt);
			//遍历并绘制每个点  
			KeyPoint += SiftFilt->nkeys;//检测到的关键点的数目
			cout << "关键点数目:" << endl << KeyPoint << endl;
			VlSiftKeypoint *pKeyPoint = SiftFilt->keys;//检测到的关键点
			for (int i = 0; i<SiftFilt->nkeys; i++)
			{
				VlSiftKeypoint TemptKeyPoint = *pKeyPoint;
				pKeyPoint++;
				circle(img_1, cvPoint(TemptKeyPoint.x, TemptKeyPoint.y), TemptKeyPoint.sigma / 2, CV_RGB(255, 0, 0));

				idx++;
				//计算并遍历每个点的方向  
				double angles[4];
				int angleCount = vl_sift_calc_keypoint_orientations(SiftFilt, angles, &TemptKeyPoint);//计算关键点的方向
				for (int j = 0; j<angleCount; j++)
				{
					double TemptAngle = angles[j];
					printf("%d: %f\n", j, TemptAngle);
					//计算每个方向的描述  
					float *Descriptors = new float[128];
					vl_sift_calc_keypoint_descriptor(SiftFilt, Descriptors, &TemptKeyPoint, TemptAngle);
					int k = 0;
					while (k<128)
					{
						printf("%d: %f", k, Descriptors[k]);
						printf("; ");
						k++;
					}

					printf("\n");
					delete[]Descriptors;
					Descriptors = NULL;
				}
			}
			//下一阶  
			if (vl_sift_process_next_octave(SiftFilt) == VL_ERR_EOF)
			{
				break;
			}
			//free(pKeyPoint);  
			KeyPoint = NULL;
		}
	}
	vl_sift_delete(SiftFilt);
	delete[]ImageData;
	ImageData = NULL;
	namedWindow("Source Image", 1);
	imshow("Source Image", img_1);

	waitKey(0);

	return 0;
}

运行效果:




运行报错:提示缺少vl.dll
   将安装vlfeat\bin\win32路径下的vl.dll拷贝到项目生成的debug文件夹再编译就OK了.

评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

微风❤水墨

你的鼓励是我最大的动力!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值