本发明公开了一种基于单基站的虚拟参考站的构建方法。它包括如下步骤，步骤一：在流动站概略位置附近选定一虚拟基准站位置，计算虚拟基准站、基准站与相同卫星的几何距离之差；步骤二：根据虚拟基准站位置、基准站位置，采用对流层延迟模型，计算虚拟基准站、基准站位置的对流层延迟；步骤三：进而计算卫星射线方向的对流层斜延迟；步骤四：进一步计算相同卫星的对流层延迟之差；步骤五：将基准站对应的卫星观测值改正上述的几何距离之差和对流层斜延迟之差，即生成虚拟参考站观测站。本发明具有能实现单基站RTK虚拟基准站的生成，提高大高差情形下的单基准站RTK精度的优点。
1.一种基于单基站的虚拟参考站的构建方法，其特征在于：考虑几何位置变化、以及位置变化引起的对流层延迟的变化，计算几何距离变化，计算对流层延迟变化，且都进行改正，分别是 ΔTi；适用于只需要一个基准站情形下的虚拟基准站生成；
步骤一：在流动站概略位置附近选定一虚拟基准站位置，计算虚拟基准站、基准站与相同卫星的几何距离之差；
步骤二：根据虚拟基准站位置、基准站位置，采用对流层延迟模型，计算虚拟基准站、基准站位置的对流层延迟；
式(1)中：Ti′为虚拟基准站对流层斜延迟；Ti为基准站对流层斜延迟；ΔTi表示虚拟基准站对流层斜延迟与基准站对流层斜延迟之差；
步骤五：将基准站对应的卫星观测值改正上述的几何距离之差和对流层斜延迟之差，即生成虚拟参考站观测站，虚拟基准站的载波相位测量与伪距测量的生成公式如下：
式(2)、式(3)中：S1表示卫星，下标i和i′分别表示基准站与虚拟基准站； 表示虚拟基准站的相位观测值； 表示虚拟基准站的伪距观测值； 为虚拟基准站、基准站与卫星的几何距离之差； 为虚拟基准站与卫星的几何距离；ΔTi为虚拟基准站对流层斜延迟与基准站对流层延迟之差；c为光速；δti′与 分别为虚拟基准站接收机钟差和卫星钟差；
Ti′为虚拟基准站对流层延迟； 为载波波长； 为模糊度参数；改正方法为：先计算ΔTi,再改正到基准站的两个观测值上，对单基站的对流层延迟改正，顾及了对流层延迟的改正，改正因基准站和流动站因高差过大引起的两个站对流层延迟不一致，适用于基准站与流动站高差较大的情况，对流层延迟的改正数ΔTi，在步骤二‑步骤四计算出来，在步骤五中进行改正。
2.根据权利要求1所述的基于单基站的虚拟参考站的构建方法，其特征在于：在步骤一中，虚拟基准站、基准站与相同卫星的几何距离之差的计算方法，包括如下步骤：假设已知，第i测站基准站的真实坐标(Xi，Yi，Zi)，第i测站虚拟基准站的真实坐标(X′i，Y′i，Z′i)，卫星的精确位置为 广播星历中已知的卫星位置为因为S2未知，卫星星历误差在信号传播路径上的投影δρ为：
3.根据权利要求2所述的基于单基站的虚拟参考站的构建方法，其特征在于：在步骤二中，对流层延迟模型采用的方式包括：在测区内部署多个地面气象观测站，进行实测的对流层延迟建模；
基于虚拟基准站和基准站高程，利用气象参数垂直递减模型得到；
在基准站和流动站配备气象观测，通过实测气象参数对流层延迟模型得到。
4.根据权利要求3所述的基于单基站的虚拟参考站的构建方法，其特征在于：在步骤二中，虚拟基准站、基准站位置的对流层延迟的计算方法，具体包括如下步骤：当选取GPT2w模型，计算第i测站基准站与第i测站虚拟基准站的对流层延迟；通过转换可得基准站纬度dlat1，基准站海拔高度h1，虚拟基准站纬度dlat2，基准站海拔高度h2；
式(8)、(9)、(10)、(11)中：ZHD1和ZWD1分别为基准站静力学延迟和湿延迟；ZHD2和ZWD2分‑1别为虚拟基准站静力学延迟和湿延迟；k′2、k3为大气折射率常数,k′2的值为16.529k·mb ，‑1k3取值3.776×105k·mb ；Tm1和Tm2为基准站和虚拟基准站处的大气加权平均温度，单位为K；gm为重力加速度；Rd为干空气比气体常数；λ1和λ2为水汽压递减率，由利用测站处的气象廓线数据拟合而来，或由GPT2w模型给出；e1和e2为基准站和虚拟基准站处的水汽压，单位为Pa；P1和P2为基准站和虚拟基准站处的大气压，单位为Pa；
其他未知量由GPT2w模型计算得，GPT2w模型可输出气压、温度、下降率、水气压以及VMF1映射函数系数；GPT2w模型利用最小二乘法估计平均值A0、年度值(A1,B1)、半年度值(A2,B2)的变化；参数r(t)由如下公式推导：式(12)中：doy代表年积日；A0、A1、B1、A2和B2均为格网化系数；
式(13)中：T0和P0为格网点参考高度的温度与气压；T和P是归化到测站高度的温度与气压；dT为温度下降率；dZ为待求高度与各网点参考高度之间的高差；Q为比湿；e为水气压；gm为重力加速度；e0为格网点水气压；λ为水汽压递减因子；dM0为干燥空气摩尔质量，其值为‑3
28.965×10 kg/mol；Rg＝8.3143J/(K·mol)为通用气体常数；
从GPT2w模型提供的格网文件中获取临近四个格网点的系数；在格网点数据归化到测站高处后，使用采用双线性内插将格网点处气象参数内插到测站处，便可得到测站处的气象参数；
式(14)中：Q0为待插值点气象参数；(x1,y1)和(x2,y2)为临近格网的左下角和右上角格网点坐标，(x,y)为待插值点坐标；Q11、Q21、Q12、Q22为GPT2w模型提供的四个格网点不同参数的值。
5.根据权利要求3或4所述的基于单基站的虚拟参考站的构建方法，其特征在于：在步骤三中，卫星射线方向的对流层斜延迟包括静力学延迟和湿延迟，具体计算方法如下：Ti＝(ZHDi+ZWDi)·mfi (15)
Ti′＝(ZHDi′+ZWDi′)·mfi ′ (16)式(15)、(16)中：Ti，Ti′为基准站与虚拟基准站位置的对流层斜延迟；mfi，mfi′为基准站与虚拟基准站位置的映射函数；ZHDi，ZHDi′为基准站与虚拟基准站位置的天顶对流层静力学延迟；ZWDi，ZWDi′为基准站与虚拟基准站位置的天顶对流层湿延迟。
一种基于单基站的虚拟参考站的构建方法
技术领域
[0001] 本发明涉及一种基于单基站的虚拟参考站的构建方法。
背景技术
[0002] 通常，在网络RTK中，为消除流动站处的相关误差（包括对流层延迟、电离层延迟、残余轨道误差等），可以在流动站附近生成一虚拟参考站，虚拟参考站的精确坐标已知，GNSS观测值通过网络RTK技术，基于周围的多个GNSS基准站计算而来。相关计算细节可参考Wanninger (1995, 1997)，姚宜斌等(2016)。虚拟参考站技术的本质是在流动站附近生成一组与流动站具有相似误差的观测值，相关误差必须通过周围站点计算得到；
[0003] 然而，也存在单基站情况下可能出现需要生成虚拟参考站的情况：（1）在短距离大高差单基站RTK作业时，需要基于单基站和相关对流层延迟信息，生成一组与流动站相同（或相近）高程的虚拟参考站。（2）在网络RTK服务中，若生成的虚拟参考站高程与流动站高程差异较大，也需要生成与流动站相同（或相近）高程的虚拟参考站。（3）在单基准RTK作业时，若基准站坐标涉密，无法以明文方式传输，需要对观测站进行加偏，生成虚拟参考站对外提供服务；
[0004] 当采用单基站为流动站提供服务时，因为只有一个基准站，无法采用常规的虚拟参考站计算方法生成虚拟参考站。此外，现有公开号为CN105929424A，一种基于GNSS基准站的虚拟化算法，其由于只对几何距离进行了改正，没有考虑对流层延迟的改正，只适用于基准站与虚拟基准站对流层延迟差异较小的情况；
[0005] 因此，开发一种适用于大高差环境的基于单基站的虚拟参考站很有必要。
发明内容
[0006] 本发明的目的是为了提供一种基于单基站的虚拟参考站的构建方法，实现单基站RTK虚拟基准站的生成，解决网络RTK或单基站RTK中需要生成虚拟参考站的问题，提高大高差情形下的采用虚拟基准站进行RTK定位的精度。
[0007] 为了实现上述目的，本发明的技术方案为：一种基于单基站的虚拟参考站的构建方法，其不同之处在于：包括如下步骤，
[0008] 步骤一：在流动站概略位置附近选定一虚拟基准站位置，计算虚拟基准站、基准站与相同卫星的几何距离之差；
[0009] 步骤二：根据虚拟基准站位置、基准站位置，采用对流层延迟模型，计算虚拟基准站、基准站位置的对流层延迟；
[0010] 步骤三：进而计算卫星射线方向的对流层斜延迟；
[0011] 步骤四：进一步计算相同卫星的对流层延迟之差；
[0012] 相同卫星的对流层延迟之差的计算公式为：
[0013] （1）
[0014] 式（1）中： 为虚拟基准站对流层斜延迟； 为基准站对流层斜延迟； 表示虚拟基准站对流层斜延迟与基准站对流层斜延迟之差；
[0015] 步骤五：将基准站对应的卫星观测值改正上述的几何距离之差和对流层斜延迟之差，即生成虚拟参考站观测站，
[0016] 虚拟基准站的载波相位测量与伪距测量的生成公式如下：
[0017] （2）
[0018] （3）
[0019] 式（2）、式（3）中： 表示卫星，下标 和 分别表示基准站与虚拟基准站； 表示虚拟基准站的相位观测值； 表示虚拟基准站的伪距观测值； 为虚拟基准站与卫星的几何距离； 为虚拟基准站、基准站与卫星的几何距离之差； 为虚拟基准站对流层斜延迟与基准站对流层延迟之差；为光速； 与 分别为虚拟基准站接收机钟差和卫星钟差； 为虚拟基准站对流层延迟； 为载波波长； 为模糊度参数；
[0020] 本发明的改正方法为：先计算 、 , 再改正到基准站的两个观测值上，本发明对单基站的对流层延迟改正，顾及了对流层延迟的改正（对流层延迟改正的作用是改正因基准站和流动站因高差过大引起的两个站对流层延迟不一致），适用于基准站与流动站高差较大的情况。本发明中的对流层延迟的改正数 ，在步骤二‑步骤四计算出来，在步骤五中进行改正。
[0021] 在上述技术方案中，在步骤一中，虚拟基准站、基准站与相同卫星的几何距离之差的计算方法，包括如下步骤：
[0022] 假设已知，第i测站基准站的真实坐标 ，第i测站虚拟基准站的真实坐标，卫星的精确位置为 ，广播星历中已知的卫星位置为；因为S2未知，卫星星历误差在信号传播路径上的投影 为：
[0023] （4）
[0024] 基准站到卫星的几何距离 为：
[0025] （5）
[0026] 虚拟基准站到卫星的几何距离 为：
[0027] （6）
[0028] 几何距离之差 为：
[0029] （7）。
[0030] 在上述技术方案中，在步骤二中，对流层延迟模型可以采用多种方式，对流层延迟模型采用的方式包括：在测区内部署多个地面气象观测站，进行实测的对流层延迟建模；
[0031] 或基于经验对流层延迟模型，如GPT2w、UNB3等；
[0032] 或基于虚拟基准站和基准站高程，利用气象参数垂直递减模型得到；
[0033] 或在基准站和流动站配备气象观测，通过实测气象参数对流层延迟模型得到；
[0034] 因为基准站之间的距离较近，本方法忽略电离层延迟的差异。
[0035] 在上述技术方案中，在步骤二中，虚拟基准站、基准站位置的对流层延迟的计算方法，具体包括如下步骤：
[0036] 当选取GPT2w模型，计算第i测站基准站（该测站基准站的真实坐标为 ）与第i测站虚拟基准站（该测站虚拟基准站的真实坐标为 ）的对流层延迟；通过转换可知基准站纬度dlat1，基准站海拔高度h1，虚拟基准站纬度dlat2，基准站海拔高度h2；
[0037] 计算得到基准站对流层延迟为：
[0038] （8）
[0039] （9）
[0040] 计算得到虚拟基准站对流层延迟为：
[0041] （10）
[0042] （11）
[0043] 式（8）、（9）、（10）、（11）中： 和 分别为基准站静力学延迟和湿延迟；和 分别为虚拟基准站静力学延迟和湿延迟； 、k3为大气折射率常数, 的值‑1 ‑1
为16.529 k•mb ，k3取值3.776×105 k•mb ； 和 为基准站和虚拟基准站处的大气加权平均温度，单位为K；gm为重力加速度；Rd为干空气比气体常数；和 为水汽压递减率，由利用测站处的气象廓线数据拟合而来，或由GPT2w模型给出；e1和e2为基准站和虚拟基准站处的水汽压，单位为Pa；P1和P2为基准站和虚拟基准站处的大气压，单位为Pa；
[0044] 其他未知量可由GPT2w模型计算得，GPT2w模型可输出气压、温度、下降率、水气压以及VMF1映射函数系数；GPT2w模型利用最小二乘法估计平均值A0、年度值（A1,B1）、半年度值（A2,B2）的变化；参数r(t)由如下公式推导：
[0045] （12）
[0046] 式（12）中： 代表年积日， 和 为格网化系数；
[0047] 各气象参数的高度改正如下：
[0048] （13）
[0049] 式（13）中：T0和P0为格网点参考高度的温度与气压；T和P是归化到测站高度的温度与气压；dT为温度下降率；dZ为待求高度与各网点参考高度之间的高差；Q为比湿；e为水气2 ‑3
压；gm=9.80665m/s 为重力加速度； 为干燥空气摩尔质量，其值为28.965×10 kg/mol；
为通用气体常数；e0为格网点水气压；λ为水汽压递减率；
[0050] 从GPT2w模型提供的格网文件中获取临近四个格网点的系数；在格网点数据归化到测站高处后，使用采用双线性内插将格网点处气象参数内插到测站处，便可得到测站处的气象参数；双线性内插公式为：
[0051] （14）
[0052] 式（14）中：Q0为待插值点气象参数；(x1, y1)和(x2, y2)为临近格网的左下角和右上角格网点坐标，(x, y)为待插值点坐标；Q11、Q21、Q12、Q22为GPT2w模型提供的四个格网点不同参数的值。
[0053] 在上述技术方案中，在步骤三中，卫星射线方向的对流层斜延迟包括静力学延迟和湿延迟，具体计算方法如下：
[0054] （15）
[0055] （16）
[0056] 式（15）、（16）中： 为基准站与虚拟基准站位置的对流层斜延迟；为基准站与虚拟基准站位置的映射函数； 为基准站与虚拟基准站位置的天顶对流层静力学延迟； 为基准站与虚拟基准站位置的天顶对流层湿延迟。
[0057] 本发明适用于短距离大高差环境；本发明适用的短距离是指最大不超过10km，一般在5km以内；本发明适用的大高差是指：高差大于100m，结合地球地表高差实际情况，一般为100m‑1000m；
[0058] 在短距离大高差情形下本发明的RTK定位精度可以提高到厘米级。
[0059] 上述*表示数学符号乘号。
[0060] 本发明具有如下优点：
[0061] (1)常规虚拟参考站技术基于多个GNSS基准站求解生成，不适用于单基准站的情况；本发明基于单基准站生成，适用于只有一个基准站对外提供服务的情况；本发明适用于单基准站（只需要一个基准站）情形下的虚拟基准站生成，且能提高单基准站RTK精度(大高差情形下)；
[0062] (2)基于单基站生成虚拟参考站需要解决两个问题虚拟参考站位置与基准站位置不同导致的卫星几何距离差异在观测值中的改正，以及虚拟参考站位置与基准站位置不同导致的空间误差（主要是对流层误差）差异在观测值中的改正；本发明考虑了几何位置变化、以及位置变化引起的对流层延迟的变化（步骤一中计算了几何距离变化，步骤四中计算了对流层延迟变化，且都在步骤五中进行了改正（分别是delta_rho， delta_Ti）），提高了单基准站RTK精度；克服了现有的虚拟化技术，只考虑了几何位置便宜，未顾及位置变化引起的对流层延迟的变化的缺陷；
[0063] (5)本发明适用于短距离大高差环境，且RTK定位精度高（本发明的适用距离最大不超过10km，一般在5km以内；本发明适用的高差大于100m，结合地表高差实际，一般为100m‑1000m；本发明的RTK定位精度为cm级），解决上述短距离大高差RTK定位精度差的问题（常规的RTK技术不适用于短距离大高差环境（如，距离小于10km、高差大于100m的环境），且高差越大，精度越差，常规的RTK技术定位精度一般在分米级）。
附图说明
[0064] 图1为本发明的流程图；
[0065] 图2为本发明应用于大高差单基站RTK的虚拟基准站服务的流程图；
[0066] 图3为本发明应用于网络RTK中虚拟基准站的高程改正的流程图；
[0067] 图4为本发明应用于单基站坐标加偏RTK服务方式的流程图。
具体实施方式
[0068] 下面结合附图详细说明本发明的实施情况，但它们并不构成对本发明的限定，仅作举例而已。同时通过说明使本发明的优点更加清楚和容易理解。
[0069] 参阅附图1可知：一种基于单基站的虚拟参考站的构建方法，包括如下步骤，[0070] 步骤一：在流动站概略位置附近选定一虚拟基准站位置，虚拟基准站位置坐标精确已知，同时，基准站坐标精确已知；计算虚拟基准站、基准站与相同卫星的几何距离之差；
[0071] 步骤二：根据虚拟基准站位置、基准站位置，采用对流层延迟模型，计算虚拟基准站、基准站位置的对流层延迟；
[0072] 步骤三：进而计算卫星射线方向的对流层斜延迟；
[0073] 步骤四：进一步计算相同卫星的对流层延迟之差；
[0074] 相同卫星的对流层延迟之差的计算公式为：
[0075] （1）
[0076] 式（1）中： 为虚拟基准站对流层斜延迟； 为基准站对流层斜延迟； 表示二者之差；
[0077] 步骤五：将基准站对应的卫星观测值（包括伪距和相位）改正上述的几何距离之差和对流层斜延迟之差，即生成虚拟参考站观测站，
[0078] 虚拟基准站的载波相位测量与伪距测量的生成公式如下：
[0079] （2）
[0080] （3）
[0081] 式（2）、式（3）中： 表示卫星，下标 和 分别表示基准站与虚拟基准站； 表示虚拟基准站的相位观测值； 表示虚拟基准站的伪距观测值； 为虚拟基准站、基准站与卫星的几何距离之差； 为虚拟基准站与卫星的几何距离； 为虚拟基准站对流层斜延迟与基准站对流层延迟之差，由公式（1）计算而来；为光速； 与 分别为虚拟基准站接收机钟差和卫星钟差； 为虚拟基准站对流层延迟； 为载波波长； 为模糊度参数。
[0082] 进一步地，在步骤一中，虚拟基准站、基准站与相同卫星的几何距离之差的计算方法，包括如下步骤：
[0083] 假设已知，第i测站基准站的真实坐标 ，第i测站虚拟基准站的真实坐标，卫星的精确位置为 ，广播星历中已知的卫星位置为；因为S2未知，卫星星历误差在信号传播路径上的投影 为：
[0084] （4）
[0085] 基准站到卫星的几何距离 为：
[0086] （5）
[0087] 虚拟基准站到卫星的几何距离 为：
[0088] （6）
[0089] 几何距离之差 为：
[0090] （7）。
[0091] 进一步地，在步骤二中，对流层延迟模型可以采用多种方式，对流层延迟模型采用的方式包括：在测区内部署多个地面气象观测站，进行实测的对流层延迟建模；
[0092] 基于经验对流层延迟模型，如GPT2w、UNB3等；
[0093] 基于虚拟基准站和基准站高程，利用气象参数垂直递减模型得到；
[0094] 在基准站和流动站配备气象观测，通过实测气象参数对流层延迟模型得到；
[0095] 因为距离较近，本方法忽略电离层延迟的差异。
[0096] 进一步地，在步骤二中，虚拟基准站、基准站位置的对流层延迟的计算方法，具体包括如下步骤：
[0097] 假设选取GPT2w模型，计算第i测站基准站 与第i测站虚拟基准站的对流层延迟；通过转换可知基准站纬度dlat1，基准站海拔高度h1，虚拟基准站纬度dlat2，基准站海拔高度h2；
[0098] 计算得到基准站对流层延迟为：
[0099] （8）
[0100] （9）
[0101] 计算得到虚拟基准站对流层延迟为：
[0102] （10）
[0103] （11）
[0104] 式（8）、（9）、（10）、（11）中： 和 分别为基准站静力学延迟和湿延迟；和 分别为虚拟基准站静力学延迟和湿延迟； 、k3为大气折射率常数, 的值‑1 ‑1
为16.529 k•mb ，k3取值3.776×105 k•mb ； 和 为基准站和虚拟基准站处的大气加权平均温度，单位为K；gm为重力加速度；Rd为干空气比气体常数；和 为水汽压递减率，由利用测站处的气象廓线数据拟合而来，或由GPT2w模型给出；e1和e2为基准站和虚拟基准站处的水汽压，单位为Pa；P1和P2为基准站和虚拟基准站处的大气压，单位为Pa；
[0105] 其他未知量可由GPT2w模型计算得，GPT2w模型可输出气压、温度、下降率、水气压以及VMF1映射函数系数。GPT2w模型利用最小二乘法估计平均值A0、年度值（A1,B1）、半年度值（A2,B2）的变化；参数r(t)由如下公式推导：
[0106] （12）
[0107] 式（12）中： 代表年积日， 和 为格网化系数；
[0108] 各气象参数的高度改正如下：
[0109] （13）
[0110] 式（13）中：T0和P0为格网点参考高度的温度与气压；T和P是归化到测站高度的温度与气压；dT为温度下降率；dZ为待求高度与各网点参考高度之间的高差；Q为比湿；e为水气2 ‑3
压；gm=9.80665m/s 为重力加速度； 为干燥空气摩尔质量，其值为28.965×10 kg/mol；
为通用气体常数；e0为格网点水气压；λ为水汽压递减率；
[0111] 从GPT2w模型提供的格网文件中获取临近四个格网点的系数；在格网点数据归化到测站高处后，使用采用双线性内插将格网点处气象参数内插到测站处，便可得到测站处的气象参数；双线性内插公式为：
[0112] （14）
[0113] 式（14）中：Q0为待插值点气象参数，如气温、气压或水汽压；(x1, y1)和(x2, y2)为临近格网的左下角和右上角格网点坐标，(x, y)为待插值点坐标；Q11、Q21、Q12、Q22为GPT2w模型提供的四个格网点不同参数的值。
[0114] 进一步地，在步骤三中，卫星射线方向的对流层斜延迟包括静力学延迟和湿延迟，具体计算方法如下：
[0115] （15）
[0116] （16）
[0117] 式（15）、（16）中： 为基准站与虚拟基准站位置的对流层斜延迟； 为基准站与虚拟基准站位置的映射函数； 为基准站与虚拟基准站位置的天顶对流层静力学延迟； 为基准站与虚拟基准站位置的天顶对流层湿延迟。实施例
[0118] 实施例1：大高差单基站RTK的虚拟基准站服务
[0119] 本实施例如图2所示，为本申请用于大高差情形下的实施例。首先，流动站接收机通过4G通信向服务端发送自身概略位置；其次，以基准站经纬度和流动站高程作为虚拟基准站位置坐标；再次，服务端利用基准站伪距和相位观测值，通过本申请方法计算虚拟基准站观测值，并通过4G通信联合虚拟基准站位置发送给流动站接收机；最后，流动站接收机接收到虚拟基准站观测值后，可以利用常规RTK定位；
[0120] 结论：本实施例采用本发明方法按照上述先后执行的步骤，使传统RTK接收机（即流动站）不做任何改动，即可解决大高差引起的对流层延迟差异过大，导致的定位精度差的问题。
[0121] 实施例2：网络RTK中虚拟基准站的高程改正
[0122] 本实施例如图3所示，为本发明所述方法应用于网络RTK服务中的高程修正的实施例；常规的网络RTK服务首先获取流动站概略位置后，忽略高程因素在流动站服务生成虚拟基准站；对于流动站与虚拟基准站存在大高差的情况，会存在较大的对流层延迟误差；因此，可利用本发明方法根据流动站发来的概略位置，顾及虚拟基准站和流动基准站的高程差异，生成新的虚拟基准站；最后，将新的虚拟基准站数据通过4G通信发送给流动站；
[0123] 结论：本实施例采用本发明方法按照上述先后执行的步骤，顾及了对流层延迟的改正，适用于基准站与流动站高差较大的情况，解决了传统网络RTK中，由于未顾及流动站高程与虚拟参考站高程差异过大引起的流动站定位精度差的问题。
[0124] 实施例3：单基站坐标加偏RTK服务方式
[0125] 本实施例如图4所示，为本发明应用于需要隐藏基准站真实坐标对外提供服务时对基准站坐标进行加偏的实施例。首先，基准站获取接收机伪距和相位观测值以及基准站坐标；其次，在基准站真实坐标上随机加上偏移数（考虑到RTK定位的适用距离和地表实际情况，水平偏移一般不超过5km，垂直偏移不超过1000m），得到虚拟基准站坐标；其次，利用本方法生成虚拟基准站处的伪距和相位虚拟观测值；最后，将生成的伪距和相位虚拟观测值及虚拟基准站坐标播发给流动站；
[0126] 结论：本实施例采用本发明方法按照上述先后执行的步骤，可以实现加密基准站坐标对外服务，不仅可以对水平坐标加偏，还可以对高程坐标加偏；
[0127] 与现有技术相比，本实施例采用本发明方法不仅可以对平面坐标进行加偏，还可以对高程进行加偏（即能实现立体加偏），彻底隐藏基准站坐标对外提供服务；克服了现有技术由于只对几何距离进行了改正，没有考虑对流层延迟的改正，只适用于基准站与虚拟基准站对流层延迟差异较小的情况，仅能对平面坐标进行加偏，无法实现立体加偏的缺陷。
[0128] 其它未说明的部分均属于现有技术。
