CN110673219A - 磁异常数据的检测方法及装置 - Google Patents

Abstract

本发明公开了一种磁异常数据的检测方法及装置，所述方法包括：构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；根据磁异常数据的模型，确定一组正交检测基底；根据正交检测基底与磁传感器在检测窗口中的接收数据，计算检测窗口的函数值；根据检测窗口的函数值，判断磁异常数据是否存在。采用本发明，通过构造磁传感器圆弧运动轨迹下，检测目标磁异常数据的正交检测基底，可以实现圆弧运动下磁异常检测，相对于标准正交基分解(OBF)方法，具有更高的信噪比，检测效果更佳。

Description

磁异常数据的检测方法及装置

技术领域

本发明涉及地球物理学技术领域，尤其涉及一种磁异常数据的检测方法及装置。

背景技术

地球磁场在一定区域内可看作是均匀的，磁性目标处于地磁场中被磁化，使得该区域均匀分布的地磁场出现扰动，产生磁异常信号，而这个磁异常信号可被装有磁传感器的飞机飞过这片区域时探测到。然而，远距离的目标信号由于磁感应强度极低往往淹没在背景磁噪声中，因而需要利用有效的航空磁异常检测技术提高检测性能。

相关技术中，航空磁异常检测技术主要以能量检测方法为主，主要从目标磁场特征和背景磁场特征两方面着手。其中，基于目标磁场特征的最具代表性的研究为BorisGinzburg等人在2002年提出的标准正交基分解(OBF)方法，该方法将磁偶极子模型分解成三个正交基底线性组合的形式，然后利用背景噪声与正交基底不相关的特性，通过接收数据与三个不同正交基底的乘积求和实现目标信号能量积累，从而提高信噪比，达到弱磁信号检测的目的。

然而，OBF方法是基于运动轨迹为直线时的假设条件，没有考虑转弯运动对检测性能的影响。但在实际检测过程中，根据不同搜索策略，飞机需要进行不同程度的转弯运动，以满足搜索需求。此时，目标模型不再满足OBF方法中的标准正交基分解，造成检测性能大幅度下降。

发明内容

本发明实施例提供一种磁异常数据的检测方法及装置，用以解决现有技术中圆弧运动过程中对于目标检测效果差的问题。

一方面，本发明实施例提出一种磁异常数据的检测方法，包括：

构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

根据所述磁异常数据的模型，确定一组正交检测基底；

根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值；

根据所述检测窗口的函数值，判断所述磁异常数据是否存在。

根据本发明的一些实施例，所述构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型包括：

以所述检测目标为原点，所述圆弧运动轨迹的圆心在所述检测目标所在平面上的投影点、与所述检测目标之间的连线为Y轴，垂直所述Y轴向上为Z轴，垂直所述Y轴与Z轴为X轴，建立三轴坐标系；

基于所述三轴坐标系，确定所述检测目标与所述磁传感器的距离矢量

所述检测目标的磁矩矢量

以及地磁矢量

基于所述所述

以及所述

根据公式1-2构建所述检测目标相对于所述磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率。

根据本发明的一些实施例，所述根据所述磁异常数据的模型，确定一组正交检测基底，包括：

将所述磁异常数据的模型分解为多个基底函数；

根据每个所述基底函数的特性，选择两个所述基底函数作为一组正交检测基底。

根据本发明的一些实施例，所述根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值，包括：

根据公式3确定所述检测窗口的长度：

其中，θ_s为所述正交检测基底中圆心角范围的起始值，θ_e为所述正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

获取所述磁传感器在所述检测窗口中的接收数据；

将所述正交检测基底分别与所述接收数据乘积累加，以获得所述检测窗口的两个函数值。

在本发明的一些实施例中，所述根据所述检测窗口的函数值，判断所述磁异常数据是否存在，包括：

计算两个所述函数值的平方和以构造检测函数；

当所述检测函数在所述检测窗口出现峰值，则确定所述检测窗口内存在所述磁异常数据。

本发明实施例还提出一种磁异常数据的检测设备，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如上所述的磁异常数据的检测方法的步骤。

本发明实施例还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有信息传递的实现程序，所述程序被处理器执行时实现如上所述的磁异常数据的检测方法的步骤。

另一方面，本发明实施例还提出一种磁异常数据的检测装置，包括：

构建模块，用于构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

分析模块，用于根据所述磁异常数据的模型，确定一组正交检测基底；

计算模块，用于根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值；

判断模块，用于根据所述检测窗口的函数值，判断所述磁异常数据是否存在。

根据本发明的一些实施例，所述构建模块，用于：

以所述检测目标为原点，所述圆弧运动轨迹的圆心在所述检测目标所在平面上的投影点、与所述检测目标之间的连线为Y轴，垂直所述Y轴向上为Z轴，垂直所述Y轴与Z轴为X轴，建立三轴坐标系；

基于所述三轴坐标系，确定所述检测目标与所述磁传感器的距离矢量

所述检测目标的磁矩矢量

以及地磁矢量

基于所述

所述

以及所述

根据公式1-2构建所述检测目标相对于所述磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率。

根据本发明的一些实施例，所述分析模块，用于：

将所述磁异常数据的模型分解为多个基底函数；

根据每个所述基底函数的特性，选择两个所述基底函数作为一组正交检测基底。

根据本发明的一些实施例，所述计算模块，用于：

根据公式3确定所述检测窗口的长度：

其中，θ_s为所述正交检测基底中圆心角范围的起始值，θ_e为所述正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

获取所述磁传感器在所述检测窗口中的接收数据；

将所述正交检测基底分别与所述接收数据乘积累加，以获得所述检测窗口的两个函数值。

在本发明的一些实施例中，所述判断模块，用于：

计算两个所述函数值的平方和以构造检测函数；

当所述检测函数在所述检测窗口中出现峰值，则确定所述检测窗口内存在所述磁异常数据。

采用本发明实施例，通过构造磁传感器圆弧运动轨迹下，检测目标磁异常数据的正交检测基底，可以实现圆弧运动下磁异常检测，相对于标准正交基分解(OBF)方法，具有更高的信噪比，检测效果更佳。

上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。

附图说明

通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：

图1是本发明实施例中磁异常数据的检测方法的流程图；

图2是本发明实施例中三轴坐标系的示意图；

图3是本发明实施例中三轴坐标系的示意图；

图4是本发明实施例中基底函数示意图，其中λ＝0.1；

图5是本发明实施例中基底函数示意图，其中λ＝0.3；

图6是本发明实施例中基底函数示意图，其中λ＝0.5；

图7是本发明实施例中基底函数示意图，其中λ＝0.7；

图8是本发明实施例中基底函数示意图，其中λ＝0.9；

图9是本发明实施例一中磁传感器在三轴坐标系下的运动轨迹；

图10是本发明实施例一中仿真得到的检测目标的磁异常数据的曲线图；

图11是本发明实施例一中磁传感器的接收数据的曲线图；

图12是本发明实施例一中采用本发明实施例的磁异常数据的检测方法(OBF-circle)和传统OBF方法所得到的检测结果对比示意图；

图13是本发明实施例二中磁传感器在三轴坐标系下的运动轨迹；

图14是本发明实施例二中仿真得到的检测目标的磁异常数据的曲线图；

图15是本发明实施例二中磁传感器的接收数据的曲线图；

图16是本发明实施例二中采用本发明实施例的磁异常数据的检测方法(OBF-circle)和传统OBF方法所得到的检测结果对比示意图；

图17是本发明实施例三中磁传感器在三轴坐标系下的运动轨迹；

图18是本发明实施例三中仿真得到的检测目标的磁异常数据的曲线图；

图19是本发明实施例三中磁传感器的接收数据的曲线图；

图20是本发明实施例三中采用本发明实施例的磁异常数据的检测方法(OBF-circle)和传统OBF方法所得到的检测结果对比示意图；

图21是本发明实施例四中磁传感器在三轴坐标系下的运动轨迹；

图22是本发明实施例四中仿真得到的检测目标的磁异常数据的曲线图；

图23是本发明实施例四中磁传感器的接收数据的曲线图；

图24是本发明实施例四中采用本发明实施例的磁异常数据的检测方法(OBF-circle)和传统OBF方法所得到的检测结果对比示意图。

具体实施方式

下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。

一方面，本发明实施例提出一种磁异常数据的检测方法，如图1所示，所示检测方法包括：

S101，构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

可以理解的是，用于检测检测目标的磁传感器的运动轨迹类似于圆弧轨迹。磁传感器在运动过程中与检测目标的相对位置是时刻发生变化的，则磁传感器所能检测到的检测目标产生的磁异常数据也是随之变化的。需要基于磁传感器的圆弧轨迹，建立检测目标所产生的磁异常数据的模型。

S102，根据磁异常数据的模型，确定一组正交检测基底；

可以理解的是，磁异常数据的模型经过分解、筛选后可以确定出两个正交检测基底。

S103，根据正交检测基底与磁传感器在检测窗口中的接收数据，计算检测窗口的函数值；

需要说明的是，磁传感器是通过采样的方式采集接收数据的，即接收数据是离散型的，这里所提到的“检测窗口”可以理解为包含有至少一个接收数据。

S104，根据检测窗口的函数值，判断磁异常数据是否存在。

可以理解的是，通过对统计出的检测窗口的函数值分析，可以确定是否存在磁异常数据，以确定是否存在检测目标。

采用本发明实施例，通过构造磁传感器圆弧运动轨迹下，检测目标磁异常数据的正交检测基底，可以实现圆弧运动下磁异常检测，相对于标准正交基分解(OBF)方法，具有更高的信噪比，检测效果更佳。

在上述实施例的基础上，进一步提出各变型实施例，在此需要说明的是，为了使描述简要，在各变型实施例中仅描述与上述实施例的不同之处。

根据本发明的一些实施例，构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型包括：

如图2及图3所示，以检测目标为原点O，圆弧运动轨迹的圆心在检测目标所在平面上的投影点O'、与检测目标(即原点O)之间的连线为Y轴，垂直Y轴向上为Z轴，垂直Y轴与Z轴为X轴，建立三轴坐标系；

需要说明的是，这里所提到的“检测目标所在平面”可以理解为检测目标所在的水平面。

基于三轴坐标系，确定检测目标与磁传感器的距离矢量

检测目标的磁矩矢量

以及地磁矢量

基于

以及

根据公式1-2构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率，其值为4π*10-7H/m。

例如，假设检测目标是静止不动的，磁传感器与目标之间的高度差也是不变的，检测目标的磁矩矢量

可表示为

地磁矢量

可表示为

基于建立的三轴坐标系，检测目标与磁传感器的距离矢量可以用公式4表示：

其中，R为圆弧运动轨迹所对应的半径，θ为磁传感器在某一时间点在圆弧运动轨迹中对应的圆心角，即磁传感器的位置点在检测目标所在平面上的投影点与O'的连线、与Y轴正方向的夹角，H为O'与O的距离，Z₀为检测目标到圆弧运动轨迹所在的水平面的垂直距离。

令k₁＝R/Z₀,k₂＝H/Z₀，则公式4可以表示为公式5：

由磁偶极子模型可知，空间任一点处的检测目标的磁感应强度可以用公式1表示：

由于检测目标的磁感应强度远小于地磁场，故磁传感器测得的检测目标的磁感应强度可以用公式2表示：

令

将公式1、5代入公式2可得磁异常数据的模型，参考公式6：

其中，θ是随时间变化的，即θ＝θ(t)，α_i为常数。

根据本发明的一些实施例，根据磁异常数据的模型，确定一组正交检测基底，包括：

将磁异常数据的模型分解为多个基底函数；

根据每个基底函数的特性，选择两个基底函数作为一组正交检测基底。

例如，假设磁异常数据的模型如公式6所示，将公式6进行分解可得，

α_i表达式如下所示：

φ_i(θ)表达式为：

其中，

φ_i(θ)则公式6中磁异常数据的模型对应的基底函数。

根据三角函数的特点，上述5个基底函数需满足如下正交关系：

由于λ在不同情况下的值不同，对应的φ_i(θ)也不同，图4-图8为5个基底函数在5个不同λ值下的函数图，图4-图8中的λ分别取0.1、0.3、0.5、0.7、0.9。从图4-图8可以看出，φ_i(θ),i＝1,3,5彼此相关，φ_j(θ),i＝2,4彼此相关(值越大，相关系数λ越接近1)。正交检测基底可以从φ_i(θ),i＝1,3,5和φ_j(θ),i＝2,4中选取。

进一步的，可以根据基于随着λ和θ变化，收敛速度最快的选取原则选择两个基底函数作为正交检测基底。例如，可以选择φ₁(θ)和φ₂(θ)作为一组正交检测基底，并令λ＝0.3，对应的，θ＝[-4π/9,4π/9]。

构造正交检测基底的目的是为了最大化正交检测基底与磁异常数据的乘积，最小化正交检测基底与噪声数据的乘积，从而提高检测信噪比。

根据本发明的一些实施例，根据正交检测基底与磁传感器在检测窗口中的接收数据，计算检测窗口的函数值，包括：

根据公式3确定检测窗口的长度：

其中，θ_s为正交检测基底中圆心角范围的起始值，θ_e为正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

例如，当正交检测基底中圆心角范围为θ＝[-4π/9,4π/9]，则

获取磁传感器在检测窗口中的接收数据；

将正交检测基底分别与接收数据乘积累加，以获得检测窗口的两个函数值。

例如，根据公式7计算检测窗口的两个函数值，

其中，m表示第m个时刻的检测窗口，x_m(n)表示第m个时刻的检测窗口中第n个接收数据。

进一步的，采样角度间隔Δθ可以根据公式8计算获得：

Δθ＝v/(f_sR) 公式8；

其中，v为磁传感器的运行速度，f_s为磁传感器的接收数据采集频率，即采样频率。

在本发明的一些实施例中，根据检测窗口的函数值，判断磁异常数据是否存在，包括：

计算两个函数值的平方和以构造检测函数；

例如，根据公式9构造检测函数，

当检测函数在检测窗口出现峰值，则确定检测窗口内存在磁异常数据。

可以理解的是，从公式9可知，E(m)与时间有关，当E(m)出现明显峰值，即可确定磁异常数据存在。

本发明实施例还提出一种磁异常数据的检测设备，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，计算机程序被处理器执行时实现如上的磁异常数据的检测方法的步骤。

本发明实施例还提出一种计算机可读存储介质，计算机可读存储介质上存储有信息传递的实现程序，程序被处理器执行时实现如上的磁异常数据的检测方法的步骤。

另一方面，本发明实施例还提出一种磁异常数据的检测装置，包括：

构建模块，用于构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

分析模块，用于根据磁异常数据的模型，确定一组正交检测基底；

计算模块，用于根据正交检测基底与磁传感器在检测窗口中的接收数据，计算检测窗口的函数值；

判断模块，用于根据检测窗口的函数值，判断磁异常数据是否存在。

根据本发明的一些实施例，构建模块，用于：

以检测目标为原点，圆弧运动轨迹的圆心在检测目标所在平面上的投影点、与检测目标之间的连线为Y轴，垂直Y轴向上为Z轴，垂直Y轴与Z轴为X轴，建立三轴坐标系；

基于三轴坐标系，确定检测目标与磁传感器的距离矢量

检测目标的磁矩矢量

以及地磁矢量

基于

以及

根据公式1-2构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率。

根据本发明的一些实施例，分析模块，用于：

将磁异常数据的模型分解为多个基底函数；

根据每个基底函数的特性，选择两个基底函数作为一组正交检测基底。

根据本发明的一些实施例，计算模块，用于：

根据公式3确定检测窗口的长度：

其中，θ_s为正交检测基底中圆心角范围的起始值，θ_e为正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

获取磁传感器在检测窗口中的接收数据；

将正交检测基底分别与接收数据乘积累加，以获得检测窗口的两个函数值。

在本发明的一些实施例中，判断模块，用于：

计算两个函数值的平方和以构造检测函数；

当检测函数在检测窗口中出现峰值，则确定检测窗口内存在磁异常数据。

采用本发明实施例，通过构造磁传感器圆弧运动轨迹下，检测目标磁异常数据的正交检测基底，可以实现圆弧运动下磁异常检测，相对于标准正交基分解(OBF)方法，具有更高的信噪比，检测效果更佳。

下面通过四个具体的实施例，进一步的说明采用本发明实施例的磁异常检测方法相对于现有技术中的标准正交基分解(OBF)方法具有更佳的检测效果。需要说明的是，下面实施例中，磁传感器承载于飞机，可以将飞机等同于磁传感器。

实施例一

飞机做180°的圆弧转弯运动，圆弧的半径R为1000m，飞行速度v为30m/s，采样率f_s为10Hz，坐标系原点O和圆弧的圆心在检测目标所在的平面上的投影点O'的距离H＝800m，检测目标到飞行平面的垂直距离Z₀＝200m，检测目标的相关参数及地磁的相关参数如表1所示。

表1检测目标的相关参数及地磁的相关参数

通过仿真，给出磁传感器(飞机)在三轴坐标系下的运动轨迹如图9所示，仿真得到的检测目标的磁异常数据的曲线图如图10所示；加入高斯白噪声作为背景噪声(信噪比为-20dB)，得到的接收数据的曲线如图11所示，图12分别显示了利用本发明实施例中磁异常数据的检测方法(实线)和传统OBF方法(虚线)进行磁异常检测得到的结果。

实施例二

飞机做180°的圆弧转弯运动，圆弧的半径R为500m，飞行速度v为30m/s，采样率f_s为10Hz，坐标系原点O和圆弧的圆心在检测目标所在的平面上的投影点O'的距离H＝100m，检测目标到飞行平面的垂直距离Z₀＝200m，检测目标的相关参数及地磁的相关参数如表2所示。

表2检测目标的相关参数及地磁的相关参数

通过仿真，给出飞机在目三轴坐标系下的运动轨迹如图13所示，仿真得到的检测目标的磁异常数据的曲线如图14所示；加入高斯白噪声作为背景噪声(信噪比为-20dB)，得到的接收数据的曲线图如图15所示，图16分别显示了利用本发明实施例中磁异常数据的检测方法(实线)和传统OBF方法(虚线)进行磁异常检测得到的结果。

实施例三

飞机做180°的圆弧转弯运动，圆弧的半径R为1000m，飞行速度v为30m/s，采样率f_s为10Hz，坐标系原点O和圆弧的圆心在检测目标所在的平面上的投影点O'的距离H＝1100m，检测目标到飞行平面的垂直距离Z₀＝300m，检测目标的相关参数及地磁的相关参数如表3所示。

表3检测目标的相关参数及地磁的相关参数

通过仿真，给出飞机在目三轴坐标系下的运动轨迹如图17所示，仿真得到的检测目标的磁异常数据的曲线图如图18所示；加入高斯白噪声作为背景噪声(信噪比为-20dB)，得到的接收数据的曲线图如图19所示，20分别显示了利用本发明实施例中磁异常数据的检测方法(实线)和传统OBF方法(虚线)进行磁异常检测得到的结果。

实施例四

飞机做90°的圆弧转弯运动，圆弧的半径R为1000m，飞行速度v为30m/s，采样率f_s为10Hz，坐标系原点O和圆弧的圆心在检测目标所在的平面上的投影点O'的距离H＝800m，检测目标到飞行平面的垂直距离Z₀＝200m，检测目标的相关参数及地磁的相关参数如表4所示。

表4检测目标的相关参数及地磁的相关参数

通过仿真，给出飞机在目三轴坐标系下的运动轨迹如图21所示，仿真得到的检测目标的磁异常数据的曲线图如图22所示；加入高斯白噪声作为背景噪声(信噪比为-20dB)，得到的接收数据的曲线图如图23所示，图24分别显示了利用本发明实施例中磁异常数据的检测方法(实线)和传统OBF方法(虚线)进行磁异常检测得到的结果。

从上述四个实施例可知，通过采用本发明实施例的磁异常数据的检测方法比传统OBF方法具有更好的检测效果。

需要说明的是，在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (12)

1.一种磁异常数据的检测方法，其特征在于，包括：

构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

根据所述磁异常数据的模型，确定一组正交检测基底；

根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值；

根据所述检测窗口的函数值，判断所述磁异常数据是否存在。

2.如权利要求1所述的方法，其特征在于，所述构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型包括：

以所述检测目标为原点，所述圆弧运动轨迹的圆心在所述检测目标所在平面上的投影点、与所述检测目标之间的连线为Y轴，垂直所述Y轴向上为Z轴，垂直所述Y轴与Z轴为X轴，建立三轴坐标系；

基于所述三轴坐标系，确定所述检测目标与所述磁传感器的距离矢量

所述检测目标的磁矩矢量以及地磁矢量

基于所述

所述

以及所述

根据公式1-2构建所述检测目标相对于所述磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率。

3.如权利要求1所述的方法，其特征在于，所述根据所述磁异常数据的模型，确定一组正交检测基底，包括：

将所述磁异常数据的模型分解为多个基底函数；

根据每个所述基底函数的特性，选择两个所述基底函数作为一组正交检测基底。

4.如权利要求1所述的方法，其特征在于，所述根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值，包括：

根据公式3确定所述检测窗口的长度：

其中，θ_s为所述正交检测基底中圆心角范围的起始值，θ_e为所述正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

获取所述磁传感器在所述检测窗口中的接收数据；

将所述正交检测基底分别与所述接收数据乘积累加，以获得所述检测窗口的两个函数值。

5.如权利要求4所述的方法，其特征在于，所述根据所述检测窗口的函数值，判断所述磁异常数据是否存在，包括：

计算两个所述函数值的平方和以构造检测函数；

当所述检测函数在所述检测窗口出现峰值，则确定所述检测窗口内存在所述磁异常数据。

6.一种磁异常数据的检测装置，其特征在于，包括：

构建模块，用于构建检测目标相对于磁传感器的圆弧运动轨迹所产生的磁异常数据的模型；

分析模块，用于根据所述磁异常数据的模型，确定一组正交检测基底；

计算模块，用于根据所述正交检测基底与所述磁传感器在检测窗口中的接收数据，计算所述检测窗口的函数值；

判断模块，用于根据所述检测窗口的函数值，判断所述磁异常数据是否存在。

7.如权利要求6所述的装置，其特征在于，所述构建模块，用于：

以所述检测目标为原点，所述圆弧运动轨迹的圆心在所述检测目标所在平面上的投影点、与所述检测目标之间的连线为Y轴，垂直所述Y轴向上为Z轴，垂直所述Y轴与Z轴为X轴，建立三轴坐标系；

基于所述三轴坐标系，确定所述检测目标与所述磁传感器的距离矢量

所述检测目标的磁矩矢量

以及地磁矢量

基于所述

所述

以及所述

根据公式1-2构建所述检测目标相对于所述磁传感器的圆弧运动轨迹所产生的磁异常数据的模型：

其中，μ₀为空气中的磁导率。

8.如权利要求6所述的装置，其特征在于，所述分析模块，用于：

将所述磁异常数据的模型分解为多个基底函数；

根据每个所述基底函数的特性，选择两个所述基底函数作为一组正交检测基底。

9.如权利要求6所述的装置，其特征在于，所述计算模块，用于：

根据公式3确定所述检测窗口的长度：

其中，θ_s为所述正交检测基底中圆心角范围的起始值，θ_e为所述正交检测基底中圆心角范围的终值，Δθ为采样角度间隔；

获取所述磁传感器在所述检测窗口中的接收数据；

将所述正交检测基底分别与所述接收数据乘积累加，以获得所述检测窗口的两个函数值。

10.如权利要求9所述的装置，其特征在于，所述判断模块，用于：

计算两个所述函数值的平方和以构造检测函数；

当所述检测函数在所述检测窗口中出现峰值，则确定所述检测窗口内存在所述磁异常数据。

11.一种磁异常数据的检测设备，其特征在于，包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至5中任一项所述的磁异常数据的检测方法的步骤。

12.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有信息传递的实现程序，所述程序被处理器执行时实现如权利要求1至5中任一项所述的磁异常数据的检测方法的步骤。

Images