[发明专利]一种浅层对地雷达信号去除地面波的方法

权利要求书

1.一种浅层对地雷达信号去除地面波的方法，其特征在于，所述地雷达信号由n个单通道雷达波形组成，并记为Ascan；所述n为大于0的自然数；所述浅层对地雷达信号去除地面波的方法，包括以下步骤：

对连续的雷达波形按等时窗进行采样，任一等时窗内的任一雷达信号Ascan含有m个采集点；所述m为大于0的自然数；

所述单通道雷达在运行时产生n个雷达信号Ascan，并标记为通道雷达波bscan；

对通道雷达波bscan的任一雷达信号Ascan采用墨西哥帽小波进行连续小波变换，并寻找通道雷达波bscan中的任一雷达信号Ascan与墨西哥帽波形最接近的位置k；

将任一雷达信号Ascan的位置k顺序排列得到位置集合序列K；

对位置集合序列K进行平滑处理，并移除地面波。

2.根据权利要求1所述的一种浅层对地雷达信号去除地面波的方法，其特征在于，对通道雷达波bscan的任一雷达信号Ascan采用墨西哥帽小波进行连续小波变换，并寻找通道雷达波bscan中的任一雷达信号Ascan与墨西哥帽波形最接近的位置k，包括以下步骤：

对通道雷达波bscan的某一雷达信号Ascan采用墨西哥帽小波进行连续小波变换，得到与该雷达信号Ascan对应的复数矩阵C；

对所述复数矩阵C中的任一复数求得与之对应的模值，并得到矩阵R；所述矩阵R的横坐标为Ascan的位置，且矩阵R的纵坐标为小波的尺度；

对矩阵R归一化至0～255中，得到矩阵norm_R；

根据矩阵norm_R求得矩阵R中最大值的位置P(row，col)；

获取位置P(row，col)在复数矩阵C中的复数，并记为c1，所述复数c1的实部为c1r；

若实部c1r大于0，则将位置P(row，col)作为该雷达信号Ascan与墨西哥帽波形最接近的位置k；

若实部c1r小于等于0，则将矩阵R中的位置P(row，col)的值置为0，并重复求得矩阵R中最大值的位置P(row，col)，直至寻找到实部c1r大于0对应的位置P(row，col)；

若为找到矩阵R中最大值的位置P(row，col)，则将位置P(row，col)置0。

3.根据权利要求1或2所述的一种浅层对地雷达信号去除地面波的方法，其特征在于，所述对位置集合序列K进行平滑处理，并移除地面波，包括以下步骤：

求得位置集合序列K的一阶导数DiffK；

从位置集合序列K的起始位置开始采用滑动窗依次取一个子序列subDiffK；所述子序列subDiffK的长度为大于1的自然数：

若subDiffK[0]＝0，则不检查该子序列，取下一个子序列检查；所述subDiffK[0]表示位置集合序列K的起始位置的一阶导数；

若subDiffK[0]＝Val，则在位置集合序列K中查询一个位置点u，其subDiffK[u]＝-val；所述Val为不等0的数；所述u的值小于子序列subDiffK的长度；将subDiffK[1]到subDiffK[u]对应的序列赋值成subDiffK[0]＝0时的序列：以移除任一雷达信号Ascan中从信号开设至计算的地面波位置的信号，完成地面波移除工作。