[發(fā)明專利]一種基于最小二乘克希霍夫偏移的探月雷達(dá)數(shù)據(jù)處理方法

在審

• 申請?zhí)枺?/td>: 201910125903.5
• 申請日: 2019-02-20
• 公開（公告）號: CN109765626A
• 公開（公告）日: 2019-05-17
• 發(fā)明（設(shè)計(jì)）人: 李靜
• 申請（專利權(quán)）人: 吉林大學(xué)
• 主分類號: G01V3/12
• 分類號: G01V3/12
• 代理公司: 成都眾恒智合專利代理事務(wù)所(普通合伙) 51239
• 代理人: 劉華平

1.一種基于最小二乘克?；舴蚱频奶皆吕走_(dá)數(shù)據(jù)處理方法，其特征在于，包括如下步驟：

(S1)讀入探月雷達(dá)數(shù)據(jù)，并對數(shù)據(jù)進(jìn)行預(yù)處理；

(S2)進(jìn)行速度譜分析，確定月球內(nèi)部主要層結(jié)構(gòu)的位置；

(S3)開展克?；舴虻奶皆吕走_(dá)數(shù)據(jù)偏移處理；

(S4)加入最小二乘迭代進(jìn)行數(shù)據(jù)偏移算法；

(S5)輸出探月雷達(dá)數(shù)據(jù)成像結(jié)果。

2.根據(jù)權(quán)利要求1所述的一種基于最小二乘克?；舴蚱频奶皆吕走_(dá)數(shù)據(jù)處理方法，其特征在于，所述步驟(S1)中數(shù)據(jù)預(yù)處理包括以下過程：

(1)數(shù)據(jù)讀?。焊鶕?jù)標(biāo)準(zhǔn)存儲格式存儲9段數(shù)據(jù)，并逐卷讀取道集及位置信息；

(2)數(shù)據(jù)拼接：將讀取到的9段數(shù)據(jù)拼接成一組完整剖面數(shù)據(jù)；

(3)道集調(diào)整：對得到的剖面數(shù)據(jù)根據(jù)連續(xù)同相軸進(jìn)行道集調(diào)整；

(4)道集擇取：將調(diào)整后的道集進(jìn)行篩選，選擇數(shù)據(jù)質(zhì)量好的道集；

(5)時(shí)間差調(diào)整；

(6)無用數(shù)據(jù)刪除：將調(diào)整后的數(shù)據(jù)信噪比低的無用數(shù)據(jù)部分刪除；

(7)帶通濾波：進(jìn)行帶通濾波，去除數(shù)據(jù)中的低頻及高頻噪聲；

(8)背景去除；

(9)自動(dòng)增益。

3.根據(jù)權(quán)利要求2所述的一種基于最小二乘克?；舴蚱频奶皆吕走_(dá)數(shù)據(jù)處理方法，其特征在于，所述步驟(S2)中進(jìn)行速度譜分析，確定月球內(nèi)部主要層結(jié)構(gòu)的位置是通過轉(zhuǎn)換數(shù)據(jù)由共偏移距道集為中心點(diǎn)數(shù)據(jù)開展動(dòng)校正，在水平迭代處理過程中開展速度譜分析，通過時(shí)頻分析技術(shù)，計(jì)算出月球內(nèi)部主要雷達(dá)信號反射層的位置。

4.根據(jù)權(quán)利要求3所述的一種基于最小二乘克?；舴蚱频奶皆吕走_(dá)數(shù)據(jù)處理方法，其特征在于，所述步驟(S3)中進(jìn)行克?；舴虻奶皆吕走_(dá)數(shù)據(jù)偏移處理的主要流程如下：

線性正演模擬算子L_～滿足

p＝L_～m (1)

其中，p表示建模數(shù)據(jù)的矢量，m表示模型反射率模型矢量，克?；舴蚱剖褂霉?1)中的正演模型算子的轉(zhuǎn)置：

m_k＝L_～^TP (2)

將公式(1)帶入公式(2)，可得

m_k＝L_～^TL_～m (3)

其中，L_～^TL_～表示Hessian矩陣，對應(yīng)L_～^TL_～的積分算子的顯式形式表達(dá)公式為：

k(x,x′)＝∫dsh_s(s)A_sxA_sx′×∫drh_r(r)A_xrA_x′rR(τ_sx+τ_xr-τ_sx′-τ_x′r) (4)

其中，x和x′分別表示與m_k和m相關(guān)的位置，τ_sx表示從源s處到x處反射器的波傳播時(shí)間，τ_xr表示從反射器到r處接收器的波傳播時(shí)間，A_sx和A_xr表示是振幅項(xiàng)，R表示時(shí)間(τ_sx+τ_xr-τ_sx′-τ_x′r)延遲的相關(guān)時(shí)間交叉系數(shù)，h_s(s)和h_r(r)分別表示源和接收器采樣的函數(shù)。