点特征提取算法实现.docx

Harris_ex实现代码：#include<opencv2/opencv.hpp>#include<iostream>using namespace std;using namespace cv;#define WINDOW_NAME1 "【程序窗口1】"#define WINDOW_NAME2 "【程序窗口2】"Mat g_srcImage, g_srcImage1, g_grayImage;int thresh = 30; //当前阈值int max_thresh = 175; //最大阈值void on_CornerHarris(int, void*); //回调函数int main() {g_srcImage = imread("Harris_example.jpg");imshow("【原始图】", g_srcImage);g_srcImage1 = g_srcImage.clone();cvtColor(g_srcImage1, g_grayImage, COLOR_BGR2GRAY);namedWindow(WINDOW_NAME1);createTrackbar("阈值：", WINDOW_NAME1, &thresh, max_thresh, on_CornerHarris);on_CornerHarris(0, 0);while (waitKey(9) != 27);return 0;}void on_CornerHarris(int, void*) {Mat dstImage, normImage, scaledImage;dstImage = Mat::zeros(g_srcImage.size(), CV_32FC1);g_srcImage1 = g_srcImage.clone();cornerHarris(g_grayImage, dstImage, 2, 3, 0.04); //进行角点检测normalize(dstImage, normImage, 0, 255, NORM_MINMAX); //归一化convertScaleAbs(normImage, scaledImage); //转换为8位无符号整形//将检测到的，且符合阈值条件的角点绘制出来for (int j = 0; j < normImage.rows; j++) {for (int i = 0; i < normImage.cols; i++) {if ((int)normImage.at<float>(j, i)>thresh + 80) {circle(g_srcImage1, Point(i, j), 5, Scalar(10, 10, 255), 2);circle(scaledImage, Point(i, j), 5, Scalar(10, 10, 255), 2);}}}imshow(WINDOW_NAME1, g_srcImage1);imshow(WINDOW_NAME2, scaledImage);}程序结果：shi-Tomasi实现代码：#include<opencv2/opencv.hpp>#include<iostream>using namespace std;using namespace cv;#define WINDOW_NAME "【Shi-Tomasi 角点检测】"Mat g_srcImage, g_grayImage;int g_maxCornerNumber = 33, g_maxTrackbarNumber = 500;RNG g_rng;void on_GoodFeaturesToTrack(int, void*);int main() {g_srcImage = imread("Shi-Tomasi.jpg");cvtColor(g_srcImage, g_grayImage, COLOR_BGR2GRAY);namedWindow(WINDOW_NAME);createTrackbar("最大角点数", WINDOW_NAME, &g_maxCornerNumber, g_maxTrackbarNumber, on_GoodFeaturesToTrack);imshow(WINDOW_NAME, g_srcImage);on_GoodFeaturesToTrack(0, 0);while (waitKey(9) != 27);return 0;}void on_GoodFeaturesToTrack(int, void*) {if (g_maxCornerNumber <= 1) g_maxCornerNumber = 1;vector<Point2f> corner