一种基于图像局部增强的山体裂缝检测方法技术

本发明专利技术公开了一种山体裂缝检测方法，包括通过图像滑动窗口对采集的图像进行分析，确定图像的每个像素点的对比度；根据该对比度，构建对比度图，对对比度图进行分析，确定对比度图的最佳分割阈值，根据该阈值，确定出待增强的区域，并对待增强区域进行图像增强处理；对处理后的图像进行分析，确定出灰度特征图和方差特征图；求取灰度特征图和反差特征图的最佳分割阈值，并求取灰度特征和方差特征的二值图，对两个二值图进行求与操作，确定出裂缝的二值图；对目标裂缝的二值图进行变换确定其细化骨架，并对骨架进行区域延伸处理，得到骨架上所有裂缝潜在区域；对骨架的像素进行遍历，得到裂缝的长度信息，并根据区域的像素个数，计算裂缝的面积。

【技术实现步骤摘要】
一种基于图像局部增强的山体裂缝检测方法
本专利技术涉及数字图像处理
，具体来说，涉及一种基于图像局部增强的山体裂缝检测方法。
技术介绍
随着自然环境的不断恶化和人为工程建设的越来越多，山体滑坡数目每年都在增多，且强度和频率每年都在增大，中国每年由山体滑坡带来的损失超过4亿美元。山体滑坡带来巨大的财产损失、人员伤亡，早期预警系统的建立非常必要。山体滑坡预警系统目前都采用监测仪器进行监测，传统的用于山体滑坡监测的方法大概有如下几种：常规大地测量法；液体静力水准测量、重力测量法；地下水位监测法；电测法、地下钻孔倾斜等。尽管这些方法都在山体滑坡预警方面起到了积极作用，但它们具有一些弊端，如受地形影响，不能连续观测，自动化程度不高，人力投入过大，数据不能实时处理等。采用基于数字图像处理技术的监测方法则可以消除上述部分弊端，如无接触式测量，智能化程度高，成本低，可以实时预警等。本专利技术基于山体滑坡前部分区域会出现裂缝的征兆，采用视频采集设备对山体进行全天候监控，采用图像处理技术对裂缝进行识别。但山体裂缝一般与背景对比度较低，造成了裂缝识别的困难。针对相关技术中的问题，目前尚未提出有效的解决方案。
技术实现思路
针对相关技术中的问题，本专利技术提出一种基于图像局部增强的山体裂缝检测方法，其解决了山体裂缝识别困难的问题，且较传统仪器监测方法具有无接触测量的优势，适用于野外实时监控环境。本专利技术的技术方案是这样实现的：一种基于图像局部增强的山体裂缝检测方法，包括：通过预先设定的图像滑动窗口，对采集的图像进行分析，确定所述滑动窗口中心点处的灰度均值，并根据所述灰度均值，确定所述图像的每个像素点的对比度；根据每个像素点的对比度，构建与所述图像对应的对比度图，并利用ostu自动分割阈值算法，对对比度图进行分析，确定对比度图的最佳分割阈值，并根据该最佳分割阈值，确定出待增强的区域，并对确定出的待增强区域进行图像增强处理；对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图；并利用上述ostu自动分割阈值算法，求取所述灰度特征图和方差特征图的最佳分割阈值，并利用该分割阈值，求取灰度特征和方差特征的二值图，并对两个二值图进行求与操作，确定出目标裂缝的二值图；采用击中击不中变换对所述目标裂缝的二值图进行变换，确定其细化骨架，并对所述细化骨架进行区域延伸处理，得到所述细化骨架上所有裂缝潜在区域；对上述区域延伸细化骨架的像素进行遍历，得到裂缝的长度信息，并且根据遍历出的区域的像素个数，计算裂缝的面积信息。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：在利用ostu自动分割阈值算法，对对比度图进行分析之前，对对比度图进行滤波处理，以滤除掉对比度较高的噪声点。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：在对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图之前，对进行图像增强处理后的图像依次进行中值滤波处理和匀光处理。优选的，所述匀光处理采用直方图均衡化技术。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：根据裂缝的面积信息，发出对应的预警信息。本专利技术的有益效果：本专利技术引入了数字图像处理技术，图像处理算法自适应提高裂缝处的对比度，解决了山体裂缝识别困难的问题，且较传统仪器监测方法具有无接触测量的优势，适用于野外实时监控环境。具体实施方式下面将结合本专利技术实施例，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本专利技术保护的范围。根据本专利技术的实施例，提供了一种基于图像局部增强的山体裂缝检测方法，包括：通过预先设定的图像滑动窗口，对采集的图像进行分析，确定所述滑动窗口中心点处的灰度均值，并根据所述灰度均值，确定所述图像的每个像素点的对比度；根据每个像素点的对比度，构建与所述图像对应的对比度图，并利用ostu自动分割阈值算法，对对比度图进行分析，确定对比度图的最佳分割阈值，并根据该最佳分割阈值，确定出待增强的区域，并对确定出的待增强区域进行图像增强处理；对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图；并利用上述ostu自动分割阈值算法，求取所述灰度特征图和方差特征图的最佳分割阈值，并利用该分割阈值，求取灰度特征和方差特征的二值图，并对两个二值图进行求与操作，确定出目标裂缝的二值图；采用击中击不中变换对所述目标裂缝的二值图进行变换，确定其细化骨架，并对所述细化骨架进行区域延伸处理，得到所述细化骨架上所有裂缝潜在区域；对上述区域延伸细化骨架的像素进行遍历，得到裂缝的长度信息，并且根据遍历出的区域的像素个数，计算裂缝的面积信息。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：在利用ostu自动分割阈值算法，对对比度图进行分析之前，对对比度图进行滤波处理，以滤除掉对比度较高的噪声点。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：在对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图之前，对进行图像增强处理后的图像依次进行中值滤波处理和匀光处理。优选的，所述匀光处理采用直方图均衡化技术。进一步的，上述基于图像局部增强的山体裂缝检测方法，还包括：根据裂缝的面积信息，发出对应的预警信息。为了方便理解本专利技术的上述技术方案，以下通过具体使用步骤和使用原理对本专利技术的上述技术方案进行详细说明。第一步图像局部增强输入图像设为f(i，j)，设定滑动窗口，窗口大小为5×5，中心点处的灰度均值计算为：式中，GM表示图像行值为i、列值为j处的灰度均值，k和l分别表示窗口内图像的行值和列值，f(k，l)为以(i，j)为中心的窗口内像素点(k，l)的灰度值，∑为求和操作符号。中心点处的对比度EM计算为：式中，p(f(k，l))表示灰度值f(k，l)在5×5窗口邻域中出现的概率估计值，计算公式如下：式中，nl，k为灰度值f(k，l)在5×5窗口邻域内出现的次数，通过遍历窗口统计即可获得。滑动窗口对整幅图像的每个像素点都进行对比度计算，由此可以得到原图对应的对比度图E(i，j)；采用3×3的模板对E(i，j)进行中值滤波，滤除掉对比度较高的噪声点。利用ostu自动分割阈值算法求取对比度图E(i，j)的最佳分割阈值d，通过分割阈值d可以求取出待增强的区域，增强算法采用如下公式：y(k，l)＝λ(x(k，l)-G)+x(k，l)式中，y(k，l)表示窗口内增强后的灰度值，x(k，l)为增强前的灰度值，λ表示动态增强系数，根据实际情况进行调整，G表示窗口内所有像素的灰度均值。对所有待增强区域都进行如上式计算。至此，完成了落石图像的局部增强。第二步图像预处理对预处理的图像首先采用5×5模板中值滤波，将滤波后的图像fL(i，j)进行匀光处理，目的是针对图像光照不均的情况，用来均衡图像的光照强度。匀光处理采用直方图均衡化技术。直方图均衡化后的图像为fZ(i，j)。第三步图像分割本方法基于裂缝目标的特征进行分割，所使用的目标特征为裂缝图像的灰度和方差。设定滑动窗口大小为3×3，窗内图像的均值GN为：式中，fZ本文档来自技高网...

【技术保护点】
一种基于图像局部增强的山体裂缝检测方法，其特征在于，包括：通过预先设定的图像滑动窗口，对采集的图像进行分析，确定所述图像的中心点处的灰度均值，并根据所述灰度均值，确定所述图像的每个像素点的对比度；根据每个像素点的对比度，构建与所述图像对应的对比度图，并利用ostu自动分割阈值算法，对对比度图进行分析，确定对比度图的最佳分割阈值，并根据该最佳分割阈值，确定出待增强的区域，并对确定出的待增强区域进行图像增强处理；对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图；并利用上述ostu自动分割阈值算法，求取所述灰度特征图和反差特征图的最佳分割阈值，并利用该分割阈值，求取灰度特征和方差特征的二值图，并对两个二值图进行求与操作，确定出目标裂缝的二值图；采用击中击不中变换对所述目标裂缝的二值图进行变换，确定其细化骨架，并对所述细化骨架进行区域延伸处理，得到所述细化骨架上所有裂缝潜在区域；对上述细化骨架的像素进行遍历，得到裂缝的长度信息，并且根据遍历出的区域的像素个数，计算裂缝的面积信息。

【技术特征摘要】
1.一种基于图像局部增强的山体裂缝检测方法，其特征在于，包括：通过预先设定的图像滑动窗口，对采集的图像进行分析，确定所述滑动窗口的中心点处的灰度均值，并根据所述灰度均值，确定所述图像的每个像素点的对比度；根据每个像素点的对比度，构建与所述图像对应的对比度图，并利用ostu自动分割阈值算法，对对比度图进行分析，确定对比度图的最佳分割阈值，并根据该最佳分割阈值，确定出待增强的区域，并对确定出的待增强区域进行图像增强处理；对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图；并利用上述ostu自动分割阈值算法，求取所述灰度特征图和方差特征图的最佳分割阈值，并利用该分割阈值，求取灰度特征和方差特征的二值图，并对两个二值图进行求与操作，确定出目标裂缝的二值图；采用击中击不中变换对所述目标裂缝的二值图进行变换，确定其细化骨架，并对所述细化骨架进行区域延伸处理，得到所述细化骨架...

【专利技术属性】
技术研发人员：焦书森，惠飞菲，宋杨，
申请(专利权)人：北京宇航时代科技发展有限公司，
类型：发明
国别省市：北京;11