CN104992429B - 一种基于图像局部增强的山体裂缝检测方法 - Google Patents

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

Info

Publication number
CN104992429B
CN104992429B CN201510195310.8A CN201510195310A CN104992429B CN 104992429 B CN104992429 B CN 104992429B CN 201510195310 A CN201510195310 A CN 201510195310A CN 104992429 B CN104992429 B CN 104992429B
Authority
CN
China
Prior art keywords
image
contrast
crack
region
feature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510195310.8A
Other languages
English (en)
Other versions
CN104992429A (zh
Inventor
焦书森
惠飞菲
宋杨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Astronavigation Age Development In Science And Technology Co Ltd
Original Assignee
Beijing Astronavigation Age Development In Science And Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Astronavigation Age Development In Science And Technology Co Ltd filed Critical Beijing Astronavigation Age Development In Science And Technology Co Ltd
Priority to CN201510195310.8A priority Critical patent/CN104992429B/zh
Publication of CN104992429A publication Critical patent/CN104992429A/zh
Application granted granted Critical
Publication of CN104992429B publication Critical patent/CN104992429B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

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

Description

一种基于图像局部增强的山体裂缝检测方法
技术领域
本发明涉及数字图像处理技术领域，具体来说，涉及一种基于图像局部增强的山体裂缝检测方法。
背景技术
随着自然环境的不断恶化和人为工程建设的越来越多，山体滑坡数目每年都在增多，且强度和频率每年都在增大，中国每年由山体滑坡带来的损失超过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(k，l)表示窗口内坐标为(k，l)处直方图均衡化后的图像灰度值。
窗内图像的方差为，
滑动窗口对每个像素点都进行上述计算，便可以得到原图对应的灰度特征图g(i，j)和方差特征图v(i，j)；
同样利用ostu自动分割阈值算法求取灰度特征图g(i，j)和方差特征图v(i，j)的最佳分割阈值，利用分割阈值可以求取灰度特征和方差特征的二值图Bg(i，j)和Bv(i，j)，对两个二值图进行求与操作，其结果则为目标裂缝的二值图B(i，j)。
式中，&表示求与操作。
第四步骨架提取及区域延伸
采用击中击不中变换对二值图B(i，j)进行变换，求取其细化骨架U(i，j)。细化的骨架会有很多断裂，利用长度和方向信息，采用如下步骤进行区域延伸
a.遍历U(i，j)的每个骨架区域，求取每个区域的长度值，去除掉小于TL＝3的区域，剩下区域为潜在区域，TL根据需求设定；
b.遍历潜在区域，利用最小二乘线性拟合计算每个潜在区域的方向α，在每个区域的头部和尾部做α±β(β＝2°)角度，r＝5为半径的扇形搜索，r为断裂误差值，β为方向偏差值；
c.将落在扇形区域的所有候选区域，按长度从小到大开始遍历，与当前潜在区域进行匹配，匹配成功则停止搜索，连接两个潜在区域，并将新的区域重新作为潜在区域，重复b，c两步，直到没有用于延伸的区域，此潜在区域不再进行搜索；
d.重复b，c步骤，直到所有潜在区域遍历完毕。
第五步裂缝参数计算
完成上述步骤后，通过遍历每个骨架的像素个数，可以得出裂缝的长度信息，遍历区域的像素个数，则可以得到裂缝的面积信息，由此可以给出相应预警信息。
至此，实现了基于图像局部增强的山体裂缝检测方法。
综上所述，本发明引入了数字图像处理技术，图像处理算法自适应提高裂缝处的对比度，解决了山体裂缝识别困难的问题，且较传统仪器监测方法具有无接触测量的优势，适用于野外实时监控环境。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种基于图像局部增强的山体裂缝检测方法，其特征在于，包括：
通过预先设定的图像滑动窗口，对采集的图像进行分析，确定所述滑动窗口的中心点处的灰度均值，并根据所述灰度均值，确定所述图像的每个像素点的对比度；
根据每个像素点的对比度，构建与所述图像对应的对比度图，并利用ostu自动分割阈值算法，对对比度图进行分析，确定对比度图的最佳分割阈值，并根据该最佳分割阈值，确定出待增强的区域，并对确定出的待增强区域进行图像增强处理；
对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图；
并利用上述ostu自动分割阈值算法，求取所述灰度特征图和方差特征图的最佳分割阈值，并利用该分割阈值，求取灰度特征和方差特征的二值图，并对两个二值图进行求与操作，确定出目标裂缝的二值图；
采用击中击不中变换对所述目标裂缝的二值图进行变换，确定其细化骨架，并对所述细化骨架进行区域延伸处理，得到所述细化骨架上所有裂缝潜在区域；
对上述区域延伸的细化骨架的像素进行遍历，得到裂缝的长度信息，并且根据遍历出的区域的像素个数，计算裂缝的面积信息。
2.根据权利要求1所述的一种基于图像局部增强的山体裂缝检测方法，其特征在于，还包括：
在利用ostu自动分割阈值算法，对对比度图进行分析之前，对对比度图进行滤波处理，以滤除掉对比度较高的噪声点。
3.根据权利要求1所述的一种基于图像局部增强的山体裂缝检测方法，其特征在于，还包括：
在对进行图像增强处理后的图像进行分析，确定出与图像对应的灰度特征图和方差特征图之前，对进行图像增强处理后的图像依次进行中值滤波处理和匀光处理。
4.根据权利要求3所述的一种基于图像局部增强的山体裂缝检测方法，其特征在于，所述匀光处理采用直方图均衡化技术。
5.根据权利要求1-4中任意一项所述的一种基于图像局部增强的山体裂缝检测方法，其特征在于，还包括：
根据裂缝的面积信息，发出对应的预警信息。
CN201510195310.8A 2015-04-23 2015-04-23 一种基于图像局部增强的山体裂缝检测方法 Active CN104992429B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510195310.8A CN104992429B (zh) 2015-04-23 2015-04-23 一种基于图像局部增强的山体裂缝检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510195310.8A CN104992429B (zh) 2015-04-23 2015-04-23 一种基于图像局部增强的山体裂缝检测方法

Publications (2)

Publication Number Publication Date
CN104992429A CN104992429A (zh) 2015-10-21
CN104992429B true CN104992429B (zh) 2017-09-12

Family

ID=54304237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510195310.8A Active CN104992429B (zh) 2015-04-23 2015-04-23 一种基于图像局部增强的山体裂缝检测方法

Country Status (1)

Country Link
CN (1) CN104992429B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105303543A (zh) * 2015-10-23 2016-02-03 努比亚技术有限公司 图像增强方法及移动终端
CN107909567B (zh) * 2017-10-31 2022-02-15 华南理工大学 数字图像的细长型连通区域提取方法
CN109655895B (zh) * 2019-01-29 2019-10-25 中国科学院地质与地球物理研究所 一种任意复杂裂缝介质格子法弹性建模方法
CN110378866B (zh) * 2019-05-22 2021-04-06 中国水利水电科学研究院 一种基于无人机巡检的渠道衬砌破损图像识别方法
CN110517255A (zh) * 2019-08-29 2019-11-29 北京理工大学 基于吸引子模型的浅裂缝检出方法
CN110751623A (zh) * 2019-09-06 2020-02-04 深圳新视智科技术有限公司 基于联合特征的缺陷检测方法、装置、设备及存储介质
CN112446871B (zh) * 2020-12-02 2022-11-15 山东大学 一种基于深度学习和OpenCV的隧道裂缝识别方法
CN114418957B (zh) * 2021-12-24 2022-11-18 广州大学 基于机器人视觉的全局和局部二值模式图像裂缝分割方法
CN116071893B (zh) * 2023-04-06 2023-07-14 湖南智慧平安科技有限公司 基于大数据与计算机视觉的预警中控指挥平台
CN116342635B (zh) * 2023-05-26 2023-08-08 山东省地质矿产勘查开发局第一地质大队（山东省第一地质矿产勘查院） 一种地质测绘中裂缝轮廓提取方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006120891A (ja) * 2004-10-22 2006-05-11 Matsushita Electric Ind Co Ltd 固体撮像素子の評価方法および評価装置
CN103870833A (zh) * 2014-03-31 2014-06-18 武汉工程大学 基于凹性测度的路面裂缝提取与测评方法
CN104021574A (zh) * 2014-07-04 2014-09-03 武汉武大卓越科技有限责任公司 路面病害自动识别方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006120891A (ja) * 2004-10-22 2006-05-11 Matsushita Electric Ind Co Ltd 固体撮像素子の評価方法および評価装置
CN103870833A (zh) * 2014-03-31 2014-06-18 武汉工程大学 基于凹性测度的路面裂缝提取与测评方法
CN104021574A (zh) * 2014-07-04 2014-09-03 武汉武大卓越科技有限责任公司 路面病害自动识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A Crack Segmentation approach using the Combination of Gray Thresholds and Fractal Feature;Yang Song et al;《Advanced Materials Research》;20120315;第487卷;第622-626页 *
结合NSCT和图像形态学的路面裂缝检测;马常霞 等;《计算机辅助设计与图像学学报》;20091230;第21卷(第12期);第1761-1767页 *
融合多特征与格式塔理论的路面裂缝检测;徐威 等;《计算机辅助设计与图形学学报》;20150131;第27卷(第1期);第147-156页 *

Also Published As

Publication number Publication date
CN104992429A (zh) 2015-10-21

Similar Documents

Publication Publication Date Title
CN104992429B (zh) 一种基于图像局部增强的山体裂缝检测方法
US8463065B2 (en) Linear feature detection method and apparatus
CN101826209B (zh) 一种基于Canny模型的三维医学图像分割方法
CN107945179A (zh) 一种基于特征融合的卷积神经网络的肺结节良恶性检测方法
CN105513077A (zh) 一种用于糖尿病性视网膜病变筛查的系统
CN108921076B (zh) 基于图像的道面裂缝病害自适应恒虚警检测方法
CN107392885A (zh) 一种基于视觉对比机制的红外弱小目标检测方法
CN107798679A (zh) 乳腺钼靶图像乳房区域分割与钙化点检测方法
WO2011038607A1 (zh) 人体检查图像处理方法和人体检查设备
CN104616006B (zh) 一种面向监控视频的胡须人脸检测方法
CN102068281A (zh) 一种占位性病变超声图像的处理方法
CN107564035B (zh) 基于重要区域识别和匹配的视频跟踪方法
CN108647693A (zh) 一种基于二元显著性特征的海面红外目标检测方法
CN105225216A (zh) 基于空间距圆标记边缘检测的虹膜预处理算法
Zhao et al. Principal curvature for infrared small target detection
CN106327451A (zh) 一种古生物化石的图像修复方法
CN113850328A (zh) 基于多视角深度学习的非小细胞肺癌亚型分类系统
CN108108740B (zh) 一种主动毫米波人体图像性别识别方法
CN106780718A (zh) 一种古生物化石的三维重建方法
CN116758421A (zh) 一种基于弱监督学习的遥感图像有向目标检测方法
CN104992139B (zh) 一种基于运动矢量与圆形度匹配的落石检测方法
CN114972272A (zh) 一种基于Grad-CAM的新冠肺炎病灶分割方法
CN104899592A (zh) 一种基于圆形模板的道路半自动提取方法及系统
CN103279960A (zh) 一种基于x射线背散射图像的人体隐藏物的图像分割方法
CN112926516A (zh) 一种鲁棒的手指静脉图像感兴趣区域提取方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant