角点检测（3）SUSAN算子

#include "stdafx.h"
#include <iostream>
#include <stdlib.h>
#include <cv.h>
#include <cxcore.h>
#include <highgui.h>
#include <math.h>

int main( int argc, char** argv )
{
    int height ,width ,step ,channels ;
         int i,j,k,same ,max,min,thresh,sum;
         uchar*data0,*data1 ;
         //char *filename="result.bmp";
        // IplImage* Img,*nimg; //声明IplImage指针
    //载入图像
       IplImage*  nimg = cvLoadImage("1.jpg");
       IplImage*  Img = cvCreateImage(cvGetSize(nimg),8,1);
		cvCvtColor( nimg, Img, CV_BGR2GRAY );
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=0;i<height;i++)
// for(j=0;j<width;j++)
    //{
   // if(data0[i*step+j]>max) max = data0[i*step+j];
    //if(data0[i*step+j]<min)   min = data0[i*step+j];
//   }
    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 = sum/37;
        if(fabs( (float)(data0[(i+OffSetY[k])*step+(j+OffSetX[k])]-data0[i*step+j]))<=thresh)
        same++;
        if(same<18)
        data1[i*step+j] = 255;
        else
        data1[i*step+j] = 0;
     }
   }
   cvNamedWindow( "Image", 1 ); //创建窗口
        cvShowImage( "Image", nimg ); //显示图像
        cvWaitKey(0); //等待按键
        cvDestroyWindow( "Image" );//销毁窗口
        cvReleaseImage( &Img ); //释放图像
        cvReleaseImage( &nimg );
        return 0;
}