CN111985108A - 一种基于三层球谐近似的gnss电离层tec建模方法 - Google Patents
一种基于三层球谐近似的gnss电离层tec建模方法 Download PDFInfo
- Publication number
- CN111985108A CN111985108A CN202010853007.3A CN202010853007A CN111985108A CN 111985108 A CN111985108 A CN 111985108A CN 202010853007 A CN202010853007 A CN 202010853007A CN 111985108 A CN111985108 A CN 111985108A
- Authority
- CN
- China
- Prior art keywords
- layer
- ionosphere
- tec
- spherical harmonic
- ipp
- 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
Images
Classifications
-
- 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
- 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/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于三层球谐近似的GNSS电离层TEC建模方法,涉及电离层建模技术领域,该方法包括:采用非差非组合PPP方法提取电离层STEC;将所述电离层在垂直方向上分为三个部分,并假设存在三个IPP及对应的三个投影函数,并在每个IPP的位置将斜TEC转换为垂直的天顶方向,获得每个电离层穿刺点IPP处的垂直电子含量VTEC和卫星与接收机传播路径上的电子总含量STEC;对每个电离层穿刺点IPP处的垂直电子含量VTEC通过球谐函数进行建模。本发明通过将电离层在垂直方向上分为三层,并构建三层球谐近似电离层TEC模型,与单层电离层TEC模型相比,本发明的电离层TEC模型精度更高,模型更加稳定。
Description
技术领域
本发明涉及电离层建模技术领域,具体涉及一种基于三层球谐近似的GNSS电离层TEC建模方法。
背景技术
电离层是近地空间环境一个重要组成部分,是地球高空大气的电离区域,位于地面上约60 km到1000 km范围内。电离层的精确建模可改善导航定位用户特别是单频用户的导航定位精度。由于电离层受太阳和地磁活动的影响,空间变化极为复杂,因此对其精确建模仍是一个挑战。单层近似的GNSS (Global Navigation Satellite System)电离层建模最为常见,考虑到电离层F2层的高电子密度,认为电离层薄层的高度在250 km到500 km之间。国际GNSS服务(IGS)各分析中心采用不同的单层近似来描述电离层总电子含量(totalelectron content, TEC)的分布,如CODE采用球谐函数展开,JPL则是在单层电离层球壳上用一个统一的三角格网来描述,ESA采用一个非线性的高斯型电离层指数模型。然而,基于单层近似的电离层TEC建模忽略了垂直领域中的电离层变化,且限制了单层模型的改进范围。因此,如何有效构建更为精确的电离层TEC模型,是一个值得研究的重要科学问题。
近年来,随着GPS、GLONASS、Galileo和北斗在内的多GNSS卫星跟踪站的数量增加,为GNSS电离层TEC建模提供了前所未有的机会。目前,常用的单层近似电离层建模是基于固定高度的单层来简化模型。但是,单层近似建模存在以下问题。其未考虑最大电离层电子密度区域的高度变化和垂直方向上电离层电子密度的不均匀分布。另外,利用投影函数将射线传播方向的TEC转换为天顶方向,且该投影函数仅取决于所观测卫星的高度和假定层的高度,其假定电离层穿刺点(ionospheric pierce point, IPP)周围的电离层电子密度对称分布。然而,在电离层空间分布不均匀的情况下,单层近似可导致1到2个TECu的误差。
发明内容
为了解决上述问题,本发明提供了一种基于三层球谐近似的GNSS电离层TEC建模方法。
本发明提供的具体技术方案为:一种基于三层球谐近似的GNSS电离层TEC建模方法,包括以下步骤:
S1:将所述电离层在垂直方向上分为三层,并采用非差非组合PPP法提取斜总电子含量;
S2:假设每层存在一个IPP及一个对应的投影函数,在每个IPP的位置将倾斜的TEC转换为垂直的天顶方向,获得每个电离层穿刺点IPP处的垂直电子含量VTEC和卫星与接收机传播路径上的电子总含量;
S3:对每个电离层穿刺点IPP处的垂直总电子含量利用球谐函数模型进行拟合,然后将三个球谐函数求和即为三层球谐近似电离层TEC模型。
进一步地,所述方法还包括在提取电离层TEC前对GNSS卫星和接收机的硬件延迟偏差进行校准的过程。
更进一步地,所述投影函数的公式为:
更进一步地,所述卫星与接收机传播路径上的电子总含量STEC通过公式(2)计算得到:
更进一步地,所述球谐函数模型的表达式见式(3):
更进一步地,所述三层球谐近似电离层TEC模型的表达式为(4):
本发明的有益效果:
本发明通过将电离层在垂直方向上分为三层,并构建三层球谐近似电离层TEC模型,与单层电离层TEC模型相比,本发明的电离层TEC模型精度更高,模型更加稳定。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将对本发明作进一步详细的说明。
附图说明:
图1是本发明实施例的电离层三层几何示意图;
图2是本发明实施例的实验区域GNSS测站分布图;
图3是本发明实施例的JFNG测站2015年172日电离层STEC;
图4是本发明实施例的年积日172天UT13:00时刻对应卫星电离层STEC对比;
图5是本发明实施例的年积日172天UT13:00时刻三层和单层电离层模型残差对比。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
参考图1至图3,一种基于三层球谐近似的GNSS电离层TEC建模方法,包括以下步骤:
S1:将所述电离层在垂直方向上分为三层,并采用非差非组合PPP法提取斜总电子含量(slant total electron content, STEC );
在提取电离层TEC时对GNSS卫星和接收机的硬件延迟偏差进行校准的过程,以便在电离层建模过程中获取“干净”的电离层TEC。
提取电离层的区域:经度范围为70°E~140°E和纬度范围为15°N~55°N进行实验分析,选取数据为中国陆态网约250个测站。如附图1的内方框区域所示。
所述将电离层在垂直方向上分为三层,即60 km-300 km,300 km-600 km,600 km-1000 km;假定电离层电子密度集中在每个部分的平均高度即200 km,450 km和800 km。如图1所示,取=200 km, =450 km, =800 km。
S2:假设每层存在一个IPP及一个对应的投影函数,在每个IPP的位置将倾斜的TEC转换为垂直的天顶方向,获得每个电离层穿刺点IPP处的垂直电子含量VTEC和卫星与接收机传播路径上的电子总含量STEC;
并假设存在三个IPP及对应的三个投影函数,所述投影函数的公式为:
在每个IPP的位置将STEC转换为垂直的天顶方向,获得每个电离层穿刺点IPP处的垂直电子含量VTEC和卫星与接收机传播路径上的电子总含量STEC;
所述卫星与接收机传播路径上的电子总含量STEC通过公式(2)计算得到:
S3:对每个电离层穿刺点IPP处的垂直电子含量VTEC利用球谐函数模型进行拟合,然后将三个球谐函数求和即为三层球谐近似电离层TEC模型。
所述球谐函数模型的表达式见式(3):
所述三层球谐近似电离层TEC模型的表达式为(4):
为验证本发明所建立的三层电离层模型的精度,选取IGS中心在中国区域建立的GNSS测站JFNG观测到的电离层总电子含量(STEC)作为本实验的参考值,用来验证三层电离层模型的精度。附图3为实验选取2015年年积日为172天的JFNG测站全天的电离层STEC数据图,在该日共观测到GPS卫星31颗,不同的图案代表不同卫星的电离层STEC。
附图4为单层电离层模型、三层电离层模型及JFNG测站在UT13:00时刻观测卫星的电离层STEC对比图,从图中可以看出,在该时刻视野中存在8颗GPS卫星,每一颗卫星相对于JFNG测站均存在一个电离层STEC,选取JFNG测站的该时刻电离层STEC为参考真值,将单层与多层模型得到的电离层斜延迟与之进行对比,实验表明三层电离层建模得到的电离层STEC精度更高。
附图5为年积日172天UT13:00时刻三层和单层电离层模型残差对比图,从图中可以看出三层电离层模型在对应时刻得到的电离层STEC与参考值的残差控制在2个TECu以内,而对比的单层电离层模型得到的电离层STEC残差则在2-3个TECu范围内,实验结果可以得到三层电离层模型精度更高,模型更加稳定。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、同替换、改进,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于三层球谐近似的GNSS电离层TEC建模方法,其特征在于,包括以下步骤:
S1:将所述电离层在垂直方向上分为三层,并采用非差非组合PPP法提取斜总电子含量;
S2:假设每层存在一个IPP及一个对应的投影函数,在每个IPP的位置将倾斜的TEC转换为垂直的天顶方向,获得每个电离层穿刺点IPP处的垂直电子含量VTEC和卫星与接收机传播路径上的电子总含量;
S3:对每个电离层穿刺点IPP处的垂直总电子含量利用球谐函数模型进行拟合,然后将三个球谐函数求和即为三层球谐近似电离层TEC模型。
2.根据权利要求1所述的基于三层球谐近似的GNSS电离层TEC建模方法,其特征在于,所述方法还包括在提取电离层TEC前对GNSS卫星和接收机的硬件延迟偏差进行校准的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010853007.3A CN111985108A (zh) | 2020-08-22 | 2020-08-22 | 一种基于三层球谐近似的gnss电离层tec建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010853007.3A CN111985108A (zh) | 2020-08-22 | 2020-08-22 | 一种基于三层球谐近似的gnss电离层tec建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111985108A true CN111985108A (zh) | 2020-11-24 |
Family
ID=73443014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010853007.3A Pending CN111985108A (zh) | 2020-08-22 | 2020-08-22 | 一种基于三层球谐近似的gnss电离层tec建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111985108A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140163938A1 (en) * | 2012-02-08 | 2014-06-12 | California Institute Of Technology | Ionospheric slant total electron content analysis using global positioning system based estimation |
CN104933316A (zh) * | 2015-07-02 | 2015-09-23 | 华东交通大学 | 一种基于双参数混合正则化的电离层掩星反演方法 |
CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
CN107085626A (zh) * | 2017-03-17 | 2017-08-22 | 东南大学 | 一种基于bp‑多项式模型融合的区域电离层垂直总电子含量建模方法 |
CN107622304A (zh) * | 2017-08-21 | 2018-01-23 | 东南大学 | 一种基于bp神经网络技术的电离层球谐函数改进方法 |
CN110275185A (zh) * | 2019-07-11 | 2019-09-24 | 武汉大学 | 基于gnss和geo卫星的电离层投影函数建模方法 |
-
2020
- 2020-08-22 CN CN202010853007.3A patent/CN111985108A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140163938A1 (en) * | 2012-02-08 | 2014-06-12 | California Institute Of Technology | Ionospheric slant total electron content analysis using global positioning system based estimation |
CN104933316A (zh) * | 2015-07-02 | 2015-09-23 | 华东交通大学 | 一种基于双参数混合正则化的电离层掩星反演方法 |
CN106405589A (zh) * | 2016-06-24 | 2017-02-15 | 西安科技大学 | 一种全球电离层格网模型的确定方法及装置 |
CN107085626A (zh) * | 2017-03-17 | 2017-08-22 | 东南大学 | 一种基于bp‑多项式模型融合的区域电离层垂直总电子含量建模方法 |
CN107622304A (zh) * | 2017-08-21 | 2018-01-23 | 东南大学 | 一种基于bp神经网络技术的电离层球谐函数改进方法 |
CN110275185A (zh) * | 2019-07-11 | 2019-09-24 | 武汉大学 | 基于gnss和geo卫星的电离层投影函数建模方法 |
Non-Patent Citations (3)
Title |
---|
ZISHEN LI 等: "Regional ionospheric TEC modeling based on a two-layer spherical harmonic approximation for real-time single-frequency PPP", 《JOURNAL OF GEODESY》 * |
刘淑琼 等: "EEMD-RBF神经网络的电离层TEC预报模型", 《测绘工程》 * |
李子申等: "国际GNSS服务组织全球电离层TEC格网精度评估与分析", 《地球物理学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yuan et al. | The BeiDou global broadcast ionospheric delay correction model (BDGIM) and its preliminary performance evaluation results | |
Nie et al. | Quality assessment of CNES real-time ionospheric products | |
Liu et al. | Real-Time Precise Point Positioning (RTPPP) with raw observations and its application in real-time regional ionospheric VTEC modeling | |
Ren et al. | Performance evaluation of real-time global ionospheric maps provided by different IGS analysis centers | |
Liu et al. | The cooperative IGS RT-GIMs: A reliable estimation of the global ionospheric electron content distribution in real time | |
Ren et al. | Global ionospheric modeling using multi-GNSS and upcoming LEO constellations: Two methods and comparison | |
Brunini et al. | A new ionosphere monitoring technology based on GPS | |
CN114690207A (zh) | 差分电离层建模方法及系统 | |
Li et al. | FY-3D and FY-3C onboard observations for differential code biases estimation | |
Chen et al. | Combination of ground-and space-based data to establish a global ionospheric grid model | |
Hoque et al. | A new climatological electron density model for supporting space weather services | |
Zus et al. | Ionospheric mapping functions based on electron density fields | |
CN110146904B (zh) | 一种适用于区域电离层tec的精确建模方法 | |
Chen et al. | A novel ionospheric mapping function modeling at regional scale using empirical orthogonal functions and GNSS data | |
She et al. | Estimation of ionospheric total electron content from a multi-GNSS station in China | |
Zhang et al. | Performance evaluation of BDS-3 ionospheric delay correction models (BDSK and BDGIM): First year for full operational capability of global service | |
Chen et al. | Assessment of BDS-3 terrestrial reference frame realized by broadcast ephemeris: comparison with GPS/Galileo | |
Gao et al. | Positioning via GEO communication satellites’ signals of opportunity | |
Chen et al. | Near real-time global ionospheric modeling based on an adaptive Kalman filter state error covariance matrix determination method | |
Bahadur et al. | Real-time single-frequency multi-GNSS positioning with ultra-rapid products | |
CN117269987A (zh) | 一种低轨增强sbas电离层监测系统及格网估计方法 | |
Angrisano et al. | Ionospheric models comparison for single-frequency GNSS positioning | |
Cheng et al. | Global monitoring of geomagnetic storm-induced ionosphere anomalies using 3-D ionospheric modeling with multi-GNSS and COSMIC measurements | |
Klimenko et al. | Longitudinal variation in the ionosphere-plasmasphere system at the minimum of solar and geomagnetic activity: Investigation of temporal and latitudinal dependences | |
CN111985108A (zh) | 一种基于三层球谐近似的gnss电离层tec建模方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201124 |