CN110275192A - 一种基于智能手机的高精度单点定位方法与装置 - Google Patents
一种基于智能手机的高精度单点定位方法与装置 Download PDFInfo
- Publication number
- CN110275192A CN110275192A CN201910428335.6A CN201910428335A CN110275192A CN 110275192 A CN110275192 A CN 110275192A CN 201910428335 A CN201910428335 A CN 201910428335A CN 110275192 A CN110275192 A CN 110275192A
- Authority
- CN
- China
- Prior art keywords
- satellite
- observation
- smart phone
- gnss
- point positioning
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 claims description 20
- 230000003068 static effect Effects 0.000 claims description 8
- 239000005433 ionosphere Substances 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 239000005436 troposphere Substances 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 230000006399 behavior Effects 0.000 claims description 2
- 230000004807 localization Effects 0.000 abstract 2
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 244000062793 Sorghum vulgare Species 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 235000019713 millet Nutrition 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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(Global Navigation 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)中,所削弱的误差影响包括对接收到观测值使用精密星历、精密钟差文件消除轨道误差和卫星钟差,使用电离层格网文件削弱电离层延迟,以及将多路径效应视作观测噪声,所形成的非组合观测值模型为:
式中,上标g、e、c分别表示GPS系统、Galileo系统和BDS系统,下标i、j、k表示第i、j、k颗卫星,左侧P、Φ分别为伪距和载波观测值,ρ为手机终端到卫星的距离,c为光速,分别为伪距和载波观测值的移动端钟差,dtrop为对流层延迟,为载波整周模糊度,∈P、∈Φ分别为伪距和载波相位观测值的残差;和分别表示Galileo和BDS系统与GPS系统之间的时间偏差。
在优选的实施方案中,步骤(3)中,根据卫星高度角确定各卫星观测值权重的定权方案具体为:
当卫星高度角小于10°时,该卫星观测值权重为0;
当卫星高度角大于10°时,该卫星观测值权重为sinE,E为卫星高度角。
在优选的实施方案中,步骤(4)中,改进的静态卡尔曼滤波具体为:
对于单个历元,根据步骤(2)所述的非组合观测值模型,一颗卫星的观测方程数为2个,设某历元观测到GPS卫星n1颗,Galileo卫星n2颗,BDS卫星n3颗,满足n1+n2+n3≥8;此时卡尔曼滤波中待估参数向量为:
式中,x、y和z为智能手机位置参数;则观测方程总数为2×(n1+n2+n3)个,待估参数为n1+n2+n3+8个,多余观测数为n1+n2+n3-8个;
卡尔曼滤波的观测方程系数阵H为:
矩阵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原始数据表
表2 GnssClock类主要GNSS原始数据
然后根据GNSS原始数据中的时间数据解算出伪距,其实质是通过信号传播时间差值计算伪距:
P=(tRx-tTx)·c (1)
其中:P表示伪距,tRx表示接收机接收到信号的时刻,tTx表示接收到卫星发射信号的时刻,c为光速。tTx可通过getReceivedSvTimeNanos()方法直接获取,数值为本GPS周内的秒数。
其中表示从1980.6.1到当前时间的总秒数,GNSSweek表示所有整GPS周秒数和,TimeNanos为接收机本地时钟值getTimeNanos(),TimeOffsetNanos为本地时钟偏移值getTimeOffsetNanos(),FullBiasNanos为接收机本地时钟总偏差getFullBiasNanos(),BiasNanos为接收机本地钟差getBiasNanos()。
第二步,在进行数据预处理削弱部分误差影响之后,根据改进的估计双钟差的PPP算法,对原始观测值形成非组合观测值模型。
前期实验将智能手机和测地型接收机置于相同地点进行同步观测,收集原始观测值,并将智能手机原始观测值文件转换为Rinex格式文件。图2为相同时间相同卫星测地型接收机的原始观测值和智能手机(华为P10手机)的观测值对比。图中,接收机载波与伪距观测值重合在一起,实际上二者之差为载波相位的整周模糊度×载波波长。测地型接收机观测值在一段时间之后会有一个大幅度的跳动,这是接收机钟跳引起的,不影响原始观测值的定位解算。而可以明显地观察出手机伪距与载波相位观测值相互之间的差值不固定,这个性质与测地型接收机是不相同的,影响手机终端载波相位数据的使用。
对同一设备不同卫星原始观测值进行分析,发现同一时刻不同卫星的伪距变化率和载波观测值变化率之间的差异是相同的,故在定位时认为手机终端伪距观测值和载波观测值具有不同的钟差。
根据以上性质改进定位方案,在参数求解时估计两个钟差参数。以GPS系统为例,单频观测方程为:
式中,上标g表示GPS系统,下标i表示第i颗卫星,为伪距和载波观测值,为手机端到卫星的距离,c为光速,分别为伪距和载波观测值的钟差,dTg为卫星钟差,为对流层延迟,为卫星轨道误差,为电离层延迟,分别为伪距和载波多路径效应延迟,为载波整周模糊度,分别伪距和载波相位观测值的残差。
实时下载IGS(International GNSS Service,国际GNSS服务组织)提供的预测精密星历、钟差文件和预测电离层格网文件,消除卫星轨道误差和钟差,削弱电离层延迟误差,并通过模型改正消除相对论效应和地球自转等误差影响,将多路径效应视作观测噪声后,简化观测值方程,得到非组合观测值模型:
式中,左侧为各个GNSS系统的伪距和载波观测值,上标e、c分别表示Galileo系统和BDS系统,下标j、k表示第j、k颗卫星,和分别表示Galileo和BDS系统与GPS系统之间的时间偏差。
第三步,确定各卫星观测值权重。本方法根据卫星高度角确定各卫星观测值权重,定权方案具体为:
式中,W为权重,E为手机端到卫星连线的高度角。
第四步,使用卡尔曼滤波进行参数求解,得到高精度定位结果。
对观测值形成第二步所示的非组合模型,此时,单个历元一颗卫星的观测方程数为2个,假设某历元观测到GPS卫星n1颗,Galileo卫星n2颗,BDS卫星n3颗,此时卡尔曼滤波中待估参数向量为:
式中,x、y和z为智能手机位置参数。
则观测方程总数为2×(n1+n2+n3)个,待估参数为n1+n2+n3+8个,多余观测数为n1+n2+n2-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:
式中,Hk为观测方程的系数阵。相对于常规的静态卡尔曼滤波,本方法中的Hk需要进行相应的修改:
Hk行数为2×(n1+n2+n3),列数为n1+n2+n2+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+n2+8列为载波相位模糊度系数,第2×i行,i+8列为1,i为卫星序号。
式中,xs、ys和zs为卫星坐标(通过星历文件获得),x0、y0和z0为智能手机近似坐标(通过伪距单点定位获得),ρ0为手机端到卫星的近似距离(通过卫星坐标和手机近似坐标计算)。
式(16)中Rk为观测噪声方差矩阵:
式中,σ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实施例各时段定位中误差
定义收敛时间为,从定位开始时刻到N和E方向定位误差均小于1m,且后续历元误差不再超过1m的时刻为止经过的时间。统计5个时段的收敛时间如表4所示,各时段均能在30s之内实现收敛,说明本方法能够应用于实时定位,提供低延迟的高精度智能手机定位结果。
表4实施例各时段收敛时间
基于相同的发明构思,本发明实施例公开的一种基于智能手机的高精度单点定位装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述计算机程序被加载至处理器时实现所述的一种基于智能手机的高精度单点定位方法。
Claims (6)
1.一种基于智能手机的高精度单点定位方法,其特征在于,包括以下步骤:
(1)获取智能手机导航定位模块的GNSS伪距和载波相位原始观测值;
(2)在进行数据预处理削弱部分误差影响之后,根据改进的估计双钟差的精密单点定位方法,对原始观测值形成非组合观测值模型;
(3)根据卫星高度角确定各卫星观测值权重;
(4)使用改进的静态卡尔曼滤波进行计算,得到高精度的单点定位结果。
2.根据权利要求1所述的基于智能手机的高精度单点定位方法,其特征在于,步骤(1)中,利用智能手机操作系统提供的基于位置服务的API直接获取部分GNSS原始数据,包括时间数据、GNSS系统类型和载波相位观测值,然后根据GNSS原始数据中的时间数据通过信号传播时间差值解算出伪距。
3.根据权利要求1所述的基于智能手机的高精度单点定位方法,其特征在于,步骤(2)中,所削弱的误差影响包括对接收到观测值使用精密星历、精密钟差文件消除轨道误差和卫星钟差,使用电离层格网文件削弱电离层延迟,以及将多路径效应视作观测噪声,所形成的非组合观测值模型为:
式中,上标g、e、c分别表示GPS系统、Galileo系统和BDS系统,下标i、j、k表示第i、j、k颗卫星,左侧P、Φ分别为伪距和载波观测值,ρ为手机终端到卫星的距离,c为光速,分别为伪距和载波观测值的移动端钟差,dtrop为对流层延迟,为载波整周模糊度,∈P、∈Φ分别为伪距和载波相位观测值的残差;和分别表示Galileo和BDS系统与GPS系统之间的时间偏差。
4.根据权利要求1所述的基于智能手机的高精度单点定位方法,其特征在于,步骤(3)中,根据卫星高度角确定各卫星观测值权重的定权方案具体为:
当卫星高度角小于10°时,该卫星观测值权重为0;
当卫星高度角大于10°时,该卫星观测值权重为sin E,E为卫星高度角。
5.根据权利要求1所述的基于智能手机的高精度单点定位方法,其特征在于,步骤(4)中,改进的静态卡尔曼滤波具体为:
对于单个历元,根据步骤(2)所述的非组合观测值模型,一颗卫星的观测方程数为2个,设某历元观测到GPS卫星n1颗,Galileo卫星n2颗,BDS卫星n3颗,满足n1+n2+n3≥8;此时卡尔曼滤波中待估参数向量为:
式中,x、y和z为智能手机位置参数;则观测方程总数为2×(n1+n2+n3)个,待估参数为n1+n2+n3+8个,多余观测数为n1+n2+n3-8个;
卡尔曼滤波的观测方程系数阵H为:
矩阵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)确定。
6.一种基于智能手机的高精度单点定位装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述计算机程序被加载至处理器时实现根据权利要求1-5任一项所述的一种基于智能手机的高精度单点定位方法。
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 true CN110275192A (zh) | 2019-09-24 |
CN110275192B 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) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110687556A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 一种适用于laas的多径误差模型化方法 |
CN111045052A (zh) * | 2019-10-14 | 2020-04-21 | 广东星舆科技有限公司 | 一种智能终端伪距差分定位及质量控制方法 |
CN111060945A (zh) * | 2020-01-13 | 2020-04-24 | 东南大学 | 一种gnss/5g紧组合融合定位方法与装置 |
CN111796313A (zh) * | 2020-06-28 | 2020-10-20 | 中国人民解放军63921部队 | 卫星定位方法及装置、电子设备、存储介质 |
WO2020233158A1 (zh) * | 2019-05-22 | 2020-11-26 | 东南大学 | 一种基于智能手机的高精度单点定位方法与装置 |
CN113156477A (zh) * | 2021-04-15 | 2021-07-23 | 中国科学院精密测量科学与技术创新研究院 | 一种基于Android智能手机的高精度RTK定位方法 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
CN113267793A (zh) * | 2021-05-26 | 2021-08-17 | 中国电子科技集团公司第五十四研究所 | 一种基于外部增强信息的gbas对流层参数生成方法 |
CN114488227A (zh) * | 2022-01-26 | 2022-05-13 | 西南交通大学 | 一种基于空间相关性的多路径误差改正方法 |
CN114640950A (zh) * | 2022-03-06 | 2022-06-17 | 南京理工大学 | 基于Android源生GPS定位API的移动设备定位方法及系统 |
CN116184441A (zh) * | 2022-12-20 | 2023-05-30 | 北京航空航天大学 | 一种精密大气改正增强的智能终端ppp瞬时亚米级定位方法 |
CN116819585A (zh) * | 2023-08-31 | 2023-09-29 | 长沙金维信息技术有限公司 | 基于非线性优化的gnss单点定位方法及导航方法 |
CN117647830A (zh) * | 2024-01-29 | 2024-03-05 | 中国测绘科学研究院 | 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11650324B2 (en) * | 2020-05-01 | 2023-05-16 | Deere & Company | Navigation apparatus and method in which measurement quantization errors are modeled as states |
CN116678421B (zh) * | 2023-06-12 | 2024-01-23 | 深圳沧穹科技有限公司 | 基于多模块ble发射装置的多源融合定位方法及系统 |
Citations (6)
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)
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 |
US9128176B2 (en) * | 2007-06-22 | 2015-09-08 | Trimble Navigation Limited | GNSS signal processing to estimate MW biases |
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 | 东南大学 | 一种基于智能手机的高精度单点定位方法与装置 |
-
2019
- 2019-05-22 CN CN201910428335.6A patent/CN110275192B/zh not_active Expired - Fee Related
-
2020
- 2020-01-16 US US17/051,428 patent/US11709281B2/en active Active
- 2020-01-16 WO PCT/CN2020/072411 patent/WO2020233158A1/zh active Application Filing
Patent Citations (6)
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)
Title |
---|
王延染等: "广域实时亚米级单点定位算法", 《山东大学学报(工学版)》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020233158A1 (zh) * | 2019-05-22 | 2020-11-26 | 东南大学 | 一种基于智能手机的高精度单点定位方法与装置 |
US11709281B2 (en) | 2019-05-22 | 2023-07-25 | Southeast University | High-precision point positioning method and device based on smartphone |
CN111045052A (zh) * | 2019-10-14 | 2020-04-21 | 广东星舆科技有限公司 | 一种智能终端伪距差分定位及质量控制方法 |
CN110687556B (zh) * | 2019-11-04 | 2021-06-22 | 中国电子科技集团公司第五十四研究所 | 一种适用于laas的多径误差模型化方法 |
CN110687556A (zh) * | 2019-11-04 | 2020-01-14 | 中国电子科技集团公司第五十四研究所 | 一种适用于laas的多径误差模型化方法 |
CN111060945A (zh) * | 2020-01-13 | 2020-04-24 | 东南大学 | 一种gnss/5g紧组合融合定位方法与装置 |
CN111796313A (zh) * | 2020-06-28 | 2020-10-20 | 中国人民解放军63921部队 | 卫星定位方法及装置、电子设备、存储介质 |
CN111796313B (zh) * | 2020-06-28 | 2023-07-21 | 中国人民解放军63921部队 | 卫星定位方法及装置、电子设备、存储介质 |
CN113156477B (zh) * | 2021-04-15 | 2022-08-26 | 中国科学院精密测量科学与技术创新研究院 | 一种基于Android智能手机的高精度RTK定位方法 |
CN113156477A (zh) * | 2021-04-15 | 2021-07-23 | 中国科学院精密测量科学与技术创新研究院 | 一种基于Android智能手机的高精度RTK定位方法 |
CN113204042A (zh) * | 2021-05-21 | 2021-08-03 | 北京交通大学 | 一种基于精密单点定位的多星座联合列车定位方法 |
CN113267793A (zh) * | 2021-05-26 | 2021-08-17 | 中国电子科技集团公司第五十四研究所 | 一种基于外部增强信息的gbas对流层参数生成方法 |
CN113267793B (zh) * | 2021-05-26 | 2022-05-06 | 中国电子科技集团公司第五十四研究所 | 一种基于外部增强信息的gbas对流层参数生成方法 |
CN114488227A (zh) * | 2022-01-26 | 2022-05-13 | 西南交通大学 | 一种基于空间相关性的多路径误差改正方法 |
CN114488227B (zh) * | 2022-01-26 | 2023-10-20 | 西南交通大学 | 一种基于空间相关性的多路径误差改正方法 |
CN114640950A (zh) * | 2022-03-06 | 2022-06-17 | 南京理工大学 | 基于Android源生GPS定位API的移动设备定位方法及系统 |
CN116184441A (zh) * | 2022-12-20 | 2023-05-30 | 北京航空航天大学 | 一种精密大气改正增强的智能终端ppp瞬时亚米级定位方法 |
CN116184441B (zh) * | 2022-12-20 | 2023-10-20 | 北京航空航天大学 | 一种精密大气改正增强的智能终端ppp瞬时亚米级定位方法 |
CN116819585A (zh) * | 2023-08-31 | 2023-09-29 | 长沙金维信息技术有限公司 | 基于非线性优化的gnss单点定位方法及导航方法 |
CN116819585B (zh) * | 2023-08-31 | 2023-12-29 | 长沙金维信息技术有限公司 | 基于非线性优化的gnss单点定位方法及导航方法 |
CN117647830A (zh) * | 2024-01-29 | 2024-03-05 | 中国测绘科学研究院 | 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法 |
CN117647830B (zh) * | 2024-01-29 | 2024-05-07 | 中国测绘科学研究院 | 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
US11709281B2 (en) | 2023-07-25 |
US20220155465A1 (en) | 2022-05-19 |
CN110275192B (zh) | 2021-01-26 |
WO2020233158A1 (zh) | 2020-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110275192A (zh) | 一种基于智能手机的高精度单点定位方法与装置 | |
CN105849589B (zh) | 全球导航卫星系统、定位终端、定位方法以及记录介质 | |
CN107949795A (zh) | 用于协作式全球导航卫星系统(gnss)诊断的方法和系统 | |
CN105044747B (zh) | 一种基于多星共视和滤波的时间同步装置及其方法 | |
CN109343089B (zh) | 定位设备的性能测试设备、测试方法及其装置 | |
CN101680944A (zh) | 部分搜索载波相位整数模糊度解算 | |
CN107807373A (zh) | 基于移动智能终端的gnss高精度定位方法 | |
CN109085619B (zh) | 多模gnss系统的定位方法及装置、存储介质、接收机 | |
CN102033236A (zh) | 一种卫星导航位置速度联合估计方法 | |
CN103235321A (zh) | Gps单站伪距定位精密授时方法 | |
CN109633703B (zh) | 一种应对遮挡场景的北斗导航无源定位方法 | |
CN106383355A (zh) | 基于高精度时间的gnss时差监测方法及系统 | |
CN103592657A (zh) | 一种基于钟差辅助的低可见星下单模raim实现方法 | |
CN110320536A (zh) | 卫星定位参数校准方法、装置、终端设备及存储介质 | |
CN105388496B (zh) | 基于gps的交通应用脆弱性检测系统及其方法 | |
CN117289312A (zh) | 一种基于最优卫星子集的智能手机伪距差分定位定权方法 | |
CN113917509B (zh) | 一种双差模糊度固定方法、设备以及可读存储介质 | |
CN111045064B (zh) | 一种cors系统数据计算的方法及装置 | |
Kashani et al. | The double difference effect of ionospheric correction latency on instantaneous ambiguity resolution in long-range RTK | |
CN114114343A (zh) | 用于gnss的测量方法、系统及rtk接收机 | |
Rocken et al. | Testing a new network RTK software system | |
He et al. | Performance evaluation for bds-2+ bds-3 combined precision point positioning | |
CN111538045A (zh) | 一种星载导航接收机在轨精度预先评估方法 | |
Zhilinskiy | GLONASS Satellite Pseudorange Errors Mitigation Using Gradient Boosting Machine | |
Wielgosz et al. | Application of predictive regional ionosphere model to medium range RTK positioning |
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 |
Granted publication date: 20210126 |
|
CF01 | Termination of patent right due to non-payment of annual fee |