[发明专利]矿区线性异常分析提取新型方法

说明书

本发明公开了矿区线性异常分析提取新型方法。初步线性特征提取；对所得线性特征进行模式分析和滤波处理，首先对所得图像进行区域分割，将彩色图像转换为灰度图像；对所得图像进行二值化，应用高斯低通滤波器去除图像噪声，对去噪后所得二值图像进行形态学打开操作，去除区域中的小对象，将细化区域取反，作为噪声，继续精细图像中的线型区域；对提取得到的噪声线型图像进行形态学膨胀操作，在区域图像上剔除掉提取到的所有短线型，得到图像中的有效线型。本发明的有益效果是对具有线性特征的异常数据进行线型提取处理取得了较好的结果。

技术领域

本发明属于地质技术领域，涉及矿区线性异常分析提取新型方法。

背景技术

对大洋构造特征进行详细解译和研究，对热液矿床成因、大小等方面预测及更广泛的研究和应用具有重要的科学意义。目前，尽管对大西洋脊尤其是北大西洋脊区的热液矿点及构造背景已有较多的发现和研究，但对大洋构造特征的掌握还不够精细和全面。对调查资料详尽的热液区域，没有详细的构造解译和研究，也就不可能有充分的进一步的解析和利用。由于地壳变动或物质喷发及填充形成的地形起伏除了具有趋势性，还具有局部区域的不均一性，传统的地形数据解译一般通过等高线或灰度色系等方式体现，这种方法因为是对区域数据的总体数值统计，很容易只体现大的地形趋势性，忽略掉小的地形起伏异常，而复杂的成矿环境中，小的地形异常恰可能是矿床或矿点的结果和标志。重力和磁力信息反映地质体的物理特性，如块的大小和密度，可以确定地质体的位置、地质体区域边界，岩性及形态特征，能够识别隐伏地质体(HarrisandPan,1991)。有线性构造的位置是地壳局部相对较薄、整体压力相对较低区域，其物质密度区别于周边非构造区，另外线性构造也是岩浆喷发的源区，并为喷发物质提供了就位空间，构造区地质体密度、大小、岩性等物理特征的差异会在重力和磁力信息中得以体现，且能够隐含地质体的新老关系。利用重力和磁力数据作为地形数据的协同要素，对地形线性特征的解译，构造的期次、方向等参数的确定具有重要意义，也进而能够更深入地解释和发现热液的控矿构造特征。构造特征的提取在于构造边界的识别和线性纹理的提取。这里运用能够指示地形突变的梯度指向累计、线性纹理提取等方法进行线性构造特征的提取。

发明内容

本发明的目的在于提供矿区线性异常分析提取新型方法。

本发明所采用的技术方案是按照以下步骤进行：

步骤1：线性特征提取；

步骤2：对所得线性特征进行模式分析和滤波处理，首先对所得图像进行区域分割，然后对线型区域进行4倍下采样；

步骤3：将彩色图像转换为灰度图像；

步骤4：水平和竖直方向采用Canny算子提取图像边框，进行图像切割；

步骤5：对所得图像进行二值化，应用高斯低通滤波器去除图像噪声，对去噪后所得二值图像进行形态学打开操作，去除区域中的小对象，采用thin，remove两种形态学操作对图像中物体进行细化，将细化区域取反，作为噪声，继续精细图像中的线型区域；

步骤6：对提取得到的噪声线型图像进行形态学膨胀操作，在区域图像上剔除掉提取到的所有短线型，得到图像中的有效线型。

进一步，步骤1中利用ArcGIS空间分析中的Hydrology方法实现线性异常的提取，将数据根据精度划分成适当的网格，生成具有空间坐标和异常属性值的栅格图层，首先利用FlowDirection方法，将准备好的栅格图层作为输入图层，计算每个单元的梯度方向，输出图层命名为Dir，则得到指示梯度方向的图层；用FlowAccumulation方法将梯度方向图层作为输入图层，计算网格单元的梯度指向累计次数，输出图层命名为Acc，利用RasterCalculator计算梯度指向次数为0的网格单元，并输出图层Acc0，将Acc0重分类，重分类图层命名为Acc0Rec，图层中值为1的网格形成的区域即为异则获得初步的线性异常特征。

本发明的有益效果是对原始数据进行线型提取取得了较好的结果。