CN110363783B - 基于Canny算子的岩体结构面迹线半自动检测方法 - Google Patents
基于Canny算子的岩体结构面迹线半自动检测方法 Download PDFInfo
- Publication number
- CN110363783B CN110363783B CN201910519868.5A CN201910519868A CN110363783B CN 110363783 B CN110363783 B CN 110363783B CN 201910519868 A CN201910519868 A CN 201910519868A CN 110363783 B CN110363783 B CN 110363783B
- Authority
- CN
- China
- Prior art keywords
- fracture
- point
- detection
- points
- segment
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
Abstract
本发明公开了一种基于Canny算子的岩体结构面迹线半自动检测方法,包括如下步骤:控制点输入;图像分割;图像获取;预处理;边缘检测;裂隙处理;迭代连接、剔除伪边缘;边缘叠加,得到完整的裂隙检测结果图。本发明将主要裂隙与次要裂隙分开处理,能够更好地减少二者之间的干扰,而且通过人机交互,能够更加准确的确定迭代连接的角度、距离阈值,更好地保证主次裂隙识别的完整性与正确性,因而具有更高的精确性。
Description
技术领域
本发明属于图像处理领域,具体涉及一种基于Canny算子的岩体结构面迹线半自动检测方法。
背景技术
结构面是岩体中由力学强度较低的部位或岩性相对软弱的夹层构成的不连续面,岩体的变形与稳定性主要取决于结构面的发育条件。所以研究岩体结构面 (rockdiscontinuity structural plane),对研究岩体的力学性质十分重要,具有十分重要的工程意义。
虽然,在前人研究中关于岩体结构面迹线检测较多,但针对岩体结构面复杂的情况全自动检测效果并不理想。目前,工程中常用的测线法或窗口统计法来测量,即通过皮尺和罗盘人工现场逐一测量结构面几何信息(迹长、倾角、间距等)。该法数据处理繁琐、工作量大,并且很多地方无法进行测量。因此,随着计算机技术的发展,运用在该领域的近景摄影测量和数字图像处理技术就应运而生。经典的灰度图像的边缘检测算法有SUSAN边缘检测算子,Canny边缘检测算子,沈俊边缘检测算子等。
但是以上传统的边缘检测算法较为简单,只适用于一些简单的图像的边缘检测,对于复杂的岩体结构面情况检测效果并不理想。所以需要对该检测方法进行改进,使其能够更好地适用于复杂图像的边缘检测。
发明内容
本发明的目的在于提供一种基于Canny算子的岩体结构面迹线半自动检测方法。
本发明解决其技术问题所采用的技术方案是:一种基于Canny算子的岩体结构面迹线半自动检测方法,具体步骤如下:
步骤1、控制点输入:对原照片进行观察,指定主要裂隙的控制点,裂隙控制点包括裂隙起始点、终点和拐点;
步骤2、图像分割:通过指定的区域将原图一分为二,一张为步骤1中的人工指定区域,另一张包括剩余裂隙,其中人工指定区域是通过控制点和偏移量确定的矩形区域;
步骤3、图像获取:将输入的两张彩色图像采用0.3倍的红基色、0.59倍的绿基色和0.11倍的蓝基色进行相加得到两张灰度图像;
步骤4、预处理:对步骤3中获取的两张灰度图像分别通过滤波器进行图像低通滤波,对低通滤波后的图像通过增强算法进行图像增强果;
步骤5、边缘检测:利用Canny算子,对步骤4处理后的两张图片分别进行边缘检测;
步骤6、裂隙处理:对步骤5处理后的两张检测结果进行细化处理,然后删除节点和单点;
步骤7、迭代连接、剔除伪边缘:对步骤6处理后的两张边缘分检测图片通过人机交互输入角度、距离阈值,将满足阈值要求的分段线进行迭代连接,同时剔除伪边缘;
步骤8、边缘叠加:对步骤7处理后的两张边缘检测图片中的迹线检测结果进行叠加,得到完整的裂隙检测结果图。
与现有技术相比,本发明的显著优点为:(1)本发明将主要裂隙与次要裂隙分开处理,能够更好地减少二者之间的干扰,而且通过人机交互,能够更加准确的确定迭代连接的角度、距离阈值,所以能够更好地保证主次裂隙识别的完整性与正确性,因而具有更高的精确性。
附图说明
图1为基于Canny算子的岩体结构面迹线半自动检测方法流程图。
图2为人工指定区域示意图。
图3为线段节点、内部点、端点和单点定义图。
图4为求两分段线最近端点示意图。
图5为最近两端点分别到另一分段线的距离示意图。
图6为最近两端点连线与另一分段线的夹角示意图。
图7为两分段线水平夹角之差示意图。
图8为两分段线相似性分析流程图。
图9为岩体露头面图像示意图。
图10为人工指定少量难以正确识别的主要裂隙图像示意图。
图11为人工指定裂隙区域提取图像示意图。
图12为提取人工指定区域后剩余的结构面图像示意图。
图13为人工指定区域主要裂隙检测结果图。
图14为Canny边缘检测、细化、删除节点和单点结果图。
图15为剩余区域裂隙检测结果图。
图16为人工指定区域与剩余区域的岩体结构面迹线叠加结果图。
具体实施方式
如图1所示,一种基于Canny算子的岩体结构面迹线半自动检测方法,将主次裂隙区域分离,实施边缘像素检测的后处理以消除伪边缘。具体步骤如下:
步骤1、控制点输入:对原照片进行观察,人工指定主要裂隙的控制点,裂隙控制点包括裂隙起始点、终点和拐点;主要裂隙是指裂隙大小超过设定阈值的裂隙。这一步骤主要是通过人工输入控制点,然后将难以正确识别的主要裂隙存在的区域确定出来,为提取该区域做准备。
步骤2、图像分割:通过人工指定的区域将原图一分为二,一张包括难以正确识别的主要裂隙,另一张包括剩余裂隙。其中人工指定区域是通过控制点和偏移量确定的矩形区域,偏移量为输入图像纵向像素值h的1/20。如图2中控制点 A、B,人工指定区域则为线段AB分别向垂直线段AB的上、下方向偏移h/20 个像素点;这一步能够将因为相互交错而导致难以正确识别的主次裂隙区域分开,从而减少二者的干扰,提高检测的精确性。
步骤3、图像获取:将输入的两张彩色图像采用0.3倍的红基色、0.59倍的绿基色和0.11倍的蓝基色进行相加得到两张灰度图像;
步骤4、预处理:对c中获取的两张灰度图像分别通过滤波器进行图像低通滤波,达到图像去噪的效果;然后对低通滤波后的图像通过增强算法,达到图像增强的效果;
步骤5、边缘检测:利用Canny算子,对d处理后的两张图片分别进行边缘检测;
步骤6、裂隙处理:对e处理后的两张检测结果进行细化处理,本发明采用形态学细化算法将裂隙细化成单像素宽度的裂隙,然后删除节点和单点。其中节点和单点的定义为:对e处理后的图像(白底黑线)每个像素点从左往右、从上往下依次进行分析,如果当前检测像素点P(i,j)为黑色时,统计其四周八个位置中为黑色像素点的个数并记为N。当N≥3时,则该点应为节点;当N=2时,则该点为分段线内部点;当N=1时,则该点为该分段线的端点;当N=0时,则该点为单点。如图3中的 (a)中N=3,则像素点P(i,j)为节点,图3中的 (b)中N=2,则像素点P(i,j)为线段内部点,图3中的 (c)中N=1,则像素点P(i,j)为线段端点,图3中的 (d)中N=0,则像素点P(i,j)为单点;
因为对图像进行边缘检测的结果宽度多为多个像素,将其细化为单像素宽度有助于减少图形冗余信息量,突出图形特征,这样可减少运算量从而缩短识别的时间和提高识别率。还有图像边缘检测及细化的结果往往存在分支(毛刺、迹线交叉)的情况,这有碍于对分段线进行拟合和参数提取。为了得到更加准确的迹线参数,就需要清除细化后图像中的节点,同时将检测的单点视为噪声进而清除。
步骤7、迭代连接、剔除伪边缘:对f处理后的两张边缘分检测图片通过人机交互输入角度、距离阈值(Dt、dt、γt、βt),将满足阈值要求的分段线进行迭代连接,同时剔除伪边缘。如图8所示,角度、距离阈值的定义如下:
①如图4,L1和L2分别为两条线段,a11、a12分别为线段L1的两端点, a21、a22分别为线段L2的两端点,d1、d2表示a11与线段L2两端点的距离, d3、d4表示a12与线段L2两端点的距离,dmin表示d1、d2、d3、d4的最小值,即图4中a12、a21之间的距离。如果dmin小于给定的判定阈值Dt,则继续往下判定。
②如图5所示,由上知a12、a21表示两条分段线距离最近的两个端点,D1、 D2分别表示a12到线段L2和a21到线段L1的距离。Sd=D1+D2,如果Sd小于给定的判定阈值dt,则继续往下判定。
③A、角度准则1:如图6所示,a12、a21连线与线段L1的夹角记为γ,如果γ小于给定的判定阈值γt,则继续往下判定。
B、角度准则2:如图7所示,α1,α2分别为L1和L2的倾斜角,β为两倾角差,β≤90°,当|α1-α2|≤90°时,β=|α1-α2|;当|α1-α2|>90°时,β=180-|α1-α2|。如果β小于给定的判定阈值βt,则继续往下判定。
④如果L1和L2两线段满足上述条件①②③,则L1和L2有很大可能属于同一条结构面的迹线。L1和L2的相似系数Sij求解如下,相似系数越大则二者属于同一结构面裂隙的可能性越大:
其中w1、w2、w3、w4、θ1、θ2均为大于0的常数。
这一步充分体现了人机交互的优越性,通过人工输入角度、距离阈值,能够更加准确地对属于同一结构面迹线的分段线进行迭代连接,伪边缘因为无法满足阈值要求无法连接变长,进而滤除。
步骤8、边缘叠加:对g处理后的两张边缘检测图片中的迹线检测结果进行叠加,得到完整的裂隙检测结果图,然后输出。
下面参照附图对本发明的示例性实施方式进行详细描述。
实施例
原图为图9,对岩体露头面中容易识别错误的主要裂隙进行人工指定控制点,如图10所示,为了便于观察,此处将属于同一裂隙的控制点进行了连接。然后通过算法将原图一分为二,分别是提取区域图11和剩余区域图12。之后分别对二图进行Canny边缘检测、细化、删除节点和单点以及迭代连接,如图14所示,二图的处理结果如图13、图15所示。最后将二者的处理结果叠加,即为最终结果图16。本发明将主次裂隙区域分离,实施边缘像素检测的后处理以消除伪边缘,能够显著提高岩体结构面迹线检测的完整性与正确性,因而具有更高的精确性。
Claims (4)
1.一种基于Canny算子的岩体结构面迹线半自动检测方法,其特征在于,包括如下步骤:
步骤1、控制点输入:对原照片进行观察,指定主要裂隙的控制点,裂隙控制点包括裂隙起始点、终点和拐点;
步骤2、图像分割:通过指定的区域将原图一分为二,一张为步骤1中的人工指定区域,另一张包括剩余裂隙,其中人工指定区域是通过控制点和偏移量确定的矩形区域;
步骤3、图像获取:将输入的两张彩色图像采用0.3倍的红基色、0.59倍的绿基色和0.11倍的蓝基色进行相加得到两张灰度图像;
步骤4、预处理:对步骤3中获取的两张灰度图像分别通过滤波器进行图像低通滤波,对低通滤波后的图像通过增强算法进行图像增强果;
步骤5、边缘检测:利用Canny算子,对步骤4处理后的两张图片分别进行边缘检测;
步骤6、裂隙处理:对步骤5处理后的两张检测结果进行细化处理,然后删除节点和单点;
步骤7、迭代连接、剔除伪边缘:对步骤6处理后的两张边缘分检测图片通过人机交互输入角度、距离阈值,将满足阈值要求的分段线进行迭代连接,同时剔除伪边缘,具体为:
(1)L1和L2分别为两条线段,a11、a12分别为线段L1的两端点,a21、a22分别为线段L2的两端点,d1、d2表示a11与线段L2两端点的距离,d3、d4表示a12与线段L2两端点的距离,dmin表示d1、d2、d3、d4的最小值,如果dmin小于给定的判定阈值Dt,则继续往下判定;
(2)设a12、a21表示两条分段线距离最近的两个端点,D1、D2分别表示a12到线段L2和a21到线段L1的距离,Sd=D1+D2,如果Sd小于给定的判定阈值dt,则继续往下判定;
(3)角度准则1:a12、a21连线与线段L1的夹角记为γ,如果γ小于给定的判定阈值γt,则继续往下判定;
角度准则2:α1,α2分别为L1和L2的倾斜角,β为两倾角差,当|α1-α2|≤90°时,β=|α1-α2|;当|α1-α2|>90°时,β=180-|α1-α2|;如果β小于给定的判定阈值βt,则继续往下判定;
(4)L1和L2的相似系数Sij求解如下:
其中w1、w2、w3、w4、θ1、θ2均为大于0的常数;
步骤8、边缘叠加:对步骤7处理后的两张边缘检测图片中的迹线检测结果进行叠加,得到完整的裂隙检测结果图。
2.根据权利要求1所述的基于Canny算子的岩体结构面迹线半自动检测方法,其特征在于,步骤2中,所述偏移量为输入图像纵向像素值h的1/20。
3.根据权利要求1所述的基于Canny算子的岩体结构面迹线半自动检测方法,其特征在于,步骤6中细化处理是用形态学细化算法将裂隙细化成单像素宽度的裂隙。
4.根据权利要求1所述的基于Canny算子的岩体结构面迹线半自动检测方法,其特征在于,步骤6中节点和单点的定义为:对步骤5处理后的图像每个像素点从左往右、从上往下依次进行分析,如果当前检测像素点为黑色时,统计其四周八个位置中为黑色像素点的个数并记为N;当N≥3时,则该点应为节点;当N=2时,则该点为分段线内部点;当N=1时,则该点为该分段线的端点;当N=0时,则该点为单点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910519868.5A CN110363783B (zh) | 2019-06-17 | 2019-06-17 | 基于Canny算子的岩体结构面迹线半自动检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910519868.5A CN110363783B (zh) | 2019-06-17 | 2019-06-17 | 基于Canny算子的岩体结构面迹线半自动检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110363783A CN110363783A (zh) | 2019-10-22 |
CN110363783B true CN110363783B (zh) | 2022-08-16 |
Family
ID=68217236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910519868.5A Active CN110363783B (zh) | 2019-06-17 | 2019-06-17 | 基于Canny算子的岩体结构面迹线半自动检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110363783B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115100224B (zh) * | 2022-06-29 | 2024-04-23 | 中国矿业大学 | 一种煤矿巷道掘进迎头交叉裂隙提取方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102620673A (zh) * | 2012-03-16 | 2012-08-01 | 同济大学 | 基于图像分析的隧道变形在线监测系统及其应用 |
CN109033538B (zh) * | 2018-06-30 | 2022-07-22 | 南京理工大学 | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 |
-
2019
- 2019-06-17 CN CN201910519868.5A patent/CN110363783B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110363783A (zh) | 2019-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108921176B (zh) | 一种基于机器视觉的指针式仪表定位与识别方法 | |
CN104751142B (zh) | 一种基于笔划特征的自然场景文本检测方法 | |
CN111260616A (zh) | 一种基于Canny算子二维阈值分割优化的绝缘子裂纹检测方法 | |
WO2016091016A1 (zh) | 一种基于胞核标记分水岭变换的粘连白细胞分割方法 | |
CN102999886B (zh) | 图像边缘检测器及标尺光栅栅线精度检测系统 | |
CN106599890B (zh) | 变电站巡检机器人数字式仪表识别算法 | |
CN111105389B (zh) | 融合Gabor滤波器与卷积神经网络的路面裂缝的检测方法 | |
CN110838117A (zh) | 一种基于孔壁图像的岩石面孔隙度识别方法 | |
CN111667477B (zh) | 磁材尺寸缺陷检测方法、装置、检测设备及可读存储介质 | |
CN104700420A (zh) | 基于Hough变换的椭圆检测方法、系统和虫卵识别方法 | |
CN107133623A (zh) | 一种基于背景差和圆心定位的指针位置精确检测方法 | |
CN111652844B (zh) | 一种基于数字图像区域增长的x射线缺陷检测方法及系统 | |
CN112307803A (zh) | 数字地质露头裂缝提取方法及装置 | |
CN110363783B (zh) | 基于Canny算子的岩体结构面迹线半自动检测方法 | |
CN113129323A (zh) | 一种基于人工智能的遥感田埂界线检测方法、系统、计算机设备及存储介质 | |
CN109187548A (zh) | 一种岩体裂隙识别方法 | |
CN104966283A (zh) | 图像分层配准方法 | |
CN106324708B (zh) | 雨量记录图纸的数字化方法、装置 | |
CN102313740B (zh) | 太阳能电池板裂痕检测方法 | |
CN116452613B (zh) | 一种地质调查中裂缝轮廓提取方法 | |
CN110298816B (zh) | 一种基于图像重生成的桥梁裂缝检测方法 | |
CN107993193B (zh) | 基于光照均衡化和改进surf算法的隧道衬砌图像拼接方法 | |
CN105160300B (zh) | 一种基于水平集分割的文本抽取方法 | |
CN115719358A (zh) | 一种x射线安检图像中直线段提取方法 | |
CN115170507A (zh) | 基于图像数据的注浆管表面缺陷检测方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |