CN104748722B - 利用卫星定位信息实时校准气压测高结果的高程定位方法 - Google Patents
利用卫星定位信息实时校准气压测高结果的高程定位方法 Download PDFInfo
- Publication number
- CN104748722B CN104748722B CN201510111831.0A CN201510111831A CN104748722B CN 104748722 B CN104748722 B CN 104748722B CN 201510111831 A CN201510111831 A CN 201510111831A CN 104748722 B CN104748722 B CN 104748722B
- Authority
- CN
- China
- Prior art keywords
- result
- satellite
- positioning
- elevation
- precision
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
- G01C5/06—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels by using barometric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Manufacturing & Machinery (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种利用卫星定位信息实时校准气压测高结果的高程定位方法,结合气压测高短时间段稳定、相对高程精度高的特点和在一些离散观测点卫星定位精度高的特点,利用卫星定位信息对气压测高结果进行实时校准,同时,通过实验和理论计算获得筛选卫星定位结果的条件公式,去除精度差的定位结果,保留精度高的结果对气压测高结果进行校准,因此可实现高精度、连续的高程定位,定位精度优于卫星定位高程精度和气压测高精度,可达到1m的精度,且该精度的实现可达到(1σ)66.7%的概率。
Description
技术领域
本发明涉及导航定位关键技术,提出了一种利用卫星定位信息实时校准气压测高结果的高程定位方法,实现室内外稳定、高精度的高程定位。
背景技术
随着卫星导航技术的发展,导航接收机的应用范围越来越广,在城市高楼间、树荫下、峡谷中等复杂环境下导航定位需求越来越多。然而导航接收机受环境影响较大,当用户在这些复杂环境下时,卫星导航信号将会被遮挡或干扰,这将导致终端定位结果不连续、定位精度降低,在高程定位振荡尤为明显。除此之外,当卫星几何结构不理想时,相比水平方向,高程方向更容易受到影响,精度下降严重。
为了解决上述问题,提出了一些高程测量方案,其中应用较广的是利用气压信息测高的方式,气压测高的原理是地球表面上随着高程的升降其气压也在不断发生变化,因此按照一定的规则,可以由当地气压信息计算得到高程结果。
然而,在气压测高定位中仍存在一些问题:气压测高有相对精度较高以及短期内较稳定的特点,而长时间工作时,受环境(风速、温度、湿度等)的影响,会发生随时间发生漂移的问题,出现连续几天测量结果均不一致的情况。
基于以上技术背景和应用背景,提出了一种利用卫星定位信息实时校准气压测高结果的高程定位方法,解决了气压测高随时间漂移的问题,可以实现室内外稳定、高精度的高程定位。
发明内容
有鉴于此,本发明的目的是提供一种利用卫星定位信息实时校准气压测高结果的高程定位方法,能够对气压测高结果进行实时校准,输出高精度的定位结果。
本发明的利用卫星定位信息实时校准气压测高结果的高程定位方法,包括如下步骤:
步骤1、利用卫星定位信息,实时获得移动目标的卫星定位结果,即待测点在大地坐标系的x,y,z轴上的坐标值,然后再分别获得所述三个坐标值的置信度Dx、Dy和Dz;
步骤2、利用待测点所在位置的实时气压信息获得移动目标的气压高程定位结果;
步骤3、针对步骤1实时获得的置信度,保留满足公式Dx+Dy+Dz≤1的置信度对应的卫星定位结果,将满足公式的卫星定位结果转换成经纬高坐标形式,获取其中的卫星高程定位结果;将实时获得的卫星高程定位结果与步骤2在同一时刻获得的气压高程定位结果求差值,然后将该各时刻的差值进行累加,累加一段设定的时间后,除以参与累加的差值个数m,得到累加结果的平均校准值bias;采用该平均校准值bias对当前时刻以后获得的气压高程定位进行校准。
较佳的,所述设定的时间选取为5分钟至30分钟。
较佳的,所述步骤2中根据实时气压信息获得移动目标的高程定位结果的方法为:
在所述移动目标上配备定位终端,采用该定位终端获取其接收范围内气象基准站获得的气象信息和位置信息,然后获得高程定位结果H:H=VBS.H+dh;其中,DisUE2BS[i]表示移动目标到第i个气象基准站的粗算距离,H为第i个气象基准站的高度;其中,N表示定位终端接收范围内基准站的个数;Tao和为常数,g为当地的重力加速度,T和P分别表示移动目标所在位置的温度和气压;BS[i].P为第i个气象基准站的气压值。
本发明具有如下有益效果:
本发明结合气压测高短时间段稳定、相对高程精度高的特点和在一些离散观测点卫星定位精度高的特点,利用卫星定位信息对气压测高结果进行实时校准,同时,通过实验和理论计算获得筛选卫星定位结果的条件公式,去除精度差的定位结果,保留精度高的结果对气压测高结果进行校准,因此可实现高精度、连续的高程定位,定位精度优于卫星定位高程精度和气压测高精度,可达到1m的精度,且该精度的实现可达到(1σ)66.7%的概率。
附图说明
图1为本发明的高程定位方法的流程图。
图2为本发明的三种测高结果的曲线对比图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明的利用卫星定位信息实时校准气压测高结果的高程定位方法,如图1所示,包括如下步骤:步骤1、利用卫星定位信息,实时获得移动目标的卫星定位结果,即待测点在大地坐标系的x,y,z轴上的坐标值,然后再分别获得所述三个坐标值的置信度Dx、Dy和Dz,具体为:
置信度计算涉及GPS定位原理,采用最小二乘方法。假设伪距观测量为L,其精度为σ。GPS定位方程组:
观测方程表达为:
L=BX
其中,X为定位结果,B为量测矩阵。
线性方程表达为:GδX=l
其中,δX=Xk-Xk-1,即上一状态与本次状态定位结果差值,l=L-BXk-1,G为雅可比矩阵,与各颗卫星相对于用户的几何位置有关。
误差估计结果:
其中,P为权阵,
σj为用户测距误差的方差,实际中,接收机根据卫星信号的强弱、卫星仰角高低和接收机跟踪环路的运行状态以及URA等指标,对不同时刻、不同卫星测量值的误差方差进行估算。误差的计算方法很多,以下将举例说明其中一种。
一般来说,随着卫星高度角的降低,导航信号的观测误差越大。因此,本文根据观测历元内各卫星的高度角计算相应观测值的方差,其中,非差观测值的方差计算方法如下:
b为常数项因子;f(e l)为观测历元内与卫星高度角e l相关的函数;常数项a和b可根据观测数据的先验信息综合计算得到或直接取经验值.
f(el)=1/sin(el)
根据误差传播定律,的误差方差为
Q=(GTPG)-1
置信度为误差方差:
DX=Q,其分量为Dx,Dy,Dz,Dt,即坐标X,Y,Z及时间t维度的置信度。
步骤2、根据气压信息计算高程定位结果:
采用基准站差分技术实现气压测高。在所述移动目标上配备定位终端,采用该定位终端获取其接收范围内气象基准站获得的气象信息和位置信息,然后定位终端结合气象基准站的气象信息和位置信息以及本地气象信息获取终端高程定位结果,具体为:
利用定位终端的粗位置计算定位终端到第i个基准站的距离为:
DisUE2BS[i]=|U.pos-BS[i].pos|
其中,U.pos为用户终端大地坐标值,BS[i].pos为第i个基准站大地坐标值。
计算定位终端的信息和气象信息:
其中,BS[i].P为第i个基准站的气压值,BS[i].T为第i个基准站的温度,BS[i].H为第i个基准站的高度,VBS.P为获得的定位终端所在处气压值,VBS.T为获得的定位终端所在处温度值,VBS.H为获得的定位终端所在处高度值。
计算高程:
其中N表示定位终端接收范围内基准站的个数;Tao=6.5,Rd=287.05,g为当地的重力加速度,T和P分别表示移动目标所在位置的温度和气压。
则H=VBS.H+dh;
步骤3、实时校准气压高程定位结果:
按照卫星导航定位结果的置信度筛选高精度卫星导航定位结果。通过以上公式推导,可以得到Dx,Dy,Dz三个方向置信度,该物理量表征三个方向的误差方差。为了获取高精度的定位结果,选择筛选条件为Dx+Dy+Dz≤d2,其中,根据方差计算特性,Dx+Dy+Dz表示为定位结果单位矢量方向误差方差,d为精度门限,当Dx+Dy+Dz≤1,即可筛选出精度为1的定位结果。以上条件为理论值,在工程实践中,通过采集大量的数据分析该筛选条件的虚警率和漏检率,可以证实该置信度筛选条件可用性。漏检率指,对于一些实测数据,定位结果满足1米的精度要求,但是其置信度不满足上述条件,此类数据没有被正确筛选出来,此类数据占所有数据的比例即为漏检率。虚警率指,对于一些实测数据,定位结果不满足1米的精度要求,但是其置信度却满足上述条件,此类数据没有被正确筛选出来,此类数据占所有数据的比例即为虚警率。上述筛选条件达到在漏检率为13%左右,虚警率在10%以内。
按照卫星导航定位结果的置信度筛选高精度卫星导航定位结果,筛选条件为Dx+Dy+Dz≤1,将高精度卫星导航定位结果经过坐标转换得到经纬高,取卫星导航高程定位结果和气压高程定位结果差值,并进行积分累加,累加时间为n分钟,根据筛选条件,累计m个可用于校准的数据,最后统计平均得到校准值bias,随后时间内将气压定位结果加上校准值bias的结果作为最终的高程结果输出。
Ho(t)=Hb(t)+bias
其中,Hg(k)为k时刻GNSS定位高程结果,Hb(k)为k时刻气压测高高程定位结果。在定位过程中,每n秒更新一次校准值,对气压定位结果进行实时校准,输出高精度定位结果。鉴于气压定位结果的短期稳定连续、长期受环境影响易漂移的特点,校准后的定位结果可达到长时间连续、稳定和高精度。典型的场景是室内外定位和卫星受遮挡情景。室内外定位过程中,在室外环境下,气压测高定位结果易漂移,校准后高程稳定、精度高,进入室内环境,校准信息依然有效,仍可以实现连续、稳定、高精度定位结果。如图2所示,为本发明的三种测高结果的曲线对比图,可发现本发明的定位方法具有更高的精度。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.利用卫星定位信息实时校准气压测高结果的高程定位方法,其特征在于,包括如下步骤:
步骤1、利用卫星定位信息,实时获得移动目标的卫星定位结果,即待测点在大地坐标系的x,y,z轴上的坐标值,然后再分别获得所述三个坐标值的置信度Dx、Dy和Dz;
步骤2、利用待测点所在位置的实时气压信息获得移动目标的气压高程定位结果;
步骤3、针对步骤1实时获得的置信度,保留满足公式Dx+Dy+Dz≤1的置信度对应的卫星定位结果,将满足公式的卫星定位结果转换成经纬高坐标形式,获取其中的卫星高程定位结果;将实时获得的卫星高程定位结果与步骤2在同一时刻获得的气压高程定位结果求差值,然后将该各时刻的差值进行累加,累加一段设定的时间后,除以参与累加的差值个数m,得到累加结果的平均校准值bias;采用该平均校准值bias对当前时刻以后获得的气压高程定位结果进行校准。
2.如权利要求1所述的利用卫星定位信息实时校准气压测高结果的高程定位方法,其特征在于,所述设定的时间选取为5分钟至30分钟。
3.如权利要求1所述的利用卫星定位信息实时校准气压测高结果的高程定位方法,其特征在于,所述步骤2中根据实时气压信息获得移动目标的高程定位结果的方法为:
在所述移动目标上配备定位终端,采用该定位终端获取其接收范围内气象基准站获得的气象信息和位置信息,然后获得高程定位结果H:H=VBS.H+dh;
其中,DisUE2BS[i]表示移动目标到第i个气象基准站的粗算距离,BS[i].H为第i个气象基准站的高度;其中,N表示定位终端接收范围内基准站的个数;Tao和Rd为常数,Tao=6.5,Rd=287.05,g为当地的重力加速度,T和P分别表示移动目标所在位置的温度和气压;BS[i].P为第i个气象基准站的气压值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510111831.0A CN104748722B (zh) | 2015-03-13 | 2015-03-13 | 利用卫星定位信息实时校准气压测高结果的高程定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510111831.0A CN104748722B (zh) | 2015-03-13 | 2015-03-13 | 利用卫星定位信息实时校准气压测高结果的高程定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104748722A CN104748722A (zh) | 2015-07-01 |
CN104748722B true CN104748722B (zh) | 2017-09-08 |
Family
ID=53588768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510111831.0A Active CN104748722B (zh) | 2015-03-13 | 2015-03-13 | 利用卫星定位信息实时校准气压测高结果的高程定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104748722B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105509770B (zh) * | 2016-01-07 | 2018-07-17 | 中国科学院嘉兴微电子与系统工程中心 | 一种gnss和mems组合导航系统中的气压计在线校正方法 |
US11761765B2 (en) | 2016-09-09 | 2023-09-19 | Nextnav, Llc | Calibrating a pressure sensor |
US10514258B2 (en) | 2016-09-09 | 2019-12-24 | Nextnav, Llc | Systems and methods for calibrating unstable sensors |
CN106840203B (zh) * | 2017-01-10 | 2020-01-17 | 南京航空航天大学 | 惯导/气压高度计/gps组合导航系统中气压高度计校正方法 |
CN106814382B (zh) * | 2017-01-11 | 2019-05-10 | 武汉大学 | 联合高度角和用户测距误差的gnss卫星观测值定权方法 |
US11248908B2 (en) | 2017-02-24 | 2022-02-15 | Here Global B.V. | Precise altitude estimation for indoor positioning |
WO2019228605A1 (en) * | 2018-05-28 | 2019-12-05 | Here Global B.V. | Multilevel altitude maps |
CN111027019B (zh) * | 2018-10-10 | 2023-09-08 | 千寻位置网络有限公司 | 统计分析高精度定位结果置信度的方法及装置 |
CN111811471A (zh) * | 2020-06-26 | 2020-10-23 | 黄珊 | 气压计校准方法、标签高度确定方法、标签、uwb定位系统 |
CN114018217B (zh) * | 2021-09-29 | 2023-09-05 | 河南省联睿智能科技研究院有限公司 | 基于二次校准的鲁棒差分绝对高度估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102620716A (zh) * | 2012-03-12 | 2012-08-01 | 无锡泽太微电子有限公司 | 高度测量方法、三维定位方法、高度测量设备及移动终端 |
US8612172B2 (en) * | 2006-09-19 | 2013-12-17 | Nokia Corporation | Relative positioning |
CN103712599A (zh) * | 2013-12-27 | 2014-04-09 | 南京信息工程大学 | 一种相对高度测量装置及方法 |
WO2014087598A1 (ja) * | 2012-12-04 | 2014-06-12 | 旭化成株式会社 | 気圧式高度計および屋内対応型気圧式高度計 |
-
2015
- 2015-03-13 CN CN201510111831.0A patent/CN104748722B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8612172B2 (en) * | 2006-09-19 | 2013-12-17 | Nokia Corporation | Relative positioning |
CN102620716A (zh) * | 2012-03-12 | 2012-08-01 | 无锡泽太微电子有限公司 | 高度测量方法、三维定位方法、高度测量设备及移动终端 |
WO2014087598A1 (ja) * | 2012-12-04 | 2014-06-12 | 旭化成株式会社 | 気圧式高度計および屋内対応型気圧式高度計 |
CN103712599A (zh) * | 2013-12-27 | 2014-04-09 | 南京信息工程大学 | 一种相对高度测量装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104748722A (zh) | 2015-07-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104748722B (zh) | 利用卫星定位信息实时校准气压测高结果的高程定位方法 | |
CN110645979B (zh) | 基于gnss/ins/uwb组合的室内外无缝定位方法 | |
CN105607093B (zh) | 一种组合导航系统及获取导航坐标的方法 | |
KR102034527B1 (ko) | 와이-파이 위치와 센서정보 결합을 통한 단말 위치 필터링 시스템 및 그 방법 | |
CN101354253B (zh) | 一种基于匹配度的地磁辅助导航方法 | |
CN110487301A (zh) | 一种雷达辅助机载捷联惯性导航系统初始对准方法 | |
CN109029454A (zh) | 一种基于卡尔曼滤波的横坐标系捷联惯导系统阻尼算法 | |
Yang et al. | Global autonomous positioning in GNSS-challenged environments: A bioinspired strategy by polarization pattern | |
CN104049269B (zh) | 一种基于激光测距和mems/gps组合导航系统的目标导航测绘方法 | |
CN106886039B (zh) | 基于城市三维地图的地面数字电视广播与导航卫星混合粒子滤波定位方法 | |
CN105527642B (zh) | 一种单星定位装置与方法 | |
CN108693547A (zh) | 一种用于水下深潜器的导航系统及精准三点定位方法 | |
CN101784908A (zh) | 卫星定位系统中的独立高度测量 | |
CN110617795A (zh) | 一种利用智能终端的传感器实现室外高程测量的方法 | |
US20240159529A1 (en) | Systems and methods for extending the spatial coverage of a reference pressure network | |
CN108151765A (zh) | 一种在线实时估计补偿磁强计误差的定位测姿方法 | |
CN105115478A (zh) | 机动经纬仪基座平台晃动测量与修正方法 | |
CN114088080A (zh) | 一种基于多传感器数据融合的定位装置及方法 | |
RU2697859C1 (ru) | Способ определения местоположения наземного подвижного объекта | |
Rózsa | Modelling tropospheric delays using the global surface meteorological parameter model GPT2 | |
US20220377499A1 (en) | Characterizing height above terrain confidence | |
WO2020044278A1 (en) | Methods for geolocation using electronic distance measurement equipment | |
CN105068104A (zh) | 一种基于惯性/双星间断伪距约束的定位方法 | |
CN115143878A (zh) | 基于方差膨胀模型的gnss/加速度计自适应融合滑坡监测方法 | |
JP2019203795A (ja) | ダイナミック座標管理装置、ダイナミック座標管理方法、プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |