CN115166788A - 一种低仰角对流层误差修正方法、系统、设备及介质 - Google Patents

一种低仰角对流层误差修正方法、系统、设备及介质 Download PDF

Info

Publication number
CN115166788A
CN115166788A CN202210805902.7A CN202210805902A CN115166788A CN 115166788 A CN115166788 A CN 115166788A CN 202210805902 A CN202210805902 A CN 202210805902A CN 115166788 A CN115166788 A CN 115166788A
Authority
CN
China
Prior art keywords
observation
carrier phase
satellite
frequency point
ionospheric
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.)
Granted
Application number
CN202210805902.7A
Other languages
English (en)
Other versions
CN115166788B (zh
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.)
Xiangtan University
Original Assignee
Xiangtan University
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 Xiangtan University filed Critical Xiangtan University
Priority to CN202210805902.7A priority Critical patent/CN115166788B/zh
Publication of CN115166788A publication Critical patent/CN115166788A/zh
Application granted granted Critical
Publication of CN115166788B publication Critical patent/CN115166788B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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中，根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值，具体包括：
根据所述目标观测数据，通过第一公式，确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，其中，所述第一公式为：
Figure BDA0003737562670000051
Figure BDA0003737562670000052
其中，PIF表示无电离层伪距组合观测值，ΦIF表示无电离层载波相位组合观测值，f1、f2分别表示所述频点A、所述频点B对应的频率值，P1、P2分别表示所述频点A、所述频点B的伪距观测值，ρ表示卫星信号发射时刻与接收机的真实几何距离，c表示真空中的光速，dtr表示接收机时钟误差，
Figure BDA0003737562670000061
表示卫星时钟误差，T表示信号在大气传播中的对流层延迟误差，Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值，λIF表示组合波长，NIF表示组合模糊度，
Figure BDA0003737562670000062
表示伪距包括的硬件延迟误差和伪距噪声，
Figure BDA0003737562670000063
表示载波相位包括的硬件延迟误差和相位噪声。
采用上述进一步方案的有益效果是：通过确定无电离层伪距组合观测值和无电离层载波相位组合观测值，以确定无电离层组合观测方程，从而消除了电离层延迟误差，提高了接收机定位精度。
进一步，所述对流层延迟误差是通过以下方式确定的：
选取卫星截止角，根据所述卫星截止角对应的卫星位置，通过Saastamoinen模型计算天顶延迟量；
通过映射函数将所述天顶延迟量投影至实际传播路径，得到对流层延迟误差，所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
采用上述进一步方案的有益效果是：通过确定对流层延迟误差，便于修正对流层误差，提高接收机位置的定位精度。
为解决上述技术问题，本发明还提供了一种低仰角对流层误差修正系统，包括：
数据获取模块，用于获取针对两个频点的第一原始观测数据，所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据，所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差，所述两个频点包括频点A和频点B，每个所述卫星位置各自对应一个卫星仰角，所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角；
第一数据处理模块，用于根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值；根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程，确定所述第一原始观测数据中产生周跳的目标历元，更新所述目标历元对应的载波相位观测值，得到目标观测数据；
第二数据处理模块，用于根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值；根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，确定无电离层组合观测方程；
第三数据处理模块，用于根据所述卫星位置和预先构建的对流层改正模型，确定对流层延迟误差；
第四数据处理模块，根据所述无电离层组合观测方程和所述对流层延迟误差，更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值，得到每个所述频点在各个历元对应的调整后的载波相位观测值；
第五数据处理模块，用于根据每个所述频点在各个历元对应的调整后的载波相位观测值，构建每个所述频点在各个历元对应的初始载波相位观测方程；
第六数据处理模块，用于根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，得到多个目标载波相位观测方程；
解算定位模块，用于根据各个所述目标载波相位观测方程，确定接收机位置。
为解决上述技术问题，本发明还提供了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如前所述的低仰角对流层误差修正方法。
为解决上述技术问题，本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如前所述的低仰角对流层误差修正方法。
附图说明
图1为本发明的方法流程示意图；
图2为本发明的系统结构示意图；
图3为本发明的电子设备结构示意图。
具体实施方式
以下对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
实施例一
为解决现有技术存在的技术问题，本实施例提供了一种低仰角对流层误差修正方法，如图1所示，包括：
步骤S1：获取针对两个频点的第一原始观测数据，所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据，所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差，所述两个频点包括频点A和频点B，每个所述卫星位置各自对应一个卫星仰角，所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角。
步骤S2：根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值；根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程，确定所述第一原始观测数据中产生周跳的目标历元，更新所述目标历元对应的载波相位观测值，得到目标观测数据。
其中，所述步骤S2中，根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值，具体包括：
根据各个所述频点的伪距观测值和各个所述频点对应的频率值，确定电离层伪距组合观测值，所述电离层伪距组合观测值为：
Figure BDA0003737562670000091
其中，Piono表示电离层伪距组合观测值，IP1表示P1码上电离层延迟；
根据各个所述频点的载波相位观测值和各个所述频点对应的频率值，确定电离层载波相位组合观测值。
所述预先构建的宽巷组合观测方程为：
Figure BDA0003737562670000092
其中，Nw表示宽巷组合模糊度，
Figure BDA0003737562670000093
表示宽巷组合观测值，λw表示组合波长。
步骤S3：根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值；根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，确定无电离层组合观测方程。
其中，所述步骤S3中，根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值，具体包括：
根据所述目标观测数据，通过第一公式，确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，其中，所述第一公式为：
Figure BDA0003737562670000094
Figure BDA0003737562670000095
其中，ΦIF表示无电离层伪距组合观测值，ΦIF表示无电离层载波相位组合观测值，f1、f2分别表示所述频点A、所述频点B对应的频率值，P1、P2分别表示所述频点A、所述频点B的伪距观测值，ρ表示卫星信号发射时刻与接收机的真实几何距离，c表示真空中的光速，
Figure BDA0003737562670000096
表示接收机时钟误差，
Figure BDA0003737562670000097
表示卫星时钟误差，T表示信号在大气传播中的对流层延迟误差，Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值，λIF表示组合波长，NIF表示组合模糊度，
Figure BDA0003737562670000101
表示伪距包括的硬件延迟误差和伪距噪声，
Figure BDA0003737562670000102
表示载波相位包括的硬件延迟误差和相位噪声。
步骤S4：根据所述卫星位置和预先构建的对流层改正模型，确定对流层延迟误差。
其中，所述对流层延迟误差是通过以下方式确定的：
选取卫星截止角，根据所述卫星截止角对应的卫星位置，通过Saastamoinen模型计算天顶延迟量；
通过映射函数将所述天顶延迟量投影至实际传播路径，得到对流层延迟误差，所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
本实施例中，选取的卫星截止角的取值大于等于5°，估计得到的天顶对流层延迟的RMS值较低，便于有效地利用低仰角的观测数据，提高对流层延迟误差的计算精度。
步骤S5：根据所述无电离层组合观测方程和所述对流层延迟误差，更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值，得到每个所述频点在各个历元对应的调整后的载波相位观测值。
步骤S6：根据每个所述频点在各个历元对应的调整后的载波相位观测值，构建每个所述频点在各个历元对应的初始载波相位观测方程。
步骤S7：根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，得到多个目标载波相位观测方程。
其中，所述步骤S7中，根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，具体包括：
对于各个所述卫星仰角中的每个所述卫星仰角，若所述卫星仰角大于预设的分界角度，确定所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重为设定值；所述设定值通过第二公式确定，所述第二公式为：
Psl＝1/sigm02
所述卫星仰角大于所述分界角度时，所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重均为Ps1，1/sigma2为设定值；
若所述卫星仰角不大于所述分界角度，根据所述卫星仰角，通过第三公式确定所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重；其中，所述第三公式为：
ps2＝cos2(E)/sigma2
所述卫星仰角不大于所述分界角度时，所述卫星仰角下每个所述频点在各个历元的载波相位观测值的权重均为Ps2，其中，E表示卫星仰角。
步骤S8：根据各个所述目标载波相位观测方程，通过加权最小二乘法解算出接收机位置。
实施例二
为解决现有技术存在的技术问题，本实施例提供了一种低仰角对流层误差修正系统，如图2所示，包括：
数据获取模块，用于获取针对两个频点的第一原始观测数据，所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据，所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差，所述两个频点包括频点A和频点B，每个所述卫星位置各自对应一个卫星仰角，所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角；
第一数据处理模块，用于根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值；根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程，确定所述第一原始观测数据中产生周跳的目标历元，更新所述目标历元对应的载波相位观测值，得到目标观测数据；
第二数据处理模块，用于根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值；根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，确定无电离层组合观测方程；
第三数据处理模块，用于根据所述卫星位置和预先构建的对流层改正模型，确定对流层延迟误差；
第四数据处理模块，根据所述无电离层组合观测方程和所述对流层延迟误差，更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值，得到每个所述频点在各个历元对应的调整后的载波相位观测值；
第五数据处理模块，用于根据每个所述频点在各个历元对应的调整后的载波相位观测值，构建每个所述频点在各个历元对应的初始载波相位观测方程；
第六数据处理模块，用于根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，得到多个目标载波相位观测方程；
解算定位模块，用于根据各个所述目标载波相位观测方程，确定接收机位置。
进一步，所述第六数据处理模块具体包括：
权重确定单元，具体用于：对于各个所述卫星仰角中的每个所述卫星仰角，若所述卫星仰角大于预设的分界角度，确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重为设定值；若所述卫星仰角不大于所述分界角度，根据所述卫星仰角，确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重。
进一步，所述第一数据处理模块具体包括：
第一观测值确定单元，用于根据各个所述频点的伪距观测值和各个所述频点对应的频率值，确定电离层伪距组合观测值；根据各个所述频点的载波相位观测值和各个所述频点对应的频率值，确定电离层载波相位组合观测值。
进一步，所述第二数据处理模块具体包括：
第二观测值确定单元，用于根据所述目标观测数据，通过第一公式，确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，其中，所述第一公式为：
Figure BDA0003737562670000131
Figure BDA0003737562670000132
其中，PIF表示无电离层伪距组合观测值，ΦIF表示无电离层载波相位组合观测值，f1、f2分别表示所述频点A、所述频点B对应的频率值，P1、P2分别表示所述频点A、所述频点B的伪距观测值，ρ表示卫星信号发射时刻与接收机的真实几何距离，c表示真空中的光速，dtr表示接收机时钟误差，
Figure BDA0003737562670000133
表示卫星时钟误差，T表示信号在大气传播中的对流层延迟误差，Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值，λIF表示组合波长，NIF表示组合模糊度，
Figure BDA0003737562670000134
表示伪距包括的硬件延迟误差和伪距噪声，
Figure BDA0003737562670000135
表示载波相位包括的硬件延迟误差和相位噪声。
进一步，所述第三数据处理模块具体包括：
天顶延迟确定单元，用于选取卫星截止角，根据所述卫星截止角对应的卫星位置，通过Saastamoinen模型计算天顶延迟量；
对流层延迟误差确定单元，用于通过映射函数将所述天顶延迟量投影至实际传播路径，得到对流层延迟误差，所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
实施例三
为解决现有技术存在的技术问题，本实施例提供了一种电子设备，如图3所示，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如实施例一所述的低仰角对流层误差修正方法。
实施例四
为解决现有技术存在的技术问题，本实施例提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如实施例一所述的低仰角对流层误差修正方法。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (8)

1.一种低仰角对流层误差修正方法，其特征在于，包括：
步骤S1：获取针对两个频点的第一原始观测数据，所述第一原始观测数据包括多个卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据，所述第二原始观测数据包括每个所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差，所述两个频点包括频点A和频点B，每个所述卫星位置各自对应一个卫星仰角，所述卫星仰角表征卫星与接收机所在处的地平线之间的夹角；
步骤S2：根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值；根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程，确定所述第一原始观测数据中产生周跳的目标历元，更新所述目标历元对应的载波相位观测值，得到目标观测数据；
步骤S3：根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值；根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，确定无电离层组合观测方程；
步骤S4：根据所述卫星位置和预先构建的对流层改正模型，确定对流层延迟误差；
步骤S5：根据所述无电离层组合观测方程和所述对流层延迟误差，更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值，得到每个所述频点在各个历元对应的调整后的载波相位观测值；
步骤S6：根据每个所述频点在各个历元对应的调整后的载波相位观测值，构建每个所述频点在各个历元对应的初始载波相位观测方程；
步骤S7：根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，得到多个目标载波相位观测方程；
步骤S8：根据各个所述目标载波相位观测方程，确定接收机位置。
2.根据权利要求1所述的方法，其特征在于，所述根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，具体包括：
对于各个所述卫星仰角中的每个所述卫星仰角，若所述卫星仰角大于预设的分界角度，确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重为设定值；
若所述卫星仰角不大于所述分界角度，根据所述卫星仰角，确定所述卫星仰角下每个所述频点在各个历元的初始载波相位观测方程的权重。
3.根据权利要求1所述的方法，其特征在于，所述步骤S2中，根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值，具体包括：
根据各个所述频点的伪距观测值和各个所述频点对应的频率值，确定电离层伪距组合观测值；
根据各个所述频点的载波相位观测值和各个所述频点对应的频率值，确定电离层载波相位组合观测值。
4.根据权利要求1所述的方法，其特征在于，所述步骤S3中，根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值，具体包括：
根据所述目标观测数据，通过第一公式，确定所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，其中，所述第一公式为：
Figure FDA0003737562660000021
Figure FDA0003737562660000022
其中，PIF表示无电离层伪距组合观测值，ΦIF表示无电离层载波相位组合观测值，f1、f2分别表示所述频点A、所述频点B对应的频率值，P1、P2分别表示所述频点A、所述频点B的伪距观测值，ρ表示卫星信号发射时刻与接收机的真实几何距离，c表示真空中的光速，dtr表示接收机时钟误差，
Figure FDA0003737562660000031
表示卫星时钟误差，T表示信号在大气传播中的对流层延迟误差，Φ1、Φ2分别表示所述频点A、所述频点B的载波相位观测值，λIF表示组合波长，NIF表示组合模糊度，
Figure FDA0003737562660000033
表示伪距包括的硬件延迟误差和伪距噪声，
Figure FDA0003737562660000032
表示载波相位包括的硬件延迟误差和相位噪声。
5.根据权利要求1所述的方法，其特征在于，所述对流层延迟误差是通过以下方式确定的：
选取卫星截止角，根据所述卫星截止角对应的卫星位置，通过Saastamoinen模型计算天顶延迟量；
通过映射函数将所述天顶延迟量投影至实际传播路径，得到对流层延迟误差，所述实际传播路径表征卫星天线与接收机之间的信号传播路径。
6.一种低仰角对流层误差修正系统，其特征在于，包括：
数据获取模块，用于获取针对两个频点的第一原始观测数据，所述第一原始观测数据包括卫星位置和两个所述频点中每个所述频点在各个卫星仰角对应的第二原始观测数据，所述第二原始观测数据包括所述频点在各个历元的伪距观测值、载波相位观测值和卫星时钟误差，所述两个频点包括频点A和频点B；
第一数据处理模块，用于根据所述第二原始观测数据，确定电离层伪距组合观测值和电离层载波相位组合观测值；根据所述电离层伪距组合观测值、电离层载波相位组合观测值和预先构建的宽巷组合观测方程，确定所述第一原始观测数据中产生周跳的目标历元，更新所述目标历元对应的载波相位观测值，得到目标观测数据；
第二数据处理模块，用于根据所述目标观测数据，确定无电离层伪距组合观测值和无电离层载波相位组合观测值；根据所述无电离层伪距组合观测值和所述无电离层载波相位组合观测值，确定无电离层组合观测方程；
第三数据处理模块，用于根据所述卫星位置和预先构建的对流层改正模型，确定对流层延迟误差；
第四数据处理模块，用于根据所述无电离层组合观测方程和所述对流层延迟误差，更新所述目标观测数据中每个所述频点在各个历元的载波相位观测值，得到每个所述频点在各个历元对应的调整后的载波相位观测值；
第五数据处理模块，用于根据每个所述频点在各个历元对应的调整后的载波相位观测值，构建每个所述频点在各个历元对应的初始载波相位观测方程；
第六数据处理模块，用于根据各个所述卫星仰角，对各个所述卫星仰角下每个所述频点在各个历元对应的初始载波相位观测方程进行定权，得到多个目标载波相位观测方程；
解算定位模块，用于根据各个所述目标载波相位观测方程，确定接收机位置。
7.一种电子设备，其特征在于，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1至5中任一项所述的低仰角对流层误差修正方法。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至5中任一项所述的低仰角对流层误差修正方法。
CN202210805902.7A 2022-07-08 2022-07-08 一种低仰角对流层误差修正方法、系统、设备及介质 Active CN115166788B (zh)

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 true CN115166788A (zh) 2022-10-11
CN115166788B 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)

* Cited by examiner, † Cited by third party
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精密卫星钟差的方法及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
戴连君;唐涛;蔡伯根;: "基于多观测方程融合的北斗卫星周跳探测方法", 北京交通大学学报, no. 02, 15 April 2013 (2013-04-15) *
李浩军;王解先;胡丛玮;: "GPS对流层延迟的历元间差分分析", 同济大学学报(自然科学版), no. 03, 15 March 2010 (2010-03-15) *
杨腾飞;高法钦;李斯凡;岳滕飞;: "基于双频数据组合的可见星电离层延迟研究", 浙江理工大学学报(自然科学版), no. 02, 10 March 2016 (2016-03-10) *

Also Published As

Publication number Publication date
CN115166788B (zh) 2024-04-26

Similar Documents

Publication Publication Date Title
CN110045407B (zh) 一种分布式伪卫星/gnss优化定位方法
CN109581452B (zh) 一种gnss参考站载波相位整周模糊度解算方法
JP7153427B2 (ja) 衛星測位システムを用いた測位方法および測位装置
CN108802782B (zh) 一种惯导辅助的北斗三频载波相位整周模糊度求解方法
JP2010528320A (ja) リアルタイムキネマティック（ｒｔｋ）測位における距離依存性誤差の軽減
CN109669196B (zh) 一种顾及基线形变的多天线gnss载波相位精密测姿方法
CN110456388B (zh) 一种星载gnss-r海面高度要素定标装置及方法
CN114966760B (zh) 一种电离层加权的非差非组合ppp-rtk技术实现方法
CN103852799A (zh) 一种基于ppp技术的船载重力测量方法
CN114859389A (zh) 一种gnss多系统抗差自适应融合rtk解算方法
CN115061167A (zh) 一种适用于短距离大高差rtk的对流层延迟改正方法
Abdullah et al. Improving ambiguity resolution rate with an accurate ionospheric differential correction
CN115980317B (zh) 基于修正相位的地基gnss-r数据土壤水分估算方法
CN115166788B (zh) 一种低仰角对流层误差修正方法、系统、设备及介质
CN109143289B (zh) 一种gnss单站位移监测方法
CN111650612A (zh) 一种基于双差模式的实时精密时间传递方法
CN116203603A (zh) 一种面向复杂环境低成本的ppp随机模型自适应调整方法
CN116466376A (zh) 数值预报模式辅助的实时ppp改善方法
CN106814376A (zh) 一种快速精确厘米级单点定位方法
CN110941002B (zh) 一种自适应抗差的序贯最小二乘精密单点定位方法
CN110764123B (zh) 一种基于glonass广播星历的伪距定位改进方法
CN115144882B (zh) 一种基于gnss和5g组合的精密单点定位方法及系统
CN110764122A (zh) 一种单频gps接收机的精密单点定位方法
Wu et al. Impact of elevation mask on multi-GNSS precise point positioning performance
CN113899957B (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