CN108845340A - 基于gnss-rtk的定位方法 - Google Patents
基于gnss-rtk的定位方法 Download PDFInfo
- Publication number
- CN108845340A CN108845340A CN201810556616.5A CN201810556616A CN108845340A CN 108845340 A CN108845340 A CN 108845340A CN 201810556616 A CN201810556616 A CN 201810556616A CN 108845340 A CN108845340 A CN 108845340A
- Authority
- CN
- China
- Prior art keywords
- solution
- fuzziness
- epoch
- observation
- float
- 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.)
- Pending
Links
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/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
-
- 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
Abstract
本发明属于卫星导航定位领域,涉及基于GNSS‑RTK的定位方法,包括以下步骤:步骤1,根据解析的原始星地观测值进行数据质量控制;步骤2,采用伪距观测值进行SPP定位,根据上一历元的RTD伪距残差和当前历元SPP定位计算的伪距残差进行RTD定位;步骤3,根据RTD定位结果作为初始值分别通过最小二乘、附加模糊度参数的Kalman滤波计算浮点解;步骤4,采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解。通过使用本发明,可以实现以下效果:能够快速固定模糊度,在复杂环境下具有较高的解算性能,提供复杂环境的准确、连续高精度位置信息;充分利用解算信息,实现模糊度快速固定,提供高精度、高可靠性的位置、速度信息。
Description
技术领域
本发明属于卫星导航定位领域,涉及基于GNSS-RTK的定位方法。
背景技术
GNSS:全球卫星导航系统,(Global Navigation Satellite System),也称为全球导航卫星系统,是能在地球表面或近地空间的任何地点为用户提供全天候的3 维坐标和速度以及时间信息的空基无线电导航定位系统。RTK(Real-time kinematic,实时动态)载波相位差分技术,是实时处理两个测量站载波相位观测 量的差分方法,将基准站采集的载波相位发给用户接收机,进行求差解算坐标。 随着人们对位置信息需求的不断增加,传统的高成本高精度定位模式已经无法满 足现实的需求,特别在建筑物、树木遮挡等复杂环境下定位时,存在精度、可靠 性差的问题。
发明内容
为解决上述问题,本发明提出定位精度高的基于GNSS-RTK的定位方法。
基于GNSS-RTK的定位方法,包括以下步骤:
步骤1,根据解析的原始星地观测值进行数据质量控制;
步骤2,采用伪距观测值进行SPP定位,根据上一历元的RTD伪距残差和 当前历元SPP定位计算的伪距残差进行RTD定位;
步骤3,根据RTD定位结果作为初始值分别通过最小二乘、附加模糊度参数 的Kalman滤波计算浮点解;
步骤4,采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解。
优选的,所述步骤1中根据解析的原始星地观测值进行数据质量控制包括以 下步骤:
通过观测数据预处理获取可靠的观测值;
利用数据探测法探测并改正粗差;
通过基于历元间二次差的周跳探测方法以及基于历元间差分的单站周跳探 测方法探测并剔除周跳。
优选的,所述通过观测数据预处理获取可靠的观测值包括:选取信噪比位于 设定范围内的观测值,通过伪距、载波精度信息和可用标识进行数据可用性判断, 并根据卫星连续跟踪、解算可用性选取观测数据。
优选的,所述利用数据探测法探测并改正粗差包括以下步骤:
构造线性化的伪距观测方程为:
V=L-A·X
其中,V:残差向量;L:观测值减计算值;A:设计矩阵;X:参数向量;
构造标准化残差作为粗差检测量:
QVV=P-1-A(ATPA)-1AT
其中,Vi:残差向量的元素;P:观测值的权阵;δ:单位权中误差;QViVi:残 差协因数矩阵QVV对角线元素;
如果标准化残差则存在粗差,标准化残差时,对权阵做 如下处理:
观测值减计算值L分为不包含粗差的l1和包含粗差向量S的l2,其维数分 别为n1和n2;则伪距观测方程变为:
设H=[0E]T,E为n2×n2单位阵,计算粗差向量S并将其改正:
PSS=HTPQVVPH。
优选的,所述通过基于历元间二次差的周跳探测方法以及基于历元间差 分的单站周跳探测方法探测并剔除周跳包括:先通过历元间二次差探测两周 以上的大周跳,然后通过基于历元间单差观测值的单站周跳探测方法探测与 剔除小周跳,所述通过历元间二次差探测两周以上的大周跳包括:
根据当前历元卫星的载波观测值与上一历元的载波观测值计算当 前历元的历元间载波观测值单差:
计算上一历元的历元间载波观测值单差:
计算当前历元的历元间载波观测值二次差:
当载波连续跟踪时,计算所有观测卫星的历元间载波观测值二次差
向量F中各元素集中分布,若某一元素超出数值 中心2周,则认为是大周跳。
优选的,所述步骤2中RTD定位包括以下步骤:
根据上一历元的RTD伪距残差和当前历元SPP定位计算的伪距残差计算得 到当前历元的观测值方差
其中,当前历元的观测值方差由以下模型进行计算:
其中,Vb:基准站SPP定位的伪距残差;Vr:流动站SPP定位的伪距残差,V0: 上一历元的RTD定位的伪距残差。
优选的,所述步骤3中根据RTD定位结果作为初始值分别通过最小二乘、 附加模糊度参数的Kalman滤波计算模糊度浮点解包括以下步骤:
构造双差观测方程:
其中,双差算子;Pi:伪距观测值;ρ:卫地距;I:电离层延迟;双差伪距观测噪声;φi:载波观测值;λi:载波波长;双差模糊度;
根据双差观测方程将流动站的位置、速度、加速度以及站间单差模糊度作为 待估参数,通过转换矩阵将站间单差模糊度转换成双差模糊度Kalman滤 波的参数矩阵和转换矩阵分别为:
Xk=[x y z vx vy vz ax ay az △N1 … △Nn]T
其中,(x,y,z):参考时刻的卫星位置;(vx,vy,vz):参考时刻的卫星运 动速度;(ax,ay,az):参考时刻的卫星运动加速度;
最小二乘法解算得到模糊度浮点解和验后残差σ,根据模糊度浮点解 和验后残差σ对Kalman滤波解算出的模糊度浮点解作如下变换:
优选的,所述步骤4中采用优化的部分模糊度固定策略将模糊度浮点解转换 为固定解,包括以下步骤:
根据浮点解方差对模糊度进行排序,模糊度固定按照浮点解精度顺序逐步进 行搜索;对双差模糊度参数和对应的方差-协方差矩阵进行降相关,得到相关性 较低的变换模糊度,降相关公式为:
其中,ZT:变换矩阵;模糊度浮点解;变换后的模糊度;模糊度 浮点解方差-协方差矩阵;变换后的模糊度浮点解方差-协方差矩阵,
矩阵元素全为整数,行列式的值为±1,变换后的搜索空间为:
其中,χ2:搜索空间;
搜索得到的整数解通过逆变换得到原始模糊度的整数解:
定义Ratio值对模糊度固定解进行检验:
其中,为次优解;为最优解;若满足上式则认为解算合格,否则保留 浮点解;
当连续若干个模糊度均固定,符合首次模糊度固定条件,保留模糊度固定解, 并设置固定解精度;对部分模糊度进行检验,若满足固定条件,将部分模糊度设 置为固定解。
通过使用本发明,可以实现以下效果:
(1)通过数据质量控制和残差处理,能够对原始数据进行有效改进,对数 据质量进行充分评估,为模糊度浮点解以及固定解提供可靠的原始观测数据和参 考依据,提高定位精度;
(2)通过最小二乘、附加模糊度参数的Kalman滤波计算模糊度浮点解,有 效提高收敛速度、防止滤波异常,提高浮点解的精度,进而缩短模糊度固定时间;
(3)充分利用解算信息,实现模糊度快速固定。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明的整体流程示意图;
图2是本发明的具体流程示意图。
具体实施方式
以下结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这 些实施例。
一种基于GNSS-RTK的定位方法,其主要思想包括:
(1)数据预处理和质量控制。选取可靠的卫星观测值、发现和修复伪距粗差、 探测和剔除载波周跳,提供可靠的伪距、载波观测值,保证模糊度的连续性,降 低或者防止固定解失锁;
(2)实时更新随机模型。RTD(实时伪距相对定位)随机模型根据SPP(伪距 单点定位)残差和上一历元的验后残差确定,浮点解双差伪距随机模型根据RTD 残差确定,保证了随机模型的准确度,提高了定位结果的精度;
(3)最小二乘原理和Kalman滤波确定浮点解。由最小二乘解算的模糊度浮点 解对Kalman滤波解进行改进,有效提高浮点解精度,提高滤波收敛效率,为进 一步固定模糊度做准备;
(4)模糊度固定。将模糊度浮点解按方差-协方差进行排序,采用逐步固定、 部分模糊度固定和保留策略,有效提高模糊度固定效率,事后滤波提高了浮点解 的精度,快速可靠地实现高精度定位信息,满足位置需求。
结合附图1,本实施例包括以下步骤:
步骤1,根据解析的原始星地观测值进行数据质量控制;
步骤2,采用伪距观测值进行SPP定位,根据上一历元的RTD伪距残差和 当前历元SPP定位计算的伪距残差进行RTD定位;
步骤3,根据RTD定位结果作为初始值分别通过最小二乘、附加模糊度参数 的Kalman滤波计算浮点解;
步骤4,采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解。
结合附图2,下面对各个步骤进行详细说明:
步骤1,根据解析的原始星地观测值进行数据质量控制。
基准站、流动站原始观测数据一般采用自定义协议或者标准协议进行编译, 通过对基准站原始观测数据和流动站原始观测数据解析得到观测值和卫星星历。 然而观测值和卫星星历存在误差较大值,同时一些观测值和卫星星历存在周跳等。
基准站差分数据一般采用RTCM标准协议进行编译。RTCM帧结构由同步码、 保留位、信息长度、数据信息和CRC组成,对于RTD需播发包含基准站位置信 息、GNSS各系统观测值的定义信息。流动站原始观测数据一般采用自定义的二 进制协议,需包含卫星观测数据以及数据精度指标。
卫星轨道信息通过历书和当地观测时间进行计算。对于CDMA(码分多址) 卫星定位系统,广播星历预报参数包括:
①时间参数:toe:星历参考时刻(秒);toc:钟差参考时刻(秒);IODE:星 历数据龄期;
②开普勒六参数:M0:toe时刻的平近点角(弧度);e:轨道偏心率;长 半轴的平方根Ω0:按参考时间计算的升交点赤经(弧度);i0:按参 考时间计算的轨道倾角(弧度);ω:近地点角距(弧度);
③轨道摄动九参数:△n:卫星平均运动速率与计算值之差(弧度);升交 点赤经变化率(弧度/秒);IDOT:轨道倾角变化率(弧度/秒);Cuc:纬度幅 角的余弦调和改正项的振幅(弧度);Cus:纬度幅角的正弦调和改正项的振幅 (弧度);Crc:轨道半径的余弦调和改正项的振幅(弧度);Crs:轨道半径的 正弦调和改正项的振幅(弧度);Cic:轨道倾角的余弦调和改正项的振幅(弧 度);Cis:轨道倾角的正弦调和改正项的振幅(弧度)。
对于FDMA(频分多址)卫星定位系统,广播星历预报参数包括:
①时间参数:toe:星历参考时刻(秒);IODE:星历数据龄期;
②位置参数:(x,y,z):参考时刻的卫星位置;(vx,vy,vz):参考时刻的卫 星运动速度;(ax,ay,az):参考时刻的卫星运动加速度。
具体的,步骤1中根据解析的原始星地观测值进行数据质量控制包括以 下步骤:
步骤1.1,观测数据预处理:通过观测数据预处理获取可靠的观测值。
首先根据观测数据的信噪比,选取信噪比位于合理范围的观测值;然后 根据数据可用标识,剔除不可用的数据;再利用数据输出的伪距、载波精度 信息,选取精度较高的数据;最后,根据卫星跟踪连续情况确定观测值是否 参与解算情况。
步骤1.2,伪距粗差探测与修复:利用数据探测法探测并改正粗差。
构造线性化的伪距观测方程为:
V=L-A·X(1)
其中:V:残差向量(米);L:观测值减计算值(米);A:设计矩阵;X: 参数向量。
构造标准化残差作为粗差检测量:
QVV=P-1-A(ATPA)-1AT(4)
其中,Vi:残差向量的元素;P:观测值的权阵;δ:单位权中误差;QViVi:残 差协因数矩阵QVV对角线元素。
如果标准化残差则存在粗差,标准化残差时,对权阵做 如下处理:
将公式(1)中L分为不包含粗差的l1和包含粗差向量S的l2,其维数分别为 n1和n2。
则(1)式变为:
设H=[0E]T,E为n2×n2单位阵,计算粗差向量S并将其改正,
PSS=HTPQVVPH(8)。
步骤1.3,载波周跳探测与剔除:先通过历元间二次差探测两周以上的大周 跳,然后通过基于历元间单差观测值的单站周跳探测方法探测与剔除小周跳。
具体的,通过历元间二次差探测两周以上的大周跳:
根据当前历元卫星的载波观测值与上一历元的载波观测值计算当前 历元的历元间载波观测值单差:
计算上一历元的历元间载波观测值单差:
计算当前历元的历元间载波观测值二次差:
当载波连续跟踪时,计算所有观测卫星的历元间载波观测值二次差向量F中各元素集中分布,若某一元素超出数值 中心2周,则认为是大周跳,数值中心通过分布集中的元素求平均值得出。
当载体处于低速运动时(运动状态不固定),载波观测值的二次差值取决于 运动状态,其数值与高次差法计算出的特征不一致,不服从于高斯白噪声分布, 但每一历元计算的二次差数值都集中分布在某个中心点附近,这是因为所有卫星 的观测值都和当前运动状态相关。根据该分布特点可确定载波是否发生周跳,具 体计算方法为:首先对载波观测值二次差值进行从小到大排序,选出该组数值的 中位数作为数值的中心点,计算所有数值到中心点的距离,超出2周的数值即认 为载波发生了大周跳。
具体的,通过基于历元间单差观测值的单站周跳探测方法探测与剔除小周跳:
构造历元间单差观测方程为:
其中,历元间单差观测值(周);ρSp:历元t的卫地距(米);ρSq:历 元t+1的卫地距(米);△t:接收机钟变化量(秒);c:光速(m/s);λ:载 波波长(米);△N:载波周跳值(周);△ε:观测噪声(米)。
待估参数包括三个坐标参数和一个接收机参数,若验后单位权中误差在毫米 量级,则不存在周跳,否则采用抗差最小二乘算法方程,逐步进行迭代计算,直 至前后两次解的坐标差值符合限差要求,等价权函数选取如下:
其中,表示标准化残差;P:权阵模型,采用高度角定权。
若标准化残差认为可能存在小周跳,然后在(9)中每次 增加一颗卫星的周跳估值作为待估参数,其余异常卫星的权值置为零,循环计算 每颗卫星周跳值,若周跳取整值不等于零则认为存在小周跳,剔除小周跳,否则, 将异常卫星归为正常卫星。
步骤2,采用伪距观测值进行SPP定位,根据上一历元的RTD伪距残差和 当前历元SPP定位计算的伪距残差进行RTD定位包括以下步骤:
步骤2.1,构造伪距双差观测方程为:
其中,双差算子;Pi:伪距观测值(米);ρ:卫地距(米);I:电离层 延迟(米);双差伪距观测噪声(米)。
根据获取的高质量观测值进行SPP定位解算,得到基准站和流动站SPP定 位计算的伪距残差。
步骤2.2,由上一历元的系统值和方差,通过Kalman滤波模型递推可得到当 前历元的状态估计,根据测量方程确定当前历元的系统值和方差。
其中,Kalman滤波模型为:
其中,I:单位矩阵;一步预测值;Pk,k-1:一步预测值方差-协方差矩阵;Tk-1:系统噪声驱动矩阵;Qk-1:系统噪声序列的方差阵;Kk:增益矩阵Lk:观 测值序列;待估参数;Hk:系数矩阵;Rk:测量噪声的方差阵;Xk:滤 波估值;Pk:滤波方差-协方差矩阵。
步骤2.3,流动站位置初始值根据滤波模型递推结果进行确定,通过Kalman 滤波模型实时更新伪距双差OMC(观测值减测量值)。设基准站和流动站SPP 定位的伪距残差为Vb和Vr,上一历元的RTD伪距残差为V0,当前历元的观测值 方差由以下模型进行计算:
在这一步骤中,首先进行SPP定位,若解算失败则不输出;若解算正确,则 进行RTD定位。根据上一历元的RTD伪距残差V0和当前历元SPP定位计算的伪 距残差Vb和Vr计算得到当前历元的观测值方差有利于提高RTD定位的精度 和稳定性。
步骤3,根据RTD定位结果作为初始值分别通过最小二乘解算、附加模糊度 参数的Kalman滤波计算模糊度浮点解包括以下步骤:
步骤3.1,构造双差观测方程:
其中,,双差算子;Pi:伪距观测值(米);ρ:卫地距(米);I:电离层 延迟(米);双差伪距观测噪声(米);φi:载波观测值(米);λi:载波 波长(米);双差模糊度(周)。
步骤3.2,根据双差观测方程将流动站的位置、速度、加速度以及站间单差 模糊度作为待估参数,通过转换矩阵将站间单差模糊度转换成双差模糊度 Kalman滤波的参数矩阵和转换矩阵分别为:
Xk=[x y z vx vy vz ax ay az △N1 … △Nn]T (15)
步骤3.3,和伪距相对定位相同,浮点解Kalman滤波观测方程同样更新双差 观测方程OMC。为提高滤波的数值稳定性,控制滤波发散,Pk矩阵改写为:
步骤3.4,根据Kalman滤波解算出的模糊度浮点解可能因为滤波收敛特征会 出现较大误差,通过最小二乘解算出的模糊度浮点解和验后残差σ,由Kalman 滤波解算出的模糊度浮点解做如下变换以提高精度:
步骤4,采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解,包 括以下步骤:
根据浮点解方差对模糊度进行排序,模糊度固定按照浮点解精度顺序逐步进 行搜索;对双差模糊度参数和对应的方差-协方差矩阵进行降相关,得到相关性 较低的变换模糊度,降相关公式为:
其中,ZT:变换矩阵;模糊度浮点解;变换后的模糊度;模糊度 浮点解方差-协方差矩阵;变换后的模糊度浮点解方差-协方差矩阵,
矩阵元素全为整数,行列式的值为±1,变换后的搜索空间为:
其中,χ2:搜索空间;
搜索得到的整数解通过逆变换得到原始模糊度的整数解:
定义Ratio值对模糊度固定解进行检验:
其中,为次优解;为最优解;若满足上式则认为解算正确,若不满足 上式则认为解算失败,保留浮点解,根据浮点解输出定位信息;
优化过程为:当连续若干个模糊度均固定,符合首次模糊度固定条件,保留 模糊度固定解,并设置固定解精度;对部分模糊度进行检验,若满足固定条件, 将部分模糊度设置为模糊度固定。该优化过程通过逐步固定并保留部分模糊度的 方法提高模糊度固定率,进一步提高系统在遮蔽、半遮蔽区域的定位的稳定性和 可靠性。
最后计算得到位置参数固定解,根据模糊度固定解、位置参数固定解输出定 位信息,位置参数固定解计算公式:
其中,位置参数浮点解;位置参数固定解;模糊度浮点解; 模糊度固定解;位置参数和模糊度参数协方差矩阵;模糊度 参数方差-协方差矩阵。
在本实施例中,对原始观测值进行数据预处理和质量控制。根据原始数据信 噪比、卫星高度角、连续跟踪情况等信息选取参与解算的卫星,对伪距粗差进行 探测和修复确保伪距数据质量,对载波观测值进行周跳探测与剔除,为高精度的 相对定位提供数据质量保障.
根据SPP和RTD解算结果对伪距数据质量进行实时评估。通过SPP计算的 伪距残差和上一历元的验后残差,确定RTD中双差伪距观测方程的随机模型, 根据RTD计算的双差伪距残差得出伪距残差,为浮点解的解算提供参考。
将站间单差模糊度作为参数通过Kalman滤波进行实时估计。Kalman滤波通 过对双差测量观测值进行实时更新提高浮点解精度,将单差模糊度作为参数进行 实时估计可以有效抵抗模糊度固定错误、避免参考卫星变化带来的复杂数据处理、 提高解算效率。
采用Kalman滤波和最小二乘原理确定浮点解。由于原始伪距观测值精度较 差,为减少收敛时间、防止滤波发散以及提高浮点解的精度,对Kalman滤波浮 点解进行合理控制,浮点解的结果根据Kalman滤波和最小二乘结果进行确定。
将浮点解转换成固定解。浮点解的精度越高,转换成固定解的成功率越高, 首先确定浮点解精度高的卫星观测值的模糊度,将其设定为固定值,通过Kalman 事后滤波提高剩余模糊度浮点解的精度,然后再逐步确定剩余的模糊度。
本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的 修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权 利要求书所定义的范围。
Claims (8)
1.基于GNSS-RTK的定位方法,其特征在于,包括以下步骤:
步骤1,根据解析的原始星地观测值进行数据质量控制;
步骤2,采用伪距观测值进行SPP定位,根据上一历元的RTD伪距残差和当前历元SPP定位计算的伪距残差进行RTD定位;
步骤3,根据RTD定位结果作为初始值分别通过最小二乘、附加模糊度参数的Kalman滤波计算浮点解;
步骤4,采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解。
2.根据权利要求1所述的基于GNSS-RTK的定位方法,其特征在于,所述步骤1中根据解析的原始星地观测值进行数据质量控制包括以下步骤:
通过观测数据预处理获取可靠的观测值;
利用数据探测法探测并改正粗差;
通过基于历元间二次差的周跳探测方法以及基于历元间差分的单站周跳探测方法探测并剔除周跳。
3.根据权利要求2所述的基于GNSS-RTK的定位方法,其特征在于,所述通过观测数据预处理获取可靠的观测值包括:选取信噪比位于设定范围内的观测值,通过伪距、载波精度信息和可用标识进行数据可用性判断,并根据卫星连续跟踪、解算可用性选取观测数据。
4.根据权利要求2所述的基于GNSS-RTK的定位方法,其特征在于,所述利用数据探测法探测并改正粗差包括以下步骤:
构造线性化的伪距观测方程为:
V=L-A·X
其中,V:残差向量;L:观测值减计算值;A:设计矩阵;X:参数向量;
构造标准化残差作为粗差检测量:
QVV=P-1-A(ATPA)-1AT
其中,Vi:残差向量的元素;P:观测值的权阵;δ:单位权中误差;残差协因数矩阵QVV对角线元素;
如果标准化残差则存在粗差,标准化残差时,对权阵做如下处理:
观测值减计算值L分为不包含粗差的l1和包含粗差向量S的l2,其维数分别为n1和n2;则伪距观测方程变为:
设H=[0 E]T,E为n2×n2单位阵,计算粗差向量S并将其改正:
PSS=HTPQVVPH。
5.根据权利要求2所述的基于GNSS-RTK的定位方法,其特征在于,所述通过基于历元间二次差的周跳探测方法以及基于历元间差分的单站周跳探测方法探测并剔除周跳包括:先通过历元间二次差探测两周以上的大周跳,然后通过基于历元间单差观测值的单站周跳探测方法探测与剔除小周跳,所述通过历元间二次差探测两周以上的大周跳包括:
根据当前历元卫星的载波观测值与上一历元的载波观测值计算当前历元的历元间载波观测值单差:
计算上一历元的历元间载波观测值单差:
计算当前历元的历元间载波观测值二次差:
当载波连续跟踪时,计算所有观测卫星的历元间载波观测值二次差向量F中各元素集中分布,若某一元素超出数值中心2周,则认为是大周跳。
6.根据权利要求1所述的基于GNSS-RTK的定位方法,其特征在于,所述步骤2中RTD定位包括以下步骤:
根据上一历元的RTD伪距残差和当前历元SPP定位计算的伪距残差计算得到当前历元的观测值方差
其中,当前历元的观测值方差由以下模型进行计算:
其中,Vb:基准站SPP定位的伪距残差;Vr:流动站SPP定位的伪距残差,V0:上一历元的RTD定位的伪距残差。
7.根据权利要求1所述的基于GNSS-RTK的定位方法,其特征在于,所述步骤3中根据RTD定位结果作为初始值分别通过最小二乘、附加模糊度参数的Kalman滤波计算浮点解包括以下步骤:
构造双差观测方程:
其中,双差算子;Pi:伪距观测值;ρ:卫地距;T:对流层延迟;I:电离层延迟;双差伪距观测噪声;φi:载波观测值;λi:载波波长;双差模糊度;双差载波观测噪声。
根据双差观测方程将流动站的位置、速度、加速度以及站间单差模糊度作为待估参数,通过转换矩阵将站间单差模糊度转换成双差模糊度Kalman滤波的参数矩阵和转换矩阵分别为:
Xk=[x y z vx vy vz ax ay az △N1 … △Nn]T
其中,(x,y,z):参考时刻的卫星位置;(vx,vy,vz):参考时刻的卫星运动速度;(ax,ay,az):参考时刻的卫星运动加速度;
最小二乘法解算得到模糊度浮点解和验后残差σ,根据模糊度浮点解和验后残差σ对Kalman滤波解算出的模糊度浮点解作如下变换:
8.根据权利要求1所述的基于GNSS-RTK的定位方法,其特征在于,所述步骤4中采用优化的部分模糊度固定策略将模糊度浮点解转换为固定解,包括以下步骤:
根据浮点解方差对模糊度进行排序,模糊度固定按照浮点解精度顺序逐步进行搜索;
对双差模糊度参数和对应的方差-协方差矩阵进行降相关,得到相关性较低的变换模糊度,降相关公式为:
其中,ZT:变换矩阵;模糊度浮点解;变换后的模糊度;模糊度浮点解方差-协方差矩阵;变换后的模糊度浮点解方差-协方差矩阵,矩阵元素全为整数,行列式的值为±1,变换后的搜索空间为:
其中,χ2:搜索空间;
搜索得到的整数解通过逆变换得到原始模糊度的整数解:
定义Ratio值对模糊度固定解进行检验:
其中,为次优解;为最优解;若满足上式则认为解算合格,否则保留浮点解;
当连续若干个模糊度均固定,符合首次模糊度固定条件,保留模糊度固定解,并设置固定解精度;对部分模糊度进行检验,若满足固定条件,将部分模糊度设置为固定解。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810556616.5A CN108845340A (zh) | 2018-06-01 | 2018-06-01 | 基于gnss-rtk的定位方法 |
US16/965,019 US11002860B2 (en) | 2018-06-01 | 2019-05-30 | GNSS-RTK-based positioning method |
PCT/CN2019/089207 WO2019228439A1 (zh) | 2018-06-01 | 2019-05-30 | 基于gnss-rtk的定位方法 |
EP19810061.2A EP3805802A4 (en) | 2018-06-01 | 2019-05-30 | POSITIONING PROCESS BASED ON GNSS-RTK |
AU2019278052A AU2019278052B2 (en) | 2018-06-01 | 2019-05-30 | GNSS-RTK-based positioning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810556616.5A CN108845340A (zh) | 2018-06-01 | 2018-06-01 | 基于gnss-rtk的定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108845340A true CN108845340A (zh) | 2018-11-20 |
Family
ID=64211179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810556616.5A Pending CN108845340A (zh) | 2018-06-01 | 2018-06-01 | 基于gnss-rtk的定位方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11002860B2 (zh) |
EP (1) | EP3805802A4 (zh) |
CN (1) | CN108845340A (zh) |
AU (1) | AU2019278052B2 (zh) |
WO (1) | WO2019228439A1 (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108196284A (zh) * | 2018-01-20 | 2018-06-22 | 中国人民解放军61540部队 | 一种进行星间单差模糊度固定的gnss网数据处理方法 |
CN110058281A (zh) * | 2019-04-29 | 2019-07-26 | 湖南国科微电子股份有限公司 | 动态定位方法及装置 |
WO2019228439A1 (zh) * | 2018-06-01 | 2019-12-05 | 浙江亚特电器有限公司 | 基于gnss-rtk的定位方法 |
CN110646820A (zh) * | 2019-09-20 | 2020-01-03 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN110737008A (zh) * | 2019-08-31 | 2020-01-31 | 广州海达安控智能科技有限公司 | 模糊度固定方法、装置及存储介质 |
CN111045052A (zh) * | 2019-10-14 | 2020-04-21 | 广东星舆科技有限公司 | 一种智能终端伪距差分定位及质量控制方法 |
CN111239787A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss动态卡尔曼滤波方法 |
CN111273328A (zh) * | 2018-12-05 | 2020-06-12 | 千寻位置网络有限公司 | 模糊度的固定方法及装置、定位系统 |
CN111665529A (zh) * | 2020-07-17 | 2020-09-15 | 广东星舆科技有限公司 | 伪距观测数据的处理方法、装置及可读介质 |
CN113267796A (zh) * | 2021-05-13 | 2021-08-17 | 中国人民解放军92859部队 | 一种双天线gnss、rtk定位及测向方法 |
CN113466903A (zh) * | 2021-08-16 | 2021-10-01 | 电子科技大学 | 一种顾及观测值系统误差的部分模糊度固定算法 |
CN115079238A (zh) * | 2022-08-23 | 2022-09-20 | 安徽交欣科技股份有限公司 | 基于rtk的公路交通智能精确定位系统和方法 |
CN115220078A (zh) * | 2022-08-24 | 2022-10-21 | 长沙金维信息技术有限公司 | 基于载波相位差分的gnss高精度定位方法及导航方法 |
CN115993623A (zh) * | 2023-03-24 | 2023-04-21 | 武汉大学 | 一种自适应选星方法、装置、设备及可读存储介质 |
CN115993620A (zh) * | 2021-10-19 | 2023-04-21 | 千寻位置网络有限公司 | 模糊度固定方法及其系统 |
CN116990841A (zh) * | 2023-06-25 | 2023-11-03 | 无锡卡尔曼导航技术有限公司南京技术中心 | 一种gnss形变监测数据质量控制方法、系统及装置 |
WO2024066028A1 (zh) * | 2022-09-29 | 2024-04-04 | 广东汇天航空航天科技有限公司 | 一种基于观测值外推的卫星定位方法和装置 |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110109166B (zh) * | 2019-04-30 | 2020-06-09 | 东南大学 | 一种快速得到高可靠性卫星定位整数解的方法 |
CN110927757B (zh) * | 2019-12-26 | 2022-03-04 | 广东星舆科技有限公司 | 卫星观测数据的质控方法、装置及定位装置 |
CN111157983B (zh) * | 2020-01-06 | 2022-02-15 | 南京鹰目电子科技有限公司 | 一种雷达目标跟踪方法 |
CN111352138B (zh) * | 2020-01-21 | 2023-02-17 | 北京眸星科技有限公司 | 卫星导航系统静态定位方法 |
CN111965676A (zh) * | 2020-07-16 | 2020-11-20 | 北京航空航天大学 | 一种加快卡尔曼滤波rtk浮点解收敛速度的方法 |
CN112924998B (zh) * | 2021-02-02 | 2024-03-19 | 武汉珈纬智能科技有限公司 | Gnss rtk实时自适应无缝换站方法、系统、终端 |
CN113050135B (zh) * | 2021-03-26 | 2022-09-13 | 同济大学 | 一种低成本gnss终端数据异常处理的方法 |
CN113124919B (zh) * | 2021-04-09 | 2021-12-07 | 中移(上海)信息通信科技有限公司 | 参考站的异常监测方法、装置及相关设备 |
CN113156477B (zh) * | 2021-04-15 | 2022-08-26 | 中国科学院精密测量科学与技术创新研究院 | 一种基于Android智能手机的高精度RTK定位方法 |
CN113296138B (zh) * | 2021-04-25 | 2022-12-13 | 北京遥测技术研究所 | 一种陆基导航高精度定位方法 |
CN113064186A (zh) * | 2021-05-18 | 2021-07-02 | 东南大学 | 基于双系统卫星差分的智能终端快速动态定位系统及方法 |
CN113568018A (zh) * | 2021-07-14 | 2021-10-29 | 广东汇天航空航天科技有限公司 | 定位信息的更新方法、装置、电子设备及存储介质 |
CN113805212B (zh) * | 2021-09-26 | 2023-11-07 | 南宁桂电电子科技研究院有限公司 | 一种自适应的gnss载波相位差分滑坡监测方法 |
CN114488238B (zh) * | 2021-11-15 | 2022-11-18 | 北京凯芯微科技有限公司 | Rtk定位方法及相关设备 |
CN114705156B (zh) * | 2022-02-28 | 2023-06-27 | 合肥工业大学 | 北斗短报文传输差分数据的电力塔筒沉降监测方法 |
CN114254568B (zh) * | 2022-02-28 | 2022-07-05 | 浙江国遥地理信息技术有限公司 | 基于人工智能决策树的gps遥感洪水预警方法 |
CN114674212B (zh) * | 2022-03-18 | 2023-11-21 | 建研地基基础工程有限责任公司 | 一种高精度基坑位移监测方法及监测系统 |
CN115201864A (zh) * | 2022-07-13 | 2022-10-18 | 涟漪位置(广州)科技有限公司 | 一种检测卫星钟差跳变的方法、装置、存储介质及设备 |
CN115390096B (zh) * | 2022-08-29 | 2023-04-25 | 浙江大学 | 基于全视角星载gnss接收系统的低轨卫星实时相对定轨方法 |
CN115343742B (zh) * | 2022-09-21 | 2024-04-02 | 北极星云空间技术股份有限公司 | 一种双星八频的gnss-rtk高维度模糊度快速解算方法 |
CN115728793A (zh) * | 2022-10-26 | 2023-03-03 | 同济大学 | 一种基于dia理论的精密单点定位粗差探测与处理方法 |
CN115902981B (zh) * | 2022-11-14 | 2024-01-30 | 中南大学 | 列车定位优化方法、系统及轨道交通车辆 |
CN115616643B (zh) * | 2022-12-16 | 2023-04-07 | 南京航空航天大学 | 一种城市区域建模辅助的定位方法 |
CN115688196B (zh) * | 2022-12-26 | 2023-07-18 | 萨科(深圳)科技有限公司 | 一种基于互联网平台订单大数据的在线数据处理方法 |
CN115877414B (zh) * | 2023-02-20 | 2023-04-28 | 广州导远电子科技有限公司 | 一种固定解坐标的检核方法、装置、存储介质及电子设备 |
CN116204756B (zh) * | 2023-04-28 | 2023-07-07 | 武汉大学 | 一种多分析中心精密站坐标产品综合方法及系统 |
CN117202342B (zh) * | 2023-08-04 | 2024-03-08 | 北京讯腾智慧科技股份有限公司 | 基于ppk的区域基准站网网络rtk定位服务评估方法及装置 |
CN116774266B (zh) * | 2023-08-23 | 2024-04-12 | 北京理工大学 | 一种基于rtk-ppp技术的大坝形变检测方法 |
CN116893433B (zh) * | 2023-09-08 | 2023-12-05 | 真点科技(北京)有限公司 | 一种实现跟踪站观测值预测的方法及装置 |
CN117289318B (zh) * | 2023-11-24 | 2024-02-20 | 北京北方联星科技有限公司 | 一种基于bds辅助的glonass频间偏差实时估计方法 |
CN117452463A (zh) * | 2023-12-22 | 2024-01-26 | 开普勒卫星科技(武汉)有限公司 | 一种适用于复杂环境下单频终端的周跳探测与修复方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030112177A1 (en) * | 2001-12-19 | 2003-06-19 | Hiroyuki Toda | Carrier-phase-based relative positioning device |
CN102096084A (zh) * | 2010-12-09 | 2011-06-15 | 东南大学 | 基于星间组合差分的精密单点定位方法 |
CN104570011A (zh) * | 2014-12-26 | 2015-04-29 | 中国电子科技集团公司第二十研究所 | 一种卫星导航相对定位装置及其载波相位周跳修复方法 |
CN106093991A (zh) * | 2016-06-21 | 2016-11-09 | 武汉大学 | 一种用于gnss定位的模糊度快速恢复方法及系统 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN107193028A (zh) * | 2017-03-29 | 2017-09-22 | 中国航空无线电电子研究所 | 基于GNSS的Kalman相对定位方法 |
CN107607973A (zh) * | 2017-09-12 | 2018-01-19 | 武汉大学 | Gnss网络rtk参考站间模糊度快速固定方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7228230B2 (en) * | 2004-11-12 | 2007-06-05 | Mitsubishi Denki Kabushiki Kaisha | System for autonomous vehicle navigation with carrier phase DGPS and laser-scanner augmentation |
CN106873017B (zh) * | 2017-02-21 | 2019-06-25 | 上海海积信息科技股份有限公司 | 模糊度固定方法 |
CN108051840B (zh) * | 2017-11-03 | 2021-07-16 | 中国航空无线电电子研究所 | 一种基于gnss的含约束ekf相对定位方法 |
US10948605B2 (en) * | 2017-12-21 | 2021-03-16 | Valeo Comfort And Driving Assistance | Broadcast and utilization of precise GNSS correction data |
CN108871179B (zh) * | 2018-05-07 | 2020-06-16 | 桂林电子科技大学 | 基于载波相位差分动静态融合的形变监测定位方法及装置 |
CN108845340A (zh) * | 2018-06-01 | 2018-11-20 | 浙江亚特电器有限公司 | 基于gnss-rtk的定位方法 |
-
2018
- 2018-06-01 CN CN201810556616.5A patent/CN108845340A/zh active Pending
-
2019
- 2019-05-30 WO PCT/CN2019/089207 patent/WO2019228439A1/zh unknown
- 2019-05-30 EP EP19810061.2A patent/EP3805802A4/en active Pending
- 2019-05-30 AU AU2019278052A patent/AU2019278052B2/en active Active
- 2019-05-30 US US16/965,019 patent/US11002860B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030112177A1 (en) * | 2001-12-19 | 2003-06-19 | Hiroyuki Toda | Carrier-phase-based relative positioning device |
CN102096084A (zh) * | 2010-12-09 | 2011-06-15 | 东南大学 | 基于星间组合差分的精密单点定位方法 |
CN104570011A (zh) * | 2014-12-26 | 2015-04-29 | 中国电子科技集团公司第二十研究所 | 一种卫星导航相对定位装置及其载波相位周跳修复方法 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN106093991A (zh) * | 2016-06-21 | 2016-11-09 | 武汉大学 | 一种用于gnss定位的模糊度快速恢复方法及系统 |
CN107193028A (zh) * | 2017-03-29 | 2017-09-22 | 中国航空无线电电子研究所 | 基于GNSS的Kalman相对定位方法 |
CN107607973A (zh) * | 2017-09-12 | 2018-01-19 | 武汉大学 | Gnss网络rtk参考站间模糊度快速固定方法及系统 |
Non-Patent Citations (3)
Title |
---|
岑新远: "粗差分类探测的选权迭代抗差估计方法", 《测绘》 * |
李盼: "GNSS精密单点定位模糊度快速固定技术和方法研究", 《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》 * |
祝会忠: "网络RTK整周模糊度快速解算", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108196284A (zh) * | 2018-01-20 | 2018-06-22 | 中国人民解放军61540部队 | 一种进行星间单差模糊度固定的gnss网数据处理方法 |
CN108196284B (zh) * | 2018-01-20 | 2021-07-27 | 中国人民解放军61540部队 | 一种进行星间单差模糊度固定的gnss网数据处理方法 |
WO2019228439A1 (zh) * | 2018-06-01 | 2019-12-05 | 浙江亚特电器有限公司 | 基于gnss-rtk的定位方法 |
US11002860B2 (en) | 2018-06-01 | 2021-05-11 | Zhejiang Yat Electrical Appliance Co., Ltd | GNSS-RTK-based positioning method |
CN111273328A (zh) * | 2018-12-05 | 2020-06-12 | 千寻位置网络有限公司 | 模糊度的固定方法及装置、定位系统 |
CN110058281A (zh) * | 2019-04-29 | 2019-07-26 | 湖南国科微电子股份有限公司 | 动态定位方法及装置 |
CN110737008B (zh) * | 2019-08-31 | 2023-12-12 | 广州市中海达测绘仪器有限公司 | 模糊度固定方法、装置及存储介质 |
CN110737008A (zh) * | 2019-08-31 | 2020-01-31 | 广州海达安控智能科技有限公司 | 模糊度固定方法、装置及存储介质 |
CN110646820A (zh) * | 2019-09-20 | 2020-01-03 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN110646820B (zh) * | 2019-09-20 | 2021-11-30 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN111045052A (zh) * | 2019-10-14 | 2020-04-21 | 广东星舆科技有限公司 | 一种智能终端伪距差分定位及质量控制方法 |
CN111045052B (zh) * | 2019-10-14 | 2022-08-09 | 广东星舆科技有限公司 | 一种智能终端伪距差分定位及质量控制方法 |
CN111239787A (zh) * | 2020-02-28 | 2020-06-05 | 同济大学 | 一种集群自主协同中的gnss动态卡尔曼滤波方法 |
CN111665529B (zh) * | 2020-07-17 | 2023-04-14 | 广东星舆科技有限公司 | 伪距观测数据的处理方法、装置及可读介质 |
CN111665529A (zh) * | 2020-07-17 | 2020-09-15 | 广东星舆科技有限公司 | 伪距观测数据的处理方法、装置及可读介质 |
CN113267796A (zh) * | 2021-05-13 | 2021-08-17 | 中国人民解放军92859部队 | 一种双天线gnss、rtk定位及测向方法 |
CN113466903A (zh) * | 2021-08-16 | 2021-10-01 | 电子科技大学 | 一种顾及观测值系统误差的部分模糊度固定算法 |
CN113466903B (zh) * | 2021-08-16 | 2022-03-25 | 电子科技大学 | 一种顾及观测值系统误差的部分模糊度固定算法 |
CN115993620A (zh) * | 2021-10-19 | 2023-04-21 | 千寻位置网络有限公司 | 模糊度固定方法及其系统 |
CN115993620B (zh) * | 2021-10-19 | 2024-03-15 | 千寻位置网络有限公司 | 模糊度固定方法及其系统 |
CN115079238B (zh) * | 2022-08-23 | 2023-10-03 | 安徽交欣科技股份有限公司 | 基于rtk的公路交通智能精确定位系统和方法 |
CN115079238A (zh) * | 2022-08-23 | 2022-09-20 | 安徽交欣科技股份有限公司 | 基于rtk的公路交通智能精确定位系统和方法 |
CN115220078A (zh) * | 2022-08-24 | 2022-10-21 | 长沙金维信息技术有限公司 | 基于载波相位差分的gnss高精度定位方法及导航方法 |
WO2024066028A1 (zh) * | 2022-09-29 | 2024-04-04 | 广东汇天航空航天科技有限公司 | 一种基于观测值外推的卫星定位方法和装置 |
CN115993623A (zh) * | 2023-03-24 | 2023-04-21 | 武汉大学 | 一种自适应选星方法、装置、设备及可读存储介质 |
CN116990841A (zh) * | 2023-06-25 | 2023-11-03 | 无锡卡尔曼导航技术有限公司南京技术中心 | 一种gnss形变监测数据质量控制方法、系统及装置 |
CN116990841B (zh) * | 2023-06-25 | 2024-01-23 | 无锡卡尔曼导航技术有限公司南京技术中心 | 一种gnss形变监测数据质量控制方法、系统及装置 |
Also Published As
Publication number | Publication date |
---|---|
US11002860B2 (en) | 2021-05-11 |
US20210072406A1 (en) | 2021-03-11 |
AU2019278052B2 (en) | 2022-04-07 |
EP3805802A1 (en) | 2021-04-14 |
AU2019278052A1 (en) | 2020-11-26 |
EP3805802A4 (en) | 2022-02-23 |
WO2019228439A1 (zh) | 2019-12-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108845340A (zh) | 基于gnss-rtk的定位方法 | |
CN107422354B (zh) | 一种模糊度固定的ppp/sins紧组合定位定姿方法 | |
CN108317949B (zh) | 一种rtk高精度差分定位形变监测系统及方法 | |
CN110275186B (zh) | Leo卫星增强的gnss电离层归一化与融合建模方法 | |
CN108169774A (zh) | 支持rtppp和rtk的多模gnss单频周跳探测与修复方法 | |
CN105158783B (zh) | 一种实时动态差分定位方法及其设备 | |
CN103837879B (zh) | 基于北斗系统民用载波相位组合实现高精度定位的方法 | |
CN109709591A (zh) | 一种面向智能终端的gnss高精度定位方法 | |
CN108196279A (zh) | 基于实时数据流的卫星钟差解算及预报方法 | |
CN109116394A (zh) | 一种适用于不同长度基线的实时动态定位方法 | |
CN108549095A (zh) | 一种区域cors网非差并行增强方法及系统 | |
CN110031879B (zh) | 模糊度域信息整合的高精度后处理定位方法及系统 | |
CN112285745B (zh) | 基于北斗三号卫星导航系统的三频模糊度固定方法及系统 | |
CN105891860A (zh) | 一种基于误差分离模式的gnss区域伪距差分增强定位方法 | |
TW201445168A (zh) | 一種接收機和衛星定位及測速方法 | |
CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
Pan et al. | A new approach for optimising GNSS positioning performance in harsh observation environments | |
CN110426717A (zh) | 一种基于三维地图辅助gnss的协同定位方法及系统、定位设备、存储介质 | |
CN110646823A (zh) | 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法 | |
Liu et al. | Glio: Tightly-coupled gnss/lidar/imu integration for continuous and drift-free state estimation of intelligent vehicles in urban areas | |
CN109613582A (zh) | 一种车载实时单频米级伪距定位方法 | |
Xu et al. | A unified cycle-slip, multipath estimation, detection and mitigation method for VIO-aided PPP in urban environments | |
Luo et al. | Relative positioning of multiple moving platforms using GPS | |
Zhao et al. | Signal quality and positioning performance of GPS/BDS-3/GLONASS/Galileo in polar regions | |
Chen et al. | Undifferenced zenith tropospheric modeling and its application in fast ambiguity recovery for long-range network RTK reference stations |
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 | ||
AD01 | Patent right deemed abandoned | ||
AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20221104 |