[发明专利]利用苹果果梗图像计算苹果旋转角度的方法

摘要

本发明公开了一种利用苹果果梗图像计算苹果旋转角度的方法。包括如下步骤：苹果目标提取、果梗侧图像的特征点提取、两幅果梗图像特征点匹配、特征点旋转角度计算、旋转角度坏点剔除和苹果旋转角度计算。本发明利用苹果果梗作为苹果旋转角度计算的依据，检测速度快，准确度高。

主权项

一种利用苹果果梗图像计算苹果旋转角度的方法，其特征在于，包括如下步骤：1)苹果目标提取对苹果果梗图像I1(x,y)和I2(x,y)进行目标提取，即去除背景图像，提取苹果图像，用最大类间方差法进行分割，提取苹果目标图像，记为Iob_1(x,y)和Iob_2(x,y)2)果梗侧图像的特征点提取对图像Iob_1(x,y)特征点的步骤如下：①建立在不同尺度σ空间内的图像，定义图像Iob_1(x,y)为S(x,y)，按式（1）定义G(x,y,σ)为可变尺度的2维高斯函数，按式（2）定义L(x,y,σ)为不同尺度σ空间的图像，有： G ( x , y , σ ) = 1 2 π σ 2 e x 2 + y 2 2 σ 2 - - - ( 1 ) L ( x , y , σ ) = G ( x , y , σ ) * S ( x , y ) - - - ( 2 ) ; ②对不同尺度σ空间的图像L(x,y,σ)按式（3）进行连续平滑和采样，记D为相邻尺度的高斯差分图像。 D ( x , y , σ ) = L ( x , y , kσ ) - L ( x , y , σ ) - - - ( 3 ) ; ③提取高斯差分图像D上的峰值点作为特征点的候选点。每一个采样点与其同尺度的8个相邻点以及上下相邻尺度对应的9×2个点共26个点比较，将局部极值点集（灰度值极大值或极小值）作为特征点候选点集；④去除极值点集中的不稳定点。利用Taylor二次级数在采样点处的尺度空间函数约束条件剔除对图像噪声敏感的低对比度的极值点；利用Hessian矩阵的轨迹和行列式去除边缘处的候选特征点；⑤确定特征点的方向。记固定尺寸σ的高斯平滑图像为L(x,y,k0σ)，利用特征点的梯度方向信息建立角度直方图，峰值即为局部梯度的主方向，分别按式（4）和式（5）计算像素梯度值m(x,y)和方向θ(x,y)； m ( x , y ) = ( L ( x + 1 , y ) - L ( x - 1 , y ) ) 2 + ( L ( x , y + 1 ) - L ( x , y - 1 ) ) 2 - - - ( 4 ) θ ( x , y ) = tan - 1 ( ( L ( x , y + 1 ) - L ( x , y - 1 ) ) ( L ( x + 1 , y ) - L ( x - 1 , y ) ) ) - - - ( 5 ) ; ⑥建立特征点的描述符：以特征点为中心点采集16×16个像素窗口，将像素窗口分为4×4的子块，则每个子块可得到8个方向的梯度大小信息和梯度方向直方图，用128维向量描述该特征点。记提取的特征点为N1(xn，yn)(n为正整数)；对图像Iob_2(x，y)重复步骤①～⑥，得到其特征点N2(xn，yn)(n为正整数)；3)两幅果梗图像特征点匹配采用k维搜索树(kd树)进行特征点的搜索匹配。以苹果果梗图像I1(x，y)中的特征点N1(xn，yn)(n为正整数)为基准点，搜索与图像I1(x，y)中的特征点最邻近的图像I2(x，y)中的特征点。记图像I1(x，y)和图像I2(x，y)中对应匹配点坐标为(xm，ym)和(xm′，ym′)；4)获取圆的中心坐标采用Canny算子提取苹果图像的边界，对苹果图像的边界点采用最小二乘圆进行拟合，得到图像Iob_1(x，y)和图像Iob_2(x，y)的圆心点坐标，分别为(x0，y0)和(x0′，y0′)；5)特征点旋转角度计算特征点的旋转角度θm计算公式为： θ m = 180 π × ( ar tan y m ′ - y 0 ′ x m ′ - x 0 ′ - ar tan y m - y 0 x m - x 0 ) - - - ( 6 ) ; 6)旋转角度坏点剔除以特征点的旋转角度θm为对象，记为： θ = { θ i | θ i ∈ θ m , i = 1,2 , · · · , m } - - - ( 7 ) 将对象分为3个类Wj(j＝1，2，3)，这3个类的聚类中心分别为：z1，z2，z3，对应的特征点数量为Mj(j＝1，2，3)；聚类过程如下：①从对象中任意选择3个对象作为初始聚类中心；②②用式(8)计算两个对象问的欧式距离： d ( θ i , θ j ) = ( θ i - θ j ) T ( θ i - θ j ) - - - ( 8 ) 将当前对象归类到d(θi，θj)取得最小值时聚类中心θj所在的类；③用式(9)分别计算3个类的聚类中心： Z j = 1 M j Σ x ∈ w j x - - - ( 9 ) ④用式(10)计算J值： J = Σ i = 1 k Σ j = 1 M j d ( θ j , z i ) - - - ( 10 ) ⑤循环步骤②～④，直到J值不再发生变化；7)苹果旋转角度计算取对象数量最多的类的聚类中心为苹果旋转角度。