一种复杂裂隙网络的裂缝提取及优化方法技术

本发明专利技术公开了一种复杂裂隙网络的裂缝提取及优化方法，首先采用直方图均衡化、高斯滤波方式对待处理的复杂裂隙网络进行处理，去除该复杂裂隙网络的裂隙图像中的随机性噪声；然后基于伽马变换的裂隙边界增强算法，增强裂隙图像的裂缝信息、矫正裂隙图像的亮度，获得复杂裂隙网络数据；建立基于形态特征的伪目标去除模型，从复杂裂隙网络数据中提取出裂隙网络；将裂隙网络作为先验形状，用基于序列图像的裂缝网络优化裂隙网络，得到完整的裂隙网络。上述方法可对复杂裂隙图进行全自动批量处理，从而得到完整准确的裂隙网络，解决传统方法带来的局限性。

【技术实现步骤摘要】
一种复杂裂隙网络的裂缝提取及优化方法
本专利技术涉及裂隙网络处理
，尤其涉及一种复杂裂隙网络的裂缝提取及优化方法。
技术介绍
在煤、矿物、石油开采过程中，地质条件复杂、煤与瓦斯突出、瓦斯爆炸等灾害不时发生，带来了巨大的经济损失和大量人员伤亡。深部煤岩中的裂隙网络是瓦斯煤层气的存储场所和运移通道，开采和掘进过程导致煤岩体中应力场发生改变，进而影响了煤岩体中裂隙的延伸、扩展及裂隙网络的重新分布。煤炭生产需要构建流体运移模型，掌握煤岩中裂隙网络以及采动应力作用下的演化规律有重要的意义。采样样品中的裂缝提取及优化用于开采过程中稳定检测富有气体含量，例如检测煤中瓦斯、油田中天然气含量。但是提取裂缝过程存在大量困难，例如：裂缝图像不明显，与背景过于相似，不易提取裂缝信息；复杂裂隙图像中存在大量随机性噪声和颗粒性噪声，使得裂缝提取十分困难，复杂裂隙网络的准确提取，直接关系到后续煤岩的深入研究。现有技术主要用人工勾画方式提取复杂裂隙图像中的裂隙网络，准确性受人主观影响大，提取效率低下；而数字图像技术从复杂裂隙图中提取裂隙网络，主要基于阈值分割、边缘检测、区域生长等方法，这些方法不适应对比度低、细节多、边界弱等特性的复杂裂隙图像，且具有不能抑制噪声、模糊图像细节、易受干扰、精确度低等缺点。
技术实现思路
本专利技术的目的是提供一种复杂裂隙网络的裂缝提取及优化方法，该方法可对复杂裂隙图进行全自动批量处理，从而得到完整准确的裂隙网络，解决传统方法带来的局限性，得到了更加理想的裂缝提取结果。本专利技术的目的是通过以下技术方案实现的：一种复杂裂隙网络的裂缝提取及优化方法，所述包括：步骤1、采用直方图均衡化、高斯滤波方式和伽马变换算法对待处理的复杂裂隙网络进行处理，实现图像中复杂裂缝与基质像素的分离；步骤2、采用自适应otsu阈值化初步提取裂缝目标，然后使用基于密度变换去除二值图像中基质材料处的噪声，保留裂缝部位，实现复杂裂缝网络的初步提取；步骤3、建立基于形态特征的伪目标去除模型，从步骤2获得的复杂裂隙网络数据中提取出清晰图像帧的裂隙网络；步骤4、建立序列图像的裂缝网络相似度优化模型，将步骤3提取出的清晰图像帧的裂隙网络作为先验形状，得到其他相邻图像的复杂裂隙网络，最终得到完整的复杂裂隙网络。由上述本专利技术提供的技术方案可以看出，上述方法可对复杂裂隙图进行全自动批量处理，从而得到完整准确的裂隙网络，解决传统方法带来的局限性，得到了更加理想的裂缝提取结果。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。图1为本专利技术实施例提供的复杂裂隙网络的裂缝提取及优化方法流程示意图；图2为本专利技术实施例所举的裂隙图像原图；图3为本专利技术实施例所述经过直方图均衡化处理后的图像；图4为本专利技术实施例所述经过高斯滤波处理后的图像；图5为本专利技术实施例所述经过二值化处理后的结果图；图6为本专利技术实施例所述经过序列图像中裂缝网络优化后的结果图。具体实施方式下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术的保护范围。下面将结合附图对本专利技术实施例作进一步地详细描述，如图1所示为本专利技术实施例提供的复杂裂隙网络的裂缝提取及优化方法流程示意图，所述方法包括：步骤1、采用直方图均衡化、高斯滤波方式和伽马变换算法对待处理的复杂裂隙网络进行处理，实现图像中复杂裂缝与基质像素的分离；1.1采用直方图均衡化处理的过程为：设待处理的复杂裂隙网络的原图像f(x,y)的像素点总数为n，有m个灰度级；用r表示原图像的灰度值，s表示经过直方图均衡化处理后图像的灰度值，nk表示灰度级为rk的像素点个数，灰度级k出现的概率如下式所示：直方图均衡化处理后各像素的灰度值表示为：举例来说，如图2所示为本专利技术实施例所举的裂隙图像的原图，如图3所示为本专利技术实施例所述经过直方图均衡化处理后的图像。经过直方图均衡化处理后，加大了裂缝与周边区域之间的差距，突出裂缝的信息，改善了煤岩序列CT图像的视觉效果，提高图像的清晰度，为图像的分割和进一步处理奠定良好的基础。1.2使用高斯滤波的方式去除随机性噪声，高斯滤波器是一种线性滤波器，二维高斯滤波器对目标像素某邻域内像素值分别赋以不同权值，对图像卷积后得到目标点平滑值。高斯权值由两个像素之间的欧式距离决定，观察正态分布曲线可以得知，在一定邻域范围内，与目标点距离越近的像素影响越大，反之影响越小。具体实现过程中，二维高斯函数公式如下：其中，σ是标准差，即高斯滤波器的尺度因子，决定高斯滤波器的范围，σ越大，滤波效果越明显。模板范围越大，滤波效果越明显，通过调节标准差σ，可以防止图像中有效的特征分量被模糊，并且平滑图像中影响图像效果的噪声以及细纹理。高斯卷积如下：其中：g[i,j]为高斯卷积模板，f[i,j]为原始输入图。举例来说，如图4所示为本专利技术实施例所述经过高斯滤波处理后的图像，通过高斯滤波处理可以抑制图像噪声，平滑图像。1.3基于伽马变换的裂隙边界增强算法，增强裂隙图像的裂缝信息、矫正裂隙图像的亮度，获得处理后的复杂裂隙网络数据；该步骤中，由于显卡或者显示器的原因会出现实际输出的图像在亮度上有偏差，而γ曲线矫正就是通过一定的方法来矫正图像的这种偏差的方法。一般情况下，当用于γ矫正的值大于1时，图像的高光部分被压缩而暗调部分被扩展，当γ矫正的值小于1时，图像的高光部分被扩展而暗调部分被压缩，γ矫正一般用于平滑的扩展暗调的细节。具体实现中，伽马变换的公式表示为：fgamma(x,y,z)＝Af′log(x,y,z)γ；其中，f′log(x,y,z)为符合要求的复杂裂隙图像数据；fgamma(x,y,z)为伽马变换输出的图像数据；A为常数，γ是可调变量，用以控制图像增强程度。步骤2、采用自适应otsu阈值化初步提取裂缝目标，然后使用基于密度变换去除二值图像中基质材料处的噪声，保留裂缝部位，实现复杂裂缝网络的初步提取；在该步骤中，采用双窗OTSU二值法对步骤1处理后的复杂裂隙网络数据进行处理，得到隐含完整裂隙网络的二值化结果，如图5所示为本专利技术实施例所述经过二值化处理后的结果图。其中，复杂裂隙网络的像素值用f(x,y)表示，经过密度过滤器去除噪声后的图像用g(x,y)表示，其中(x,y)为像素点坐标，过滤器的尺寸选取为a×a，则有：gx,y＝1&i＝x-a2x+a2j＝y-a2y+a2fi,ja2>70％0&否则&&fx,y＝1；再利用轮廓查找算法获得所有轮廓，用八链码原理将每个轮廓边界编码为一组被称为链码的代码，用0～7表示8个方向，边界链码组合用于形状面积和周长的计算。步骤3、建立基于形态特征的伪目标去除模型，从步骤2获得的复杂裂隙网络数据中提取出清晰图像帧的裂隙网络；3.1所述形态特征包括：(1)面积S：直接用区域的像素之和作为区域的面积；(2)本文档来自技高网...

【技术保护点】
1.一种复杂裂隙网络的裂缝提取及优化方法，其特征在于，所述方法包括：步骤1、采用直方图均衡化、高斯滤波方式和伽马变换算法对待处理的复杂裂隙网络进行处理，实现图像中复杂裂缝与基质像素的分离；步骤2、采用自适应otsu阈值化初步提取裂缝目标，然后使用基于密度变换去除二值图像中基质材料处的噪声，保留裂缝部位，实现复杂裂缝网络的初步提取；步骤3、建立基于形态特征的伪目标去除模型，从步骤2获得的复杂裂隙网络数据中提取出清晰图像帧的裂隙网络；步骤4、用建立序列图像的裂缝网络相似度优化模型，将步骤3提取出的清晰图像帧的裂隙网络作为先验形状，得到其他相邻图像的复杂裂隙网络，最终得到完整的复杂裂隙网络。

【技术特征摘要】
1.一种复杂裂隙网络的裂缝提取及优化方法，其特征在于，所述方法包括：步骤1、采用直方图均衡化、高斯滤波方式和伽马变换算法对待处理的复杂裂隙网络进行处理，实现图像中复杂裂缝与基质像素的分离；步骤2、采用自适应otsu阈值化初步提取裂缝目标，然后使用基于密度变换去除二值图像中基质材料处的噪声，保留裂缝部位，实现复杂裂缝网络的初步提取；步骤3、建立基于形态特征的伪目标去除模型，从步骤2获得的复杂裂隙网络数据中提取出清晰图像帧的裂隙网络；步骤4、用建立序列图像的裂缝网络相似度优化模型，将步骤3提取出的清晰图像帧的裂隙网络作为先验形状，得到其他相邻图像的复杂裂隙网络，最终得到完整的复杂裂隙网络。2.根据权利要求1所述复杂裂隙网络的裂缝提取及优化方法，其特征在于，所述步骤2的过程为：采用双窗OTSU二值法对复杂裂隙网络数据进行处理，得到隐含完整裂隙网络的二值化结果；其中，复杂裂隙网络的像素值用f(x,y)表示，经过密度过滤器去除噪声后的图像用g(x,y)表示，其中(x,y)为像素点坐标，过滤器的尺寸选取为a×a，则有：gx,y＝1&i＝x-a2x+a2j＝y-a2y+a2fi,ja2>70％0&否则&&fx,y＝1；利用轮廓查找算法获得所有轮廓，用八链码原理将每个轮廓边界编码为一组被称为链码的代码，用0～7表示8个方向，边界链码组合用于形状面积和周长的计算。3.根据权利要求1所述复杂裂隙网络的裂缝提取及优化方法，其特征在于，在步骤3中，所述形态特征包括：(1)面积S：直接用区域的像素之和作为区...

【专利技术属性】
技术研发人员：张国英，许文祥，
申请(专利权)人：中国矿业大学北京，
类型：发明
国别省市：北京,11