一种自适应阈值的Harris角点检测方法技术

本发明专利技术提供了一种自适应阈值的Harris角点检测方法，计算每个像素点的R值，通过使用自适应迭代阈值算法的策略，求取针对图像来讲最合适的阈值，使用阈值及临近角点剔除策略提取角点。本发明专利技术避免人工对阈值的设置，能使角点分布的更加均匀，提高角点检测的精确性，避免了角点聚簇以及伪角点的产生，在不影响Harris角点检测算法计算方便与稳定的前提下，避免的原算法因为阈值选取而带来的不便，同时有效的避免了聚簇现象和伪角点的出现。本发明专利技术在阈值自适应选取方面达到了较好的效果。

【技术实现步骤摘要】

【技术保护点】
一种自适应阈值的Harris角点检测方法，其特征在于包含如下步骤：（1）对原图像I中的每个像素点进行梯度运算，得到M矩阵中的各个元素，M=Ix2IxIyIxIyIy2,其中，Ix为图像在水平方向的梯度图像，Iy为图像在竖直方向的梯度图像，对原图像I求卷积求得，：Ix=∂I∂x=I⊗(-1,0,1)Iy=∂I∂y=I⊗(-1,0,1)T其中，I表示原图像；表示卷积；即由水平方向梯度图像的各个像素点的灰度值的平方值组成的图像；即计算垂直方向梯度图像的各个像素点的灰度值的平方值组成的图像；（2）对M矩阵中的各个元素分别进行高斯平滑滤波，得到新的元素图像I“；（3）利用公式R={Ix2×Iy2-(IxIy)2}-0.04{Ix2+Iy2}计算步骤(2)中图像I“中各个像素点的R值，组成矩阵MatrixR；（4）设置迭代终止值K＝1及迭代初值T0，T0取步骤(3)矩阵MatrixR中所有元素的最大值与最小值的算数平均数，即：T0=Rmax+Rmin2其中：Rmax为原图像中各个像素点计算出的R值中最大的R值，即矩阵MatrixR中最大的元素值。Rmin为原图像中各个像素点计算的R值中最小的R值，即矩阵MatrixR中最小的元素值。（5）设置迭代中间变量初值T＝0，并将T0赋给T；（6）使用T对MatrixR中的各个元素进行分类，分为两个区域G1和G2，G1由MatrixR中所有大于T的元素组成，G2由MatrixR中所有小于等于T的元素组成。（7）分别G1和G2两个区域中元素的算数平均值μ1和μ2。（8）计算μ1和μ2的算数平均值，并赋给μ，即：μ=μ1+μ22（9）设置迭代输出变量初值Tout＝0，求μ与T的差值的绝对值，若差值的绝对值小于等于步骤（4）设定的迭代终止值K，则停止迭代，并将T赋给阈值Tout；若差值绝对值大于K，将μ的值赋给T，然后重复步骤（6）至步骤（9）；（10）使用阈值Tout提取角点：找出矩阵MatrixR中元素值大于该元素周边其他8个元素的值且大于阈值Tout的元素，其中矩阵边缘的元素不参与角点提取，找到这些元素所对应的坐标，在原图像中找出坐标所对应的像素点，坐标对应的像素点就是角点，在原图像I中将角点标记出来。FDA00003609800500013.jpg,FDA00003609800500014.jpg,FDA00003609800500015.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：谢红梅，王志成，李会方，冯晓毅，吴俊，彭进业，蒋晓悦，何贵青，王保平，王珺，李永恒，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：