VLFeat——SIFT图像特征提取（VC++实现）【转载】

由于VLFeat说明文档中只提供了matlab的调用SIFT函数的实例，没有VC++的调用说明，本文是根据实验，在VS2010环境下，结合opencv开源库计算图像的SIFT的描述子，具体实现细节如下所示：

1.配置VLFeat环境已经在之前的一篇博文VLFeat-----mean sift开源库 中介绍了详细的配置过程，可以参考；

2.包含头文件：

  extern "C"{
#include <vl/generic.h>
#include <vl/stringop.h>
#include <vl/pgm.h>
#include <vl/sift.h>
#include <vl/getopt_long.h>
};

3.实现过程如下：

#include <stdio.h>
#include <tchar.h>
#include <opencv2/opencv.hpp>
#include <stdio.h>
#include <string.h>
#include <ctype.h>
#pragma comment(lib, "opencv_core230d.lib")
#pragma comment(lib, "opencv_imgproc230d.lib")
#pragma comment(lib, "opencv_highgui230d.lib")
#pragma comment(lib, "opencv_objdetect230d.lib")

using namespace cv;
using namespace std;

extern "C"{
#include <vl/generic.h>
#include <vl/stringop.h>
#include <vl/pgm.h>
#include <vl/sift.h>
#include <vl/getopt_long.h>
};

int _tmain(int argc, _TCHAR* argv[])
{
    VL_PRINT ("Hello world!\n") ;
char *ImagePath="a.jpg";
    IplImage *Image=cvLoadImage(ImagePath,0);
//  int min=0;
//  min=Image->width>Image->height?Image->height:Image->width;
    int noctaves=4,nlevels=2,o_min=0;
// noctaves=(int)(log(min)/log(2));
    vl_sift_pix *ImageData=new vl_sift_pix[Image->height*Image->width];
    unsigned char *Pixel;
for (int i=0;i<Image->height;i++)
    {
for (int j=0;j<Image->width;j++)
        {
            Pixel=(unsigned char*)(Image->imageData+i*Image->width+j);
            ImageData[i*Image->width+j]=*(Pixel);
        }
    }
    VlSiftFilt *SiftFilt=NULL;
    SiftFilt=vl_sift_new(Image->width,Image->height,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;
            VlSiftKeypoint *pKeyPoint=SiftFilt->keys;
for (int i=0;i<SiftFilt->nkeys;i++)
            {
                VlSiftKeypoint TemptKeyPoint=*pKeyPoint;
                pKeyPoint++;
                cvDrawCircle(Image,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;
    cvNamedWindow("Source Image",1);
    cvShowImage("Source Image",Image);
    cvWaitKey(0);
    cvReleaseImage(&Image);
    cvDestroyAllWindows();
return 0;
}