DoG角点检测（C++&openCV2.4.6版本）

#include <opencv2/opencv.hpp> #include <iostream> using namespace std; using namespace cv; int main() { Mat ori_img=imread("G:\\nianfudong\\zhixin.jpg"); Mat gray_img; cvtColor(ori_img,gray_img,CV_RGB2GRAY); gray_img.convertTo(gray_img,CV_32F); Mat gauss1, gauss2; GaussianBlur(gray_img,gauss1,Size(5,5),0.3,0.3); GaussianBlur(gray_img,gauss2,Size(5,5),0.4,0.4); Mat DoG1,DoG2,DoG3; DoG1=gauss1-gauss2; GaussianBlur(gray_img,gauss1,Size(5,5),0.6,0.6); GaussianBlur(gray_img,gauss2,Size(5,5),0.7,0.7); DoG2=gauss1-gauss2; GaussianBlur(gray_img,gauss1,Size(5,5),0.7,0.7); GaussianBlur(gray_img,gauss2,Size(5,5),0.8,0.8); DoG3=gauss1-gauss2; for (int j=1;j<gray_img.rows-1;j++) { for (int i=1;i<gray_img.cols-1;i++) { if (DoG2.at<float>(j,i)<DoG2.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)<DoG2.at<float>(j-1,i)&& DoG2.at<float>(j,i)<DoG2.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)<DoG2.at<float>(j,i-1)&&DoG2.at<float>(j,i)<DoG2.at<float>(j,i+1)&& DoG2.at<float>(j,i)<DoG2.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)<DoG2.at<float>(j+1,i)&&DoG2.at<float>(j,i)<DoG2.at<float>(j+1,i+1) &&DoG2.at<float>(j,i)<DoG1.at<float>(j,i)&&DoG2.at<float>(j,i)<DoG1.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)<DoG1.at<float>(j-1,i)&& DoG2.at<float>(j,i)<DoG1.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)<DoG1.at<float>(j,i-1)&&DoG2.at<float>(j,i)<DoG1.at<float>(j,i+1)&& DoG2.at<float>(j,i)<DoG1.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)<DoG1.at<float>(j+1,i)&&DoG2.at<float>(j,i)<DoG1.at<float>(j+1,i+1) &&DoG2.at<float>(j,i)<DoG3.at<float>(j,i)&&DoG2.at<float>(j,i)<DoG3.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)<DoG3.at<float>(j-1,i)&& DoG2.at<float>(j,i)<DoG3.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)<DoG3.at<float>(j,i-1)&&DoG2.at<float>(j,i)<DoG3.at<float>(j,i+1)&& DoG2.at<float>(j,i)<DoG3.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)<DoG3.at<float>(j+1,i)&&DoG2.at<float>(j,i)<DoG3.at<float>(j+1,i+1)) { if (DoG2.at<float>(j,i)<-3) { circle(ori_img,Point(i,j),3,CV_RGB(0,0,255)); } } else if (DoG2.at<float>(j,i)>DoG2.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)>DoG2.at<float>(j-1,i)&& DoG2.at<float>(j,i)>DoG2.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)>DoG2.at<float>(j,i-1)&&DoG2.at<float>(j,i)>DoG2.at<float>(j,i+1)&& DoG2.at<float>(j,i)>DoG2.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)>DoG2.at<float>(j+1,i)&&DoG2.at<float>(j,i)>DoG2.at<float>(j+1,i+1) &&DoG2.at<float>(j,i)>DoG1.at<float>(j,i)&&DoG2.at<float>(j,i)>DoG1.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)>DoG1.at<float>(j-1,i)&& DoG2.at<float>(j,i)>DoG1.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)>DoG1.at<float>(j,i-1)&&DoG2.at<float>(j,i)>DoG1.at<float>(j,i+1)&& DoG2.at<float>(j,i)>DoG1.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)>DoG1.at<float>(j+1,i)&&DoG2.at<float>(j,i)>DoG1.at<float>(j+1,i+1) &&DoG2.at<float>(j,i)>DoG3.at<float>(j,i)&&DoG2.at<float>(j,i)>DoG3.at<float>(j-1,i-1)&&DoG2.at<float>(j,i)>DoG3.at<float>(j-1,i)&& DoG2.at<float>(j,i)>DoG3.at<float>(j-1,i+1)&&DoG2.at<float>(j,i)>DoG3.at<float>(j,i-1)&&DoG2.at<float>(j,i)>DoG3.at<float>(j,i+1)&& DoG2.at<float>(j,i)>DoG3.at<float>(j+1,i-1)&&DoG2.at<float>(j,i)>DoG3.at<float>(j+1,i)&&DoG2.at<float>(j,i)>DoG3.at<float>(j+1,i+1)) { if (DoG2.at<float>(j,i)>3) { circle(ori_img,Point(i,j),3,CV_RGB(255,0,0)); } } } } imshow("corner",ori_img); waitKey(0); }