CN113759402A - Gnss形变监测方法、*** - Google Patents
Gnss形变监测方法、*** Download PDFInfo
- Publication number
- CN113759402A CN113759402A CN202010496023.1A CN202010496023A CN113759402A CN 113759402 A CN113759402 A CN 113759402A CN 202010496023 A CN202010496023 A CN 202010496023A CN 113759402 A CN113759402 A CN 113759402A
- Authority
- CN
- China
- Prior art keywords
- difference
- double
- satellite
- observation equation
- error
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 241001061260 Emmelichthys struhsakeri Species 0.000 claims description 7
- 238000007667 floating Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000001934 delay Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000005433 ionosphere Substances 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000005436 troposphere Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 description 1
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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/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
-
- 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
-
- 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
本申请公开了GNSS形变监测方法、***。该方法包括:根据双差观测方程计算卫星的双差残差;将所述双差残差转换为单差残差;计算每颗卫星的重复周期;根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;在基线解算过程中消除所述每颗卫星的多路径误差。
Description
技术领域
本说明书一般涉及卫星导航技术领域,尤其涉及一种GNSS形变监测方法、***。
背景技术
目前基于GNSS的监测方案都是通过监测站数据与基准站数据(或虚拟的VRS数据)组成基线,进行基线解算获得监测点的坐标信息,从而确定建构筑物是否发生形变。首先,地面接收的卫星数据不但包含接收机位置等有用信息,还包含随机观测误差和各种***误差。其中***误差,例如电离层误差、对流层误差、轨道误差、卫星钟差和接收机钟差等,可通过双频观测、双差观测及精确模型修正得到消除或大大削弱。但是与站址环境相关的多路径效应不具空间相关性,无法采用差分模式消除,也无法通过线性模型估计来分离。
发明内容
本说明书提供了一种GNSS形变监测方法,消除卫星导航定位信号传播过程中的多路径误差。
本申请的一实施例中公开了一种GNSS形变监测方法,包括:
根据双差观测方程计算卫星的双差残差;
将所述双差残差转换为单差残差;
计算每颗卫星的重复周期;
根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
在基线解算过程中消除所述每颗卫星的多路径误差。
在一个优选例中,所述根据双差观测方程计算卫星的双差残差的步骤,进一步包括:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,分别表示如下:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
根据最小二乘原理解算获得所述待估参数的估计值为坐标参数浮点解,为双差模糊度参数浮点解,利用lambda方法对所述双差模糊度浮点解进行固定,获得双差模糊度参数固定解将所述双差模糊度参数固定解代入所述双差观测方程,根据最小二乘原理可获得坐标参数固定解将所述坐标参数固定解所述双差模糊度参数固定解代入所述双差观测方程计算获得所述双差残差
在一个优选例中,所述将卫星的双差残差转换为单差残差的步骤,进一步包括:
利用矩阵求逆将所述双差残差转换到仅包括多路径误差和随机噪声的单差残差。
在一个优选例中,所述计算每颗卫星的重复周期的步骤,进一步包括:
根据以下公式计算所述每颗卫星的重复周期:
Ta1=86400-2(2π/n)
Ta2=86400-(2π/n)
Ta3=86400·7-13·(2π/n)
其中为地球引力常数,a为卫星轨道的长半轴,n为卫星转动的角速度,Δn为平均角速度的摄动参数,Ta1为GPS卫星重复周期时间偏移,Ta2为BDS GEO/IGSO卫星重复周期时间偏移,Ta3为BDS MEO卫星重复周期时间偏移。
在一个优选例中,所述根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差的步骤,进一步包括:对相邻前几个周期的单差残差进行平均处理。
在一个优选例中,进一步包括:选取作为网平差时的观测值的基线向量;形成网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。
本申请的另一实施例中公开了一种GNSS形变监测***,包括:
残差获取模块,被配置为根据双差观测方程计算卫星的双差残差;
残差转换模块,被配置为将所述双差残差转换为单差残差;
周期计算模块,被配置为计算每颗卫星的重复周期;
合并计算模块,被配置为根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
误差提取模块,被配置为对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
解算模块,被配置为在基线解算过程中消除所述每颗卫星的多路径误差。
在一个优选例中,所述残差获取模块进一步被配置为:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,分别表示如下:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
根据最小二乘原理解算获得所述待估参数的估计值为坐标参数浮点解,为双差模糊度参数浮点解,利用lambda方法对所述双差模糊度浮点解进行固定,获得双差模糊度参数固定解将所述双差模糊度参数固定解代入所述双差观测方程,根据最小二乘原理可获得坐标参数固定解将所述坐标参数固定解所述双差模糊度参数固定解代入所述双差观测方程计算获得所述双差残差
在一个优选例中,所述残差转换模块进一步被配置为:
利用矩阵求逆将所述双差残差转换到仅包括多路径误差和随机噪声的单差残差。
在一个优选例中,所述周期计算模块根据以下公式计算所述每颗卫星的重复周期:
Ta1=86400-2(2π/n)
Ta2=86400-(2π/n)
Ta3=86400·7-13·(2π/n)
其中为地球引力常数,a为卫星轨道的长半轴,n为卫星转动的角速度,Δn为平均角速度的摄动参数,Ta1为GPS卫星重复周期时间偏移,Ta2为BDS GEO/IGSO卫星重复周期时间偏移,Ta3为BDS MEO卫星重复周期时间偏移。
在一个优选例中,所述合并计算模块对相邻前几个周期的单差残差进行平均处理。
在一个优选例中,进一步包括基线网平差计算模块,被配置为:选取作为网平差时的观测值的基线向量;形成网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。
本申请的另一实施例还公开了一种GNSS形变监测***包括:
存储器,用于存储计算机可执行指令;以及
处理器,与所述存储器耦合,用于在执行所述计算机可执行指令时实现如前文描述的方法中的步骤。
本申请的另一实施例还公开了一种计算机可读存储介质所述计算机可读存储介质中存储有计算机可执行指令,所述计算机可执行指令被处理器执行时实现如前文描述方法中的步骤。
相对于现有技术,本申请具有以下有益效果:
本申请可有效的削弱多路径效应对基线解算的影响,提升单条基线的解算精度,同时采用GNSS基线网平差方案,消除误差基线对监测网的影响,改善监测网质量,提升整体监测精度,适用于大坝桥梁、危房、铁塔等复杂场景的毫米级形变监测。
本说明书中记载了大量的技术特征,分布在各个技术方案中,如果要罗列出本申请所有可能的技术特征的组合(即技术方案)的话,会使得说明书过于冗长。为了避免这个问题,本说明书上述发明内容中公开的各个技术特征、在下文各个实施方式和例子中公开的各技术特征、以及附图中公开的各个技术特征,都可以自由地互相组合,从而构成各种新的技术方案(这些技术方案均应该视为在本说明书中已经记载),除非这种技术特征的组合在技术上是不可行的。例如,在一个例子中公开了特征A+B+C,在另一个例子中公开了特征A+B+D+E,而特征C和D是起到相同作用的等同技术手段,技术上只要择一使用即可,不可能同时采用,特征E技术上可以与特征C相组合,则,A+B+C+D的方案因技术不可行而应当不被视为已经记载,而A+B+C+E的方案应当视为已经被记载。
附图说明
参考以下附图描述本申请的非限制性和非穷举性实施例,其中除非另有说明,否则相同的附图标记在各个视图中指代相同的部分。
图1是根据本说明书一实施例中GNSS形变监测方法的流程图。
图2是根据本说明书一实施例中GNSS形变监测***的示意图。
具体实施方式
现有的GNSS形变监测方案有以下两种。一是新建GNSS形变监测基准站:在建构筑物被测量点上安装监测基站,并在建构筑物附近相对无沉降区设置基准站,将所述监测基站与基准站数据通过无线或有线传输***传输给服务***进行分析解算获得形变信息。二是在监测点接收增强信号:接收机在监测点接收卫星电文信号及低轨卫星导航定位辅助增强***信号或导航定位及形变监测定位辅助增强***信号,如NRTK中VRS信号,进行分析解算获得形变信息。
上述两种监测方案都是通过监测站数据与基准站数据(或虚拟的VRS数据)组成基线,进行基线解算获得监测点的坐标信息,从而确定建构筑物是否发生形变。地面接收的卫星数据不但包含接收机位置等有用信息,还包含随机观测误差和各种***误差,如电离层误差、对流层误差、轨道误差、卫星钟差和接收机钟差等。
方案一中监测站与基准站距离往往较短,短基线情况下双差观测可有效消除或大大削弱电离层误差、对流层误差、轨道误差、卫星钟差和接收机钟差,此时估计参数中将只考虑坐标参数项及模糊度参数项,通过模糊度搜索固定,可获得监测点坐标进而获得形变量。
方案二中是利用增强***如CORS***,通过综合各个参考站的观测信息,建立精确的误差模型来修正距离相关误差,在用户站(监测站)附近产生一个实时上不存在的虚拟基准站VRS,然后,利用虚拟基准站VRS和监测站的观测值进行差分,精确确定监测站的坐标继而获得形变信息。
方案一的具体方案如下:
卫星导航定位信号从播发、传播、直至接收机天线捕获阶段,会受到各种类型误差的影响。其中卫星端的误差包括硬件延迟、卫星钟差、卫星轨道误差等,传播过程中的电离层延迟、对流层延迟,以及由反射物等引起的多路径效应,接收机端的误差包括硬件延迟、接收机钟差等。因此载波相位和伪距非差观测距离方程可表示成如下形式(以米为单位):
上式中各参数的意义如下:P和L分别为伪距和载波相位的非差观测值,为卫星到测站的几何距离,Oi为轨道误差,c为真空中光速,δtu和δti分别为接收机和卫星钟差,为对流层延迟,为电离层延迟,Dpi和Dpu分别为卫星端和接收机端伪距硬件延迟,Dli和Dlu分别为卫星端和接收机端载波相位硬件延迟,λu为当前频点的载波相位波长,整周模糊度,和分布为伪距与载波相位多路径误差,和分别为伪距与载波相位残余误差。
假设在当前历元下,u,r同时观测到了卫星i,将非差观测方程站间做差可得到单差观测方程,公式如下:
上式中短基线情况下站间单差卫星端的轨道误差和钟误差基本消除,同时会削弱具有强相关性的大部分电离层延迟和对流程延迟误差,多路径误差由于测站之间不具有相关性,因此站间单差后多路径误差仍然存在,此外站间单差接收机端的钟差和硬件延迟没有被消除。
假设某一历元,u,r测站同时观测到卫星i,j,则可以根据上式组成双差观测方程,具体表现形式如下:
上式中的多路径误差项无法忽略,此时估计参数中有坐标参数项、多路径误差项及模糊度参数项。并且,当前GNSS形变监测方案多采用单基线模式,即一个监测站对应一个基准站,以一条基线解算结果获取此监测点形变量,受观测环境影响单基线结果往往无法保证结果的可靠性和稳定性。本申请中提出一种顾及多路径误差的GNSS形变监测方案。
在以下的叙述中,为了使读者更好地理解本申请而提出了许多技术细节。但是,本领域的普通技术人员可以理解,即使没有这些技术细节和基于以下各实施方式的种种变化和修改,也可以实现本申请所要求保护的技术方案。为使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请的实施方式作进一步地详细描述。
本申请的一实施例中公开了一种GNSS形变监测方法,图1是本实施例中的GNSS形变监测方法的流程图,该方法包括:
步骤101,根据双差观测方程计算卫星的双差残差。
在一个优选例中,根据双差观测方程计算卫星的双差残差的步骤101进一步包括:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,具体的,
测站u观测到了卫星i,载波相位和伪距非差观测距离方程可表示成如下形式(以米为单位):
其中,P和L分别为伪距和载波相位的非差观测值,为卫星到测站的几何距离,Oi为轨道误差,c为真空中光速,δtu和δti分别为接收机和卫星钟差,为对流层延迟,为电离层延迟,Dpi和Dpu分别为卫星端和接收机端伪距硬件延迟,Dli和Dlu分别为卫星端和接收机端载波相位硬件延迟,λu为当前频点的载波相位波长,整周模糊度,和分布为伪距与载波相位多路径误差,和分别为伪距与载波相位残余误差。
测站u,r同时观测到了卫星i,卫星i的单差观测值采用如下公式表示:
测站u,r同时观测到卫星i,j,卫星i,j的双差观测值采用如下公式表示:
根据上述公式可得:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
根据最小二乘原理解算获得所述待估参数的估计值为坐标参数浮点解,为双差模糊度参数浮点解,再利用lambda方法对双差模糊度参数浮点解进行固定,获得双差模糊度参数固定解将双差模糊度参数固定解代入所述双差观测方程,根据最小二乘原理可获得坐标参数固定解将坐标参数固定解双差模糊度参数固定解代入所述双差观测方程计算获得所述双差残差值。
步骤102,将卫星的双差残差转换为单差残差。
式中wi表示i卫星在基线解算中的权因子,假设观测卫星足够多,观测值中的多路径误差可视为随机噪声,单差残差加权平均值理论上应该为零,即上式中利用矩阵求逆可将双差残差转换到只剩下多路径误差和随机噪声的单差残差。
步骤103,计算每颗卫星的重复周期。
在一个优选例中,根据以下公式计算所述每颗卫星的重复周期:
Ta1=86400-2(2π/n)
Ta2=86400-(2π/n)
Ta3=86400·7-13·(2π/n)
其中为地球引力常数,a为卫星轨道的长半轴,n为卫星转动的角速度,Δn为平均角速度的摄动参数,Ta1为GPS卫星重复周期时间偏移,Ta2为BDS GEO/IGSO卫星重复周期时间偏移,Ta3为BDS MEO卫星重复周期时间偏移。
例如,BDS卫星、GEO卫星由于其是地球同步轨道卫星,相对地球固定不动,其运行周期为一个恒星日,且IGSO卫星运行周期也为一个恒星日,因此GEO卫星、IGSO卫星的重复周期可认为是一个恒星日。对于MEO卫星其运行周期不具半个恒星日或者一个恒星日的特性,而根据其运行周期和恒星日的最小公倍数中可以分析出,当其运行7个恒星日周期后,MEO卫星将围绕地球运行13圈,所以MEO卫星其重复周期为7个恒星日。
步骤104,根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差。
在一个优选例中,所述根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差的步骤,进一步包括:对相邻前几个周期的单差残差进行平均处理。
本实施例中,GNSS形变监测观测条件较差,由于遮挡等原因,数据存在缺失,同时在数据解算过程中往往剔除数据质量较差的数据,导致最后获得的卫星单差残差存在数据缺失,为保证一个周期内卫星残差数据的完整性,获取卫星相邻3个周期数据进行合并均值化处理,获得较为完整的单个周期残差数据。
步骤105,对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差。
步骤106,在基线解算过程中消除所述每颗卫星的多路径误差,也就是在观测方程中减去计算的多路径误差。
本实施例中的GNSS形变监测方法利用相邻前几个周期的数据提取多路径误差,用于当前数据改正,可实现多路径误差的实时改正。同时为确保建立的多路径模型完整性,采用了相邻的前3个周期的卫星单差残差提取多路径误差,提升基线解算精度;
在一个优选例中,所述GNSS形变监测***进一步包括:选取作为网平差时的观测值的基线向量;形成无约束网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。本实施例中采用GNSS基线网平差方案,消除误差基线对监测网的影响,改善监测网质量,提升整体监测精度,适用于大坝桥梁、危房、铁塔等复杂场景的毫米级形变监测。
本申请的另一实施例中公开了一种GNSS形变监测***,图2示出了本实施例中GNSS形变监测***的框图,该***包括:
残差获取模块201,被配置为根据双差观测方程计算卫星的双差残差;
残差转换模块202,被配置为将所述双差残差转换为单差残差;
周期计算模块203,被配置为计算每颗卫星的重复周期;
合并计算模块204,被配置为根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
误差提取模块205,被配置为对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
解算模块206,被配置为在基线解算过程中消除所述每颗卫星的多路径误差。
在一个优选例中,所述残差获取模块101进一步被配置为:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,具体的,
测站u观测到了卫星i,载波相位和伪距非差观测距离方程可表示成如下形式(以米为单位):
其中,P和L分别为伪距和载波相位的非差观测值,为卫星到测站的几何距离,Oi为轨道误差,c为真空中光速,δtu和δti分别为接收机和卫星钟差,为对流层延迟,为电离层延迟,Dpi和Dpu分别为卫星端和接收机端伪距硬件延迟,Dli和Dlu分别为卫星端和接收机端载波相位硬件延迟,λu为当前频点的载波相位波长,整周模糊度,和分布为伪距与载波相位多路径误差,和分别为伪距与载波相位残余误差。
测站u,r同时观测到了卫星i,卫星i的单差观测值采用如下公式表示:
u,r测站同时观测到卫星i,j,卫星i,j的双差观测值采用如下公式表示:
根据上述公式可得:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
根据最小二乘原理解算获得所述待估参数的估计值为坐标参数浮点解,为双差模糊度参数浮点解,再利用lambda方法对双差模糊度参数浮点解进行固定,获得双差模糊度参数固定解将双差模糊度参数固定解代入所述双差观测方程,根据最小二乘原理可获得坐标参数固定解将坐标参数固定解双差模糊度参数固定解代入所述双差观测方程计算获得所述双差残差值;
式中wi表示i卫星在基线解算中的权因子,假设观测卫星足够多,观测值中的多路径误差可视为随机噪声,单差残差加权平均值理论上应该为零,即上式中利用矩阵求逆可将双差残差转换到只剩下多路径误差和随机噪声的单差残差。
在一个优选例中,所述周期计算模块根据以下公式计算所述每颗卫星的重复周期:
Ta1=86400-2(2π/n)
Ta2=86400-(2π/n)
Ta3=86400·7-13·(2π/n)
在一个优选例中,所述合并计算模块对相邻前几个周期的单差残差进行平均处理。
在一个优选例中,进一步包括基线网平差计算模块,被配置为:选取作为网平差时的观测值的基线向量;形成无约束网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。
第一实施方式是与本实施方式相对应的***实施方式,第一实施方式中的技术细节可以应用于本实施方式,本实施方式中的技术细节也可以应用于第一实施方式。
需要说明的是,本领域技术人员应当理解,上述GNSS形变监测***的实施方式中所示的各模块的实现功能可参照前述GNSS形变监测方法的相关描述而理解。上述在GNSS形变监测***的实施方式中所示的各模块的功能可通过运行于处理器上的程序(可执行指令)而实现,也可通过具体的逻辑电路而实现。本说明书实施例上述GNSS形变监测***如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本说明书实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本说明书各个实施例所述方法的全部或部分。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。这样,本说明书实施例不限制于任何特定的硬件和软件结合。
相应地,本说明书实施方式还提供一种计算机可读存储介质,其中存储有计算机可执行指令,该计算机可执行指令被处理器执行时实现本说明书的各方法实施方式。计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括但不限于,相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读存储介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
此外,本说明书实施方式还提供一种GNSS形变监测***,其中包括用于存储计算机可执行指令的存储器,以及,处理器;该处理器用于在执行该存储器中的计算机可执行指令时实现上述各方法实施方式中的步骤。
在一个实施例中,该计算机可执行指令可以用于:
根据双差观测方程计算卫星的双差残差;
将卫星的双差残差转换为单差残差;
计算每颗卫星的重复周期;
根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
在基线解算过程中消除所述每颗卫星的多路径误差。
在一个实施例中,其中,该处理器可以是中央处理单元(Central ProcessingUnit,简称“CPU”),还可以是其他通用处理器、数字信号处理器(Digital SignalProcessor,简称“DSP”)、专用集成电路(Application Specific Integrated Circuit,简称“ASIC”)等。前述的存储器可以是只读存储器(read-only memory,简称“ROM”)、随机存取存储器(random access memory,简称“RAM”)、快闪存储器(Flash)、硬盘或者固态硬盘等。本发明各实施方式所公开的方法的步骤可以直接体现为硬件处理器执行完成,或者用处理器中的硬件及软件模块组合执行完成。在一个实施例中,该海面监控***的***还包括总线和通信接口。处理器、存储器和通信接口都通过总线相互连接。通信接口可以是无线通信接口也可以是有线通信接口,用于使得处理器能够与其他的***通信。
需要说明的是,在本专利的申请文件中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。本专利的申请文件中,如果提到根据某要素执行某行为,则是指至少根据该要素执行该行为的意思,其中包括了两种情况:仅根据该要素执行该行为、和根据该要素和其它要素执行该行为。多个、多次、多种等表达包括2个、2次、2种以及2个以上、2次以上、2种以上。
在本说明书提及的所有文献都被认为是整体性地包括在本说明书的公开内容中,以便在必要时可以作为修改的依据。此外应理解,以上所述仅为本说明书的较佳实施例而已,并非用于限定本说明书的保护范围。凡在本说明书一个或多个实施例的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本说明书一个或多个实施例的保护范围之内。
在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描述的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
Claims (14)
1.一种GNSS形变监测方法,其特征在于,包括:
根据双差观测方程计算卫星的双差残差;
将所述双差残差转换为单差残差;
计算每颗卫星的重复周期;
根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
在基线解算过程中消除所述每颗卫星的多路径误差。
2.如权利要求1所述的GNSS形变监测方法,其特征在于,所述根据双差观测方程计算卫星的双差残差的步骤,进一步包括:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,分别表示如下:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
5.如权利要求1所述的GNSS形变监测方法,其特征在于,所述根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差的步骤,进一步包括:对相邻前几个周期的单差残差进行平均处理。
6.如权利要求1所述的GNSS形变监测方法,其特征在于,进一步包括:选取作为网平差时的观测值的基线向量;形成无约束网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。
7.一种GNSS形变监测***,其特征在于,包括:
残差获取模块,被配置为根据双差观测方程计算卫星的双差残差;
残差转换模块,被配置为将所述双差残差转换为单差残差;
周期计算模块,被配置为计算每颗卫星的重复周期;
合并计算模块,被配置为根据每颗卫星相邻前几个周期的单差残差合并计算一个周期的单差残差;
误差提取模块,被配置为对所述合并计算的单差残差进行经验模态分析,获取每颗卫星的多路径误差;
解算模块,被配置为在基线解算过程中消除所述每颗卫星的多路径误差。
8.如权利要求7所述的GNSS形变监测***,其特征在于,所述残差获取模块进一步被配置为:
构建双差观测方程,所述双差观测方程包括伪距双差观测方程和载波相位双差观测方程,分别表示如下:
对所述双差观测方程进行线性化处理得到以下公式:
V=A·x+B·y-L
其中,V为双差观测值残差,A、B分别为坐标参数及双差模糊度参数系数阵,L为观测值与计算值之差;
11.如权利要求7所述的GNSS形变监测***,其特征在于,所述合并计算模块对相邻前几个周期的单差残差进行平均处理。
12.如权利要求7所述的GNSS形变监测***,其特征在于,进一步包括基线网平差计算模块,被配置为:选取作为网平差时的观测值的基线向量;形成无约束网平差的函数模型;求解所述函数模型,得到所述基线向量的平差值、改正数及相应的精度统计值;判断所述观测值是否存在粗差;若是,将所述观测值剔除或降权,若不是,输出最终结果。
13.一种GNSS形变监测***,其特征在于,包括:
存储器,用于存储计算机可执行指令;
处理器,与所述存储器耦合,用于在执行所述计算机可执行指令时实现如权利要求1至6中任意一项所述的方法中的步骤。
14.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机可执行指令,所述计算机可执行指令被处理器执行时实现如权利要求1至6中任意一项所述方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010496023.1A CN113759402A (zh) | 2020-06-03 | 2020-06-03 | Gnss形变监测方法、*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010496023.1A CN113759402A (zh) | 2020-06-03 | 2020-06-03 | Gnss形变监测方法、*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113759402A true CN113759402A (zh) | 2021-12-07 |
Family
ID=78783385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010496023.1A Pending CN113759402A (zh) | 2020-06-03 | 2020-06-03 | Gnss形变监测方法、*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113759402A (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101218354B1 (ko) * | 2012-03-05 | 2013-01-03 | 주식회사 이제이텍 | Gnss 기반 교량 모니터링 과대오차 처리 방법 |
CN104880717A (zh) * | 2015-05-27 | 2015-09-02 | 中国科学院嘉兴微电子与***工程中心 | 一种卫星测量值多路径误差检测装置及其算法 |
KR20160030644A (ko) * | 2014-09-11 | 2016-03-21 | 창신대학교 산학협력단 | 다중경로 오차 보정 방법 |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
CN107314763A (zh) * | 2017-07-18 | 2017-11-03 | 上海海洋大学 | 一种基于制约函数非线性估计的卫星影像区域网平差方法 |
CN109738917A (zh) * | 2018-12-30 | 2019-05-10 | 广州海达安控智能科技有限公司 | 一种北斗变形监测中的多路径误差削弱方法及装置 |
KR101981225B1 (ko) * | 2018-12-03 | 2019-05-22 | 한국해양과학기술원 | Gnss 측정치에 포함된 특이신호요소 정밀추정 방법 |
CN110006332A (zh) * | 2019-04-19 | 2019-07-12 | 中国矿业大学(北京) | 一种地震前兆形变gnss监测网的三角形形变监测方法 |
CN111007537A (zh) * | 2019-10-21 | 2020-04-14 | 桂林电子科技大学 | 一种卫星信号模拟器的时间同步***、方法和卫星信号模拟器 |
-
2020
- 2020-06-03 CN CN202010496023.1A patent/CN113759402A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101218354B1 (ko) * | 2012-03-05 | 2013-01-03 | 주식회사 이제이텍 | Gnss 기반 교량 모니터링 과대오차 처리 방법 |
KR20160030644A (ko) * | 2014-09-11 | 2016-03-21 | 창신대학교 산학협력단 | 다중경로 오차 보정 방법 |
CN104880717A (zh) * | 2015-05-27 | 2015-09-02 | 中国科学院嘉兴微电子与***工程中心 | 一种卫星测量值多路径误差检测装置及其算法 |
CN106569242A (zh) * | 2016-11-04 | 2017-04-19 | 中国科学院测量与地球物理研究所 | 固定参考星的gnss单差处理方法 |
CN107229061A (zh) * | 2017-07-18 | 2017-10-03 | 武汉大学 | 一种基于低轨卫星的星地差分实时精密定位方法 |
CN107314763A (zh) * | 2017-07-18 | 2017-11-03 | 上海海洋大学 | 一种基于制约函数非线性估计的卫星影像区域网平差方法 |
KR101981225B1 (ko) * | 2018-12-03 | 2019-05-22 | 한국해양과학기술원 | Gnss 측정치에 포함된 특이신호요소 정밀추정 방법 |
CN109738917A (zh) * | 2018-12-30 | 2019-05-10 | 广州海达安控智能科技有限公司 | 一种北斗变形监测中的多路径误差削弱方法及装置 |
CN110006332A (zh) * | 2019-04-19 | 2019-07-12 | 中国矿业大学(北京) | 一种地震前兆形变gnss监测网的三角形形变监测方法 |
CN111007537A (zh) * | 2019-10-21 | 2020-04-14 | 桂林电子科技大学 | 一种卫星信号模拟器的时间同步***、方法和卫星信号模拟器 |
Non-Patent Citations (3)
Title |
---|
李仁安;: "基于单差短基线的北斗观测值域多路径误差滤波算法", 工程勘察, no. 12 * |
石强;戴吾蛟;曾凡河;张超;: "BDS多路径效应特征及其对静态基线解精度的影响", 大地测量与地球动力学, no. 10 * |
袁兵等: "顾及多路径误差改正的 GNSS 大坝形变监测研究", 导航定位与授时, vol. 3, no. 1, pages 53 - 57 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646538B (zh) | 一种基于单差滤波的变形监测gnss信号多路径改正方法 | |
CN109738917B (zh) | 一种北斗变形监测中的多路径误差削弱方法及装置 | |
CN111045034B (zh) | 基于广播星历的gnss多***实时精密时间传递方法及*** | |
AU2008260578B2 (en) | Distance dependant error mitigation in real-time kinematic (RTK) positioning | |
CN104502935B (zh) | 一种基于非差非组合模型的网络rtk模糊度解算方法 | |
CN108828640B (zh) | 一种卫星导航定位观测值定权方法及装置 | |
US10338232B2 (en) | Navigation satellite wide-lane bias determination system and method | |
CN108196284B (zh) | 一种进行星间单差模糊度固定的gnss网数据处理方法 | |
WO2011046866A1 (en) | System and method for compensating for faulty measurements | |
CN109143298B (zh) | 北斗和gps观测值周跳探测与修复方法、设备及存储设备 | |
CN107966722B (zh) | 一种gnss钟差解算方法 | |
CN112068161B (zh) | 一种多路径误差削减方法及装置 | |
CN113253314B (zh) | 一种低轨卫星间时间同步方法及*** | |
CN111290004A (zh) | 伪距差分定位方法、装置、电子设备及存储介质 | |
CN111522032B (zh) | 一种北斗三号***用户完好性处理的优化方法及优化装置 | |
Matviichuk et al. | Estimating ocean tide loading displacements with GPS and GLONASS | |
CN111929705A (zh) | 基于动态基准站的定位方法、装置、电子设备及存储介质 | |
CN111123315A (zh) | 非差非组合ppp模型的优化方法及装置、定位*** | |
Zeng et al. | GPS triple-frequency undifferenced and uncombined precise orbit determination with the consideration of receiver time-variant bias | |
Yin et al. | A novel cycle slips detection model for the high precision positioning | |
CN115993620B (zh) | 模糊度固定方法及其*** | |
CN114002722A (zh) | Gnss定位方法及其装置 | |
Kuang et al. | Galileo real-time orbit determination with multi-frequency raw observations | |
CN113759402A (zh) | Gnss形变监测方法、*** | |
Yang et al. | Real-time kinematic GPS positioning for centimeter level ocean surface monitoring |
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 |