本发明涉及图像处理领域，尤其涉及一种水泥刻槽路面图像降噪增强与裂缝特征提取方法。
背景技术：
一直以来，公路交通与我们的生活息息相关，它不仅保障了我们正常的工作、学习，同时也为各类出行提供了方便，极大的推进了社会的发展。截至2017年7月，我国高速公路通车里程达13.1万公里，但随着公路总里程数目的增长，路面防护工作变得愈来愈重要。早期建设的公路，由于持续受到来自自然环境及车辆负载等因素的影响，不可避免地会受到损伤，因此，为避免损伤扩大化，及时制定出相应的路面防护策略显得尤为重要。
道路破损程度作为公路健康状况的重要指标，其中最常见的破损类型是路面裂缝，及时检测出路面裂缝并加以修补是进行路面防护的一项重要工作。水泥混凝土路面作为公路路面的一种抗滑构造，其通过添加刻槽增大道路表面的摩擦系数。刻槽的存在虽然极大地降低了交通意外的发生，但与此同时也对路面裂缝的检测工作造成了诸多不便。早期因路面检测技术未臻完善，普遍使用人工方法进行检测。在交通量较小的公路，可采用人工手段进行检测，但是在车辆来往频繁的公路采用此方法则会引起交通堵塞。此外，加之人的主观意识，人工检测的质量很难得到保证，为此，还需耗费一定的人力物力对检测人员进行道路健康状况相关知识的专业培训。
自进入20世纪以来，随着计算机技术的飞速发展，计算机图形图像处理技术、人工智能以及各种存储技术等在路面裂缝检测方面得到了普遍应用。因此，研究水泥刻槽路面裂缝检测方法，快速统计路面病害，将大大加快路面防护自动化技术的进步，同时对公路管理部门制定科学养护决策也具有重大意义。
现阶段裂缝提取方法在整个裂缝的检测工作中或多或少都会存在以下问题，如在图像对比度增强方面：部分方法采用传统的局部自适应对比度增强算法虽然能够使图像的细节信息更加突出，但并未解决图像亮度不均带来的裂缝识别精度偏低问题；在干扰噪声的去除方面：部分方法采用平移不变shearlet变换去噪方法虽然能在一定程度上解决伪gibbs效应带来的图像边缘模糊问题，但是对低信噪比路面图像仍然无法完整的保留图像中裂缝的边缘信息；在路面刻槽的去除方面：霍夫变换只能去除部分裂缝干扰，且会去除裂缝中的横向片段；快速傅里叶变换在弱化刻槽边缘的同时也弱或了裂缝的边缘信息，使结果图像产生模糊现象。
对此，本发明方法针对水泥刻槽路面裂缝识别过程中由外部因素导致路面对比度过低、由路面材质引起路面斑点以及刻槽自身干扰，提出了一种解决方案。
技术实现要素：
本发明要解决的技术问题在于针对现有技术中的缺陷，提供一种水泥刻槽路面图像降噪增强与裂缝特征提取方法。
本发明解决其技术问题所采用的技术方案是：一种水泥刻槽路面图像降噪增强与裂缝特征提取方法，包括以下步骤：
步骤s100，对原始水泥路面图像进行灰度化处理，然后采用改进的局部自适应对比度增强算法处理灰度化后的图像，以增强图像对比度；
步骤s200，对增强对比度后的图像采用平移不变shearlet变换，针对其中的高频子带应用改良后的p-m扩散模型进行处理，对低频子带采用双边滤波进行处理；
步骤s300，应用基于单向全变分utv模型建立的水泥刻槽路面图像平滑模型，对步骤s200得到的图像进行处理，获得去除刻槽后的路面图像；
步骤s400，采用形态学滤波对去除刻槽后的路面二值化图像进行后处理；
步骤s500，结合连通域标记法、投影法和矩形框法，完成裂缝长度、宽度、倾角、面积的提取。
进一步的，步骤s100所述采用改进的局部自适应对比度增强算法处理灰度化后的图像的实现方式如下，
设图像大小为m×n，滑动窗口的大小为m×n，起始坐标为(x0,y0)，(i,j)为相应像素点位置，f(i,j)是灰度图像(i,j)处的灰度值，favg为全局均值，wavg为滑动窗口内的局部均值，δf为全局标准差，δw为滑动窗口内的局部标准差，则有：
若滑动窗口内局部均值低于全局均值，则增大窗口内像素点的灰度值；若滑动窗口内局部标准差低于全局标准差，则增强图像对比度。其处理过程可表示为：
其中：f”为处理后的图像；(i,j)为相应像素点位置；f(i,j)是灰度图像(i,j)处的灰度值；wavg为滑动窗口内的局部均值；favg为全局均值；δw为滑动窗口内的局部标准差；δf为全局标准差；η,λ,λ1,λ2,λ3均为常数。
进一步的，步骤s200对高频子带应用改良后的p-m扩散模型进行处理的实现方式如下，
记高频子带图像f(x)的sist系数为ci,j，i和j分别代表尺度和方向，高频子带图像的方向局部化过程可表示为：
其中：ψi,j(x)为剪切波函数，其公式如下：
ψi,j(x)＝a-3/4ψ(e-1f-1(x-m))
其中，a∈r+表示尺度参数；e＝(a,0；0,a1/2)为各项异性膨胀矩阵；m∈r2表示平移参数；f＝(1,s；0,1)为剪切波矩阵，s∈r表示剪切波参数；
改进的sist域的p-m去噪模型如下所示：
其中：f0(x)为输入源图像；
其中：λ为常数；η为扩散门限，用于区别噪声和边缘；k为图像梯度阈值。
进一步的，步骤s300所述水泥刻槽路面图像平滑模型的表达式如下，
在utv模型的基础上引入保真项和稀疏优化项，其能量泛函表示如下：
其中：f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
针对上述能量泛函具有不可微和不可分离性，采用admm算法对其进行优化，能量泛函的增广拉格朗日函数可表示如下：
其中：a,b,c为引入的辅助变量，其表达式为
进一步的，a,b,c,g的求解更新方式如下，
因a,b,c,g互不耦合，故可将admm的每次迭代过程转换为下式4个简单子问题的求解，每个子问题的参数以一种交替有序的方式进行更新：
(1)a相关子问题：即l1范数优化问题，使用软阈值收缩算子求解，其迭代公式如下所示：
其中：k表示迭代次数，为正整数；s1为软阈值函数，其表达式如下所示：
其中：t为预先选定的阈值；
(2)b相关子问题：与a相关子问题求解过程相同，其迭代公式为：
其中：k表示迭代次数，为正整数，s1为软阈值函数；
(3)c相关子问题：即l0惩罚项优化问题，使用硬阈值算法求解，其迭代公式如下：
其中：k表示迭代次数，为正整数；s0为硬阈值函数，其表达式如下所示：
其中：t为预先选定的阈值；
(4)g相关子问题：即二次最小化问题，其迭代公式如下，可使用快速傅里叶变换fft求解：
其中：k表示迭代次数，为正整数；
上述(1)～(4)求解过程中的每一步拉格朗日算子可使用下式迭代求解：
其中：k表示迭代次数，为正整数。
进一步的，步骤s400中对去除刻槽后的路面二值化图像先进行开运算，然后对开运算后的结果进行闭运算完成后处理。
进一步的，步骤s500的具体实现方式为，首先采用连通域标记算法统计步骤s400处理后的二值图像中连通域的个数，并据此判断裂缝是网状还是线状；其次，采用投影法判断线性裂缝的方向，对于无法判定方向的线性裂缝采用矩形框法进行二次判定；最后基于矩形框法完成裂缝长度、宽度、倾角、面积信息的提取。
进一步的，基于矩形框法完成裂缝长度、宽度、倾角、面积信息的提取，具体实现方式如下，
(1)线性裂缝长度计算：设利用矩形框法得到的裂缝外接矩形的边长分别为x和y，则裂缝的实际长度可由下式计算：
其中：l表示裂缝的实际长度；α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm；
(2)线性裂缝宽度计算：利用裂缝外接矩形内裂缝的像素总数进行计算，其公式如下：
其中：w表示裂缝的实际宽度；α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm；t1为裂缝外接矩形内裂缝的像素总数；l为裂缝的实际长度；
(3)斜向裂缝倾角计算：根据裂缝外接矩形的两条边进行估算，设斜向裂缝的倾角为θ，则π/6≤θ≤π/3，且
其中：x、y表示裂缝外接矩形的边长；
(4)网状裂缝面积计算：在裂缝外接矩形的基础上采用滑动窗口的方式计算网状裂缝的面积。设窗口的大小为n×n，图像大小为m×n，以滑动窗口的边长为步长对裂缝外接矩形内的子图像进行扫描，并规定：滑动窗口扫描到裂缝称为命中，命中次数用λ表示；对于滑动窗口，若窗口内像素值为0的个数超过一半，则认为该滑动窗口命中裂缝，λ增加1，最后依据命中次数计算网状裂缝的面积s，计算公式如下：
s＝λα2n2
其中：α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm。
进一步的，步骤s100中采用加权均值法取原始水泥图像像素点r,g,b三个颜色分量的均值作为灰度值，其公式如下：
其中：gray表示加权均值法得到的灰度值；r,g,b为原始水泥图像像素点的三种颜色分量；ω1,ω2,ω3为加权系数。
本发明具有以下主要优点：
本发明方法采用改进的局部自适应对比度增强算法，解决了因图像亮度不均、对比度低带来的裂缝识别精度偏低问题；采用改进p-m模型的平移不变shearlet变换去噪算法，最大程度保留了图像边缘信息，提高了图像清晰度，减小了图像失真；基于utv模型提出了一种用于消除刻槽的路面图像平滑算法，不仅将刻槽完整的从图像中去除，而且尽可能清晰地保留了裂缝的边缘信息；采用形态学滤波对去除刻槽后的路面二值化图像进行了后处理，去除了残留的少量孤立噪声点和孤立片段，剔除了裂缝边缘毛刺，修复了断裂的裂缝，填补了缝空洞；提出了一种新的裂缝类型判定方法及裂缝特征计算方法，实现了裂缝特征的数字化描述。此外，本发明方法系统全面、计算量小、易于应用。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明实施例的操作流程图；
图2是传统对比度增强方法与本发明实施例方法处理结果对比图，其中，图2a是本发明实施例采用的第一图像未处理原始路面灰度图像；图2b是采用灰度拉伸法处理第一图像后的图像；图2c是采用直方图均衡法处理第一图像后的图像；图2d是采用经典局部自适应对比度增强算法处理第一图像后的图像；图2e是采用本发明方法改进的局部自适应对比度增强算法处理第一图像后的图像；
图3是传统去噪方法与本发明实施例方法处理结果对比图，其中，图3a是采用本发明方法改进的局部自适应对比度增强算法处理第一图像后的图像；图3b是采用shearlet软阈值法处理图3a后的图像；图3c是采用二维连续小波变换法处理图3a后的图像；图3d是采用传统p-m扩散去噪算法处理图3a后的图像；图3e是采用本发明方法改进p-m模型的sist去噪算法处理图3a后的图像；
图4是传统刻槽去除方法与本发明实施例方法处理结果对比图，其中，图4a是采用霍夫变换刻槽去除法处理图3e后的二值图像；图4b是采用快速傅里叶变换刻槽去除法处理图3e后的二值图像；图4c是采用陷波滤波刻槽去除法处理图3e后的二值图像；图4d是采用本发明方法水泥刻槽路面图像平滑算法处理图3e后的二值图像；
图5是本发明实施例后处理之后的图像，其中，图5a是采用本发明方法水泥刻槽路面图像平滑算法处理图3e后的二值图像；图5b是采用开运算处理图5a后的二值图像；图5c是采用闭运算处理图5a后的二值图像。
图6是利用矩形框法得到的裂缝外接矩形框图。
具体实施方式
下面结合附图和实施例对本发明的技术方案作进一步说明。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
如图1所示，一种水泥刻槽路面图像降噪增强与裂缝特征提取方法，包括以下步骤：
步骤s100、对原始水泥路面图像灰度化处理后，采用改进的局部自适应对比度增强算法，解决因图像亮度不均、对比度低带来的裂缝识别精度偏低问题。
本实施例中首先需要对原始水泥图像灰度化处理，即采用加权均值法取原始水泥图像像素点r,g,b三个颜色分量的均值作为灰度值，其公式如下：
其中：gray表示加权均值法得到的灰度值；r,g,b为原始水泥图像像素点的三种颜色分量；ω1,ω2,ω3为加权系数。本发明取gray＝0.2989r+0.5870g+0.1140b；
其次，沿用经典局部自适应对比度增强算法借助滑动窗口的方式，改变窗口内子图像的处理方式，达到增强图像细节的同时均匀化图像亮度：
设图像大小为m×n，滑动窗口的大小为m×n，起始坐标为(x0,y0)，(i,j)为相应像素点位置，f(i,j)是灰度图像(i,j)处的灰度值，favg为全局均值，wavg为滑动窗口内的局部均值，δf为全局标准差，δw为滑动窗口内的局部标准差，则有：
若滑动窗口内局部均值低于全局均值，则增大窗口内像素点的灰度值；若滑动窗口内局部标准差低于全局标准差，则增强图像对比度。其处理过程可表示为：
其中：f”为处理后的图像；(i,j)为相应像素点位置；wavg为滑动窗口内的局部均值；favg为全局均值；δw为滑动窗口内的局部标准差；δf为全局标准差；η,λ,λ1,λ2,λ3均为常数。λ的取值与经典局部自适应对比度增强算法的取值区间[1,2]相同；η,λ1,λ2,λ3的取值应按实际需要依据均值和标准差的局部与全局对比得到，且λ1,λ2均小于1，η,λ,λ3均大于1；引入λ2δf的目的是为了避免对局部标准差为零的平滑区域进行增强。
步骤s200、对增强对比度后的图像采用平移不变shearlet变换(sist)，并针对其中的高频子带应用改良后的p-m扩散模型进行处理，对低频子带采用双边滤波进行处理，最大程度保留图像边缘信息，提高图像清晰度，减小图像失真：
记高频子带图像f(x)的sist系数为ci,j，i和j分别代表尺度和方向，高频子带图像的方向局部化过程可表示为：
其中：ψi,j(x)为剪切波函数，其公式如下：
ψi,j(x)＝a-3/4ψ(e-1f-1(x-m))
其中，a∈r+表示尺度参数；e＝(a,0；0,a1/2)为各项异性膨胀矩阵；m∈r2表示平移参数；f＝(1,s；0,1)为剪切波矩阵，s∈r表示剪切波参数。
改进的sist域的p-m去噪模型如下所示：
其中：f0(x)为输入源图像；
其中：λ为常数，取0.5；η为扩散门限，用于区别噪声和边缘，取19；k为图像梯度阈值，取58.93。
步骤s300、对去噪后的图像应用基于单向全变分utv模型建立的水泥刻槽路面图像平滑模型去除影响裂缝特征提取的路面刻槽类干扰，尽可能清晰地保留裂缝边缘信息。
本实施例中为了避免经典变分去噪后产生“阶梯效应”导致图像失真，同时为了去除路面图像中残留的孤立噪声点，在utv模型的基础上引入了保真项和稀疏优化项，建立了一种用于消除刻槽的水泥刻槽路面图像平滑模型，其能量泛函表示如下：
其中：f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
针对上述能量泛函具有不可微和不可分离性，在此采用admm算法对其进行优化，能量泛函的增广拉格朗日函数可表示如下：
其中：f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
其中：f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
因a,b,c,g互不耦合，故将admm的每次迭代过程转换为下式4个简单子问题的求解，每个子问题的参数将以一种交替有序的方式进行更新：
其中：f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
每个子问题的求解过程阐述如下：
(1)a相关子问题：即l1范数优化问题，可使用软阈值收缩算子求解，其迭代公式如下所示：
其中：k表示迭代次数，为正整数；g表示经过步骤s200去噪处理后得到的图像；
其中：t为预先选定的阈值。
(2)b相关子问题：与a相关子问题求解过程相同，其迭代公式为：
其中：k表示迭代次数，为正整数；f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
(3)c相关子问题：即l0惩罚项优化问题，可使用硬阈值算法求解，其迭代公式如下：
其中：k表示迭代次数，为正整数；g表示经过步骤s200去噪处理后得到的图像；ω代表图像变换域转换方式，本发明采用曲波变换；λ3为正则化参数，其取值范围为[0.001,0.01]；p3为拉格朗日乘子；ρ3为惩罚系数，用于计算每一步的拉格朗日乘子，其值为ρ3＝100λ1，λ1为正则化参数，其取值范围为[0.1,1]；s0为硬阈值函数，其表达式如下所示：
其中：t为预先选定的阈值。
(4)g相关子问题：即二次最小化问题，其迭代公式如下，可使用快速傅里叶变换fft求解：
其中：k表示迭代次数，为正整数；f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
上述(1)～(4)求解过程中的每一步拉格朗日算子可使用下式迭代求解：
其中：k表示迭代次数，为正整数；f表示带有裂缝的原始水泥路面图像；g表示经过步骤s200去噪处理后得到的图像；
本实施例中路面图像平滑模型的优化步骤是一个迭代处理过程，故以每一步去除的干扰||(f-gk)-(f-gk-1)||与当前残留干扰||(f-gk)||的比值
步骤s400、采用形态学滤波对去除刻槽后的路面二值化图像进行后处理，去除残留的少量孤立噪声点和孤立片段，剔除裂缝边缘毛刺，修复断裂的裂缝，填补缝空洞。
本实施例中首先对路面二值化图像进行图像腐蚀处理，处理过程表示如下：
其中：f代表目标图像；θ表示与操作，w代表结构元素图像，本发明采用矩形化的结构元素；x代表w的平移距离；
然后，对处理后的二值图像进行图像膨胀处理，处理过程表示如下：
其中：f代表目标图像；
图像先经过腐蚀然后膨胀的过程称为开运算，反之则称为闭运算。为了完成二值图像的后处理，本实施例先进行开运算，然后对开运算后的结果进行闭运算，其处理过程表示如下：
其中：f代表目标图像；w代表结构元素图像，在此采用矩形化的结构元素。
步骤s500、提出一种新的裂缝类型判定方法及裂缝特征计算方法，实现裂缝特征的数字化描述，其具体做法为：首先采用连通域标记算法统计步骤s400处理后的二值图像中连通域的个数，并据此判断裂缝是网状还是线状；其次，采用投影法判断线性裂缝的方向，对于无法判定方向的线性裂缝采用矩形框法进行二次判定；最后基于矩形框法完成裂缝长度、宽度、倾角、面积等数字化信息的提取。(1)线性裂缝长度计算：设利用矩形框法得到的裂缝外接矩形的边长分别为x和y，则裂缝的实际长度可由下式计算：
其中：l表示裂缝的实际长度；α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm。
(2)线性裂缝宽度计算：利用裂缝外接矩形内裂缝的像素总数进行计算，其公式如下：
其中：w表示裂缝的实际宽度；α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm；t1为裂缝外接矩形内裂缝的像素总数；l为裂缝的实际长度。
(3)斜向裂缝倾角计算：根据裂缝外接矩形的两条边进行估算，设斜向裂缝的倾角为θ，则π/6≤θ≤π/3，且
其中：x、y表示裂缝外接矩形的边长。
(4)网状裂缝面积计算：在裂缝外接矩形的基础上采用滑动窗口的方式计算网状裂缝的面积。设窗口的大小为n×n，图像大小为m×n，以滑动窗口的边长为步长对裂缝外接矩形内的子图像进行扫描，并规定：滑动窗口扫描到裂缝称为命中，命中次数用λ表示。对于滑动窗口，若窗口内像素值为0的个数超过一半，则认为该滑动窗口命中裂缝，λ增加1，最后依据命中次数计算网状裂缝的面积s，计算公式如下：
s＝λα2n2
其中：α为裂缝图像中单个像素的长度所对应的路面裂缝的实际长度，单位为mm。
对于滑动窗口的边长n需要满足如下条件：
1)n＜2w，其中w表示裂缝的实际宽度；
2)n不能被m和n整除；
3)n的取值不应当过小。
上述过程即为一种水泥刻槽路面图像降噪增强与裂缝特征提取方法的整体技术方案，随后进行实验测试。
步骤s100的实验结果如下：图2a是未处理的原始路面灰度图像。图2b是灰度拉伸后的图像。图2c是直方图均衡化后的图像。图2d是经典局部自适应对比度增强处理后的图像。图2e是采用本发明方法改进的局部自适应对比度增强算法处理后的图像。从图2可以看出，灰度拉伸法虽然对比增强效果较为明显，但同时也加剧了图像亮度不均的程度；直方图均衡法与灰度拉伸法相比，解决了亮度不均的问题，但对比度较低；局部自适应对比度增强算法处理后的图像对比度得到了明显增强，但仍存在图像局部亮度偏暗的问题；相较而言，本发明改进的算法不仅增强了原始灰度图像的对比度，同时也解决了图像亮度不均的问题。
步骤s200的实验结果如下：图3a为经过步骤s100处理后得到的路面灰度图像。图3b是shearlet软阈值去噪结果图像。图3c是二维连续小波变换去噪结果图像。图3d是p-m扩散去噪结果图像。图3e是采用本发明方法改进p-m模型的sist去噪结果图像。从图3可以看出，shearlet阈值法虽然能够去除部分路面斑点干扰，但是在去噪的同时弱化了图像的细节信息，导致图像模糊，而且处理后的图像仍残留部分噪声；二维连续小波去噪的效果要优于shearlet阈值法，但仍然会弱化图像的边缘信息；与前面两种方法相比，p-m扩散去噪比较清晰的保留了裂缝的边缘信息，并去除了部分随机噪声，但是由于扩散过程未考虑图像边缘部分的处理，因此处理后的图像带有黑边，而且去噪效果不理想；相较而言，本发明改进的算法不仅去除了绝大部分路面斑点造成的随机噪声干扰，而且较为完整的保留了裂缝的边缘信息。
步骤s300的实验结果如下：为了更直观的从视觉上看出刻槽去除效果的差异，特将刻槽去除后的路面图像采用统一阈值转换为二值图像。图4a是霍夫变换去除刻槽后的二值图像。图4b是快速傅里叶变换去除刻槽后的二值图像。图4c是陷波滤波去除刻槽后的二值图像。图4d是采用本发明方法水泥刻槽路面图像平滑算法去除刻槽后的二值图像。从图4可以看出，霍夫变换在去除直线的同时会造成部分裂缝断裂或缺失，尤其对包含横向裂缝的路面图像处理效果较差；fft对刻槽的去除效果最不明显，处理后的图像仍然保留有大量的刻槽干扰；陷波滤波从频域角度对刻槽这种周期性干扰进行去除，虽然取得了一定的效果，但仍不能完全除去图像中的刻槽干扰。而且，上述三种方法均未考虑残留噪声的因素，因此二值化图像中会有较多的孤立片段。相比之下，本发明水泥刻槽路面图像平滑算法不仅去除了部分孤立噪声点的干扰，同时也较为清晰地保留了裂缝的边缘，但是二值化的路面图像中仍然存在少量的孤立噪声点和孤立片段，并且裂缝边缘存在毛刺，对此有必要进行相应的后处理。
步骤s400的实验结果如下：图5a是采用本发明方法水泥刻槽路面图像平滑算法去除刻槽后的二值图像。图5b是开运算后的二值图像。图5c是闭运算后的二值图像。从图5可以看出，开运算后的路面二值图像去除了大部分的孤立噪声点和孤立片段，并去除了裂缝边缘的毛刺，但造成了裂缝的断裂，经过闭运算膨胀处理后，基本修复了裂缝中断裂的部分，而且对裂缝骨架中的空洞进行了填充，处理后的裂缝边缘比原始二值图像中的裂缝边缘更加平滑。
步骤s500的实验结果如下：图6是利用矩形框法得到的裂缝外接矩形框图，通过该图可完成裂缝长度、宽度、倾角等数字化信息的提取。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。