一种GNSS水汽站气压与温度无仪器获取方法及系统与流程

本发明涉及gnss气象学应用领域，尤其涉及一种gnss水汽站气压与温度无仪器获取方法和一种gnss水汽站气压与温度无仪器获取系统。

背景技术：

水汽分布及其变化是灾害性天气形成和演变过程中的重要影响因子，是气象学和天气预报研究的重要问题。gnss解算水汽技术以其高时空分辨率、高精度、全天候、近实时等优点成为新一代大气遥感技术中最有发展前景的方法之一，被广泛应用于暴雨灾害的预警和气候变化的监测。

在地基gps反演水汽的实际应用中，站点温度和气压的精确测定至关重要。精确确定站点气压，可以依据干延迟模型计算准确的干延迟，有利于从天顶总延迟中剥离出准确的天顶湿延迟；精确确定温度有利于计算加权平均温度模型，提高天顶湿延迟到水汽可降水量pwv的换算精度。因我国地域辽阔、经济和社会发展的差异较大，并不是每个gps台站附近都安装有气象要素观测仪。如何精确确定这些台站附近的气压和温度，成为地基gps气象站网布设的一个重要的问题。

技术实现要素：

本发明旨在至少解决现有技术或相关技术中存在的技术问题之一。

为此，本发明的一个目的在于提供一种gnss水汽站气压与温度无仪器获取方法，其能够为gnss水汽站提供高精度的气压与温度，尤其是在未布设温度和气压传感器的gnss气象站网中，提高gnss水汽解算中gnss水汽站的利用率。

本发明的另一个目的在于提供一种gnss水汽站气压与温度无仪器获取系统，其能够为gnss水汽站提供高精度的气压与温度，尤其是在未布设温度和气压传感器的gnss气象站网中，提高gnss水汽解算中gnss水汽站的利用率。

为实现上述目的，本发明第一方面的技术方案提供了一种gnss水汽站气压与温度无仪器获取方法，包括以下步骤：获取待求gnss水汽站的三维大地坐标；根据预设条件检索对应待求gnss水汽站的三维大地坐标的其他gnss水汽站的周围气温数据和周围气压数据；根据待求gnss水汽站与其他gnss水汽站的海拔高度差对周围气温数据进行线性高度订正和对周围气压数据进行非线性高度订正，得到气温订正值和气压订正值；计算待求gnss水汽站与其他gnss水汽站的距离，并根据距离求出待求gnss水汽站与其他gnss水汽站的权函数；根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据。

在该技术方案中，该方法主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要其周边区域内有气温传感器和气压传感器的gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

在上述技术方案中，优选地，预设条件包括其他gnss水汽站的数量为4～6个，其他gnss水汽站与待求gnss水汽站的最大距离小于或等于300km。

在该技术方案中，该方法主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要位于待求gnss水汽站周边300km以内有气温传感器和气压传感器的4-6个其他gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

在上述任一技术方案中，优选地，气温订正值的表达式为：δt＝δh×0.0065，δt为气温订正值，δh为海拔高度差；气压订正值的表达式为：δp为气压订正值，p1为周围气压数据，ti为周围气温数据，δh为海拔高度差。

在上述任一技术方案中，优选地，计算待求gnss水汽站与其他gnss水汽站的距离的表达式为：x＝(n+h)*cosb*cosl；y＝(n+h)*cosb*sinl；δs为距离；n为椭球面卯酉圈的曲率半径，e为椭球面的第一偏心率，(b，l，h)为其他gnss气象站的三维大地坐标；权函数的表达式为：其中，ωi为其他gnss水汽站中第i个gnss气象站的权函数，μ的值为2，n为gnss水汽站的总数。

在上述任一技术方案中，优选地，根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据的表达式为：

其中，p、t为待求gnss水汽站的气压数据和气温数据，δti′和δpi′分别为气温订正值和气压订正值，p′和ti′分别为第i个其他gnss气水汽站的周围气压数据和周围气温数据。

本发明第二方面的技术方案提供了一种gnss水汽站气压与温度无仪器获取系统，包括：获取模块，被设置为用于获取待求gnss水汽站的三维大地坐标；检索模块，被设置为用于根据预设条件检索对应待求gnss水汽站的三维大地坐标的其他gnss水汽站的周围气温数据和周围气压数据；订正模块，被设置为用于根据待求gnss水汽站与其他gnss水汽站的海拔高度差对周围气温数据进行线性高度订正和对周围气压数据进行非线性高度订正，得到气温订正值和气压订正值；权函数构建模块，被设置为用于计算待求gnss水汽站与其他gnss水汽站的距离，并根据距离求出待求gnss水汽站与其他gnss水汽站的权函数；计算模块，被设置为用于根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据。

在该技术方案中，该系统主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要其周边区域内有气温传感器和气压传感器的gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

在上述技术方案中，优选地，预设条件包括其他gnss水汽站的数量为4～6个，其他gnss水汽站与待求gnss水汽站的最大距离小于或等于300km。

在该技术方案中，该系统主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要位于待求gnss水汽站周边300km以内有气温传感器和气压传感器的4-6个其他gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

在上述任一技术方案中，优选地，气温订正值的表达式为：δt＝δh×0.0065，δt为气温订正值，δh为海拔高度差；气压订正值的表达式为：δp为气压订正值，p1为周围气压数据，ti为周围气温数据，δh为海拔高度差。

在上述任一技术方案中，优选地，计算待求gnss水汽站与其他gnss水汽站的距离的表达式为：x＝(n+h)*cosb*cosl；y＝(n+h)*cosb*sinl；δs为距离；n为椭球面卯酉圈的曲率半径，e为椭球面的第一偏心率，(b，l，h)为其他gnss气象站的三维大地坐标；权函数的表达式为：其中，ωi为其他gnss水汽站中第i个gnss气象站的权函数，μ的值为2，n为gnss水汽站的总数。

在上述任一技术方案中，优选地，根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据的表达式为：

其中，p、t为待求gnss水汽站的气压数据和气温数据，δt′i和δp′i分别为气温订正值和气压订正值，p′和t′i分别为第i个其他gnss气水汽站的周围气压数据和周围气温数据。

附图说明

本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1示出了本发明实施例所涉及获取方法的流程框图；

图2示出了本发明实施例所涉及获取系统的结构框图；

图3示出了本发明具体实施例所涉及具体计算步骤的原理图；

图4示出了本发明具体实施例所涉及气压残差图；

图5示出了本发明具体实施例所涉及气温残差图。

具体实施方式

为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。

在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不限于下面公开的具体实施例的限制。

下面参照图1至图5描述根据本发明一些实施例的gnss水汽站气压与温度无仪器获取方法及系统。

如图1所示，按照本发明一个实施例的gnss水汽站气压与温度无仪器获取方法，包括以下步骤：

s100，获取待求gnss水汽站的三维大地坐标；

s200，根据预设条件检索对应待求gnss水汽站的三维大地坐标的其他gnss水汽站的周围气温数据和周围气压数据；

s300，根据待求gnss水汽站与其他gnss水汽站的海拔高度差对周围气温数据进行线性高度订正和对周围气压数据进行非线性高度订正，得到气温订正值和气压订正值；

s400，计算待求gnss水汽站与其他gnss水汽站的距离，并根据距离求出待求gnss水汽站与其他gnss水汽站的权函数；

s500，根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据。

在该实施例中，该方法主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要其周边区域内有气温传感器和气压传感器的gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

如图2所示，按照本发明另一个实施例的gnss水汽站气压与温度无仪器获取系统1000，包括：

获取模块100，被设置为用于获取待求gnss水汽站的三维大地坐标；

检索模块200，被设置为用于根据预设条件检索对应待求gnss水汽站的三维大地坐标的其他gnss水汽站的周围气温数据和周围气压数据；

订正模块300，被设置为用于根据待求gnss水汽站与其他gnss水汽站的海拔高度差对周围气温数据进行线性高度订正和对周围气压数据进行非线性高度订正，得到气温订正值和气压订正值；

权函数构建模块400，被设置为用于计算待求gnss水汽站与其他gnss水汽站的距离，并根据距离求出待求gnss水汽站与其他gnss水汽站的权函数；

计算模块500，被设置为用于根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据。

在该实施例中，该系统主要针对未安装气温传感器和气压传感器的gnss水汽站，仅需要其周边区域内有气温传感器和气压传感器的gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

具体而言，按照本发明上述任一实施例所涉及的gnss水汽站气压与温度无仪器获取方法或系统，预设条件包括其他gnss水汽站的数量为4～6个，其他gnss水汽站与待求gnss水汽站的最大距离小于或等于300km。

在该实施例中，针对未安装气温传感器和气压传感器的gnss水汽站，仅需要位于待求gnss水汽站周边300km以内有气温传感器和气压传感器的4-6个其他gnss水汽站，就能求得gnss水汽反演所必需的气压和温度，使用简单，结果精度高，对于提高gnss水汽反演站网的密度具有重要意义。

在上述任一实施例中，优选地：

气压订正值的表达式为：

δp为气压订正值，p1为周围气压数据，ti为周围气温数据，δh为海拔高度差；

气温订正值的表达式为：

δt＝δh×0.0065(2)；

δt为气温订正值，δh为海拔高度差。

在上述任一实施例中，优选地：

计算待求gnss水汽站与其他gnss水汽站的距离的表达式为：

δs为距离；n为椭球面卯酉圈的曲率半径，e为椭球面的第一偏心率，(b，l，h)为其他gnss气象站的三维大地坐标；

权函数的表达式为：

其中，ωi为其他gnss水汽站中第i个gnss气象站的权函数，μ的值为2，n为gnss水汽站的总数。

在上述任一实施例中，优选地：

根据权函数、气温订正值和气压订正值计算待求gnss水汽站的气压数据和气温数据的表达式为：

其中，p、t为待求gnss水汽站的气压数据和气温数据，δt′i和δp′i分别为气温订正值和气压订正值，p′和t′i分别为第i个其他gnss气水汽站的周围气压数据和周围气温数据。

具体实施例

以某gnss水汽站为例，具体计算步骤如下：

步骤1，参见图3，a为待求gnss水汽站，其大地坐标为(ba，la，ha)，b1～b5为a周围的5个有气象数据(气温和气压)的其他gnss水汽站，其大地坐标分别为(b1，l1，h1)～(b5，l5，h5)(因坐标保密需要，真实值不列出)。b1～b5五个站点某一时刻的气压和温度分别为p1～p5，t1～t5。

步骤2，采用上式(1)和(2)，对b1～b5五个站点的气压和温度进行高度订正，其中δh为ha与h1～h5的差值，求出b1～b5五个站点的气压和温度的气压订正值和气温订正值，分别为△p1～△p5，△t1～△t5。

步骤3，采用上式(3)，将gnss水汽站的大地坐标转换成wgs84下的空间直角坐标，并求出b1～b5与a的空间距离，分别为s1～s5。

采用上式(4)，求出b1～b5五个站点的权函数，分别为ω1～ω5。

步骤4，采用上式(5)，求出待求站点a的气压和温度：

pa＝ω1(p1+△p1)+ω2(p2+△p2)+ω3(p3+△p3)+ω4(p4+△p4)+ω5(p5+△p5)；

ta＝ω1(t1+△t1)+ω2(t2+△t2)+ω3(t3+△t3)+ω4(t4+△t4)+ω5(t5+△t5)。

以该gnss水汽站90天数据为例，利用本发明方法求得的气压和温度，其均方根误差分别为0.68hpa、0.74k。气压和温度的残差参见附图4和附图5。上述均方根误差值和附图，都说明本发明提出的一种gnss水汽站气压与温度无仪器获取方法和系统具有精度高的优点。

在本发明中，术语“第一”、“第二”、“第三”仅用于描述的目的，而不能理解为指示或暗示相对重要性；术语“多个”则指两个或两个以上，除非另有明确的限定。术语“安装”、“相连”、“连接”、“固定”等术语均应做广义理解，例如，“连接”可以是固定连接，也可以是可拆卸连接，或一体地连接；“相连”可以是直接相连，也可以通过中间媒介间接相连。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。

本发明的描述中，需要理解的是，术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或单元必须具有特定的方向、以特定的方位构造和操作，因此，不能理解为对本发明的限制。

在本说明书的描述中，术语“一个实施例”、“一些实施例”、“具体实施例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或实例。而且，描述的具体特征、结构、材料或特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。