1.一种基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,包括:步骤1,获取卫星星历、基准站原始观测值及基准站坐标数据;步骤2,根据基准站坐标,划定长方形区域,计算出长方形区域左下角最小经纬度值BL min 和右上角最大经纬度值BL max ,并按间隔inval依次划分出格网点;步骤3,将区域内基准站点按照Delaunay三角网元形式进行剖分,划分为一个个三角网元;步骤4,计算其中一个三角网元内基准站间模糊度,先进行宽巷模糊度 固定,再利用无电离层组合模型固定宽巷及L1模糊度 恢复其整数特性,最后利用已固定的宽巷及L1模糊度确定L2模糊度 大小;步骤5,计算基准站间空间相关误差,当基准站间模糊度固定后,实时计算生成各基准站间双差相位观测值 双差伪距观测值 双差电离层延迟 双差对流层延迟 和轨道等综合误差 步骤6,通过线性内插方式内插出各三角网元内主参考站与格网点处虚拟基准站间双差电离层残差 双差对流层残差 和双差轨道残差 在此基础上,通过主参考站伪距观测值和载波相位观测值计算得到虚拟基准站伪距观测值P(VRS)和载波相位观测值 步骤7,重复步骤6,依次计算该三角网元内所有格网点虚拟基准站的伪距观测值和载波相位观测值,并将该格网点虚拟基准站的伪距观测值、载波相位观测值与格网点行列号进行一一映射;步骤8,重复步骤4-7,依次计算所有格网点处虚拟基准站伪距观测值和载波相位观测值;步骤9,当有用户进行服务请求时,首先验证其身份,验证失败则直接结束,验证通过时,根据用户上传的经纬度信息(B i ,L i ),与区域最小经纬度BL min (B min ,L min )作差并除以间隔inval,计算出离用户位置最近的格网点位置(H i ,R i );步骤10,将格网点(H i ,R i )所映射的差分改正数产品发送给用户进行定位。
2.根据权利要求1所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,最大行列号(H max ,R max )计算公式为:其中,MIN()为取最小值函数,MAX()为取最大值函数,int[]为取整函数,B i 表示纬度,L i 表示经度。
3.根据权利要求2所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,网格点经纬度GW (i,j) ,通过下式计算:GW (i,j) =(MIN(B i )+i*inval,MIN(L j )+j*inval)(0≤i≤H max ,0≤j≤R max )其中i表示行号、j表示列号。
4.根据权利要求1所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,步骤4具体包括以下步骤:第一步,利用MW组合固定宽巷模糊度,得到MW组合双差观测方程;第二步,利用无电离层组合模型通过自适应卡尔曼滤波估算得到L1模糊度 浮点解;第三步,通过LAMBDA算法对L1模糊度进行搜索固定,对L1模糊度固定后,利用已知的宽巷模糊度,得到L2模糊度 固定解。
5.根据权利要求4所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,L1波段模糊度计算公式为:式中 为L1双差模糊度, 为宽巷模糊度,λ n 、λ w 分别为窄巷、宽巷波长, 分别为双差星地距、双差轨道误差、双差对流层延迟、双差多路径效应;f 1 、f 2 为L1、L2频率,φ 1 、φ 2 为L1、L2载波观测值。
6.根据权利要求5所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于, 由式(9)计算,实时计算生成基准站间双差对流层延迟和轨道误差时合在一起计算,由式(10)计算:N 1 、N 2 分别为L1、L2模糊度,c为光速, 分别为L1、L2频点的载波相观测值,λ 1 、λ 2 分别为L1、L2波长。
7.根据权利要求5所述的基于行列号机制的网络RTK用户请求快速匹配方法,其特征在于,离用户位置最近的格网点位置(H i ,R i )为:其中(B i ,L i )为用户上传的经纬度信息,(B min ,L min )为区域最小经纬度,int[]为取整函数。