CN116609799A - 一种厘米级斜向电离层延迟产品的生成方法及装置 - Google Patents

一种厘米级斜向电离层延迟产品的生成方法及装置 Download PDF

Info

Publication number
CN116609799A
CN116609799A CN202310892383.7A CN202310892383A CN116609799A CN 116609799 A CN116609799 A CN 116609799A CN 202310892383 A CN202310892383 A CN 202310892383A CN 116609799 A CN116609799 A CN 116609799A
Authority
CN
China
Prior art keywords
satellite
bias
oblique
delay
reference station
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
CN202310892383.7A
Other languages
English (en)
Other versions
CN116609799B (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.)
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 CN202310892383.7A priority Critical patent/CN116609799B/zh
Publication of CN116609799A publication Critical patent/CN116609799A/zh
Application granted granted Critical
Publication of CN116609799B publication Critical patent/CN116609799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S19/072—Ionosphere corrections
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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

Abstract

本发明涉及一种厘米级斜向电离层延迟产品的生成方法及装置，方法包括：利用非组合PPP技术提取出厘米级斜向电离层延迟信息；基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品。可避免电离层薄层假设以及投影函数误差带来的精度损失，从而提高电离层延迟信息的精度。

Description

一种厘米级斜向电离层延迟产品的生成方法及装置
技术领域
本发明涉及卫星导航定位技术领域，特别涉及一种厘米级斜向电离层延迟产品的生成方法及装置。
背景技术
作为全球导航卫星系统（GNSS）的一个主要误差源，电离层延迟的校正精度会直接影响到GNSS服务的可用性、准确性、可靠性和完整性。
目前为了满足GNSS的应用需求，相关技术涉及两类电离层改正方式：一类是根据电离层资料构建出的经验模型（IRI模型、Bent模型）改正，这类改正方式由于参数设置复杂，改正精度有限，所以不太适合直接应用于GNSS用户。另一类是利用电离层格网产品加以改正，利用若干格网点内插得到穿刺点处的垂向电离层延迟（VTEC），再通过投影函数获取斜向电离层延迟。但这类方式中电离层改正精度较低，仅能达到5TECU左右。但研究发现，限制其精度的主要因素包括：电离层建模是基于电离层薄层假设建立的所有GNSS卫星的平均电离层延迟模型，无法表征每个单独卫星的准确电离层延迟；而且采用投影函数将斜向电离层延迟（STEC）转换到VTEC进行建模也将存在较大误差，无法有效表征STEC的各向异性。
发明内容
本发明实施例提供一种厘米级斜向电离层延迟产品的生成方法及装置，可解决解决相关技术中电离层产品基于电离层薄层假设以及投影函数所带来的精度低、可用性差的问题。
一方面，本发明实施例提供了一种厘米级斜向电离层延迟产品的生成方法，其特征在于，所述方法包括步骤：
利用非组合PPP技术提取出厘米级斜向电离层延迟信息；
基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品。
一些实施例中，所述利用非组合PPP技术提取出厘米级斜向电离层延迟信息，包括步骤：
收集地面基准站的观测数据并结合实时卫星轨道和钟差产品，利用非组合PPP模型进行解算；
利用FCB产品恢复模糊度的整周特性并对宽项和窄项进行级联固定完成模糊度解算；
将所述级联固定后的宽项和窄项的模糊度约束进所述非组合PPP模型并更新滤波参数后获取所述厘米级斜向电离层延迟信息。
一些实施例中，所述厘米级斜向电离层延迟信息基于第一公式获取，所述第一公式包括：
，
，
其中，为所述厘米级斜向电离层延迟信息，上标s和n分别表示GNSS系统和卫星PRN号，下标m表示基站，/>为第一频点上的纯净的斜向电离层延迟，分别为接收机端和卫星端第一频点和第二频点间的伪距硬件延迟，分别为GNSS第一频点和第二频点的信号频率，/>分别为接收机端和卫星端FCB产品的误差，c为光速。
一些实施例中，所述基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型，包括步骤：
将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式；
将预设时间间隔内基准站m到所有卫星的电离层延迟信息当作虚拟观测值组建误差方程并根据所述误差方程求解二阶多项式系数；
基于所述关于基准站到卫星间高度角和方位角的二阶多项式以及所述二阶多项式系数确定所述斜向电离层延迟模型。
一些实施例中，所述将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式，包括步骤：
基于第二公式生成所述二阶多项式，所述第二公式包括：
，
其中，是t时刻基准站m到卫星n的斜向电离层延迟，/>是待估参数，是t时刻基准站m到卫星n间的高度角和方位角。
一些实施例中，所述斜向电离层延迟模型的残差基于第三公式获取，所述第三公式包括：
，
其中，为所述斜向电离层延迟模型的残差，/>为厘米级斜向电离层延迟信息。
一些实施例中，所述将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品，包括步骤：
挑选区域格网点并将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上；
按照指定格式播发斜向电离层延迟信息产品，并利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟。
一些实施例中，所述将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上，包括步骤：
基于第四公式进行插值，所述第四公式包括：
，
其中,是n卫星在/>格网点出的残差，m是地面基准站个数，/>是/>基准站对n卫星的模型化残差,/>是/>格网点到/>基准站的平面距离。
一些实施例中，所述利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟，包括步骤：
利用播发的m基准站的坐标参数和广播星历计算的t时刻卫星n位置计算出m基准站到n卫星的高度角/>和方位角/>；
将所述和所述/>带入所述第二公式中获取t时刻基准站m到卫星n间的模型化斜向电离层延迟/>再利用反距离加权算法插值出t时刻用户站usr到卫星n间的模型化斜向电离层延迟/>；
基于第五公式获取用户对应的斜向电离层延迟残差，所述第五公式包括：
其中，分别是T1和T2时距离GNSS用户最近的4个格网点处的电离层延迟模型残差改正数，/>是GNSS用户到/>格网点的水平距离；
将所述模型化斜向电离层延迟和所述斜向电离层延迟残差/>相加后获得所述用户到卫星的厘米级斜向电离层延迟。
第二方面，本发明实施例提供了一种厘米级斜向电离层延迟产品的生成装置，其包括：
厘米级斜向电离层信息获取模块，其用于采用非组合PPP技术提取出厘米级斜向电离层延迟信息；
模型构建模块，其用于根据高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
插值与产品生产模块，其用于将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品。
本发明实施例提供了一种厘米级斜向电离层延迟产品的生成方法及装置，将距离适当的基准站获得的电离层斜向延迟（STEC）直接进行插值，避免电离层薄层假设以及投影函数误差带来的精度损失，从而提高电离层延迟信息的精度。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种厘米级斜向电离层延迟产品的生成方法的流程示意图；
图2为本发明实施例提供的一种厘米级斜向电离层延迟产品的生成流程示意图；
图3为本发明实施例提供的格网点处残差的插值示意图；
图4为本发明实施例提供的用户获取高精度斜向电离层延迟示意图；
图5为本发明实施例提供的一种厘米级斜向电离层延迟产品的生成装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供了一种厘米级斜向电离层延迟产品的生成方法，包括步骤：
S100:利用非组合PPP技术提取出厘米级斜向电离层延迟信息；
S200:基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
S300:将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品。
需要说明的是，S100之前，可接入IGS（International GNSS Service）实时轨道和钟差产品以及FCB(Fractional-cycle bias)产品并遍历地面基准站数据。
可以理解的是，本发明实施例提供了一种厘米级斜向电离层延迟产品的生成方法，采用非组合PPP-AR方法提取厘米级的斜向电离层延迟，然后将各基准站到卫星间的斜向电离层延迟信息模型参数化，用户使用时只需利用这些模型化电离层延迟参数恢复出基准站到某颗卫星间的斜向电离层延迟，并结合区域格网点改正数，通过反距离加权平均算法进行插值即可获取用户到这颗卫星的斜向电离层延迟。本发明实施例将距离适当的基准站获得的电离层斜向延迟（STEC）直接进行插值，避免电离层薄层假设以及投影函数误差带来的精度损失，从而提高电离层延迟信息的精度。
一些实施例中，S100包括步骤：
S110:收集地面基准站的观测数据并结合实时卫星轨道和钟差产品，利用非组合PPP模型进行解算；
S120:利用FCB产品恢复模糊度的整周特性并对宽项和窄项进行级联固定完成模糊度解算；
S130:将所述级联固定后的宽项和窄项的模糊度约束进所述非组合PPP模型并更新滤波参数后获取所述厘米级斜向电离层延迟信息。
需要说明的是，S110的解算结果包括浮点解的电离层延迟、模糊度等参数；S130将固定后的模糊度约束进非组合PPP模型中获得固定解的电离层延迟、坐标等信息。
可优选地，考虑到FCB产品存在误差且伪距硬件延迟和电离层延迟无法分离，本发明实施例中厘米级斜向电离层延迟信息基于第一公式获取，且第一公式包括：
，
，
其中，为所述厘米级斜向电离层延迟信息，上标s和n分别表示GNSS系统和卫星PRN号，下标m表示基站，/>为第一频点上的纯净的斜向电离层延迟，分别为接收机端和卫星端第一频点和第二频点间的伪距硬件延迟，分别为GNSS第一频点和第二频点的信号频率，/>分别为接收机端和卫星端FCB产品的误差，c为光速。
一些实施例中，S200包括步骤：
S210:将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式；
S220:将预设时间间隔内基准站m到所有卫星的电离层延迟信息当作虚拟观测值组建误差方程并根据所述误差方程求解二阶多项式系数；
S230:基于所述关于基准站到卫星间高度角和方位角的二阶多项式以及所述二阶多项式系数确定所述斜向电离层延迟模型。
可优选地，S210中基于第二公式生成所述二阶多项式，且第二公式包括：
，
其中，是t时刻基准站m到卫星n的斜向电离层延迟，/>是待估参数，是t时刻基准站m到卫星n间的高度角和方位角。
可优选地，S220中预设时间间隔设为30分钟，那么按照30分钟的时间间隔，将30分钟内基准站到所有卫星的电离层延迟信息当作虚拟观测值组建误差方程，通过最小二乘平差求解二阶多项式系数为：
。
可优选地，考虑到模型化的电离层和真实电离层/>存在差异，S300中斜向电离层延迟模型的残差基于第三公式获取，且第三公式包括：
，
其中，为所述斜向电离层延迟模型的残差，/>为厘米级斜向电离层延迟信息。
如图3和图4所示，一些实施例中，S300包括步骤：
S310:挑选区域格网点并将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上；
S320:按照指定格式播发斜向电离层延迟信息产品，并利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟。
可优选地，S310中基于第四公式进行插值，所述第四公式包括：
，
其中,是n卫星在/>格网点出的残差，m是地面基准站个数，/>是/>基准站对n卫星的模型化残差,/>是/>格网点到/>基准站的平面距离。
需要说明的是，最终播发的斜向电离层延迟信息包括各基准站的多项式系数（基于第二公式）以及格网点残差值。
可以理解的是，将残差提供给用户时，首先挑选区域格网点（推荐以为间隔划分），之后将卫星n的模型化残差/>插值到每个区域格网点上。
本实施例中，将斜向电离层延迟信息产品按照指定形式播发，用户使用时，只需利用播发的产品恢复出基准站到某颗卫星间的斜向电离层延迟，并结合区域格网点改正数，通过反距离加权平均算法进行插值即可获取用户到这颗卫星的斜向电离层延迟。
一个具体的实施例中，按照指定格式播发斜向电离层延迟信息产品时，每30分钟播发一组用于恢复各基准站到所有卫星间斜向电离层延迟信息的多项式系数，每5min播发一组用于恢复模型化残差的格网改正数，具体播发信息如表1所示：
其中，是参与数据处理的基准站点个数，/>是参与数据处理的卫星个数，是所划分格网点的个数。
一些实施例中，S320中利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟，包括步骤：
S321:利用播发的m基准站的坐标参数和广播星历计算的t时刻卫星n位置/>计算出m基准站到n卫星的高度角/>和方位角/>；
S322:将所述和所述/>带入所述第二公式中获取t时刻基准站m到卫星n间的模型化斜向电离层延迟/>再利用反距离加权算法插值出t时刻用户站usr到卫星n间的模型化斜向电离层延迟/>；
S323:基于所述第五公式获取用户对应的斜向电离层延迟残差，所述第五公式包括：
其中，分别是T1和T2时距离GNSS用户最近的4个格网点处的电离层延迟模型残差改正数，/>是GNSS用户到/>格网点的水平距离；
S324:将所述模型化斜向电离层延迟和所述斜向电离层延迟残差/>相加后获得所述用户到卫星的厘米级斜向电离层延迟。
可以理解的是，如图2所示，最终的用户相对于卫星n的厘米级斜向电离层延迟为：/>.....，由于恢复出的斜向电离层延迟包含伪距硬件延迟和FCB误差的影响，所以用户使用本文产品时GNSS模型需采用相同的非组合PPP模型。
如图5所示，本发明实施例还提供一种厘米级斜向电离层延迟产品的生成装置，其包括：
厘米级斜向电离层信息获取模块，其用于采用非组合PPP技术提取出厘米级斜向电离层延迟信息；
模型构建模块，其用于根据高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
插值与产品生产模块，其用于将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品。
一些实施例中，厘米级斜向电离层信息获取模块还用于：
收集地面基准站的观测数据并结合实时卫星轨道和钟差产品，利用非组合PPP模型进行解算；
利用FCB产品恢复模糊度的整周特性并对宽项和窄项进行级联固定完成模糊度解算；
将所述级联固定后的宽项和窄项的模糊度约束进所述非组合PPP模型并更新滤波参数后获取所述厘米级斜向电离层延迟信息。
可优选地，厘米级斜向电离层延迟信息基于第一公式获取，且第一公式包括：
，
，
其中，为所述厘米级斜向电离层延迟信息，上标s和n分别表示GNSS系统和卫星PRN号，下标m表示基站，/>为第一频点上的纯净的斜向电离层延迟，分别为接收机端和卫星端第一频点和第二频点间的伪距硬件延迟，分别为GNSS第一频点和第二频点的信号频率，/>分别为接收机端和卫星端FCB产品的误差，c为光速。
一些实施例中，模型构建模块还用于：
将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式；
将预设时间间隔内基准站m到所有卫星的电离层延迟信息当作虚拟观测值组建误差方程并根据所述误差方程求解二阶多项式系数；
基于所述关于基准站到卫星间高度角和方位角的二阶多项式以及所述二阶多项式系数确定所述斜向电离层延迟模型。
可优选地，模型构建模块基于第二公式生成所述二阶多项式，且第二公式包括：
，
其中，是t时刻基准站m到卫星n的斜向电离层延迟，/>是待估参数，是t时刻基准站m到卫星n间的高度角和方位角。
可优选地，S300中斜向电离层延迟模型的残差基于第三公式获取，且第三公式包括：
，
其中，为所述斜向电离层延迟模型的残差，/>为厘米级斜向电离层延迟信息。
一些实施例中，插值与产品生产模块还用于：
挑选区域格网点并将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上；
按照指定格式播发斜向电离层延迟信息产品，并利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟。
可优选地，插值与产品生产模块基于第四公式进行插值，且第四公式包括：
，
其中,是n卫星在/>格网点出的残差，m是地面基准站个数，/>是/>基准站对n卫星的模型化残差,/>是/>格网点到/>基准站的平面距离。
一些实施例中，插值与产品生产模块还用于：
利用播发的m基准站的坐标参数和广播星历计算的t时刻卫星n位置计算出m基准站到n卫星的高度角/>和方位角/>；
将所述和所述/>带入所述第二公式中获取t时刻基准站m到卫星n间的模型化斜向电离层延迟/>再利用反距离加权算法插值出t时刻用户站usr到卫星n间的模型化斜向电离层延迟/>；
基于所述第五公式获取用户对应的斜向电离层延迟残差，所述第五公式包括：
其中，分别是T1和T2时距离GNSS用户最近的4个格网点处的电离层延迟模型残差改正数，/>是GNSS用户到/>格网点的水平距离；
将所述模型化斜向电离层延迟和所述斜向电离层延迟残差/>相加后获得所述用户到卫星的厘米级斜向电离层延迟。
本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统、装置中的功能模块/单元可以被实施为软件、固件、硬件及其适当的组合。在硬件实施方式中，在以上描述中提及的功能模块/单元之间的划分不一定对应于物理组件的划分；例如，一个物理组件可以具有多个功能，或者一个功能或步骤可以由若干物理组件合作执行。某些物理组件或所有物理组件可以被实施为由处理器，如中央处理器、数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读存储介质上，计算机可读存储介质可以包括计算机可读存储介质（或非暂时性介质）和通信介质（或暂时性介质）。
需要说明的是，在本发明中，诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式，使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所申请的原理和新颖特点相一致的最宽的范围。

Claims (9)

1.一种厘米级斜向电离层延迟产品的生成方法，其特征在于，所述方法包括步骤：
利用非组合PPP技术提取出厘米级斜向电离层延迟信息；
基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品;
所述利用非组合PPP技术提取出厘米级斜向电离层延迟信息，包括步骤：
收集地面基准站的观测数据并结合实时卫星轨道和钟差产品，利用非组合PPP模型进行解算；
利用FCB产品恢复模糊度的整周特性并对宽项和窄项进行级联固定完成模糊度解算；
将所述级联固定后的宽项和窄项的模糊度约束进所述非组合PPP模型并更新滤波参数后获取所述厘米级斜向电离层延迟信息。
2.如权利要求1所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述厘米级斜向电离层延迟信息基于第一公式获取，所述第一公式包括：
，
，
其中，为所述厘米级斜向电离层延迟信息，上标s和n分别表示GNSS系统和卫星PRN号，下标m表示基站，/>为第一频点上的纯净的斜向电离层延迟，/>分别为接收机端和卫星端第一频点和第二频点间的伪距硬件延迟，/>分别为GNSS第一频点和第二频点的信号频率，/>分别为接收机端和卫星端FCB产品的误差，c为光速。
3.如权利要求1所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述基于高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型，包括步骤：
将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式；
将预设时间间隔内基准站m到所有卫星的电离层延迟信息当作虚拟观测值组建误差方程并根据所述误差方程求解二阶多项式系数；
基于所述关于基准站到卫星间高度角和方位角的二阶多项式以及所述二阶多项式系数确定所述斜向电离层延迟模型。
4.如权利要求3所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述将t时刻地面基准站网中各基准站到卫星间的斜向电离层延迟模型参数化为关于基准站到卫星间高度角和方位角的二阶多项式，包括步骤：
基于第二公式生成所述二阶多项式，所述第二公式包括：
，
其中，是t时刻基准站m到卫星n的斜向电离层延迟，/>是待估参数，是t时刻基准站m到卫星n间的高度角和方位角。
5.如权利要求4所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述斜向电离层延迟模型的残差基于第三公式获取，所述第三公式包括：
，
其中，为所述斜向电离层延迟模型的残差，/>为厘米级斜向电离层延迟信息。
6.如权利要求5所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品，包括步骤：
挑选区域格网点并将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上；
按照指定格式播发斜向电离层延迟信息产品，并利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟。
7.如权利要求6所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述将卫星的斜向电离层延迟模型的残差插值到每个区域格网点上，包括步骤：
基于第四公式进行插值，所述第四公式包括：
，
其中,是n卫星在/>格网点出的残差，m是地面基准站个数，/>是/>基准站对n卫星的模型化残差,/>是/>格网点到/>基准站的平面距离。
8.如权利要求7所述的厘米级斜向电离层延迟产品的生成方法,其特征在于，所述利用所述播发的斜向电离层延迟信息产品获取用户到卫星的厘米级斜向电离层延迟，包括步骤：
利用播发的m基准站的坐标参数和广播星历计算的t时刻卫星n位置计算出m基准站到n卫星的高度角/>和方位角/>；
将所述和所述/>带入所述第二公式中获取t时刻基准站m到卫星n间的模型化斜向电离层延迟/>再利用反距离加权算法插值出t时刻用户站usr到卫星n间的模型化斜向电离层延迟/>；
基于第五公式获取用户对应的斜向电离层延迟残差，所述第五公式包括：
其中，分别是T1和T2时距离GNSS用户最近的4个格网点处的电离层延迟模型残差改正数，/>是GNSS用户到/>格网点的水平距离；
将所述模型化斜向电离层延迟和所述斜向电离层延迟残差/>相加后获得所述用户到卫星的厘米级斜向电离层延迟。
9.一种厘米级斜向电离层延迟产品的生成装置，其包括：
厘米级斜向电离层信息获取模块，其用于采用非组合PPP技术提取出厘米级斜向电离层延迟信息；
模型构建模块，其用于根据高度角和方位角的二阶多项式函数构建出基准站到卫星间的斜向电离层延迟模型；
插值与产品生产模块，其用于将所述斜向电离层延迟模型的残差插值到格网点上获取所述厘米级斜向电离层延迟产品;
所述利用非组合PPP技术提取出厘米级斜向电离层延迟信息，包括：
收集地面基准站的观测数据并结合实时卫星轨道和钟差产品，利用非组合PPP模型进行解算；
利用FCB产品恢复模糊度的整周特性并对宽项和窄项进行级联固定完成模糊度解算；
将所述级联固定后的宽项和窄项的模糊度约束进所述非组合PPP模型并更新滤波参数后获取所述厘米级斜向电离层延迟信息。
CN202310892383.7A 2023-07-20 2023-07-20 一种厘米级斜向电离层延迟产品的生成方法及装置 Active CN116609799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310892383.7A CN116609799B (zh) 2023-07-20 2023-07-20 一种厘米级斜向电离层延迟产品的生成方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310892383.7A CN116609799B (zh) 2023-07-20 2023-07-20 一种厘米级斜向电离层延迟产品的生成方法及装置

Publications (2)

Publication Number Publication Date
CN116609799A true CN116609799A (zh) 2023-08-18
CN116609799B CN116609799B (zh) 2023-10-20

Family

ID=87685784

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310892383.7A Active CN116609799B (zh) 2023-07-20 2023-07-20 一种厘米级斜向电离层延迟产品的生成方法及装置

Country Status (1)

Country Link
CN (1) CN116609799B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060017610A1 (en) * 2004-07-23 2006-01-26 California Institute Of Technology Generating high precision ionospheric ground-truth measurements
US20160377730A1 (en) * 2014-12-16 2016-12-29 Trimble Navigation Limited Navigation satellite system positioning involving the generation of correction information
CN106873009A (zh) * 2017-01-12 2017-06-20 中国科学院测量与地球物理研究所 利用非组合ppp辅助的长基线模糊度解算方法
CN108254762A (zh) * 2016-12-28 2018-07-06 千寻位置网络有限公司 伪距差分定位方法及系统
KR101925624B1 (ko) * 2018-01-24 2018-12-05 한국항공대학교산학협력단 지역전리층지도 생성 장치 및 방법
CN112835082A (zh) * 2021-01-05 2021-05-25 广州星际互联科技有限公司 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质
WO2021146775A1 (en) * 2020-01-23 2021-07-29 Ied Foundation Pty Ltd Systems and methods for processing gnss data streams for determination of hardware and atmosphere-delays
US20220171076A1 (en) * 2020-11-30 2022-06-02 Qualcomm Incorporated Global navigation satellite system precise positioning engine with estimated ionosphere
CN116203608A (zh) * 2023-03-09 2023-06-02 广州市中海达测绘仪器有限公司 一种实时精密单点定位的电离层时延补偿方法
CN116203598A (zh) * 2023-02-28 2023-06-02 广州市中海达测绘仪器有限公司 基于地基与星基增强融合的电离层建模方法、装置和介质
CN116299623A (zh) * 2023-05-12 2023-06-23 武汉大学 一种城市复杂场景下的ppp和ins紧组合方法与系统

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060017610A1 (en) * 2004-07-23 2006-01-26 California Institute Of Technology Generating high precision ionospheric ground-truth measurements
US20160377730A1 (en) * 2014-12-16 2016-12-29 Trimble Navigation Limited Navigation satellite system positioning involving the generation of correction information
CN108254762A (zh) * 2016-12-28 2018-07-06 千寻位置网络有限公司 伪距差分定位方法及系统
CN106873009A (zh) * 2017-01-12 2017-06-20 中国科学院测量与地球物理研究所 利用非组合ppp辅助的长基线模糊度解算方法
KR101925624B1 (ko) * 2018-01-24 2018-12-05 한국항공대학교산학협력단 지역전리층지도 생성 장치 및 방법
WO2021146775A1 (en) * 2020-01-23 2021-07-29 Ied Foundation Pty Ltd Systems and methods for processing gnss data streams for determination of hardware and atmosphere-delays
US20220171076A1 (en) * 2020-11-30 2022-06-02 Qualcomm Incorporated Global navigation satellite system precise positioning engine with estimated ionosphere
CN112835082A (zh) * 2021-01-05 2021-05-25 广州星际互联科技有限公司 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质
CN116203598A (zh) * 2023-02-28 2023-06-02 广州市中海达测绘仪器有限公司 基于地基与星基增强融合的电离层建模方法、装置和介质
CN116203608A (zh) * 2023-03-09 2023-06-02 广州市中海达测绘仪器有限公司 一种实时精密单点定位的电离层时延补偿方法
CN116299623A (zh) * 2023-05-12 2023-06-23 武汉大学 一种城市复杂场景下的ppp和ins紧组合方法与系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KE SU 等: "A Novel GNSS Single-Frequency PPP Approach to Estimate the Ionospheric TEC and Satellite Pseudorange Observable-Specific Signal Bias", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol. 60, pages 1 - 12, XP011900624, DOI: 10.1109/TGRS.2021.3126397 *
宋伟伟 等: "不同纬度区域电离层增强PPP-RTK性能分析", 武汉大学学报.信息科学, vol. 46, no. 12, pages 1832 - 1842 *
汪登辉 等: "基于区域CORS的电离层总电子含量提取及其在单频PPP中的应用", 东南大学学报(自然科学版), vol. 43, no. 2, pages 388 - 393 *

Also Published As

Publication number Publication date
CN116609799B (zh) 2023-10-20

Similar Documents

Publication Publication Date Title
CN109764879B (zh) 一种卫星定轨方法、装置及电子设备
JP7054270B2 (ja) 精密単独測位とリアルタイムキネマティックを組み合わせる（ｐｐｐ－ｒｔｋ）の測位方法及び装置
CN111045034B (zh) 基于广播星历的gnss多系统实时精密时间传递方法及系统
CN107861131B (zh) 一种斜路径电离层延迟的获取方法及系统
US20050024263A1 (en) Method for generating clock corrections for a wide-area or global differential GPS system
CN111596321B (zh) 利用非差改正数的多gnss多路径误差恒星日滤波方法及系统
CN110007326B (zh) 一种用于星基增强系统的双频测距误差参数生成方法
CN112327340B (zh) 终端定位精度评估方法、装置、设备以及介质
CN111856534B (zh) 智能终端的双模gnss载波精密单点定位方法及系统
CN114721012B (zh) 一种电离层延迟效应改正方法及装置
CN114966760B (zh) 一种电离层加权的非差非组合ppp-rtk技术实现方法
CN110109158A (zh) 基于gps、glonass和bds多系统的事后超快速rtk定位算法
CN112835082A (zh) 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质
CN112505735B (zh) 对终端进行定位的方法、装置和存储介质
CN108535749B (zh) 基于cors的定位增强方法及系统、定位系统
CN110824505B (zh) Gnss卫星接收机的偏差估计方法及系统、定位方法及终端
CN114879239B (zh) 一种增强瞬时ppp固定解的区域三频整数钟差估计方法
CN116148909B (zh) 多频多模非组合精密单点定位瞬时窄巷模糊度固定方法
CN114114334B (zh) Glonass频间偏差标定方法及rtk方法
CN114928428A (zh) 一种基于增强手段的精密时间传递方法和系统
CN115327593B (zh) 基于无人机的定位方法、系统及存储介质
CN116609799B (zh) 一种厘米级斜向电离层延迟产品的生成方法及装置
CN108254762B (zh) 伪距差分定位方法及系统
CN111123315A (zh) 非差非组合ppp模型的优化方法及装置、定位系统
CN109856656B (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