[发明专利]一种探地雷达层位自动追踪方法及系统

权利要求书

1.一种探地雷达层位自动追踪方法，包括：

从探地雷达二维剖面图中确定前一道种子点在相应道的第一位置和在预设的相关系数计算范围内的第二位置；

基于所述第一位置确定当前道的映射点，并在预设的滑动范围的限制下，将以所述映射点为中心且以所述计算范围为半径所对应的计算区域分别向上、下滑动，确定每个计算区域与所述前一道种子点在所述计算范围内对应的区域间的相关系数，得到相应的系数集合；

计算所有相关系数的方差，判断集合内的最大相关系数是否超过预设的系数阈值并且所述方差是否超过预设的相关系数方差阈值，若满足，则将所述最大相关系数对应的计算区域作为最优相关区域，进入下一步；

基于所述第二位置确定所述最优相关区域内的同位置数据点，基于此，利用预设的搜索范围，确定当前道的种子点，以追踪下一道层位。

2.根据权利要求1所述的方法，其特征在于，所述方法，还包括：按照预设的特征点检验条件，对所述当前道的种子点进行有效性验证，其中，

判断所述当前道的种子点的层位深度数据与所述前一道种子点的层位深度数据的差值是否小于预设的深度差值阈值，若满足，则进入下一步；

判断所述当前道的种子点的层位深度数据是否处于预设的深度范围阈值之内，若处于，则进入下一步；

确定以所述当前道的种子点为中心点且以所述搜索范围为半径所对应的区域内的最大振幅极值点和最小振幅极值点，并计算两个极值点的第一极值深度差和第一极值振幅差；

判断所述第一极值深度差是否大于预设的极值深度差阈值，并且判断所述第一极值振幅差是否处于极值幅度差范围阈值之内，若满足，则确定当前种子点有效。

3.根据权利要求1或2所述的方法，其特征在于，所述特征点的类型选自谷值、峰值和零交叉点中的一种。

4.根据权利要求1或2所述的方法，其特征在于，所述方法还包括：对完成层位追踪的所有道的层位线进行层位置信度评价，其中，

按照预设间距，将所有道分成若干组；

计算每组的层位特征指标；

对所述每组的层位特征指标进行加权打分处理，得到关于层位特征的置信度。

5.根据权利要求4所述的方法，其特征在于，在计算每组的层位特征指标步骤中，包括：

计算每组中的不连续道所占比例；和/或

按照预设的第一数量，以同组内每道种子点为中心，计算每道种子点向上所述第一数量个数据点的探地雷达数据、以及向下所述第一数量个数据点的探地雷达数据的第一方差，并计算同组内所有道的所述第一方差的均值；和/或

按照预设的第二数量，以同组内每道种子点为中心，计算每道种子点向上所述第二数量个数据点的探地雷达数据、以及向下所述第二数量个数据点的探地雷达数据的能量平均，并计算同组内所有道的所述能量平均的均值。

6.根据权利要求4所述的方法，其特征在于，所述方法还包括：将属于不同置信度评价区间的组所对应的层位线，按照预设的展示类型，以不同形式显示在所述二维剖面图上，其中，所述展示类型选自线宽、颜色和线型差异中的一种。

7.根据权利要求1或2所述的方法，其特征在于，所述方法，还包括：从探地雷达二维剖面图中确定前一道的种子点，其中，

对所述二维剖面图进行去干扰预处理，所述预处理包括背景去噪、滤波和增益处理；

从经过预处理的所述二维剖面图中，选取前一道的初始特征点；

以所述初始特征点为中心在所述搜索范围内，选取相应的种子点。

8.根据权利要求1或2所述的方法，其特征在于，在计算所有相关系数的方差，判断集合内的最大相关系数是否超过预设的系数阈值、并且所述方差是否超过预设的相关系数方差阈值步骤中，若不满足任一标准，则停止追踪，并标记所述当前道与前一道为不连续道。