[发明专利]一种全自动的航空电磁数据调平方法及装置

权利要求书

1.一种全自动的航空电磁数据调平方法，其特征在于：包括以下步骤，

1).对待测区域进行航测并记录航测数据D，航测所得的测区数据中包含多条测线数据；

2).录入航测数据，将航测数据中异常区域数据剔除，剔除后得到测区数据记为

3).所述的D_{ext_gru}经过变分模态分解后，变分模态分解后将分解得到2个模态分量u₁(t)和u₂(t)；

4).求解调平误差D_{ext_error}，模态分量u₁(t)代入下式

得到调平误差D_{ext_error}

5).求解调平结果

D_result＝D-D_error，得到各测区数据调平后调平航测数据结果；

输出测区数据调平结果D_result

6).依据步骤5)中输出的调平航测数据结果D_result进行成图，得到调平后航空电磁数据图。

2.根据权利要求1所述的一种全自动的航空电磁数据调平方法，其特征在于：步骤2)中，所述的将航测数据中异常区域数据剔除，还包括，

在测线数据的两端数据分别进行镜像周期延拓，并利用格拉布斯准则进行异常区域数据的剔除处理；

其中数据延拓过程中，设定各测线两端的延拓测点数为：

此时各测线包含测点为：

各测线数据经镜像周期延拓后，测区数据表示为：

3.根据权利要求2所述的一种全自动的航空电磁数据调平方法，其特征在于：数据延拓过程中计算测点数据与该测线数据均值的残差与格拉布斯准则，当测点数据与该测线数据均值的残差满足时记为异常测点数据，并剔除且由该测线的中值数据代替：

为该测线数据的均值，g(num，a)为临界值系数，取决于测量次数num和显著性水平a，为测线数据的标准差；

测线数据剔除所述的异常测点数据后，表示为，

4.根据权利要求3所述的一种全自动的航空电磁数据调平方法，其特征在于：

并依次对测区各测线数据进行异常区域数据的剔除处理，剔除处理后的测区数据记为

5.根据权利要求1所述的一种全自动的航空电磁数据调平方法，其特征在于：步骤3)中所述的D_{ext_gru}经过变分模态分解，还包括

①变分模态分量u(t)、变分模态分量的中心频率ω和拉格朗日乘子λ，代入变分模态分量计算式循环更新迭代，在第m次时候得到：

其中，δ(t)为狄拉克函数，*表示卷积；

②变换到谱域并积分可得到第m次迭代的变分模态分量：

③变分模态分量的中心频率计算公式如下：

④变换到谱域，从而得到第m次迭代的变分模态分量的中心频率

⑤拉格朗日乘子更新迭代公式如下：

⑥设置迭代中止条件，符合中止条件时迭代终止。

对所得到的对变分模态分量和进行傅里叶反变换，得到时域的分解结果u₁(t)和u₂(t)，至此完成测区数据D_{ext_gru}的变分模态分解。

6.根据权利要求5所述的一种全自动的航空电磁数据调平方法，其特征在于：所述的变换到谱域具体的为，

将步骤①中得到的式变换为谱域式：

sgn为符号函数。

7.根据权利要求6所述的一种全自动的航空电磁数据调平方法，其特征在于：所述的积分为构建埃尔米特矩阵，并将所述的谱域式的两项结果写为非负频率上的半空间积分，同时令ω＝ω+ω₁，得到：