[发明专利]基于实时网格的网络RTK电离层延迟内插方法及系统

权利要求书

1.一种基于实时网格的网络RTK电离层延迟内插方法，其特征在于，包括以下步骤：

(1)构建Delaunay三角网，选择网络最优独立基线，构成参考站网，计算参考站网实时双差电离层延迟；

(2)计算参考站网的双差电离层延迟对应的所有电离层穿刺点的经纬度；

(3)利用所有电离层穿刺点的经纬度确定实时最优网格；

(4)基于所述实时最优网格，利用实时解算的参考站网双差电离层延迟以及电离层穿刺点经纬度信息，建立双差电离层延迟与周围网格点VTEC的函数关系，然后计算网格点的VTEC，即将双差电离层延迟还原到网格点的VTEC；

(5)选择距流动站最近的参考站为主参考站，根据计算得到的流动站与主参考站间双差电离层延迟对应的电离层穿刺点经纬度，建立流动站与主参考站间双差电离层延迟与网格点VTEC的函数关系，然后计算流动站与主参考站间的双差电离层延迟，即重构流动站与主参考站间的双差电离层延迟。

2.根据权利要求1所述的基于实时网格的网络RTK电离层延迟内插方法，其特征在于，所述步骤(1)中，参考站网实时双差电离层延迟Δ▽I的计算公式如下：

其中，λ₁为载波φ₁的波长，λ₂为载波φ₂的波长，为载波φ₁的双差观测值，为载波φ₂的双差观测值，Δ▽N₁表示φ₁上的双差整周模糊度，Δ▽N₂表示φ₂上的双差整周模糊度。

3.根据权利要求2所述的基于实时网格的网络RTK电离层延迟内插方法，其特征在于，所述步骤(2)中，电离层穿刺点的经纬度的计算过程包括：

已知当前历元的卫星s坐标和测站接收机p坐标，则当前历元电离层穿刺点的大地经纬度l_IPP和B_IPP由下式计算得到：

其中，(l_p，B_p)为测站接收机p的经纬度，A为卫星方位角，为测站接收机p与电离层穿刺点的地心张角；

当前时刻地心张角的计算公式如下：

其中，R为地球半径,H为电离层薄层高度，E为卫星高度角；

当前时刻的卫星高度角E和方位角A的计算公式如下：

其中，(X，Y，Z)为当前历元站心坐标系下的卫星坐标，由下式转换得到：

其中，(ΔX，ΔY，ΔZ)为测站接收机p与卫星s的地心坐标差。

4.根据权利要求1-3任意一项所述的基于实时网格的网络RTK电离层延迟内插方法，其特征在于，所述步骤(3)中，实时最优网格的确定包括以下步骤：

1)根据所有的电离层穿刺点的经纬度划定网格覆盖范围；

2)按照不同的网格大小，划分出所有实时网格，得到一个网格序列；

3)按照优先级将所有的网格排序，得到实时最优网格。