[發明專利]一種淺層對地雷達信號去除地面波的方法

在審

• 申請號: 202011248058.X
• 申請日: 2020-11-10
• 公開（公告）號: CN112162279A
• 公開（公告）日: 2021-01-01
• 發明（設計）人: 黃敏;余世杰;楊輝;桂仲成
• 申請（專利權）人: 上海圭目機器人有限公司
• 主分類號: G01S13/88
• 分類號: G01S13/88;G01S7/02;G01S7/36
• 代理公司: 成都佳劃信知識產權代理有限公司 51266
• 代理人: 幸偉山

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中從信號開設至計算的地面波位置的信號，完成地面波移除工作。