一种探测地表采动裂缝的方法技术

本发明专利技术涉及一种探测地表采动裂缝的方法，包括以下步骤：首先，建立地表采动裂缝的解译标志；然后，获取研究区的遥感影像与采空区矢量图，对遥感影像进行预处理及分类，得到矿井口位置，作为佐证采空区和工作面的逻辑信息；将遥感影像、采空区矢量图与工作面矢量图叠加，生成缓冲区图，划定地表采动裂缝存在范围并通过裁剪提取该区域；根据解译标志对该区域进行高通滤波和方向滤波，弱化非地表采动裂缝；再利用区域生长法得到疑似地表采动裂缝，对其进行概率分级，获取最终地表采动裂缝，采集其坐标位置。该方法可精准识别地表采动裂缝，为土地复垦、生态修复提供重要基础数据，也可为损害赔偿、裂缝治理提供科学依据和技术支撑。

A Method for Detecting Surface Mining Cracks

【技术实现步骤摘要】
一种探测地表采动裂缝的方法
本专利技术涉及矿山生态环境监测的
，尤其是涉及一种探测地表采动裂缝的方法。
技术介绍
地表采动裂缝是地裂缝中一类较为特殊的类型，其形成主要是由地下采矿引起，易诱发次生地质灾害，改变微地形微地貌，扰动生态系统，不但制约了煤炭工业与经济的发展，而且对矿区建设和人民生命财产以及生态环境也造成了严重的危害。对煤自燃区域生态环境进行监测及修复，需对矿区的地表采动裂缝进行探测。目前，对于地裂缝的监测主要是利用地裂缝监测装置对其进行实地直接测量。如，《一种井工开采地表裂缝动态监测设备及方法》，专利技术公开专利号为CN109084726A；《一种SAA监测地裂缝的方法及系统》，专利技术公开专利号为CN109059827A；《一种采动裂缝监测装置》，专利技术公开专利号为CN204188107U；《一种探测浅层地裂缝的方法》，专利技术公开专利号为CN106507971B；《一种智能拉线裂缝监测仪》，专利技术公开专利号为CN109373855A。这些专利技术适合小范围地裂缝监测，对于大范围监测适用性有限。针对因采矿引起的地表采动裂缝的大范围精准监测，目前仍缺乏有效解决办法。
技术实现思路
针对上述问题，本专利技术提供了一种探测地表采动裂缝的方法，利用遥感技术具有覆盖范围广、获取便捷且成本较低的特点，将其应用于地表采动裂缝的探测方面，为满足上述需求提供了一种可行性方案，对保证矿区建设和人民生命财产、修复矿区生态环境提供了有力的技术支撑。本专利技术采用的技术方案是一种探测地表采动裂缝的方法，包括以下步骤：步骤1：建立地表采动裂缝影像解译标志；(1)裂缝发育方向基本与工作面方向平行；(2)裂缝呈现直线型，且分叉较少；(3)裂缝长宽比数值在30至60之间；(4)地缝对整块地物产生切割效果；步骤2：获取研究区的遥感影像、采空区与工作面矢量图，对影像进行预处理，得到预处理后影像；利用支持向量机的方法对预处理后影像进行分类，得到分类后影像，再对分类后影像进行矢量化；利用采空区与工作面矢量图生成缓冲区范围图，将分类后影像与缓冲区范围图叠加，重合区域即为地表采动裂缝存在区域,提取其重合区域图像，记为f。其过程主要包括步骤2-1至步骤2-3。步骤2-1：对遥感影像进行预处理，其过程主要包括坐标变换、辐射标定和影像融合；步骤2-2：利用支持向量机方法对预处理后遥感影像进行分类，分为矿井厂房、阴影、裸土、积雪四类训练样本，对分类后影像进行矢量化，得到矿井厂房信息，以确保采空区与工作面位置正确；步骤2-3：利用采空区与工作面矢量图作为面图层，缓冲区范围设为100米，生成缓冲区范围图，叠加分类后遥感影像与缓冲区范围图，重合部分即为地表采动裂缝存在区域，通过裁剪提取该区域影像，记为f；步骤3：对地表采动裂缝存在区域影像f进行高通滤波和方向滤波处理；首先，通过高通滤波器实现对线状地物的提取，再通过方向滤波弱化非地表采动裂缝，得到滤波后影像，该影像记为g，相关公式如下：RH为高通滤波结果，R为原始影像，GH为高通滤波器，RD为方向滤波别结果，GD为方向滤波器；步骤4：分别对滤波后影像g中的非线状地物区域和线状地物区域进行区域生长，两者求差获取疑似地表采动裂缝Rc，主要步骤及相关公式如下：步骤4-1：获取像素点的颜色特征和纹理特征作为区域生长的种子点特征信息，假设彩色图像为f(x，y)，x＝1，2，...，M；y＝1，2，...，N，Pxy为该图像像素点；其中Col(i，j)、HH(i，j)分别为像素点(i，j)颜色特征和纹理强度信息；H(i，j)、S(i，j)、V(i，j)依次表示图像中(i，j)像素归一化处理后的H(色调)、S(饱和度)、V(明度)三个分量值；w1、w2、w3表示三个权重的系数；步骤4-2：求取像素(i，j)和某一个区域颜色和纹理均值的距离；DCol(i，j；avg)＝[w1[H(i，j)-Havg]2+w2[S(i，j)-Savg]2+w3[V(i，j)-Vavg]2]1/2(公式4)DGra(i，j；avg)＝HH(i，j)-HHavg(公式5)上式中，DCol(i，j；avg)、DGra(i，j；avg)分别为像素(i，j)和某一个区域颜色和纹理均值的距离；HHavg、Savg、Vavg分别为已加入到生长区域中像素点的H、S、V均值；HHavg表示已经加入到生长区域中像素点的纹理强度信息均值；步骤4-3：获取区域分割阈值条件，DF(i，j；avg)＜T(公式6)T表示区域分割的阈值条件。步骤4-4：获取疑似地表采动裂缝结果Rc，Rc＝|RG1-RG2|(公式7)RG1是第一次非线状地物区域生长的结果，RG2是第二次线状地物区域生长结果，RC是两次区域生长求差结果，即疑似地表采动裂缝；步骤5：将疑似地表采动裂缝Rc进行概率等级划分，根据地表采动裂缝解译标志，将结果Rc分成四个等级，即地表采动裂缝的存在概率等级；符合其中一个条件划分为Ⅳ级地表采动裂缝，符合两个条件划分为Ⅲ级地表采动裂缝，符合三个条件划分为Ⅱ级地表采动裂缝，四个判断条件全部满足则划分为Ⅰ级地表采动裂缝，本专利技术取Ⅰ级地表采动裂缝为最终地表采动裂缝识别结果；步骤6：进行地表采动裂缝地理坐标采集，识别出来的地表采动裂缝为不规则曲线，采集3个或3个以上坐标点作为其地理坐标。本专利技术的技术效果在于：(1)利用地表采动裂缝特征及其与开采工作面、采空区位置相对应关系，实现了对地表采动裂缝的精准探测，降低了地表裂隙、人工陡坎及雨水裂缝等易错分地物对探测结果的影响，保证了地表采动裂缝探测结果的可靠性。(2)利用遥感影像进行地表采动裂缝的精准探测，可实现地表采动裂缝的大范围精准监测，并且提高了工作效率，降低了生产成本。(3)本专利技术最终提供地表采动裂缝探测的地图，可以为煤炭工业与经济的发展，对地表采动裂缝区域生态环境的监测及精准修复提供有力依据。附图说明图1为本专利技术的技术流程图。图2为矿井厂房矢量化图。图3为实施例缓冲区范围图。图4为实施例的地表采动裂缝探测结果图。图5为实施例的地表采动裂缝与真彩色影像叠加图。图6为实施例的地表采动裂缝上点位坐标获取图。具体实施方式以下结合具体实施例对本专利技术作进一步详细描述。在本实施例中，采用高空间分辨率WorldView-2卫星数据作为探测数据源，并结合图1进行详细描述，具体步骤如下：步骤1：建立地表采动裂缝影像解译标志。(1)裂缝发育方向基本与工作面方向平行；(2)裂缝呈现直线型，且分叉较少；(3)裂缝长宽比数值在30至60之间；(4)裂缝对整块地物产生切割效果。步骤2：选取山西省忻州市保德县康家塔村、贾家梁村、路家沟村所在的保德矿区作为研究区域，获取研究区的遥感影像、采空区与工作面矢量图，对影像进行预处理，得到预处理后影像。利用支持向量机的方法对预处理后影像进行分类，得到分类后影像，再对分类后影像进行矢量化，图2为矿井厂房矢量图。利用采空区与工作面矢量图生成缓冲区范围图，本实施例中，研究区缓冲区范围图如图3所示，再将分类后影像与缓冲区范围图叠加，重合区域即为地表采动裂缝存在区域,提取其重合区域图像，记为f。其过程主要包括步骤2-1至步骤2-3。步骤2-1：对遥感影像进行预处理，其过程主要包括坐标变换、辐射标定和影像融合。本实施例中利用ENVI软件建本文档来自技高网...

【技术保护点】
1.一种探测地表采动裂缝的方法，其特征在于，包括以下步骤：步骤1：建立地表采动裂缝影像解译标志；步骤2：获取研究区的遥感影像、采空区与工作面矢量图，对影像进行预处理，得到预处理后影像；利用支持向量机的方法对预处理后影像进行分类，得到分类后影像，再对分类后影像进行矢量化；利用采空区与工作面矢量图生成缓冲区范围图，将分类后影像与缓冲区范围图叠加，重合区域即为地表采动裂缝存在区域,提取其重合区域图像，记为f；步骤3：对地表采动裂缝存在区域影像f进行高通滤波和方向滤波处理；首先，通过高通滤波器实现对线状地物的提取，再通过方向滤波弱化非地表采动裂缝，得到滤波后影像，该影像记为g，相关公式如下：

【技术特征摘要】
1.一种探测地表采动裂缝的方法，其特征在于，包括以下步骤：步骤1：建立地表采动裂缝影像解译标志；步骤2：获取研究区的遥感影像、采空区与工作面矢量图，对影像进行预处理，得到预处理后影像；利用支持向量机的方法对预处理后影像进行分类，得到分类后影像，再对分类后影像进行矢量化；利用采空区与工作面矢量图生成缓冲区范围图，将分类后影像与缓冲区范围图叠加，重合区域即为地表采动裂缝存在区域,提取其重合区域图像，记为f；步骤3：对地表采动裂缝存在区域影像f进行高通滤波和方向滤波处理；首先，通过高通滤波器实现对线状地物的提取，再通过方向滤波弱化非地表采动裂缝，得到滤波后影像，该影像记为g，相关公式如下：RH为高通滤波结果，R为原始影像，GH为高通滤波器，RD为方向滤波结果，GD为方向滤波器；步骤4：分别对滤波后影像g中的非线状地物区域和线状地物区域进行区域生长，两者求差获取疑似地表采动裂缝Rc，主要步骤及相关公式如下：步骤4-1：获取像素点的颜色特征和纹理特征作为区域生长的种子点特征信息，假设彩色图像为f(x，y)，x＝1，2，…，M；y＝1，2，…，N，Pxy为该图像像素点；其中Col(i，j)、HH(i，j)分别为像素点(i，j)颜色特征和纹理强度信息；H(i，j)、S(i，j)、V(i，j)依次表示图像中(i，j)像素归一化处理后的H(色调)、S(饱和度)、V(明度)三个分量值；w1、w2、w3表示三个权重的系数；步骤4-2：求取像素(i，j)和某一个区域颜色和纹理均值的距离；DCol(i，j；avg)＝[w1[H(i，j)-Havg]2+w2[S(i，j)-Savg]2+w3[V(i，j)-Vavg]2]1/2(公式4)DGra(i，j；avg)＝HH(i，j)-HHavg(公式5)上式中，DCol(i，j；avg)、DGra(i，j；avg)分别为像素(i，j)和某一个区域颜色和纹理均值的距离；HH...

【专利技术属性】
技术研发人员：朱笑笑，杨永均，侯湖平，张绍良，张晋纶，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏,32