OpenCV 角点检测(一) Moravec
2017-01-21 22:54
465 查看
什么是角点
角点一般反应的是图像中局部最大值或最小值的孤立点,以角点为中心点的窗口向图像中的任意方向滑动,都会引起较大的灰度变化。
Moravec角点
Moravec角点检测算子于1981年提出,是最早的角点检测算法之一。Moravec角点检测算子通过计算中心像素点的0°,45°,90°,135°四个方向的灰度差的平方和,并取其中最小的值作为中心点的兴趣值。
其中(u,v)为中心点(x,y)的0°,45°,90°,135°方向的像素点坐标。
Moravec角点检测算子相比于其他角点检测算法具有如下特点:
(1).非均匀性响应,很容易受到临边特性干扰的影响,对倾斜边缘的检测响应很强。
(2).不管像素点与中心点的距离是多少,都赋予相同的权重,导致 Moravec角点检测算子对噪声很敏感
opencv实现
原图:
角点检测结果:
角点一般反应的是图像中局部最大值或最小值的孤立点,以角点为中心点的窗口向图像中的任意方向滑动,都会引起较大的灰度变化。
Moravec角点
Moravec角点检测算子于1981年提出,是最早的角点检测算法之一。Moravec角点检测算子通过计算中心像素点的0°,45°,90°,135°四个方向的灰度差的平方和,并取其中最小的值作为中心点的兴趣值。
其中(u,v)为中心点(x,y)的0°,45°,90°,135°方向的像素点坐标。
Moravec角点检测算子相比于其他角点检测算法具有如下特点:
(1).非均匀性响应,很容易受到临边特性干扰的影响,对倾斜边缘的检测响应很强。
(2).不管像素点与中心点的距离是多少,都赋予相同的权重,导致 Moravec角点检测算子对噪声很敏感
opencv实现
// Moravec角点检测 #include <opencv2/opencv.hpp> #include <opencv2/highgui/highgui.hpp> #include <opencv2/imgproc/imgproc.hpp> using namespace std; using namespace cv; cv::Mat MoravecCorners(cv::Mat SrcImage, int kSize, int threshold) { cv::Mat MorImage = SrcImage.clone(); // 获取初始化参数信息 int r = kSize / 2; const int nRows = SrcImage.rows; const int nCols =SrcImage.cols; int nConut = 0; CvPoint *pPoint = new CvPoint[nRows*nCols]; // 图像遍历 for (int i = r; i < SrcImage.rows-r; i++) { for (int j = r; j<SrcImage.cols-r; j++) { int wV1, wV2, wV3, wV4; wV1 = wV2 = wV3 = wV4 = 0; // 计算水平方向窗内兴趣值 for (int k = -r; k < r; k++) wV1 += (SrcImage.at<uchar>(i,j+k)- SrcImage.at<uchar>(i,j+k+1))*(SrcImage.at <uchar>(i,j+k)-SrcImage.at<uchar>(i,j+k+1)); // 计算垂直方向窗内兴趣值 for (int k = -r; k < r; k++) wV2 += (SrcImage.at<uchar>(i+k,j)- SrcImage.at<uchar>(i+k+1,j))*(SrcImage.at <uchar>(i+k,j)-SrcImage.at<uchar>(i+k+1,j)); // 计算45度方向窗内兴趣值 for (int k = -r; k < r; k++) wV3 += (SrcImage.at<uchar>(i+k,j+k)- SrcImage.at<uchar>(i+k+1,j+k+1))*(SrcImage.at <uchar>(i+k,j+k)-SrcImage.at<uchar>(i+k+1,j+k+1)); // 计算135度方向窗内兴趣值 for (int k = -r; k < r; k++) wV4 += (SrcImage.at<uchar>(i+k,j-k)- SrcImage.at<uchar>(i+k+1,j-k-1))*(SrcImage.at <uchar>(i+k,j-k)-SrcImage.at<uchar>(i+k+1,j-k-1)); // 取其中的最小值作为该像素点的最终兴趣值 int value = min(min(wV1,wV2), min(wV3,wV4)); // 若兴趣值大于阈值,则将点的坐标存入数组中 if (value > threshold) { pPoint[nConut] = cvPoint(j,i); nConut++; } } } //绘制兴趣点 for (int i = 0; i < nConut; i++) cv::circle(MorImage, pPoint[i], 5, cv::Scalar(255,0,0)); return MorImage; } int main() { cv::Mat SrcImage = imread("1.jpg",0); if (!SrcImage.data) return -1; //阈值设置为10000 cv::Mat MorImage = MoravecCorners(SrcImage, 5,10000); cv::imshow("srcImage", SrcImage); cv::imshow("MorMat",MorImage); cv::waitKey(0); return 0; }
原图:
角点检测结果:
相关文章推荐
- 架构故障
- mysql高可用之PXC(Percona XtraDB Cluster)
- 权限管理架构
- Linux 文件 操作
- kafka2:常用topic命令
- shell常用操作积累
- CentOS下安装C/C++开发工具包的最佳方式
- mha的搭建步骤(一主一从架构)
- OpenGrok安装
- 基础学习笔记之opencv(24):imwrite函数的使用
- shell脚本中8种字符串截取方法_转自脚本之家
- 解决在Linuxmint18.1环境下fcitx无法显示候选框的问题
- Shell高亮显示
- nginx笔记
- linux内核(kernel)版本号的意义
- 2016的个人总结
- CentOS安装MySQL与登录退出
- php学习笔记:thinkphp 配置在nginx的坑
- OpenCV打开摄像头
- HDU 5938 - Four Operation(贪心)