[发明专利]一种自适应阈值的Harris角点检测方法

权利要求书

1.一种自适应阈值的Harris角点检测方法，其特征在于包含如下步骤：

（1）对原图像I中的每个像素点进行梯度运算，得到M矩阵中的各个元素，M=Ix2IxIyIxIyIy2,]]>其中，I_x为图像在水平方向的梯度图像，I_y为图像在竖直方向的梯度图像，对原图像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值，组成矩阵Matrix_R；

（4）设置迭代终止值K＝1及迭代初值T₀，T₀取步骤(3)矩阵Matrix_R中所有元素的最大值与最小值的算数平均数，即：

T0=Rmax+Rmin2]]>

其中：R_max为原图像中各个像素点计算出的R值中最大的R值，即矩阵Matrix_R中最大的元素值。R_min为原图像中各个像素点计算的R值中最小的R值，即矩阵Matrix_R中最小的元素值。

（5）设置迭代中间变量初值T＝0，并将T₀赋给T；

（6）使用T对Matrix_R中的各个元素进行分类，分为两个区域G₁和G₂，G₁由Matrix_R中所有大于T的元素组成，G₂由Matrix_R中所有小于等于T的元素组成。

（7）分别G₁和G₂两个区域中元素的算数平均值μ₁和μ₂。

（8）计算μ₁和μ₂的算数平均值，并赋给μ，即：

μ=μ1+μ22]]>

（9）设置迭代输出变量初值T_out＝0，求μ与T的差值的绝对值，若差值的绝对值小于等于步骤（4）设定的迭代终止值K，则停止迭代，并将T赋给阈值T_out；若差值绝对值大于K，将μ的值赋给T，然后重复步骤（6）至步骤（9）；

（10）使用阈值T_out提取角点：找出矩阵Matrix_R中元素值大于该元素周边其他8个元素的值且大于阈值T_out的元素，其中矩阵边缘的元素不参与角点提取，找到这些元素所对应的坐标，在原图像中找出坐标所对应的像素点，坐标对应的像素点就是角点，在原图像I中将角点标记出来。