CN111142133B - 基于多个连续运行参考站的后处理定位方法和系统 - Google Patents
基于多个连续运行参考站的后处理定位方法和系统 Download PDFInfo
- Publication number
- CN111142133B CN111142133B CN201811303873.4A CN201811303873A CN111142133B CN 111142133 B CN111142133 B CN 111142133B CN 201811303873 A CN201811303873 A CN 201811303873A CN 111142133 B CN111142133 B CN 111142133B
- Authority
- CN
- China
- Prior art keywords
- data
- station
- continuously operating
- operating reference
- reference stations
- 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
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/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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- 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
本发明提供了一种基于多个连续运行参考站的后处理定位方法及系统，所述方法包括以下步骤：服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；对连续运行参考站数据进行处理，计算参考站差分改正数；基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。本发明解决了数据中断RTK长时间初始化不能获得高精度位置的问题，以及观测条件不好时解算结果不准确的问题。
Description
技术领域
本发明涉及定位技术领域，具体涉及一种基于多个连续运行参考站的后处理定位方法及系统。
背景技术
RTK(Real-time kinematic，实时动态)技术利用参考站的GNSS(GlobalNavigation Satellite System，全球导航卫星系统)观测值和用户接收到的GNSS观测值进行差分，解算得到用户端的位置。用户端在接收卫星信号的同时接收网络播发的参考站数据，能够快速实时的获得厘米级的定位精度。常规的RTK技术随着用户距离参考站距离的增大，其定位精度也会有所下降。
网络RTK是利用设置在区域内的连续运行参考站接收GNSS数据计算各个基线上的电离层，对流层误差，然后根据各个基线上的误差值内插得到区域内各个点上的电离层，对流层等误差，进而对参考站的观测值进行改正得到某一点上的观测值。再将该观测值发送到用户端进行常规的RTK解算。
RTK服务中只能按照时间顺序的先后进行滤波，解算用户的位置。实际应用中由于用户所处的环境动态的变化，在部分场景无法得到高精度的位置，因此现有技术存在以下技术问题：
1、传统RTK定位技术适用于距离参考站15km以内的场景，随着距离的增加受到大气延迟误差的影响定位精度将迅速减弱。
2、实时场景下用户流动站接收参考站差分改正数但不会保存，使得历史数据无法回溯。
3、实时情况下RTK解算仅能够按照时间的前后顺序进行，单向的滤波会出现发生周跳后需要再次初始化RTK。
4.本发明利用多个连续运行参考站对流动站数据解算，扩展RTK定位范围的同时利用历史信息提高了定位的精度。
发明内容
本发明公开了一种基于GNSS动态数据后处理的方法，利用采集到的用户数据和多个参考站数据，解决了数据中断RTK长时间初始化不能获得高精度位置的问题，以及观测条件不好时解算结果不准确的问题。
本发明采用的技术方案如下：
一种基于多个连续运行参考站的后处理定位方法，所述方法包括以下步骤：
服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
对连续运行参考站数据进行处理，计算参考站差分改正数；
基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。
进一步地，所述用户流动站数据和连续运行参考站数据是指GNSS的观测数据，包含每一历元的时间和每一颗卫星的载波相位与伪距的观测值。
进一步地，所述用户流动站数据是在动态或静态下采集，所述连续运行参考站数据是在静态下采集。
进一步地，所述对流动站原始数据进行分割具体包括以下步骤：
从用户流动站数据的第一个时刻起向空的经纬度序列中加入流动站的坐标；
分别计算最大纬度最小纬度之差对应的距离长度ΔB和最大经度最小经度之差对应的距离长度ΔD；
根据ΔB或者ΔD决定是否对流动站原始数据进行分割。
进一步地，如果ΔB或者ΔD大于阈值，则对流动站原始数据进行分割，将经纬度序列对应的时间段作为一个格网进行划分，并分别记录存储经度，纬度的平均值作为一个格网点的坐标等待后续使用。
进一步地，对于未划分的数据再次进行划分，直至最后一个时间段数据均符合流动站原始数据分割的规则。
进一步地，所述对连续运行参考站数据进行处理具体步骤如下：根据各个格网的时间段与格网点的坐标，选取对应的连续运行参考站接收GNSS观测数据和广播星历数据，计算对应的连续运行参考站差分改正数。
进一步地，所述选取对应的连续运行参考站的具体步骤如下：
以格网的位置为起始，搜索符合要求的所有连续运行参考站；
对各个连续运行参考站的数据进行数据质量评估，保留评估合格的连续运行参考站；
计算上一步中各个连续运行参考站对应于格网中心的方位角；
分别在0°至90°、90°至180°、180°至270°、270°至360°的四个区域内找到各区域内距离中心点最近的4个连续运行参考站，然后在剩余的所有连续运行参考站中找到2个距离格网中心最近的连续运行参考站；
如果一个区域缺少连续运行参考站，则就近选择1个连续运行参考站作为补充。
进一步地，根据数据集对每一个时间段的数据进行精准位置解算，所述数据集包括分割完成后的用户流动站数据、连续运行参考站差分改正数和广播星历数据，所述精准位置解算步骤如下：
对每一个时间段的数据进行正向滤波解算，从第一个历元开始利用正向滤波的方法解算各个历元的位置和模糊度参数；
对每一个时间段的数据进行反向滤波解算，从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数；
正向滤波解算和反向滤波解算完成后，对各个时间段内同一历元的正向滤波结果和反向滤波结果进行联合处理。
进一步地，所述正向滤波解算过程中第k个历元时在获得浮点解后进行模糊度固定及检核，如果固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，然后继续进行第k+1历元的解算，直至当前时间段数据解算完毕，进入下一时间段的正向滤波解算。
进一步地，所述反向滤波解算过程中第k个历元时在获得浮点解后进行模糊度固定及检核，如果固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，然后继续进行第k-1历元的解算，直至当前时间段数据解算完毕，进入下一时间段的反向滤波解算。
本发明还提供了一种基于多个连续运行参考站的后处理定位系统，所述系统包括：
数据接收与存储单元，用于服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
流动站数据预处理单元，用于对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
参考站数据处理单元，用于对连续运行参考站数据进行处理，计算参考站差分改正数；
精准位置解算单元，用于基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算；
结果合并与存储单元，用于将解算结果按照时间先后合并位置结果并存储。
进一步地，所述精准位置解算步骤如下：
对每一个时间段的数据进行正向滤波解算，从第一个历元开始利用正向滤波的方法解算各个历元的位置和模糊度参数；
对每一个时间段的数据进行反向滤波解算，从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数；
正向滤波解算和反向滤波解算完成后，对各个时间段内同一历元的正向滤波结果和反向滤波结果进行联合处理。
本发明还提供了一种存储器，所述存储器存储有计算机程序，所述计算机程序执行如下步骤：
服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
对连续运行参考站数据进行处理，计算参考站差分改正数；
基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。
本发明的有益效果在于，利用多个连续运行参考站对流动站数据解算，扩展RTK定位范围的同时利用历史信息提高了定位的精度。
附图说明
图1为本发明数据处理方法流程图；
图2为本发明流动站数据预处理流程图；
图3为本发明精准位置解算流程图；
图4为本发明联合处理流程图；
图5为本发明数据处理系统结构图。
具体实施方式
本发明提出利用流动站和多个参考站数据处理方案，根据流动站的位置动态生成对应虚拟参考站的观测值，并对流动站数据进行对应分割，根据各时段的数据采用正反向滤波提高流动站定位精度。下文中，结合附图和实施例对本发明作进一步阐述。
实施例一：
本发明提供了一种基于多个连续运行参考站的后处理定位方法，数据处理如图1所示，具体步骤如下：
用户将接收到的数据上传到服务器，同时服务器接收GNSS连续运行参考站的数据进行存储。用户和连续运行参考站的数据是指GNSS的观测数据，它包含了每一历元的时间和每一颗卫星的载波相位与伪距的观测值。所不同的是用户数据是在动态或静态下采集，而连续运行参考站一般静态不动，且参考站的坐标已知。
流动站数据预处理的目的是为了得到流动站在每一历元的概略位置并根据概略位置将流动站原始数据进行分割，如图2所示。首先获得流动站数据对应的广播星历，计算每一历元流动站单点的位置。然后对流动站原始数据进行分割，分割的规则为：从用户数据的第一个时刻起向空的经纬度序列中加入流动站的坐标，并分别计算最大最小纬度之差对应的距离长度ΔB和最大最小经度之差对应的距离长度ΔD。其中
ΔB＝(Bmax-Bmin)·h (1)
式(1)中Bmax，Bmin分别为纬度的最大值和最小值，h＝πR/180表示地球表面每一度所对应的长度，R＝6378137.0m为地球半径。式(2)中
ΔL＝Lmax-Lmin，Lmax，Lmin分别为经度的最大值和最小值。当ΔB或ΔD大于10km(阈值)时将经纬度序列对应的时间段作为一个格网进行划分，并分别记录存储经度，纬度的平均值作为一个格网点的坐标等待后续使用。后续未划分的数据再次利用上述划分方法对数据进行划分，直至最后一段数据均符合分割的规则。分割流动站数据建立格网点的优势是每段数据中的每一时间点的位置与格网点的距离保持在较短距离，进而保证定位精度。
参考站数据处理主要是根据上述过程记录的各个格网点时间段与格网点坐标，选取对应的参考站并计算对应的参考站差分改正数。根据各个格网的时间段和格网中心的位置进行参考站的选取，选取的规则为：(1)以格网的位置为起始搜索距离在70km以内的所有参考站；(2)对各个参考站的数据进行数据质量评估，保留评估合格的参考站；(3)计算上一步中各参考站站对应于格网中心的方位角；(4)分别在0°至90°、90°至180°、180°至270°、270°至360°的四个区域内找到各区域内距离中心点最近的4个站，然后在剩余的所有站点中找到两个距离格网中心最近的2个站；(6)若某一区域缺少站，则就近选择1个站作为补充。6个站选定后根据该格网点记录的对应时间，利用选定的参考站接收到的GNSS数据和广播星历数据，计算该格网中心对应的参考站差分改正数。其余格网的数据采用相同的策略计算每一格网点对应时间段内的参考站差分改正数。
将流动站各时段的GNSS数据分割完毕，并得到对应的参考站差分改正数和广播星历等数据后，对每一段数据进行精准位置解算，如图3所示。精准位置解算具体步骤如下：
(1)对每一段数据进行正向滤波，正向滤波对当前时段的数据从第一个历元开始利用滤波的方法解算各个历元的位置和模糊度参数，解算过程中第k个历元时在获得浮点解后可进行模糊度固定及检核，若固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，然后继续进行第k+1历元的解算，直至当前时段数据解算完毕，进入下一时段的正向滤波解算。
(2)对每一段数据进行反向滤波，每一段数据从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数，解算过程中第k个历元时在获得浮点解后可进行模糊度固定及检核，若固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，与正向滤波不同的是第k个历元解算完成后进行第k-1历元的解算，直至当前时段数据解算完毕进入下一时段的反向滤波解算。正向滤波在流动站数据周跳，遮挡发生后需要进行初始化，而反向滤波按照时间倒序进行能够提高这些时段的解算精度。
(3)正反向滤波解算完成后，对各个数据段内同一历元的正反向滤波结果进行联合处理，如图4所示。若正反向滤波模糊度均已固定，直接采用式(3)和式(4)对位置参数进行处理。
Qc＝Qa-Qa(Qa+Qb)-1Qa (4)
其中x，Q分别表示参数的估值及其方差协方差矩阵，下标a表示正向滤波结果，b表示反向滤波结果，c表示正反向滤波联合处理结果。若正反向滤波存在模糊度不固定的情况，对位置和模糊度参数采用式(3)和式(4)的方法进行联合处理并尝试模糊度固定及检核，即得到第k历元的结果，然后进行第k+1历元的联合处理。联合处理利用了正向和反向的信息，融合后能够显著提高解算结果的精度。
每一段数据均处理完成后，对所有历元的数据结果按照时间先后合并位置结果并存储。
实施例二：
本发明提供了一种基于多个连续运行参考站的后处理定位系统，如图5所示，包括：
数据接收与存储单元，用于服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
流动站数据预处理单元，用于对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
参考站数据处理单元，用于对连续运行参考站数据进行处理，计算参考站差分改正数；
精准位置解算单元，用于基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算；
结果合并与存储单元，用于将解算结果按照时间先后合并位置结果并存储。
进一步地，所述精准位置解算步骤如下：
对每一个时间段的数据进行正向滤波解算，从第一个历元开始利用正向滤波的方法解算各个历元的位置和模糊度参数；
对每一个时间段的数据进行反向滤波解算，从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数；
正向滤波解算和反向滤波解算完成后，对各个时间段内同一历元的正向滤波结果和反向滤波结果进行联合处理。
实施例三：
本发明还提供了一种存储器，所述存储器存储有计算机程序，所述计算机程序执行如下步骤：
服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
对连续运行参考站数据进行处理，计算参考站差分改正数；
基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
Claims (12)
1.一种基于多个连续运行参考站的后处理定位方法，其特征在于，
所述方法包括以下步骤：
服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
对连续运行参考站数据进行处理，计算参考站差分改正数，其中，所述对连续运行参考站数据进行处理具体步骤如下：根据各个格网的时间段与格网点的坐标，选取对应的连续运行参考站接收GNSS观测数据和广播星历数据，计算对应的连续运行参考站差分改正数，其中，以格网的位置为起始，搜索符合要求的所有连续运行参考站；对各个连续运行参考站的数据进行数据质量评估，保留评估合格的连续运行参考站；计算各个连续运行参考站对应于格网中心的方位角；分别在0°至90°、90°至180°、180°至270°、270°至360°的四个区域内找到各区域内距离中心点最近的4个连续运行参考站，然后在剩余的所有连续运行参考站中找到2个距离格网中心最近的连续运行参考站；如果一个区域缺少连续运行参考站，则就近选择1个连续运行参考站作为补充；
基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。
2.如权利要求1所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
所述用户流动站数据和连续运行参考站数据是指GNSS的观测数据，包含每一历元的时间和每一颗卫星的载波相位与伪距的观测值。
3.如权利要求2所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
所述用户流动站数据是在动态或静态下采集，所述连续运行参考站数据是在静态下采集。
4.如权利要求2所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
所述对流动站原始数据进行分割具体包括以下步骤：
从用户流动站数据的第一个时刻起向空的经纬度序列中加入流动站的坐标；
分别计算最大纬度最小纬度之差对应的距离长度ΔB和最大经度最小经度之差对应的距离长度ΔD；
根据ΔB或者ΔD决定是否对流动站原始数据进行分割。
5.如权利要求4所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
如果ΔB或者ΔD大于阈值，则对流动站原始数据进行分割，将经纬度序列对应的时间段作为一个格网进行划分，并分别记录存储经度，纬度的平均值作为一个格网点的坐标等待后续使用。
6.如权利要求5所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
对于未划分的数据再次进行划分，直至最后一个时间段数据均符合流动站原始数据分割的规则。
7.如权利要求6所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
根据数据集对每一个时间段的数据进行精准位置解算，所述数据集包括分割完成后的用户流动站数据、连续运行参考站差分改正数和广播星历数据，所述精准位置解算步骤如下：
对每一个时间段的数据进行正向滤波解算，从第一个历元开始利用正向滤波的方法解算各个历元的位置和模糊度参数；
对每一个时间段的数据进行反向滤波解算，从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数；
正向滤波解算和反向滤波解算完成后，对各个时间段内同一历元的正向滤波结果和反向滤波结果进行联合处理。
8.如权利要求7所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
所述正向滤波解算过程中第k个历元时在获得浮点解后进行模糊度固定及检核，如果固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，然后继续进行第k+1历元的解算，直至当前时间段数据解算完毕，进入下一时间段的正向滤波解算。
9.如权利要求7所述的一种基于多个连续运行参考站的后处理定位方法，其特征在于，
反向滤波解算过程中第k个历元时在获得浮点解后进行模糊度固定及检核，如果固定成功则存储固定后的模糊度和位置参数，否则存储模糊度浮点解和位置参数，然后继续进行第k-1历元的解算，直至当前时间段数据解算完毕，进入下一时间段的反向滤波解算。
10.一种基于多个连续运行参考站的后处理定位系统，其特征在于，
所述系统包括：
数据接收与存储单元，用于服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
流动站数据预处理单元，用于对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
参考站数据处理单元，用于对连续运行参考站数据进行处理，计算参考站差分改正数，其中，所述对连续运行参考站数据进行处理具体步骤如下：根据各个格网的时间段与格网点的坐标，选取对应的连续运行参考站接收GNSS观测数据和广播星历数据，计算对应的连续运行参考站差分改正数，其中，以格网的位置为起始，搜索符合要求的所有连续运行参考站；对各个连续运行参考站的数据进行数据质量评估，保留评估合格的连续运行参考站；计算各个连续运行参考站对应于格网中心的方位角；分别在0°至90°、90°至180°、180°至270°、270°至360°的四个区域内找到各区域内距离中心点最近的4个连续运行参考站，然后在剩余的所有连续运行参考站中找到2个距离格网中心最近的连续运行参考站；如果一个区域缺少连续运行参考站，则就近选择1个连续运行参考站作为补充；
精准位置解算单元，用于基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算；
结果合并与存储单元，用于将解算结果按照时间先后合并位置结果并存储。
11.如权利要求10所述的一种基于多个连续运行参考站的后处理定位系统，其特征在于，
所述精准位置解算步骤如下：
对每一个时间段的数据进行正向滤波解算，从第一个历元开始利用正向滤波的方法解算各个历元的位置和模糊度参数；
对每一个时间段的数据进行反向滤波解算，从最后一个历元开始利用反向滤波的方法解算各个历元的位置和模糊度参数；
正向滤波解算和反向滤波解算完成后，对各个时间段内同一历元的正向滤波结果和反向滤波结果进行联合处理。
12.一种存储器，所述存储器存储有计算机程序，其特征在于，
所述计算机程序执行如下步骤：
服务器接收上传的用户流动站数据和连续运行参考站数据并进行存储；
对用户流动站数据进行预处理，得到流动站在每一历元的概略位置，并根据概略位置对流动站原始数据进行分割；
对连续运行参考站数据进行处理，计算参考站差分改正数，其中，所述对连续运行参考站数据进行处理具体步骤如下：根据各个格网的时间段与格网点的坐标，选取对应的连续运行参考站接收GNSS观测数据和广播星历数据，计算对应的连续运行参考站差分改正数，其中，以格网的位置为起始，搜索符合要求的所有连续运行参考站；对各个连续运行参考站的数据进行数据质量评估，保留评估合格的连续运行参考站；计算各个连续运行参考站对应于格网中心的方位角；分别在0°至90°、90°至180°、180°至270°、270°至360°的四个区域内找到各区域内距离中心点最近的4个连续运行参考站，然后在剩余的所有连续运行参考站中找到2个距离格网中心最近的连续运行参考站；如果一个区域缺少连续运行参考站，则就近选择1个连续运行参考站作为补充；
基于流动站数据预处理结果和连续运行参考站数据处理结果对数据进行精准位置解算，将解算结果按照时间先后合并位置结果并存储。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811303873.4A CN111142133B (zh) | 2018-11-02 | 2018-11-02 | 基于多个连续运行参考站的后处理定位方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811303873.4A CN111142133B (zh) | 2018-11-02 | 2018-11-02 | 基于多个连续运行参考站的后处理定位方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111142133A CN111142133A (zh) | 2020-05-12 |
| CN111142133B true CN111142133B (zh) | 2022-06-24 |
Family
ID=70515497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811303873.4A Active CN111142133B (zh) | 2018-11-02 | 2018-11-02 | 基于多个连续运行参考站的后处理定位方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111142133B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111596330B (zh) * | 2020-05-22 | 2022-11-11 | 深圳思凯微电子有限公司 | 定位方法、装置、终端和存储介质 |
| CN116106948A (zh) * | 2021-11-09 | 2023-05-12 | 千寻位置网络（浙江）有限公司 | 网络rtk抗电离层干扰定位方法及相关设备 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101661096A (zh) * | 2008-08-29 | 2010-03-03 | 夏晓清 | 基于三角形逼近算法的虚拟参考站的生成方法及系统 |
| CN101614808B (zh) * | 2009-07-30 | 2012-04-25 | 北京北斗星通导航技术股份有限公司 | 一种航迹测试系统 |
| CN102054388B (zh) * | 2009-11-03 | 2013-01-23 | 厦门雅迅网络股份有限公司 | 一种接收和显示地图的方法 |
| CN102118198B (zh) * | 2009-12-31 | 2016-01-20 | 中兴通讯股份有限公司 | 一种终端反馈宽带信道质量指示的系统及方法 |
| CN103065472A (zh) * | 2012-12-24 | 2013-04-24 | 中国科学院深圳先进技术研究院 | 实时路况分析方法及系统 |
| CN104007457A (zh) * | 2014-05-28 | 2014-08-27 | 武汉市测绘研究院 | 一种基于cors系统的无参数实时定位方法 |
| CN106019336B (zh) * | 2015-08-28 | 2019-07-23 | 千寻位置网络有限公司 | 全球导航卫星系统的差分中继方法及其装置 |
| CN106814379B (zh) * | 2017-01-19 | 2019-07-05 | 湖南北云科技有限公司 | 一种云rtk定位方法及系统 |
| CN107422351A (zh) * | 2017-08-02 | 2017-12-01 | 湖南省测绘科技研究所 | 一种基于虚拟网格的gnss分米级差分定位方法 |
-
2018
- 2018-11-02 CN CN201811303873.4A patent/CN111142133B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111142133A (zh) | 2020-05-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
| CN112230252B (zh) | 终端定位方法、装置、计算机设备和存储介质 | |
| Wen et al. | 3D LiDAR aided GNSS NLOS mitigation in urban canyons | |
| JP5794646B2 (ja) | 衛星測位システム、測位端末、測位方法、及びプログラム | |
| CN108828626B (zh) | 基于实时网格的网络rtk电离层延迟内插方法及系统 | |
| CN108845345B (zh) | 利用gnss测速原理的双天线定向测姿的方法 | |
| CN111142133B (zh) | 基于多个连续运行参考站的后处理定位方法和系统 | |
| CN114562992A (zh) | 一种基于因子图及场景约束的多径环境组合导航方法 | |
| CN116125514A (zh) | 基于北斗ppp-rtk虚拟观测值地灾监测方法、装置、终端及介质 | |
| CN115343745A (zh) | 三维激光雷达辅助的高精度卫星定位方法 | |
| CN115932923A (zh) | 一种基于v2v的共享gnss车辆增强协作定位方法 | |
| Bisnath et al. | Innovation: Examining precise point positioning now and in the future | |
| CN116381760B (zh) | Gnss rtk/ins紧耦合定位方法、装置及介质 | |
| CN110850447B (zh) | 对列车定位单元的定位精度进行综合评估的方法 | |
| CN111736183A (zh) | 一种联合bds2/bds3的精密单点定位方法和装置 | |
| CN112731493B (zh) | 一种精度高的静态伪距单点定位方法及装置 | |
| CN115932920A (zh) | 对流层延迟的插值方法 | |
| CN113534206B (zh) | 基于北斗地基增强系统的接入虚拟参考站快速选择方法 | |
| Marques et al. | Shoreline monitoring by gnss-ppp aiming to attendance the law 14.258/2010 from Pernambuco state, Brazil | |
| CN112578423A (zh) | 整周模糊度确定方法、装置及设备 | |
| CN110208841B (zh) | 一种改进的面向非重叠频率的gnss紧组合方法 | |
| CN116125506A (zh) | 一种多基准站的载波相位的差分定位方法及装置 | |
| CN113933872A (zh) | 多系统差分定位方法及其系统 | |
| US20240118435A1 (en) | Post-processing of global navigation satellite system (gnss) data | |
| CN117289318B (zh) | 一种基于bds辅助的glonass频间偏差实时估计方法 |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 200438 9 / F, 10 / F, 11 / F, 12 / F, 38 Lane 1688, Guoquan North Road, Yangpu District, Shanghai Applicant after: QIANXUN SPATIAL INTELLIGENCE Inc. Address before: Room j165, 1st floor, building 64, 1436 Jungong Road, Yangpu District, Shanghai, 200433 Applicant before: QIANXUN SPATIAL INTELLIGENCE Inc. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |