Harris角点检测源码

第一种

#include "opencv2/highgui/highgui.hpp"#include "opencv2/imgproc/imgproc.hpp"#include <iostream>#include <stdio.h>#include <stdlib.h>using namespace cv;using namespace std;/// Global variablesMat src, src_gray;int thresh = 200;int max_thresh = 255;char* source_window = "Source image";char* corners_window = "Corners detected";/// Function headervoid cornerHarris_demo( int, void* );/** @function main */int main( int argc, char** argv ){  /// Load source image and convert it to gray  src = imread( argv[1], 1 );  cvtColor( src, src_gray, CV_BGR2GRAY );  /// Create a window and a trackbar  namedWindow( source_window, CV_WINDOW_AUTOSIZE );  createTrackbar( "Threshold: ", source_window, &thresh, max_thresh, cornerHarris_demo );  imshow( source_window, src );  cornerHarris_demo( 0, 0 );  waitKey(0);  return(0);}/** @function cornerHarris_demo */void cornerHarris_demo( int, void* ){  Mat dst, dst_norm, dst_norm_scaled;  dst = Mat::zeros( src.size(), CV_32FC1 );  /// Detector parameters  int blockSize = 2;  int apertureSize = 3;  double k = 0.04;  /// Detecting corners  cornerHarris( src_gray, dst, blockSize, apertureSize, k, BORDER_DEFAULT );  /// Normalizing  normalize( dst, dst_norm, 0, 255, NORM_MINMAX, CV_32FC1, Mat() );  convertScaleAbs( dst_norm, dst_norm_scaled );  /// Drawing a circle around corners  for( int j = 0; j < dst_norm.rows ; j++ )     { for( int i = 0; i < dst_norm.cols; i++ )          {            if( (int) dst_norm.at<float>(j,i) > thresh )              {               circle( dst_norm_scaled, Point( i, j ), 5,  Scalar(0), 2, 8, 0 );              }          }     }  /// Showing the result  namedWindow( corners_window, CV_WINDOW_AUTOSIZE );  imshow( corners_window, dst_norm_scaled );}

第二种

#include "cv.h"#include "highgui.h"#include<iostream>int main(int argc,char** argv){     IplImage* pImg;     IplImage* pHarrisImg;     IplImage* grayImage;     IplImage* dst8;    double minVal=0.0, maxVal=0.0;     double scale, shift;     double min=0, max=255;     if((pImg=cvLoadImage("/home/dz/study/7class/a.jpg",1))!=NULL)     {         cvNamedWindow("source",1);         cvShowImage("source",pImg);         pHarrisImg=cvCreateImage(cvGetSize(pImg),IPL_DEPTH_32F,1);         //there we should define IPL_DEPTH_32F rather than IPL_DEPTH_8U         grayImage=cvCreateImage(cvGetSize(pImg),IPL_DEPTH_8U,1);         dst8=cvCreateImage(cvGetSize(pImg),IPL_DEPTH_8U,1);//this is for the result image         grayImage->origin=pImg->origin;  //there make sure the same  origin between grayImage and pImg         cvCvtColor(pImg,grayImage,CV_BGR2GRAY);//cause harris need gray scale image,we should convert RGB 2 gray         int block_size=7;         //do harris algorithm         cvCornerHarris(grayImage,pHarrisImg,block_size,3,0.04);         //convert scale so that we see the clear image         cvMinMaxLoc(pHarrisImg,&minVal,&maxVal,NULL,NULL,0);          std::cout<<minVal<<std::endl;          std::cout<<maxVal<<std::endl;         scale=(max-min)/(maxVal-minVal);         shift=-minVal*scale+min;         cvConvertScale(pHarrisImg,dst8,scale,shift);         cvNamedWindow("Harris",1);         cvShowImage("Harris",dst8);         cvWaitKey(0);         cvDestroyWindow("source");         cvDestroyWindow("Harris");         cvReleaseImage(&dst8);         cvReleaseImage(&pHarrisImg);         return 0;     }     return 1;}