[发明专利]一种自适应探地雷达数据去噪方法

权利要求书

1.一种自适应探地雷达数据去噪方法，包括如下步骤：

(1)对于原始探地雷达数据，截取其前k纳秒内的数据并确定数据的标准道，k为大于0的实数；

(2)对于任一扫描道，计算该扫描道与标准道在前k纳秒内的互相关函数并取其最大值；

(3)根据步骤(2)遍历所有扫描道，得到由各扫描道互相关函数最大值组成的初至曲线L，进而计算出初至曲线L的均方根高D；

(4)根据以下公式计算出与相干噪声对应的奇异值个数N；

N＝ROUND{0.2634D+1.3086}

其中：ROUND{}为向上取整函数；

(5)对原始探地雷达数据进行奇异值分解，剔除奇异值分解后的前N个数据成分，将剩余数据成分重新合成后即完成去噪过程。

2.根据权利要求1所述的自适应探地雷达数据去噪方法，其特征在于：所述步骤(1)中的原始探地雷达数据为n×m大小的数据矩阵，n和m分别对应为时间和空间上的采样点数；所述扫描道即对应探测对象表面上任一采样点的雷达数据，所述标准道即对应探测对象表面上某一采样点的雷达数据，该采样点为随机或人为指定。

3.根据权利要求2所述的自适应探地雷达数据去噪方法，其特征在于：所述步骤(2)中的互相关函数表达式如下：

Rj(τ)=Σt=-∞∞x0(t)xj*(t-τ)]]>

其中：R_j(τ)为第j扫描道与标准道在前k纳秒内的互相关函数，x₀(t)为t时刻的标准道，t和τ均为时间采样序号且1≤τ≤K，K为前k纳秒内雷达数据的采样点数，x_j^*(t-τ)为x_j(t-τ)的共轭复数，x_j(t-τ)为t-τ时刻的第j扫描道，j为自然数且1≤j≤m。

4.根据权利要求2所述的自适应探地雷达数据去噪方法，其特征在于：所述步骤(3)中根据以下公式计算初至曲线L的均方根高D：

D=Σj=1m(d(j)-d‾)2m]]>

其中：d(j)为第j扫描道与标准道在前k纳秒内的互相关函数最大值，为初至曲线L的平均值。

5.根据权利要求1所述的自适应探地雷达数据去噪方法，其特征在于：k＝5。