本发明涉及机器视觉测量技术领域，尤其是涉及一种基于图像处理技术的水准泡气泡长度测量方法。
背景技术：
在众多工业制造生产活动中，人们常常需要对物体表面进行水平度和垂直度的检测，因此通常会借助一些带有水准泡的计量器具(例如：水平尺等)来完成此类精密测量。目前使用最为广泛的计量器具是水平尺。水平尺主要由水准泡和尺身两部分组成，其中水准泡是其检测功能的核心，它内部有刻线和气泡。在水平尺制造过程中，水准泡气泡长度是影响水准泡质量的因素之一，因此在将水准泡装配到水平尺身前需要对其中的气泡长度进行测量。
目前对水准泡气泡长度的测量仍是采用传统的人工测量方式，但由于水准泡气泡形状不是规则的矩形，导致人工测量气泡长度的方式耗时较长，工人容易疲劳且测量结果的精度较低。而且，当今世界对于这些信息的处理方式也要求要更加地智能化，尽量减少人工的干预，提高自动化程度和工作效率。但到目前为止，国内外对水准泡气泡长度的视觉测量方法的研究成果非常少。
技术实现要素：
本发明针对现有技术的不足，提供了一种基于图像处理技术的水准泡气泡长度测量方法，该方法对获取的水准泡灰度图像进行边缘检测，再对边缘检测后的图像进行形态学闭运算；对形态学闭运算操作后的水准泡图像进行轮廓检测，计算各轮廓所包围区域的面积、各轮廓最小外接矩形的长宽比以及各轮廓的占空比；通过对各轮廓所包围区域的面积、最小外接矩形的长宽比以及占空比的限制，提取水准泡中气泡的候选轮廓，并将候选轮廓进行填充放置于新图像中；对得到的气泡候选轮廓填充图像进行形态学开运算操作，将气泡区域与其它干扰区域分隔开；再对图像进行区域筛选，选取面积最大的区域并获取其最小外接矩形，通过外接矩形的长度和所得图像标定结果计算水准泡气泡物理长度。该方法可用于机器替代人工对水准泡气泡长度进行自动测量，解决了人工测量费时费力，效率低下的问题。
本发明所采用的技术方案是，一种基于图像处理技术的水准泡气泡长度测量方法，具体按照以下步骤实施：
步骤1：读取水准泡原始灰度图像并对其进行边缘检测；
步骤2：对边缘检测后的图像进行形态学闭运算操作；
步骤3：对形态学闭运算操作后的水准泡图像进行轮廓检测，计算各轮廓所包围区域的面积、各轮廓最小外接矩形的长宽比以及各轮廓的占空比；
步骤4：通过对各轮廓所包围区域的面积、最小外接矩形的长宽比以及占空比的限制，提取水准泡中气泡的候选轮廓，并将候选轮廓进行填充放置于新图像中；
步骤5：对步骤4得到的气泡候选轮廓填充图像进行形态学开运算操作，将气泡区域与其它干扰区域分隔开；
步骤6：对步骤5得到的图像进行区域筛选，选取面积最大的区域并获取其最小外接矩形；
步骤7：计算步骤6中的最大面积区域的最小外接矩形的长度；
步骤8：根据步骤7得到的外接矩形的长度和所得图像标定结果计算水准泡气泡物理长度。
本发明的有益效果是：通过图像处理技术对水准泡图像进行边缘检测、形态学闭运算、轮廓检测，根据气泡特征获取气泡候选轮廓，然后通过轮廓填充和形态学开运算将气泡区域和其它干扰区域分离，最后获取气泡最小外接矩形长度并结合图像标定结果得到气泡物理长度。本发明方法能有效解决人工测量费时费力的问题，可有效提高测量效率。
附图说明
图1是本发明方法的步骤流程图；
图2是本发明方法获取的水准泡原始灰度图像；
图3是本发明方法获取的水准泡边缘检测结果图像；
图4是本发明方法获取的形态学闭运算操作结果图像；
图5是本发明方法提取的候选气泡轮廓填充结果图像；
图6是本发明方法获取的形态学开运算操作结果图像；
图7是本发明方法获取的水准泡气泡最小外接矩形标记图像。
具体实施方式
以下结合附图对本发明作进一步说明。
如图1所示，本发明方法的步骤是：
步骤1：读取水准泡原始灰度图像并对其进行边缘检测
读取的水准泡原始灰度图像isrc如图2所示，利用canny边缘检测对原始灰度图像isrc进行边缘检测，得到水准泡边缘检测结果图像icanny，如图3所示。
步骤2：对边缘检测后的图像进行形态学闭运算操作
采用尺寸为5×5的椭圆形结构元素对水准泡边缘检测结果图像icanny进行形态学闭运算处理，得到形态学闭运算操作结果图像iclose如图4所示。
步骤3：对形态学闭运算操作后的水准泡图像进行轮廓检测，计算各轮廓所包围区域的面积、各轮廓最小外接矩形的长宽比以及各轮廓的占空比
对形态学闭运算操作结果图像iclose进行轮廓检测，计算每个轮廓ci的包围区域面积si、每个轮廓最小外接矩形的宽高比ki以及每个轮廓的占空比ri，即轮廓包围区域面积与其外接矩形面积之比，其中i＝1,2,3,…,n，n为在图像iclose中检测到的轮廓总数。
轮廓ci的最小外接矩形宽高比ki通过式(1)求得：
ki＝max(wi，hi)/min(wi,hi)(1)
其中wi和hi分别为轮廓ci最小外接矩形的宽度和高度，max()表示取较大值，min()表示取较小值。
步骤4：通过对各轮廓所包围区域的面积、最小外接矩形的宽高比以及占空比的限制，提取水准泡中气泡的候选轮廓，并将候选轮廓进行填充放置于新图像中
对步骤3得到的各轮廓包围区域面积si、最小外接矩形的宽高比ki以及占空比ri进行限制，当轮廓包围面积si在某个固定范围[sl,sh]内，最小外接矩形的宽高比ki小于某一阈值kl，占空比ri大于某一阈值rh时，即；
sl<si<sh(2)
ki<kl(3)
ri＞rh(4)
则判定该轮廓为水准泡气泡的候选轮廓；所有轮廓判定后，将气泡的候选轮廓进行内部填充，即将轮廓内部像素的灰度值设置为255，并放置于一幅背景像素值均为0的新图像中，得到候选气泡轮廓填充结果图像icb，如图5所示。
步骤5：对步骤4得到的气泡候选轮廓填充图像进行形态学开运算操作，将气泡区域与其它干扰区域分隔开
采用寸为7×7的椭圆形结构元素对候选气泡轮廓填充图像icb进行形态学开运算处理，将气泡区域与其它可能存在的干扰区域分隔开，得到形态学开运算操作结果图像iopen，如图6所示。
步骤6：对步骤5得到的图像进行区域筛选，选取面积最大的区域并获取其最小外接矩形
对形态学开运算操作结果图像iopen进行轮廓检测，选取面积最大的区域amax，此区域即为水准泡气泡区域，获取此区域的最小外接矩形rmer。将最小外接矩形rmer绘制于水准泡原始灰度图像中，结果如图7所示。
获取面积最大区域时，可以用区域像素个数作为面积计算，也可以用区域对应的最小外接矩形面积计算。
步骤7：计算步骤6中的最大面积区域的最小外接矩形的长度
根据步骤6中得到的最小外接矩形rmer，按式(5)计算最小接外接矩形rmer的长度li：
li＝max(wr,hr)(5)
其中，wr和hr分别为最小外接矩形rmer的宽度和高度，max()表示取较大值。
步骤8：根据步骤7得到的外接矩形的长度和所得图像标定结果计算水准泡气泡物理长度
根据步骤7得到的像素平面的水准泡气泡最小外接矩形的长度li，单位为pixel，和对图像进行标定获得的标定结果res，单位为mm/pixel，可计算得到水准泡气泡物理长度lr，单位为mm，即：
lr＝li×res(6)。
1.一种基于图像处理技术的水准泡气泡长度测量方法，其特征在于：
具体按照以下步骤实施：
步骤1：读取水准泡原始灰度图像并对其进行边缘检测；
步骤2：对边缘检测后的图像进行形态学闭运算操作；
步骤3：对形态学闭运算操作后的水准泡图像进行轮廓检测，计算各轮廓所包围区域的面积、各轮廓最小外接矩形的宽高比以及各轮廓的占空比；
步骤4：通过对各轮廓所包围区域的面积、最小外接矩形的宽高比以及占空比的限制，提取水准泡中气泡的候选轮廓，并将候选轮廓进行填充放置于新图像中；
步骤5：对步骤4得到的气泡候选轮廓填充图像进行形态学开运算操作，将气泡区域与其它干扰区域分隔开；
步骤6：对步骤5得到的图像进行区域筛选，选取面积最大的区域并获取其最小外接矩形；
步骤7：计算步骤6中的最大面积区域的最小外接矩形的长度；
步骤8：根据步骤7得到的外接矩形的长度和所得图像标定结果计算水准泡气泡物理长度；
所述的步骤3具体是：
对形态学闭运算操作后获得的图像iclose进行轮廓检测，计算每个轮廓ci的包围区域面积si、每个轮廓最小外接矩形的宽高比ki以及每个轮廓的占空比ri，即轮廓包围区域面积与其外接矩形面积之比，其中i＝1,2,3,…,n；
所述的步骤4具体是：
对步骤3得到的各轮廓包围区域面积si、最小外接矩形的宽高比ki以及占空比ri进行限制，当轮廓包围面积si在某个固定范围[sl,sh]内，最小外接矩形的宽高比ki小于某一阈值kl，占空比ri大于某一阈值rh时，即；
sl<si<sh
ki<kl
ri＞rh
则判定该轮廓为水准泡气泡的候选轮廓；所有轮廓判定后，将气泡的候选轮廓进行内部填充并放置于一幅新图像icb中；
所述的步骤8具体是：
根据步骤7得到的像素平面的水准泡气泡最小外接矩形的长度li，单位为pixel，和对图像进行标定获得的标定结果res，单位为mm/pixel，可计算得到水准泡气泡物理长度度lr，单位为mm，即：
lr＝li×res。
2.根据权利要求1所述的基于图像处理技术的水准泡气泡长度测量方法，其特征在于：
所述的步骤2具体是：采用尺寸为5×5的椭圆形结构元素对边缘检测后的图像进行形态学闭运算处理。
3.根据权利要求1所述的基于图像处理技术的水准泡气泡长度测量方法，其特征在于：
所述的步骤5具体是：采用尺寸为7×7的椭圆形结构元素对气泡候选轮廓填充图像进行形态学开运算处理，将气泡区域与其它干扰区域分隔开。
4.根据权利要求1所述的基于图像处理技术的水准泡气泡长度测量方法，其特征在于：
所述的步骤6中获取面积最大区域时，可以用区域像素个数作为面积计算，也可以用区域对应的最小外接矩形面积计算。
技术总结