[发明专利]一种基于抛物线的区域对流层湿延迟计算方法

权利要求书

1.一种基于抛物线的区域对流层湿延迟计算方法，其特征在于，包括：

S1：采集目标区域中每个探空站点的探空数据、年积日和空间位置数据，根据采集到的探空数据计算每个探空站点的对流层湿延迟真值T_ZWD；所述空间位置数据为纬度值、海拔高度值；

S2：建立对流层湿延迟抛物线函数模型：

其中，Q_ZWD为对流层湿延迟计算值，lat为纬度值，h₀为海拔高度值，doy为年积日，A,B,C和A′,B′,C′均为模型系数；

S3：将S1中获取的探空站点的对流层湿延迟真值、年积日、纬度值、海拔高度值代入对流层湿延迟抛物线函数模型中，确定各项系数，得到目标区域内对流层湿延迟抛物线函数模型；

S4：获取目标区域内待测地的纬度值、海拔高度值和年积日，根据S3确定的目标区域内对流层湿延迟抛物线函数模型得到待测地的对流层湿延迟计算值。

2.根据权利要求1所述的区域对流层湿延迟计算方法，其特征在于，所述步骤S1中计算第n个探空站点的对流层湿延迟真值T_ZWD,n的步骤包括：

S1.1：计算第n个探空站点的对流层总延迟S_ZTD,n：

其中，h_0,n为第n个探空站点的海拔高度，N(h_0,n)第n个探空站点地面处的GNSS无线电信号折射率，N(11000)为距地面11km高度处的折射率，c₁为地面处的折射率衰减系数，c₂是距地面11km高处的折射率衰减系数；

S1.2：计算第n个探空站点的大气延迟A_ZHD,n：

其中P_S,n、f(lat_n,h_0,n)、lat_n分别为第n个探空站点的地面气压、地球自转所引起的重力加速度变化的修正、纬度；

S1.3：第n个探空站点的对流层湿延迟真值T_ZWD,n为：

T_ZWD,n＝S_ZTD,n-A_ZHD,n。

3.根据权利要求1所述的区域对流层湿延迟计算方法，其特征在于，所述步骤S3中采用最小二乘法确定各项系数。

4.根据权利要求1所述的区域对流层湿延迟计算方法，其特征在于，中国区域内纬度为lat，海拔高度为h₀处的对流层湿延迟抛物线函数模型：

其中的系数为：

5.根据权利要求2所述的区域对流层湿延迟计算方法，其特征在于，所述步骤S1.2中，第n个探空站点的地球自转所引起的重力加速度变化的修正为：

f(lat_n,h_0,n)＝1-0.00266cos2lat_n-0.00028h_0,n。