win7 regedit打不开:Surf算法学习心得(二)——源码简析

说明：作为初学者，我对于源代码也只是简单的分析，开始和(一)中一样都叫做源码分析，后来感觉自己分析的质量不太好，还是都改为源码简析吧，结合起(一)及后面的心得来看估计效果会好点，呵呵。只是希望对于即将要学习Surf算法的人有一定的帮助就行！对于一些介绍得不对的地方，也希望各位大虾能过指出，相互交流，共同进步！

Surf算法源代码分析

surf算法源代码分为两种文件，学过C/C++的都知道，在此不多说。头文件主要包括：imload.h、ipoint.h、image.h、fasthessian.h、surf.h、surflib.h，其中每个文件用于声明一个特定的相应类，下面大体进行简单介绍。ImLoad.h——声明类ImLoad，主要封装了对图像的读取和保存函数。Image *readImage(const char *filename);                             //从图像文件中读取图像void saveImage(const char *filename, Image* im);           //将图像保存到文件中ipoint.h——声明类Ipoint，主要定义关键点的相应属性。Ipoint(){                                                                                                //构造函数ivec = NULL;ori = 0.0;};~Ipoint(){                                                                                             //析构函数if (ivec)delete [] ivec;};void allocIvec(const int si){                                                         //内存空间分配ivec = new double[si];};double x, y;                                                                                          //特征点在图像中的坐标double scale;                                                                                       //检测范围double strength;                                                                                //特征点的强度double ori;                                                                                           //特征点主方向int laplace;                                                                                           //Laplacian相关的值double *ivec;                                                                                      //特征点描述器(局部特征)image.h——声明类Image，主要定义图像的相关属性。Image(const int w, const int h);                                                 //带参数的构造函数~Image();                                                                                            //析构函数Image(double **pixels, int w, int h);                                      //根据已存在的像素数组构造图像Image(Image *im, bool doubleImSize=false);                   //构造积分图像void setFrame(unsigned char *im);                                       //通过单一的帧到(预初始化)结构void setFrame(Image *im);Image *HalfImage();                                                                    //将图片长宽个变为原来的1/2double getHessian(int *x);                                                        //获得特定点的Hessian检测值int getTrace(int *x);                                                                    //获得Hessian矩阵的迹(线性代数中学过迹)double **getPixels() const;                                                      //获得指向图像像素的指针double getPix(const int x, const int y) const {                 //获得某一指定位置的像素值return _pixels[y][x];}double &getPix(const int x, const int y) {                          //重载getPix并返回引用return _pixels[y][x];}void setPix(const int x, const int y, const double val) {  //设置指定位置的像素值_pixels[y][x] = val;}int getWidth();int getHeight();void setWidth(int wi);void setHeight(int hi);void allocPixels(int w, int h);                                                  //为图像像素分配二维内存空间double *_buf;                                                                                //指向图像实际缓冲区的指针double **_pixels;                                                                         //指向图像像素二维数组的指针int _height, _width;int _orihi;                                                                                        //原始图像高度bool _ref;                                                                                         //对图像是否为引用的一种标志fasthessian.h——声明类FastHessian，主要定义算法中的快速Hessian检测子方法~FastHessian();//带参数的构造方法FastHessian(Image *im, std::vector< Ipoint >& ip, double thres = 0.2, bool doub = false,short int initMasksize = 9, short int samplingStep = 2,short int octaves = 4);void setIimage( Image *iim );                                                  //传入积分图像void getInterestPoints();                                                           //检测图像的所有特征点//在特定位置创建新的点，并在一定范围内void makeIpoint(double x, double y, double scale, double strength=0);void allocateOctave();                                                                //为Octave分配内存空间//Fast non-maximum-suppressionvoid findMaximum(int *borders, int o, int octave);void interpFeature(int s, int row, int col, Image *map,int o, int octave, int movesRemain,int *borders);int fitQuadrat(int s, int r, int c, double &res);Image *_Iimage;Image **_scaleLevel;                                                               //Octavesint _vas[9];                                                                                    //向量变量double _threshold;                                                                    //检测特征点时的阈值bool _doubled;                                                                            //图像是否放大std::vector< Ipoint >& _ipts;                                               //从外部传进来的指向特征点向量的引用short int _initLobe;                                               //在某一方向第二次求导时的初始lobe值，默认为3short int _maxScales;                                                              //Number scalesshort int _maxOctaves;                                                          //Number octavesshort int _sampling;                                                                 //The sampling stepint _width;                                                                                   //积分图像的宽int _height;double _offset[3];                                                                      //二次拟合的结果surf.h——声明surf，主要用于定义Surf中关键点相应的描述器Surf();Surf(Image *im, bool dbl=false, bool usurf=false,bool ext=false, int insi=4);~Surf();int getVectLength();                                                              //获得特征描述器向量的长度void setIpoint(Ipoint *ipt);                                               //为一个需要计算的描述器设置相应点void assignOrientation();                                                   //定向分配重现void makeDescriptor();                                                      //计算不变图像特征void createVector(double scale,                                    //创建向量double row, double col);void createUprightVector(double scale,double row, double col);void AddSample(int r, int c, double rpos,                   //向向量添加样本double cpos, double rx, double cx, int step);void AddUprightSample(int r, int c, double rpos,double cpos, double rx, double cx, int step);void PlaceInIndex(double mag1, int ori1,                  //为向量中样本设置索引值double mag2, int ori2, double rx, double cx);//! Normalise descriptor vector for illumination invariance for Lambertian surfacesvoid normalise();void createLookups();                                                       //创建查找表surflib.h——声明surf算法要用到的库函数。//针对整个图像定义关键点及其相应描述器（关键点附加的详细信息（局部特征））//其中关键点信息保存在向量ipts当中inline void surfDetDes(Image *im, std::vector< Ipoint >& ipts,double thres = 4.0, bool doubleImageSize = false,int initLobe = 3, int samplingStep = 2, int octaves = 4,bool upright = false, bool extended = false, int indexSize = 4)//针对一个给定的特征点，计算相应的描述器（局部特征）inline void surfDes(Image *im, std::vector< Ipoint >& ipts,bool doubleImageSize = false,bool upright = false, bool extended = false, int indexSize = 4)

编译运行

cmd下进入可执行文件目录，输入可执行文件名，得到如下提示：可以在相应的提示下进行运行，如：注意：输入文件 -i 参数后面的文件必须是PGM格式的图像文件，可以自行网上下载，有个“人脸pgm图片库”可以拿来使用，输出不限，如本程序中是output.txt运行完成后，打开文件output.txt可以看到文件中的如下数据：这只是一个简单的示例，直接运行可执行文件名出现的帮助里面还有很多选项，别的选项也就是与Surf算法源代码中的那些参数对应的，大家应该都懂的，要学习的可以试一下，这样就能更深入的了解Surf算法。另外源文件中有文件README的说明 ，里面有关于数据格式以及数据输入输出格式的说明，有兴趣的朋友可以自行研究下