CN106407560A - 表征大气各向异性的对流层映射函数模型的构建方法 - Google Patents
表征大气各向异性的对流层映射函数模型的构建方法 Download PDFInfo
- Publication number
- CN106407560A CN106407560A CN201610831005.8A CN201610831005A CN106407560A CN 106407560 A CN106407560 A CN 106407560A CN 201610831005 A CN201610831005 A CN 201610831005A CN 106407560 A CN106407560 A CN 106407560A
- Authority
- CN
- China
- Prior art keywords
- troposphere
- mapping function
- function model
- delay
- azimuth
- 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
Links
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
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种表征大气各向异性的对流层映射函数模型的构建方法,包括步骤:S1收集测站周围空间的气象数据;S2根据气象数据,利用射线追踪法计算测站的对流层天顶延迟及不同方位角和不同高度角下的对流层斜延迟;S3分别计算各对流层斜延迟与对流层天顶延迟的比值;S4以比值为观测值,构建误差方程;S5采用最小二乘原理求解对流层映射函数模型系数的最优解。本发明将卫星信号的高度角和方位角同时作为自变量融入对流层映射函数模型中,能完整、准确地表达对流层延迟在不同的高度角和方位角上的差异,将所构建对流层映射函数模型用于卫星导航定位、水汽探测等,可进一步提高卫星导航定位、水汽探测的精度。
Description
技术领域
本发明属于卫星导航、水汽遥感及无线电传播技术领域,尤其涉及一种表征大气各向异性的对流层映射函数模型的构建方法。
背景技术
对流层延迟是卫星导航信号传播过程中的一个重要误差源,其数值在天顶方向约2.3m,到地平方向可达30m。对流层延迟在射电波段具有非色散性,无法通过多个频率信号的组合予以消除。通常将任意高度角方向的对流层斜延迟表达为天顶延迟与映射函数的乘积。现有映射函数模型均基于大气层球对称这一假设,绝大部分采用三阶连分式的表达形式,将映射函数模型表达为高度角的函数,见式(1):
式(1)中,m(ε)表示以高度角ε为自变量的映射函数模型;ai、bi、ci均为映射函数模型的系数,i=h,nh,h和nh分别表示流体静力学分量和非流体静力学分量。
实际上,大气状态在不同方位上是有差异的,极端天气下尤其明显,因此其映射函数也应表达为方位角的函数。然而三阶连分式无法将方位角纳入其中。为此,现有方法采用附加梯度参数估计的方法进行补偿,然而梯度只考虑了南-北、东-西方向的两个参数,不足以弥补现有映射函数的缺陷。因而亟需一种新的映射函数模型,能够完整、准确地表达对流层延迟在各个高度角和方位角上的差异,从而进一步提高卫星导航定位、水汽探测的精度。
发明内容
针对现有技术存在的问题,本发明将卫星信号的高度角和方位角同时作为自变量纳入对流层映射函数模型,提供了一种表征大气各向异性的对流层映射函数模型的构建方法。
为解决上述技术问题,本发明采用如下技术方案:
一种表征大气各向异性的对流层映射函数模型的构建方法,包括步骤:
S1收集测站周围空间的气象数据;
S2根据气象数据,利用射线追踪法计算测站的对流层天顶延迟ZD及不同方位角和不同高度角下的对流层斜延迟SD(φ,θ),SD(φ,θ)即方位角φ和高度角θ下的对流层斜延迟;
S3分别计算各对流层斜延迟SD(φ,θ)与对流层天顶延迟ZD的比值L(φ,θ);
S4以比值L(φ,θ)为观测值,构建误差方程V=MF(φ,θ)-L(φ,θ),其中,V表示观测值L(φ,θ)的改正数;MF(φ,θ)表示对流层映射函数模型;
所述的对流层映射函数模型其中:
和为系数;
N为满足的最小整数,α为模型分辨率,α根据实际的精度要求取值;
MF(φ,θ)表示以方位角φ和高度角θ为自变量的映射函数模型;
r为满足条件的最大整数;
S5采用最小二乘原理求解对流层映射函数模型系数和的最优解。
上述气象数据来源于数值天气模型。
上述对流层天顶延迟ZD和对流层斜延迟SD(φ,θ)为流体静力学延迟、非流体静力学延迟或总延迟。
步骤S2中,以为步长,在方位角和高度角的范围取一系列不同的方位角和高度角,然后计算各不同方位角和不同高度角下的对流层斜延迟SD(φ,θ);其中,方位角的范围为0°~360°,高度角的范围为0°~90°,α为模型分辨率,q为满足q≥2的整数。
和现有技术相比,本发明具有如下优点和有益效果:
(1)将卫星信号的高度角和方位角同时作为自变量融入对流层映射函数模型中,有效克服了传统对流层映射函数模型精度不够的缺陷,为解决大气各向异性问题提供了新的途径。
(2)能完整、准确地表达对流层延迟在不同的高度角和方位角上的差异,将所构建对流层映射函数模型用于卫星导航定位、水汽探测等,可进一步提高卫星导航定位、水汽探测的精度。
附图说明
图1为本发明所提出对流层映射函数模型的原理示意图。
具体实施方式
为便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
本发明提供了一种表征大气各向异性的对流层映射函数模型,其原理见图1,其数学形式见公式(1):
式(1)中:
MF(φ,θ)表示以方位角φ和高度角θ为自变量的映射函数模型;
和为系数,其值待定;
N为满足的最小整数,α为模型分辨率,其单位为度,模型分辨率α根据实际需求进行取值,当对模型精度要求越高,模型分辨率α取值越小,模型分辨率α一般在0°~90°范围内取值;
r为满足条件的最大整数。
本发明提供的表征大气各向异性的对流层映射函数模型的构建方法,具体步骤如下:
第一步,收集测站周围空间的气象数据,所述的气象数据包括气温、气压和相对湿度。具体实施中,气象数据可来源于数值天气模型。
第二步,根据气象数据,利用射线追踪法计算测站的对流层天顶延迟ZD以及不同方位角和不同高度角下的对流层斜延迟SD(φ,θ),SD(φ,θ)即方位角φ和高度角θ下的对流层斜延迟。
方位角的取值范围为0°~360°,高度角的取值范围为0°~90°,以为步长,在方位角和高度角的范围取一系列的方位角和高度角。α为模型分辨率,q为满足q≥2的整数。
本发明中,对流层天顶延迟ZD和对流层斜延迟SD(φ,θ)可以为流体静力学延迟(即干延迟)、非流体静力学延迟(即湿延迟)或总延迟。
第二步,计算各对流层斜延迟SD(φ,θ)与对流层天顶延迟的比值L(φ,θ):
第三步,以比值L(φ,θ)为观测值,构建误差方程V=MF(φ,θ)-L(φ,θ),其中,V表示观测值L(φ,θ)的改正数,MF(φ,θ)为以方位角φ和高度角θ为自变量的映射函数模型,见公式(1)。
第四步,采用最小二乘原理求解映射函数模型系数和的最优解。
在已知对流层天顶延迟的情况下,采用本发明所构建的对流层映射函数模型可准确计算出卫星信号在任意高度和方位上的对流层延迟误差。在对流层天顶延迟未知或仅有模型初始值的情况下,采用本发明所构建的对流层映射函数模型,联立多个历元的GNSS相位观测方程进行参数估计可精确估算出对流层天顶延迟。
Claims (4)
1.表征大气各向异性的对流层映射函数模型的构建方法,其特征是,包括步骤:
S1收集测站周围空间的气象数据;
S2根据气象数据,利用射线追踪法计算测站的对流层天顶延迟ZD及不同方位角和不同高度角下的对流层斜延迟SD(φ,θ),SD(φ,θ)即方位角φ和高度角θ下的对流层斜延迟;
S3分别计算各对流层斜延迟SD(φ,θ)与对流层天顶延迟ZD的比值L(φ,θ);
S4以比值L(φ,θ)为观测值,构建误差方程V=MF(φ,θ)-L(φ,θ),其中,V表示观测值L(φ,θ)的改正数;MF(φ,θ)表示对流层映射函数模型;
所述的对流层映射函数模型其中:
和为系数;
N为满足的最小整数,α为模型分辨率,α根据实际的精度要求取值;
MF(φ,θ)表示以方位角φ和高度角θ为自变量的映射函数模型;
r为满足条件的最大整数;
S5采用最小二乘原理求解对流层映射函数模型系数和的最优解。
2.如权利要求1所述的表征大气各向异性的对流层映射函数模型的构建方法,其特征是:
所述的气象数据来源于数值天气模型。
3.如权利要求1所述的表征大气各向异性的对流层映射函数模型的构建方法,其特征是:
所述的对流层天顶延迟ZD和对流层斜延迟SD(φ,θ)为流体静力学延迟、非流体静力学延迟或总延迟。
4.如权利要求1所述的表征大气各向异性的对流层映射函数模型的构建方法,其特征是:
步骤S2中,以为步长,在方位角和高度角的范围取一系列不同的方位角和高度角,然后计算各不同方位角和不同高度角下的对流层斜延迟SD(φ,θ);其中,方位角的范围为0°~360°,高度角的范围为0°~90°,α为模型分辨率,q为满足q≥2的整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610831005.8A CN106407560B (zh) | 2016-09-19 | 2016-09-19 | 表征大气各向异性的对流层映射函数模型的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610831005.8A CN106407560B (zh) | 2016-09-19 | 2016-09-19 | 表征大气各向异性的对流层映射函数模型的构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106407560A true CN106407560A (zh) | 2017-02-15 |
CN106407560B CN106407560B (zh) | 2019-03-19 |
Family
ID=57996987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610831005.8A Expired - Fee Related CN106407560B (zh) | 2016-09-19 | 2016-09-19 | 表征大气各向异性的对流层映射函数模型的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106407560B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031877A (zh) * | 2019-04-17 | 2019-07-19 | 山东科技大学 | 一种基于grnn模型的区域nwp对流层延迟改正方法 |
CN111273320A (zh) * | 2020-02-27 | 2020-06-12 | 东南大学 | 一种顾及对流层残余延迟的gnss随机模型建立方法 |
CN112130173A (zh) * | 2020-09-29 | 2020-12-25 | 长安大学 | 一种用于gnss反射高程测量的大气改正方法 |
CN112949294A (zh) * | 2021-02-05 | 2021-06-11 | 国家基础地理信息中心 | 湿延迟数据文本的生成方法、装置、设备及存储介质 |
CN113093241A (zh) * | 2021-03-12 | 2021-07-09 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113093242A (zh) * | 2021-03-17 | 2021-07-09 | 山东科技大学 | 一种基于球谐展开的gnss单点定位方法 |
CN114910939A (zh) * | 2022-07-13 | 2022-08-16 | 武汉大学 | 短距离大高差rtk中对流层延迟实测气象改正方法 |
CN116068670B (zh) * | 2023-03-30 | 2023-06-06 | 中国科学院精密测量科学与技术创新研究院 | 一种无地面网络区适用的北斗水汽场实时重构方法及装置 |
CN116361714A (zh) * | 2023-06-01 | 2023-06-30 | 山东科技大学 | 一种顾及非各向同性的水平方向对流层延迟分类方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101802643A (zh) * | 2007-09-18 | 2010-08-11 | 莱卡地球系统公开股份有限公司 | 对卫星辅助定位系统的基于网络的校正进行精确度估计的方法 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
-
2016
- 2016-09-19 CN CN201610831005.8A patent/CN106407560B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101802643A (zh) * | 2007-09-18 | 2010-08-11 | 莱卡地球系统公开股份有限公司 | 对卫星辅助定位系统的基于网络的校正进行精确度估计的方法 |
CN102928850A (zh) * | 2012-11-22 | 2013-02-13 | 桂林电子科技大学 | 一种广域电离层误差改正新方法 |
Non-Patent Citations (3)
Title |
---|
DI ZHANG 等: "Quantitative assessment of meteorological and tropospheric Zenith Hydrostatic Delay models", 《ADVANCES IN SPACE RESEARCH》 * |
郭际明 等: "利用射线追踪法分析三种典型对流层映射函数在中国区域的精度", 《武汉大学学报(信息科学版)》 * |
高兴国 等: "GNSS对流层延时映射函数影响分析比较研究", 《武汉大学学报(信息科技版)》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031877B (zh) * | 2019-04-17 | 2020-05-26 | 山东科技大学 | 一种基于grnn模型的区域nwp对流层延迟改正方法 |
CN110031877A (zh) * | 2019-04-17 | 2019-07-19 | 山东科技大学 | 一种基于grnn模型的区域nwp对流层延迟改正方法 |
CN111273320A (zh) * | 2020-02-27 | 2020-06-12 | 东南大学 | 一种顾及对流层残余延迟的gnss随机模型建立方法 |
CN111273320B (zh) * | 2020-02-27 | 2020-10-02 | 东南大学 | 一种顾及对流层残余延迟的gnss随机模型建立方法 |
CN112130173A (zh) * | 2020-09-29 | 2020-12-25 | 长安大学 | 一种用于gnss反射高程测量的大气改正方法 |
CN112949294B (zh) * | 2021-02-05 | 2022-09-30 | 国家基础地理信息中心 | 湿延迟数据文本的生成方法、装置、设备及存储介质 |
CN112949294A (zh) * | 2021-02-05 | 2021-06-11 | 国家基础地理信息中心 | 湿延迟数据文本的生成方法、装置、设备及存储介质 |
CN113093241A (zh) * | 2021-03-12 | 2021-07-09 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113093241B (zh) * | 2021-03-12 | 2024-02-02 | 东南大学 | 一种考虑高度角的单测站对流层斜延迟计算方法 |
CN113093242B (zh) * | 2021-03-17 | 2022-03-11 | 山东科技大学 | 一种基于球谐展开的gnss单点定位方法 |
CN113093242A (zh) * | 2021-03-17 | 2021-07-09 | 山东科技大学 | 一种基于球谐展开的gnss单点定位方法 |
CN114910939A (zh) * | 2022-07-13 | 2022-08-16 | 武汉大学 | 短距离大高差rtk中对流层延迟实测气象改正方法 |
CN114910939B (zh) * | 2022-07-13 | 2022-11-04 | 武汉大学 | 短距离大高差rtk中对流层延迟实测气象改正方法 |
CN116068670B (zh) * | 2023-03-30 | 2023-06-06 | 中国科学院精密测量科学与技术创新研究院 | 一种无地面网络区适用的北斗水汽场实时重构方法及装置 |
CN116361714A (zh) * | 2023-06-01 | 2023-06-30 | 山东科技大学 | 一种顾及非各向同性的水平方向对流层延迟分类方法 |
CN116361714B (zh) * | 2023-06-01 | 2023-08-04 | 山东科技大学 | 一种顾及非各向同性的水平方向对流层延迟分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106407560B (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106407560B (zh) | 表征大气各向异性的对流层映射函数模型的构建方法 | |
WO2022048694A1 (zh) | 一种基于球谐展开的gnss单点定位方法 | |
CN105629263B (zh) | 一种对流层大气延迟误差估计改正方法和改正系统 | |
Alber et al. | Obtaining single path phase delays from GPS double differences | |
CN102034238B (zh) | 基于光学成像测头和视觉图结构的多摄像机系统标定方法 | |
CN110487301A (zh) | 一种雷达辅助机载捷联惯性导航系统初始对准方法 | |
CN104111467B (zh) | 一种基于北斗三频宽巷组合的网络rtk瞬时定位方法 | |
CN104748722B (zh) | 利用卫星定位信息实时校准气压测高结果的高程定位方法 | |
CN109752696B (zh) | 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法 | |
CN104457684A (zh) | 免固定设站点的全站仪三维变形监测方法 | |
CN102278970B (zh) | 一种基于全站仪角距差分的定位变形监测技术 | |
CN104535080B (zh) | 大方位失准角下基于误差四元数的传递对准方法 | |
CN109959898B (zh) | 一种座底式水声被动定位基阵自校准方法 | |
CN109597040A (zh) | 一种星载sar影像无场几何定标方法 | |
CN110045340A (zh) | 平面相控阵天线中信标天线的坐标测算方法 | |
CN108919304A (zh) | 一种基于参考平面的移动测量系统中pos误差补偿方法 | |
CN105572639B (zh) | 一种室内超声波差分定位方法 | |
CN107316280A (zh) | 离岛卫星影像rpc模型高精度几何定位方法 | |
US20240159529A1 (en) | Systems and methods for extending the spatial coverage of a reference pressure network | |
CN102207380A (zh) | 一种高精度横轴倾斜误差补偿方法 | |
CN115597586A (zh) | 一种基于大气偏振模式对称性的航向角提取方法 | |
CN102155956B (zh) | 一种高精度垂直角横轴倾斜误差补偿方法 | |
CN104994580B (zh) | 一种室内定位方法 | |
CN102426352A (zh) | 风廓线雷达风计算方法 | |
CN107037463A (zh) | 一种地基长波授时系统与gnss系统的组合定位授时方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190319 Termination date: 20190919 |