叶型孔图像的轮廓型线提取与分割方法及系统、设备、存储介质

摘要

本发明公开了一种叶型孔图像的轮廓型线提取与分割方法及系统、设备、存储介质，所述叶型孔图像的轮廓型线提取与分割方法通过对叶型孔图像进行滤波降噪、区域划分、目标区域确定、像素级边缘提取、边缘平滑、两次轮廓型线分割与合并、轮廓型线提取等步骤，获取了完整的叶型孔轮廓型线并进一步将其分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，为后续的叶型孔轮廓型线参数计算奠定基础，保证检测的准确性。其中，通过第一次轮廓型线分割可以将较大直径的圆弧段有效地分割出来，而通过第二次轮廓型线分割可以将较小直径的圆弧段分割出来，并且可以细化较大直径的圆弧段末端，可以实现更加精确的轮廓型线分割，并达到较高的算法效率。

权利要求

1.一种叶型孔图像的轮廓型线提取与分割方法，其特征在于，包括以下内容：对采集到的叶型孔图像进行滤波降噪处理；

将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

提取出目标区域的像素级边缘；

对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

2.如权利要求1所述的叶型孔图像的轮廓型线提取与分割方法，其特征在于，所述对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线的过程具体包括以下内容：采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

3.如权利要求2所述的叶型孔图像的轮廓型线提取与分割方法，其特征在于，所述采用圆弧方程对相邻的两条轮廓段进行第一次拟合以得到一条圆弧段的过程具体包括以下内容：设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数目，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

2 2 2

(u‑u0) +(v‑v0) ＝R，u∈[u1,un]且v∈[v1,vn]其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

4.如权利要求1所述的叶型孔图像的轮廓型线提取与分割方法，其特征在于，所述对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线的过程具体包括以下内容：以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

5.一种叶型孔图像的轮廓型线提取与分割系统，其特征在于，包括：滤波模块，用于对采集到的叶型孔图像进行滤波降噪处理；

区域划分模块，用于将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

区域分解模块，用于将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

边缘提取模块，用于提取出目标区域的像素级边缘；

边缘处理模块，用于对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

轮廓型线分割模块，用于对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

轮廓型线提取模块，用于对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

6.如权利要求5所述的叶型孔图像的轮廓型线提取与分割系统，其特征在于，所述轮廓型线分割模块具体包括：第一分割单元，用于采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

第一拟合单元，用于采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

第二分割单元，用于采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

第二拟合单元，用于采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

7.如权利要求6所述的叶型孔图像的轮廓型线提取与分割系统，其特征在于，第一拟合单元和第二拟合单元采用圆弧方程对相邻两条轮廓段进行拟合得到圆弧段的过程包括以下内容：设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数目，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

2 2 2

(u‑u0) +(v‑v0) ＝R，u∈[u1,un]且v∈[v1,vn]其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

8.如权利要求5所述的叶型孔图像的轮廓型线提取与分割系统，其特征在于，所述轮廓型线提取模块包括：前缘曲线提取单元，用于以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

后缘曲线提取单元，以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

叶盆和叶背曲线提取单元，用于将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

9.一种设备，其特征在于，包括处理器和存储器，所述存储器中存储有计算机程序，所述处理器通过调用所述存储器中存储的所述计算机程序，用于执行如权利要求1～4任一项所述的方法的步骤。

10.一种计算机可读取的存储介质，用于存储对叶型孔图像进行轮廓型线提取与分割的计算机程序，其特征在于，所述计算机程序在计算机上运行时执行如权利要求1～4任一项所述的方法的步骤。

说明书

叶型孔图像的轮廓型线提取与分割方法及系统、设备、存储

介质

技术领域

[0001] 本发明涉及图像处理技术领域，特别地，涉及一种叶型孔图像的轮廓型线提取与分割方法及系统、设备、计算机可读取的存储介质。

背景技术

[0002] 在中小型航空发动机中的静子安装环类零件上，沿其圆柱壁面的周向方向均匀分布着多个安装和固定导流叶片用的叶型孔特征，这些叶型孔在整台发动机中有数百个之多，而且形状复杂、规格较多。在发动机的装配过程中，需要将导流叶片插装并钎焊到相应的叶型孔中，因而叶型孔的轮廓型线参数就决定了叶片能否插装到位以及叶片钎焊后的位置和姿态，从而影响到发动机装配的可行性、一致性和互换性等。

[0003] 当前，随着我国航空发动机性能的不断提升，对于批量叶型孔的轮廓度、位置度和均布度等轮廓型线参数的测量精度和检测效率提出了越来越高的要求，迫切需要研制出新型、高效的自动化测量设备。近年来，随着测量技术及相关学科的飞速发展，基于图像技术的机器视觉测量作为一项高新技术，已经广泛应用于航空工业的各个领域中。视觉测量是以图像作为检测手段，将机器视觉与图像处理技术应用于几何特征的测量和定位，不仅具有非接触式测量方法的优点，还具有成本低、操作简便、机动灵活和信息丰富等独特优点，因而被广泛应用于工业现场，这就为叶型孔特征的高精高效测量及相关设备研发提供了一项技术解决方案和思路。

[0004] 通常说来，叶型孔的轮廓型线是由四段不同参数和类型的弧线段组成的，即前缘、后缘、叶盆和叶背。其中，前缘和后缘为两个圆弧或椭圆弧曲线，而叶盆和叶背则为两条平面自由曲线。在基于机器视觉测量原理并采用背向照明方式采集到被测叶型孔的背向照明图像后，为了最终解算出和获取到该叶型孔的轮廓型线参数，从而判断该叶型孔的加工质量，还需要通过一系列后续的图像处理过程对叶型孔图像进行处理，从而提取出完整的叶型孔轮廓型线，并进一步将其分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，以便于后续对于各段曲线的分别处理。然而，由于叶型孔轮廓型线的复杂性及非标准性，尚未有有效且可靠的解决方法。

发明内容

[0005] 本发明提供了一种叶型孔图像的轮廓型线提取与分割方法及系统、设备、计算机可读取的存储介质，以解决现有技术存在的上述缺点。

[0006] 根据本发明的一个方面，提供一种叶型孔图像的轮廓型线提取与分割方法，包括以下内容：

[0007] 对采集到的叶型孔图像进行滤波降噪处理；

[0008] 将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

[0009] 将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

[0010] 提取出目标区域的像素级边缘；

[0011] 对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

[0012] 对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

[0013] 对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

[0014] 进一步地，所述对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线的过程具体包括以下内容：

[0015] 采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

[0016] 采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

[0017] 采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

[0018] 采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

[0019] 进一步地，所述采用圆弧方程对相邻的两条轮廓段进行第一次拟合以得到一条圆弧段的过程具体包括以下内容：

[0020] 设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数目，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

[0021] (u‑u0)2+(v‑v0)2＝R2，u∈[u1,un]且v∈[v1,vn]

[0022] 其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：

[0023]

[0024] 基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

[0025] 进一步地，所述对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线的过程具体包括以下内容：

[0026] 以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

[0027] 以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

[0028] 将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

[0029] 另外，本发明还提供一种叶型孔图像的轮廓型线提取与分割系统，包括：

[0030] 滤波模块，用于对采集到的叶型孔图像进行滤波降噪处理；

[0031] 区域划分模块，用于将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

[0032] 区域分解模块，用于将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

[0033] 边缘提取模块，用于提取出目标区域的像素级边缘；

[0034] 边缘处理模块，用于对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

[0035] 轮廓型线分割模块，用于对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

[0036] 轮廓型线提取模块，用于对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

[0037] 进一步地，所述轮廓型线分割模块具体包括：

[0038] 第一分割单元，用于采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

[0039] 第一拟合单元，用于采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

[0040] 第二分割单元，用于采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

[0041] 第二拟合单元，用于采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

[0042] 进一步地，第一拟合单元和第二拟合单元采用圆弧方程对相邻两条轮廓段进行拟合得到圆弧段的过程包括以下内容：

[0043] 设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数目，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

[0044] (u‑u0)2+(v‑v0)2＝R2，u∈[u1,un]且v∈[v1,vn]

[0045] 其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：

[0046]

[0047] 基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

[0048] 进一步地，所述轮廓型线提取模块包括：

[0049] 前缘曲线提取单元，用于以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

[0050] 后缘曲线提取单元，以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

[0051] 叶盆和叶背曲线提取单元，用于将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

[0052] 另外，本发明还提供一种设备，包括处理器和存储器，所述存储器中存储有计算机程序，所述处理器通过调用所述存储器中存储的所述计算机程序，用于执行如上所述的方法的步骤。

[0053] 另外，本发明还提供一种计算机可读取的存储介质，用于存储对叶型孔图像进行轮廓型线提取与分割的计算机程序，所述计算机程序在计算机上运行时执行如上所述的方法的步骤。

[0054] 本发明具有以下效果：

[0055] 本发明的叶型孔图像的轮廓型线提取与分割方法，通过对叶型孔图像进行滤波降噪、区域划分、目标区域确定、像素级边缘提取、边缘平滑、两次轮廓型线分割与合并、轮廓型线提取等步骤，获取了完整的叶型孔轮廓型线并进一步将其分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，作为复杂及非标准的异形孔图像数据转换为数字化评价的前提技术，为后续的叶型孔轮廓型线参数计算奠定基础，保证检测的准确性。本发明采用两次轮廓型线分割和合并过程来获取最终的分割结果，通过第一次轮廓型线分割可以将较大直径的圆弧段有效地分割出来，而通过第二次轮廓型线分割可以将较小直径的圆弧段分割出来，并且可以细化较大直径的圆弧段末端，可以实现更加精确的轮廓型线分割，并达到较高的算法效率。

[0056] 另外，本发明的叶型孔图像的轮廓型线提取与分割系统、设备、计算机可读取的存储介质同样具有上述优点。

[0057] 除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照图，对本发明作进一步详细的说明。

附图说明

[0058] 构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

[0059] 图1是本发明优选实施例的叶型孔图像的轮廓型线提取与分割方法的流程示意图。

[0060] 图2是本发明优选实施例中的叶型孔的背向照明图像的示意图。

[0061] 图3是本发明优选实施例中的从前景区域中分解出的目标区域图像的示意图。

[0062] 图4是本发明优选实施例中提取出的像素级边缘图像的示意图。

[0063] 图5是图1中步骤S6的子流程示意图。

[0064] 图6是本发明优选实施例中的第一次轮廓型线分割结果图像的示意图。

[0065] 图7是本发明优选实施例中的第一次轮廓段合并结果图像的示意图。

[0066] 图8是图1中步骤S7的子流程示意图。

[0067] 图9是本发明优选实施例中的前缘曲线和后缘曲线的最小外接矩形图像的示意图。

[0068] 图10是本发明优选实施例中的轮廓型线分割结果图像的示意图。

[0069] 图11是本发明另一实施例中的叶型孔图像的轮廓型线提取与分割系统的模块结构示意图。

具体实施方式

[0070] 以下结合附图对本发明的实施例进行详细说明，但是本发明可以由下述所限定和覆盖的多种不同方式实施。

[0071] 如图1所示，本发明的优选实施例提供一种叶型孔图像的轮廓型线提取与分割方法，包括以下内容：

[0072] 步骤S1：对采集到的叶型孔图像进行滤波降噪处理；

[0073] 步骤S2：将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

[0074] 步骤S3：将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

[0075] 步骤S4：提取出目标区域的像素级边缘；

[0076] 步骤S5：对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

[0077] 步骤S6：对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

[0078] 步骤S7：对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

[0079] 可以理解，本实施例的叶型孔图像的轮廓型线提取与分割方法，通过对叶型孔图像进行滤波降噪、区域划分、目标区域确定、像素级边缘提取、边缘平滑、两次轮廓型线分割与合并、轮廓型线提取等步骤，获取了完整的叶型孔轮廓型线并进一步将其分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，作为复杂及非标准的异形孔图像数据转换为数字化评价的前提技术，为后续的叶型孔轮廓型线参数计算奠定基础，保证检测的准确性。本发明采用两次轮廓型线分割和合并过程来获取最终的分割结果，通过第一次轮廓型线分割可以将较大直径的圆弧段有效地分割出来，而通过第二次轮廓型线分割可以将较小直径的圆弧段分割出来，并且可以细化较大直径的圆弧段末端，可以实现更加精确的轮廓型线分割，并达到较高的算法效率。

[0080] 可以理解，具体采用背向照明方式采集被测叶型孔的背向照明图像，被测叶型孔的背向照明图像如图2所示，在叶型孔的背向照明图像中，遮光的部分为黑色，透光的部分为白色，从而使叶型孔图像中的像素黑白分明而易于后续处理和分析。在所述步骤S1中，对采集到的被测叶型孔的背向照明图像进行中值滤波处理，以降低原始图像中的噪声干扰。

[0081] 可以理解，在所述步骤S2中，具体通过最大类间方差法确定图像阈值，并应用该阈值将经过中值滤波的叶型孔图像划分为前景区域和背景区域，其中，将灰度值高于阈值的像素所构成的区域设定为前景区域，而将灰度值低于阈值的像素所构成的区域设定为背景区域，前景区域即为叶型孔区域，即图2中白色区域。其中，最大类间方差法属于现有算法，故在此不再赘述。

[0082] 可以理解，在所述步骤S3中，将前景区域进一步分解为多个互不相连的子区域，并选择出其中面积最大的子区域作为后续处理的目标区域，该目标区域是叶型孔图像中的一个完整的叶型孔区域，目标区域图像如图3所示。

[0083] 可以理解，在所述步骤S4中，具体采用基于图像灰度梯度的边缘检测算法提取出目标区域的像素级边缘，其中，具体的边缘检测算法属于现有技术，故在此不再赘述。从目标区域图像中提取出来的像素级边缘图像如图4所示。

[0084] 可以理解，在所述步骤S5中，将像素级边缘上的每个像素均投影到该像素的局部回归线上，从而对像素级边缘进行平滑处理，以抑制像素级边缘上的异常值，获取到叶型孔的亚像素级边缘作为该叶型孔的完整轮廓型线。

[0085] 可以理解，如图5所示，所述步骤S6具体包括以下内容：

[0086] 步骤S61：采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

[0087] 步骤S62：采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

[0088] 步骤S63：采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

[0089] 步骤S64：采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

[0090] 具体地，先采用Ramer算法(ramer‑douglas‑peucker，拉默‑道格拉斯‑普克算法，又简称RDP算法)并设定第一距离阈值D1，对叶型孔的完整轮廓型线进行第一次分割，将其分割为若干条轮廓段，每条轮廓段由多个数据点组成，然后将每条轮廓段用一条直线段来替代，该直线段由被替代的轮廓段的起点与终点直接连接而成，每条直线段与相应轮廓段之间的最大距离均小于D1。其中，具体的Ramer算法对轮廓型线进行分割的具体过程属于现有技术，故在此不再赘述。其中，第一次轮廓型线分割结果的图像如图6所示。

[0091] 然后进行第一次轮廓段合并，采用圆弧方程对相连两条轮廓段的数据点进行第一次拟合，得到一条拟合的圆弧段。具体的拟合过程如下：

[0092] 设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数量，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

[0093] (u‑u0)2+(v‑v0)2＝R2，u∈[u1,un]且v∈[v1,vn]

[0094] 其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：

[0095]

[0096] 基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

[0097] 如果拟合得到的圆弧段与这两条轮廓段之间的最大距离不大于相应的两条直线段与这两条轮廓段之间的最大距离，则采用拟合得到的圆弧段来替代这两条轮廓段，由此完成了两条轮廓段的合并；如果拟合圆弧段与这两条轮廓段之间的最大距离大于相应的两条直线段与这两条轮廓段之间的最大距离，则依然采用这两条直线段来替代相应的两条轮廓段。不断迭代计算一直到不再发生任何改变为止，迭代过程结束后，整个叶型孔轮廓型线被分割为多条轮廓段，每条轮廓段用直线段或者圆弧段来替代。其中，第一次轮廓段合并结果的图像如图7所示。

[0098] 接着，进行第二次轮廓型线分割，再次采用Ramer算法并设定第二距离阈值D2(D2＜D1)，对上述迭代过程中得到的用直线段替代的那些轮廓段再次进行分割，将这些轮廓段再次分割为新的轮廓段，并将每条新的轮廓段用一条直线段来替代，该直线段由被替代的轮廓段的起点与终点直接连接而成，每条直线段与相应轮廓段之间的最大距离均小于D2。

[0099] 最后，进行第二次轮廓段合并，再次采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，具体的迭代拟合过程与第一次相同，故在此不再赘述，最终得到由多条圆弧段拟合而成的轮廓型线。

[0100] 可以理解，基于叶型孔轮廓型线的特点和设计要求，前缘曲线和后缘曲线是连接叶盆曲线和叶背曲线的两段圆弧或椭圆弧，而叶盆曲线和叶背曲线则为两段自由曲线。前缘曲线和后缘曲线同时与叶盆曲线和叶背曲线相切，而且前缘曲线和后缘曲线的曲率半径较小，而叶盆曲线和叶背曲线的曲率半径较大，两类曲线之间的曲率差异明显。而且，对于前缘曲线和后缘曲线而言，应用圆弧段替代轮廓段的误差会比应用若干直线段进行替代的误差要小得多，因而在轮廓段合并过程中会达到良好的合并效果，并通过迭代将多个直线段逐步合并为一个圆弧段，从而将完整的前缘曲线和后缘曲线分割出来。因此，基于前缘曲线和后缘曲线自身的圆弧或椭圆弧特性，通过两次轮廓型线分割和合并后，所得到的多条圆弧段中必然有前缘轮廓段和后缘轮廓段。通过第一次轮廓型线分割可以将较大直径的圆弧段有效地分割出来，而通过第二次轮廓型线分割可以将较小直径的圆弧段分割出来，并且可以细化较大直径的圆弧段末端，可以实现更加精确的轮廓型线分割，并达到较高的算法效率。

[0101] 可以理解，如图8所示，所述步骤S7具体包括以下内容：

[0102] 步骤S71：以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

[0103] 步骤S72：以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

[0104] 步骤S73：将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

[0105] 具体地，如图6所示，在叶型孔图像中，坐标系的原点在图像的左上方，u轴(即列坐标)正方向朝右，v轴(即行坐标)正方向朝下。因此，前缘曲线的行坐标值是最大的，而后缘曲线的列坐标值是最大的。如图9所示，通过以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线，以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线。然后，将前缘曲线和后缘曲线从所有轮廓段中分离出来，再将剩余的相邻轮廓段进行联合，由于叶背段位于叶盆段的上方，可以基于每个轮廓段的数据点的坐标信息可以得到其属于叶盆段还是叶背段，从而提取出叶盆曲线和叶背曲线。通过上述步骤，最终将提取到的完整的叶型孔轮廓型线分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，从而为后续的叶型孔轮廓型线参数计算奠定坚实基础。其中，将轮廓型线分割结果图像如图10所示。

[0106] 另外，如图11所示，本发明的另一实施例还提供一种叶型孔图像的轮廓型线提取与分割系统，优选采用上述实施例的方法，所述系统包括：

[0107] 滤波模块，用于对采集到的叶型孔图像进行滤波降噪处理；

[0108] 区域划分模块，用于将经过滤波降噪处理后的叶型孔图像划分为前景区域和背景区域，前景区域为叶型孔区域；

[0109] 区域分解模块，用于将前景区域分解为多个互不相连的子区域，并选择面积最大的子区域作为目标区域，所述目标区域为一个完整的叶型孔；

[0110] 边缘提取模块，用于提取出目标区域的像素级边缘；

[0111] 边缘处理模块，用于对像素级边缘进行平滑处理，平滑处理获得的叶型孔的亚像素级边缘作为叶型孔的完整轮廓型线；

[0112] 轮廓型线分割模块，用于对叶型孔的完整轮廓型线进行两次轮廓型线分割和合并，得到由多条圆弧段拟合而成的轮廓型线；

[0113] 轮廓型线提取模块，用于对轮廓型线中的多条圆弧段进行整理和排序，分别提取出前缘曲线、后缘曲线、叶盆曲线和叶背曲线。

[0114] 可以理解，本实施例的叶型孔图像的轮廓型线提取与分割系统，通过对叶型孔图像进行滤波降噪、区域划分、目标区域确定、像素级边缘提取、边缘平滑、两次轮廓型线分割与合并、轮廓型线提取等步骤，获取了完整的叶型孔轮廓型线并进一步将其分割为前缘曲线、后缘曲线、叶盆曲线和叶背曲线，作为复杂及非标准的异形孔图像数据转换为数字化评价的前提技术，为后续的叶型孔轮廓型线参数计算奠定基础，保证检测的准确性。本发明采用两次轮廓型线分割和合并过程来获取最终的分割结果，通过第一次轮廓型线分割可以将较大直径的圆弧段有效地分割出来，而通过第二次轮廓型线分割可以将较小直径的圆弧段分割出来，并且可以细化较大直径的圆弧段末端，可以实现更加精确的轮廓型线分割，并达到较高的算法效率。

[0115] 可以理解，所述轮廓型线分割模块具体包括：

[0116] 第一分割单元，用于采用Ramer算法并设定第一距离阈值对叶型孔的完整轮廓型线进行第一次分割，将其分割成若干条轮廓段，并将每条轮廓段用一条直线段来替代；

[0117] 第一拟合单元，用于采用圆弧方程对相邻两条轮廓段的数据点进行第一次拟合以得到一条圆弧段，若拟合得到的圆弧段与这两条轮廓段之间的最大距离小于等于相应的两条直线段与这两条轮廓段之间的最大距离，则采用圆弧段替代相邻的两条轮廓段，否则，仍然采用两条直线段来替代相应的两条轮廓段，不断迭代，整个叶型孔轮廓型线被分割成多条轮廓段，每条轮廓段用直线段或者圆弧段替代；

[0118] 第二分割单元，用于采用Ramer算法并设定第二距离值对第一次拟合过程中用直线段替代的那些轮廓段进行第二次分割，将其分割为新的轮廓段，并将每条新的轮廓段用直线段来替代，其中，第二距离值小于第一距离值；

[0119] 第二拟合单元，用于采用圆弧方程对第二次分割得到的新的轮廓段进行第二次迭代拟合，最终得到由多条圆弧段拟合而成的轮廓型线。

[0120] 其中，第一拟合单元和第二拟合单元采用圆弧方程对相邻两条轮廓段进行拟合得到圆弧段的过程包括以下内容：

[0121] 设相邻两条轮廓段上的数据点的坐标分别为(u1,v1)、(u2,v2)、…、(un,vn)，其中，(u1,v1)和(un,vn)分别为起点坐标和终点坐标，n为两条轮廓段上的数据点数目，则在二维平面坐标系o‑uv中，圆弧方程的表达式为：

[0122] (u‑u0)2+(v‑v0)2＝R2，u∈[u1,un]且v∈[v1,vn]

[0123] 其中，u0、v0和R为未知数，(u0,v0)为圆弧段所在圆周的圆心坐标，R为圆弧所在圆周的半径，由数据点(u1,v1)、(u2,v2)、…、(un,vn)和圆弧方程的表达式构成如下目标函数f：

[0124]

[0125] 基于最小二乘法求解出使目标函数f取得极小值的u0、v0和R，从而得到由相邻的两条轮廓段上的数据点拟合而成的圆弧段。

[0126] 另外，所述轮廓型线提取模块包括：

[0127] 前缘曲线提取单元，用于以各个圆弧段的最小外接矩形左上角的行坐标按升序排列，则排在最后的圆弧段即为前缘曲线；

[0128] 后缘曲线提取单元，以各个轮廓段的最小外接矩形右上角的列坐标按升序排列，则排在最后的圆弧段为后缘曲线；

[0129] 叶盆和叶背曲线提取单元，用于将前缘曲线和后缘曲线从轮廓型线中分离出来，将剩余的相邻轮廓段进行联合，得到叶盆曲线和叶背曲线。

[0130] 可以理解，本实施例的系统中的各个模块和单元分别与上述方法实施例中的各个步骤相对应，故各个模块和单元的具体工作过程在此不再赘述，参考上述方法实施例即可。

[0131] 另外，本发明的另一实施例还提供一种设备，包括处理器和存储器，所述存储器中存储有计算机程序，所述处理器通过调用所述存储器中存储的所述计算机程序，用于执行如上所述的方法的步骤。

[0132] 另外，本发明的另一实施例还提供一种计算机可读取的存储介质，用于存储对叶型孔图像进行轮廓型线提取与分割的计算机程序，所述计算机程序在计算机上运行时执行如上所述的方法的步骤。

[0133] 一般计算机可读取存储介质的形式包括：软盘(floppy disk)、可挠性盘片(flexible disk)、硬盘、磁带、任何其与的磁性介质、CD‑ROM、任何其余的光学介质、打孔卡片(punch cards)、纸带(paper tape)、任何其余的带有洞的图案的物理介质、随机存取存储器(RAM)、可编程只读存储器(PROM)、可抹除可编程只读存储器(EPROM)、快闪可抹除可编程只读存储器(FLASH‑EPROM)、其余任何存储器芯片或卡匣、或任何其余可让计算机读取的介质。指令可进一步被一传输介质所传送或接收。传输介质这一术语可包含任何有形或无形的介质，其可用来存储、编码或承载用来给机器执行的指令，并且包含数字或模拟通信信号或其与促进上述指令的通信的无形介质。传输介质包含同轴电缆、铜线以及光纤，其包含了用来传输计算机数据信号的总线的导线。

[0134] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。