CN114777766A - 一种基于磁场梯度张量的目标定位方法及装置 - Google Patents
一种基于磁场梯度张量的目标定位方法及装置 Download PDFInfo
- Publication number
- CN114777766A CN114777766A CN202210433815.3A CN202210433815A CN114777766A CN 114777766 A CN114777766 A CN 114777766A CN 202210433815 A CN202210433815 A CN 202210433815A CN 114777766 A CN114777766 A CN 114777766A
- Authority
- CN
- China
- Prior art keywords
- observation
- magnetic field
- field gradient
- target body
- positioning
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 239000013598 vector Substances 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 12
- 230000004807 localization Effects 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 2
- 238000004891 communication Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
- G01C21/08—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/022—Measuring gradient
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Automation & Control Theory (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种基于磁场梯度张量的目标定位方法及装置,其包括:获取观测站记录的观测数据;基于记录的观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于距离矢量以及在T时刻的观测点P的观测位置确定在T时刻目标体的位置。其中,获取的观测数据包括观测点P的两个相邻观测位置在同一时刻T或同一观测点P在T时刻的两个相邻时刻记录的磁场矢量及磁场梯度张量,或包括:观测点P及其两个相邻观测位置在同一时刻T或同一观测点P在T时刻及其两个相邻时刻记录的磁场梯度张量。其中,推导出的定位公式均是消除了地球背景磁场的目标定位公式,解决了现有技术中因地球背景磁场产生的定位误差问题,实现目标体定位。
Description
技术领域
本发明属于地磁学技术领域,具体涉及一种基于磁场梯度张量的目标定位方法及装置。
背景技术
在地磁学中,磁梯度张量除用于勘探地球地质结构构造,寻找油气、矿产资源以外,还被广泛应用于水下磁性目标体定位。现有的磁定位方法一般是在建立在忽略背景场影响的重力或磁异常及其梯度张量的基础上实现。然而,忽略背景场影响的磁定位方法的定位精度必然受到影响,有待进一步提升。
实际上,空间里面一点上的总重力场或总磁场T一般包括三个成分:地球磁场Be、地下物质体引起的异常场Ba以及未知场Bu:
T=B0+Ba+Bu
由于地下目标异常体通常淹没在相对较大的地球磁场Be(中纬度地区约为50000nT),直接测量异常场Ba通常具有挑战性,必须从测量结果中消除地球磁场Be。但地球磁场相对均匀,与人造的磁性目标体相比,其梯度通常较小,其相对人造磁性目标体所产生的梯度场来说可以忽略不计,未知场Bu通常也是被忽略不计,即
因此,为了实现基于磁场梯度张量的目标定位,怎么从强大的区域背景场中提取异常场是十分重要的,也是本领域的技术难点。为此,实有必要设计一种消除地球背景磁场影响的基于磁场梯度张量的目标定位方法,以解决现有技术中因地球背景磁场产生的定位误差问题。
发明内容
本发明的目的是解决现有技术中因地球背景磁场产生的定位误差问题,进而提供一种基于磁场梯度张量的目标定位方法及装置,用以实现目标体定位。本发明抛开了传统面积性数据采集方式,其观测磁场梯度张量,并利用同一观测站在相邻时刻或相邻观测站在同一时刻记录两组或三组磁场梯度张量以及采用差分的手段推导出消除了地球背景磁场的目标定位公式,实现目标体定位。且所述方法适用于静止目标体定位,也适用于运动目标体定位;适用于基于固定基站的目标定位以及适用于基于运动平台的目标定位。
一方面,本发明提供的一种基于磁场梯度张量的目标定位方法,包括以下步骤:
步骤S1:获取观测站记录的观测数据;
步骤S2:基于步骤S1记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置;
定位公式如下:
ri=[Gi-1-Gi+1]-1·[n(Ti+1-Ti-1)-Gi+1·d]
或
ri=[Gi-1-Gi+1]-1·[nGi-Gi+1]·d
其中,ri表示在T时刻所述观测点P与所述目标体之间的距离矢量,Ti-1、Ti+1为观测点P的两个相邻观测位置i-1、i+1在同一时刻T记录的磁场矢量或者为观测点P在T时刻的两个相邻时刻i-1、i+1记录的磁场矢量;Gi以及Gi-1、Gi+1依次为所述观测点P的对应观测位置i及两个相邻观测位置i-1、i+1在同一时刻T记录的磁场梯度张量或者为所述观测点P在T时刻及两个相邻时刻i-1、i+1记录的磁场梯度张量;
d=ri+1-ri-1表示两个观测位置i+1、i-1或两个时刻i+1、i-1对应的观测点与目标体之间的距离矢量的差值,n为目标体的构造指数。对应观测位置时,d=ri+1-ri-1表示观测点P的两个相邻观测位置i+1、i-1在时刻T对应的观测点与目标体之间的距离矢量的差值。对应观测时刻时,d=ri+1-ri-1表示两个时刻i+1、i-1对应的观测点P与目标体之间的距离矢量的差值。
可选地,在任意时刻任意观测点与目标体的距离矢量R表示为:
R=[xs-x* ys-y* zs-z*]
式中,x*、y*、z*分别表示观测点在xyz轴上的坐标,xs、ys、zs分别表示目标体在xyz轴上的坐标。
可选地,所述目标体的构造指数n的取值范围为[1,3]。
可选地,所述目标体为球体、或长方体、或二度体时,所述构造指数的取值对应为3,2,1。
可选地,所述目标体位于太空、或地表、或地表以上、或地表以下、或水面、或水面以上、或水面以下。
可选地,所述观测站为固定基站或者位于运动平台上;所述目标体为静止目标体或运动目标体。
第二方面,本发明提供一种基于所述目标定位方法的磁场梯度张量测量装置,所述磁场梯度张量测量装置包括8个磁传感器,所述8个磁传感器分别分布在正方体的8个顶角;
其中,所述正方体上非相邻的两个侧面的中心点及其正方体的中心点构成三个相邻观测位置时,磁传感器F1、F2、F3和F4组成一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi-1;磁传感器F5、F6、F7和F8组成另一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi+1;磁传感器F1-F8组成立方体磁场梯度测量系统,用于测定磁场梯度张量Gi。
第三方面,本发明提供一种基于所述目标定位方法的装置,其包括:
观测模块,用于记录观测数据;
定位模块,用于基于记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。
第四方面,本发明提供一种电子终端,其包括:一个或多个处理器和存储了一个或多个计算机程序的存储器;
所述计算机程序被处理器调用以实现:前述一种基于磁场梯度张量的目标定位方法的步骤。
第五方面,本发明提供一种可读存储介质,其存储了计算机程序,所述计算机程序被处理器调用以实现:前述一种基于磁场梯度张量的目标定位方法的步骤。
有益效果
本发明技术方案提供的上述一种基于磁场梯度张量的目标定位方法及其装置巧妙地消除了地球背景磁场的影响,克服了现有技术中难以从强大的区域背景场中提取异常场,从而导致定位精度有待提高的技术难点。其中,本发明抛开了传统面积性数据采集方式,其观测磁场梯度张量,并利用同一观测站在相邻时刻或相邻观测站在同一时刻记录两组或三组磁场梯度张量并采用差分的手段推导出消除了地球背景磁场的定位公式,再基于所述定位公式实现运动目标体的定位。解决了现有技术中因地球背景磁场产生的定位误差问题,实现了可靠、精准的定位追踪。
其中,本发明所述目标定位方法适用于静止目标体也适用于运动目标体,从而利用所述目标定位方法可以实时地定位监测运动目标体的位置,亦可以求出运动目标体的运动轨迹。本发明所述目标定位方法也适用于静止平台和运动平台,应用空间极为宽广。
尤其是,本发明推导出的定位公式中,仅需要测量三个连续的位场梯度张量也可实现定位,除了实现消除背景场影响功能之外,还不再需要测量位场矢量,进一步满足应用需求。
附图说明
图1为本发明实施例提供的观测装置示意图。其中F1、F2、F3和F4组成正方形磁场梯度测量系统Gi-1,F5、F6、F7和F8组成正方形磁场梯度测量系统Gi+1,F1-F8组成立方体磁场梯度测量系统Gi。
图2为本发明实施例提供的一种基于磁场梯度张量的目标定位方法的流程示意图。
图3为基于磁场梯度张量的差分定位与追踪方法在x-y平面结果。图中实线1表示三点差分定位的误差结果,图中实线2表示两点差分定位的误差结果。
具体实施方式
本发明提供的一种基于磁场梯度张量的目标定位方法及装置,其目的是消除地球背景磁场的影响,利用观测或计算得到的磁场梯度张量实现目标定位。其中,所述目标定位方法不仅仅适用于静止目标体,还适用于运动目标体。对应的观测站可以是静止基站,也可以是基于运动平台的移动观测站。下面将结合实施例对本发明做进一步的说明。
本发明提供的目标定位方法是基于差分磁场梯度张量来实现的,其原理如下:
采用运动平台探测定位相对静止的物体时,设运动平台第i个测量位置相对于磁性目标体的位置参数为ri。假设区域背景场在短距离或短时间内是相对均匀的,根据欧拉反褶积原理,由区域背景场B0、测量的场B和梯度张量G,可得:
Gr=-n(B-B0) (1)
其中,r=(xs-x ys-y zs-z)T,n为构造指数。为了消除区域背景场B0的影响,在相邻的两个测量点上或者两次测量上得到Bi-1、Bi+1和Gi-1、Gi+1,则公式可以写成如下形式:
Gi-1ri-1-Gi+1ri+1=-n(Bi-1-Bi+1) (2)
设两个点之间的距离矢量为d=(dx dy dz)T,将ri+1=ri-1+d代入公式,则有:
Gi-1ri-1-Gi+1(ri-1+d)=-n(Bi-1-Bi+1 )(3)
将公式的右端项移动到左边,则有:
(Gi-1-Gi+1)ri-1+n(Bi-1-Bi+1)=Gi+1d (4)
进一步可得到如下形式:
假设地磁场或地球重力场分布相对均匀,即Bi+1-Bi-1≈Ti+1-Ti-1,根据上述差分定位公式,利用第i-1和第i+1个测量位置点测量的磁场矢量Ti-1、Ti+1和位场梯度张量Gi-1、Gi+1得到ri:
ri=[Gi-1-Gi+1]-1·[n(Ti+1-Ti-1)-Gi+1·d] (6)
其中,d=ri+1-ri-1表示两个观测位置i+1、i-1或两个时刻i+1、i-1对应的观测点与目标体之间的距离矢量的差值。n为目标体的构造指数,构造指数的取值范围为:[1,3]。本实施中设定所述目标体为球体、或长方体、或二度体时,所述构造指数的取值对应为3,2,1。
考虑传感器的采样频率和运动平台的运动速度,对于常规探测距离满足平台采样间隔位移远远小于其相对目标的位移,即(ri+1-ri)<<ri。由于磁梯度张量Gi为磁场矢量Ti的一阶空间导数,因此三个连续测量点存在以下关系式:
Gi·(ri+1-ri-1)=Ti+1-Ti-1 (7)
将等式代入等式,则可以得到三点差分目标定位与追踪计算公式:
ri=[Gi-1-Gi+1]-1·[nGi-Gi+1]·d (8)
从等式(8)可以看出,该定位算法仅仅需要测量三个连续的位场梯度张量,而不需要测量位场矢量。且有效地消除了背景场的影响。但是需要先验给出目标体的构造指数。需要指出的是,上述推理出的公式(7)、(8)也可以应用到固定观测基站定位和追踪运动目标体。
基于上述原理性陈述,如图2所示,本发明提供的一种基于磁场梯度张量的目标定位方法,包括以下步骤:
步骤S1:获取观测站记录的观测数据。
本实施例以T时刻的观测点P为例进行简述,若获取的观测数据包括:所述观测点P的两个相邻观测位置在同一时刻T或同一观测点P在T时刻的两个相邻时刻记录的磁场矢量及磁场梯度张量,则选用对应公式(7)计算出T时刻观测点P与目标体之间的距离矢量。若获取的观测数据包括:观测点P及其两个相邻观测位置在同一时刻T或同一观测点P在T时刻及其两个相邻时刻记录的磁场梯度张量,则选用对应公式(8)计算出T时刻观测点P与目标体之间的距离矢量。
应当理解,不论是采用公式(7)或公式(8)计算距离矢量,均需要测量磁场梯度张量。如图1所示,本实施例中设定的磁场梯度张量测量装置包括8个磁传感器,所述8个磁传感器分别分布在正方体的8个顶角。其中,正方体的两个正对的端面分别设有4个磁传感器,因此正方体的两个正对的端面中心点及其正方体的中心点构成三个相邻观测位置时,磁传感器F1、F2、F3和F4组成一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi-1;磁传感器F5、F6、F7和F8组成另一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi+1;磁传感器F1-F8组成立方体磁场梯度测量系统,用于测定磁场梯度张量Gi。
本实施例中,由图1中所示F1、F2、F3和F4磁传感器的数据计算出磁场梯度张量Gi-1,对应公式为:
由图1中所示F5、F6、F7和F8磁传感器的数据计算出磁场梯度张量Gi+1,对应公式为:
式中,L为相邻两个磁传感器之间的距离,Bjx、Bjy、Bjz分别表示第j个磁传感器测量得到x,y,z方向的磁场强度,j∈{1,2,3...8}。
Gi由F1-F8组成立方体磁场梯度测量系统获得,其中,梯度张量的各分量表示如下:
其中,Bxx,Bxy,Bxz,Byy,Byz,Byx,Bzz,Bzx,Bzy均为磁场梯度张量Gi的元素。
应当理解,上述磁场梯度张量测量装置为本实施例的举例说明,本发明并不局限于采用该唯一的方式来测量磁场梯度张量,其他可行的实施例中,可以采取其他技术手段来测定磁场梯度张量。
步骤S2:基于步骤S1记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。
如前文所述,本发明根据获取的观测数据采用公式(7)或者公式(8)计算出距离矢量,进而基于已知的观测点位置计算出目标体的位置。其中,观测点与目标体的距离矢量表示如下:
R=[xs-x* ys-y* zs-z*] (12)
式中,x*、y*、z*分别表示观测点在xyz轴上的坐标,xs、ys、zs分别表示目标体在xyz轴上的坐标。
另一方面,本发明还提供一种基于所述目标定位方法的装置,其包括:观测模块以及定位模块。
其中,观测模块用于记录观测数据;定位模块用于基于记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。应当理解,观测点P以及T时刻均是指代任一观测点位置以及任一观测时刻。
需要说明的是,上述观测模块以及定位模块均可以以硬件方式来实现,也可以以软件方式来实现。譬如,观测模块可以理解为获取观测数据的传感设备或者基于传感设备采集的数据计算相关观测数据的软件功能模块。定位模块可以理解为具备数据计算以及处理功能的处理器等硬件设备,也可以理解为具备数据计算以及处理功能的软件功能模块。即上述单元模块的具体实现过程参照方法内容,本发明在此不进行具体的赘述,且上述功能模块单元的划分仅仅是一种逻辑功能的划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。同时,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
另一方面,本发明还提供一种电子终端,其包括:一个或多个处理器和存储了一个或多个计算机程序的存储器;
所述计算机程序被处理器调用以实现:前述一种基于磁场梯度张量的目标定位方法的步骤。譬如具体执行:
步骤S1:获取观测站记录的观测数据;
步骤S2:基于步骤S1记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。
该终端还包括:通信接口,用于与外界设备进行通信,进行数据交互传输。譬如与外部观测站的传感设备通讯,获取观测数据。
其中,存储器可能包含高速RAM存储器,也可能还包括非易失性除颤器,例如至少一个磁盘存储器。
如果存储器、处理器和通信接口独立实现,则存储器、处理器和通信接口可以通过总线相互连接并完成相互间的通信。所述总线可以是工业标准体系结构总线,外部设备互联总线或扩展工业标准体系结构总线等。所述总线可以分为地址总线、数据总线、控制总线等。
可选的,在具体实现上,如果存储器、处理器和通信接口集成在一块芯片上,则存储器、处理器即通信接口可以通过内部接口完成相互之间的通信。
各个步骤的具体实现过程请参照前述方法的阐述。
应当理解,在本发明实施例中,所称处理器可以是中央处理单元(CentralProcessing Unit,CPU),该处理器还可以是其他通用处理器、数字信号处理器(DigitalSignal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器可以包括只读存储器和随机存取存储器,并向处理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如,存储器还可以存储设备类型的信息。
另一方面,本发明还提供一种可读存储介质,其存储了计算机程序,所述计算机程序被处理器调用以实现:前述一种基于磁场梯度张量的目标定位方法的步骤。譬如具体执行:
步骤S1:获取观测站记录的观测数据;
步骤S2:基于步骤S1记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。
各个步骤的具体实现过程请参照前述方法的阐述。
所述可读存储介质为计算机可读存储介质,其可以是前述任一实施例所述的控制器的内部存储单元,例如控制器的硬盘或内存。所述可读存储介质也可以是所述控制器的外部存储设备,例如所述控制器上配备的插接式硬盘,智能存储卡(SmartMediaCard,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述可读存储介质还可以既包括所述控制器的内部存储单元也包括外部存储设备。所述可读存储介质用于存储所述计算机程序以及所述控制器所需的其他程序和数据。所述可读存储介质还可以用于暂时地存储已经输出或者将要输出的数据。
基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分,或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的可读存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
为了验证本发明所述目标定位方法,进行如下仿真实验来验证方法的可行性:
仿真实验利用运动平台定位相对静止目标,运动平台搭载磁梯度张量系统。设磁性目标位于空间直角坐标系的原点。设置磁性目标的磁矩m为(5×107,1×106,1×106)A·m2,运动平台的起始坐标为(100m,-100m,100m),距离矢量d=(6m,8m,0m),平台运动速度为50m/s,构造指数为3;仿真实验中磁力仪的测量精度为0.01nT。其定位误差结果如图3。
从图3可以看出,在测量精度为0.01nT的条件下,本发明的目标定位算法的定位误差没有超过2.5m,说明本本发明的目标定位算法在该精度条件下可以准确定位目标体的位置,并且不受背景场的影响。
此外,在消除背景场影响下,基于等式进行的三点定位误差整体上要大于基于等式进行的两点定位误差。引起这一误差的一个原因是由于磁梯度张量Gi为磁场矢量Ti的一阶空间导数,求导使得噪声放大,导致磁梯度张量Gi受噪声的影响大于磁场矢量Ti。
需要强调的是,本发明所述的实例是说明性的,而不是限定性的,因此本发明不限于具体实施方式中所述的实例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,不脱离本发明宗旨和范围的,不论是修改还是替换,同样属于本发明的保护范围。
Claims (10)
1.一种基于磁场梯度张量的目标定位方法,其特征在于:包括以下步骤:
步骤S1:获取观测站记录的观测数据;
步骤S2:基于步骤S1记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体之间的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置;
定位公式如下:
ri=[Gi-1-Gi+1]-1·[n(Ti+1-Ti-1)-Gi+1·d]
或
ri=[Gi-1-Gi+1]-1·[nGi-Gi+1]·d
其中,ri表示在T时刻所述观测点P与所述目标体之间的距离矢量,Ti-1、Ti+1为观测点P的两个相邻观测位置i-1、i+1在同一时刻T记录的磁场矢量或者为观测点P在T时刻的两个相邻时刻i-1、i+1记录的磁场矢量;Gi以及Gi-1、Gi+1依次为所述观测点P的对应观测位置i及两个相邻观测位置i-1、i+1在同一时刻T记录的磁场梯度张量或者为所述观测点P在T时刻及两个相邻时刻i-1、i+1记录的磁场梯度张量;
d=ri+1-ri-1表示两个观测位置i+1、i-1或两个时刻i+1、i-1对应的观测点与目标体之间的距离矢量的差值,n为目标体的构造指数。
2.根据权利要求1所述的目标定位方法,其特征在于:在任意时刻任意观测点与目标体的距离矢量表示为:
R=[xs-x* ys-y* zs-z*]
式中,x*、y*、z*分别表示观测点在xyz轴上的坐标,xs、ys、zs分别表示目标体在xyz轴上的坐标。
3.根据权利要求1所述的目标定位方法,其特征在于:所述目标体的构造指数n的取值范围为[1,3]。
4.根据权利要求3所述的目标定位方法,其特征在于:所述目标体为球体、或长方体、或二度体时,所述构造指数的取值对应为3,2,1。
5.根据权利要求1所述的目标定位方法,其特征在于:所述目标体位于太空、或地表、或地表以上、或地表以下、或水面、或水面以上、或水面以下。
6.根据权利要求1所述的目标定位方法,其特征在于:所述观测站为固定基站或者位于运动平台上;所述目标体为静止目标体或运动目标体。
7.一种基于权利要求1-6任一项所述的目标定位方法的磁场梯度张量测量装置,其特征在于:所述磁场梯度张量测量装置包括8个磁传感器,所述8个磁传感器分别分布在正方体的8个顶角;
其中,所述正方体上非相邻的两个侧面的中心点及其正方体的中心点构成三个相邻观测位置时,磁传感器F1、F2、F3和F4组成一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi-1;磁传感器F5、F6、F7和F8组成另一个正方形磁场梯度测量系统,用于测定磁场梯度张量Gi+1;磁传感器F1-F8组成立方体磁场梯度测量系统,用于测定磁场梯度张量Gi。
8.一种基于权利要求1-6任一项所述的目标定位方法的装置,其特征在于:包括:
观测模块,用于记录观测数据;
定位模块,用于基于记录的所述观测数据,并利用定位公式计算出在T时刻观测点P与目标体的距离矢量,再基于所述距离矢量以及在T时刻的所述观测点P的观测位置确定在T时刻所述目标体的位置。
9.一种电子终端,其特征在于:包括:一个或多个处理器和存储了一个或多个计算机程序的存储器;
所述计算机程序被处理器调用以实现:
权利要求1-6任一项所述目标定位方法的步骤。
10.一种可读存储介质,其特征在于:存储了计算机程序,所述计算机程序被处理器调用以实现:
权利要求1-6任一项所述目标定位方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210433815.3A CN114777766B (zh) | 2022-04-24 | 2022-04-24 | 一种基于磁场梯度张量的目标定位方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210433815.3A CN114777766B (zh) | 2022-04-24 | 2022-04-24 | 一种基于磁场梯度张量的目标定位方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114777766A true CN114777766A (zh) | 2022-07-22 |
CN114777766B CN114777766B (zh) | 2023-10-31 |
Family
ID=82432769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210433815.3A Active CN114777766B (zh) | 2022-04-24 | 2022-04-24 | 一种基于磁场梯度张量的目标定位方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114777766B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170074660A1 (en) * | 2015-09-16 | 2017-03-16 | Raytheon Company | Magnetic Field Gradient Navigation Aid |
CN106990440A (zh) * | 2017-02-22 | 2017-07-28 | 清华大学 | 一种基于两个探测位置处磁场空间梯度的潜艇定位方法 |
CN111522835A (zh) * | 2020-04-02 | 2020-08-11 | 天津大学 | 一种基于数据库特征匹配的多磁性目标位置探测方法 |
CN111721301A (zh) * | 2020-07-13 | 2020-09-29 | 中南大学 | 一种基于重力矢量及其梯度的差分定位方法与装置 |
CN112684511A (zh) * | 2021-01-06 | 2021-04-20 | 中国人民解放军海军工程大学 | 基于两点磁梯度全张量的线性定位方法 |
-
2022
- 2022-04-24 CN CN202210433815.3A patent/CN114777766B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170074660A1 (en) * | 2015-09-16 | 2017-03-16 | Raytheon Company | Magnetic Field Gradient Navigation Aid |
CN106990440A (zh) * | 2017-02-22 | 2017-07-28 | 清华大学 | 一种基于两个探测位置处磁场空间梯度的潜艇定位方法 |
CN111522835A (zh) * | 2020-04-02 | 2020-08-11 | 天津大学 | 一种基于数据库特征匹配的多磁性目标位置探测方法 |
CN111721301A (zh) * | 2020-07-13 | 2020-09-29 | 中南大学 | 一种基于重力矢量及其梯度的差分定位方法与装置 |
CN112684511A (zh) * | 2021-01-06 | 2021-04-20 | 中国人民解放军海军工程大学 | 基于两点磁梯度全张量的线性定位方法 |
Non-Patent Citations (2)
Title |
---|
常帅 等: "基于磁信标的水下SLAM 方法", 水下无人系统学报, vol. 27, no. 3 * |
张光 等: "载体平动条件下的磁梯度张量定位方法", 华中科技大学学报(自然科学版), vol. 41, no. 1 * |
Also Published As
Publication number | Publication date |
---|---|
CN114777766B (zh) | 2023-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6999021B2 (en) | Method and apparatus for detecting, mapping and locating underground utilities | |
CN108931241B (zh) | 基于磁总场梯度和磁张量模量梯度的磁目标单点定位方法 | |
US9355453B2 (en) | Three-dimensional measurement apparatus, model generation apparatus, processing method thereof, and non-transitory computer-readable storage medium | |
US9703001B2 (en) | Method for automatic recognition of a mobile magnetic object | |
CN109425906B (zh) | 一种磁异常探测矢量磁目标识别方法 | |
CN103874904B (zh) | 偏移估计装置、偏移估计方法、偏移估计程序以及信息处理装置 | |
CN106643792A (zh) | 惯性测量单元和地磁传感器整体标定装置及标定方法 | |
KR20180131033A (ko) | 카메라와 레이더의 캘리브레이션 장치 및 방법 | |
CN111458573A (zh) | 导线电力参数测量方法、装置、系统和计算机设备 | |
US11946986B2 (en) | Magnetic detector, detection method, and non-transitory computer readable storage medium | |
JP2014149211A (ja) | オフセット推定装置及びプログラム | |
CN104848861A (zh) | 一种基于图像消失点识别技术的移动设备姿态测量方法 | |
US20160370441A1 (en) | Magnetic field imaging system | |
WO2019013673A1 (ru) | Магнитный дефектоскоп для диагностики подземных стальных трубопроводов | |
CN111650539A (zh) | 导线运行参数检测方法、系统、装置和计算机设备 | |
CN109633540B (zh) | 一种磁源的实时定位系统及实时定位方法 | |
Zhao et al. | A detection system for pipeline direction based on shielded geomagnetic field | |
Wang et al. | From model to algorithms: Distributed magnetic sensor system for vehicle tracking | |
CN114449439B (zh) | 一种地下管廊空间的定位方法及装置 | |
CN109633541B (zh) | 一种磁源定位装置及磁源定位方法 | |
Bikmaev et al. | Improving the accuracy of supporting mobile objects with the use of the algorithm of complex processing of signals with a monocular camera and LiDAR | |
CN109633539B (zh) | 一种磁源的静态定位装置及静态定位方法 | |
CN114777766A (zh) | 一种基于磁场梯度张量的目标定位方法及装置 | |
CN112344925A (zh) | 基于归一化磁矩的水下目标相关性检测定位方法 | |
CN114266835A (zh) | 一种非量测相机的变形监测控制方法与系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |