CN110275192B - 一种基于智能手机的高精度单点定位方法与装置 - Google Patents

一种基于智能手机的高精度单点定位方法与装置 Download PDF

Info

Publication number
CN110275192B
CN110275192B CN201910428335.6A CN201910428335A CN110275192B CN 110275192 B CN110275192 B CN 110275192B CN 201910428335 A CN201910428335 A CN 201910428335A CN 110275192 B CN110275192 B CN 110275192B
Authority
CN
China
Prior art keywords
satellite
observation
smart phone
positioning
observation value
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.)
Expired - Fee Related
Application number
CN201910428335.6A
Other languages
English (en)
Other versions
CN110275192A (zh
Inventor
高成发
陈波
刘永胜
孙璞玉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201910428335.6A priority Critical patent/CN110275192B/zh
Publication of CN110275192A publication Critical patent/CN110275192A/zh
Priority to US17/051,428 priority patent/US11709281B2/en
Priority to PCT/CN2020/072411 priority patent/WO2020233158A1/zh
Application granted granted Critical
Publication of CN110275192B publication Critical patent/CN110275192B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/425—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between signals derived from different satellite radio beacon positioning systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/51—Relative positioning
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • G01S19/072—Ionosphere corrections
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/14—Receivers specially adapted for specific applications
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/243—Demodulation of navigation message
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/33—Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种基于智能手机的高精度单点定位方法与装置，本发明方法改进了常规的PPP非组合定位模型，无需GNSS参考站，仅使用智能手机接收的原始GNSS观测值进行高精度定位，属于卫星定位技术领域。本发明的定位方法包括以下步骤：获取智能手机原始GNSS伪距和载波相位等观测值；在进行数据预处理削弱部分误差影响之后，根据改进的估计双钟差的精密单点定位方法，对原始观测值形成非组合模型；根据卫星高度角确定各卫星观测值权重；通过改进的卡尔曼滤波方法进行滤波定位，得到高精度单点定位结果。使用本发明提出的定位方法，能够对智能手机实现平面优于0.9m，高程优于1.7m的定位精度，且收敛时间小于30s。

Description

一种基于智能手机的高精度单点定位方法与装置
技术领域
本发明涉及全球导航卫星系统(Global Navigation Satellite System，GNSS)卫星定位方法，具体涉及一种基于智能手机原始GNSS观测值的高精度单点定位方法与装置。
背景技术
智能手机是如今人们生活中不可或缺的工具，智能手机中的GNSS(GlobalNavigation Satellite System，全球导航卫星系统)模块已经极大地改善了现代人类生活。在GNSS导航定位技术的发展中，导航或定位精度一直是制约其进一步应用于人类生产生活并发挥巨大作用的关键问题，对于智能手机的GNSS导航定位也是如此。由于一直以来智能手机中的GNSS导航定位模块都封装于手机操作系统之中，研究人员只能获取定位最终结果来进行应用层面的开发，因此相对于传统测地型接收机，智能手机原始观测值的分析与定位算法研究都非常缺乏。
从2016年开始，谷歌公司在Android N操作系统上提供了访问GNSS原始观测数据的接口，学者们开始着手于对智能手机接收的GNSS数据进行质量评估，并对定位性能进行分析。目前国内外学者关于利用原始GNSS观测值实现智能手机高精度定位的研究，主要分为三个方向：1)移动终端GNSS原始观测值数据质量分析；2)综合利用各种类型的GNSS数据，通过滤波的方式提高定位精度；3)利用伪距或载波差分的方式，提升移动终端定位精度。这些研究都是将现有的导航定位增强手段应用到智能手机上来，并且大部分的研究对象都是Android平板电脑，对于普通大众使用的智能手机，目前的定位性能较差。更重要的是，目前还没有根据智能手机原始观测值性质改进现有定位方案的先例。
发明内容
发明目的：基于以上背景技术，针对测试发现的智能手机伪距观测值和载波观测值变化率不一致的性质，改进PPP(Precise Point Positioning，精密单点定位)算法，实现高精度的智能手机实时单点定位。
技术方案：为实现上述发明目的，本发明采用的技术方案为：
一种基于智能手机的高精度单点定位方法，包括以下步骤：
(1)获取智能手机导航定位模块的GNSS伪距和载波相位原始观测值；
(2)在进行数据预处理削弱部分误差影响之后，根据改进的估计双钟差的精密单点定位方法，对原始观测值形成非组合观测值模型；
(3)根据卫星高度角确定各卫星观测值权重；
(4)使用改进的静态卡尔曼滤波进行计算，得到高精度的单点定位结果。
在优选的实施方案中，步骤(1)中，利用智能手机操作系统提供的基于位置服务的API直接获取部分GNSS原始数据，包括时间数据、GNSS系统类型和载波相位观测值，然后根据GNSS原始数据中的时间数据通过信号传播时间差值解算出伪距。
在优选的实施方案中，步骤(2)中，所削弱的误差影响包括对接收到观测值使用精密星历、精密钟差文件消除轨道误差和卫星钟差，使用电离层格网文件削弱电离层延迟，以及将多路径效应视作观测噪声，所形成的非组合观测值模型为：
Figure GDA0002779780660000021
Figure GDA0002779780660000022
Figure GDA0002779780660000023
Figure GDA0002779780660000024
Figure GDA0002779780660000025
Figure GDA0002779780660000026
式中，上标g、e、c分别表示GPS系统、Galileo系统和BDS系统，下标i、j、k表示第i、j、k颗卫星，左侧P、Φ分别为伪距和载波观测值，ρ为手机终端到卫星的距离，c为光速，
Figure GDA0002779780660000027
分别为伪距和载波观测值的移动端钟差，dtrop为对流层延迟，
Figure GDA0002779780660000028
为载波整周模糊度，∈P、∈Φ分别为伪距和载波相位观测值的残差；
Figure GDA0002779780660000029
和
Figure GDA00027797806600000210
分别表示Galileo和BDS系统与GPS系统之间的时间偏差。
在优选的实施方案中，步骤(3)中，根据卫星高度角确定各卫星观测值权重的定权方案具体为：
当卫星高度角小于10°时，该卫星观测值权重为0；
当卫星高度角大于等于10°时，该卫星观测值权重为sinE，E为卫星高度角。
在优选的实施方案中，步骤(4)中，改进的静态卡尔曼滤波具体为：
对于单个历元，根据步骤(2)所述的非组合观测值模型，一颗卫星的观测方程数为2个，设某历元观测到GPS卫星n1颗，Galileo卫星n2颗，BDS卫星n3颗，满足n1+n2+n3≥8；此时卡尔曼滤波中待估参数向量为：
Figure GDA0002779780660000031
式中，x、y和z为智能手机位置参数；则观测方程总数为2×(n1+n2+n3)个，待估参数为n1+n2+n3+8个，多余观测数为n1+n2+n3-8个；
卡尔曼滤波的观测方程系数阵H为：
Figure GDA0002779780660000032
矩阵H中，行数为2×(n1+n2+n3)，列数为n1+n2+n3-8，下标1至n表示卫星序号，n＝2×(n1+n2+n3)，上标表示GNSS系统类型，g、e、c分别表示GPS系统、Galileo系统和BDS系统；奇数行对应伪距观测值，偶数行对应载波相位观测值；前3列α、β、γ为卫星-移动端的方向余弦，第4、5列为伪距观测值和载波观测值的钟差系数，奇数行第4列为1，偶数行第5列为1；第6、7列为系统间偏差系数，当观测值来自Galileo卫星时，第6列为1，当观测值来自BDS卫星时，第7列为1；MF为对流层湿延迟投影系数；第9列之后为载波相位模糊度系数，奇数行为0，偶数行第i+8列为1，i为卫星序号；滤波过程中各卫星观测值的权重由步骤(3)确定。
基于相同的发明构思，本发明提供的一种基于智能手机的高精度单点定位装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的一种基于智能手机的高精度单点定位方法。
有益效果：本发明提出一种基于智能手机的高精度单点定位方法，在现有的定位算法上，根据手机端原始GNSS观测值独特性质，优化了非组合PPP观测值模型，能够在不需要参考站的情况下，实现普通智能手机的亚米级定位精度，且收敛速度较快。基于手机端的高精度定位技术，可以向大众用户会提供更好的城市定位、车载导航和旅游交通等用户体验，在正在到来智能互联时代，更高精度的位置服务意味着难以预估的应用前景，物联网、自动驾驶和智慧城市等现代人类生活发展方向中，无一不包含着低成本的位置服务，因此本发明具有重大意义。
附图说明
图1是基于智能手机的高精度单点定位方法的流程图；
图2是同步观测的测地型接收机和智能手机的观测值对比图；
图3是本发明实施例的定位结果图。
具体实施方式
下面结合附图对本发明的技术方案作进一步详细说明。
图1为本发明实施例公开的基于智能手机的高精度单点定位方法流程图。本实施例中，首先通过Android系统(7.0及以上版本)提供的基于位置服务的API，获取GNSS原始伪距和载波观测值等原始观测值。在进行数据预处理削弱部分误差影响之后，根据改进的估计双钟差的精密单点定位算法，对原始观测值形成非组合PPP模型。然后采用卫星高度角定位模型，确定各卫星观测值权重。最后，使用改进的静态卡尔曼滤波方式，得到高精度、快速收敛的实时单点定位结果。以下叙述详细过程。
第一步，通过Android系统提供的基于位置服务的API，获取GNSS伪距和载波相位等原始观测值。首先利用Android系统开发中的android.location包中的GnssMeasurement类与GnssClock类直接获取部分GNSS原始数据，包括时间数据、GNSS系统类型、载波相位(表1第1行)、伪距率、载噪比，GnssMeasurement类和GnssClock类包含的原始数据内容如表1、表2所示。
表1 GnssMeasurement类主要GNSS原始数据表
Figure GDA0002779780660000041
表2 GnssClock类主要GNSS原始数据
Figure GDA0002779780660000051
然后根据GNSS原始数据中的时间数据解算出伪距，其实质是通过信号传播时间差值计算伪距：
P＝(tRx-tTx)·c (1)
其中：P表示伪距，tRx表示接收机接收到信号的时刻，tTx表示接收到卫星发射信号的时刻，c为光速。tTx可通过getReceivedSvTimeNanos()方法直接获取，数值为本GPS周内的秒数。
Figure GDA0002779780660000052
Figure GDA0002779780660000053
其中
Figure GDA0002779780660000054
表示从1980.6.1到当前时间的总秒数，GNSSweek表示所有整GPS周秒数和，TimeNanos为接收机本地时钟值getTimeNanos()，TimeOffsetNanos为本地时钟偏移值getTimeOffsetNanos()，FullBiasNanos为接收机本地时钟总偏差getFullBiasNanos()，BiasNanos为接收机本地钟差getBiasNanos()。
第二步，在进行数据预处理削弱部分误差影响之后，根据改进的估计双钟差的PPP算法，对原始观测值形成非组合观测值模型。
前期实验将智能手机和测地型接收机置于相同地点进行同步观测，收集原始观测值，并将智能手机原始观测值文件转换为Rinex格式文件。图2为相同时间相同卫星测地型接收机的原始观测值和智能手机(华为P10手机)的观测值对比。图2中，接收机载波与伪距观测值重合在一起，实际上二者之差为载波相位的整周模糊度×载波波长。测地型接收机观测值在一段时间之后会有一个大幅度的跳动，这是接收机钟跳引起的，不影响原始观测值的定位解算。而可以明显地观察出手机伪距与载波相位观测值相互之间的差值不固定，这个性质与测地型接收机是不相同的，影响手机终端载波相位数据的使用。
对同一设备不同卫星原始观测值进行分析，发现同一时刻不同卫星的伪距变化率和载波观测值变化率之间的差异是相同的，故在定位时认为手机终端伪距观测值和载波观测值具有不同的钟差。
根据以上性质改进定位方案，在参数求解时估计两个钟差参数。以GPS系统为例，单频观测方程为：
Figure GDA0002779780660000061
Figure GDA0002779780660000062
式中，上标g表示GPS系统，下标i表示第i颗卫星，Pi g、
Figure GDA0002779780660000063
为伪距和载波观测值，
Figure GDA0002779780660000064
为手机端到卫星的距离，c为光速，
Figure GDA0002779780660000065
分别为伪距和载波观测值的钟差，dTg为卫星钟差，
Figure GDA0002779780660000066
为对流层延迟，
Figure GDA0002779780660000067
为卫星轨道误差，
Figure GDA0002779780660000068
为电离层延迟，
Figure GDA0002779780660000069
分别为伪距和载波多路径效应延迟，
Figure GDA00027797806600000610
为载波整周模糊度，
Figure GDA00027797806600000611
分别伪距和载波相位观测值的残差。
实时下载IGS(International GNSS Service,国际GNSS服务组织)提供的预测精密星历、钟差文件和预测电离层格网文件，消除卫星轨道误差和钟差，削弱电离层延迟误差，并通过模型改正消除相对论效应和地球自转等误差影响，将多路径效应视作观测噪声后，简化观测值方程，得到非组合观测值模型：
Figure GDA00027797806600000612
Figure GDA00027797806600000613
Figure GDA00027797806600000614
Figure GDA00027797806600000615
Figure GDA00027797806600000616
Figure GDA00027797806600000617
式中，左侧为各个GNSS系统的伪距和载波观测值，上标e、c分别表示Galileo系统和BDS系统，下标j、k表示第j、k颗卫星，
Figure GDA00027797806600000618
和
Figure GDA00027797806600000619
分别表示Galileo和BDS系统与GPS系统之间的时间偏差。
第三步，确定各卫星观测值权重。本方法根据卫星高度角确定各卫星观测值权重，定权方案具体为：
Figure GDA0002779780660000071
式中，W为权重，E为手机端到卫星连线的高度角。
第四步，使用卡尔曼滤波进行参数求解，得到高精度定位结果。
对观测值形成第二步所示的非组合模型，此时，单个历元一颗卫星的观测方程数为2个，假设某历元观测到GPS卫星n1颗，Galileo卫星n2颗，BDS卫星n3颗，此时卡尔曼滤波中待估参数向量为：
Figure GDA0002779780660000072
式中，x、y和z为智能手机位置参数。
则观测方程总数为2×(n1+n2+n3)个，待估参数为n1+n2+n3+8个，多余观测数为n1+n2+n3-8个。
本发明采用的卡尔曼滤波，具体步骤为：
(1)输入系统状态初始值初值X0，及其方差P0。X即为式(13)中的待估参数向量，下标0表示初始时刻，即设置待估参数的初值，此时前3个参数为智能手机近似坐标(通过伪距单点定位获得)，之后的参数分别为伪距钟差、载波钟差、系统间偏差、对流层湿延迟和n1+n2+n3颗卫星的模糊度，全设为0。
(2)计算进一步预测值Xk,k-1及预测值的方差-协方差阵Pk,k-1(下标k为历元数)：
Xk，k-1＝φk，k-1Xk-1 (14)
Pk，k-1＝φk，k-1Pk-1φT k，k-1+Γk，k-1Qk-1ΓT k，k-1 (15)
式中，Фk,k-1为状态转移矩阵(因为卫星数量会改变，需要将上个历元的n阶Xk-1变为这个历元m阶Xk)。Γk,k-1为系统噪声驱动阵，Qk-1为系统误差(模型的误差)为正定矩阵，每个历元都要进行输入，可以直接将Γk,k-1Qk-1ΓT k,k-1看为一个整体，定为m阶对称矩阵。在首次定位是跟P0一样，在之后历元的定位中，三个位置参数和n个模糊度参数真值应该不变，所以过程噪声为0(系统噪声对应的方差)，手机钟差、系统间偏差及对流层延迟参数的随机模型采用随机游走过程模型模拟。
(3)计算增益矩阵Kk：
Figure GDA0002779780660000073
式中，Hk为观测方程的系数阵。相对于常规的静态卡尔曼滤波，本方法中的Hk需要进行相应的修改：
Figure GDA0002779780660000081
Hk行数为2×(n1+n2+n3)，列数为n1+n2+n3+8，式中n＝2×(n1+n2+n3)；下标表示卫星序号，上标表示GNSS系统类型；奇数行对应伪距观测值，偶数行对应载波相位观测值；前3列，α、β、γ为卫星-移动端的方向余弦，计算公式见式(18)；第4、5列为伪距观测值和载波观测值的钟差系数，奇数行第4列为1，偶数行第5列为1；第6、7列为系统间偏差系数，当观测值来自Galileo卫星时，第6列为1，当观测值来自BDS卫星时，第7列为1；第8列，MF为对流层湿延迟投影系数(对流层干延迟通过Hopfield模型改正，投影函数为Niell模型)。第9～n1+n2+n3+8列为载波相位模糊度系数，第2×i行，i+8列为1，i为卫星序号。
Figure GDA0002779780660000082
式中，xs、ys和zs为卫星坐标(通过星历文件获得)，x0、y0和z0为智能手机近似坐标(通过伪距单点定位获得)，ρ0为手机端到卫星的近似距离(通过卫星坐标和手机近似坐标计算)。
式(16)中Rk为观测噪声方差矩阵：
Figure GDA0002779780660000083
式中，σP和σΦ分别为伪距观测值和载波观测值的天顶方向标准差，W为各卫星权重，即式(12)计算的各卫星观测值权重。卡尔曼滤波过程中，观测值的方差设置非常重要，设置不当容易造成滤波发散，严重影响定位结果。本发明中，观测值天顶方向标准差根据先前实验测得的智能手机原始GNSS观测值数据质量来确定，实施例所使用的小米8手机，伪距观测值和载波相位观测值天顶方向标准差分别设置为3.0m和0.2m。
(4)计算当前历元滤波参数估值向量Xk和对应的方差-协方差阵Pk：
Xk＝Xk,k-1+Kk(Lk-HkXk,k-1) (20)
Pk＝(I-KkHk)Pk,k-1(I-KkHk)T+KkRkKk T (21)
式中，Lk为观测值向量，即式(6)～(11)中各GNSS系统的伪距和载波观测值，I为单位阵。
计算完成后，Xk向量中前3项即为当前历元计算的智能手机精确坐标。若当前历元不是最后一个历元，则返回步骤(2)进行循环滤波计算，最终得到多历元滤波解。
实验例：
实验时间为2018.10.19，使用小米8智能手机，测试地点为东南大学校园内的GE01控制点，进行了5次观测，每次6分钟左右，采样间隔1s。此测站精确坐标已经通过网络RTK的方式事先得到。
使用本发明方法对小米8手机进行定位测试，结果如图3所示。在观测一段时间之后，5个时段数据的N、E方向误差都能收敛到1m以下。
统计5个时段的定位中误差，如表3所示。平均的平面中误差为0.81m，而平均高程中误差为1.65m，这个结果优于目前使用普通智能手机能够获得的最高定位精度。
表3实施例各时段定位中误差
Figure GDA0002779780660000091
定义收敛时间为，从定位开始时刻到N和E方向定位误差均小于1m，且后续历元误差不再超过1m的时刻为止经过的时间。统计5个时段的收敛时间如表4所示，各时段均能在30s之内实现收敛，说明本方法能够应用于实时定位，提供低延迟的高精度智能手机定位结果。
表4实施例各时段收敛时间
Figure GDA0002779780660000101
基于相同的发明构思，本发明实施例公开的一种基于智能手机的高精度单点定位装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的一种基于智能手机的高精度单点定位方法。

Claims (5)

1.一种基于智能手机的高精度单点定位方法，其特征在于，包括以下步骤：
(1)获取智能手机导航定位模块的GNSS伪距和载波相位原始观测值；
(2)在进行数据预处理削弱部分误差影响之后，根据改进的估计双钟差的精密单点定位方法，对原始观测值形成非组合观测值模型；所削弱的误差影响包括对接收到观测值使用精密星历、精密钟差文件消除轨道误差和卫星钟差，使用电离层格网文件削弱电离层延迟，以及将多路径效应视作观测噪声，所形成的非组合观测值模型为：
Figure FDA0002788799220000011
Figure FDA0002788799220000012
Figure FDA0002788799220000013
Figure FDA0002788799220000014
Figure FDA0002788799220000015
Figure FDA0002788799220000016
式中，上标g、e、c分别表示GPS系统、Galileo系统和BDS系统，下标i、j、k表示第i、j、k颗卫星，左侧P、Φ分别为伪距和载波观测值，ρ为手机终端到卫星的距离，c为光速，
Figure FDA0002788799220000017
分别为伪距和载波观测值的移动端钟差，dtrop为对流层延迟，
Figure FDA0002788799220000018
为载波整周模糊度，∈P、∈Φ分别为伪距和载波相位观测值的残差；
Figure FDA0002788799220000019
和
Figure FDA00027887992200000110
分别表示Galileo和BDS系统与GPS系统之间的时间偏差；
(3)根据卫星高度角确定各卫星观测值权重；
(4)使用改进的静态卡尔曼滤波进行计算，得到高精度的单点定位结果。
2.根据权利要求1所述的基于智能手机的高精度单点定位方法，其特征在于，步骤(1)中，利用智能手机操作系统提供的基于位置服务的API直接获取部分GNSS原始数据，包括时间数据、GNSS系统类型和载波相位观测值，然后根据GNSS原始数据中的时间数据通过信号传播时间差值解算出伪距。
3.根据权利要求1所述的基于智能手机的高精度单点定位方法，其特征在于，步骤(3)中，根据卫星高度角确定各卫星观测值权重的定权方案具体为：
当卫星高度角小于10°时，该卫星观测值权重为0；
当卫星高度角大于等于10°时，该卫星观测值权重为sin E，E为卫星高度角。
4.根据权利要求1所述的基于智能手机的高精度单点定位方法，其特征在于，步骤(4)中，改进的静态卡尔曼滤波具体为：
对于单个历元，根据步骤(2)所述的非组合观测值模型，一颗卫星的观测方程数为2个，设某历元观测到GPS卫星n1颗，Galileo卫星n2颗，BDS卫星n3颗，满足n1+n2+n3≥8；此时卡尔曼滤波中待估参数向量为：
Figure FDA0002788799220000021
式中，x、y和z为智能手机位置参数；则观测方程总数为2×(n1+n2+n3)个，待估参数为n1+n2+n3+8个，多余观测数为n1+n2+n3-8个；
卡尔曼滤波的观测方程系数阵H为：
Figure FDA0002788799220000022
矩阵H中，行数为2×(n1+n2+n3)，列数为n1+n2+n3-8，下标1至n表示卫星序号，n＝2×(n1+n2+n3)，上标表示GNSS系统类型，g、e、c分别表示GPS系统、Galileo系统和BDS系统；奇数行对应伪距观测值，偶数行对应载波相位观测值；前3列α、β、γ为卫星-移动端的方向余弦，第4、5列为伪距观测值和载波观测值的钟差系数，奇数行第4列为1，偶数行第5列为1；第6、7列为系统间偏差系数，当观测值来自Galileo卫星时，第6列为1，当观测值来自BDS卫星时，第7列为1；MF为对流层湿延迟投影系数；第9列之后为载波相位模糊度系数，奇数行为0，偶数行第i+8列为1，i为卫星序号；滤波过程中各卫星观测值的权重由步骤(3)确定。
5.一种基于智能手机的高精度单点定位装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述计算机程序被加载至处理器时实现根据权利要求1-4任一项所述的一种基于智能手机的高精度单点定位方法。
CN201910428335.6A 2019-05-22 2019-05-22 一种基于智能手机的高精度单点定位方法与装置 Expired - Fee Related CN110275192B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201910428335.6A CN110275192B (zh) 2019-05-22 2019-05-22 一种基于智能手机的高精度单点定位方法与装置
US17/051,428 US11709281B2 (en) 2019-05-22 2020-01-16 High-precision point positioning method and device based on smartphone
PCT/CN2020/072411 WO2020233158A1 (zh) 2019-05-22 2020-01-16 一种基于智能手机的高精度单点定位方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910428335.6A CN110275192B (zh) 2019-05-22 2019-05-22 一种基于智能手机的高精度单点定位方法与装置

Publications (2)

Publication Number Publication Date
CN110275192A CN110275192A (zh) 2019-09-24
CN110275192B true CN110275192B (zh) 2021-01-26

Family

ID=67959071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910428335.6A Expired - Fee Related CN110275192B (zh) 2019-05-22 2019-05-22 一种基于智能手机的高精度单点定位方法与装置

Country Status (3)

Country Link
US (1) US11709281B2 (zh)
CN (1) CN110275192B (zh)
WO (1) WO2020233158A1 (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110275192B (zh) * 2019-05-22 2021-01-26 东南大学 一种基于智能手机的高精度单点定位方法与装置
CN111045052B (zh) * 2019-10-14 2022-08-09 广东星舆科技有限公司 一种智能终端伪距差分定位及质量控制方法
CN110687556B (zh) * 2019-11-04 2021-06-22 中国电子科技集团公司第五十四研究所 一种适用于laas的多径误差模型化方法
CN111060945B (zh) * 2020-01-13 2021-10-19 东南大学 一种gnss/5g紧组合融合定位方法与装置
US11650324B2 (en) * 2020-05-01 2023-05-16 Deere & Company Navigation apparatus and method in which measurement quantization errors are modeled as states
CN111796313B (zh) * 2020-06-28 2023-07-21 中国人民解放军63921部队 卫星定位方法及装置、电子设备、存储介质
CN113156477B (zh) * 2021-04-15 2022-08-26 中国科学院精密测量科学与技术创新研究院 一种基于Android智能手机的高精度RTK定位方法
CN113204042B (zh) * 2021-05-21 2022-11-18 北京交通大学 一种基于精密单点定位的多星座联合列车定位方法
CN113267793B (zh) * 2021-05-26 2022-05-06 中国电子科技集团公司第五十四研究所 一种基于外部增强信息的gbas对流层参数生成方法
CN114488227B (zh) * 2022-01-26 2023-10-20 西南交通大学 一种基于空间相关性的多路径误差改正方法
CN114640950A (zh) * 2022-03-06 2022-06-17 南京理工大学 基于Android源生GPS定位API的移动设备定位方法及系统
CN116184441B (zh) * 2022-12-20 2023-10-20 北京航空航天大学 一种精密大气改正增强的智能终端ppp瞬时亚米级定位方法
CN116678421B (zh) * 2023-06-12 2024-01-23 深圳沧穹科技有限公司 基于多模块ble发射装置的多源融合定位方法及系统
CN116819585B (zh) * 2023-08-31 2023-12-29 长沙金维信息技术有限公司 基于非线性优化的gnss单点定位方法及导航方法
CN117647830B (zh) * 2024-01-29 2024-05-07 中国测绘科学研究院 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104236579A (zh) * 2014-10-09 2014-12-24 武汉大学 一种基于Android内核层实现高精度卫星导航定位方法
CN104714244A (zh) * 2015-03-31 2015-06-17 东南大学 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法
CN106324629A (zh) * 2016-11-07 2017-01-11 中国电子科技集团公司第二十研究所 一种bds_gps_glonass融合精密单点定位方法
CN107807373A (zh) * 2017-10-17 2018-03-16 东南大学 基于移动智能终端的gnss高精度定位方法
CN109343090A (zh) * 2018-07-10 2019-02-15 东南大学 一种嵌入式gps/ bds/glonass实时精密单点定位设备
CN109709591A (zh) * 2018-12-07 2019-05-03 中国科学院光电研究院 一种面向智能终端的gnss高精度定位方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9651667B2 (en) * 2007-06-22 2017-05-16 Trimble Inc. Combined cycle slip indicators for regionally augmented GNSS
CN102576080B (zh) * 2009-09-19 2013-12-18 天宝导航有限公司 用以估计相位分级的时钟的gnss信号处理
CN104597465B (zh) 2015-01-23 2017-03-15 河海大学 一种提高gps与glonass组合精密单点定位收敛速度的方法
CN107356947B (zh) * 2017-05-31 2019-06-18 中国科学院测量与地球物理研究所 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN108363079B (zh) * 2018-01-30 2022-06-10 上海交通大学 一种面向便携式智能设备的gnss伪距双差定位方法及系统
CN110275192B (zh) * 2019-05-22 2021-01-26 东南大学 一种基于智能手机的高精度单点定位方法与装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104236579A (zh) * 2014-10-09 2014-12-24 武汉大学 一种基于Android内核层实现高精度卫星导航定位方法
CN104714244A (zh) * 2015-03-31 2015-06-17 东南大学 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法
CN106324629A (zh) * 2016-11-07 2017-01-11 中国电子科技集团公司第二十研究所 一种bds_gps_glonass融合精密单点定位方法
CN107807373A (zh) * 2017-10-17 2018-03-16 东南大学 基于移动智能终端的gnss高精度定位方法
CN109343090A (zh) * 2018-07-10 2019-02-15 东南大学 一种嵌入式gps/ bds/glonass实时精密单点定位设备
CN109709591A (zh) * 2018-12-07 2019-05-03 中国科学院光电研究院 一种面向智能终端的gnss高精度定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
广域实时亚米级单点定位算法;王延染等;《山东大学学报（工学版）》;20180228;第48卷(第1期);第31-35页 *

Also Published As

Publication number Publication date
CN110275192A (zh) 2019-09-24
US20220155465A1 (en) 2022-05-19
WO2020233158A1 (zh) 2020-11-26
US11709281B2 (en) 2023-07-25

Similar Documents

Publication Publication Date Title
CN110275192B (zh) 一种基于智能手机的高精度单点定位方法与装置
Cai et al. Precise point positioning with quad-constellations: GPS, BeiDou, GLONASS and Galileo
Li et al. Review of PPP–RTK: Achievements, challenges, and opportunities
CN105044747B (zh) 一种基于多星共视和滤波的时间同步装置及其方法
CN106842268A (zh) 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法
CN107765275A (zh) 广域差分定位方法、装置、终端及计算机可读存储介质
CN104898145B (zh) 一种基于半周模糊度的模糊度固定方法和系统
CN107807373A (zh) 基于移动智能终端的gnss高精度定位方法
Bahadur A study on the real-time code-based GNSS positioning with Android smartphones
CN113109840A (zh) 一种基于gnss接收机的电离层tec实时测量方法
CN113358017A (zh) 一种多站协同处理的gnss高精度变形监测方法
CN110146904B (zh) 一种适用于区域电离层tec的精确建模方法
Banville et al. Defining the basis of an integer-levelling procedure for estimating slant total electron content
Geng et al. A Robust Android Gnss Rtk Positioning Scheme Using Factor Graph Optimization
Li et al. A grid-based ionospheric weighted method for PPP-RTK with diverse network scales and ionospheric activity levels
CN104833993B (zh) 一种基于卫星和接收机频间偏差之和的北斗定位方法
Bahadur Real-time single-frequency precise positioning with Galileo satellites
CN104991265A (zh) 一种北斗卫星导航系统用户统一性定位方法
WO2023216766A1 (zh) 一种定位方法以及相关装置
CN117289312A (zh) 一种基于最优卫星子集的智能手机伪距差分定位定权方法
CN109444931B (zh) 一种静态伪距单点定位的方法及装置
Kuang et al. Galileo real-time orbit determination with multi-frequency raw observations
CN115308781A (zh) 基于bdgim辅助的相位平滑伪距高精度时间传递方法
Kashani et al. Towards instantaneous network-based RTK GPS over 100 km distance
CN112987048A (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210126