OpenCV

Stella981
• 阅读 628

opencv-自定义harris角点检测

关于harris角点检测的原理以及matlab版本,请移步https://www.cnblogs.com/klitech/p/5779600.html

小白初学,这里采用opencv实现之,把自己遇到的疑问一一表述出来,以备后用。

疑问1. CV_32FC1,CV_32FC(6)是什么意思?

   CV_<bit_depth>(S|U|F)C<number_of_channels>

    • bit_depth,可取的值为8,16,32,64.它表示一个像素所占用的bite数
    • S|U|F, S--- signed int,  U--- unsigned int,  F---float
    • C<number_of_channels> 通道数,灰度图像取值1,rgb图像取值3

疑问2. Mat::zeros的两种初始化

  • Mat::zeros(int rows, int cols, int type)
  • Mat::zeros(Size size,int type)

第一种形式,返回特定尺寸与类型的零矩阵,比如  Mat A = Mat::zeros(3,3,CV_32FC1);

第二种形式,程序中采用的方式,Mat::zeros(src.size(), CV_32FC(6));

我的理解是这两种形式实质一样

疑问3. cornerEigenValsAndVecs()使用方法

   函数原型, cornerEigenValsAndVecs( InputArray src, OutputArray dst,

int blockSize, int ksize,

int borderType = BORDER_DEFAULT );

  • src 图像类型应该为单通道,或者float
  • dst 图像类型应该为CV_32FC(6),包含2个特征值,以及对应的2个2维向量,总计6个结果。
  • blocksize 邻域大小
  • ksize 函数采用sobel算子
  • borderType 取默认BORDER_DEFAULT

函数调用参看后面的程序。

 1 #include <opencv2/opencv.hpp>
 2 #include <iostream>
 3 #include <math.h>
 4 using namespace cv;
 5 using namespace std;
 6 Mat src, gray_src;
 7 Mat harris_dst, harrisRspImg;
 8 double harris_min_rsp;
 9 double harris_max_rsp;
10 int qualityLevel = 30;
11 const char* harris_win = "Custom Harris Corners Dector";
12 int max_count = 100;
13 void CustomHarris_Demo(int, void *);
14 
15 int main()
16 {
17 
18     src = imread("D:/1.png");
19     if (src.empty())
20     {
21         cout << "could not load image..." << endl;
22         return -1;
23     }
24     namedWindow("input_image", CV_WINDOW_AUTOSIZE);
25     imshow("input_image", src);
26     cvtColor(src, gray_src, COLOR_BGR2GRAY);
27     // 计算特征值
28     int blockSize = 3;
29     int ksize = 3;
30     double k = 0.04;
31 
32     harris_dst = Mat::zeros(src.size(), CV_32FC(6)); //6通道
33     harrisRspImg = Mat::zeros(src.size(), CV_32FC1);
34     cornerEigenValsAndVecs(gray_src, harris_dst, blockSize, ksize, 4);
35     //计算响应
36     for (int row = 0; row < harris_dst.rows; row++)
37     {
38         for (int col = 0; col < harris_dst.cols; col++)
39         {
40             double lambda1 = harris_dst.at<Vec6f>(row, col)[0];
41             double lambda2 = harris_dst.at<Vec6f>(row, col)[1];
42             harrisRspImg.at<float>(row, col) = lambda1 * lambda2 - k*pow((lambda1 + lambda2), 2);
43         }
44     }
45     minMaxLoc(harrisRspImg, &harris_min_rsp, &harris_max_rsp, 0, 0, Mat());//求最大最小响应
46     namedWindow(harris_win, CV_WINDOW_AUTOSIZE);
47     createTrackbar("Quality Value", harris_win, &qualityLevel, max_count, CustomHarris_Demo);
48     CustomHarris_Demo(0, 0);
49     waitKey(0);
50     return 0;
51 }
52 void CustomHarris_Demo(int, void*) {
53     if (qualityLevel < 10) {
54         qualityLevel = 10;
55     }
56     Mat resultImg = src.clone();
57     float t = harris_min_rsp + (((double)qualityLevel) / max_count)*(harris_max_rsp - harris_min_rsp);
58     for (int row = 0; row < src.rows; row++) {
59         for (int col = 0; col < src.cols; col++) {
60             float v = harrisRspImg.at<float>(row, col);
61             if (v > t) {
62                 circle(resultImg, Point(col, row), 2, Scalar(0, 0, 255), 2, 8, 0);
63             }
64         }
65     }
66 
67     imshow(harris_win, resultImg);
68 }
点赞
收藏
评论区
推荐文章
blmius blmius
3年前
MySQL:[Err] 1292 - Incorrect datetime value: ‘0000-00-00 00:00:00‘ for column ‘CREATE_TIME‘ at row 1
文章目录问题用navicat导入数据时,报错:原因这是因为当前的MySQL不支持datetime为0的情况。解决修改sql\mode:sql\mode:SQLMode定义了MySQL应支持的SQL语法、数据校验等,这样可以更容易地在不同的环境中使用MySQL。全局s
皕杰报表之UUID
​在我们用皕杰报表工具设计填报报表时,如何在新增行里自动增加id呢?能新增整数排序id吗?目前可以在新增行里自动增加id,但只能用uuid函数增加UUID编码,不能新增整数排序id。uuid函数说明:获取一个UUID,可以在填报表中用来创建数据ID语法:uuid()或uuid(sep)参数说明:sep布尔值,生成的uuid中是否包含分隔符'',缺省为
待兔 待兔
3个月前
手写Java HashMap源码
HashMap的使用教程HashMap的使用教程HashMap的使用教程HashMap的使用教程HashMap的使用教程22
Jacquelyn38 Jacquelyn38
3年前
2020年前端实用代码段,为你的工作保驾护航
有空的时候,自己总结了几个代码段,在开发中也经常使用,谢谢。1、使用解构获取json数据let jsonData  id: 1,status: "OK",data: 'a', 'b';let  id, status, data: number   jsonData;console.log(id, status, number )
Stella981 Stella981
3年前
Opencv中Mat矩阵相乘——点乘、dot、mul运算详解
Opencv中Mat矩阵相乘——点乘、dot、mul运算详解2016年09月02日00:00:36 \牧野(https://www.oschina.net/action/GoToLink?urlhttps%3A%2F%2Fme.csdn.net%2Fdcrmg) 阅读数:59593
Wesley13 Wesley13
3年前
mysql设置时区
mysql设置时区mysql\_query("SETtime\_zone'8:00'")ordie('时区设置失败,请联系管理员!');中国在东8区所以加8方法二:selectcount(user\_id)asdevice,CONVERT\_TZ(FROM\_UNIXTIME(reg\_time),'08:00','0
Wesley13 Wesley13
3年前
00:Java简单了解
浅谈Java之概述Java是SUN(StanfordUniversityNetwork),斯坦福大学网络公司)1995年推出的一门高级编程语言。Java是一种面向Internet的编程语言。随着Java技术在web方面的不断成熟,已经成为Web应用程序的首选开发语言。Java是简单易学,完全面向对象,安全可靠,与平台无关的编程语言。
Stella981 Stella981
3年前
Django中Admin中的一些参数配置
设置在列表中显示的字段,id为django模型默认的主键list_display('id','name','sex','profession','email','qq','phone','status','create_time')设置在列表可编辑字段list_editable
Wesley13 Wesley13
3年前
MySQL部分从库上面因为大量的临时表tmp_table造成慢查询
背景描述Time:20190124T00:08:14.70572408:00User@Host:@Id:Schema:sentrymetaLast_errno:0Killed:0Query_time:0.315758Lock_
Python进阶者 Python进阶者
9个月前
Excel中这日期老是出来00:00:00,怎么用Pandas把这个去除
大家好,我是皮皮。一、前言前几天在Python白银交流群【上海新年人】问了一个Pandas数据筛选的问题。问题如下:这日期老是出来00:00:00,怎么把这个去除。二、实现过程后来【论草莓如何成为冻干莓】给了一个思路和代码如下:pd.toexcel之前把这