CN103115624A - 一种基于地磁匹配的地磁日变修正方法 - Google Patents
一种基于地磁匹配的地磁日变修正方法 Download PDFInfo
- Publication number
- CN103115624A CN103115624A CN201310028198XA CN201310028198A CN103115624A CN 103115624 A CN103115624 A CN 103115624A CN 201310028198X A CN201310028198X A CN 201310028198XA CN 201310028198 A CN201310028198 A CN 201310028198A CN 103115624 A CN103115624 A CN 103115624A
- Authority
- CN
- China
- Prior art keywords
- matching
- geomagnetic
- terrestrial magnetism
- geomagnetism
- component
- 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
Abstract
本发明公开了一种基于地磁匹配的地磁日变修正方法。针对地磁匹配导航定位中的地磁日变效应,本发明提出基于FMI方法与多维特征量匹配对地磁日变效应进行修正,包括如下步骤:1)选取地磁匹配特征量;2)基于FMI方法从地磁匹配特征量的实测数据中提取各地磁匹配特征量的日变量;3)通过地磁匹配特征量的实测数据减去与其对应的日变量,建立地磁实时图;4)基于多维特征量匹配对地磁实时图和地磁基准图进行匹配;5)输出匹配结果。本发明能较为有效地消除地磁匹配导航定位中的地磁日变效应,能显著提高地磁匹配导航定位精度。
Description
技术领域
本发明涉及一种基于地磁匹配的地磁日变修正方法,属于高精度地磁匹配导航的关键技术领域。
背景技术
地磁导航是当前国际导航领域的一大研究热点,具有高度自主性、隐蔽性、无辐射、能耗低等诸多优点,能弥补惯导系统位置误差长期累积的不足,也能克服GPS等卫星有源导航系统易被侦测和受干扰等缺陷,已成为隐蔽性运动载体无源导航的重要手段,广泛应用于潜艇、导弹和航天器辅助导航等领域,展现出巨大的军事价值和应用前景。然而,地磁匹配导航定位精度不高严重制约了其在相关领域的工程应用。
地磁日变是影响地磁匹配导航精度的重要因素之一。在地磁匹配导航研究中,许多学者指出了地磁日变的重要性,但没有开展深入的研究。目前的研究表明地磁日变对地磁导航的影响是不容忽视的,但存在以下问题:
①多数研究对地磁日变影响只进行了简单的探讨;
②部分研究直接将日变以噪声的形式加入到仿真中;
③部分研究完全忽略地磁日变效应。
传统日变场的分析方法有日均值法和子夜均值法。其中日均值法模糊了日变白天变化大、夜间变化小的特点。子夜均值法以子夜均值为基线值,但基线值的选择具有很大的人为影响因素,容易给匹配定位带来新的误差源。FMI(Finish Meteorological Institute,芬兰气象局)方法是一种基于三天的地磁场X、Y分量分钟值观测数据为数据源计算中间一天K指数的计算方法,此方法由该气象局提出。与前述两者相比,FMI方法是对观测数据的直接拟合,不需选取基线值,减少了人为误差,所提取的日变更为客观、真实,具有较大的适用性。
发明内容
本发明所要解决的技术问题是提供一种能用于地磁匹配导航定位中的地磁日变修正方法,基于FMI方法提取地磁日变成分,在多维地磁特征量匹配中消除地磁日变效应,从而提高地磁匹配导航定位精度。
本发明为解决上述技术问题采用以下技术方案:
一种基于地磁匹配的地磁日变修正方法,包括以下步骤:
步骤1,选取地磁匹配特征量,所述地磁匹配特征量为三个地磁分量,包括北向分量X、东向分量Y、垂直分量Z;
步骤2,基于FMI方法从地磁匹配特征量的实测数据中提取各地磁匹配特征量的日变量;
步骤3,通过地磁匹配特征量的实测数据减去与其对应的日变量,建立地磁实时图;
步骤4,基于多维特征量匹配对地磁实时图和地磁基准图进行匹配;
步骤5,输出匹配结果。
作为本发明的基于地磁匹配的地磁日变修正方法的进一步优化方案,所述步骤2涉及的实测数据是由磁传感器采集的;所述基于FMI方法从实测数据中提取日变量包括以下步骤:
步骤2.1,根据地磁匹配特征量的扰动幅度a分别确定与各地磁匹配特征量相对应的K指数的初始值KX0、KY0、KZ0,取K0=max(KX0,KY0,KZ0),所述扰动幅度a为每3小时时段内地磁匹配特征量的实测数据的最大值和最小值间的差值;
步骤2.2,地磁匹配特征量的时均值由该小时及其前后各(m+n)分钟的数据平均值来计算,然后利用5次谐波函数对时均值进行拟合,生成SR曲线;所述n=K33;m与地方时有关:地方时3点到6点和18点到21点,m取60;地方时7点到17点,m取0;其余时间,m取90;
步骤2.3,分别计算地磁匹配特征量的实测数据与SR曲线的残差序列,由残差序列的3小时扰动幅度a分别确定与各地磁匹配特征量相对应的K指数的中间值KXm、KYm、KZm,取Km=max(KXm,KYm,KZm);
步骤2.4,重复步骤2.2和步骤2.3的计算一次,再次得到K指数的中间值
Km=max(KXm,KYm,KZm),此Km即为K指数的最终值;
步骤2.5,由K指数的最终值,代入步骤2.2计算,生成最终SR曲线;
步骤2.6,计算最终SR曲线与地球主磁场和磁异常间的差值,得到地磁匹配特征量的日变量。
作为本发明的基于地磁匹配的地磁日变修正方法的进一步优化方案,所述步骤4涉及的地磁基准图是指基于地磁模型建立的基准数据库,或基于测量资料绘制的地磁图。
作为本发明的基于地磁匹配的地磁日变修正方法的进一步优化方案,步骤4涉及的多维特征量匹配是指通过3个地磁匹配特征量按照先后顺序分别进行匹配,即以北向分量X、东向分量Y、垂直分量Z这三个独立特征量进行匹配,具体包含以下步骤:
步骤4.1,以北向分量X作为第一匹配特征量,根据子夜均值附近北向分量X的日变大小设定阈值λX,以北向分量X的实时数据BX在地磁基准数据库Rx,y,z={Bx,By,Bz}中搜索满足给定地磁匹配算法条件的匹配点,设该匹配结果为点集Ry,z={B′y,B′z};其中,Bx,By,Bz分别表示三个地磁分量上的地磁基准数据,B′y,B′z分别表示在基准数据库中搜索出的东向分量、垂直分量上的地磁基准数据;
步骤4.2,以东向分量Y作为第二匹配特征量,根据子夜均值附近东向分量Y的日变大小设定阈值λY,以东向分量Y的实时数据BY在集合Ry,z中搜索满足给定地磁匹配算法条件的匹配点,设该匹配结果为点集Rz={B″z},其中,B″z表示在集合Ry,z中进一步搜索出的垂直分量上的地磁基准数据;
步骤4.3,以垂直分量Z作为第三匹配特征量,以垂直分量Z的实时数据BZ在集合Rz中搜索满足给定地磁匹配算法条件的匹配点,该点即为最佳匹配位置。
作为本发明的基于地磁匹配的地磁日变修正方法的进一步优化方案,所述地磁匹配算法为常用的相关分析匹配算法,包括平均平方差算法MSD、平均绝对差算法MAD、平方差算法SD、绝对差算法AD、积相关算法PROD。
作为本发明的基于地磁匹配的地磁日变修正方法的进一步优化方案,所述地磁匹配算法为平均绝对差算法MAD,具体如下:
其中,P为实时数据序列,Q为基准数据序列,实时数据序列P和基准数据序列Q的长度相同,N表示实时数据序列P或基准数据序列Q的数据序列长度,pi表示实时数据序列P中第i个值,qi表示基准数据序列Q中第i个值,i、N均为正整数。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
本发明将FMI方法应用到地磁匹配导航中,基于该方法提取地磁静日期间地球变化磁场的最主要成分——地磁日变,并建立地磁实时图,通过多维地磁特征量进行匹配,消除地磁匹配中的地磁日变效应,从而提高地磁匹配导航定位精度。
同时,本发明能有效地消除地磁匹配中的日变效应,显著地减小了匹配误差。
附图说明
图1是日变修正方法流程图。
图2是基于FMI方法提取日变量的流程图。
图3是多维特征量匹配流程图。
图4是2005年7月25日广州台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图5是2005年7月25日兰州台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图6是2005年7月25日北京十三陵台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图7是2009年7月2日广州台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图8是2009年7月2日兰州台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图9是2009年7月2日北京十三陵台站X、Y、Z分量的观测值曲线、拟合曲线和日变曲线。
图10是2005年7月25日广州、兰州和北京十三陵三台站日变修正前后的匹配误差图。
图11是2009年7月2日广州、兰州和北京十三陵三台站日变修正前后的匹配误差图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本实施例是基于地磁台站观测数据、FMI方法和多维特征量匹配进行的。以分布在不同纬度上的三个地磁台站(广州GZH、兰州LZH和北京十三陵BMT)在太阳活动高年(2005年)和太阳活动低年(2009年)的地磁观测数据分钟值模拟地磁实测数据。三个台站的地理位置见下表:
地磁台站 | 纬度(°) | 经度(°) | 海拔(m) |
广州 | 23.093N | 113.343E | 11 |
兰州 | 36.087N | 103.845E | 1560 |
北京十三陵 | 40.3N | 116.2E | 183 |
对实测数据进行预处理:所选用的地磁台站数据已对环境干扰进行了处理,但需要进行地方时校正,并剔除坏值。
参照图1,基于地磁匹配的地磁日变修正方法的分析步骤如下:
步骤1,选取地磁匹配特征量:地磁场有七个分量,分别为北向分量X、东向分量Y、垂直分量Z、磁偏角D、水平强度H、磁倾角I和总强度F。其中,X、Y、Z三个分量是独立的。本实施例选取X、Y、Z三分量作为地磁匹配特征量。
步骤2,基于FMI方法从地磁匹配特征量的实测数据中提取各地磁匹配特征量的日变量:参照图2所示提取流程,具体步骤为:
步骤2.1,分别从各地磁台站提取X、Y、Z三分量的分钟值观测数据;根据观测数据的3小时扰动幅度a(3小时时段内的最大值和最小值间的差值)分别确定各匹配特征量对应的K指数(描述单个地磁台站3小时时段内地磁扰动强度的指数)的初始值KX0、KY0、KZ0,取K0=max(KX0,KY0,KZ0);所述K指数从0到9共分10级,由3小时扰动幅度a确定:a<3nT,K取0;3≤a<6nT,K取1;6≤a<12nT,K取2;12≤a<24nT,K取3;24≤a<40nT,K取4;40≤a<70nT,K取5;70≤a<120nT,K取6;120≤a<200nT,K取7;200≤a<300nT,K取8;a≥300nT,K取9。
步骤2.2,地磁X、Y、Z三分量的时均值由该小时及其前后各(m+n)分钟的数据平均值来计算,然后利用5次谐波函数(对时间序列x=[x1,x2,...,XN],x包含N个数据,其5次谐波函数为 其中, t=1,2,....,N,j=1,2,...,5)对时均值进行拟合,生成SR曲线(Solar Regular variations,规则日变化);所述n=K33,m与地方时有关:地方时3点到6点和18点到21点,m取60;地方时7点到17点,m取0;其余时间,m取90。
步骤2.3,分别计算地磁X、Y、Z三分量的观测值与SR曲线的残差序列,由残差序列的3小时扰动幅度a分别确定与各地磁匹配特征量相对应的K指数的中间值KXm、KYm、KZm,取Km=max(KXm,KYm,KZm);
步骤2.4,重复步骤2.2和步骤2.3的计算一次,再次得到K指数的中间值
Km=max(KXm,KYm,KZm),此时的Km即为K指数的最终值;
步骤2.5,由K指数的最终值,代入步骤2.2计算,生成最终SR曲线;
步骤2.6,计算最终SR曲线与地球主磁场和磁异常间的差值,得到地磁三分量X、Y、Z的日变量。
参照图4至图9,分别为2005年7月25日和2009年7月2日广州、兰州和北京十三陵三个地磁台站X、Y、Z三分量的日变量处理结果。
步骤3,通过地磁X、Y、Z三分量的观测数据分别减去对应的日变量建立地磁实时图。
步骤4,基于多维特征量匹配对地磁实时图和地磁基准图进行匹配。
本实施例基于多维特征量匹配对地磁实时图和基准图进行匹配,具体如下:
采用IGRF11(the11th generation International Geomagnetic Reference Field,第11代国际参考地磁场)+磁异常,分别以广州、兰州和北京十三陵三个台站所在位置为中心建立2005年7月25日和2009年7月2日的地磁基准图,位置参数采用经纬度,取网格大小为1000×1000,网格精度为0.001°×0.001°。
参照图3,本实施例所涉及的多维特征量匹配具体是指:
(1)选择MAD(平均绝对差)作为匹配算法:
其中,P为实时数据序列,Q为基准数据序列,实时数据序列P和基准数据序列Q的长度相同,N表示实时数据序列P或基准数据序列Q的数据序列长度,pi表示实时数据序列P中第i个值,qi表示基准数据序列Q中第i个值,i、N均为正整数。
本发明适用的相关分析匹配算法还有很多,常规的比如:MSD(Mean Square Deviation,平均平方差)、MAD(Mean Absolute Deviation,平均绝对差)、SD(Square Deviation,平方差)、AD(Absolute Deviation,绝对差)、PROD(Product Correlation,积相关)等,都属于常用的匹配算法,在这里就不一一赘述。
(2)匹配过程:
①以北向分量X作为第一匹配特征量,设定阈值λX(地磁日变具有在子夜附近几乎不变的特点,因此,阈值的设定参考子夜均值附近北向分量X的日变大小),以北向分量X的实时数据BX在基准地磁数据库Rx,y,z={Bx,By,Bz}中搜索使BX满足MAD(BX,Bx)<λX的匹配点,设该匹配结果为点集Ry,z={B′y,B′z};其中,Bx,By,Bz分别表示X、Y、Z三个地磁分量上的地磁基准数据,B′y,B′z分别表示在基准数据库中搜索出的东向分量、垂直分量上的地磁基准数据;
②以东向分量Y作为第二匹配特征量,设定阈值λY(地磁日变具有在子夜附近几乎不变的特点,因此,阈值的设定参考子夜均值附近东向分量Y的日变大小),以东向分量Y的实时数据BY在集合Ry,z中搜索使BY满足MAD(BY,B′y)<λY的匹配点,设该匹配结果为点集Rz={B″z};其中,B″z表示在集合Ry,z中进一步搜索出的垂直分量上的地磁基准数据;
③以垂直分量Z作为第三匹配特征量,以垂直分量Z的实时数据BZ在集合Rz中搜索使BZ满足min(MAD(BZ,B″z))的点,该点即为最佳匹配位置。图10和图11分别是2005年7月25日和2009年7月2日广州、兰州和北京十三陵三个地磁台站日变修正前后的匹配误差图。参阅图10和图11可知,本发明在地方时6点至18点之间修正效果尤为明显,匹配定位精度能提高一个数量级。
(1)综合来看,地磁匹配中的日变效应显著,随太阳活动性、地方时、纬度而分布。日变修正前后的匹配误差曲线具有明显差异:修正前,匹配误差波动大,变化剧烈;修正后,匹配误差曲线波动小,变化较平缓。尤其在地方时6点到18点之间,日变修正后匹配定位误差明显减小,匹配定位精度得到显著提高。在图10中,广州台站白天14时,未进行日变修正前经度误差达到0.4°,而日变修正后经度误差仅为0.05°左右,匹配定位精度提高了约一个数量级。
(2)从图10和图11对比来看,日变修正效果在太阳活动低年要优于太阳活动高年。这是因为,在太阳活动高年,地磁活动剧烈,除地磁日变外,还存在一些较大的扰动,如地磁脉动、地磁钩扰等。图10中可以看出,在太阳活动高年,在低纬度上的日变效应尤为显著。经日变修正后,广州台站纬度方向上误差仍可达到0.1°,经度方向上达0.2°。而图11中,在太阳活动低年,经日变修正后,三个台站纬度方向上误差控制在0.005°左右,经度方向上则控制在0.05°左右。
(3)三个台站地理位置不同,但修正效果都比较明显,说明该方法具有较大的适用性。
Claims (6)
1.一种基于地磁匹配的地磁日变修正方法,其特征在于,包括以下步骤:
步骤1,选取地磁匹配特征量,所述地磁匹配特征量为三个地磁分量,包括北向分量X、东向分量Y、垂直分量Z;
步骤2,基于FMI方法从地磁匹配特征量的实测数据中提取各地磁匹配特征量的日变量;
步骤3,通过地磁匹配特征量的实测数据减去与其对应的日变量,建立地磁实时图;
步骤4,基于多维特征量匹配对地磁实时图和地磁基准图进行匹配;
步骤5,输出匹配结果。
2.根据权利要求1所述的基于地磁匹配的地磁日变修正方法,其特征在于,所述步骤2涉及的实测数据是由磁传感器采集的;所述基于FMI方法从实测数据中提取日变量包括以下步骤:
步骤2.1,根据地磁匹配特征量的扰动幅度a分别确定与各地磁匹配特征量相对应的K指数的初始值KX0、KY0、KZ0,取K0=max(KX0,KY0,KZ0),所述扰动幅度a为每3小时时段内地磁匹配特征量的实测数据的最大值和最小值间的差值;
步骤2.2,地磁匹配特征量的时均值由该小时及其前后各(m+n)分钟的数据平均值来计算,然后利用5次谐波函数对时均值进行拟合,生成SR曲线;所述n=K33;m与地方时有关:地方时3点到6点和18点到21点,m取60;地方时7点到17点,m取0;其余时间,m取90;
步骤2.3,分别计算地磁匹配特征量的实测数据与SR曲线的残差序列,由残差序列的3小时扰动幅度a分别确定与各地磁匹配特征量相对应的K指数的中间值KXm、KYm、KZm,取Km=max(KXm,KYm,KZm);
步骤2.4,重复步骤2.2和步骤2.3的计算一次,再次得到K指数的中间值
Km=max(KXm,KYm,KZm),此Km即为K指数的最终值;
步骤2.5,由K指数的最终值,代入步骤2.2计算,生成最终SR曲线;
步骤2.6,计算最终SR曲线与地球主磁场和磁异常间的差值,得到地磁匹配特征量的日变量。
3.根据权利要求1所述的基于地磁匹配的地磁日变修正方法,其特征在于,所述步骤4涉及的地磁基准图是指基于地磁模型建立的基准数据库,或基于测量资料绘制的地磁图。
4.根据权利要求1所述的基于地磁匹配的地磁日变修正方法,其特征在于,步骤4涉及的多维特征量匹配是指通过3个地磁匹配特征量按照先后顺序分别进行匹配,即以北向分量X、东向分量Y、垂直分量Z这三个独立特征量进行匹配,具体包含以下步骤:
步骤4.1,以北向分量X作为第一匹配特征量,根据子夜均值附近北向分量X的日变大小设定阈值λX,以北向分量X的实时数据BX在地磁基准数据库Rx,y,z={Bx,By,Bz}中搜索满足给定地磁匹配算法条件的匹配点,设该匹配结果为点集Ry,z={B′y,B′z};其中,Bx,By,Bz分别表示三个地磁分量上的地磁基准数据,B′y,B′z分别表示在基准数据库中搜索出的东向分量、垂直分量上的地磁基准数据;
步骤4.2,以东向分量Y作为第二匹配特征量,根据子夜均值附近东向分量Y的日变大小设定阈值λY,以东向分量Y的实时数据BY在集合Ry,z中搜索满足给定地磁匹配算法条件的匹配点,设该匹配结果为点集Rz={B″z},其中,B″z表示在集合Ry,z中进一步搜索出的垂直分量上的地磁基准数据;
步骤4.3,以垂直分量Z作为第三匹配特征量,以垂直分量Z的实时数据BZ在集合Rz中搜索满足给定地磁匹配算法条件的匹配点,该点即为最佳匹配位置。
5.根据权利要求4所述的基于地磁匹配的地磁日变修正方法,其特征在于,所述地磁匹配算法为常用的相关分析匹配算法,包括平均平方差算法MSD、平均绝对差算法MAD、平方差算法SD、绝对差算法AD、积相关算法PROD。
6.根据权利要求4或5所述的基于地磁匹配的地磁日变修正方法,其特征在于,所述地磁匹配算法为平均绝对差算法MAD,具体如下:
其中,P为实时数据序列,Q为基准数据序列,实时数据序列P和基准数据序列Q的长度相同,N表示实时数据序列P或基准数据序列Q的数据序列长度,pi表示实时数据序列P中第i个值,qi表示基准数据序列Q中第i个值,i、N均为正整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310028198.XA CN103115624B (zh) | 2013-01-24 | 2013-01-24 | 一种基于地磁匹配的地磁日变修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310028198.XA CN103115624B (zh) | 2013-01-24 | 2013-01-24 | 一种基于地磁匹配的地磁日变修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103115624A true CN103115624A (zh) | 2013-05-22 |
CN103115624B CN103115624B (zh) | 2014-12-10 |
Family
ID=48414044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310028198.XA Expired - Fee Related CN103115624B (zh) | 2013-01-24 | 2013-01-24 | 一种基于地磁匹配的地磁日变修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103115624B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104807462A (zh) * | 2015-04-30 | 2015-07-29 | 中测高科(北京)测绘工程技术有限责任公司 | 室内地磁导航基准图生成方法及系统 |
CN105043387A (zh) * | 2015-06-26 | 2015-11-11 | 武汉科技大学 | 基于惯导辅助地磁的个人室内定位系统 |
CN105758403A (zh) * | 2016-01-06 | 2016-07-13 | 东南大学 | 一种基于Vague集模糊推理的地磁图适配性的综合评价方法 |
CN106197405A (zh) * | 2016-08-01 | 2016-12-07 | 哈尔滨工业大学 | 地磁日变影响下的惯性地磁匹配定位方法 |
CN106767789A (zh) * | 2017-01-12 | 2017-05-31 | 南京航空航天大学 | 一种基于自适应卡尔曼滤波的行人航向最优融合方法 |
CN106908058A (zh) * | 2017-03-03 | 2017-06-30 | 哈尔滨工程大学 | 一种确定地磁定位阵列孔径的方法 |
CN106918350A (zh) * | 2017-02-27 | 2017-07-04 | 中国科学院国家空间科学中心 | 一种应用于地磁导航中的地磁场模型误差补偿方法 |
WO2017174843A1 (es) * | 2016-04-05 | 2017-10-12 | Universidad De Alcala | Dispositivo y procedimiento de obtención de la perturbación geomagnética local a latitudes medias |
CN107462241A (zh) * | 2017-06-27 | 2017-12-12 | 中测高科(北京)测绘工程技术有限责任公司 | 结合图像的地磁基础图的测量方法、装置和系统 |
CN107543541A (zh) * | 2017-08-25 | 2018-01-05 | 北京麦钉艾特科技有限公司 | 一种适合室内自由运动载体的地磁定位方法 |
CN109073380A (zh) * | 2016-05-04 | 2018-12-21 | 比利时皇家气象研究中心 | 自校准且自主的地磁观测台 |
CN111895995A (zh) * | 2020-06-03 | 2020-11-06 | 东南大学 | 基于pso的飞行器编队多维地磁匹配导航方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1530023A2 (en) * | 2003-11-10 | 2005-05-11 | The Foundation: The Research Institute of Electric and Magnetic Materials | Magnetic north detecting device and magnetic north detecting method |
US20060293844A1 (en) * | 2005-06-20 | 2006-12-28 | Denso Corporation | Vehicle controller |
CN2872384Y (zh) * | 2006-02-28 | 2007-02-21 | 国家海洋局第一海洋研究所 | 自持式浮标地磁日变站 |
CN101354253A (zh) * | 2008-09-19 | 2009-01-28 | 北京航空航天大学 | 一种基于匹配度的地磁辅助导航方法 |
CN101354252A (zh) * | 2008-09-19 | 2009-01-28 | 北京航空航天大学 | 一种基于多尺度估计的地磁辅助导航算法 |
CN101520328A (zh) * | 2009-04-01 | 2009-09-02 | 西北工业大学 | 地磁场线图匹配自主导航方法 |
CN102322858A (zh) * | 2011-08-22 | 2012-01-18 | 南京航空航天大学 | 用于地磁/捷联惯导组合导航系统的地磁匹配导航方法 |
CN102445201A (zh) * | 2011-09-28 | 2012-05-09 | 东北林业大学 | 用于水下载体的地磁异常特征点匹配导航方法 |
-
2013
- 2013-01-24 CN CN201310028198.XA patent/CN103115624B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1530023A2 (en) * | 2003-11-10 | 2005-05-11 | The Foundation: The Research Institute of Electric and Magnetic Materials | Magnetic north detecting device and magnetic north detecting method |
US20060293844A1 (en) * | 2005-06-20 | 2006-12-28 | Denso Corporation | Vehicle controller |
CN2872384Y (zh) * | 2006-02-28 | 2007-02-21 | 国家海洋局第一海洋研究所 | 自持式浮标地磁日变站 |
CN101354253A (zh) * | 2008-09-19 | 2009-01-28 | 北京航空航天大学 | 一种基于匹配度的地磁辅助导航方法 |
CN101354252A (zh) * | 2008-09-19 | 2009-01-28 | 北京航空航天大学 | 一种基于多尺度估计的地磁辅助导航算法 |
CN101520328A (zh) * | 2009-04-01 | 2009-09-02 | 西北工业大学 | 地磁场线图匹配自主导航方法 |
CN102322858A (zh) * | 2011-08-22 | 2012-01-18 | 南京航空航天大学 | 用于地磁/捷联惯导组合导航系统的地磁匹配导航方法 |
CN102445201A (zh) * | 2011-09-28 | 2012-05-09 | 东北林业大学 | 用于水下载体的地磁异常特征点匹配导航方法 |
Non-Patent Citations (2)
Title |
---|
王建军: "中国地区地磁场日变化时空分布特征初步研究", 《中国优秀硕士学位论文全文数据库基础科学辑》, no. 10, 15 October 2008 (2008-10-15) * |
邓翠婷等: "地磁匹配导航算法综述", 《科学技术与工程》, vol. 12, no. 24, 31 August 2012 (2012-08-31) * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104807462B (zh) * | 2015-04-30 | 2018-09-07 | 中测高科(北京)测绘工程技术有限责任公司 | 室内地磁导航基准图生成方法及系统 |
CN104807462A (zh) * | 2015-04-30 | 2015-07-29 | 中测高科(北京)测绘工程技术有限责任公司 | 室内地磁导航基准图生成方法及系统 |
CN105043387A (zh) * | 2015-06-26 | 2015-11-11 | 武汉科技大学 | 基于惯导辅助地磁的个人室内定位系统 |
CN105758403A (zh) * | 2016-01-06 | 2016-07-13 | 东南大学 | 一种基于Vague集模糊推理的地磁图适配性的综合评价方法 |
CN105758403B (zh) * | 2016-01-06 | 2018-12-14 | 东南大学 | 一种基于Vague集模糊推理的地磁图适配性的综合评价方法 |
WO2017174843A1 (es) * | 2016-04-05 | 2017-10-12 | Universidad De Alcala | Dispositivo y procedimiento de obtención de la perturbación geomagnética local a latitudes medias |
CN109073380A (zh) * | 2016-05-04 | 2018-12-21 | 比利时皇家气象研究中心 | 自校准且自主的地磁观测台 |
CN109073380B (zh) * | 2016-05-04 | 2020-11-10 | 比利时皇家气象研究中心 | 自校准且自主的地磁观测台 |
CN106197405A (zh) * | 2016-08-01 | 2016-12-07 | 哈尔滨工业大学 | 地磁日变影响下的惯性地磁匹配定位方法 |
CN106197405B (zh) * | 2016-08-01 | 2019-02-26 | 哈尔滨工业大学 | 地磁日变影响下的惯性地磁匹配定位方法 |
CN106767789B (zh) * | 2017-01-12 | 2019-12-24 | 南京航空航天大学 | 一种基于自适应卡尔曼滤波的行人航向最优融合方法 |
CN106767789A (zh) * | 2017-01-12 | 2017-05-31 | 南京航空航天大学 | 一种基于自适应卡尔曼滤波的行人航向最优融合方法 |
CN106918350B (zh) * | 2017-02-27 | 2019-08-20 | 中国科学院国家空间科学中心 | 一种应用于地磁导航中的地磁场模型误差补偿方法 |
CN106918350A (zh) * | 2017-02-27 | 2017-07-04 | 中国科学院国家空间科学中心 | 一种应用于地磁导航中的地磁场模型误差补偿方法 |
CN106908058A (zh) * | 2017-03-03 | 2017-06-30 | 哈尔滨工程大学 | 一种确定地磁定位阵列孔径的方法 |
CN106908058B (zh) * | 2017-03-03 | 2019-12-10 | 哈尔滨工程大学 | 一种确定地磁定位阵列孔径的方法 |
CN107462241A (zh) * | 2017-06-27 | 2017-12-12 | 中测高科(北京)测绘工程技术有限责任公司 | 结合图像的地磁基础图的测量方法、装置和系统 |
CN107462241B (zh) * | 2017-06-27 | 2021-08-17 | 中测高科(北京)测绘工程技术有限责任公司 | 结合图像的地磁基础图的测量方法、装置和系统 |
CN107543541B (zh) * | 2017-08-25 | 2019-11-08 | 北京麦钉艾特科技有限公司 | 一种适合室内自由运动载体的地磁定位方法 |
CN107543541A (zh) * | 2017-08-25 | 2018-01-05 | 北京麦钉艾特科技有限公司 | 一种适合室内自由运动载体的地磁定位方法 |
CN111895995A (zh) * | 2020-06-03 | 2020-11-06 | 东南大学 | 基于pso的飞行器编队多维地磁匹配导航方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103115624B (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103115624B (zh) | 一种基于地磁匹配的地磁日变修正方法 | |
Melton et al. | Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP) | |
CN102708307B (zh) | 一种应用于城市的植被指数构造方法 | |
Molteni et al. | The new ECMWF seasonal forecast system (System 4) | |
CN106909722A (zh) | 一种近地面气温的大面积精准反演方法 | |
CN104021309B (zh) | 基于方位偏转法的测风数据风向订正方法 | |
CN106197460B (zh) | 一种应用gps出行数据进行出行目的地预测的方法 | |
CN104239706A (zh) | 一种地面观测气温时空数据集的制备方法 | |
CN102592181A (zh) | 一种作物种植面积统计数据空间分布的优化方法 | |
CN102445201A (zh) | 用于水下载体的地磁异常特征点匹配导航方法 | |
CN104035138A (zh) | 一种全球及局部海洋扰动重力的精确快速计算方法 | |
Iacono et al. | Seasonal variability of the Tyrrhenian Sea surface geostrophic circulation as assessed by altimeter data | |
CN110378576A (zh) | 城市化植被效应有效距离的定量化检测方法 | |
CN104936287A (zh) | 基于矩阵补全的传感网室内指纹定位方法 | |
Jing et al. | Spatial variability analysis of soil nutrients based on GIS and geostatistics: a case study of Yisa Township, Yunnan, China | |
CN104406698A (zh) | 一种城市热岛空间分布的评价方法 | |
Yang et al. | Characteristics of vegetation cover change in Xilin Gol League based on MODIS EVI data | |
De Paula et al. | Characteristics of the ionospheric F-region plasma irregularities over Brazilian longitudinal sector | |
CN108957502A (zh) | 基于深度学习lstm的gnss多系统多路径误差实时削弱方法 | |
CN104657935A (zh) | 一种近地面气温反演方法 | |
CN105184048A (zh) | 一种基于多源数据的逐年土地利用空间数据的制备方法 | |
CN101793977A (zh) | 一种水文地质参数的估计方法 | |
CN109916396A (zh) | 一种基于多维地磁信息的室内定位方法 | |
Newitt et al. | Location of the north magnetic pole in April 2007 | |
Araujo-Pradere et al. | Consistent features of TEC changes during ionospheric storms |
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: 20141210 Termination date: 20210124 |
|
CF01 | Termination of patent right due to non-payment of annual fee |