CN106125110B - 基于分区改正的星基增强***定位精度提高方法 - Google Patents
基于分区改正的星基增强***定位精度提高方法 Download PDFInfo
- Publication number
- CN106125110B CN106125110B CN201610407689.9A CN201610407689A CN106125110B CN 106125110 B CN106125110 B CN 106125110B CN 201610407689 A CN201610407689 A CN 201610407689A CN 106125110 B CN106125110 B CN 106125110B
- Authority
- CN
- China
- Prior art keywords
- observation
- data
- station
- correction
- residual
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 122
- 230000003416 augmentation Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000006872 improvement Effects 0.000 title claims abstract description 7
- 230000004807 localization Effects 0.000 title abstract 3
- 230000009467 reduction Effects 0.000 claims abstract description 7
- 238000005192 partition Methods 0.000 claims description 60
- 230000001934 delay Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 9
- 238000013461 design Methods 0.000 abstract description 7
- 238000013316 zoning Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 12
- 239000005433 ionosphere Substances 0.000 description 9
- 239000005436 troposphere Substances 0.000 description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000011541 reaction mixture Substances 0.000 description 2
- 241000284466 Antarctothoa delta Species 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 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/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/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
-
- 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/10—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
- G01S19/11—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals wherein the cooperating elements are pseudolites or satellite radio beacon positioning system signal repeaters
-
- 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
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
本发明提供一种基于分区改正的星基增强***定位精度提高方法,包括,将一待观测区域划分为复数个观测区,并在每一观测区内设置复数个观测站;通过处理观测数据获取每一观测站的残差数据;根据观测残差数据获取各观测站的钟差均值数据;根据残差数据获取各观测区的模糊度归算值数据;根据残差数据和模糊度归算值数据获取各观测区的分区改正数数据;通过网络或者卫星链路向用户提供分区改正数数据调用服务。本发明的基于分区改正的星基增强***定位精度提高方法,通过对观测站进行广域的分区设计,并且基于观测站计算用于网络和星基播发的分区改正数,精化了用户定位中的误差改正,从而能够提高用户定位精度,具有精度高、适用范围广的优点。
Description
技术领域
本发明涉及广域星基增强领域,尤其涉及一种基于分区改正的星基增强***定位精度提高方法。
背景技术
全球导航卫星定位***(Global Navigation Satellite System,GNSS)提供的服务包括基本导航以及星基增强。其中星基增强为了提高***实时服务的精度,基本思想是对卫星轨道误差、卫星钟差及电离层延迟等主要误差源进行区分,建立每一种误差源的模型,用于修正基本导航中这些参数的误差。计算得到的这些改正数称为广域差分改正数,卫星电文上注***则将广域差分改正数通过卫星链路广播给用户。
目前的导航增强***主要由美国的WAAS(Wide Area Augmentation System,广域增强***)***,欧洲的EGNOS(European Geostationary Navigation Overlay Service,欧洲地球同步导航覆盖服务)***,日本的MSAS(Multi-functional SatelliteAugmentation System,多功能卫星增强***)***以及俄罗斯的SDCM(DifferentialCorrections and Monitoring,差分修正与监测)***等。这些***独立于GPS或GLONASS(GLOBAL NAVIGATION SATELLITE SYSTEM,全球卫星导航***)的运行控制***,基于地面测站的观测数据计算各种差分改正数,用于改正导航电文中的误差。
对导航***卫星的观测值包括了伪距和相位,理论上,相位观测值的精度比伪距要高100倍左右。相位观测数据包含了卫星和地面测站的测距信息以及未知的整周模糊度。实时处理模式下,模糊度的连续处理存在较长的收敛时间,在数据中断或者周跳的情况下,模糊度需要重新收敛。考虑到以上相位数据处理的复杂性,GPS/GLONASS星基增强***通常采用伪距观测值计算以上轨道、钟差以及格网电离层等广域差分改正数。我国区域卫星导航***采用了卫星等效钟差的算法。等效钟差基本思想如下所述:首先对若干测站位置已知的接收机伪距观测数据进行各项公共误差改正,得到新的距离观测量;在此基础上,计算距离观测量和测站卫星理论距离的差值组成残差观测量;将此残差观测量认为完全是测站、卫星钟差的误差,求解卫星钟差改正数。由于以上计算过程忽略了卫星轨道误差的影响,因此被称为等效钟差改正数。
星基增强技术在卫星导航定位***广播电文的基础上,计算更为精确的差分改正信息,并且通过卫星链路向用户进行广播。我国区域卫星导航***采用了基本导航以及星基增强一体化的设计,其中星基增强***在北斗服务区内面向授权用户提供更高精度的服务。
北斗***播发的差分信息包括等效钟差以及电离层格网改正数。其中等效钟差主要用于改正卫星钟差的快变误差,电离层格网改正数则是为了进一步提高单频接收机用户的电离层改正的精度。北斗授权用户采用***播发的广域差分改正数后的定位精度约为5米,相对于基本导航10米的服务精度有明显改善。
但现有星基增强***的处理技术存在以下问题:(1)只采用伪距观测量,广域差分改正数受伪距精度的影响,对定位精度的提高存在限制;(2)若采用相位观测量,由于相位观测值中存在模糊度以及周跳,导致相应的数据处理存在收敛时间较长的问题;(3)空间传播段的改正(主要是电离层)也是基于伪距观测值,对定位精度的提高存在限制。
由于目前星基增强处理只采用伪距观测数据,因此其服务精度仅限于米级,其定位精度最高只能达到1m左右,无法满足广域用户更高精度的定位需求。
发明内容
针对上述现有技术中的不足,本发明提供一种基于分区改正的星基增强***定位精度提高方法,通过对观测站进行广域的分区设计,并且基于观测站计算用于网络和星基播发的分区改正数,精化了用户定位中的误差改正,从而能够提高用户定位精度,具有精度高、适用范围广的优点。
为了实现上述目的,本发明提供一种基于分区改正的星基增强***定位精度提高方法,包括步骤:
S1:将一待观测区域划分为复数个观测区,并在每一所述观测区内设置复数个观测站;
S2:所述观测站观测至少一观测目标并获得观测数据;
S3:根据所述观测数据计算获取每一观测站的残差数据,所述残差数据包括伪距观测残差数据和相位观测残差数据;
S4:根据所述伪距观测残差数据计算获取各所述观测站的钟差均值数据;
S5:根据所述相位观测残差数据计算获取各所述观测区的模糊度归算值数据;
S6:根据所述残差数据和所述模糊度归算值数据计算获取各所述观测区的分区改正数数据,所述分区改正数数据包括伪距分区改正数数据和相位分区改正数数据;
S7:通过网络或者卫星链路向用户站提供分区改正数数据调用服务,所述用户站通过调用所述分区改正数数据获取所述用户站的三维坐标。
进一步地,所述步骤S3中,通过所述观测数据并利用公式(1)计算一第i观测站对一第j观测目标的伪距观测残差和所述第i观测站对所述第j观测目标的相位观测残差获取所述残差数据:
其中,i和j为大于零的自然数;f为频率,为所述第i观测站与所述第j观测目标间的距离观测误差,c为光速,dti为第i观测站的观测站钟差,dtj为第j观测目标的目标钟差,Δdtj为第j观测目标的目标钟差误差影响造成的观测误差;为模糊度参数,为与频率有关的电离层模型延迟改正误差,为所述第i观测站与所述第j观测目标间的对流层路径延迟,为所述第i观测站与所述第j观测目标间的对流层模型延迟改正误差,为所述第i观测站与所述第j观测目标间的第一剩余误差,为所述第i观测站与所述第j观测目标间的第二剩余误差。
进一步地,在所述步骤S4中,根据所述伪距观测残差数据并利用公式(2)获取所述第i观测站的钟差均值
其中,n为所述第i观测站观测到的所述观测目标的总数。
进一步地,在所述步骤S5中,根据所述相位观测残差数据并利用公式(3)获取各所述观测区的模糊度归算值数据
其中,为扣除所述观测站钟差后的相位残差改正数,m为当前观测区内的所述观测站个数,为所述第i观测站的钟差均值,t为历元编号。
进一步地,在所述步骤S6中,根据所述残差数据和所述模糊度归算值数据并利用公式(4)获取各所述伪距分区改正数数据ΔPj(f)和相位分区改正数数据ΔLj(f)|t:
其中,t为历元编号,为所述观测站的钟差均值。
本发明由于采用了以上技术方案,使其具有以下有益效果:
本发明中通过特定算法获取的分区改正数数据在发送给转发服务器后可向用户提供调用服务,用户根据需要通过对分区改正数数据的调用和使用,精化了用户定位中的误差改正,从而提高用户实现星基增强***定位精度的提高。
附图说明
图1为本发明实施例的一种基于分区改正的星基增强***定位精度提高方法的流程图;
图2为本发明实施例的复数个观测区的划分示意图。
具体实施方式
下面根据附图1-2,给出本发明的较佳实施例,并予以详细描述,使能更好地理解本发明的功能、特点。
请参阅图1,本发明的一种基于分区改正的星基增强***定位精度提高方法,包括步骤:
S1:将一待观测区域划分为复数个观测区,并在每一观测区内设置复数个观测站;
S2:观测站观测至少一观测目标获得观测数据;
S3:通过处理观测数据获取每一观测站的残差数据,残差数据包括伪距观测残差数据和相位观测残差数据;
具体地,在步骤S3中,通过观测数据并利用公式(1)计算一第i观测站对一第j观测目标的伪距观测残差和第i观测站对第j观测目标的相位观测残差获取残差数据:
其中,i和j为大于零的自然数;f为频率,为第i观测站与第j观测目标间的距离观测误差,c为光速,dti为第i观测站的观测站钟差,dtj为第j观测目标的目标钟差,为模糊度参数,为与频率有关的电离层模型延迟改正误差,为第i观测站与第j观测目标间的对流层路径延迟,为第i观测站与第j观测目标间的对流层模型延迟改正误差,为第i观测站与第j观测目标间的第一剩余误差,为第i观测站与第j观测目标间的第二剩余误差。
S4:根据伪距观测残差数据获取各观测站的钟差均值数据;
具体地,在步骤S4中,根据伪距观测残差数据并利用公式(2)获取第i观测站的钟差均值
其中,n为第i观测站观测到的观测目标的总数。
S5:根据相位观测残差数据获取各观测区的模糊度归算值数据;
具体地,在步骤S5中,根据相位观测残差数据并利用公式(3)获取各观测区的模糊度归算值数据
其中,为扣除观测站钟差后的相位残差改正数,m为当前观测区内的观测站个数,为第i观测站的钟差均值,t为历元编号。
S6:根据残差数据和模糊度归算值数据获取各观测区的分区改正数数据,分区改正数数据包括伪距分区改正数数据和相位分区改正数数据;
具体地,在步骤S6中,根据残差数据和模糊度归算值数据并利用公式(4)获取各伪距分区改正数数据ΔPj(f)和相位分区改正数数据ΔLj(f)|t:
其中,t为历元编号,为观测站的钟差均值。
S7:将分区改正数数据发送给一转发服务器并通过网络或者卫星链路向用户站提供调用服务,用户站通过对分区改正数数据的调用提高了其三维位置求取的精度。
下面对本发明的一种基于分区改正的星基增强***定位精度提高方法,做进一步描述。
本发明的内容包括星基增强相位分区改正数的计算模型以及适用于星基增强***的分区设计。
(1)首先,对于任意历元观测站i对观测目标(如:卫星)j在频点f的伪距、相位观测值为:
公式(5)中,为观测站i对观测目标j的伪距观测值,为观测站i对观测目标j的相位观测值;f为频点标识;为观测目标j与观测站i间的理论几何距离;c为光速;dti为第i观测站的观测站钟差,dtj为第j观测目标的目标钟差,为观测站i对观测目标j的相位观测值的模糊度,为与频率有关的电离层模型延迟改正误差,其对伪距、相位观测分别起到延迟和缩短观测时间的效果;为基于理论模型的对流层斜路径延迟;为第i观测站与第j观测目标间的第一剩余误差,为第i观测站与第j观测目标间的第二剩余误差,包含了多路径误差等噪声信息。
在公式(5)中,观测站已知坐标和广播星历求得的观测目标轨道、钟差,并利用实测的气象数据以及经验模型对大气延迟(电离层和对流层)进行改正,可求得第i观测站对第j观测目标的伪距观测残差和第i观测站对第j观测目标的相位观测残差
公式(1)中,Δdtj分别为观测站i受观测目标j轨道、钟差误差影响造成的观测误差。分别为电离层、对流层模型改正的剩余误差。观测残差中观测站钟差dti和模糊度参数是最大的项。为观测站上的残差改正,可作为局域/广域的增强信息提供给用户使用。
在相同历元,一用户站u对观测目标j在频点f的伪距、相位观测值为:
公式(6)中变量含义与公式(5)相同,只是测站由观测站i变成了用户站u。
其中,表示用户站对流层斜路径延迟;表示用户站与第j观测目标间的第一剩余误差;
表示用户站的电离层模型改正误差;
表示用户站u对观测目标j的相位观测值的模糊度;
dtu表示用户站钟差。
利用广播星历求得观测目标j轨道、钟差代入公式(6),并利用实测的气象数据以及经验模型对大气延迟(电离层和对流层)进行改正,得到:
公式(7)中,为基于广播星历计算的星地距离,为用户站u受卫星轨道误差影响造成的观测误差。对公式(7)引入该历元观测站i上的残差改正写为:
将公式(1)代入公式(8),得到:
其中为伪距剩余误差及观测噪声,为相位剩余误差及观测噪声。
公式(9)中,卫星钟差误差c·Δdtj得到了消除;用户站距离观测站如果小于2000公里,则的影响也在毫米级,可以忽略;观测站的站钟差c·dti可以被用户站站钟c·dtu完全吸收,成为新的站钟在观测站和用户站都没有发生周跳的情况下将重新组合为新的模糊度参数通过以上分析,重新组合部分项,公式(9)可重新写为:
与常规的观测方程相比,通过采用观测站提供的残差,公式(10)消除了观测目标钟差的影响;此外,公式(10)中对于离观测站千公里远的测站,电离层、对流层误差影响最大的公共部分得到了消除,从而降低了这两部分的影响。因此,基于观测站提供的残差改正,能够精化用户定位中的误差改正,提高用户定位精度。
(2)关于计算不同观测站钟差均值以及不同观测站公共卫星模糊度归算。
由于相位改正数中包含部分模糊度残余项,对同一颗卫星,不同历元的所能观测到的测站数量不一样,如果不进行公共卫星模糊度归算,这时相位综合改正数中所包含的模糊度残余项不一样,导致相位综合改正数不连续。
以上基于单个观测站能够提供相应的伪距相位改正信息。由于观测站残差中包含的观测站钟差能够被用户站钟差吸收,因此为了减少参数的位数,可在公式(1)残差中扣除钟差均值,形成新的残差改正数:
公式(11)中,分别为扣除观测站钟差后的伪距相位残差改正数;为观测站的钟差均值,可通过公式(1)中的伪距残差进行计算:
公式(2)中,n为观测站i观测到的观测目标总数。
单个观测站可能存在故障等问题,从而服务的连续性将受到影响。在一定的区域内可以建立多个观测站,可将单站模型拓展到多个观测站。参照公式(10),在任意历元,用户站u接收观测站k的改正数之后,对观测目标j在频点f的伪距、相位观测值为:
公式12中,
表示用户站新的站钟;
表示组合的新的模糊度参数;
表示基于广播星历计算的星地距离;
表示用户站电离层模型改正的剩余误差;
表示观测站k电离层模型改正的剩余误差;
表示用户站对流层模型改正的剩余误差;
表示观测站k对流层模型改正的剩余误差;
δ′表示伪距剩余误差及观测噪声;
ε′表示相位剩余误差机观测噪声。
比较公式(10)、公式(12),可见除了电离层、对流层改正信息有所差异之外,钟差、模糊度信息是组合了不同观测站相关项后的结果,存在差异。站钟由于是每个历元进行求解,因此不同观测站引起的钟差跳变会直接影响用户站钟差,而不会影响用户定位结果。模糊度信息在相邻历元的变化将导致用户的周跳,从而造成用户定位结果的跳变。因此在一个观测站存在故障需要切换至另外一个观测站时,为保证用户定位的平稳、连续,需要保证用户模糊度参数的连续,要求不同观测站对相同观测目标的模糊度参数保持一致。
根据公式(1),将两个观测站的残差作差得到:
公式(13)中下标i,k代表两个观测站间作差。公式(13)中要保证不同观测站对同一观测目标模糊度参数的一致性,则两个观测站对所有观测目标的模糊度必须相同或者相差一个相同的常数而实际上不同观测站对不同观测目标模糊度的差异不是常数,因此必须进行多观测站对相同观测目标模糊度的归算。
模糊度的归算是为了保证相位分区综合改正数的连续性,也即:需要保证其前后历元中包含的模糊度相同。采用对于某个观测目标,如果前后历元用于计算其综合改正数的观测站个数产生变化,则取前后历元公共的观测站数据,计算其改正数(扣除了站钟)均值的变化:
公式(14)中,为相邻历元综合改正数的变化,也即模糊度归算值;m为相同观测站的个数,下标t代表历元编号。观测站个数发生变化后,前面历元对该观测目标加上以上模糊度归算值就保证了模糊度参数的连续性。
(3)关于分区改正数数据的获取。
在进行了观测站钟差以及模糊度归算之后,可对相同分区内不同观测站的残差改正数进行综合,形成每个分区内每颗卫星的改正数,即为综合分区改正数。
多观测站综合分区改正数包括伪距改正数的综合和相位改正数的综合。即分区改正数数据包括伪距分区改正数数据和相位分区改正数数据。由于不同观测站扣除钟差均值后的伪距改正数的差异主要是指观测噪声的差异,因此可以直接对其进行综合平均。多站相位分区改正数综合时,通过模糊度归算,可直接进行叠加:
公式(15)中下标t,t-1代表历元编号,其中初始历元相位分区综合改正数为:
获得综合分区改正数后,通过网络或者卫星链路将其发送给用户使用。
(4)关于观测区的划分
将本发明提出的基于测站残差改正的分区改正模型可应用于北斗卫星导航***星基增强***。可依据现有观测站的分布以及***服务指标需求,对北斗服务区进行分区设计。
对于重点服务区域,如我国东南沿海地区,其分区大小可以适当密一些,应尽量保证每个分区内至少有两个观测站,每个观测站应分尽量分布均匀或在分区中心,另外后期也可以根据需要新增若干观测站。
例如,请参见图2,本实施例的一种观测区的分区设计,每个分区覆盖的范围各不相同,并且每个分区含有若干个观测站。其具体分区坐标,可见表1。
表1分区格网点坐标
以华东地区分区为例,假设在华东地区设置南京,上海,杭州三个观测站,在每个历元分别计算这三个观测站对观测到的每个观测目标(如卫星)的分区改正数。计算步骤如下:
1)每个观测站残差计算
以1号卫星为例,根据公式(1),在不同历元1、2、3、4时各观测站的伪距、相位残差计算结果请见表2。
表2单观测站伪距相位分区改正数
注:--表示没有数据或数据被剔除。
2)计算不同观测站钟差均值
根据公式(2),计算不同观测站的钟差均值并扣除,计算结果请参见表3。
表3扣除不同观测站站钟的钟差均值
3)多观测站公共卫星模糊度归算
根据公式(3),对于不同观测站的相位改正数,进行卫星模糊度归算,计算结果请参见表4。
历元 | 模糊度归算值 |
1 | -- |
2 | -0.010 |
3 | 0.0415 |
4 | 0.0155 |
表4多观测站卫星模糊度归算值
4)分区改正数综合
按照公式(13),对分区改正数进行综合,得到伪距分区改正数和相位分区改正数,计算结果请参见见表5。
历元 | 伪距分区改正数 | 相位分区改正数 |
1 | 0.9867 | 0.5043 |
2 | 1.0567 | 0.4943 |
3 | 1.0100 | 0.5358 |
4 | 1.1465 | 0.5513 |
表5分区综合改正数
本发明的通过对待观测区域进行广域的分区设计,并且基于观测站计算用于网络和星基播发的分区改正数。本发明的改正数数据精化了用户定位中的误差改正,从而能够提高用户定位精度。本发明的计算模型适用于任意频点及不同频点组合的分区改正数计算;也可以在不同频点计算之后进行组合获取不同频点组合的改正模型。同时,模型处理中也可以事先考虑其它已有等效钟差、格网电离层等各种广域差分改正数,进行叠加求取。
本发明通过计算同一分区内的误差改正数,修正了公共误差,精化了广域差分误差改正的精度,能够提高用户定位精度;本发明提出的模型可以扩展到基于相位观测的用户;本发明将待观测区域进行分区计算,各分区内的误差具有相似性,因此大大扩展了观测站与用户站的距离,实现了用户广域差分定位。
以上记载的,仅为本发明的较佳实施例,并非用以限定本发明的范围,本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰,皆落入本发明专利的权利要求保护范围。
Claims (5)
1.一种基于分区改正的星基增强***定位精度提高方法,包括步骤:
S1:将一待观测区域划分为复数个观测区,并在每一所述观测区内设置复数个观测站;
S2:所述观测站观测至少一观测目标并获得观测数据;
S3:根据所述观测数据计算获取每一观测站的残差数据,所述残差数据包括伪距观测残差数据和相位观测残差数据;
S4:根据所述伪距观测残差数据计算获取各所述观测站的钟差均值数据;
S5:根据所述相位观测残差数据计算获取各所述观测区的模糊度归算值数据;
S6:根据所述残差数据和所述模糊度归算值数据计算获取各所述观测区的分区改正数数据,所述分区改正数数据包括伪距分区改正数数据和相位分区改正数数据;
S7:通过网络或者卫星链路向用户站提供分区改正数数据调用服务,所述用户站通过调用所述分区改正数数据获取所述用户站的三维坐标。
2.根据权利要求1所述的基于分区改正的星基增强***定位精度提高方法,其特征在于,所述步骤S3中,通过所述观测数据并利用公式(1)计算一第i观测站对一第j观测目标的伪距观测残差和所述第i观测站对所述第j观测目标的相位观测残差获取所述残差数据:
其中,i和j为大于零的自然数;f为频率,为所述第i观测站与所述第j观测目标间的距离观测误差,c为光速,dti为第i观测站的观测站钟差,dtj为第j观测目标的目标钟差,Δdtj为第j观测目标的目标钟差误差影响造成的观测误差;为模糊度参数,为与频率有关的电离层模型延迟改正误差,为所述第i观测站与所述第j观测目标间的对流层路径延迟,为所述第i观测站与所述第j观测目标间的对流层模型延迟改正误差,为所述第i观测站与所述第j观测目标间的第一剩余误差,为所述第i观测站与所述第j观测目标间的第二剩余误差。
3.根据权利要求2所述的基于分区改正的星基增强***定位精度提高方法,其特征在于,在所述步骤S4中,根据所述伪距观测残差数据并利用公式(2)获取所述第i观测站的钟差均值
其中,n为所述第i观测站观测到的所述观测目标的总数。
4.根据权利要求3所述的基于分区改正的星基增强***定位精度提高方法,其特征在于,在所述步骤S5中,根据所述相位观测残差数据并利用公式(3)获取各所述观测区的模糊度归算值数据
其中,为扣除所述观测站钟差后的相位残差改正数,m为当前观测区内的所述观测站个数,为所述第i观测站的钟差均值,t为历元编号。
5.根据权利要求4所述的基于分区改正的星基增强***定位精度提高方法,其特征在于,在所述步骤S6中,根据所述残差数据和所述模糊度归算值数据并利用公式(4)获取各所述伪距分区改正数数据ΔPj(f)和相位分区改正数数据ΔLj(f)|t:
其中,t为历元编号,为所述观测站的钟差均值。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610407689.9A CN106125110B (zh) | 2016-06-12 | 2016-06-12 | 基于分区改正的星基增强***定位精度提高方法 |
US16/309,438 US20200096648A1 (en) | 2016-06-12 | 2017-06-12 | Zone correction-based method for improving the positioning accuracy in a satellite-based augmentation system |
PCT/CN2017/087922 WO2017215563A1 (zh) | 2016-06-12 | 2017-06-12 | 基于分区改正的星基增强***定位精度提高方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610407689.9A CN106125110B (zh) | 2016-06-12 | 2016-06-12 | 基于分区改正的星基增强***定位精度提高方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106125110A CN106125110A (zh) | 2016-11-16 |
CN106125110B true CN106125110B (zh) | 2018-08-24 |
Family
ID=57270067
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610407689.9A Active CN106125110B (zh) | 2016-06-12 | 2016-06-12 | 基于分区改正的星基增强***定位精度提高方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200096648A1 (zh) |
CN (1) | CN106125110B (zh) |
WO (1) | WO2017215563A1 (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106125110B (zh) * | 2016-06-12 | 2018-08-24 | 中国科学院上海天文台 | 基于分区改正的星基增强***定位精度提高方法 |
CN108780155B (zh) * | 2017-01-11 | 2021-03-09 | 华为技术有限公司 | 一种定位方法、终端及服务器 |
DE102017206275A1 (de) * | 2017-04-12 | 2018-10-18 | Robert Bosch Gmbh | Verfahren zum Betreiben eines Korrekturdienstsystems und Korrekturdienstsystem |
CN108376477B (zh) * | 2017-05-02 | 2021-04-16 | 山东大学 | 车路协同公交精准数据融合中心***及工作方法 |
CN108089207A (zh) * | 2017-06-06 | 2018-05-29 | 中国科学院光电研究院 | 一种基于单差电离层建模的nrtk增强定位方法 |
CN108089214B (zh) * | 2017-12-20 | 2021-06-15 | 北京卫星导航中心 | 一种卫星定位方法和卫星定位*** |
CN108415046B (zh) * | 2017-12-20 | 2021-01-05 | 中国科学院上海天文台 | 一种接收机导航定位的方法以及接收机 |
CN108563108A (zh) * | 2017-12-23 | 2018-09-21 | 航天恒星科技有限公司 | 一种gnss卫星钟健康状态监测方法 |
CN108318905B (zh) * | 2018-01-05 | 2021-12-17 | 北京北方联星科技有限公司 | 在智能手机上实现亚米级定位的方法及亚米定位智能手机 |
CN108761501A (zh) * | 2018-07-03 | 2018-11-06 | 千寻位置网络有限公司 | Ppp差分改正数的播发方法及***、定位终端 |
CN109613572B (zh) * | 2018-08-07 | 2023-03-31 | 中国科学院上海天文台 | 一种基于分区切换的星基增强用户定位方法 |
CN110007326B (zh) * | 2019-04-15 | 2022-06-21 | 中国电子科技集团公司第二十研究所 | 一种用于星基增强***的双频测距误差参数生成方法 |
CN113093237B (zh) * | 2020-01-09 | 2024-06-07 | 中移(上海)信息通信科技有限公司 | Ssr轨钟改正数质量因子实时评估方法、装置、设备及介质 |
CN111505742A (zh) * | 2020-04-29 | 2020-08-07 | 中国科学院国家空间科学中心 | 一种gnss电离层掩星数据气候研究的参数网格化方法及*** |
JP7183485B2 (ja) * | 2020-09-04 | 2022-12-05 | 三菱電機株式会社 | 測位補強信号配信装置および測位補強信号配信方法 |
CN113267792B (zh) * | 2021-05-26 | 2022-07-19 | 中国电子科技集团公司第五十四研究所 | 一种基于监测站网的电离层异常监测预警方法 |
CN116299586B (zh) * | 2022-12-07 | 2023-11-03 | 中国科学院上海天文台 | 基于广播星历的精密单点定位方法、接收机、设备和介质 |
CN115688196B (zh) * | 2022-12-26 | 2023-07-18 | 萨科(深圳)科技有限公司 | 一种基于互联网平台订单大数据的在线数据处理方法 |
CN116050483B (zh) * | 2023-01-10 | 2024-02-13 | 辽宁工程技术大学 | 一种ssa-bp神经网络与灰色模型卫星钟差预报方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7548196B2 (en) * | 2005-02-15 | 2009-06-16 | Fagan John E | Navigation system using external monitoring |
US7415354B2 (en) * | 2006-04-28 | 2008-08-19 | L-3 Communications Corporation | System and method for GPS acquisition using advanced tight coupling |
JP2010019728A (ja) * | 2008-07-11 | 2010-01-28 | Denso Corp | 測位方法および測位装置 |
US8174437B2 (en) * | 2009-07-29 | 2012-05-08 | Hemisphere Gps Llc | System and method for augmenting DGNSS with internally-generated differential correction |
CN102788985A (zh) * | 2012-08-31 | 2012-11-21 | 中国科学院国家授时中心 | 广域加权伪距差分对定位误差的改正方法 |
US20140077991A1 (en) * | 2012-09-20 | 2014-03-20 | California Institute Of Technology | Enhanced broadcast ephemeris for high accuracy assisted gps positioning |
CN104730549A (zh) * | 2013-12-24 | 2015-06-24 | 深圳艾科创新微电子有限公司 | 用于北斗卫星导航***的定位方法、装置和*** |
CN104407368A (zh) * | 2014-11-03 | 2015-03-11 | 中国人民解放军63961部队 | 基于北斗卫星导航***的高精度差分信息生成方法 |
CN104749595A (zh) * | 2015-03-19 | 2015-07-01 | 上海北伽导航科技有限公司 | 一种基于高精度地基增强***的定位服务方法及*** |
CN105044741B (zh) * | 2015-06-29 | 2019-03-08 | 中国科学院上海天文台 | 一种伪距相位综合广域差分改正值的求解方法 |
CN106125110B (zh) * | 2016-06-12 | 2018-08-24 | 中国科学院上海天文台 | 基于分区改正的星基增强***定位精度提高方法 |
-
2016
- 2016-06-12 CN CN201610407689.9A patent/CN106125110B/zh active Active
-
2017
- 2017-06-12 WO PCT/CN2017/087922 patent/WO2017215563A1/zh active Application Filing
- 2017-06-12 US US16/309,438 patent/US20200096648A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
WO2017215563A1 (zh) | 2017-12-21 |
CN106125110A (zh) | 2016-11-16 |
US20200096648A1 (en) | 2020-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106125110B (zh) | 基于分区改正的星基增强***定位精度提高方法 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
AU2018426707B2 (en) | Precise point position and real-time kinematic (PPP-RTK) positioning method and device | |
EP3035080B1 (en) | Navigation satellite system positioning involving the generation of correction information | |
CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
CN111190203B (zh) | 基于非差观测模型和消秩亏的ppp-rtk定位方法 | |
Geng et al. | Triple-frequency GPS precise point positioning with rapid ambiguity resolution | |
US6799116B2 (en) | GPS correction methods, apparatus and signals | |
CN108196272A (zh) | 一种基于实时精密单点定位的卫星导航定位装置及方法 | |
KR101360918B1 (ko) | 무선항법 신호 처리방법 | |
CN110007320B (zh) | 网络rtk解算方法 | |
CN105353391A (zh) | 一种支持多类型定位终端的多网融合定位增强***及方法 | |
Cao et al. | The wide-area difference system for the regional satellite navigation system of COMPASS | |
CN105044741B (zh) | 一种伪距相位综合广域差分改正值的求解方法 | |
CN107861131A (zh) | 一种斜路径电离层延迟的获取方法及*** | |
CN111983641B (zh) | 一种用于实时生成北斗星基增强***完好性参数的方法 | |
CN104316943B (zh) | 一种伪距离和多普勒组合差分定位***及方法 | |
CN114236587A (zh) | 基于北斗地基增强的网络rtk解算方法及存储介质 | |
CN113703021A (zh) | 一种基于码伪距的秒级实时高精度定位方法与*** | |
Landau et al. | Trimble’s RTK and DGPS solutions in comparison with precise point positioning | |
Yuan et al. | The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS | |
CN104898144A (zh) | 一种北斗***参考站整周模糊度单历元确定方法 | |
Liu et al. | Real-time multi-constellation precise point positioning with integer ambiguity resolution | |
Liu et al. | Beidou wide-area augmentation system clock error correction and performance verification | |
CN109613572B (zh) | 一种基于分区切换的星基增强用户定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |