CN105301619A - 大规模gnss网数据整体快速处理方法及系统 - Google Patents

大规模gnss网数据整体快速处理方法及系统 Download PDF

Info

Publication number
CN105301619A
CN105301619A CN201510875063.6A CN201510875063A CN105301619A CN 105301619 A CN105301619 A CN 105301619A CN 201510875063 A CN201510875063 A CN 201510875063A CN 105301619 A CN105301619 A CN 105301619A
Authority
CN
China
Prior art keywords
survey station
poor
float ambiguities
upd
network data
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
Application number
CN201510875063.6A
Other languages
English (en)
Inventor
姜卫平
陈华
葛茂荣
周晓慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201510875063.6A priority Critical patent/CN105301619A/zh
Publication of CN105301619A publication Critical patent/CN105301619A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier 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

本发明公开了一种大规模GNSS网数据整体快速处理方法及系统,包括:步骤1,采用整网解法对参考站网数据进行定轨;步骤2,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;步骤3,对GNSS网中测站逐一进行PPP定位和模糊度固定;步骤4,将相位观测值转换成高精度距离观测值;步骤5,根据高精度距离观测值对GNSS网进行整网解算。本发明中,模糊度在单站解算过程中已被固定,整网解算时则不需要估计或仅估计未固定的模糊度参数,如此可大幅度减少待估参数,提高计算效率,使得大规模GNSS网数据可同时解算;同时,本发明可改善数据连续性,能提供稳定可靠且高精度的整网解结果。

Description

大规模GNSS网数据整体快速处理方法及系统
技术领域
本发明属于卫星导航定位与应用技术领域,尤其涉及一种大规模GNSS网数据整体快速处理方法及系统。
背景技术
随着卫星导航定位卫星系统(GNSS)的不断推广应用以及新卫星系统的建设和发展,地面GNSS跟踪站也越来越多,如IGS(Beutleretal.,1994b;Dowetal.,2009)全球跟踪网测站已达460个(Neilanetal.,2013),美国国家CORS网已有近2000个(SnayandSoler,2008),日本地震监测网测站数已超过1300多个(Sagiyaetal.,2004),中国CORS站也已远超过1000个(Dangetal.,2011)。将所有的测站、所有卫星的数据在观测值层面统一处理(下称整网解),不仅能充分利用各个系统之间的优势,获取更稳定、可靠、统一、兼容的多系统GNSS产品(etal.,2011;Montenbrucketal.,2013),还可为参考框架的建立和维持(Altamimietal.,2008)、验潮站海平面变化分析(etal.,2009)以及冰后回弹(LarsonandvanDam,2012)等提供不可或缺的完整的方差和协方差阵。
然而,由于待解参数过多,目前的整网解方法计算负担重,计算效率低。同时,这导致几乎目前所有的GNSS软件都难以甚至无法同时解算含有数百甚至上千测站的网,如国际著名的GAMIT软件不能同时解算超过100个测站的网。对于包含数百个站的GNSS网,当前常采用的解决办法是将整网划分为若干个子网,每个子网单独解算,最后将各子网的结果合并,得到一个整网解。但是为了将各个子网联合在一起,各子网间必须有一定数目的公共观测站,由于这些公共站被使用两次或两次以上,它最终得到的方差协方差阵与直接整网解的结果不完全一致,这将可能降低后续分析的精度和可靠性。因此研究一种高效率的GNSS数据处理方法,解决大规模GNSS网数据处理效率低的难题是十分有意义的。
发明内容
针对现有技术存在的不足,本发明提供了一种大规模GNSS网数据整体快速处理方法及系统,可高效、高精度地处理包含数百个甚至上千个测站的大规模GNSS网数据。
为解决上述技术问题,本发明采用如下的技术方案:
一种大规模GNSS网数据整体快速处理方法,包括:
步骤1,从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
步骤3,以步骤1获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用步骤2所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
步骤4,将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
步骤5,利用高精度距离观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度。
另一种大规模GNSS网数据整体快速处理方法,包括:
步骤1,从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
步骤3,以步骤1获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用步骤2所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
步骤4,将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
步骤5,利用高精度距离观测值和接收机获得的伪距观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的测站的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度;同时,在高精度距离观测值上估计窄巷UPD。
上述参考站网包含分布于全球的80~100个测站。
上述步骤3中所述的模糊度固定进一步为:
3.1采用判定函数将可准确固定的非差浮点模糊度固定;
3.2更新余下的非差浮点模糊度,更新后执行子步骤3.1,直至无可准确固定的非差浮点模糊度;
3.3将固定了的非差浮点模糊度转化为单频率上的模糊度,即整数模糊度。
与上述一种大规模GNSS网数据整体快速处理方法对应的系统,包括:
定轨模块,用来从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
UPD解算模块,用来以定轨模块获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
PPP定位及模糊度固定模块,用来以定轨模块获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用UPD解算模块所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
高精度距离观测值生成模块,用来将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
整网解算模块,用来利用高精度距离观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度。
与上述另一种大规模GNSS网数据整体快速处理方法对应的系统,包括:
定轨模块,用来从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
UPD解算模块,用来以定轨模块获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
PPP定位及模糊度固定模块,用来以定轨模块获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用UPD解算模块所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
高精度距离观测值生成模块,用来将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
整网解算模块,用来利用高精度距离观测值和接收机获得的伪距观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的测站的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度;同时,在高精度距离观测值上估计窄巷UPD。
和现有技术相比,本发明具有如下优点和有益效果:
1、本发明中,模糊度在单站解算过程中已被固定,整网解算时则不需要估计或仅估计未固定的模糊度参数,如此可大幅度减少待估参数,提高计算效率,使得大规模GNSS网数据可同时整体解算。
2、本发明可改善数据连续性,进而能提供稳定可靠且高精度的整网解结果。
附图说明
图1为本发明方法的具体流程图;
图2为实施例选取的基准站分布图;
图3为本发明方法与传统方法处理GNSS网数据的处理时间对比图;
图4为本发明方法与传统方法得到的轨道的相邻天之间重复轨道的RMS值对比图。
具体实施方式
下面将结合附图对本发明具体实施方式做进一步说明。
见图1,本发明包括如下步骤:
步骤1,选取参考站网,采用整网解法对参考站网数据进行定轨,得卫星的精密轨道、精密钟差和非差浮点模糊度,。
为确定高精度的卫星轨道和钟差,需选取全球均匀分布的测站作为参考站,从而构建参考站网,参考站网包含参考站数目不能太少,否则难以获得高精度的轨道;也不能太多,否则计算效率低,耗时长,部分软件甚至无法解算。作为优选,参考站网应包含分布于全球的80~100个测站。
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD。
宽巷UPD即宽巷的小数偏差部分,窄巷UPD即窄巷的小数偏差部分。每颗卫星每天估计一个宽巷UPD值,每颗卫星每隔5~~15分钟估计一个窄巷UPD值。
本步骤可采用本领域内的现有技术实现,为便于理解,下面将提供宽巷UPD和窄巷UPD的估计过程。
宽巷UPD的计算中,各参考站上MW组合的小数部分是观测值,卫星端和接收机端的宽巷UPD是待估值。计算出宽巷UPD后,可利用宽巷UPD进行宽巷模糊度固定,得到整周宽巷模糊度。
窄巷UPD的估计是建立于宽巷UPD计算的基础上,首先利用步骤1中得到的非差浮点模糊度和整周宽巷模糊度计算窄巷模糊度,将窄巷模糊度小数部分分离出来作为观测值来计算卫星端和接收机端的窄巷UPD。
将计算的卫星端和接收机端的宽巷UPD和窄巷UPD分别保存为宽巷UPD文件和窄巷UPD文件。
步骤3,以步骤1获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位(精密单点定位),得各测站的非差浮点模糊度,利用步骤2所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度。
对测站进行PPP定位可得测站的非差浮点模糊度,利用卫星端的宽巷UPD和窄巷UPD对测站的非差浮点模糊度施加UPD改正,恢复非差浮点模糊度的整周特性,采用判定函数对非差浮点模糊度进行固定,得测站的整数模糊度。
为了固定更多的非差浮点模糊度,模糊度固定需要采用序贯迭代法。先采用判定函数将可准确固定的非差浮点模糊度固定,然后更新余下的非差浮点模糊度;更新后,再采用判定函数将可准确固定的非差浮点模糊度固定;重复上述步骤直至无可准确固定的非差浮点模糊度。之后,将固定了的非差浮点模糊度转化为单频率上的模糊度,即整数模糊度,并将其存储于测站相关的文件中。
本步骤是在单测站上完成,可利用多个CPU并行处理,以加快处理速度。
步骤4,将测站的整数模糊度叠加到该测站的相位观测值,得该测站的高精度距离观测值。
本步骤用来完成测站的相位观测值到高精度距离观测值的转换,相位观测值由接收机端获得。
转换过程中,若模糊度固定了,则将该测站的相位观测值转换为高精度距离观测值;若模糊度未固定,则不将该测站的相位观测值转换为高精度距离观测值。将高精度距离观测值保存于新RINEX文件。新RINEX文件中,仅模糊度固定了的测站的相位观测值转换成了高精度距离观测值,标记转换成高精度距离观测值的相位观测值。
对GNSS网中所有测站逐测站分别执行步骤3~4。
步骤5,整网解算。
利用新RINEX文件中保存的高精度距离观测值重新整网解算,获得测站的位置、精密轨道和精密钟差等参数。解算过程中,不估计模糊度已固定的测站的非差浮点模糊度,仅估计步骤3中模糊度未固定的测站的非差浮点模糊度,从而可大幅度减少待估参数,大大提高计算效率。
若同时使用高精度距离观测值和伪距观测值进行GNSS网整网解算,为保证两种观测值的钟差定义一致,高精度距离观测值上还需额外同时估计窄巷UPD参数。当仅使用高精度距离观测值时,窄巷UPD参数将完全被钟差参数吸收,此时则无需额外估计窄巷UPD参数。
具体实施时,整网解算可以采用两种模式,一是仅利用高精度距离观测值进行整网解算,解算公式见式(1)。这种模式下,卫星的精密钟差和窄巷UPD无法分离,窄巷UPD参数将完全被钟差参数吸收,不再需要进行估计。另一种是同时采用高精度距离观测值和伪距观测值进行整网解算,伪距观测值由接收机端获得。为保证两种观测值的精密钟差定义一致,需要同时在高精度距离观测值上额外估计窄巷UPD,见公式(2)。
L ~ c = f 1 2 f 1 2 - f 2 2 L ~ 1 - f 2 2 f 1 2 - f 2 2 L ~ 2 = ρ + m T + dt R ′ - dt S ′ dt R ′ = dt R + λ 1 f 1 f 1 + f 2 δb n r dt S ′ = dt S - λ 1 f 1 f 1 + f 2 δb n s - - - ( 1 )
L ~ c = f 1 2 f 1 2 - f 2 2 L ~ 1 - f 2 2 f 1 2 - f 2 2 L ~ 2 = ρ + dt R - dt S + m T + λ 1 f 1 f 1 + f 2 ( δb n r + δb n s ) P c = f 1 2 f 1 2 - f 2 2 P 1 - f 2 2 f 1 2 - f 2 2 P 2 = ρ + dt R - dt S + m T - - - ( 2 )
式(1)~(2)中:
表示无电离层组合的高精度距离观测值;
分别表示频率f1和f2上的高精度距离观测值,高精度距离观测值记载于新RINEX文件中;
ρ表示接收机和卫星间的几何距离,根据卫星和测站的位置获得;
dtR表示接收机钟差;
dtS表示卫星钟差;
T表示天顶对流层延迟;
m表示映射系数,由对流层映射函数给出;
λ1表示频率f1的波长,为GNSS系统提供的常量;
δbnr分别表示接收机端和卫星端的UPD;
dtR'、dtS'分别为重新定义的接收机钟差和卫星钟差,包含了接收机和卫星端的UPD值;
Pc表示无电离层的伪距观测值;
P1和P2分别表示频率f1和f2上的伪距观测值,由接收机获得。
实施例
见图2,本实施例选取了全球分布的约460个IGS基准站从2012年201天至289天的数据,其中约100个全球均匀分布的测站被选为参考站,见图2中三角形标记,被用来计算初始精密轨道、初始精密钟差和初始非差浮点模糊度。本实施例采用一台配有CPUi7(2.6GHZ)处理器、16GB内存、MACOS10.8系统的电脑上进行数据处理。
数据处理流程如下:
步骤1,利用参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度。
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,利用参考站网数据计算卫星端的宽巷UPD和窄巷UPD,并将宽巷UPD和窄巷UPD分别保存为宽巷UPD文件和窄巷UPD文件。
步骤3,以步骤1获得的精密轨道和精密钟差为已知值,并固定;对所有的460个测站逐测站进行PPP定位得各测站的非差浮点模糊度;采用步骤2所得宽巷UPD和窄巷UPD对各测站的非差浮点模糊度施加UPD改正,恢复模糊度整周特性,并加以固定。将固定的非差浮点模糊度,即测站的整数模糊度保存于与测站相关的文件。
步骤4,利用模糊度已固定的测站的整数模糊度,将RINEX格式的原始数据文件中保存的相位观测值转换为高精度距离观测值,并保存为新RINEX文件。转换过程中,仅将模糊度已固定的相位观测值转换为高精度距离观测值,模糊度未固定的相位观测值保持不变。
重复步骤3~4,直到所有模糊度可固定的测站均生成对应的新RINEX文件。
步骤5,对新RINEX文件中数据进行整网解算。
可采用两种处理模式进行整网解算。模式一,同时采用高精度距离观测值和伪距观测值,并同时估计卫星的轨道、钟差、对流层、测站坐标、UPD、地球自转参数等。模式二,仅采用高精度距离观测值,同时估计卫星的轨道、钟差、对流层、测站坐标和地球自转参数等。
为比较分析,同时采用目前常用方法和上述两种模式分别对含有100、150、200、250、300、350、400和460个测站的GNSS网数据进行处理,并记录处理时间。
本发明主要优势有两点:1)计算效率高,可同时解算测站数目多;2)能改善数据连续性,进而提高轨道结果精度。为验证这两点优势,将数据处理结果进行两种比较,一是比较处理效率,即将目前常用方法和本发明方法分别处理含有100、150、200、250、300、350、400和460个测站的GNSS网数据所需时间进行比较,见图3;二是将步骤1中利用约100个测站得到的轨道、步骤5中采用模式二处理同样约100个测站轨道以及步骤5中采用模式二处理所有460个测站得到的轨道进行比较,以分析利用本发明方法得到的轨道的精度。
(1)数据处理效率比较分析
GNSS网数据处理过程中,参数估计和残差编辑迭代进行,这种迭代耗时占据了大部分数据处理时间。一般而言,基于验后残差的编辑有三种模式:网模式(如GAMIT软件)、基线模式(如Bernese软件)和单站模式(如EPOS和PANDA软件)。当处理含数百个测站的GNSS网时,这种迭代将大大增加计算时间,尤其是在含有少数观测数据质量差的测站时,这种迭代次数还需要增加。因此为了提高效率,基于PPP的单站清理技术被用来处理大规模GPS网数据(Zhangetal.,2007),如IGSTIGA工程数据处理等。在本发明方法中,同样可以理解为其是通过PPP来进行数据清理。相对现有整网解而言,额外的步骤仅为UPD估计、PPP模糊度固定以及新RINEX产生三个步骤。UPD估计对于24天解数据仅需要约1分钟,而对于一个采样间隔为30s的数据,其他步骤仅需5s。更重要的是,这些步骤是基于单站处理,因而可以使用多台计算机同时并行处理。因此在本实验分析中,仅比较单次参数估计过程的时间来分析计算效率是合理的。
见图3,其中正方形标记表示目前常用方法,即传统方法;圆形标记表示本发明方法模式一,三角形标记表示本发明方法模式二,从图3可以看出,当GNSS网中测站数量较大时,本发明方法所需时间明显少于传统方法,当处理460个测站的数据时,使用本发明方法模式一和模式二的处理时间分别为16和14分钟,而传统方法处理时间则为82分钟,本发明方法数据处理效率高出常用方法近6倍。此外,本发明方法需要的计算时间随着测站数目的增长近似呈线性增长,而传统方法近似呈指数增长,这意味着当处理更多测站时,本发明方法的优势更加明显。
(2)轨道精度比较
相邻天之间重叠轨道弧段偏差的RMS值可作为衡量轨道精度的可靠标准。图4给出了采用传统方法处理全球约100个测站得到的轨道、采用本发明方法处理同样的100个测站得到的轨道以及采用本发明方法处理460个测站得到的轨道的相邻天之间的重叠轨道弧段偏差的RMS值。从图4可以看出,处理同样100个测站,采用本发明方法的结果明显优于传统方法的结果,尤其是那些在传统方法结果中较差卫的RMS减少显著,如G06、G09、G15和G27卫星的提高幅度分别为19%、15%、14%和26%。总体来说,所有卫星的平均重叠轨道弧段RMS值减小了约9.8%,从27.6mm减小为24.8mm。而采用本发明方法对460个测站数据处理得到的结果可进一步减少RMS到23.2mm,这也说明参与解算的测站越多,结果精度越好。
从实验可以看出,本发明方法不仅实现了大规模GNSS的整网解,还大幅度提高了大规模GNSS网的数据处理效率,同时还提高轨道的精度。

Claims (6)

1.一种大规模GNSS网数据整体快速处理方法,其特征是,包括:
步骤1,从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
步骤3,以步骤1获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用步骤2所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
步骤4,将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
步骤5,利用高精度距离观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度。
2.一种大规模GNSS网数据整体快速处理方法,其特征是,包括:
步骤1,从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
步骤2,以步骤1获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
步骤3,以步骤1获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用步骤2所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
步骤4,将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
步骤5,利用高精度距离观测值和接收机获得的伪距观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的测站的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度;同时,在高精度距离观测值上估计窄巷UPD。
3.如权利要求1或2所述的大规模GNSS网数据整体快速处理方法,其特征是:
所述的参考站网包含分布于全球的80~100个测站。
4.如权利要求1或2所述的大规模GNSS网数据整体快速处理方法,其特征是:
上述步骤3中所述的模糊度固定进一步为:
3.1采用判定函数将可准确固定的非差浮点模糊度固定;
3.2更新余下的非差浮点模糊度,更新后执行子步骤3.1,直至无可准确固定的非差浮点模糊度;
3.3将固定了的非差浮点模糊度转化为单频率上的模糊度,即整数模糊度。
5.一种大规模GNSS网数据整体快速处理系统,其特征是,包括:
定轨模块,用来从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
UPD解算模块,用来以定轨模块获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
PPP定位及模糊度固定模块,用来以定轨模块获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用UPD解算模块所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
高精度距离观测值生成模块,用来将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
整网解算模块,用来利用高精度距离观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度。
6.一种大规模GNSS网数据整体快速处理系统,其特征是,包括:
定轨模块,用来从测站中选取参考站,得参考站网,采用整网解法对参考站网数据进行定轨,获得卫星的精密轨道、精密钟差和非差浮点模糊度;
UPD解算模块,用来以定轨模块获得的非差浮点模糊度为初始非差浮点模糊度,根据参考站网数据估计卫星端的宽巷UPD和窄巷UPD;
PPP定位及模糊度固定模块,用来以定轨模块获得的精密轨道和精密钟差作为已知值,并固定;对GNSS网中测站逐一进行PPP定位,得各测站的非差浮点模糊度,利用UPD解算模块所得宽巷UPD和窄巷UPD对测站的非差浮点模糊度进行改正,得改正后的浮点非差模糊度,对改正后的浮点非差模糊度进行模糊度固定,得测站的整数模糊度;
高精度距离观测值生成模块,用来将模糊度已固定的整数模糊度叠加到该测站的相位观测值上,得该测站的高精度距离观测值,相位观测值由接收机端获得;
整网解算模块,用来利用高精度距离观测值和接收机获得的伪距观测值对GNSS网进行整网解算,解算中不估计模糊度已固定的测站的非差浮点模糊度,但需估计模糊度未固定的测站的非差浮点模糊度;同时,在高精度距离观测值上估计窄巷UPD。
CN201510875063.6A 2015-12-02 2015-12-02 大规模gnss网数据整体快速处理方法及系统 Pending CN105301619A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510875063.6A CN105301619A (zh) 2015-12-02 2015-12-02 大规模gnss网数据整体快速处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510875063.6A CN105301619A (zh) 2015-12-02 2015-12-02 大规模gnss网数据整体快速处理方法及系统

Publications (1)

Publication Number Publication Date
CN105301619A true CN105301619A (zh) 2016-02-03

Family

ID=55199095

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510875063.6A Pending CN105301619A (zh) 2015-12-02 2015-12-02 大规模gnss网数据整体快速处理方法及系统

Country Status (1)

Country Link
CN (1) CN105301619A (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107422354A (zh) * 2017-09-25 2017-12-01 武汉大学 一种模糊度固定的ppp/sins紧组合定位定姿方法
CN109581441A (zh) * 2018-12-18 2019-04-05 武汉大学 基于站间相关性空间结构函数构建的gnss成像方法
CN110398762A (zh) * 2019-07-15 2019-11-01 广州中海达卫星导航技术股份有限公司 实时钟差估计中的模糊度固定方法、装置、设备及介质
CN110874346A (zh) * 2018-08-30 2020-03-10 Sap欧洲公司 浮点值的压缩方案
CN112394376A (zh) * 2020-11-20 2021-02-23 中国人民解放军战略支援部队信息工程大学 一种大规模gnss网观测数据非差整网并行处理方法
CN112596088A (zh) * 2020-11-26 2021-04-02 成都网感科技有限公司 应用于土地测量的高精度定位方法、装置及存储介质
CN114879239A (zh) * 2022-07-11 2022-08-09 开普勒卫星科技(武汉)有限公司 一种增强瞬时ppp固定解的区域三频整数钟差估计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102353969A (zh) * 2011-09-02 2012-02-15 东南大学 精密单点定位技术中相位偏差的估计方法
WO2014065664A1 (en) * 2012-10-25 2014-05-01 Fugro N.V. Ppp-rtk method and system for gnss signal based position determination
CN104483690A (zh) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 一种gnss三频精密单点定位模糊度固定方法
CN104730548A (zh) * 2015-04-08 2015-06-24 武汉大学 一种基于导航卫星精密轨道和精密钟差的精密定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102353969A (zh) * 2011-09-02 2012-02-15 东南大学 精密单点定位技术中相位偏差的估计方法
WO2014065664A1 (en) * 2012-10-25 2014-05-01 Fugro N.V. Ppp-rtk method and system for gnss signal based position determination
CN104483690A (zh) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 一种gnss三频精密单点定位模糊度固定方法
CN104730548A (zh) * 2015-04-08 2015-06-24 武汉大学 一种基于导航卫星精密轨道和精密钟差的精密定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HUA CHEN ET.AL: ""An enhanced strategy for GNSS data processing of massive networks"", 《JOURNAL OF GEODESY》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107422354A (zh) * 2017-09-25 2017-12-01 武汉大学 一种模糊度固定的ppp/sins紧组合定位定姿方法
CN107422354B (zh) * 2017-09-25 2019-06-25 武汉大学 一种模糊度固定的ppp/sins紧组合定位定姿方法
CN110874346A (zh) * 2018-08-30 2020-03-10 Sap欧洲公司 浮点值的压缩方案
CN110874346B (zh) * 2018-08-30 2023-09-29 Sap欧洲公司 浮点值的压缩方案
CN109581441A (zh) * 2018-12-18 2019-04-05 武汉大学 基于站间相关性空间结构函数构建的gnss成像方法
CN109581441B (zh) * 2018-12-18 2019-11-08 武汉大学 基于站间相关性空间结构函数构建的gnss成像方法
CN110398762A (zh) * 2019-07-15 2019-11-01 广州中海达卫星导航技术股份有限公司 实时钟差估计中的模糊度固定方法、装置、设备及介质
CN112394376A (zh) * 2020-11-20 2021-02-23 中国人民解放军战略支援部队信息工程大学 一种大规模gnss网观测数据非差整网并行处理方法
CN112596088A (zh) * 2020-11-26 2021-04-02 成都网感科技有限公司 应用于土地测量的高精度定位方法、装置及存储介质
CN112596088B (zh) * 2020-11-26 2024-02-02 成都网感科技有限公司 应用于土地测量的高精度定位方法、装置及存储介质
CN114879239A (zh) * 2022-07-11 2022-08-09 开普勒卫星科技(武汉)有限公司 一种增强瞬时ppp固定解的区域三频整数钟差估计方法

Similar Documents

Publication Publication Date Title
CN105301619A (zh) 大规模gnss网数据整体快速处理方法及系统
CN102288978B (zh) 一种cors基站周跳探测与修复方法
CN107356947A (zh) 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN103728643B (zh) 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法
CN103176188B (zh) 一种区域地基增强ppp-rtk模糊度单历元固定方法
CN101403790B (zh) 单频gps接收机的精密单点定位方法
CN107942346B (zh) 一种高精度gnss电离层tec观测值提取方法
CN103323867B (zh) 一种北斗二代geo卫星精密定轨系统偏差改进方法
CN105301617B (zh) 一种卫星导航系统中的整周模糊度有效性检验方法
CN104502935A (zh) 一种基于非差非组合模型的网络rtk模糊度解算方法
CN104375157B (zh) 短基线下惯导辅助的北斗单频整周模糊度求解方法
CN103760572A (zh) 一种基于区域cors的单频ppp电离层加权方法
CN103941272A (zh) Gps、glonass和bds联合解算的定位方法
CN102998681A (zh) 一种卫星导航系统的高频钟差估计方法
CN105158782A (zh) 一种bds和gps观测信息融合的宽巷模糊度解算方法
CN103293550B (zh) 利用单频gnss接收机的实时高精度地震形变监测方法
RU2010146231A (ru) Способ и устройство для определения положения с помощью гибридных данных об орбите sps
CN106850177A (zh) 一种基于精密单点定位的实时时间传递方法
CN104614741A (zh) 一种不受glonass码频间偏差影响的实时精密卫星钟差估计方法
CN109212562A (zh) 一种生成载波伪距观测量的方法
CN104111467A (zh) 一种基于北斗三频宽巷组合的网络rtk瞬时定位方法
CN105182388A (zh) 一种快速收敛的精密单点定位方法
CN105204048A (zh) 一种北斗兼容gps/glonass的网络rtk参考站间模糊度快速解算方法
CN112394376B (zh) 一种大规模gnss网观测数据非差整网并行处理方法
CN103197335A (zh) 采用改进正则化方法抑制dgps整周模糊度病态性的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160203