CN102621584A - 航磁矢量数据处理方法 - Google Patents
航磁矢量数据处理方法 Download PDFInfo
- Publication number
- CN102621584A CN102621584A CN2011104160881A CN201110416088A CN102621584A CN 102621584 A CN102621584 A CN 102621584A CN 2011104160881 A CN2011104160881 A CN 2011104160881A CN 201110416088 A CN201110416088 A CN 201110416088A CN 102621584 A CN102621584 A CN 102621584A
- Authority
- CN
- China
- Prior art keywords
- data
- magnetometer
- inertial navigation
- sin
- cos
- 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
Images
Landscapes
- Measuring Magnetic Variables (AREA)
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、在一个磁场均匀恒定的空间,转动探头,依次记录下等概率的在各
Bi″=Tinv×Bi′为折算到标准正交坐标系中的磁场,Tinv=T-1,
步骤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为测量点数,则校正模型经过整理得到:
B2-Bt(n)2=2*(c11*Btx(n)2+c22*Bty(n)2+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)2,Bty(n)2,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)
经过函数B2-2*S(n)后得到输出Bt1(n)2,更新算法可修改权向量w,减少输出Bt1(n)2与实际测量磁场Bt(n)2之差e(n),以逐渐优化权向量w,
其中Bt1(n)2=B2-2*W(n)*w(n)T, (5)
式中w(n)表示权向量w第n步迭代结果;
偏差表达式:e(n)=Bt1(n)2-Bt(n)2; (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,则地理坐标系到载体坐标系的转换矩阵为:
通过地理坐标系与载体坐标系之间的姿态变换矩阵,将测量的空间任意姿态磁场三分量数据(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、在一个磁场均匀恒定的空间,转动探头,依次记录下等概率的在各
个方向上的分量值Bi′=(Bxi′,Byi′,Bzi′)T,i从0~N,N为总数据点数且N>>3;
Bi″=Tinv×Bi′为折算到标准正交坐标系中的磁场,Tinv=T-1,
步骤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为测量点数,则校正模型经过整理得到:
B2-Bt(n)2=2*(c11*Btx(n)2+c22*Bty(n)2+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)2,Bty(n)2,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)
经过函数B2-2*S(n)后得到输出Bt1(n)2,更新算法可修改权向量w,减少输出Bt1(n)2与实际测量磁场Bt(n)2之差e(n),以逐渐优化权向量w,
其中Bt1(n)2=B2-2*W(n)*w(n)T, (5)
式中w(n)表示权向量w第n步迭代结果;
偏差表达式:e(n)=Bt1(n)2-Bt(n)2; (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,则地理坐标系到载体坐标系的转换矩阵为:
通过地理坐标系与载体坐标系之间的姿态变换矩阵,将测量的空间任意姿态磁场三分量数据(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。
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为测量点数,则校正模型经过整理得到:
B2-Bt(n)2=2*(c11*Btx(n)2+c22*Bty(n)2+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)2,Bty(n)2,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)
经过函数B2-2*S(n)后得到输出Bt1(n)2,更新算法可修改权向量w,减少输出Bt1(n)2与实际测量磁场Bt(n)2之差e(n),以逐渐优化权向量w,
其中Bt1(n)2=B2-2*W(n)*w(n)T, (5)
式中w(n)表示权向量w第n步迭代结果;
偏差表达式:e(n)=Bt1(n)2-Bt(n)2; (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,则地理坐标系到载体坐标系的转换矩阵为:
通过地理坐标系与载体坐标系之间的姿态变换矩阵,将测量的空间任意姿态磁场三分量数据(X分量、Y分量、Z分量)投影到地理坐标系下得到磁场矢量数据,依该数据绘出矢量曲线图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110416088.1A CN102621584B (zh) | 2011-12-14 | 2011-12-14 | 航磁矢量数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110416088.1A CN102621584B (zh) | 2011-12-14 | 2011-12-14 | 航磁矢量数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102621584A true CN102621584A (zh) | 2012-08-01 |
CN102621584B CN102621584B (zh) | 2015-07-22 |
Family
ID=46561609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110416088.1A Expired - Fee Related CN102621584B (zh) | 2011-12-14 | 2011-12-14 | 航磁矢量数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102621584B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424780A (zh) * | 2013-08-27 | 2013-12-04 | 中国航空无线电电子研究所 | 基于感应线圈的载机磁环境补偿方法 |
CN103926626A (zh) * | 2014-04-30 | 2014-07-16 | 中色地科矿产勘查股份有限公司 | 航磁数据高度矫正方法及*** |
CN104199115A (zh) * | 2014-09-18 | 2014-12-10 | 中国人民解放军国防科学技术大学 | 一种基于线性模型的地磁矢量测量误差综合补偿方法 |
CN104535062A (zh) * | 2015-01-20 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 基于磁梯度张量和地磁矢量测量的运动式定位方法 |
CN104678340A (zh) * | 2015-02-27 | 2015-06-03 | 清华大学 | 一种磁强计测量误差纠正方法及*** |
CN105068131A (zh) * | 2015-08-03 | 2015-11-18 | 中国科学院电子学研究所 | 一种航磁数据调平方法 |
CN105572749A (zh) * | 2016-01-25 | 2016-05-11 | 中国地质大学(武汉) | 地面三分量磁力定向方法及地面三分量磁力定向勘探装置 |
CN104459828B (zh) * | 2014-12-08 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 基于绕轴法的地磁矢量***非对准校正方法 |
CN107024724A (zh) * | 2017-06-22 | 2017-08-08 | 吉林大学 | 反卷积快速优化直升机航空电磁探测早期数据的方法 |
CN107493531A (zh) * | 2017-08-04 | 2017-12-19 | 歌尔科技有限公司 | 一种头部姿态检测方法、装置和耳机 |
CN110308490A (zh) * | 2019-06-13 | 2019-10-08 | 南京理工大学 | 基于磁传感器阵列的旋转运动磁异常目标探测装置及方法 |
CN110417489A (zh) * | 2018-04-26 | 2019-11-05 | 杭州海康机器人技术有限公司 | 一种无人机定向天线的校准方法及装置 |
CN111257940A (zh) * | 2020-02-14 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
US11118913B2 (en) | 2016-10-19 | 2021-09-14 | Huawei Technologies Co., Ltd. | Vehicle positioning correction method and mobile device |
CN118050814A (zh) * | 2024-04-16 | 2024-05-17 | 山东省地质矿产勘查开发局第五地质大队(山东省第五地质矿产勘查院) | 低空光学定姿三分量磁测***及测定磁场矢量数据的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6427122B1 (en) * | 2000-12-23 | 2002-07-30 | American Gnc Corporation | Positioning and data integrating method and system thereof |
CN1644456A (zh) * | 2005-01-21 | 2005-07-27 | 清华大学 | 用于微小卫星的微型组合姿态测量*** |
CN101556155A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 小卫星姿态确定***及方法 |
CN101644579A (zh) * | 2009-08-17 | 2010-02-10 | 中国人民解放军国防科学技术大学 | 一种空中微小抛撒物运动参数测量方法及*** |
US20100312509A1 (en) * | 2009-06-05 | 2010-12-09 | Apple Inc. | Calibration techniques for an electronic compass in a portable device |
WO2011152141A1 (ja) * | 2010-05-31 | 2011-12-08 | オリンパスメディカルシステムズ株式会社 | 内視鏡形状検出装置及び内視鏡の挿入部の形状検出方法 |
-
2011
- 2011-12-14 CN CN201110416088.1A patent/CN102621584B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6427122B1 (en) * | 2000-12-23 | 2002-07-30 | American Gnc Corporation | Positioning and data integrating method and system thereof |
CN1644456A (zh) * | 2005-01-21 | 2005-07-27 | 清华大学 | 用于微小卫星的微型组合姿态测量*** |
CN101556155A (zh) * | 2009-05-20 | 2009-10-14 | 上海微小卫星工程中心 | 小卫星姿态确定***及方法 |
US20100312509A1 (en) * | 2009-06-05 | 2010-12-09 | Apple Inc. | Calibration techniques for an electronic compass in a portable device |
CN101644579A (zh) * | 2009-08-17 | 2010-02-10 | 中国人民解放军国防科学技术大学 | 一种空中微小抛撒物运动参数测量方法及*** |
WO2011152141A1 (ja) * | 2010-05-31 | 2011-12-08 | オリンパスメディカルシステムズ株式会社 | 内視鏡形状検出装置及び内視鏡の挿入部の形状検出方法 |
Non-Patent Citations (6)
Title |
---|
寇义民: "地磁导航关键技术研究", 《中国博士学位论文全文数据库》 * |
李玎,卜雄洙: "基于磁/惯性传感器旋转弹体定姿的Kalman滤波器设计", 《中国惯性技术学报》 * |
林春生等: "三轴磁强计正交误差分析与校正", 《探测与控制学报》 * |
贺飞: "高动态环境下MIMU/地磁组合姿测技术研究", 《中国优秀硕士学位论文全文数据库》 * |
赵捍东: "脉冲发送机提供控制力的火箭弹道修正理论及技术研究", 《中国博士学位论文全文数据库》 * |
闫辉等: "一种海洋环境地磁场三分量的测量方法", 《海军工程大学学报》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103424780B (zh) * | 2013-08-27 | 2016-05-18 | 中国航空无线电电子研究所 | 基于感应线圈的载机磁环境补偿方法 |
CN103424780A (zh) * | 2013-08-27 | 2013-12-04 | 中国航空无线电电子研究所 | 基于感应线圈的载机磁环境补偿方法 |
CN103926626A (zh) * | 2014-04-30 | 2014-07-16 | 中色地科矿产勘查股份有限公司 | 航磁数据高度矫正方法及*** |
CN104199115A (zh) * | 2014-09-18 | 2014-12-10 | 中国人民解放军国防科学技术大学 | 一种基于线性模型的地磁矢量测量误差综合补偿方法 |
CN104199115B (zh) * | 2014-09-18 | 2016-09-21 | 中国人民解放军国防科学技术大学 | 一种基于线性模型的地磁矢量测量误差综合补偿方法 |
CN104459828B (zh) * | 2014-12-08 | 2017-05-31 | 中国人民解放军国防科学技术大学 | 基于绕轴法的地磁矢量***非对准校正方法 |
CN104535062A (zh) * | 2015-01-20 | 2015-04-22 | 中国人民解放军国防科学技术大学 | 基于磁梯度张量和地磁矢量测量的运动式定位方法 |
CN104535062B (zh) * | 2015-01-20 | 2017-09-15 | 中国人民解放军国防科学技术大学 | 基于磁梯度张量和地磁矢量测量的运动式定位方法 |
CN104678340A (zh) * | 2015-02-27 | 2015-06-03 | 清华大学 | 一种磁强计测量误差纠正方法及*** |
CN104678340B (zh) * | 2015-02-27 | 2017-09-22 | 清华大学 | 一种磁强计测量误差纠正方法及*** |
CN105068131A (zh) * | 2015-08-03 | 2015-11-18 | 中国科学院电子学研究所 | 一种航磁数据调平方法 |
CN105572749A (zh) * | 2016-01-25 | 2016-05-11 | 中国地质大学(武汉) | 地面三分量磁力定向方法及地面三分量磁力定向勘探装置 |
US11118913B2 (en) | 2016-10-19 | 2021-09-14 | Huawei Technologies Co., Ltd. | Vehicle positioning correction method and mobile device |
CN107024724A (zh) * | 2017-06-22 | 2017-08-08 | 吉林大学 | 反卷积快速优化直升机航空电磁探测早期数据的方法 |
CN107493531A (zh) * | 2017-08-04 | 2017-12-19 | 歌尔科技有限公司 | 一种头部姿态检测方法、装置和耳机 |
CN110417489A (zh) * | 2018-04-26 | 2019-11-05 | 杭州海康机器人技术有限公司 | 一种无人机定向天线的校准方法及装置 |
CN110417489B (zh) * | 2018-04-26 | 2021-09-17 | 杭州海康机器人技术有限公司 | 一种无人机定向天线的校准方法及装置 |
CN110308490A (zh) * | 2019-06-13 | 2019-10-08 | 南京理工大学 | 基于磁传感器阵列的旋转运动磁异常目标探测装置及方法 |
CN111257940A (zh) * | 2020-02-14 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
CN111257940B (zh) * | 2020-02-14 | 2020-09-04 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
CN118050814A (zh) * | 2024-04-16 | 2024-05-17 | 山东省地质矿产勘查开发局第五地质大队(山东省第五地质矿产勘查院) | 低空光学定姿三分量磁测***及测定磁场矢量数据的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102621584B (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102621584A (zh) | 航磁矢量数据处理方法 | |
CN103630137B (zh) | 一种用于导航***的姿态及航向角的校正方法 | |
CN109556632A (zh) | 一种基于卡尔曼滤波的ins/gnss/偏振/地磁组合导航对准方法 | |
CN103630139B (zh) | 一种基于地磁梯度张量测量的水下载体全姿态确定方法 | |
CN1740746B (zh) | 微小型动态载体姿态测量装置及其测量方法 | |
Li et al. | Magnetic sensors for navigation applications: an overview | |
CN102520455B (zh) | 航空地磁矢量检测装置 | |
CN109541704B (zh) | 一种三轴磁通门航磁测量***及矫正补偿方法 | |
CN102393201B (zh) | 航空遥感用位置和姿态测量***(pos)动态杆臂补偿方法 | |
CN103363992B (zh) | 基于梯度下降的四旋翼无人机姿态航向参考***解算方法 | |
CN101393022B (zh) | 有磁环境的数字磁罗盘标定方法 | |
CN105509737B (zh) | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 | |
JP4412381B2 (ja) | 方位検出装置 | |
CN109556631A (zh) | 一种基于最小二乘的ins/gnss/偏振/地磁组合导航***对准方法 | |
CN103926627B (zh) | 水下载体地磁三分量测量方法 | |
CN104613983A (zh) | 一种应用于微小型无人机的整机磁强计校准方法 | |
CN102901977A (zh) | 一种飞行器的初始姿态角的确定方法 | |
CN110779496B (zh) | 三维地图构建***、方法、设备和存储介质 | |
CN110887472B (zh) | 一种偏振-地磁信息深度融合全自主姿态解算方法 | |
CN105737858A (zh) | 一种机载惯导***姿态参数校准方法与装置 | |
CN109507706B (zh) | 一种gps信号丢失的预测定位方法 | |
CN109612460B (zh) | 一种基于静止修正的垂线偏差测量方法 | |
US11408735B2 (en) | Positioning system and positioning method | |
CN112334736A (zh) | 用于校准物体磁力计的方法 | |
CN103278165A (zh) | 基于剩磁标定的磁测及星光备份的自主导航方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150722 Termination date: 20151214 |
|
EXPY | Termination of patent right or utility model |