CN110414060A - 一种基于四阶谱矩的位场边界识别方法 - Google Patents
一种基于四阶谱矩的位场边界识别方法 Download PDFInfo
- Publication number
- CN110414060A CN110414060A CN201910579203.3A CN201910579203A CN110414060A CN 110414060 A CN110414060 A CN 110414060A CN 201910579203 A CN201910579203 A CN 201910579203A CN 110414060 A CN110414060 A CN 110414060A
- Authority
- CN
- China
- Prior art keywords
- quadravalence
- potential field
- point
- data point
- matrix
- 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 53
- 230000003595 spectral effect Effects 0.000 title claims abstract description 42
- 238000007373 indentation Methods 0.000 claims abstract description 22
- 238000000605 extraction Methods 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 11
- 102000002274 Matrix Metalloproteinases Human genes 0.000 claims description 9
- 108010000684 Matrix Metalloproteinases Proteins 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 238000007796 conventional method Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000005415 magnetization Effects 0.000 description 2
- 241001046947 Ectropis obliqua Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于四阶谱矩的位场边界识别方法，首先通过滑动窗口计算出位场每一个数据点的四阶谱矩元，然后根据本发明定义的四阶曲率弧刻痕系数公式计算其特征提取结果，最后由该系数的高值指示出场源的边界位置。本发明提出的基于四阶谱矩的位场边界识别方法，是一种适用于提取弧形边界的方法，克服了常规传统方法和表面二阶谱矩不能有效识别弧形边界的缺陷。重磁异常的仿真试验结果表明，本发明提出的基于四阶谱矩的位场边界识别方法能够均衡收敛地识别弧形边界，且对深度的敏感性较低，可用于提取位场中的弱曲率弧刻痕信号。
Description
技术领域
本发明属于地球物理学位场数据处理技术领域，具体涉及一种基于四阶谱矩的位场边界识别方法。
背景技术
揭示地下地质体的边缘位置是位场数据解释中的重要内容。各地质体边缘是具有一定密度或磁性差异的，在地质体的边缘附近，重磁异常的变化率较大，因此可以利用这一特点对其进行边界识别。目前，存在三大类位场源边界识别方法，它们分别是数值类计算方法、数理统计类计算方法及其它计算方法。其中，数值类计算方法是现阶段研究最多的一类方法，且其应用也最为广泛，数理统计类计算方法次之，其他计算方法应用较少。
数值类计算方法大多基于位场异常不同方向、不同阶次的导数及其组合，如总水平导数(TotalHorizontal Derivative,THDR)、垂向导数(Vertical Derivative,VDR)、解析信号振幅(Analytical Signal Amplitude,ASM)、倾斜角(Tilt Angle,TA)和Theta map。数理统计类方法主要包括两种，小子域滤波和归一化标准差(Normalized StandardDeviation,NSTD)，该类方法通过选取一定尺寸大小的滑动窗口来进行计算。上述方法均属于传统的边界识别方法，其中，目前比较常用的方法有以下几种：
总水平导数法(THDR)是一种应用广泛的位场边界识别方法，计算公式为
其中，u表示位场异常。
解析信号振幅(ASM)是由Nabighian提出的，该方法也叫总梯度模量法，由三个方向上的一阶导数组合得到，表达式为
式中，THDR代表总水平导数，而VDR则表示垂向导数，u表示位场异常；ASM计算结果的极大值表示地质体的位置，零值刻画的是地质体边界信息。
倾斜角法(TA)，即斜导数法，是一种具有均衡意义的场源边界识别方法，它在1994年由Miller和Singh提出。倾斜角法实质上是各一阶导数组合的比值，具体计算公式如下，
总水平导数的倾斜角(TAHG)法是Ferreira于2013年提出的，在总水平导数的结果上计算倾斜角，以达到边界增强的效果。其计算式子如下
Theta map方法是由Wijns等于2005年首次提出，它实质上是THDR与ASM的比值，表达式如下，
归一化标准差(NSTD)方法属于数理统计类边界识别方法，由Cooper和Cowan于2008年提出，它是垂向坐标方向的一阶导数的标准差与三个坐标方向一阶导数标准差的比值，计算式子如下
其中，u代表位场异常，σ表示标准差。这种方法是利用一定大小的滑动窗口来计算窗口中心点的归一化标准差值，其计算结果的极大值位置代表地质体的边界位置。
上述常规的边界识别方法可以识别位场的线形边界，但对于弧形边界提取效果较差。谱矩方法以随机过程理论为基础，早先被广泛应用于工程学的表面形貌识别，因其能够细致地描述表面的几何特征而被引入地球物理学领域，用于位场表面数据的解释。目前，孙艳云等基于表面二阶谱矩，定义了位场的脊形化系数Λ2，计算公式如下
式中，Δ2和M2均为二阶谱矩的统计不变量，其表达式如下
M2＝m20+m02 (9)
m20、m11、m02是二阶谱矩元，分别为x和y方向上一阶导数的方差或协方差。由于二阶谱矩表征位场的斜率信息，故该系数主要提取的是地质体的线形边界，而不能有效识别弧形边界。
发明内容
针对传统方法以及表面二阶谱矩识别场源边界的不足，考虑表面四阶谱矩在重磁异常数据处理中的应用。由于四阶谱矩主要描述位场表面与曲率相关的几何特征，因此可利用其提取弧形的边界信息。本发明提出一种基于四阶谱矩的边界识别方法，理论试验结果表明，该方法可以进一步完善弧形边界信息的识别，同时对深度的敏感性较低，可用于提取位场中的弱曲率弧刻痕信号。
本发明所采用的技术方案是：一种基于四阶谱矩的位场边界识别方法，其特征在于，包括以下步骤：
步骤1：在x、y方向上分别以间隔Δx、Δy等间距采样，采样点数分别为M和N，得到研究区域的位场异常数据矩阵A；
步骤2：采用滑动窗口圈定窗口中心点及其邻近数据点的异常值，计算窗口面元的四阶谱矩，将结果定位于面元中心点；
步骤3：逐点移动滑动窗口，得到矩阵A上每个数据点对应的四阶谱矩元m40、m31、m22、m13、m04的值；
步骤4：根据矩阵A中每一个数据点处四阶谱矩元的值计算该点处统计不变量M4、Δ4的值：
M4＝m40+2m22+m04
步骤5：计算矩阵A中每一个数据点处的四阶曲率弧刻痕系数值：
得到矩阵B，即为四阶曲率弧刻痕系数对位场数据的提取结果，其高值指示场源的边界位置。
作为优选，步骤1中，对研究区域的位场异常进行测量，在x、y方向上分别以间隔Δx、Δy等间距采样，采样点数分别为M和N，得到N×M的位场异常数据矩阵A；选取尺寸为N1×M1的滑动窗口，沿x方向包含M1个数据点，沿y方向包含N1个数据点，M1、N1均为奇数。
作为优选，步骤2中，将N×M的矩阵A的列向外延拓(M1-1)/2列，行向外延拓(N1-1)/2行，得到大小为[N+(N1-1)]×[M+(M1-1)]的新矩阵对矩阵中每个数据点求二阶导数，得到在矩阵中，选定延拓前第一行第一列的数据点作为滑动窗口的中心点，利用窗口圈定窗口中心点及其邻近数据点的异常值，形成大小为N1×M1的局部面元，通过如下公式计算其四阶谱矩元的值：
其中，zxx(xi,yk)、zxy(xi,yk)、zyy(xi,yk)分别为数据点(xi,yk)处位场异常的二阶导数；将上述各谱矩元的值赋给窗口中心的数据点。
作为优选，步骤3中，在矩阵中，按延拓前的数据点逐点移动滑动窗口，计算出每一个数据点处m40、m31、m22、m13、m04的值。
相对于现有技术，本发明的有益效果是：本发明提出的基于四阶谱矩的位场边界识别方法能够均衡、收敛地识别弧形边界，且不产生干扰刻痕。同时，场源埋深的深浅对提取结果影响较小，即该方法对深度的敏感性较低，可用于提取位场中的弱曲率弧刻痕信号。
附图说明
图1为本发明实施例的流程图；
图2为本发明实施例的球体模型平面图；
图3为本发明实施例的不同埋深下球体模型的理论重力异常和四阶曲率弧刻痕系数Λ4的提取结果：(a)埋深10km的理论重力异常；(b)埋深10km时Λ4的提取结果；(c)埋深12km的理论重力异常；(d)埋深12km时Λ4的提取结果；(e)埋深14km的理论重力异常；(f)埋深14km时Λ4的提取结果；
图4为本发明实施例的常规边界识别方法对埋深10km的球体模型理论重力异常的提取结果：(a)总水平导数(THDR)；(b)总水平导数的倾斜角(TAHG)；(c)Theta map；(d)归一化标准差(NSTD)；
图5为本发明实施例的不同埋深下球体模型的理论磁异常和四阶曲率弧刻痕系数Λ4的提取结果：(a)埋深10km的理论磁异常；(b)埋深10km时Λ4的提取结果；(c)埋深12km的理论磁异常；(d)埋深12km时Λ4的提取结果；(e)埋深14km的理论磁异常；(f)埋深14km时Λ4的提取结果；
图6为本发明实施例的常规边界识别方法对埋深10km的球体模型理论磁异常的提取结果：(a)总水平导数(THDR)；(b)总水平导数的倾斜角(TAHG)；(c)Theta map；(d)归一化标准差(NSTD)。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。
请见图1，本发明提供的一种基于四阶谱矩的位场边界识别方法，包括以下步骤：
步骤1：对研究区域的位场异常进行测量，在x、y方向上分别以间隔Δx、Δy等间距采样，采样点数分别为M和N，得到N×M的位场异常数据矩阵A；选取尺寸为N1×M1的滑动窗口提取网格点处的位场异常值，即沿x方向包含M1个数据点，沿y方向包含N1个数据点，M1、N1均为奇数；
步骤2：将N×M的矩阵A的列向外延拓(M1-1)/2列，行向外延拓(N1-1)/2行，得到大小为[N+(N1-1)]×[M+(M1-1)]的新矩阵对矩阵中每个数据点求二阶导数，得到在矩阵中，选定延拓前第一行第一列的数据点作为滑动窗口的中心点，利用窗口圈定中心点及其邻近数据点的异常值，形成大小为N1×M1的局部面元。令u(xi,yk)为滑动窗口中数据点(xi,yk)处的异常值，通过如下公式计算其四阶谱矩元的值：
式中，uxx(xi,yk)、uxy(xi,yk)、uyy(xi,yk)分别为数据点(xi,yk)处位场异常的二阶导数；将上述各谱矩元的值赋给窗口中心的数据点；
步骤3：在矩阵中，按延拓前的数据点逐点移动滑动窗口，计算出每一个数据点处m40、m31、m22、m13、m04的值；
步骤4：根据矩阵A中每一个数据点处四阶谱矩元的值计算该点处统计不变量M4、Δ4的值：
M4＝m40+2m22+m04
步骤5：计算矩阵A中每一个数据点处的四阶曲率弧刻痕系数值：
得到矩阵B，即为四阶曲率弧刻痕系数对位场数据的提取结果，其高值指示场源的边界位置。
本发明提出的基于四阶谱矩的位场边界识别方法，是一种适用于提取弧形边界的方法，克服了常规传统方法和表面二阶谱矩不能有效识别弧形边界的缺陷。重磁异常的仿真试验结果表明，本发明提出的基于四阶谱矩的位场边界识别方法能够均衡收敛地识别弧形边界，且对深度的敏感性较低，可用于提取位场中的弱曲率弧刻痕信号。
为了验证是发明提供的一种基于四阶谱矩的边界识别方法，现将该方法应用于仿真试验数据。由于球体在地面的位场异常具有完整的曲率弧刻痕特征，是应用四阶谱矩来提取位场中曲率弧刻痕信息的理想模型，故设计平面图如图2所示的球体仿真模型，两球体半径均为8km，中心埋深均为10km，球心坐标为(-20km，-20km)的球体剩余密度为0.8g/cm3，磁化强度为0.8A/m，球心坐标为(20km，20km)的球体剩余密度0.2g/cm3，磁化强度为0.3A/m，x轴与y轴方向点距均为0.5km。
在理论试验部分，首先测试本发明提出的四阶曲率弧刻痕系数对深度的敏感性，然后与当前一些常规的位场边界识别方法进行比较，包括总水平导数(THDR)、总水平导数的倾斜角(TAHG)、Thetamap法和归一化标准差(NSTD)，以验证本发明提出的边界识别方法的有效性。以下给出两组试验分别检验四阶曲率弧刻痕系数在重、磁异常数据中的边界识别效果。
重、磁异常仿真试验中选取的滑动窗口大小分别为5×5和3×3。图3和图5分别为该球体模型在中心埋深10km、12km、14km时的理论重、磁异常及对应四阶曲率弧刻痕系数Λ4的提取结果。分析其提取结果可知，Λ4提取的曲率弧分布信息受场源强度影响较小，能稳定、有效地提取不同深度地质体的弧形边界信息，对深度的敏感性较低。图4和图6分别为几种常规边界识别方法对重、磁异常的提取结果，可见四种传统方法都能提取出两个球状场源的边界，但均存在边界发散的问题，除总水平导数外的另外三种方法的提取结果中存在干扰刻痕，这些缺点使得实际应用中难以准确分辨出地质体的边界位置。对比图3(b)和图5(b)易知，四阶曲率弧刻痕系数Λ4相较于传统方法能更准确、收敛地刻画球状重磁异常体的边界，且没有干扰刻痕。从实验结果可以看出，本实施例的边界识别方法在提取场源边界时明显优于常规的传统方法。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。
Claims (4)
1.一种基于四阶谱矩的位场边界识别方法，其特征在于，包括以下步骤：
步骤1：在x、y方向上分别以间隔Δx、Δy等间距采样，采样点数分别为M和N，得到研究区域的位场异常数据矩阵A；
步骤2：采用滑动窗口圈定窗口中心点及其邻近数据点的异常值，计算窗口面元的四阶谱矩，将结果定位于面元中心点；
步骤3：逐点移动滑动窗口，得到矩阵A上每个数据点对应的四阶谱矩元m40、m31、m22、m13、m04的值；
步骤4：根据矩阵A中每一个数据点处四阶谱矩元的值计算该点处统计不变量M4、Δ4的值：
M4＝m40+2m22+m04
步骤5：计算矩阵A中每一个数据点处的四阶曲率弧刻痕系数值：
得到矩阵B，即为四阶曲率弧刻痕系数对位场数据的提取结果，其高值指示场源的边界位置。
2.根据权利要求1所述的基于四阶谱矩的位场边界识别方法，其特征在于：步骤1中，对研究区域的位场异常进行测量，在x、y方向上分别以间隔Δx、Δy等间距采样，采样点数分别为M和N，得到N×M的位场异常数据矩阵A；选取尺寸为N1×M1的滑动窗口，沿x方向包含M1个数据点，沿y方向包含N1个数据点，M1、N1均为奇数。
3.根据权利要求1所述的基于四阶谱矩的位场边界识别方法，其特征在于：步骤2中，将N×M的矩阵A的列向外延拓(M1-1)/2列，行向外延拓(N1-1)/2行，得到大小为[N+(N1-1)]×[M+(M1-1)]的新矩阵对矩阵中每个数据点求二阶导数，得到在矩阵中，选定延拓前第一行第一列的数据点作为滑动窗口的中心点，利用窗口圈定窗口中心点及其邻近数据点的异常值，形成大小为N1×M1的局部面元，通过如下公式计算其四阶谱矩元的值：
其中，zxx(xi,yk)、zxy(xi,yk)、zyy(xi,yk)分别为数据点(xi,yk)处位场异常的二阶导数；将上述各谱矩元的值赋给窗口中心的数据点。
4.根据权利要求2所述的基于四阶谱矩的位场边界识别方法，其特征在于：步骤3中，在矩阵中，按延拓前的数据点逐点移动滑动窗口，计算出每一个数据点处m40、m31、m22、m13、m04的值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910579203.3A CN110414060B (zh) | 2019-06-28 | 2019-06-28 | 一种基于四阶谱矩的位场边界识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910579203.3A CN110414060B (zh) | 2019-06-28 | 2019-06-28 | 一种基于四阶谱矩的位场边界识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110414060A true CN110414060A (zh) | 2019-11-05 |
| CN110414060B CN110414060B (zh) | 2023-01-03 |
Family
ID=68358873
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910579203.3A Expired - Fee Related CN110414060B (zh) | 2019-06-28 | 2019-06-28 | 一种基于四阶谱矩的位场边界识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110414060B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111580174A (zh) * | 2020-05-29 | 2020-08-25 | 中国地质科学院 | 一种基于帕德近似的重磁数据向下延拓方法 |
| CN112305623A (zh) * | 2020-08-20 | 2021-02-02 | 中国地质科学院地球物理地球化学勘查研究所 | 一种基于谱融合的位场特征获取方法和装置 |
| CN113886753A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Tilt法航磁边界检测方法、装置 |
| CN116609844A (zh) * | 2023-03-13 | 2023-08-18 | 中国自然资源航空物探遥感中心 | 一种用于提取地球物理三维断裂信息的方法和计算设备 |
| CN116774303A (zh) * | 2023-06-26 | 2023-09-19 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101144860A (zh) * | 2007-10-16 | 2008-03-19 | 哈尔滨工业大学 | 基于选择性核主成份分析的高光谱图像异常点的检测方法 |
| CA2978500A1 (en) * | 2015-03-04 | 2016-09-09 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field data |
| CN109407161A (zh) * | 2018-09-21 | 2019-03-01 | 中国国土资源航空物探遥感中心 | 用于提取地球物理磁异常场边界的磁场刻痕分析方法 |
-
2019
- 2019-06-28 CN CN201910579203.3A patent/CN110414060B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101144860A (zh) * | 2007-10-16 | 2008-03-19 | 哈尔滨工业大学 | 基于选择性核主成份分析的高光谱图像异常点的检测方法 |
| CA2978500A1 (en) * | 2015-03-04 | 2016-09-09 | Institute of Mineral Resources, Chinese Academy of Geological Sciences | Method for automatically extracting tectonic framework of potential field data |
| CN109407161A (zh) * | 2018-09-21 | 2019-03-01 | 中国国土资源航空物探遥感中心 | 用于提取地球物理磁异常场边界的磁场刻痕分析方法 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111580174A (zh) * | 2020-05-29 | 2020-08-25 | 中国地质科学院 | 一种基于帕德近似的重磁数据向下延拓方法 |
| CN111580174B (zh) * | 2020-05-29 | 2023-08-11 | 中国地质科学院 | 一种基于帕德近似的重磁数据向下延拓方法 |
| CN112305623A (zh) * | 2020-08-20 | 2021-02-02 | 中国地质科学院地球物理地球化学勘查研究所 | 一种基于谱融合的位场特征获取方法和装置 |
| CN113886753A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Tilt法航磁边界检测方法、装置 |
| CN113886753B (zh) * | 2021-10-09 | 2022-04-12 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Tilt法航磁边界检测方法、装置 |
| CN116609844A (zh) * | 2023-03-13 | 2023-08-18 | 中国自然资源航空物探遥感中心 | 一种用于提取地球物理三维断裂信息的方法和计算设备 |
| CN116609844B (zh) * | 2023-03-13 | 2024-03-12 | 中国自然资源航空物探遥感中心 | 一种用于提取地球物理三维断裂信息的方法和计算设备 |
| CN116774303A (zh) * | 2023-06-26 | 2023-09-19 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
| CN116774303B (zh) * | 2023-06-26 | 2024-05-07 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110414060B (zh) | 2023-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110414060A (zh) | 一种基于四阶谱矩的位场边界识别方法 | |
| CN106289777B (zh) | 一种基于几何度量的多工况滚动轴承性能评估方法 | |
| CN104867126B (zh) | 基于点对约束和三角形网的有变化区域的合成孔径雷达图像配准方法 | |
| CN110007350B (zh) | 一种磁探测方法盲区的分析方法 | |
| CN106651942A (zh) | 基于特征点的三维旋转运动检测与旋转轴定位方法 | |
| CN108010116A (zh) | 点云特征点检测方法和点云特征提取方法 | |
| CN106023298A (zh) | 基于局部泊松曲面重建的点云刚性配准方法 | |
| CN104299260A (zh) | 一种基于sift和lbp的点云配准的接触网三维重建方法 | |
| CN103852795A (zh) | 一种水下小目标的磁异常信号提取方法 | |
| CN106500594B (zh) | 融合反射强度和几何特征的铁路轨道半自动检测方法 | |
| CN108664983A (zh) | 一种尺度与特征强度自适应的surf特征点匹配方法 | |
| CN108830899A (zh) | 一种基于几何特征的快速自适应角点结合的圆心检测方法 | |
| CN105469402B (zh) | 基于空间形状上下文特征的汽车零件识别方法 | |
| CN109579827A (zh) | 一种基于弧形阵列的磁性目标探测和定位方法 | |
| CN107632319B (zh) | 基于gst的碳酸盐岩非均质储层中溶洞识别标定方法 | |
| Liu et al. | The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies | |
| CN106127147B (zh) | 一种基于三维数据的人脸深度纹理修复方法 | |
| Kruglov | The algorithm of the roundwood volume measurement via photogrammetry | |
| CN110688767A (zh) | 一种评价岩体裂隙网络模型精度的综合差异度的方法 | |
| CN110458111A (zh) | 基于LightGBM的车载激光点云电力线的快速提取方法 | |
| CN115588178A (zh) | 一种高精地图要素自动化提取的方法 | |
| CN112946760B (zh) | 基于正则化方法的未爆弹三维立体成像方法、装置及*** | |
| CN104766085B (zh) | 一种多尺度图形识别方法 | |
| Liu et al. | Detection of arbitrary triangle | |
| CN112505752B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230103 |