CN120028810B 一种基于gfs和gnss的大高差地区ztd预报方法 （南京邮电大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN120028810B(45)授权公告日2025.07.01

(21)申请号202510503201.1

(22)申请日2025.04.22

(65)同一申请的已公布的文献号申请公布号CN120028810A

(43)申请公布日2025.05.23

(73)专利权人南京邮电大学

地址210003江苏省南京市鼓楼区新模范

马路66号

(72)发明人杨柳王一帆杨薪俊朱丹彤潘志文

(74)专利代理机构南京业腾知识产权代理事务所(特殊普通合伙)32321

专利代理师缪友益

G01S19/37(2010.01)

GO1W1/10(2006.01)

(56)对比文件

CN117669171A,2024.03.08

WO2021094753A1,2021.05.20审查员曾思萌

(51)Int.CI.

GO1S19/03(2010.01)权利要求书2页说明书7页附图1页

(54)发明名称

一种基于GFS和GNSS的大高差地区ZTD预报方法

(57)摘要

CN120028810B本发明提出了一种基于GFS和GNSS的大高差地区ZTD预报方法，首先基于历史GFS预报数据获取GNSS测站位置处的GFS湍流混合延迟和GNSS湍流混合延迟。其次，使用XGBoost方法构建GFS湍流混合延迟和GNSS湍流混合延迟之间的非线性映射关系。最后，使用训练好的XGBoost方法校正预报时刻的GFS数据，实现了大高差地区ZTD的高精度预报。本发明的方法不需要实时GNSSZTD的参与，仅需要时间、位置和GFS预报数据，就能够通过构建好的XGBoost模型获取高精度的

CN120028810B

计算研究区域内所有

计算研究区域内所有GNSS测站位置处的GNSSZTD和GFS分层ZTD的历史数据

基于GFS分层ZTD获取GNSS测站位置处的GNSS湍

流混合延迟和GFS湍流混合延迟

基于XGBoost构建GFS湍流混合延迟到GNSS湍流混合延迟的映射模型

通过预报时刻的GFS垂直分层延迟以及训练好的XGBoost模型预报指定位置及高度的ZTD

CN120028810B权利要求书1/2页

2

1.一种基于GFS和GNSS的大高差地区ZTD预报方法，其特征在于，包括以下步骤：

S1,计算研究区域内所有GNSS测站位置处的GNSSZTD和GFS分层ZTD的历史数据；

S2,基于GFS分层ZTD获取GNSS测站位置处的GNSS湍流混合延迟和GFS湍流混合延迟，具体为：

S21,利用指数函数模型描述ZTD的垂直变化：

ZTD=TD?+ZTDe+ε。

式中，ZTD、TD和ZTD?分别是α位置处的对流层延迟、湍流混合延迟和高度为0m处的垂直分层延迟，e。是α位置处未建模的残差，β是模型系数，ZTDe为垂直分层延迟分量；基于GFS分层ZTD利用指数函数模型，使用最小二乘法估计湍流混合延迟TD。和模型系数β,从而获取GNSS测站位置处GFS垂直分层延迟和GFS湍流混合延迟；

S22,将GNSSZTD与对应位置的处GFS垂直分层延迟作差，得到GNSS湍流混合延迟；

S3,基于XGBoost构建GFS湍流混合延迟到GNSS湍流混合延迟的映射模型，通过预报时刻的GFS垂直分层延迟以及训练好的XGBoost模型预报指定位置及高度的ZTD。

2.根据权利要求1所述的一种基于GFS和GNSS的大高差地区ZTD预报方法，其特征在于，步骤S1具体为：

S11,通过精密单点定位PPP计算所有GNSS测站位置处的GNSSZTD;

S12,计算研究区域内历史GFS预报数据中所有格网点位置处10km以下所有气压层的GFSZTD值；

S13,根据GNSS测站位置最近的四个GFS格网点的GFSZTD值利用反距离加权法计算GNSS测站位置处10km以下所有气压层的GFSZTD值，得到GFS分层ZTD。

3.根据权利要求2所述的一种基于GFS和GNSS的大高差地区ZTD预报方法，其特征在于，S12中GFSZTD值的计算方法如下：

使用Saastamoinen模型计算GFS最高气压层处的ZTDto