本发明涉及疲劳端口形貌分析领域，具体涉及一种基于漫水填充算法的疲劳断口形貌特征提取方法。
背景技术：
材料的疲劳断口分析是将样品断口表面腐蚀清洗后分别进行宏观和微观的分析。宏观分析是用肉眼、放大镜和实体显微镜对断裂零件进行直接观察与分析；微观分析是对整体断口成像分析的方法主要为使用透射电镜和扫描电镜。然而，宏观分析的放大倍数通常为100倍以下，很难用肉眼判定、分析出裂纹的微观构成；而通过电镜得到的断口图像仅经过放大和灰度处理，最低放大倍数1000倍，分辨率最高只可达0.01um且大都通过人为的观察与断定。
疲劳裂纹分为表面和内部萌生，对于某些复杂的断口形貌如塑性变形的累积导致难以分辨出裂纹源与塑性累积区域的差别、存在的凹陷容易被误认为是裂纹源、裂纹从萌生区形核在一个疲劳断口形貌中常常存在多种扩展机制，使得裂纹从源头出发沿着各个方向扩展，有时难以观察得到裂纹的扩展走势等，这些都难以通过肉眼直观的观察得到和快速确定。
疲劳断口解理面的数量确定需要人为计量、面积大小计算需要先提取解理面轮廓再通过计算机积分计算、占比的统计分布需要提取整个裂纹扩展低速区和各个解理面轮廓计算面积大小后得到解理面的占比，这些过程需要频繁的步骤和计算。
技术实现要素：
针对现有技术中的上述不足，本发明提供的一种基于漫水填充算法的疲劳断口形貌特征提取方法提供了一种快速获取疲劳断口形貌特征的方法。
为了达到上述发明目的，本发明采用的技术方案为：
提供一种基于漫水填充算法的疲劳断口形貌特征提取方法，其包括以下步骤：
s1、采集疲劳断口图像并采用matlab进行预处理，得到预处理后的图像；
s2、对预处理后的图像进行漫水填充分割编号，得到序列区域；
s3、绘制序列区域的图像轮廓，得到具有闭合轮廓的图像；获取序列区域的图像特征值，图像特征值包括熵、灰度平均值和灰度方差；
s4、计算具有闭合轮廓的图像的区域面积；
s5、根据区域面积计算各个轮廓的强度因子；
s6、基于图像特征值和强度因子进行线性拟合，得到拟合曲线；
s7、根据具有闭合轮廓的图像和拟合曲线获取疲劳断口形貌特征。
进一步地，步骤s1的具体方法为：
通过数码相机采集疲劳断口图像，并采用matlab进行灰度化、去噪和二值化处理，得到预处理后的图像；
或，通过扫描电镜或透射电镜进行图像采集与灰度处理，并采用matlab进行去噪和二值化处理，得到预处理后的图像。
进一步地，步骤s2的具体方法包括以下子步骤：
s2-1、通过vs软件调用floodfill算子进行空域、渐变或固定范围、渐变或浮动范围、操作标志符低八位的四连通域范围，以及浮动范围、操作标志符低八位的八连通域范围漫水填充得到预处理后的图像的连通区域图像；
s2-2、以连通区域图像的初始像素点为起点，根据灰度的梯度变化依次进行区域分割和序列编号，得到序列区域。
进一步地，步骤s3中绘制序列区域的图像轮廓，得到具有闭合轮廓的图像的具体方法为：
采用vs软件的mask函数查找序列区域边界处的灰度值，根据该灰度值连通边界像素点形成闭合轮廓，得到具有闭合轮廓的图像。
进一步地，步骤s3中获取序列区域的熵的具体方法为：
在vs系列软件中加载opencv库函数，对序列区域运用全局阈值法搜索每个序列区域内全部像素点，调用doubleentropy函数计算得到每个序列区域图像的熵。
进一步地，步骤s5的具体方法为：
根据公式：
计算强度因子δk；其中n为常数，当裂纹从表面萌生时该值为0.6，从内部萌生时该值为0.5；π为圆周率；areα为区域面积。
进一步地，步骤s6的具体方法为：
将各个具有闭合轮廓的图像的图像特征值作比，得到各个具有闭合轮廓的图像的对比度；采用oringin软件调用图像特征值，分别以熵、灰度平均值、灰度方差和各个具有闭合轮廓的图像的对比度作为自变量，以强度因子作为因变量，进行一次函数拟合，分别得到与熵对应的拟合曲线、与灰度平均值对应的拟合曲线、与灰度方差对应的拟合曲线，以及与各个具有闭合轮廓的图像的对比度对应的拟合曲线。
进一步地，步骤s7中疲劳断口形貌特征包括图像特征和数值特征；其中图像特征包括裂纹萌生源位置、塑性累积区域、裂纹走势和受力变形程度；数值特征包括由线性拟合得到的图像特征值与强度因子的大小关系和统计分布。
本发明的有益效果为：
1、本发明利用图像像素灰度值的敏感性，结合图像采集系统，利用漫水填充算子以及oringin数据拟合可定量获得不同疲劳寿命下裂纹的断口形貌，轮廓提取中由裂纹走势可以反推得到裂纹源头；区域面积计算可以得到不同裂纹形貌的统计占比；由不同的区域面积可以得到与之对应的裂纹强度因子。克服了以往直观的分析由采集系统得出的灰度图像肉眼辨别的困难和不准确性以及区域面积、强度因子复杂计算的问题，使分析的结论更快捷准确。
2、本发明中的线性拟合部分可以表征出断后形貌各个连通区域的组织受力应变响应情况，具有普遍应用性。
3、本发明程序设计函数调用与操作较为简便，适用性强。
附图说明
图1为本发明的流程示意图；
图2为实施例中的原始图像；
图3为操作标志符的低八位八连通的区域提取效果示意图；
图4为绘制序列区域的图像轮廓的示意图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，该基于漫水填充算法的疲劳断口形貌特征提取方法包括以下步骤：
s1、采集疲劳断口图像并采用matlab进行预处理，得到预处理后的图像；
s2、对预处理后的图像进行漫水填充分割编号，得到序列区域；
s3、绘制序列区域的图像轮廓，得到具有闭合轮廓的图像；获取序列区域的图像特征值，图像特征值包括熵、灰度平均值和灰度方差；
s4、计算具有闭合轮廓的图像的区域面积；
s5、根据区域面积计算各个轮廓的强度因子；
s6、基于图像特征值和强度因子进行线性拟合，得到拟合曲线；
s7、根据具有闭合轮廓的图像和拟合曲线获取疲劳断口形貌特征。
步骤s1的具体方法为：
通过数码相机采集疲劳断口图像，并采用matlab进行灰度化、去噪和二值化处理，得到预处理后的图像；
或，通过扫描电镜或透射电镜进行图像采集与灰度处理，并采用matlab进行去噪和二值化处理，得到预处理后的图像。
步骤s2的具体方法包括以下子步骤：
s2-1、通过vs软件调用floodfill算子进行空域、渐变或固定范围、渐变或浮动范围、操作标志符低八位的四连通域范围，以及浮动范围、操作标志符低八位的八连通域范围漫水填充得到预处理后的图像的连通区域图像；
s2-2、以连通区域图像的初始像素点为起点，根据灰度的梯度变化依次进行区域分割和序列编号，得到序列区域。
步骤s3中绘制序列区域的图像轮廓，得到具有闭合轮廓的图像的具体方法为：采用vs软件的mask函数查找序列区域边界处的灰度值，根据该灰度值连通边界像素点形成闭合轮廓，得到具有闭合轮廓的图像。
步骤s3中获取序列区域的熵的具体方法为：在vs系列软件中加载opencv库函数，对序列区域运用全局阈值法搜索每个序列区域内全部像素点，调用doubleentropy函数计算得到每个序列区域图像的熵。
步骤s5的具体方法为：根据公式：
计算强度因子δk；其中n为常数，当裂纹从表面萌生时该值为0.6，从内部萌生时该值为0.5；π为圆周率；areα为区域面积。
步骤s6的具体方法为：将各个具有闭合轮廓的图像的图像特征值作比，得到各个具有闭合轮廓的图像的对比度；采用oringin软件调用图像特征值，分别以熵、灰度平均值、灰度方差和各个具有闭合轮廓的图像的对比度作为自变量，以强度因子作为因变量，进行一次函数拟合，分别得到与熵对应的拟合曲线、与灰度平均值对应的拟合曲线、与灰度方差对应的拟合曲线，以及与各个具有闭合轮廓的图像的对比度对应的拟合曲线。
步骤s7中疲劳断口形貌特征包括图像特征和数值特征；其中图像特征包括裂纹萌生源位置、塑性累积区域、裂纹走势和受力变形程度；数值特征包括由线性拟合得到的图像特征值与强度因子的大小关系和统计分布。
在本发明的一个实施例中，原始图像如图2所示，对原始图像进行预处理后，采用操作标志符的低八位八连通的区域提取的连通区域图像如图3所示，对该连通区域图像分割和序列编号得到的序列区域进行轮廓绘制，得到的具有闭合轮廓的图像如图4所示。在具有闭合轮廓的图像进行后续操作，可以实现快速提取疲劳断开形貌特征。
在具体实施过程中，断口形貌分析：主要包含图像特征和数值特征两个部分。图像特征为：裂纹萌生源位置确定、塑性累积区域确定、裂纹走势确定以及帮助肉眼识别各种特殊形貌如“小平面”“鱼眼”“暗区”“光区”“萌生区、扩展区、瞬断区”等的功能。通过分析对提取的特征区域根据灰度值的不同可以进行序列编号，同一编号下的区域结合特征值数据反应组成成分相近，受力变形程度相当；对轮廓的提取可以分析出裂纹走向、裂纹累积处以及裂纹萌生源；数值特征主要功能为由线性拟合解释图像特征值与强度因子的大小关系和统计分部：1、由拟合曲线的统计分布可以确定熵反应的图像灰度混乱程度、灰度平均值、方差、图像区域对比度与强度因子的大小关系。2、由图像区域对比度与强度因子的大小关系进而可以根据该强度因子值所在的序列区域于整个图像上的分布做到对强度因子在整个断面的位置和大小分布。
综上所述，本发明利用图像像素灰度值的敏感性，结合图像采集系统，利用漫水填充算子以及oringin数据拟合可定量获得不同疲劳寿命下裂纹的断口形貌，轮廓提取中由裂纹走势可以反推得到裂纹源头；区域面积计算可以得到不同裂纹形貌的统计占比；由不同的区域面积可以得到与之对应的裂纹强度因子。克服了以往直观的分析由采集系统得出的灰度图像肉眼辨别的困难和不准确性以及区域面积、强度因子复杂计算的问题，使分析的结论更快捷准确。
1.一种基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，包括以下步骤：
s1、采集疲劳断口图像并采用matlab进行预处理，得到预处理后的图像；
s2、对预处理后的图像进行漫水填充分割编号，得到序列区域；
s3、绘制序列区域的图像轮廓，得到具有闭合轮廓的图像；获取序列区域的图像特征值，图像特征值包括熵、灰度平均值和灰度方差；
s4、计算具有闭合轮廓的图像的区域面积；
s5、根据区域面积计算各个轮廓的强度因子；
s6、基于图像特征值和强度因子进行线性拟合，得到拟合曲线；
s7、根据具有闭合轮廓的图像和拟合曲线获取疲劳断口形貌特征。
2.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s1的具体方法为：
通过数码相机采集疲劳断口图像，并采用matlab进行灰度化、去噪和二值化处理，得到预处理后的图像；
或，通过扫描电镜或透射电镜进行图像采集与灰度处理，并采用matlab进行去噪和二值化处理，得到预处理后的图像。
3.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s2的具体方法包括以下子步骤：
s2-1、通过vs软件调用floodfill算子进行空域、渐变或固定范围、渐变或浮动范围、操作标志符低八位的四连通域范围，以及浮动范围、操作标志符低八位的八连通域范围漫水填充得到预处理后的图像的连通区域图像；
s2-2、以连通区域图像的初始像素点为起点，根据灰度的梯度变化依次进行区域分割和序列编号，得到序列区域。
4.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s3中绘制序列区域的图像轮廓，得到具有闭合轮廓的图像的具体方法为：
采用vs软件的mask函数查找序列区域边界处的灰度值，根据该灰度值连通边界像素点形成闭合轮廓，得到具有闭合轮廓的图像。
5.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s3中获取序列区域的熵的具体方法为：
在vs系列软件中加载opencv库函数，对序列区域运用全局阈值法搜索每个序列区域内全部像素点，调用doubleentropy函数计算得到每个序列区域图像的熵。
6.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s5的具体方法为：
根据公式：
计算强度因子δk；其中n为常数，当裂纹从表面萌生时该值为0.6，从内部萌生时该值为0.5；π为圆周率；areα为区域面积。
7.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s6的具体方法为：
将各个具有闭合轮廓的图像的图像特征值作比，得到各个具有闭合轮廓的图像的对比度；采用oringin软件调用图像特征值，分别以熵、灰度平均值、灰度方差和各个具有闭合轮廓的图像的对比度作为自变量，以强度因子作为因变量，进行一次函数拟合，分别得到与熵对应的拟合曲线、与灰度平均值对应的拟合曲线、与灰度方差对应的拟合曲线，以及与各个具有闭合轮廓的图像的对比度对应的拟合曲线。
8.根据权利要求1所述的基于漫水填充算法的疲劳断口形貌特征提取方法，其特征在于，所述步骤s7中疲劳断口形貌特征包括图像特征和数值特征；其中图像特征包括裂纹萌生源位置、塑性累积区域、裂纹走势和受力变形程度；数值特征包括由线性拟合得到的图像特征值与强度因子的大小关系和统计分布。
技术总结