CN110927761B - 一种gnss阵列式联合高精度定位算法 - Google Patents
一种gnss阵列式联合高精度定位算法 Download PDFInfo
- Publication number
- CN110927761B CN110927761B CN201911171580.XA CN201911171580A CN110927761B CN 110927761 B CN110927761 B CN 110927761B CN 201911171580 A CN201911171580 A CN 201911171580A CN 110927761 B CN110927761 B CN 110927761B
- Authority
- CN
- China
- Prior art keywords
- observation
- equation
- receiver
- representing
- point
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000006073 displacement reaction Methods 0.000 claims description 21
- 238000003339 best practice Methods 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000005433 ionosphere Substances 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 239000005436 troposphere Substances 0.000 claims description 3
- 230000003068 static effect Effects 0.000 abstract description 4
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002344 surface layer Substances 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/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
本发明公开了一种GNSS阵列式联合高精度定位算法,是基于GNSS定位原理,参考局域差分技术,提出的一种适合GNSS定位系统的静态多点位GNSS阵列式联合定位算法,本算法首先使用载波相位观测法,得出单点位位置信息;其次通过点间双差方法,得出某一点位的多个位置信息,然后根据统计学原理,得出多个位置信息的有效值,并使用所有点位的有效值进行数学模型的建立,构建一个平面基准网络。本算法有效的解决了多点位联合解算的难题,并有效降低了数据解算难度,于此同时,本算法还很好的解决了实际使用中无法设立基准站或基准站基线过长无法精确解算的问题。
Description
技术领域
本发明涉及计算机数据处理算法领域,尤其涉及一种GNSS阵列式联合高精度定位算法。
背景技术
中国北斗卫星导航系统(BDS)是中国自行研制的全球卫星导航系统。是继美国全球定位系统(GPS)、俄罗斯格洛纳斯卫星导航系统(GLONASS)、欧洲伽利略卫星导航系统(GSNS)之后第四个成熟的卫星导航系统。
北斗高精度定位是指利用BDS系统“静态基线相对定位技术”获取监测点的高精度坐标(毫米级),并利用数学运算,得出观测位置的位移变化量,其在地质灾害、桥梁、矿山等有着广泛应用价值。
目前采用的高精度定位技术,一般均采用“静态相对定位”,此种监测方式需要建立基准观测站,并人为假设其为静止状态,通过此基准站的位置信息,通过差分运算,得出与观测站间的相关信息,但是大部分应用场景,无法建立基准站。此种方法就给实际应用带来了困难。
发明内容
本发明所要解决的技术问题是克服现有技术中存在的不足,提供一种GNSS阵列式联合高精度定位算法。
本发明是通过以下技术方案予以实现:一种GNSS阵列式联合高精度定位算法,包括以下步骤:步骤(1),使用载波相位观测得出的单点位状态量进行多点位相对差分运算;步骤(2),根据统计学原理得出最佳观测值;步骤(3),基于多点位的最佳观测值构建平面基准网络。
进一步的,:所述步骤(1)中的载波相位观测方程为:
式中,表示接收机与卫星间的实际空间距离,dtu(t)和dti(t)分别表示接收机和卫星的钟差,/>为电离层延迟,/>为对流程延迟,ε表示测量噪声,c表示光速,t表示观测历元;
同理,可得出在t+1时刻的载波相位方程如下:
将式(1-1)与式(1-2)做差,得出此观测站对单一卫星的单差载波相位方程,如下:
详细方程为如下:
由于在相邻历元观测同一颗卫星,所以钟差可以消除;同时,相邻历元的时间范围内,可以认定对流层与电离层的变化微弱,所以此三个误差因素可以通过该单差方式消除;此外,由于观测时间短,且为同一颗卫星,可以认定观测噪声相同,同样可以消除;式(1-4)可简化为:
进一步的,将任意两个观测站的单差载波相位观测方程做差,得出基于两个相邻历元t和(t+1)任意两个观测站对同一颗卫星的双差方程为:
进一步的,所述任意一点u,通过方程(1-6)解算,得出(m-1)个三维坐标矩阵:
进一步的,将式(1-7)使用统计学方法,得出其最佳实践值:
进一步的,将式(1-8)所有点位进行平面拟合,得出平面拟合曲线,算法方程如下:
式中,S表示得出的拟合曲面方程,n表示所要拟合的所有点位数,带下标的X/Y/Z分别为式(1-8)得出的最佳实践值,带下标的a为设定的参数,其取得遵循如下方程:AX+BY+CZ+D=0,(C≠0)(1-10)
如上可得:
Z=-AX/C-BY/C-D/C (1-11)
我们假设参数:
可得Z=aoX+a1Y+a2(1-12)
为得出最佳拟合曲线,需满足如下约束条件:
由数学运算可得:
解上述方程组可得ao,a1,a2,即可得出式(1-12)的三维曲面拟合。
进一步的,将式(1-8)得出的点位最佳实践值与三维曲面拟合值进行基准库定义,得出所有观测点位的原始基准值并将原始基准值定义为基准曲面,以此为参考进行后续测量。
进一步的,将式(1-8)得出的点位最佳实践值与基准库进行逐一求差,得出点位三维位移量,方程如下:
ΔUn(t)=Un(t)-Un (1-15)
其中ΔUn(t)表示第n个点位的位移变化量,Un(t)为t时刻观测值,Un为基准值。
进一步的,将点位三维位移量带入到曲面拟合方程,并进行重拟合,得出新的拟合曲面,将基准曲面与拟合曲面通过数学运算,得出点位矩阵中非观测点位以外所有观测范围内的变化量。
本发明的有益效果是:本发明采用短基线下载波相位测量方程进行数据解算,并根据历元间平差及站间平差的双差方式,计算出多个观测值,在使用统计学原理得出每个观测点的基准观测值,再利用最小二乘法,得出曲面拟合方程,作为标准观测曲面。通过反复运算,得出观测时刻的位移量。可以解决短基线情况下,无基准站的问题,以提高短基线下的定位精度。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
为了解决上述问题,本发明采用的技术方案是:提供一种GNSS阵列式联合高精度定位算法,采用短基线下载波相位测量方程进行数据解算,并根据历元间平差及站间平差的双差方式,计算出多个观测值,在使用统计学原理得出每个观测点的基准观测值,再利用最小二乘法,得出曲面拟合方程,作为标准观测曲面。通过反复运算,得出观测时刻的位移量,
其中载波相位观测方程为:
式中,表示接收机与卫星间的实际空间距离,dtu(t)和dti(t)分别表示接收机和卫星的钟差,/>为电离层延迟,/>为对流程延迟,ε表示测量噪声,c表示光速,t表示观测历元。
同理,可得出在t+1时刻的载波相位方程如下:
将式(1-1)与式(1-2)做差,得出此观测站对单一卫星的单差载波相位方程,如下:
详细方程为如下:
由于在相邻历元观测同一颗卫星,所以钟差可以消除。同时,相邻历元的时间范围内,可以认定对流层与电离层的变化微弱,所以此三个误差因素可以通过该单差方式消除。此外,由于观测时间短,且为同一颗卫星,可以认定观测噪声相同,同样可以消除。式(1-4)可简化为:
将所述任意两个观测站的单差载波相位观测方程做差,得出基于两个相邻历元t和(t+1)任意两个观测站对同一颗卫星的双差方程为:
周跳探测与模糊度求解分别使用伪距与相位组合法及最小二乘模糊度降相关法。
通过上述解算,任意一点u,得出(m-1)个三维坐标矩阵:
式(1-7)使用统计学方法,得出其最佳实践值:
将式(1-8)所有点位进行平面拟合,得出平面拟合曲线,算法方程如下:
式中,S表示得出的拟合曲面方程,n表示所要拟合的所有点位数,带下标的X/Y/Z分别为式(1-8)得出的最佳实践值,带下标的a为设定的参数,其取得遵循如下方程:
AX+BY+CZ+D=0,(C≠0) (1-10)
如上可得:
Z=-AX/C-BY/C-D/C, (1-11)
进行人为参数设定:
可得Z=aoX+a1Y+a2 (1-12)
为得出最佳拟合曲线,需满足如下约束条件:
由数学运算可得:
对上述方程组求解可得ao,a1,a2,即可得出式(1-12)的三维曲面拟合。
将上述得出的点位最佳实践值与三维曲面拟合值进行基准库定义,得出所有观测点位的原始基准值,以此为参考进行后续测量。
使用相同算法得出任意观测时刻的点位最佳实践值,并与基准库进行逐一求差,得出点位三维位移量,方程如下:
ΔUn(t)=Un(t)-Un (1-15)
其中ΔUn(t)表示第n个点位的位移变化量,Un(t)为t时刻观测值,Un为基准值。
将点位三维位移量带入到曲面拟合方程,并进行重拟合,得出新的拟合曲面,将基准曲面与拟合曲面通过数学运算,得出点位矩阵中非观测点位以外所有观测范围内的变化量。
以一个滑坡监测为例:
当一个山体坡面表层出现松动时,或发生降水时,山体坡面标称会形成滑坡体,滑坡体会沿坡面径向方向位移,此位移可分解为水平位移和高程位移。位移变化量由小变大,最终导致滑坡产生。我们可以通过对滑坡体的水平位移分量和高程位移分量的变化进行滑坡体的监测。
设在一个滑坡体上部署了n个监测点,
根据上述算法,得出的观测点n的水平监测量与高程监测量为:
水平:Hn(0),垂直:Vn(0)(mm),
得出第n点的观测值矩阵:Wn(0)=[Hn(0),Vn(0)]
将该矩阵作为基准矩阵。
观测时刻t的观测值矩阵:Wn(t)=[Hn(t),Vn(t)]
将t时刻的观测值与初始基准阵列做差运算,即可得出该点位的位移变化量:
ΔWn(t)=Wn(t)-Wn(0),
即:ΔWn(t)=[ΔHn(t),ΔVn(t)]
将Hn(0),Vn(0)带入式(1-9)、(1-14)即可得出标准拟合曲面。
综上所述,本发明用短基线下载波相位测量方程进行数据解算,并根据历元间平差及站间平差的双差方式,计算出多个观测值,在使用统计学原理得出每个观测点的基准观测值,再利用最小二乘法,得出曲面拟合方程,作为标准观测曲面,通过反复运算,得出观测时刻的位移量,并能测算出监测范围内非观测点的可能位移变化量,从而达到短基线情况下,无基准站的问题,以此来提高短基线下的定位精度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种GNSS阵列式联合高精度定位算法,其特征在于:包括以下步骤:
步骤(1),使用载波相位观测得出的单点位状态量进行多点位相对差分运算;步骤(2),根据统计学原理得出最佳观测值;步骤(3),基于多点位的最佳观测值构建平面基准网络;
步骤(1)中的载波相位观测方程为:
式中,表示接收机与卫星间的实际空间距离,/>和/>分别表示接收机和
卫星的钟差,为电离层延迟,/>为流程延迟,/>表示测量噪声,/>表示光速,/>表示观测历元,/>表示由地面站观测到的第/>颗卫星,/>表示地面标记为/>的接收机,/>表示在/>时刻接收机/>接收到的第/>颗卫星的导航信号的观测值,/>表示接收机/>接收到的第/>颗卫星的导航信号的观测值的整周模糊度;
同理,可得出在时刻的载波相位方程如下:
将式(1-1)与式(1-2)做差,得出此观测站对单一卫星的单差载波相位方程,如下:
详细方程为如下:
由于在相邻历元观测同一颗卫星,所以钟差可以消除;同时,相邻历元的时间范围内,可以认定对流层与电离层的变化微弱,所以此三个误差因素可以通过该单差方式消除;此外,由于观测时间短,且为同一颗卫星,可以认定观测噪声相同,同样可以消除;式(1-4)可简化为:
2.根据权利要求1所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将任意两个观测站的单差载波相位观测方程做差,得出基于两个相邻历元和(/>任意两个观测站对同一颗卫星的双差方程为:
式中,表示地面标记为/>的接收机,/>表示与接收机关联的监测站,/>表示由地面站观测到的第/>颗卫星,/>表示相邻两个历元之间间隔/>时间内对应的同一卫星的接收机与监测站的观测值,/>表示相邻两个历元之间间隔/>时间内对应的同一卫星的接收机与监测站的实际空间距离,/>表示接收机与监测站对应的统一卫星的整周模糊度。
3.根据权利要求2所述的一种GNSS阵列式联合高精度定位算法,其特征在于:对于任意一点,通过方程(1-6)解算,得出/>个三维坐标矩阵:
式中,表示与接收机关联的监测站,/>表示包括接收机在内共有的点位数量,/>表示接收机/>与其余除/>点外的/>个点位做解算得出的个矩阵。
4.根据权利要求3所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将式(1-7)使用统计学方法,得出其最佳实践值:
式中,表示与接收机关联的监测站。
5.根据权利要求4所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将式(1-8)所有点位进行平面拟合,得出平面拟合曲线,算法方程如下:
式中,表示得出的拟合曲面方程,/>表示所要拟合的所有点位数,带下标的/>分别为式(1-8)得出的最佳实践值,带下标的/>为设定的参数,其取得遵循如下方程:
如上可得:
我们假设参数:
可得
为得出最佳拟合曲线,需满足如下约束条件:
其中/>
由数学运算可得:
解上述方程组可得,即可得出式(1-12)的三维曲面拟合。
6.根据权利要求5所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将式(1-8)得出的点位最佳实践值与三维曲面拟合值进行基准库定义,得出所有观测点位的原始基准值并将原始基准值定义为基准曲面,以此为参考进行后续测量。
7.根据权利要求6所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将式(1-8)得出的点位最佳实践值与基准库进行逐一求差,得出点位三维位移量,方程如下:
其中,表示第/>个点位的位移变化量,/>为/>时刻观测值,/>为基准值。
8.根据权利要求7所述的一种GNSS阵列式联合高精度定位算法,其特征在于:将点位三维位移量带入到曲面拟合方程,并进行重拟合,得出新的拟合曲面,将基准曲面与拟合曲面通过数学运算,得出点位矩阵中非观测点位以外所有观测范围内的变化量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911171580.XA CN110927761B (zh) | 2019-11-22 | 2019-11-22 | 一种gnss阵列式联合高精度定位算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911171580.XA CN110927761B (zh) | 2019-11-22 | 2019-11-22 | 一种gnss阵列式联合高精度定位算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110927761A CN110927761A (zh) | 2020-03-27 |
CN110927761B true CN110927761B (zh) | 2023-09-05 |
Family
ID=69851900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911171580.XA Active CN110927761B (zh) | 2019-11-22 | 2019-11-22 | 一种gnss阵列式联合高精度定位算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110927761B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298151A (zh) * | 2011-07-20 | 2011-12-28 | 东南大学 | 一种gnss网络差分定位系统中的误差改正方法 |
CN104898144A (zh) * | 2015-06-09 | 2015-09-09 | 辽宁工程技术大学 | 一种北斗系统参考站整周模糊度单历元确定方法 |
CN105549048A (zh) * | 2016-01-29 | 2016-05-04 | 申研 | 基于三角形接收机阵列的高精度卫星定位装置及其方法 |
CN108061911A (zh) * | 2018-01-16 | 2018-05-22 | 东南大学 | 一种glonass载波单差残差估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105589087B (zh) * | 2016-01-29 | 2017-06-23 | 申研 | 一种高精度实时卫星定位装置及其方法 |
-
2019
- 2019-11-22 CN CN201911171580.XA patent/CN110927761B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102298151A (zh) * | 2011-07-20 | 2011-12-28 | 东南大学 | 一种gnss网络差分定位系统中的误差改正方法 |
CN104898144A (zh) * | 2015-06-09 | 2015-09-09 | 辽宁工程技术大学 | 一种北斗系统参考站整周模糊度单历元确定方法 |
CN105549048A (zh) * | 2016-01-29 | 2016-05-04 | 申研 | 基于三角形接收机阵列的高精度卫星定位装置及其方法 |
CN108061911A (zh) * | 2018-01-16 | 2018-05-22 | 东南大学 | 一种glonass载波单差残差估计方法 |
Non-Patent Citations (1)
Title |
---|
张帆 等.无人机集群作战应用卫星导航差分定位技术研究.现代导航.2019,(第3期),157-162. * |
Also Published As
Publication number | Publication date |
---|---|
CN110927761A (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2017128871A1 (zh) | 一种高精度实时卫星定位装置及其方法 | |
CN104483691B (zh) | 一种gnss组合精密单点定位方法 | |
CN116540303B (zh) | 一种历元间差分观测方程建立方法和地震位移计算方法 | |
Maciuk | GPS-only, GLONASS-only and combined GPS+ GLONASS absolute positioning under different sky view conditions | |
CN111044972B (zh) | 基于gnss精密时间同步的飞行器时差定位方法及系统 | |
CN103543454B (zh) | 一种嵌入在移动通讯网中的卫星定轨系统 | |
CN113358017B (zh) | 一种多站协同处理的gnss高精度变形监测方法 | |
CN109613582B (zh) | 一种车载实时单频米级伪距定位方法 | |
CN108732601A (zh) | 基于差分卫星导航的垂直返回运载器着陆段导航方法 | |
CN116009042A (zh) | 一种单站载波历元间差分实时探测相对形变的方法及系统 | |
CN110927761B (zh) | 一种gnss阵列式联合高精度定位算法 | |
CN103760582A (zh) | 一种遮挡环境下卫星双差观测结构的优化方法 | |
CN116858290B (zh) | 基于大型无人机的深远海海面高度观测定标方法和系统 | |
CN109977499B (zh) | 基于位置约束的北斗三频/静力水准仪斜拉桥监测方法 | |
CN117492035A (zh) | 一种用于解算ppp-rtk服务完好性信息的方法 | |
Karadeniz et al. | Comparison of Network-RTK and PPP Technique in terms of Position Accuracy | |
CN116047563A (zh) | 一种基于海洋北斗定位增强的水面高精度定位及评估方法 | |
CN115561793A (zh) | 一种基于并行计算的实时北斗相位小数偏差快速估计方法 | |
CN110082800B (zh) | 一种差分定位方法 | |
CN113126133A (zh) | 基于bds或gps的中长基线多频rtk定位的快速收敛方法 | |
CN114137588A (zh) | 一种基于伪距观测噪声的伪距差分定位分区定权方法 | |
CN113899957B (zh) | 基于双差相位观测值的室外天线绝对标定方法及系统 | |
CN115144882B (zh) | 一种基于gnss和5g组合的精密单点定位方法及系统 | |
Zhang et al. | Research on virtual reference station (VRS) technology and differential correction information | |
CN110646817A (zh) | 计算定位误差的方法和高精度定位方法 |
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 |