harris角点检测

1、课程实验报告2017 - 2018学年第 一学期课程名称： 计算机视觉及应用 实验名称： harris角点检测 班 级： 学生姓名: 学号: 实验日期: 2017.12.6 地点: 指导教师: 成绩评定: 批改日期: 实验目的及要求本实验旨在清楚计算机视觉中角点的概念，学习推导并使用一种经典的角点检测算法，利用合适的编程仿真软件，用程序实现该角点检测算法并用此程序成功提取出示例图片上的角点。实验仪器设备实验设备为一台装有win10系统的PC，qt+opencv软件。实验原理 角点的概念：在现实世界中，角点对应于物体的拐角，道路的十字路口、丁字路口等。从图像分析的角度来说，角点可以是两个边缘的角

2、点，角点也是领域内具有两个主方向的特征点，角点附近区域的像素点无论在梯度方向上还是其梯度幅值上有着较大变化。Harris角点检测算法基本思想：使用一个固定窗口在图像上进行任意方向上的滑动，比较窗口中的像素灰度在窗口滑动前与滑动后的变化程度。如果在任意方向上的滑动都有着较大灰度值的变化，那么我们可以认为该窗口中存在角点。实验内容1.理解Harris角点检测算法的实现原理；2.从理论上推导Harris角点检测算法；3.将Harris角点算法转变为代码的形式，通过程序验证Harris角点算法的可行性；4.自定一张示例图片，用自己编写的程序提取出该示例图片中的角点。实验步骤及方法1.Harris角点检

3、测算法的理论推导：首先建立一个数学模型：当窗口发生移动时，滑动前与滑动后的窗口中的像素点灰度变化用以下方式描述：Eu，v=u,v(ix+u,y+v-ix,y)2其中ix+u,y+v，ix,y分别表示平移前后灰度值的变化。在上式中添加权重wx,y，则描述其变化函数为：Eu，v=u,vwx,y(ix+u,y+v-ix,y)2对ix+u,y+v进行一阶泰勒展开：ix+u,y+v=ix,y+ixx,yu+iyx,yv+o(x2+y2)因为o(x2+y2)对变化贡献不大，故忽略此项不计，则上式可写为：Eu，v=u,vwx,yixx,yu+iyx,yv2将上式平方项展开并写成矩阵形式，即：ixx,yu+i

4、yx,yv2=ixx,y2u2+iyx,y2v2+2ixx,yiyx,yuv =u，vixx,yixx,yixx,yiyx,yiyx,yixx,yiyx,yiyx,y(uv)对于变化函数Eu，v，其本质为自相关函数，对其进行变换为：Eu，v=u，vu,vwix2u,vu,vwiyx,yixx,yu,vwiyx,yixx,yu,vwiy2u,v(uv)令A=u,vwix2u,v，B=u,vwiyx,yixx,y，C=u,vwiy2u,v，则上式可写为：Eu，v=u，vABBCuv=MM是一个实对称阵，则其必定可以相似对角化且存在：PMP-1=PMPT，存在一个变换矩阵P，使得x=Px，ABBC=

5、P1002P-1= P1002PT。求出其特征值及特征向量分别为：1,2=（A+C）±(A-C)2+(2B)22V1,2=12A-C±A-C2+2B22,1TV1和V2正交，进一步化简E：E=u，vx1y1x2y21002x1x2y1y2uv=(u',v') 1002(u'v')= 1u'2+2v'2其中u',v'=u，vx1y1x2y2，由此看出E表征一个椭圆，而x1y1x2y2矩阵表示对u，v旋转。画出Eu，v=1的等高线L，即：E=u'2(1/1)2+v'2(1/2)2=1现在我们可以通过

6、1和2去衡量变化量。此时有三种情况：假设12不考虑v的前提下，u只需变化很小的量便可以使得E=1，同理21，不考虑u的前提下，v只需变化很小的量便可以使得E=1，此时便可以说明在某一个方向上的变化很大，也就是说此时检测到的是图像的边缘。假设1，2都很很小，那么u，v需很大的变化量才可以使得E=1，此时便可以说明在任意方向上的变化不明显，也就是说此时检测到的是图像平坦区域（灰度值变化不大的区域）。假设1，2都很大，那么u，v只需很小的变化量就可以使得E=1，此时便可以说明在任意方向上的变化剧烈，也就是说此时检测到的是图像角点。为了衡量这三种变化，找出其特征值与其的对应关系，建立数学表达式：R=1

7、2-k(1+2)2当R<0，检测到的为图像边缘；当|R|非常小的时候，检测到的为图像中的平坦区域；当R>0，检测到的为图像角点。2.使用qt+opencv编写Harris角点检测程序，具体代码见实验数据；3.使用编写的代码检测一张示例图片的角点。实验数据OPENCV程序代码：#include <opencv2/opencv.hpp>#include "opencv2/highgui/highgui.hpp"#include "opencv2/imgproc/imgproc.hppusing namespace cv;using namesp

8、ace std;Mat g_srcImage, g_srcImage1,g_grayImage;int thresh = 30; /当前阈值int max_thresh = 175; /最大阈值void on_CornerHarris( int, void* );/回调函数int main( int argc, char* argv ) g_srcImage = imread( "H:/picture/a.jpg"); if(!g_srcImage.data ) printf("no picture n"); return false; imshow(&

9、quot;原始图",g_srcImage); g_srcImage1=g_srcImage.clone( ); cvtColor( g_srcImage1, g_grayImage, CV_BGR2GRAY ); on_CornerHarris( 0, 0 ); waitKey(0); return(0);void on_CornerHarris( int, void* ) Mat dstImage;/目标图 Mat normImage;/归一化后的图 Mat scaledImage;/线性变换后的八位无符号整型的图 dstImage = Mat:zeros( g_srcImage.

10、size(), CV_32FC1 ); g_srcImage1=g_srcImage.clone( ); cornerHarris( g_grayImage, dstImage, 2, 3, 0.04, BORDER_DEFAULT ); normalize( dstImage, normImage, 0, 255, NORM_MINMAX, CV_32FC1, Mat() ); convertScaleAbs( normImage, scaledImage ); for( int j = 0; j < normImage.rows ; j+ ) for( int i = 0; i &l

11、t; normImage.cols; i+ ) if( (int) normImage.at<float>(j,i) > thresh+80 ) circle( g_srcImage1, Point( i, j ), 5, Scalar(10,10,255), 2, 8, 0 ); circle( scaledImage, Point( i, j ), 5, Scalar(0,10,255), 2, 8, 0 ); imshow("g_srcImage1" , g_srcImage1 ); imshow( "scaledImage", scaledImage );实验数据分析及处理示例图片角点检测情况： 图一 图二图一为示例图像，是最大阈值为160的时候harris的角点检测，其中红色小点即为检测到的角点，图二中是最大阈值为130的时候harris的角点检测，其中红色小点即为检测到的角点实验结果分析通过用示例