CN115166788B - 一种低仰角对流层误差修正方法、系统、设备及介质 - Google Patents
一种低仰角对流层误差修正方法、系统、设备及介质 Download PDFInfo
- Publication number
- CN115166788B CN115166788B CN202210805902.7A CN202210805902A CN115166788B CN 115166788 B CN115166788 B CN 115166788B CN 202210805902 A CN202210805902 A CN 202210805902A CN 115166788 B CN115166788 B CN 115166788B
- Authority
- CN
- China
- Prior art keywords
- observation
- carrier phase
- satellite
- ionosphere
- combined
- 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
- 239000005436 troposphere Substances 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012937 correction Methods 0.000 title claims abstract description 25
- 239000005433 ionosphere Substances 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000004590 computer program Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000001934 delay Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012876 topography 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/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/04—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing carrier phase data
-
- 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
- G01S19/072—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/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
-
- 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
本发明公开了一种低仰角对流层误差修正方法、系统、设备及介质,涉及卫星定位数据处理技术领域,方法包括:获取针对两个频点的第一原始观测数据,第一原始观测数据包括多个卫星位置和两个频点中每个频点在各个卫星仰角对应的第二原始观测数据,基于第二原始观测数据,探测及修复周跳,得到目标观测数据;基于目标观测数据,确定无电离层组合观测方程,消除电离层延迟误差;根据卫星位置,确定对流层延迟误差;根据无电离层组合观测方程和对流层延迟误差,更新目标观测数据中的每个载波相位观测值,构建多个初始载波相位观测方程,根据各个卫星仰角,对各个初始载波相位观测方程进行定权,根据定权后的各个初始载波相位观测方程,确定接收机位置。
Description
技术领域
本发明涉及卫星定位数据处理技术领域,具体涉及一种低仰角对流层误差修正方法、系统、设备及介质。
背景技术
卫星信号为电磁波信号,由卫星端发出,穿过地球表面的大气层被接收机获取,其中大气层中主要影响信号传播的有电离层和对流层。大气延迟是卫星定位中的主要误差,是影响GPS定位精度的主要因素之一。如果不能有效地消除或削弱大气延迟,最终的定位偏差将达到上百米。
GPS定位过程中,无线电信号经过中性大气层时会发生对流层延迟,在天顶方向,对流层延迟量约为2.4m,当仰角较低时,其误差可达20m。由于对流层延迟与信号频率无关,因此无法像电离层延迟一样采用多频组合观测值进行消除,需要建立相应的对流层延迟修正模型。在测站高程不变的情况下,卫星仰角越低,对流层传播延迟和弯曲延迟会逐渐增大,在观测数据中会存在一个较大的延迟误差。经研究表明,卫星仰角在接近0°时,对流层弯曲延迟约为0.2m,对流层传播延迟将达到25m以上。对流层传播路径上的斜延迟很难直接获取,因此一般将对流层斜延迟投影到天顶对流层延迟,天顶对流层延迟的精度则影响着定位的精度。目前,常用的天顶方向对流层改正模型有Saastamoinen模型和Hopfield模型,但单独使用这两种模型中的任一个以修正天顶对流层延迟,最终的定位精度都不是很高。
GPS非差相位精密单点定位是利用GPS卫星精密星历及由一定方法确定的精密卫星钟差,以单台双频GPS接收机采集的非差相位数据作为主要观测值来进行单点定位计算,其精度可达分米级甚至厘米级。对于GPS非差相位精密单点定位,许多误差(如接收机钟差、卫星钟差、对流层延迟、电离层延迟)不能直接去掉,只能通过模型和设未知数进行消除,因此,采用好的随机模型对估算对流层天顶延迟非常关键。卫星仰角与多路径效应和噪声密切相关,传统的随机模型有sine模型、sine平方型,这类模型认为信号在低仰角方向的噪声比较大,天顶方向的噪声比较小,但忽略了信号在低仰角时包含的信息多,而随机模型遗漏的信息都会影响天顶延迟的结果和精度。
北斗卫星导航系统偏差(简称系统偏差)具有随仰角呈非线性变化的特性,对非线性模型的拟合主要采用多项式拟合,采用多项式拟合非线性模型的方法大致为:首先选择多项式基函数,然后利用最小二乘法计算多项式系数,最后通过单个多项式对系统偏差进行拟合。然而,在大量数据的情况下,使用单个多项式对系统偏差进行拟合的精度较低,且不同仰角下观测数据的质量不同(如低仰角卫星受到的多路径效应以及观测噪声的影响较大),更是降低了采用单个多项式对系统偏差进行拟合的精度。
实际上,北斗卫星导航系统定位过程中,系统偏差占主要部分,系统偏差包括多路径效应、电离层影响、对流层影响等,这些误差对于不同测站、不同观测方向上的卫星来说都是不等的,因此每颗卫星的测量精度也不同。然而,利用现有的非线性模型进行拟合,以修复系统偏差过程中,对所有的观测数据都是等权看待的,导致在进行最小二乘法拟合时,模型对接收机位置的定位精度会降低。
目前,我国正在全力建设北斗卫星导航系统,但卫星数量有限,且存在地形因素导致接收机受遮挡严重,以及由于地区原因观测到的北斗卫星大多为低仰角卫星,在这种情况下,低仰角北斗卫星定位成为了一大热点,因此,如何更有效地利用低仰角卫星并在某种程度上缓和低仰角下延迟误差大的问题,以提高对接收机位置的定位精度,是一个亟待解决的问题。
发明内容
本发明所要解决的技术问题为:现有的利用非线性模型进行拟合,以修正北斗卫星导航系统偏差过程中,对所有的观测数据都是等权看待的,导致模型对接收机位置的定位精度低。为解决该技术问题,本发明提供了一种低仰角对流层误差修正方法、系统、设备及介质。
本发明解决上述技术问题的技术方案如下:
一种低仰角对流层误差修正方法,包括:
步骤S1:获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B,每个所述卫星位置各自对应一个卫星仰角,所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角;
步骤S2:根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据;
步骤S3:根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程;
步骤S4:根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差;
步骤S5:根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值;
步骤S6:根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程;
步骤S7:根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程;
步骤S8:根据各个所述目标载波相位观测方程,确定接收机位置。
采用本发明的有益效果是:通过先探测第一观测数据中产生周跳的目标历元,更新目标历元对应的载波相位观测值,以得到一个准确的载波相位观测值,然后通过确定无电离层组合观测方程和对流层延迟误差,以消除目标观测数据中各个载波相位观测值包含的电离层延迟误差和对流层延迟误差,得到精度高的载波相位观测值,再根据精度高的载波相位观测值构建初始载波相位观测方程,并根据各个精度高的载波相位观测值对应的卫星仰角,对各个精度高的载波相位观测值对应的初始载波相位观测方程赋予权重,得到目标载波观测方程,从而解算接收机位置,解决了现有的利用非线性模型进行拟合,以修正北斗卫星导航系统偏差过程中,对所有的观测数据等权看待,导致模型对接收机位置的定位精度低的问题,通过对各个卫星仰角下目标观测数据中的载波相位观测值赋予权重,以修正系统偏差,从而提高了对接收机的定位精度。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,具体包括:
对于各个所述卫星仰角中的每个所述卫星仰角,若所述卫星仰角大于预设的分界角度,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重为设定值;
若所述卫星仰角不大于所述分界角度,根据所述卫星仰角,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重。
采用上述进一步方案的有益效果是:通过对各个卫星仰角下初始载波相位观测方程赋予权重,根据各个赋予了权重的目标载波相位观测方程以修正系统偏差,解算接收机位置,避免了对所有的观测数据等权看待,从而提高了对接收机定位精度。
进一步,所述步骤S2中,根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值,具体包括:
根据各个所述频点的伪距观测值和各个所述频点对应的频率值,确定电离层伪距组合观测值;
根据各个所述频点的载波相位观测值和各个所述频点对应的频率值,确定电离层载波相位组合观测值。
采用上述进一步方案的有益效果是:通过确定电离层伪距组合观测值和电离层载波相位组合观测值,再基于电离层伪距组合观测值和电离层载波相位组合观测值探测和修复第一原始观测数据中产生周跳的目标历元,并对目标历元对应的载波相位观测值进行修复,从而得到一个准确的载波相位观测值,为后续解算接收机位置,提高接收机位置精度奠定了基础。
进一步,所述步骤S3中,根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值,具体包括:
根据所述目标观测数据,通过第一公式,确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,其中,所述第一公式为:
其中,PIF表示无电离层伪距组合观测值,ΦIF表示无电离层载波相位组合观测值,f1、f2分别表示所述频点A、所述频点B对应的频率值,P1、P2分别表示所述频点A、所述频点B的伪距观测值,ρ表示卫星信号发射时刻与接收机的真实几何距离,c表示真空中的光速,dtr表示接收机时钟误差,表示卫星时钟误差,T表示信号在大气传播中的对流层延迟误差,Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值,λIF表示组合波长,NIF表示组合模糊度,/>表示伪距包括的硬件延迟误差和伪距噪声,/>表示载波相位包括的硬件延迟误差和相位噪声。
采用上述进一步方案的有益效果是:通过确定无电离层伪距组合观测值和无电离层载波相位组合观测值,以确定无电离层组合观测方程,从而消除了电离层延迟误差,提高了接收机定位精度。
进一步,所述对流层延迟误差是通过以下方式确定的:
选取卫星截止角,根据所述卫星截止角对应的卫星位置,通过Saastamoinen模型计算天顶延迟量;
通过映射函数将所述天顶延迟量投影至实际传播路径,得到对流层延迟误差,所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
采用上述进一步方案的有益效果是:通过确定对流层延迟误差,便于修正对流层误差,提高接收机位置的定位精度。
为解决上述技术问题,本发明还提供了一种低仰角对流层误差修正系统,包括:
数据获取模块,用于获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B,每个所述卫星位置各自对应一个卫星仰角,所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角;
第一数据处理模块,用于根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据;
第二数据处理模块,用于根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程;
第三数据处理模块,用于根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差;
第四数据处理模块,根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值;
第五数据处理模块,用于根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程;
第六数据处理模块,用于根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程;
解算定位模块,用于根据各个所述目标载波相位观测方程,确定接收机位置。
为解决上述技术问题,本发明还提供了一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如前所述的低仰角对流层误差修正方法。
为解决上述技术问题,本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如前所述的低仰角对流层误差修正方法。
附图说明
图1为本发明的方法流程示意图;
图2为本发明的系统结构示意图;
图3为本发明的电子设备结构示意图。
具体实施方式
以下对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
实施例一
为解决现有技术存在的技术问题,本实施例提供了一种低仰角对流层误差修正方法,如图1所示,包括:
步骤S1:获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B,每个所述卫星位置各自对应一个卫星仰角,所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角。
步骤S2:根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据。
其中,所述步骤S2中,根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值,具体包括:
根据各个所述频点的伪距观测值和各个所述频点对应的频率值,确定电离层伪距组合观测值,所述电离层伪距组合观测值为:
其中,Piono表示电离层伪距组合观测值,IP1表示P1码上电离层延迟;
根据各个所述频点的载波相位观测值和各个所述频点对应的频率值,确定电离层载波相位组合观测值。
所述预先构建的宽巷组合观测方程为:
其中,Nw表示宽巷组合模糊度,表示宽巷组合观测值,λw表示组合波长。
步骤S3:根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程。
其中,所述步骤S3中,根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值,具体包括:
根据所述目标观测数据,通过第一公式,确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,其中,所述第一公式为:
其中,ΦIF表示无电离层伪距组合观测值,ΦIF表示无电离层载波相位组合观测值,f1、f2分别表示所述频点A、所述频点B对应的频率值,P1、P2分别表示所述频点A、所述频点B的伪距观测值,ρ表示卫星信号发射时刻与接收机的真实几何距离,c表示真空中的光速,表示接收机时钟误差,/>表示卫星时钟误差,T表示信号在大气传播中的对流层延迟误差,Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值,λIF表示组合波长,NIF表示组合模糊度,/>表示伪距包括的硬件延迟误差和伪距噪声,/>表示载波相位包括的硬件延迟误差和相位噪声。
步骤S4:根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差。
其中,所述对流层延迟误差是通过以下方式确定的:
选取卫星截止角,根据所述卫星截止角对应的卫星位置,通过Saastamoinen模型计算天顶延迟量;
通过映射函数将所述天顶延迟量投影至实际传播路径,得到对流层延迟误差,所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
本实施例中,选取的卫星截止角的取值大于等于5°,估计得到的天顶对流层延迟的RMS值较低,便于有效地利用低仰角的观测数据,提高对流层延迟误差的计算精度。
步骤S5:根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值。
步骤S6:根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程。
步骤S7:根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程。
其中,所述步骤S7中,根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,具体包括:
对于各个所述卫星仰角中的每个所述卫星仰角,若所述卫星仰角大于预设的分界角度,确定所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重为设定值;所述设定值通过第二公式确定,所述第二公式为:
Psl=1/sigm02
所述卫星仰角大于所述分界角度时,所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重均为Ps1,1/sigma2为设定值;
若所述卫星仰角不大于所述分界角度,根据所述卫星仰角,通过第三公式确定所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重;其中,所述第三公式为:
ps2=cos2(E)/sigma2
所述卫星仰角不大于所述分界角度时,所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重均为Ps2,其中,E表示卫星仰角。
步骤S8:根据各个所述目标载波相位观测方程,通过加权最小二乘法解算出接收机位置。
实施例二
为解决现有技术存在的技术问题,本实施例提供了一种低仰角对流层误差修正系统,如图2所示,包括:
数据获取模块,用于获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B,每个所述卫星位置各自对应一个卫星仰角,所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角;
第一数据处理模块,用于根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据;
第二数据处理模块,用于根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程;
第三数据处理模块,用于根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差;
第四数据处理模块,根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值;
第五数据处理模块,用于根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程;
第六数据处理模块,用于根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程;
解算定位模块,用于根据各个所述目标载波相位观测方程,确定接收机位置。
进一步,所述第六数据处理模块具体包括:
权重确定单元,具体用于:对于各个所述卫星仰角中的每个所述卫星仰角,若所述卫星仰角大于预设的分界角度,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重为设定值;若所述卫星仰角不大于所述分界角度,根据所述卫星仰角,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重。
进一步,所述第一数据处理模块具体包括:
第一观测值确定单元,用于根据各个所述频点的伪距观测值和各个所述频点对应的频率值,确定电离层伪距组合观测值;根据各个所述频点的载波相位观测值和各个所述频点对应的频率值,确定电离层载波相位组合观测值。
进一步,所述第二数据处理模块具体包括:
第二观测值确定单元,用于根据所述目标观测数据,通过第一公式,确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,其中,所述第一公式为:
其中,PIF表示无电离层伪距组合观测值,ΦIF表示无电离层载波相位组合观测值,f1、f2分别表示所述频点A、所述频点B对应的频率值,P1、P2分别表示所述频点A、所述频点B的伪距观测值,ρ表示卫星信号发射时刻与接收机的真实几何距离,c表示真空中的光速,dtr表示接收机时钟误差,表示卫星时钟误差,T表示信号在大气传播中的对流层延迟误差,Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值,λIF表示组合波长,NIF表示组合模糊度,/>表示伪距包括的硬件延迟误差和伪距噪声,/>表示载波相位包括的硬件延迟误差和相位噪声。
进一步,所述第三数据处理模块具体包括:
天顶延迟确定单元,用于选取卫星截止角,根据所述卫星截止角对应的卫星位置,通过Saastamoinen模型计算天顶延迟量;
对流层延迟误差确定单元,用于通过映射函数将所述天顶延迟量投影至实际传播路径,得到对流层延迟误差,所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
实施例三
为解决现有技术存在的技术问题,本实施例提供了一种电子设备,如图3所示,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如实施例一所述的低仰角对流层误差修正方法。
实施例四
为解决现有技术存在的技术问题,本实施例提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如实施例一所述的低仰角对流层误差修正方法。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (8)
1.一种低仰角对流层误差修正方法,其特征在于,包括:
步骤S1:获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B,每个所述卫星位置各自对应一个卫星仰角,所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角;
步骤S2:根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据;
步骤S3:根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程;
步骤S4:根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差;
步骤S5:根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值;
步骤S6:根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程;
步骤S7:根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程;
步骤S8:根据各个所述目标载波相位观测方程,确定接收机位置。
2.根据权利要求1所述的方法,其特征在于,所述根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,具体包括:
对于各个所述卫星仰角中的每个所述卫星仰角,若所述卫星仰角大于预设的分界角度,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重为设定值;
若所述卫星仰角不大于所述分界角度,根据所述卫星仰角,确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重。
3.根据权利要求1所述的方法,其特征在于,所述步骤S2中,根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值,具体包括:
根据各个所述频点的伪距观测值和各个所述频点对应的频率值,确定电离层伪距组合观测值;
根据各个所述频点的载波相位观测值和各个所述频点对应的频率值,确定电离层载波相位组合观测值。
4.根据权利要求1所述的方法,其特征在于,所述步骤S3中,根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值,具体包括:
根据所述目标观测数据,通过第一公式,确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,其中,所述第一公式为:
其中,PIF表示无电离层伪距组合观测值,ΦIF表示无电离层载波相位组合观测值,f1、f2分别表示所述频点A、所述频点B对应的频率值,P1、P2分别表示所述频点A、所述频点B的伪距观测值,ρ表示卫星信号发射时刻与接收机的真实几何距离,c表示真空中的光速,dtr表示接收机时钟误差,表示卫星时钟误差,T表示信号在大气传播中的对流层延迟误差,Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值,λIF表示组合波长,NIF表示组合模糊度,/>表示伪距包括的硬件延迟误差和伪距噪声,/>表示载波相位包括的硬件延迟误差和相位噪声。
5.根据权利要求1所述的方法,其特征在于,所述对流层延迟误差是通过以下方式确定的:
选取卫星截止角,根据所述卫星截止角对应的卫星位置,通过Saastamoinen模型计算天顶延迟量;
通过映射函数将所述天顶延迟量投影至实际传播路径,得到对流层延迟误差,所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
6.一种低仰角对流层误差修正系统,其特征在于,包括:
数据获取模块,用于获取针对两个频点的第一原始观测数据,所述第一原始观测数据包括卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据,所述第二原始观测数据包括所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差,所述两个频点包括频点A和频点B;
第一数据处理模块,用于根据所述第二原始观测数据,确定电离层伪距组合观测值和电离层载波相位组合观测值;根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程,确定所述第一原始观测数据中产生周跳的目标历元,更新所述目标历元对应的载波相位观测值,得到目标观测数据;
第二数据处理模块,用于根据所述目标观测数据,确定无电离层伪距组合观测值和无电离层载波相位组合观测值;根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值,确定无电离层组合观测方程;
第三数据处理模块,用于根据所述卫星位置和预先构建的对流层改正模型,确定对流层延迟误差;
第四数据处理模块,用于根据所述无电离层组合观测方程和所述对流层延迟误差,更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值,得到每个所述频点在各个历元对应的调整后的载波相位观测值;
第五数据处理模块,用于根据每个所述频点在各个历元对应的调整后的载波相位观测值,构建每个所述频点在各个历元对应的初始载波相位观测方程;
第六数据处理模块,用于根据各个所述卫星仰角,对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权,得到多个目标载波相位观测方程;
解算定位模块,用于根据各个所述目标载波相位观测方程,确定接收机位置。
7.一种电子设备,其特征在于,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1至5中任一项所述的低仰角对流层误差修正方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至5中任一项所述的低仰角对流层误差修正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210805902.7A CN115166788B (zh) | 2022-07-08 | 2022-07-08 | 一种低仰角对流层误差修正方法、系统、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210805902.7A CN115166788B (zh) | 2022-07-08 | 2022-07-08 | 一种低仰角对流层误差修正方法、系统、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115166788A CN115166788A (zh) | 2022-10-11 |
CN115166788B true CN115166788B (zh) | 2024-04-26 |
Family
ID=83493477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210805902.7A Active CN115166788B (zh) | 2022-07-08 | 2022-07-08 | 一种低仰角对流层误差修正方法、系统、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115166788B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009008104A1 (ja) * | 2007-07-10 | 2009-01-15 | Electronic Navigation Research Institute, An Independent Administrative Institution | 天頂対流圏遅延量の算出方法及び衛星測位信号の対流圏遅延の補正方法 |
CN107153209A (zh) * | 2017-07-06 | 2017-09-12 | 武汉大学 | 一种短弧段低轨导航卫星实时精密定轨方法 |
WO2018034580A1 (en) * | 2016-08-15 | 2018-02-22 | Limited Liability Company "Topcon Positioning Systems" | Adaptive gnss positioning method capable of adjusting to ionospheric delay |
CN109782313A (zh) * | 2019-01-17 | 2019-05-21 | 上海卫星工程研究所 | 地面精密单点定位数据处理方法 |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
-
2022
- 2022-07-08 CN CN202210805902.7A patent/CN115166788B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009008104A1 (ja) * | 2007-07-10 | 2009-01-15 | Electronic Navigation Research Institute, An Independent Administrative Institution | 天頂対流圏遅延量の算出方法及び衛星測位信号の対流圏遅延の補正方法 |
WO2018034580A1 (en) * | 2016-08-15 | 2018-02-22 | Limited Liability Company "Topcon Positioning Systems" | Adaptive gnss positioning method capable of adjusting to ionospheric delay |
CN107153209A (zh) * | 2017-07-06 | 2017-09-12 | 武汉大学 | 一种短弧段低轨导航卫星实时精密定轨方法 |
CN109782313A (zh) * | 2019-01-17 | 2019-05-21 | 上海卫星工程研究所 | 地面精密单点定位数据处理方法 |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
Non-Patent Citations (3)
Title |
---|
GPS对流层延迟的历元间差分分析;李浩军;王解先;胡丛玮;;同济大学学报(自然科学版);20100315(03);全文 * |
基于双频数据组合的可见星电离层延迟研究;杨腾飞;高法钦;李斯凡;岳滕飞;;浙江理工大学学报(自然科学版);20160310(02);全文 * |
基于多观测方程融合的北斗卫星周跳探测方法;戴连君;唐涛;蔡伯根;;北京交通大学学报;20130415(02);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115166788A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN110045407B (zh) | 一种分布式伪卫星/gnss优化定位方法 | |
CA2681918A1 (en) | Distance dependant error mitigation in real-time kinematic (rtk) positioning | |
CN107607971B (zh) | 基于gnss共视时间比对算法的时间频率传递方法及接收机 | |
CN110456388B (zh) | 一种星载gnss-r海面高度要素定标装置及方法 | |
CN111551971B (zh) | 一种支持异频gnss信号伪距差分定位的方法 | |
CN104535976A (zh) | 一种相控阵传感器的卫星标校方法 | |
CN114994727B (zh) | 一种实现高精度时间校准及卫星定位的设备 | |
CN105510945A (zh) | 一种应用于卫导着陆外场检测的ppp定位方法 | |
CN103852799A (zh) | 一种基于ppp技术的船载重力测量方法 | |
CN115980317B (zh) | 基于修正相位的地基gnss-r数据土壤水分估算方法 | |
CN115166788B (zh) | 一种低仰角对流层误差修正方法、系统、设备及介质 | |
CN106814376B (zh) | 一种快速精确厘米级单点定位方法 | |
CN114779301B (zh) | 一种基于广播星历的卫星导航实时精密单点定位方法 | |
CN116466376A (zh) | 数值预报模式辅助的实时ppp改善方法 | |
CN113267793B (zh) | 一种基于外部增强信息的gbas对流层参数生成方法 | |
CN115421172A (zh) | 一种基于实时与准实时结合的北斗变形监测方法 | |
CN115144882B (zh) | 一种基于gnss和5g组合的精密单点定位方法及系统 | |
CN110764122A (zh) | 一种单频gps接收机的精密单点定位方法 | |
Wu et al. | Impact of elevation mask on multi-GNSS precise point positioning performance | |
CN113899957B (zh) | 基于双差相位观测值的室外天线绝对标定方法及系统 | |
CN116953746B (zh) | 基于单台相位失真的卫星导航天线的定向方法及装置 | |
CN113009525B (zh) | 一种建立实时对流层格网产品的方法 | |
CN115712130B (zh) | 一种伪距校正方法、装置、电子设备及存储介质 | |
CN111399019B (zh) | 一种gnss中卫星定位方法和装置 |
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 |