角点检测 susan角点检测

程序参见：

code：

#include "stdafx.h"#include "cv.h"#include "highgui.h"#define max_corners 100int main( int argc, char** argv ){int cornerCount=max_corners;CvPoint2D32f corners[max_corners];double qualityLevel = 0.05;double minDistance = 5;IplImage *srcImage = 0, *grayImage = 0, *corners1 = 0, *corners2 = 0;int i;CvScalar color = CV_RGB(255,0,0);cvNamedWindow( "image", 2); //创建显示窗口//加载一副图片srcImage = cvLoadImage("33.jpg", 1);grayImage = cvCreateImage(cvGetSize(srcImage), IPL_DEPTH_8U, 1);//将原图灰度化cvCvtColor(srcImage, grayImage, CV_BGR2GRAY);corners1= cvCreateImage(cvGetSize(srcImage), IPL_DEPTH_32F, 1);corners2= cvCreateImage(cvGetSize(srcImage),IPL_DEPTH_32F, 1);cvGoodFeaturesToTrack (grayImage, corners1, corners2, corners,&cornerCount, qualityLevel, minDistance, 0);printf("num corners found: %d\n", cornerCount);if(cornerCount>0){for (i=0; i <cornerCount;++i){cvCircle(srcImage, cvPoint((int)(corners[i].x), (int)(corners[i].y)), 6,color, 2, CV_AA, 0);}} cvShowImage("image", srcImage);cvWaitKey(0);cvReleaseImage(&srcImage);cvReleaseImage(&grayImage);cvReleaseImage(&corners1);cvReleaseImage(&corners2);return 0; }

另外一个：

忘记在哪里下的了。。感谢原作者：

　　　　// eyes_susan.cpp : 定义控制台应用程序的入口点。//可以检测眼睛动态匹配轮廓的susan角点检测算法，准确率98#include "stdafx.h"#include<math.h>#include<stdlib.h>#include<stdio.h>#include "cv.h"#include "highgui.h"int main( int argc, char** argv ){    int height ,width ,step ,channels ;         int i,j,k,same ,max,min,sum; float thresh;         uchar *data0,*data1 ;         //char *filename="result.bmp";         IplImage* Img, *nimg; //声明IplImage指针    //载入图像        Img = cvLoadImage( "yy.jpg",0);        cvNamedWindow( "Image0", 2); //创建窗口        cvShowImage( "Image0", Img ); //显示图像        nimg = cvCreateImage(cvGetSize(Img),8,1);        height    = Img->height;width     = Img->width;step = Img->widthStep/sizeof(uchar);channels = Img->nChannels;        data0   = (uchar*)Img->imageData;        data1 =    (uchar*)nimg->imageData;printf("Processing a %d X %d image with %d channels\n",width,height,channels);int OffSetX[37] = { -1, 0, 1,-2,-1, 0, 1, 2,-3,-2,-1, 0, 1, 2, 3,-3,-2,-1, 0, 1, 2, 3,-3,-2,-1, 0, 1, 2, 3,-2,-1, 0, 1, 2,-1, 0, 1 };int OffSetY[37] = { -3,-3,-3,-2,-2,-2,-2,-2,-1,-1,-1,-1,-1,-1,-1,0, 0, 0, 0, 0, 0, 0,1, 1, 1, 1, 1, 1, 1,2, 2, 2, 2, 2,3, 3, 3 };    max = min = data0[0];    for(i=3;i<height-3;i++)    for(j=3;j<width-3;j++)   {        same =0;sum = 0;        for(k=0;k<37;k++){        sum+=data0[(i+OffSetY[k])*step+(j+OffSetX[k])];        thresh = (float)sum/37;float data_fabs;data_fabs= (float)(data0[(i+OffSetY[k])*step+(j+OffSetX[k])]-data0[i*step+j]);if(fabs( data_fabs)<=thresh)        same++;}        if(same<18)        nimg->imageData[i*step+j] = 255;        else        nimg->imageData[i*step+j]  = 0;printf("same = %d\n", same);    }        cvNamedWindow( "Image", 2); //创建窗口        cvShowImage( "Image", nimg ); //显示图像        cvWaitKey(0); //等待按键        cvDestroyWindow( "Image" );//销毁窗口        cvReleaseImage( &Img ); //释放图像        cvReleaseImage( &nimg );        return 0;}