CN108492260A - 基于张量投票耦合霍夫变换的地质线性体提取方法 - Google Patents
基于张量投票耦合霍夫变换的地质线性体提取方法 Download PDFInfo
- Publication number
- CN108492260A CN108492260A CN201810121330.4A CN201810121330A CN108492260A CN 108492260 A CN108492260 A CN 108492260A CN 201810121330 A CN201810121330 A CN 201810121330A CN 108492260 A CN108492260 A CN 108492260A
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing image
- image
- coordinate system
- lineament
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000009466 transformation Effects 0.000 title claims abstract description 23
- 230000002708 enhancing effect Effects 0.000 claims abstract description 18
- 238000001914 filtration Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 16
- 239000002131 composite material Substances 0.000 claims abstract description 15
- 230000008878 coupling Effects 0.000 claims abstract description 9
- 238000010168 coupling process Methods 0.000 claims abstract description 9
- 238000005859 coupling reaction Methods 0.000 claims abstract description 9
- 238000003708 edge detection Methods 0.000 claims abstract description 9
- 238000002203 pretreatment Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000000605 extraction Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000012892 rational function Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000006870 function Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000003707 image sharpening Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G06T5/73—
-
- 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/143—Segmentation; Edge detection involving probabilistic approaches, e.g. Markov random field [MRF] modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- 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/10032—Satellite or aerial image; Remote sensing
-
- 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/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral 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
-
- 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/20048—Transform domain processing
- G06T2207/20061—Hough transform
-
- 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/20076—Probabilistic image processing
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Software Systems (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提供了一种基于张量投票耦合霍夫变换的地质线性体提取方法，包括以下步骤：对遥感影像进行预处理，得到预处理后的遥感影像；在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，得到由最优波段组合的遥感影像，其中N为大于等于3的自然数；利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，增强线性化边缘信息；对增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点；将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体。本发明对解译专家的知识和经验依赖更少，极大程度上缩短了处理的时间，节省了大量的人力，具有较强的实用价值和推广意义。
Description
技术领域
本发明属于遥感地质及影像处理技术领域，具体涉及一种基于张量投票耦合霍夫变换的地质线性体提取方法。
背景技术
断裂带、断层等地质构造属于地质薄弱区，由于侵蚀等作用形成线状地貌。在地质构造的作用下，这些线状地貌在遥感影像上往往表现出明显的线状分布，称为地质线性体，其控制着地下流体(矿液、地下水和油气等)的运移和矿产资源的空间赋存，其方位与密度对分析区域构造运动趋势与活动程度有着深远的科学意义和实用价值。
从上世纪八十年代以来，国内外不少遥感地质工作者利用专家知识经验与图像处理的方法，提取了遥感影像与数字高程模型(Digital elevation model,DEM)的地质线性体，分析了区域地质构造趋势及影响程度。如Jansson等(2005)使用Landsat 7ETM+和数字高程模型提取地质线性体，并绘制了Wales东北部冰川地貌；吴婧等(2011)使用Canny边缘检测和霍夫变换提取了断裂构造信息，并使用ENVI+IDL等编程语言实现；袁小祥等(2011)利用假彩色合成、主成分变换、缨帽变换、波段比值、地貌渲染等方法，突出了地质线性体在多源遥感数据上的反差，提取了活动断层；Alaa,A.M等(2011)根据DEM数据山体阴影增强可以线性特征的特点，使用线索追踪算法提取山体阴影渲染后的影像上的地质线性体，并用B-spline曲线整合表达，最终评价了区域地质构造环境状况；Yusof等(2011)分析了高速公路周围滑坡灾害分布与地质线性体密度之间的关系；刘智荣等(2012)利用反差增强、彩色合成、方向滤波与图像融合等图像增强处理算法，提取与分析了银川活动断层信息及分布；Bahiru等(2016)使用Landsat ETM+和SRTMDEM数据提取并绘制了乌干达地区地质线性构造分布，研究了该区域金矿的分布，对矿产预测有非常重要的研究意义。
上述工作都不同程度地取得了一些成就，主要存在以下不足：
1)目视解译结果的正确性依赖于解译专家的经验和知识背景，耗时耗力且效率低。
2)计算机解译的精度与处理速度、数据源的分辨率等因素有关，影像越大，处理速度就慢；分辨率过高，则容易受到道路、土地利用界线等线性特征影响，而产生大量的错误线性边缘，也就产生了过多的噪声。过度依赖参数设置，导致方普适性差。
发明内容
针对现有技术存在的不足，本发明的目的在于，提供一种基于张量投票耦合霍夫变换的地质线性体提取方法，解决现有技术中依赖于解译专家的经验和知识背景，耗时耗力且效率低、处理速度慢、噪声大及普适性差的问题。
为了解决上述技术问题，本发明采用如下技术方案予以实现：
一种基于张量投票耦合霍夫变换的地质线性体提取方法，包括以下步骤：
步骤1，对遥感影像进行预处理，得到预处理后的遥感影像；
步骤2，在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，得到由最优波段组合的遥感影像，N为大于等于3的自然数；
步骤3，利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，增强线性化边缘信息；
步骤4，对步骤3增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点；
步骤5，将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体；
所述图像坐标系为以遥感影像的任一角作为原点，遥感影像的水平方向为x轴，遥感影像的竖直方向为y轴；
所述参数坐标系为ρ＝x cosθ+y sinθ，其中θ和ρ为边缘点对应的参数坐标系下的极坐标。
进一步地，所述步骤1中对遥感影像进行预处理，包括：
从遥感影像库中任选一遥感影像作为当前遥感影像，该当前遥感影像上含云量低于5％，且具备有理函数模型及统计影像灰度信息；所述影像灰度信息包括每个波段的灰度方差和标准差；
对该当前遥感影像进行辐射定标、大气校正、影像裁剪处理。
进一步地，所述步骤2中在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，包括：
选择预处理后的遥感影像的多光谱波段中波段指数OIF最大时所对应的3个波段作为最优波段；
其中，通过式(1)得到波段指数OIF：
式(1)中，Si为第i个波段的标准差；Rij为第i个波段和第j个波段的相关系数；i≠j，i＝1,2,...,N,j＝1,2,...,N；
D(i)、D(j)分别为第i个波段和第j个波段的方差；Cov(i,j)为第i个波段和第j个波段之间的协方差。
进一步地，所述步骤3中利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，包括：
通过式(2)对由最优波段组合的遥感影像进行高斯高通滤波，得到滤波后的遥感图像H(u,v)：
式(2)中，D(u,v)是由最优波段组合的遥感影像的频率域中点(u,v)与频率矩形中心的距离，D0为常数。
进一步地，所述步骤4中对步骤3增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点，包括：
步骤41，任选增强线性化边缘信息的遥感影像中的一个像素点作为当前像素点，对当前像素点作拉普拉斯算子卷积操作，得到张量矩阵T；
其中，I为增强线性化边缘信息的遥感影像，和别为遥感影像I沿x和y方向的二阶导数；其中，x、y是以遥感影像I的任一角作为原点，遥感影像I的水平方向为x轴，遥感影像I的竖直方向为y轴建立的图像坐标系中的x轴和y轴；
步骤42，将张量矩阵T进行矩阵谱分解，得到棒状分量和球状分量；
其中， 为棒状分量，为球状分量；
步骤43，当(λ1-λ2)＞λ2时，所述当前像素点为边缘点；否则，所述当前像素点为非边缘点；
步骤44，重复步骤41至步骤43，直至增强线性化边缘信息的遥感影像的所有像素点都作为当前像素点为止，得到遥感影像中的所有边缘点。
进一步地，所述步骤5中将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体，包括：
步骤51，遍历参数坐标系，寻找局部最大值的点，将该局部最大值的点作为峰值点，设峰值点的坐标为(ρ，θ)，其中(ρ，θ)为遥感影像中地质线性体的斜率和截距；
步骤52，将参数坐标系下峰值点所对应的坐标转换到图像坐标系中，按照边缘点的方向及端点距离连接成线，得到地质线性体影像，即完成地质线性体的提取。本发明与现有技术相比，具有如下技术效果：
1.本发明综合了波段选择、图像增强、边界检测、线性提取等算法与规则，提供了基于张量投票耦合霍夫变换的地质线性体提取方法，相对于单纯的目视解译等方法，对解译专家的知识和经验依赖更少，极大程度上缩短了处理的时间，节省了大量的人力，具有较强的实用价值和推广意义；
2.较Canny边缘检测算法，基于张量投票的边缘检测算法既可以在边缘检测的基础上进行边界检测，也可以在灰度图像上直接用二维圆形投票域进行张量投票过程，然后进行投票解释，并根据边界特征的显著性提取边界，具有鲁棒性；
3.本发明可以处理多源遥感数据，方法的参数设置更为均衡，具有较好的普适性，对区域构造演化及板块运动方向有较强的指示作用。
附图说明
图1为本发明的基于张量投票耦合霍夫变换的地质线性体提取方法的总体流程图；
图2为本发明实例的研究区7-5-4波段合成遥感影像；
图3为本发明实例的研究区高斯高通滤波后增强影像；
图4为本发明张量投票的示意图；
图5为本发明实例的研究区张量投票后边缘检测影像；
图6为本发明霍夫变换的示意图；
图7为本发明地质线性体提取影像。
以下结合附图对本发明的具体内容作进一步详细解释说明。
具体实施方式
本发明中的遥感影像为Landsat 8OLI冬季多光谱遥感影像。
以下给出本发明的具体实施例，需要说明的是本发明并不局限于以下具体实施例，凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
如图1所示，本实施例提供了基于张量投票耦合霍夫变换的地质线性体提取方法，包括以下步骤：
步骤1，对遥感影像进行预处理，得到预处理后的遥感影像；
具体地，本实施例从遥感影像库中任选一遥感影像作为当前遥感影像，该当前遥感影像上含云量低于5％，且具备有理函数模型及统计影像灰度信息；所述影像灰度信息包括每个波段的灰度方差和标准差；
对该当前遥感影像进行辐射定标、大气校正、影像裁剪处理。
其中，本实施例中的遥感影像取Landsat 8OLI冬季多光谱遥感影像，因为冬季影像植被覆盖较少，线性边缘在影像上突出明显。影像上云含量应低于5％，具备有理函数模型(Rational Polynomial Coefficients，RPC)及统计影像灰度信息，以便与地质资料比对，验证预处理的精度；
影像灰度统计信息主要包括，每个波段的灰度方差Di、标准差Si，统计信息如下表所示，用于后续分析。
表1 Landsat 8OLI多光谱遥感影像
步骤2，在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，得到由最优波段组合的遥感影像，N为大于等于3的自然数；
具体地，本实施例选择预处理后的遥感影像的多光谱波段中波段指数OIF最大时所对应的3个波段作为最优波段；
其中，通过式(1)得到波段指数OIF：
式(1)中，Si为第i个波段的标准差；Rij为第i个波段和第j个波段的相关系数；i≠j，i＝1,2,...,N,j＝1,2,...,N；
D(i)、D(j)分别为第i个波段和第j个波段的方差；Cov(i,j)为第i个波段和第j个波段之间的协方差。
本实施例中利用Landsat 8OLI的遥感影像的6个多光谱波段，计算波段相关系数矩阵，如表2所示。
表2研究区Landsat 8OLI波段间相关系数矩阵
其中，相关系数越大说明，波段间冗余信息越多，对影像的波段之间的相关系数进行分析，不利于突出不同地质体的信息。表2中Band 5与其他波段的相关系数普遍较小，因此影像合成可优先考虑Band 5。Band 7为短波红外波段，对岩石、特定矿物反应敏感，用于区分主要岩石类型、岩石水热蚀变，探测相关粘土矿物等，其次优先考虑Band 7。
遥感影像的波段组合如表3所示，从表3中可以看出，7-5-4组合方式与7-5-6的组合方式是最大的，但Band分别与Band 4、Band 6的相关系数为0.447和0.417，因此应选择7-5-4波段合成方式。最终确定7、5、4波段分别赋予红、绿、蓝合成得到的影像，如图2，增强各种地物间的对比度。即图2为由最优波段组合的遥感影像。
表3 Landsat 8OLI波段组合和相应的OIF指数
此外，去除多余波段，目的在于抵制其余波段对合成影像的干扰，如使用全部波段，任意两个波段的相关系数一旦过高，必会造成数据冗余。
步骤3，利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，增强线性化边缘信息；
具体地，通过式(2)对由最优波段组合的遥感影像进行高斯高通滤波，得到滤波后的遥感图像H(u,v)：
式(2)中，D(u,v)是由最优波段组合的遥感影像的频率域中点(u,v)与频率矩形中心的距离，D0为常数。本实施例中D0为截止频率，具体取20。
本实施例利用高斯高通滤波对遥感影像线性边缘的锐化作用，处理步骤2中得到的由最优波段组合的遥感影像，突出线性边缘，获取地质线性体提取的增强线性化边缘信息的遥感影像。
其优点在于：高通滤波是一种典型的图像锐化增强算子，采用高通滤波的方法可以让高频分量顺利通过、低频分量受到抑制，也就是增强边缘特征，抑制非边缘等噪声。而高斯高通滤波就是高通滤波的一种，增强了图像的对比度，使其具有更为明显的区分信息。
如图3所示为经过高斯高通滤波后的遥感影像，山谷线和山脊线等地质地貌线在影像上被突出，呈现明显的线性地貌。
步骤4，对步骤3增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点；
进一步的，步骤4的边界点矢量和叠加的显著性特征，定义如下：
其中，DF(L，k)显著性函数，L是曲线长度，k是曲率，σ是投票邻域范围，c是控制曲率衰减度的系数。
具体地，步骤41，任选增强线性化边缘信息的遥感影像中的一个像素点作为当前像素点，对当前像素点作拉普拉斯算子卷积操作，得到张量矩阵T；
其中，I为增强线性化边缘信息的遥感影像，和分别为遥感影像I沿x和y方向的二阶导数；其中，x、y是以遥感影像I的任一角作为原点，遥感影像I的水平方向为x轴，遥感影像I的竖直方向为y轴建立的图像坐标系中的x轴和y轴；
本实施例利用拉普拉斯算子计算张量矩阵T中的二阶导数，算子如下：
进行奇异值分解，可得：
步骤42，将张量矩阵T进行矩阵谱分解，得到棒状分量和球状分量；
其中， 为棒状分量，为球状分量；
步骤43，当(λ1-λ2)＞λ2时，所述当前像素点为边缘点；否则，所述当前像素点为非边缘点；
具体地，当λ1≈λ2时，该当前像素点位于一个区域内部点或者交叉处；当λ1和λ2取值都非常小，该当前像素点判断为非边缘点。
步骤44，重复步骤41至步骤43，直至增强线性化边缘信息的遥感影像的所有像素点都作为当前像素点为止，得到遥感影像中的所有边缘点。
如图4为张量投票的示意图，其中，线性边缘点由于其灰度突变性，经过8邻域范围点对其投票，矢量均有显著变化，具体表现为显著性函数大于0，因此将其视为边缘点且被保存；非线性边缘点由于邻域范围点为同一物质，经过邻域投票，矢量和无显著变化，具体表现为显著性函数约等于0；因此将其视为非边缘点且被删除。
本实施例的步骤4利用张量表示的待提取特征经过稀疏和多尺度密集投票，叠加矢量和，非地质线性体边界点得到不同方向的投票，而矢量相互作用而抵消，边界点由于只得到某一侧的投票而得到增强，最后通过投票解释过程获取区域的边界；利用显著性函数计算边界点叠加后矢量和的显著特征，实现张量投票方法的边缘检测，得到二值黑白影像，如图5所示。
步骤5，将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体；
本实施例使用霍夫变换将边缘点从图像坐标系转换到参数坐标系中，转而求参数坐标系的峰值，记录峰值对应的位置，提取地质线性体。
图像坐标系为以遥感影像的任一角作为原点，遥感影像的水平方向为x轴，遥感影像的竖直方向为y轴；
参数坐标系为ρ＝x cosθ+y sinθ，其中θ和ρ为边缘点对应的参数坐标系下的极坐标。
具体地，步骤51，遍历参数坐标系，寻找局部最大值的点，将该局部最大值的点作为峰值点，设峰值点的坐标为(ρ，θ)，其中(ρ，θ)为遥感影像中地质线性体的斜率和截距；
本实施例使用霍夫变换将边缘点从图像坐标系转换到参数坐标系，原始的边缘点坐标转换到了以(ρ，θ)为基准的参数坐标系，不断累加。遍历ρ-θ空间，寻找出局部最大值(极值)的点，称为峰值点。
图6从左至右分别为图像坐标系下的像素点通过霍夫变换转换到参数坐标下的示意图，再通过不断累加得到图6最右边中的一个峰值点。
步骤52，将参数坐标系下峰值点所对应的坐标转换到图像坐标系中，按照边缘点的方向及端点距离连接成线，得到地质线性体影像，即完成地质线性体的提取。如图7所示为叠加霍夫变换后的该区域地质线性体分布图。
Claims (6)
1.一种基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，包括以下步骤：
步骤1，对遥感影像进行预处理，得到预处理后的遥感影像；
步骤2，在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，得到由最优波段组合的遥感影像，N为大于等于3的自然数；
步骤3，利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，增强线性化边缘信息；
步骤4，对步骤3增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点；
步骤5，将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体；
所述图像坐标系为以遥感影像的任一角作为原点，遥感影像的水平方向为x轴，遥感影像的竖直方向为y轴；
所述参数坐标系为ρ＝x cosθ+y sinθ，其中θ和ρ为边缘点对应的参数坐标系下的极坐标。
2.根据权利要求1所述的基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，所述步骤1中对遥感影像进行预处理，包括：
从遥感影像库中任选一遥感影像作为当前遥感影像，该当前遥感影像上含云量低于5％，且具备有理函数模型及统计影像灰度信息；所述影像灰度信息包括每个波段的灰度方差和标准差；
对该当前遥感影像进行辐射定标、大气校正、影像裁剪处理。
3.根据权利要求1或2所述的基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，所述步骤2中在预处理后的遥感影像的N个多光谱波段中选择3个最优波段，包括：
选择预处理后的遥感影像的多光谱波段中波段指数OIF最大时所对应的3个波段作为最优波段；
其中，通过式(1)得到波段指数OIF：
式(1)中，Si为第i个波段的标准差；Rij为第i个波段和第j个波段的相关系数；i≠j，i＝1,2,...,N,j＝1,2,...,N；
D(i)、D(j)分别为第i个波段和第j个波段的方差；Cov(i,j)为第i个波段和第j个波段之间的协方差。
4.根据权利要求1或2所述的基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，所述步骤3中利用高斯高通滤波对由最优波段组合的遥感影像进行锐化处理，包括：
通过式(2)对由最优波段组合的遥感影像进行高斯高通滤波，得到滤波后的遥感图像H(u,v)：
式(2)中，D(u,v)是由最优波段组合的遥感影像的频率域中点(u,v)与频率矩形中心的距离，D0为常数。
5.根据权利要求1或2所述的基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，所述步骤4中对步骤3增强线性化边缘信息的遥感影像中进行边缘检测，得到遥感影像中的所有边缘点，包括：
步骤41，任选增强线性化边缘信息的遥感影像中的一个像素点作为当前像素点，对当前像素点作拉普拉斯算子卷积操作，得到张量矩阵T；
其中，I为增强线性化边缘信息的遥感影像，和分别为遥感影像I沿x和y方向的二阶导数；其中，x、y是以遥感影像I的任一角作为原点，遥感影像I的水平方向为x轴，遥感影像I的竖直方向为y轴建立的图像坐标系中的x轴和y轴；
步骤42，将张量矩阵T进行矩阵谱分解，得到棒状分量和球状分量；
其中， 为棒状分量，为球状分量；
步骤43，当(λ1-λ2)＞λ2时，所述当前像素点为边缘点；否则，所述当前像素点为非边缘点；
步骤44，重复步骤41至步骤43，直至增强线性化边缘信息的遥感影像的所有像素点都作为当前像素点为止，得到遥感影像中的所有边缘点。
6.根据权利要求1或2所述的基于张量投票耦合霍夫变换的地质线性体提取方法，其特征在于，所述步骤5中将遥感影像中的所有边缘点从图像坐标系转换到参数坐标系中，从参数坐标系中提取地质线性体，包括：
步骤51，遍历参数坐标系，寻找局部最大值的点，将该局部最大值的点作为峰值点，设峰值点的坐标为(ρ，θ)，其中(ρ，θ)为遥感影像中地质线性体的斜率和截距；
步骤52，将参数坐标系下峰值点所对应的坐标转换到图像坐标系中，按照边缘点的方向及端点距离连接成线，得到地质线性体影像，即完成地质线性体的提取。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810121330.4A CN108492260B (zh) | 2018-02-07 | 2018-02-07 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
| US16/627,900 US11430128B2 (en) | 2018-02-07 | 2018-09-17 | Geological linear body extraction method based on tensor voting coupled with Hough transformation |
| PCT/CN2018/105966 WO2019153746A1 (zh) | 2018-02-07 | 2018-09-17 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810121330.4A CN108492260B (zh) | 2018-02-07 | 2018-02-07 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108492260A true CN108492260A (zh) | 2018-09-04 |
| CN108492260B CN108492260B (zh) | 2019-01-08 |
Family
ID=63344647
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810121330.4A Active CN108492260B (zh) | 2018-02-07 | 2018-02-07 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US11430128B2 (zh) |
| CN (1) | CN108492260B (zh) |
| WO (1) | WO2019153746A1 (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109657717A (zh) * | 2018-12-17 | 2019-04-19 | 罗世彬 | 一种基于多尺度密集结构特征提取的异源图像匹配方法 |
| CN109902557A (zh) * | 2019-01-14 | 2019-06-18 | 云南大学 | 蚀变特征矿物的蚀变异常信息的提取方法 |
| WO2019153746A1 (zh) * | 2018-02-07 | 2019-08-15 | 长安大学 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
| CN112652010A (zh) * | 2020-12-25 | 2021-04-13 | 安徽工业大学 | 一种基于图像增强和霍夫变换的球团粒径检测方法 |
| CN116228599A (zh) * | 2023-05-08 | 2023-06-06 | 山东省地质测绘院 | 基于遥感影像增强的矿产资源动态监测方法 |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113269056B (zh) * | 2021-05-06 | 2023-11-07 | 四川省亚通工程咨询有限公司 | 基于卫星影像的公路勘察设计方法 |
| CN113408462B (zh) * | 2021-06-29 | 2023-05-02 | 西南交通大学 | 基于卷积神经网络与类别热力图的滑坡遥感信息提取方法 |
| CN115526924B (zh) * | 2022-11-29 | 2023-06-23 | 湖南长理尚洋科技有限公司 | 一种蒙特卡洛模拟的水文环境建模方法及系统 |
| CN115661669B (zh) * | 2022-12-13 | 2023-03-21 | 山东省土地发展集团有限公司 | 一种基于视频监控的耕地违法占用监测方法及系统 |
| CN115797502B (zh) * | 2023-02-08 | 2023-04-25 | 航天宏图信息技术股份有限公司 | 一种基于dem数据的地性线提取方法和装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156975A (zh) * | 2011-04-22 | 2011-08-17 | 西安电子科技大学 | 基于支撑值变换和多尺度冗余字典学习的自然图像去噪方法 |
| KR101580594B1 (ko) * | 2015-07-06 | 2015-12-28 | 국방과학연구소 | 정규격자 기반 허프변환을 이용한 고해상도 위성영상의 선형 추출 방법 |
| CN106295604A (zh) * | 2016-08-19 | 2017-01-04 | 厦门大学 | 基于联合滤波的遥感影像路网提取技术 |
| CN107301661A (zh) * | 2017-07-10 | 2017-10-27 | 中国科学院遥感与数字地球研究所 | 基于边缘点特征的高分辨率遥感图像配准方法 |
| CN107610103A (zh) * | 2017-08-30 | 2018-01-19 | 西安石油大学 | 一种基于x射线成像的钢管螺旋弧焊缝区域提取方法 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6912322B2 (en) * | 2000-12-14 | 2005-06-28 | Itt Manufacturing Enterprises Inc. | Adaptive process for removing streaks in multi-band digital images |
| US20130066666A1 (en) * | 2010-01-22 | 2013-03-14 | Monsanto Technology Llc | Enhancing Performance of Crops Within An Area of Interest |
| CN103488968B (zh) * | 2012-06-14 | 2016-12-21 | 株式会社日立制作所 | 遥感图像的混合像素物质构成精细化分解装置及方法 |
| WO2014183259A1 (zh) * | 2013-05-14 | 2014-11-20 | 中国科学院自动化研究所 | 一种全色−多光谱遥感图像融合方法 |
| ES2530687B1 (es) * | 2013-09-04 | 2016-08-19 | Shot & Shop. S.L. | Método implementado por ordenador para recuperación de imágenes por contenido y programa de ordenador del mismo |
| JP6660313B2 (ja) * | 2014-06-30 | 2020-03-11 | ベンタナ メディカル システムズ， インコーポレイテッド | 画像解析を用いた核のエッジの検出 |
| US10303916B2 (en) * | 2015-07-31 | 2019-05-28 | Fotonation Limited | Image processing apparatus |
| US10373339B2 (en) * | 2015-12-14 | 2019-08-06 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Hyperspectral scene analysis via structure from motion |
| US10217236B2 (en) * | 2016-04-08 | 2019-02-26 | Orbital Insight, Inc. | Remote determination of containers in geographical region |
| CN106294705B (zh) * | 2016-08-08 | 2017-12-15 | 长安大学 | 一种批量遥感影像预处理方法 |
| US10154624B2 (en) * | 2016-08-08 | 2018-12-18 | The Climate Corporation | Estimating nitrogen content using hyperspectral and multispectral images |
| CN107563964B (zh) * | 2017-08-22 | 2020-09-04 | 长光卫星技术有限公司 | 大面阵亚米级夜景遥感影像的快速拼接方法 |
| CN108492260B (zh) * | 2018-02-07 | 2019-01-08 | 长安大学 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
-
2018
- 2018-02-07 CN CN201810121330.4A patent/CN108492260B/zh active Active
- 2018-09-17 WO PCT/CN2018/105966 patent/WO2019153746A1/zh active Application Filing
- 2018-09-17 US US16/627,900 patent/US11430128B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156975A (zh) * | 2011-04-22 | 2011-08-17 | 西安电子科技大学 | 基于支撑值变换和多尺度冗余字典学习的自然图像去噪方法 |
| KR101580594B1 (ko) * | 2015-07-06 | 2015-12-28 | 국방과학연구소 | 정규격자 기반 허프변환을 이용한 고해상도 위성영상의 선형 추출 방법 |
| CN106295604A (zh) * | 2016-08-19 | 2017-01-04 | 厦门大学 | 基于联合滤波的遥感影像路网提取技术 |
| CN107301661A (zh) * | 2017-07-10 | 2017-10-27 | 中国科学院遥感与数字地球研究所 | 基于边缘点特征的高分辨率遥感图像配准方法 |
| CN107610103A (zh) * | 2017-08-30 | 2018-01-19 | 西安石油大学 | 一种基于x射线成像的钢管螺旋弧焊缝区域提取方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘志恒: "一种资源三号卫星影像的云量判读方法与应用", 《遥感信息》 * |
| 高景昌: "遥感图像线性构造信息微机机助提取", 《遥感学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019153746A1 (zh) * | 2018-02-07 | 2019-08-15 | 长安大学 | 基于张量投票耦合霍夫变换的地质线性体提取方法 |
| US11430128B2 (en) | 2018-02-07 | 2022-08-30 | Chang'an University | Geological linear body extraction method based on tensor voting coupled with Hough transformation |
| CN109657717A (zh) * | 2018-12-17 | 2019-04-19 | 罗世彬 | 一种基于多尺度密集结构特征提取的异源图像匹配方法 |
| CN109902557A (zh) * | 2019-01-14 | 2019-06-18 | 云南大学 | 蚀变特征矿物的蚀变异常信息的提取方法 |
| CN112652010A (zh) * | 2020-12-25 | 2021-04-13 | 安徽工业大学 | 一种基于图像增强和霍夫变换的球团粒径检测方法 |
| CN116228599A (zh) * | 2023-05-08 | 2023-06-06 | 山东省地质测绘院 | 基于遥感影像增强的矿产资源动态监测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20210287376A1 (en) | 2021-09-16 |
| WO2019153746A1 (zh) | 2019-08-15 |
| US11430128B2 (en) | 2022-08-30 |
| US20220164960A9 (en) | 2022-05-26 |
| CN108492260B (zh) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108492260B (zh) | 基于张量投票耦合霍夫变换的地质线性体提取方法 | |
| Han et al. | Extraction and analysis of geological lineaments combining a DEM and remote sensing images from the northern Baoji loess area | |
| CN102609726B (zh) | 利用面向对象技术融合高空间和高时间分辨率数据的遥感图像分类方法 | |
| AU2020101520A4 (en) | Method for identifying exposed and buried fault structures in granite-type uranium province | |
| Wang et al. | Estimation of soil salt content using machine learning techniques based on remote-sensing fractional derivatives, a case study in the Ebinur Lake Wetland National Nature Reserve, Northwest China | |
| CN107765323A (zh) | 一种基于aster卫星数据的成矿预测方法 | |
| CN108171210B (zh) | 一种植被覆盖区蚀变遥感异常信息的提取方法及系统 | |
| Lasaponara et al. | Image enhancement, feature extraction and geospatial analysis in an archaeological perspective | |
| CN109872389A (zh) | 一种基于三维地表模型的遥感地质构造解译方法 | |
| CN110008908A (zh) | 一种基于高分遥感影像的草原围栏提取方法 | |
| CN105513060A (zh) | 一种视觉感知启发的高分辨率遥感图像分割方法 | |
| Ni et al. | Semi-automatic extraction and mapping of dyke swarms based on multi-resolution remote sensing images: Applied to the dykes in the Kuluketage region in the northeastern Tarim Block | |
| CN114386497A (zh) | 面向铀成矿构造的航空高光谱和伽玛能谱数据融合方法 | |
| Pal et al. | Utilization of Landsat ETM+ data for mineral-occurrences mapping over Dalma and Dhanjori, Jharkhand, India: an Advanced Spectral Analysis approach | |
| Mohammadpour et al. | Automatic lineament extraction method in mineral exploration using CANNY algorithm and hough transform | |
| CN114463623A (zh) | 一种基于多尺度遥感影像的耕地变化检测方法及装置 | |
| Salehi et al. | Mapping ultramafic complexes using airborne imaging spectroscopy and spaceborne data in Arctic regions with abundant lichen cover, a case study from the Niaqornarssuit complex in South West Greenland | |
| Chen et al. | Improving the accuracy of hydrothermal alteration mapping based on image fusion of ASTER and Sentinel-2A data: a case study of Pulang Cu deposit, Southwest China | |
| Bektas Balcik et al. | Determination of magnitude and direction of land use/land cover changes in Terkos Water Basin, Istanbul | |
| Brandmeier | Remote sensing of Carhuarazo volcanic complex using ASTER imagery in Southern Peru to detect alteration zones and volcanic structures–a combined approach of image processing in ENVI and ArcGIS/ArcScene | |
| Sui et al. | Processing of multitemporal data and change detection | |
| Hashim et al. | Land use land cover analysis with pixel-based classification approach | |
| Nanare et al. | Remote sensing satellite image analysis for deforestation in yavatmal district, maharashtra, india | |
| Pisek et al. | Assessment of a global leaf area index product from SPOT-4 VEGETATION data over selected sites in Canada | |
| Li et al. | Multisource remotely sensed wetland information exploration using interactive visualization methods |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20180904 Assignee: Xi'an tianbang surveying and Mapping Technology Co.,Ltd. Assignor: CHANG'AN University Contract record no.: X2022980004777 Denomination of invention: Geological linear volume extraction method based on tensor voting coupled Hough transform Granted publication date: 20190108 License type: Common License Record date: 20220425 |