[发明专利]一种工件表面轮廓的二维激光实时探测方法

权利要求书

1.一种工件表面轮廓的二维激光实时探测方法，其特征在于，包括以下步骤：

1)从2D激光传感器中获取到一帧轮廓点数据，每个点用如下符号表示：

R＝{r₁，...，r_l}；

2)对获取的轮廓点数据做预处理，预处理包括零值轮廓点去除和重复轮廓点去除；

3)对预处理过后的轮廓点数据进行线段提取，得到轮廓点数据的线段拟合结果；

4)对线段提取的结果进行类内点删除操作，去除错误分类的边缘点，保证整个点集的直线拟合参数的正确性；

5)轮廓形状检测，获得轮廓点数据线段提取的结果之后，把所要提取的目标形状检测并定位出来；

所述对预处理过后的轮廓点数据进行线段提取的方法，包括以下步骤：

(1)使用自适应阈值的迭代中断点拟合IEPF算法对轮廓点数据进行过分割，通过点集R的起点r₁和终点r_l确定直线L，在该点集中搜索到该直线L距离最大的点r_j，验证准则为若该点到该直线L的距离d_j大于设定的阈值threshold，则以该点为断点，将该点集分割为R'和R″两个点集，之后对每个分割点集执行相同的分割过程，直到所有点集都满足验证准则；

得到过分割之后的点集集合V：

V＝{q₁，...，q_n}，q_k＝{r_i，...，r_j}，k∈n

其中q_k为第k个过分割点集，r为该点集所包含的点；

(2)建立无向图G＝[v，ε]，节点v为过分割之后的n个点集V＝{q₁，...，q_n}，无向图的边ε则表征两个节点之间的关系，分为主边和次边，定义ε表征两个点集能进一步融合的概率，其计算函数如下：

$F(q_i,q_j)=\frac{f_r(q_i,q_j)+\mathrm{\α}*f_k(q_i,q_j)+\mathrm{\β}*f_s(q_i,q_j)+\mathrm{\γ}*f_d(q_i,q_j)}{1+\mathrm{\α}+\mathrm{\β}+\mathrm{\γ}}$

其中q_i，q_j为任意两个节点,α，β，γ为常量；

f_r(q_i，q_j)为粗IEPF融合项：

其中length(q_i)为q_i点集所含轮廓点的长度，D_mid(q_i)为q_i集合远端点到IEPF拟合直线的距离，D_MAX为IEPF融合阈值,D_mid为D_mid(q_i)和D_mid(q_j)中较小的那个值；

f_k(q_i，q_j)为kmeans聚类项：

$f_k(q_i,q_j)=\left\{\begin{array}{cc}1& q_i\∈S_{i}and{q}_j\∈S_i\\ \frac{K\left(q_i\right)-K_{mean}\left(S_i\right)}{K\left(q_j\right)-K_{mean}\left(S_i\right)}& q_i\∈S_{i}and{q}_j\∉S_i\end{array}\right.$

其中S_i为kmeans聚类集合S中的第i类，K(q_i)为q_i集合直线拟合斜率，K_mean(S_i)为kmeans聚类后集合的斜率均值；

f_s(q_i，q_j)为锯齿边缘项:

$f_s(q_i,q_j)=\left\{\begin{array}{cc}\sqrt{\frac{{\left({\mathrm{length}}_{\min }\right(q_1,\mathrm{...},q_n\left)\right)}^2}{length\left(q_i\right)*length\left(q_j\right)}}& \begin{array}{c}j=i+1and{k}_{q_i}*k_{q_j}\<0and\\ k_{q_i}*k_{q_{i-1}}\<0and{k}_{q_j}*k_{q_{j+1}}\<0\end{array}\\ 0& else\end{array}\right.$

其中q_i和q_j为点集p的左右临近点集，length_min(q₁，...，q_n)为所有点集中的最短点集的长度，k为该点集的平均斜率；

f_d(q_i，q_j)为点集距离项:

$f_d(q_i,q_j)=\exp (-\frac{1}{2\mathrm{\σ}}||c_{qi}-c_{qj}|\left|\right)$

其中为点集q_i，q_j中距离最近的两个点的坐标；

(3)创建分割无向图的能量函数，使分割无向图的能量函数最小化则得到过分割的点集之间的融合关系；

首先定义一下主边和次边的概念，认为节点V＝{q₁，...，q_n}是有顺序的，当一条边表示的q_i，q_j满足j＝i+1时，则该边为主边，否则该边为次边；实际上就是紧挨的点集组成的边即为主边；

在构造能量函数之前要先对无向图做几步预处理：

1)在主边中，找到极大值的边作为生长种子边：

Seed＝{c_i|sign(F_main(c_i)′)′＝-2}，c_i∈ε

2)在主边中，找到函数值小于0.375的边作为截断边，作为种子边生长的最大范围：

Break＝{c_j|F_main(c_j)＜0.375}，c_j∈ε

所有点集被截断边分为N类，故能量函数设为：

$E\left(M\right)={\min }_{M\≥N}-\log \left(\frac{{\Σ}_{k=1}^M\frac{{\Σ}_{i=1}^{m_k}{\Σ}_{k=i+1}^{m_k}\frac{F(q_i,q_j)}{j-i}}{m_k-1}}{M}\right)$

其中M为融合之后得到的集合的数量，m_k为融合之后第k类中包含的点集数目，m_k-1则为第k类包含的主边数目，j-i衡量的是次边的权重程度。