CN102621584A - 航磁矢量数据处理方法 - Google Patents

Abstract

本发明涉及一种航磁矢量数据的处理方法。磁矢量数据处理流程为：惯导传感器与磁力仪装置连接，上位机接收磁力仪和惯性导航同一时刻采集的数据→磁力仪数据与惯性导航数据同步处理并保存→磁力仪数据不正交度校正→磁力仪数据转向差校正→磁矢量姿态变换→作磁矢量曲线图。通过对磁力仪探头采集精度影响的因素进行单独的校正，有针对性的消除对数据采集精度的影响；经二次校正后的磁力仪精度得到较大提高；通过对磁力仪数据包与惯性导航数据包进行数据同步处理，确保二者采集信息匹配；通过惯性导航姿态信息对磁矢量进行姿态变换，使磁场数据转换为基于固定地理坐标系下的磁矢量，为航磁测量提供可靠的矢量数据。

Description

航磁矢量数据处理方法

技术领域

本发明涉及一种地球物理磁法勘探数据的处理，尤其是三分量磁通门磁力仪并配以惯性导航系统与GPS全球定位系统进行航空磁矢量数据的处理方法。

背景技术

航空磁矢量数据测量系统中，固定在飞机上的磁力仪传感器实时测量空间矢量磁场三分量值，磁力仪传感器基于载体坐标系；同时惯性导航系统实时记录飞行器的姿态信息(翻滚角、俯仰角、方位角)，系统基于地理坐标系，且惯性导航系统与磁力仪传感器刚性固定，磁力仪数据与惯性导航数据由GPS秒脉冲同步。通过坐标体系之间的坐标变换，将基于实时改变的载体坐标系下的磁场测量值转化为基于固定地理坐标系下的磁场矢量值，以便于数据分析。

现有的磁矢量数据处理方法中，主要的是对磁力仪传感器测量误差进行单方面的校正。探测与控制学报，2005，27(2)，林春生等《三轴磁强计正交误差分析与校正》，对三轴磁强计中由于三轴间不正交和灵敏度不一致所引起的测量误差进行了详细分析和理论计算，提出了三轴磁强计的正交变换，灵敏度一致性校正及消除零点漂移的方法并进行了实验验证。其仅仅对磁力仪传感器进行了数据处理，未涉及到磁力仪数据与惯性导航数据的同步，也未涉及到坐标体系间的姿态变换。

海军工程大学学报，2005，17(6)，闫辉等《一种海洋环境地磁场三分量的测量方法》，提出了一种利用安装在舰船上的三分量磁传感器进行地磁场测量的方法。从理论上解决了在任意舰船姿态下，消除测量中舰船磁场干扰的问题。其首先通过舰船上的三分量磁传感器测量环境中的磁场值，然后同步采集舰船姿态信号和航向信号，接着将舰船磁场三分量值变换成地磁场值。其不足之处在于仅仅进行了理论分析仿真而未现场试验，对航磁三分量磁传感器的三轴间不正交和灵敏度不一致所引起的测量误差不能有效的校正，而且该方法是用于海洋环境地磁场数据处理，尚不能适用于航空磁场数据处理。

发明内容：

本发明的目的就在于针对上述现有技术的不足，提供一种磁矢量数据处理方法。

本发明的目的是通过以下技术方案实现的：

磁矢量数据处理方法，磁矢量数据处理流程为：首将惯导传感器与磁力仪装置连接，其次磁力仪探头三通道灵敏度标定、磁力仪探头零点漂移校正，然后是磁力仪装置数据采集，上位机接收磁力仪和惯导传感器同一时刻采集的数据→磁力仪数据与惯导传感器数据同步处理并保存→磁力仪数据不正交度校正→磁力仪数据转向差校正→磁矢量姿态变换→作磁矢量曲线图。

磁力仪数据与惯性导航数据同步处理并保存，包括以下顺序和步骤：

步骤1、磁力仪数据包接收和惯性导航数据包接收，磁力仪装置每一秒钟上传一个数据包，该数据包包含20组采集数据，每组包含X、Y、Z三个分量的数据；惯性导航装置每50ms上传一个数据包，该数据包包含经度、纬度、时间和姿态偏转角的采集数据；

步骤2、上位机软件分别将当前接收的磁力仪数据包和惯性导航数据包进行解压，并以磁力仪数据包中的采集时间为基准，匹配惯性导航数据包中的时间；

步骤3、若时间信息一致且惯性导航数据包为该秒的第一个数据包，则将磁力仪数据与惯性导航数据重新打包为一个新的数据包；若时间信息一致且惯性导航数据包不是该秒的第一个数据包，则等待并接收下一组磁力仪数据包，并跳转到步骤2；

步骤4、判断该秒惯性导航数据包数目是否等于20，若惯性导航数据包数目＜20，则等待新的导航数据包上传，解压数据惯性导航数据，并将数据打包到步骤3中的数据包之后；

步骤5、重复步骤4，直到该秒的20组惯性导航数据全部接收并打包完毕，然后将完整的数据包保存到文件；

步骤6、若两者时间信息不一致，则通过判断两者时间超前/滞后量，来移动软件当前处理的数据包；

步骤7、若当前的磁力仪数据时间滞后于惯性导航数据时间，则以磁力仪数据时间为基准，等待并接收下一组惯性导航数据包；

步骤8、判断接收的惯性导航数据包时间，并重复步骤7，直到二者时间信息一致，并执行步骤3；

步骤9、若当前的磁力仪数据时间超前于惯性导航数据时间，则等待并接收下一组磁力仪数据包；

步骤10、判断二者的时间信息，并重复步骤9，直到二者时间信息一致，执行步骤3。

传感器三轴不正交度校正的步骤如下：

步骤1、在一个磁场均匀恒定的空间，转动探头，依次记录下等概率的在各

个方向上的分量值B_i′＝(B_xi′，B_yi′，B_zi′)^T，i从0～N，N为总数据点数且N＞＞3；步骤2、寻找参数α、β、γ使得目标函数

取得最小值；其中，

B_i″＝Tinv×B_i′为折算到标准正交坐标系中的磁场，T_inv＝T^-1，

$T=\left[\begin{array}{ccc}1& 0& 0\\ \sin \left(\mathrm{\α}\right)& \cos \left(\mathrm{\α}\right)& 0\\ \cos \left(\mathrm{\γ}\right)\sin \left(\mathrm{\β}\right)& \cos \left(\mathrm{\γ}\right)\cos \left(\mathrm{\β}\right)& \sin \left(\mathrm{\γ}\right)\end{array}\right];$

步骤3、重复步骤1与步骤2进行多次独立实验并比较多次实验的反演结果。

磁力仪数据转向差校正：

步骤1：在正交坐标系OX1Y1Z1与实际坐标系OXYZ转换中，先将0Z与0Z1重合，再将YOZ与Y1OZ1共面，其中0X与X1OY1的夹角为c，与X1OZ1的夹角为a，OY与OY1夹角为b。则实际坐标系与正交坐标系之间的关系矩阵为A1＝[cos(a)*cos(c)sin(c)sin(a)；0cos(b)sin(b)；0 0 1]；令B1为实际磁分量，B为正交磁分量，则B1＝A1*B；

步骤2：对于磁力仪三轴灵敏度校正，以OZ1轴灵敏度为标准，设OX1和OY1两轴灵敏度分别为Kx和Ky。则灵敏度关系矩阵为A2＝[Kx 0 0；0 Ky 0；0 0 Kz]；令B2为实际磁分量，则B2＝A2*B1；

步骤3：零点漂移校正，设磁力仪三轴的零点漂移分别为Bx0、By0、Bz0，令B3为实际磁场三分量输出，则B3＝B2+B0，其中B0＝[Bx0，By0，Bz0]；

步骤4：结合上述三个误差分析，得到磁力仪实际输出B3与理想输出B的关系为B3＝f(B)＝A2*A1*B+B0＝C*B+B0；令磁力仪误差校正模型为B4＝g(B3)，若g(B3)为f(B)的逆映射，就能使输出B4在数值上与B一致，达到校正的目的；

步骤5：在处于稳定的磁场环境中，磁力仪的空间姿态会不断改变，但其输出的总磁场值不变，即B3的模值不变，则校正模型的目标转变为：使校正结果的模值与理想值一致；令Bt(n)＝[Btx(n)，Bty(n)，Btz(n)]为实际测量的磁场数据；

n为测量点数，则校正模型经过整理得到：

B²-Bt(n)²＝2*(c11*Btx(n)²+c22*Bty(n)²+c12*Btx(n)*Bty(n)+c13*Btx(n)*Btz(n)+c23*Bty(n)*Btz(n)-Bx0*Btx(n)-By0*Bty(n)-Bz0*Btz(n))；          (1)

其中c11、c22、c12、c13、c23、Bx0、By0、Bz0为各项的系数，根据实际测量的数据，求解出上述系数，即实现校正；令

W(n)＝[Btx(n)²，Bty(n)²，Btx(n)*Bty(n)，Btx(n)*Btz(n)，Bty(n)*Btz(n)，-Btx(n)，-Bty(n)，-Btz(n)]；                                      (2)

权向量w(n)＝[c11，c22，c12，c13，c23，Bx0，By0，Bz0]；       (3)

则2*S(n)＝B2-Bt(n)2＝2*W(n)*w(n)T；                          (4)

经过函数B²-2*S(n)后得到输出Bt1(n)²，更新算法可修改权向量w，减少输出Bt1(n)²与实际测量磁场Bt(n)²之差e(n)，以逐渐优化权向量w，

其中Bt1(n)²＝B²-2*W(n)*w(n)T，                (5)

式中w(n)表示权向量w第n步迭代结果；

偏差表达式：e(n)＝Bt1(n)²-Bt(n)²；            (6)

参数调整表达式w(n+1)＝w(n)-d*e(n)*W(n)，      (7)

式中d为迭代参数，本测试预制为0.00001，测试中反复迭代至偏差e(n)不再减小，此时得到的权值向量w即为所求系数，将所求系数代入到校正矩阵公式(1)中，即对磁力仪转向差进行校正。

磁矢量姿态变换：

步骤1：提取磁力仪装置校正后的磁场矢量数据：X分量、Y分量、Z分量，该矢量数据为磁力仪传感器测量的矢量空间某一位置地球磁场的各个分量值，基于载体坐标系；

步骤2：提取惯性导航姿态数据：翻滚角(Roll)、俯仰角(Pitch)、方位角(Heading)，该姿态数据为惯导传感器旋转时的欧拉角度，基准坐标系为地理坐标系；惯导传感器与磁力仪传感器刚性固定在无磁性平台上，此时翻滚角(Roll)、俯仰角(Pitch)和方位角(Heading)同样是磁力仪传感器旋转时的欧拉角度；

步骤3：令导航系统姿态角度分别是：方位角JH，俯仰角JP，翻滚角JR，则地理坐标系到载体坐标系的转换矩阵为：

$S=\left[\begin{array}{ccc}(\cos \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right))& (-\sin \left(\mathrm{JP}\right))\\ (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)-\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)+\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\\ (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)+\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)-\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\end{array}\right]$

通过地理坐标系与载体坐标系之间的姿态变换矩阵，将测量的空间任意姿态磁场三分量数据(X分量、Y分量、Z分量)投影到地理坐标系下得到磁场矢量数据，依该数据绘出矢量曲线图。

有益效果：本发明通过对磁力仪探头采集精度影响的因素进行单独的校正，有针对性的消除特定影响因素对数据采集精度的影响；经二次校正后的磁力仪精度得到较大提高；通过对磁力仪数据包与惯性导航数据包进行数据同步处理，确保二者采集信息匹配；通过惯性导航姿态信息对磁矢量进行姿态变换，使磁场数据转换为基于固定地理坐标系下的磁矢量，为航磁测量提供可靠的矢量数据。

附图说明：

附图1磁矢量数据处理方法流程图

附图2磁力仪数据与惯性导航数据同步并保存处理流程图

附图3不正交度校正中两坐标系示意图

附图4惯导传感器与磁力仪探头安装图及坐标系转换图

具体实施方式：

下面结合附图和实施例作进一步的详细说明：

磁矢量数据处理方法，磁矢量数据处理流程为：首将惯导传感器与磁力仪装置连接，其次磁力仪探头三通道灵敏度标定、磁力仪探头零点漂移校正，然后是磁力仪装置数据采集，上位机接收磁力仪和惯性导航同一时刻采集的数据→磁力仪数据与惯性导航数据同步处理并保存→磁力仪数据不正交度校正→磁力仪数据转向差校正→磁矢量姿态变换→作曲线图。

磁力仪数据与惯性导航数据同步处理并保存，包括以下顺序和步骤：

步骤1、磁力仪数据包接收和惯性导航数据包接收，磁力仪装置每一秒钟上传一个数据包，该数据包包含20组采集数据，每组包含X、Y、Z三个分量的数据；惯性导航装置每50ms上传一个数据包，该数据包包含经度、纬度、时间和姿态偏转角的采集数据；

步骤2、上位机软件分别将当前接收的磁力仪数据包和惯性导航数据包进行解压，并以磁力仪数据包中的采集时间为基准，匹配惯性导航数据包中的时间；

步骤3、若时间信息一致且惯性导航数据包为该秒的第一个数据包，则将磁力仪数据与惯性导航数据重新打包为一个新的数据包；若时间信息一致且惯性导航数据包不是该秒的第一个数据包，则等待并接收下一组磁力仪数据包，并跳转到步骤2；

步骤4、判断该秒惯性导航数据包数目是否等于20，若惯性导航数据包数目＜20，则等待新的导航数据包上传，解压数据惯性导航数据，并将数据打包到步骤3中的数据包之后；

步骤5、重复步骤4，直到该秒的20组惯性导航数据全部接收并打包完毕，然后将完整的数据包保存到文件；

步骤6、若两者时间信息不一致，则通过判断两者时间超前/滞后量，来移动软件当前处理的数据包；

步骤7、若当前的磁力仪数据时间滞后于惯性导航数据时间，则以磁力仪数据时间为基准，等待并接收下一组惯性导航数据包；

步骤8、判断接收的惯性导航数据包时间，并重复步骤7，直到二者时间信息一致，并执行步骤3；

步骤9、若当前的磁力仪数据时间超前于惯性导航数据时间，则等待并接收下一组磁力仪数据包；

步骤10、判断二者的时间信息，并重复步骤9，直到二者时间信息一致，执行步骤3。

传感器三轴不正交度校正的步骤如下：

步骤1、在一个磁场均匀恒定的空间，转动探头，依次记录下等概率的在各

个方向上的分量值B_i′＝(B_xi′，B_yi′，B_zi′)^T，i从0～N，N为总数据点数且N＞＞3；

步骤2、寻找参数α、β、γ使得目标函数

取得最小值；其中，

B_i″＝Tinv×B_i′为折算到标准正交坐标系中的磁场，T_inv＝T^-1，

$T=\left[\begin{array}{ccc}1& 0& 0\\ \sin \left(\mathrm{\α}\right)& \cos \left(\mathrm{\α}\right)& 0\\ \cos \left(\mathrm{\γ}\right)\sin \left(\mathrm{\β}\right)& \cos \left(\mathrm{\γ}\right)\cos \left(\mathrm{\β}\right)& \sin \left(\mathrm{\γ}\right)\end{array}\right];$

步骤3、重复步骤1与步骤2进行多次独立实验并比较多次实验的反演结果。

磁力仪数据转向差校正：

步骤1：在正交坐标系OX1Y1Z1与实际坐标系OXYZ转换中，先将0Z与0Z1重合，再将YOZ与Y1OZ1共面，其中0X与X1OY1的夹角为c，与X1OZ1的夹角为a，OY与OY1夹角为b。则实际坐标系与正交坐标系之间的关系矩阵为A1＝[cos(a)*cos(c)sin(c)sin(a)；0 cos(b)sin(b)；0 0 1]；令B1为实际磁分量，B为正交磁分量，则B1＝A1*B；

步骤2：对于磁力仪三轴灵敏度校正，以OZ1轴灵敏度为标准，设OX1和OY1两轴灵敏度分别为Kx和Ky。则灵敏度关系矩阵为A2＝[Kx 0 0；0 Ky 0；0 0 Kz]；令B2为实际磁分量，则B2＝A2*B1；

步骤3：零点漂移校正，设磁力仪三轴的零点漂移分别为Bx0、By0、Bz0，令B3为实际磁场三分量输出，则B3＝B2+B0，其中B0＝[Bx0，By0，Bz0]；

步骤4：结合上述三个误差分析，得到磁力仪实际输出B3与理想输出B的关系为B3＝f(B)＝A2*A1*B+B0＝C*B+B0；令磁力仪误差校正模型为B4＝g(B3)，若g(B3)为f(B)的逆映射，就能使输出B4在数值上与B一致，达到校正的目的；

步骤5：在处于稳定的磁场环境中，磁力仪的空间姿态会不断改变，但其输出的总磁场值不变，即B3的模值不变，则校正模型的目标转变为：使校正结果的模值与理想值一致；令Bt(n)＝[Btx(n)，Bty(n)，Btz(n)]为实际测量的磁场数据；

n为测量点数，则校正模型经过整理得到：

B²-Bt(n)²＝2*(c11*Btx(n)²+c22*Bty(n)²+c12*Btx(n)*Bty(n)+c13*Btx(n)*Btz(n)+c23*Bty(n)*Btz(n)-Bx0*Btx(n)-By0*Bty(n)-Bz0*Btz(n))；(1)

其中c11、c22、c12、c13、c23、Bx0、By0、Bz0为各项的系数，根据实际测量的数据，求解出上述系数，即实现校正；令

W(n)＝[Btx(n)²，Bty(n)²，Btx(n)*Bty(n)，Btx(n)*Btz(n)，Bty(n)*Btz(n)，-Btx(n)，-Bty(n)，-Btz(n)]；                            (2)

权向量w(n)＝[c11，c22，c12，c13，c23，Bx0，By0，Bz0]；(3)

则2*S(n)＝B2-Bt(n)2＝2*W(n)*w(n)T；            (4)

经过函数B²-2*S(n)后得到输出Bt1(n)²，更新算法可修改权向量w，减少输出Bt1(n)²与实际测量磁场Bt(n)²之差e(n)，以逐渐优化权向量w，

其中Bt1(n)²＝B²-2*W(n)*w(n)T，                 (5)

式中w(n)表示权向量w第n步迭代结果；

偏差表达式：e(n)＝Bt1(n)²-Bt(n)²；             (6)

参数调整表达式w(n+1)＝w(n)-d*e(n)*W(n)，       (7)

式中d为迭代参数，本测试预制为0.00001，测试中反复迭代至偏差e(n)不再减小，此时得到的权值向量w即为所求系数，将所求系数代入到校正矩阵公式(1)中，即对磁力仪转向差进行校正。

磁矢量姿态变换：

步骤1：提取磁力仪装置校正后的磁场矢量数据：X分量、Y分量、Z分量，该矢量数据为磁力仪传感器测量的矢量空间某一位置地球磁场的各个分量值，基于载体坐标系；

步骤2：提取惯性导航姿态数据：翻滚角(Roll)、俯仰角(Pitch)、方位角(Heading)，该姿态数据为惯导传感器探头旋转时的欧拉角度，基准坐标系为地理坐标系；惯导传感器与磁力仪传感器刚性固定在无磁性平台上，此时翻滚角(Roll)、俯仰角(Pitch)和方位角(Heading)同样是磁力仪传感器旋转时的欧拉角度；

步骤3：令导航系统姿态角度分别是：方位角JH，俯仰角JP，翻滚角JR，则地理坐标系到载体坐标系的转换矩阵为：

$S=\left[\begin{array}{ccc}(\cos \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right))& (-\sin \left(\mathrm{JP}\right))\\ (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)-\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)+\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\\ (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)+\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)-\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\end{array}\right]$

通过地理坐标系与载体坐标系之间的姿态变换矩阵，将测量的空间任意姿态磁场三分量数据(X分量、Y分量、Z分量)投影到地理坐标系下得到磁场矢量数据，依该数据绘出矢量曲线图。

Claims (5)

1.一种航磁矢量数据处理方法，其特征在于，磁矢量数据处理流程为：首将惯导传感器与磁力仪装置连接，其次磁力仪探头三通道灵敏度标定、磁力仪探头零点漂移校正，然后是磁力仪装置数据采集，上位机接收磁力仪和惯导传感器同一时刻采集的数据→磁力仪数据与惯性导航数据同步处理并保存→磁力仪数据不正交度校正→磁力仪数据转向差校正→磁矢量姿态变换→作曲线图。

2.按照权利要求1所述的航磁矢量数据处理方法，其特征在于，磁力仪数据与惯性导航数据同步处理并保存，包括以下顺序和步骤：

步骤1、磁力仪数据包接收和惯性导航数据包接收，磁力仪装置每一秒钟上传一个数据包，该数据包包含20组采集数据，每组包含X、Y、Z三个分量的数据；惯性导航装置每50ms上传一个数据包，该数据包包含经度、纬度、时间和姿态偏转角的采集数据；

步骤2、上位机软件分别将当前接收的磁力仪数据包和惯性导航数据包进行解压，并以磁力仪数据包中的采集时间为基准，匹配惯性导航数据包中的时间；

步骤3、若时间信息一致且惯性导航数据包为该秒的第一个数据包，则将磁力仪数据与惯性导航数据重新打包为一个新的数据包；若时间信息一致且惯性导航数据包不是该秒的第一个数据包，则等待并接收下一组磁力仪数据包，并跳转到步骤2；

步骤4、判断该秒惯性导航数据包数目是否等于20，若惯性导航数据包数目＜20，则等待新的导航数据包上传，解压数据惯性导航数据，并将数据打包到步骤3中的数据包之后；

步骤5、重复步骤4，直到该秒的20组惯性导航数据全部接收并打包完毕，然后将完整的数据包保存到文件；

步骤6、若两者时间信息不一致，则通过判断两者时间超前/滞后量，来移动软件当前处理的数据包；

步骤7、若当前的磁力仪数据时间滞后于惯性导航数据时间，则以磁力仪数据时间为基准，等待并接收下一组惯性导航数据包；

步骤8、判断接收的惯性导航数据包时间，并重复步骤7，直到二者时间信息一致，并执行步骤3；

步骤9、若当前的磁力仪数据时间超前于惯性导航数据时间，则等待并接收下一组磁力仪数据包；

步骤10、判断二者的时间信息，并重复步骤9，直到二者时间信息一致，执行步骤3。

3.按照权利要求1所述的航磁矢量数据处理方法，其特征在于，传感器三轴不正交度校正的步骤如下：

步骤1、在一个磁场均匀恒定的空间，转动探头，依次记录下等概率的在各

个方向上的分量值B_i′＝(B_xi′，B_yi′，B_zi′)^T，i从0～N，N为总数据点数且N＞＞3；

步骤2、寻找参数α、β、γ使得目标函数

取得最小值；其中，B_i″＝Tinv×B_i′为折算到标准正交坐标系中的磁场，T_inv＝T^-1，

$T=\left[\begin{array}{ccc}1& 0& 0\\ \sin \left(\mathrm{\α}\right)& \cos \left(\mathrm{\α}\right)& 0\\ \cos \left(\mathrm{\γ}\right)\sin \left(\mathrm{\β}\right)& \cos \left(\mathrm{\γ}\right)\cos \left(\mathrm{\β}\right)& \sin \left(\mathrm{\γ}\right)\end{array}\right];$

步骤3、重复步骤1与步骤2进行多次独立实验并比较多次实验的反演结果。

4.按照权利要求1所述的航磁矢量数据处理方法，其特征在于，磁力仪数据转向差校正：

步骤1：在正交坐标系OX1Y1Z1与实际坐标系OXYZ转换中，先将0Z与0Z1重合，再将YOZ与Y1OZ1共面，其中0X与X1OY1的夹角为c，与X1OZ1的夹角为a，OY与OY1夹角为b。则实际坐标系与正交坐标系之间的关系矩阵为A1＝[cos(a)*cos(c)sin(c)sin(a)；0cos(b)sin(b)；001]；令B1为实际磁分量，B为正交磁分量，则B1＝A1*B；

步骤2：对于磁力仪三轴灵敏度校正，以OZ1轴灵敏度为标准，设OX1和OY1两轴灵敏度分别为Kx和Ky。则灵敏度关系矩阵为A2＝[Kx 0 0；0 Ky 0；0 0 Kz]；令B2为实际磁分量，则B2＝A2*B1；

步骤3：零点漂移校正，设磁力仪三轴的零点漂移分别为Bx0、By0、Bz0，令B3为实际磁场三分量输出，则B3＝B2+B0，其中B0＝[Bx0，By0，Bz0]；

步骤4：结合上述三个误差分析，得到磁力仪实际输出B3与理想输出B的关系为B3＝f(B)＝A2*A1*B+B0＝C*B+B0；令磁力仪误差校正模型为B4＝g(B3)，若g(B3)为f(B)的逆映射，就能使输出B4在数值上与B一致，达到校正的目的；

步骤5：在处于稳定的磁场环境中，磁力仪的空间姿态会不断改变，但其输出的总磁场值不变，即B3的模值不变，则校正模型的目标转变为：使校正结果的模值与理想值一致；令Bt(n)＝[Btx(n)，Bty(n)，Btz(n)]为实际测量的磁场数据；

n为测量点数，则校正模型经过整理得到：

B²-Bt(n)²＝2*(c11*Btx(n)²+c22*Bty(n)²+c12*Btx(n)*Bty(n)+c13*Btx(n)*Btz(n)+c23*Bty(n)*Btz(n)-Bx0*Btx(n)-By0*Bty(n)-Bz0*Btz(n))；(1)

其中c11、c22、c12、c13、c23、Bx0、By0、Bz0为各项的系数，根据实际测量的数据，求解出上述系数，即实现校正；令

W(n)＝[Btx(n)²，Bty(n)²，Btx(n)*Bty(n)，Btx(n)*Btz(n)，Bty(n)*Btz(n)，-Btx(n)，-Bty(n)，-Btz(n)]；                               (2)

权向量w(n)＝[c11，c22，c12，c13，c23，Bx0，By0，Bz0]；(3)

则2*S(n)＝B2-Bt(n)2＝2*W(n)*w(n)T；                   (4)

经过函数B²-2*S(n)后得到输出Bt1(n)²，更新算法可修改权向量w，减少输出Bt1(n)²与实际测量磁场Bt(n)²之差e(n)，以逐渐优化权向量w，

其中Bt1(n)²＝B²-2*W(n)*w(n)T，                (5)

式中w(n)表示权向量w第n步迭代结果；

偏差表达式：e(n)＝Bt1(n)²-Bt(n)²；            (6)

参数调整表达式w(n+1)＝w(n)-d*e(n)*W(n)，      (7)

式中d为迭代参数，本测试预制为0.00001，测试中反复迭代至偏差e(n)不再减小，此时得到的权值向量w即为所求系数，将所求系数代入到校正矩阵公式(1)中，即对磁力仪转向差进行校正。

5.按照权利要求1所述的磁矢量数据处理方法，其特征在于，磁矢量姿态变换：

步骤1：提取磁力仪装置校正后的磁场矢量数据：X分量、Y分量、Z分量，该矢量数据为磁力仪传感器测量的矢量空间某一位置地球磁场的各个分量值，基于载体坐标系；

步骤2：提取惯性导航姿态数据：翻滚角(Roll)、俯仰角(Pitch)、方位角(Heading)，该姿态数据为惯导传感器旋转时的欧拉角度，基准坐标系为地理坐标系；惯导传感器与磁力仪传感器刚性固定在无磁性平台上，此时翻滚角(Roll)、俯仰角(Pitch)和方位角(Heading)同样是磁力仪传感器旋转时的欧拉角度；

步骤3：令导航系统姿态角度分别是：方位角JH，俯仰角JP，翻滚角JR，则地理坐标系到载体坐标系的转换矩阵为：

$S=\left[\begin{array}{ccc}(\cos \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right))& (-\sin \left(\mathrm{JP}\right))\\ (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)-\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)+\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\\ (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\cos \left(\mathrm{JH}\right)+\sin \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\sin \left(\mathrm{JP}\right)\·\sin \left(\mathrm{JH}\right)-\sin \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JH}\right))& (\cos \left(\mathrm{JR}\right)\·\cos \left(\mathrm{JP}\right))\end{array}\right]$

通过地理坐标系与载体坐标系之间的姿态变换矩阵，将测量的空间任意姿态磁场三分量数据(X分量、Y分量、Z分量)投影到地理坐标系下得到磁场矢量数据，依该数据绘出矢量曲线图。

Images