CN101449179A - 从区域、广域或全球载波相位差分导航（wadgps）向本地实时动态（rtk）导航系统过渡时增加位置信息可靠性的方法 - Google Patents

Abstract

本发明包括一种组合使用本地定位系统、本地RTK系统与区域、广域或全球差分载波相位定位系统(WADGPS)的方法，以避免与所述本地定位系统、所述RTK和所述WADGPS导航技术在单独使用时可能存在的缺点。所述方法包括：根据来自WADGPS的信息确定物体的第一位置，及根据来自本地定位/RTK定位系统的信息确定物体的第二位置；之后，比较通过所述WADGPS确定的位置与通过所述本地定位/RTK定位系统确定的位置。当所述WADGPS位置与本地定位/RTK位置相差大于一预定阈值时，所述WADGPS被用于导航物体，而当所述WADGPS位置与本地定位/RTK位置相差小于所述预定阈值时，所述本地定位/RTK位置被用于导航物体。

Description

从区域、广域或全球载波相位差分导航(WADGPS)向本地实时动态(RTK)导航系统过渡时增加位置信息可靠性的方法

交叉引用相关案子

本申请是于2006年1月31日递交的申请号为11/345,124申请的部分继续申请，所述申请号为11/345,124的申请是2004年1月13日递交的申请号为10/757,340申请的部分继续申请，所述两个申请的内容通过参考被引用到本申请中。

技术领域

本发明一般涉及与使用卫星定位和导航相关联的技术，并且更具体地说，涉及在区域、广域或全球载波相位定位和/或导航系统中求解载波浮点模糊度。

背景技术

全球定位系统(GPS)使用太空中的卫星定位地球上的物体。利用GPS，信号从卫星到达GPS接收机，并用于确定GPS接收机的位置。目前，对应于具有锁定GPS卫星信号的每个相关器信道的两种类型的GPS测量，可提供用于民用GPS接收机。所述两种类型的GPS测量是伪距和两个载波信号L1和L2的综合载波相位，这两个载波信号L1和L2分别具有频率1.5754GHz和1.2276GHz或波长0.1903米和0.2442米。伪距测量(或代码测量)是所有类型的GPS接收机可进行的基本的GPS可观察量。它利用调制到载波信号上的C/A或P代码。测量记录相关代码从卫星传播到接收机所费的视时，即根据接收机时钟信号到达接收机的时间减去根据卫星时钟信号离开卫星的时间。载波相位测量是通过综合信号到达接收机时信号的重建载波而获得。因此，载波相位测量也是如根据卫星时钟信号离开卫星的时间和根据接收机时钟信号到达接收机的时间所确定的转接时间差的量度。然而，由于在接收机开始跟踪信号的载波相位时在卫星与接收机之间转接中的初始整周数通常未知，因此，转接时间差可能有多个载波周期的误差，即，在载波相位测量中有整周模糊度(whole-cycle ambiguity)。

在GPS测量可用时，在GPS接收机与多个卫星中的每一卫星之间的距离可通过将信号传播时间乘以光速而计算得出。这些距离通常称为伪距(假距离)，这是因为接收机时钟一般具有大的时间误差，导致测量的距离有常见的偏差。作为普通导航计算的一部分求出由接收机时钟误差引起的这种常见的偏差和接收机的位置坐标。各种其它因素也可导致计算的距离中出现误差或噪声，包括星历表误差、卫星时钟计时误差、大气效应、接收机噪声和多径误差。利用孤立GPS导航，其中具有GPS接收机的用户获得相对于视野中的多个卫星的代码和/或载波相位距离而不咨询任何参考站，用户在减少距离中出现的误差或噪声方面极为受限。

为了消除或减少这些误差，典型地在GPS应用中使用差分操作。差分GPS(DGPS)操作典型地包括基准参考GPS接收机、用户(或导航)GPS接收机及在用户与参考接收机之间的通信链路。参考接收机放在已知位置，并且该已知位置用于生成与一些或所有上述误差因素相关联的校正。所述校正被提供给所述用户接收机，然后所述用户接收机使用所述校正适当地校正其计算的位置。所述校正的形式可以是对在参考站点确定的参考接收机位置的校正或者对特定GPS卫星时钟和/或轨道的校正。使用载波相位测量的差分操作经常称为实时动态(RTK)定位/导航操作。

差分GPS(DGPS)的基本概念是利用GPS测量中固有的误差的空间和时间相关，以消除在伪距和/或载波相位测量中由这些误差因素导致的噪声因素。然而，虽然作为在伪距或载波相位测量上的偏差出现的GPS卫星时钟计时误差在参考接收机与用户接收机之间优选是相关的，但大多数其它误差因素是不相关的，或者在广域应用中，即当在参考接收机与用户接收机之间的距离变大时，所述相关会消除。

为克服在广域应用中DGPS系统的不准确性，开发了各种区域、广域或全球DGPS(以下称为广域DGPS或WADGPS)技术。WADGPS包括与计算中心通信的多个参考站的网络。基于参考站的已知位置和由它们进行的测量，误差校正在中心被计算出。所述计算出的误差校正随后通过诸如卫星、电话或无线电的通信链路传送到用户。通过使用多个参考站，WADGPS提供更准确估计的所述误差校正。

因此，已开发多种不同的技术以使用GPS载波相位测量获得高准确度的差分导航。所述的RTK技术具有约1厘米的典型准确度。然而为获得该准确度，必须确定所述差分载波相位测量中的整周模糊度。当在用户接收机与参考接收机之间的距离(基线距离)短时，RTK技术非常有利，因为在这种情况下可准确和快速地求解整周模糊度。另一方面，当所述基线距离超过几十千米时，可能无法确定所述整周模糊度，普通的RTK准确度无法实现。所述RTK技术的另一个限制是它需要在参考接收机与导航接收机之间维持本地无线电链路。

采用载波相位差分方法的WADGPS技术也可实现很高的导航准确度。所述WADGPS差分技术的特征还在于可靠的长距离低频率通信链路或可靠的卫星通信链路。因此，校正一般可传递到导航接收机而无大的中断。然而，所述WADGPS技术通常将整周模糊度视为实值(非整数)变量，并求出“浮点模糊度”，在获得经过大的卫星几何变化的时间间隔的测量数据前，所述浮点模糊度的定义通常很差。因此，在WADGPS应用中，求出“浮点模糊度”经常需要长达1或2小时时间间隔以产生导航位置的少于10厘米的准确度。

采用本地定位系统能够也可以获得准确度定位(<1cm)。常规的本地定位系统利用有源或无源元件，包括声波和激光测距系统，例如，基于信号的飞行时间和/或多普勒频移。声波系统典型地使用地标和/或应答信标来测量网络内的装置的距离，所述装置中的一些被固定以形成本地坐标系。不幸地，由于声传播穿过空气的特性，声波系统只能够测量距离具有1厘米或更高准确度，并只能是相对短的距离。基于激光的本地定位系统利用装置与诸如棱镜的一个或多个反射物体之间的角度和距离的测量，以三角测量或三边测量所述装置的位置。基于激光的本地定位系统的有效的操作范围也典型地被限制在相对短的距离(约1000-10,000米)。

发明内容

本发明包括一种用于组合使用RTK和WADGPS导航技术的方法，以便可通过另一技术的优点补充每个技术的弱点。WADGPS技术的主要不足之处在于导航接收机花费长的经过时间(经常超过1小时)确定将载波相位测量转换成准确的距离测量所需的浮点模糊度值。RTK技术的主要不足之处在于它需要在用户GPS接收机与参考GPS接收机之间的实时(通常为直线对传)数据链路，并且整周模糊度只可当在参考GPS接收机与用户GPS接收机之间的间距相对较短时确定。

通过使用根据本发明一个实施例的用于组合使用RTK和WADGPS导航技术的方法，可消除这些各自的不足之处。该方法包括使用用户接收机的已知位置在WADGPS系统中初始化浮点模糊度值。当用户接收机已固定时，用户接收机的已知位置可以是测量的位置或从先前操作获得的位置。当用户接收机在移动时，可使用RTK系统获得所述已知位置。

因此，在组合操作中，当用于本地定位系统和/或RTK导航的通信链路可用时，可使用所述本地定位系统和/或所述RTK系统获得所述用户接收机的位置、速度和时间(PVT)输出，而WADGPS系统在后台运行，并且其输出不断被初始化以与来自所述RTK系统的输出一致。当用于所述本地定位系统的通信链路丢失时，可使用所述RTK系统和/或所述WADGPS系统获得所述用户接收机的所述PVT输出，在所述本地定位系统操作的同时该系统已被初始化。当用于所述RTK导航的通信链路丢失时，或者当用户接收机偏离所述RTK系统中的参考站太远时，可使用所述WADGPS系统获得所述用户接收机的所述PVT输出，在所述RTK操作的同时该系统已被初始化。初始化避免了当用户GPS接收机的位置未知时求出浮点模糊度值所需的一般15分钟到2个小时的“捕捉(Pull-in)”时间。这在本地定位系统和/或RTK系统不可用或不准确时从WADGPS系统提供了很准确的PVT解，并且使WADGPS技术对实时高准确度定位和导航用途更实用。

附图说明

图1是根据本发明一个实施例的WADGPS系统、本地定位系统和本地RTK系统的组合的方框图；

图2是藕合到用户GPS接收机的计算机系统方框图；

图3A是表示用于组合使用WADGPS系统、本地RTK系统和/或本地定位系统的方法的流程图；

图3B是表示用于使用本地RTK系统更新接收机位置的方法的流程图；

图4A是表示使用WADGPS系统和本地RTK系统两者的组合操作的进程流的流程图；

图4B是表示使用本地定位系统、本地RTK系统和/或WADGPS系统的进程流的流程图；

图4C是表示使用具有本地定位系统和/或本地RTK系统的WADGPS系统的进程流的流程图；

图5是表示可使用组合操作的情况的图。

具体实施方案

图1表示根据本发明一个实施例的广域或全球差分GPS(WADGPS)系统100。如图1所示，WADGPS系统100包括参考站120的网络和一个或多个处理中心105，每个参考站具有GPS接收机122。参考站120不断将原始GPS可观察量提供到中心105进行处理。这些可观察量包括GPS代码和载波相位测量、星历表及在参考站120根据从多个卫星110接收的信号获得的其它信息。对于广域DGPS系统，参考站120位于跨诸如大陆的广域101的已知位置，或对于全球DGPS网络，参考站120位于跨全球的已知位置。中心105是处理GPS可观察量和计算DGPS校正的设施。如果提供了多个独立的中心，则优选是它们在地理上分开且平行操作。

所述WADGPS系统100可由一个或多个用户(或用户装置或物体)140利用，每个用户具有用于定位和/或导航用途的用户GPS接收机142。在本发明的一个实施例中，用户140通过RTK无线电链路与附近的参考站120相关联，使得用户接收机142和附近的参考站120形成本地RTK系统150。在一些实施例中，所述用户140也可以与具有一个或多个地标176的本地定位系统174相关联。所述一个或多个地标176可以是有源的或无源的。每一所述一个或多个地标176可以具有GPS接收机122。

系统100还包括常规数据链路(未表示)，用于为将GPS可观察量从参考站120发送到中心105并且为将计算得出的校正从中心105广播到参考站120和用户140提供可靠的传输机制。大陆WADGPS系统通常具有大约3-10个参考接收机，而全球WADGPS系统通常具有大约20-100个参考接收机，以将数据反馈至所述中心105。在本发明的一个实施例中，所述GPS可观察量从参考站120经因特网发送到中心105，并且也经因特网将计算出的校正从所述中心发送到一个或多个地面站(未表示)以上行传输到一个或多个卫星(未表示)，然后，卫星广播计算出的校正以由参考站120和用户接收机142接收。

在本发明的一个实施例中，用户或物体140还配有耦合到用户GPS接收机142的计算机系统144。如图2所示，计算机系统144包括中央处理器(CPU)146、存储器148、一个或多个输入端口154、一个或多个辅助接收机155、一个或多个输出端口156及(任选地)用户接口158，通过一个或多个通信总线152彼此藕合。所述存储器148可包括高速随机存取存储器，并可包括非易失性海量储存器，诸如一个或多个磁盘存储装置或闪存装置。

存储器148优选存储操作系统162、GPS应用程序164和数据库170。如下面更详细描述的一样，GPS应用程序164可包括用于执行用于组合使用本地定位系统174、本地RTK系统150和/或WADGPS系统160的方法300的程序166。存储在存储器148中的操作系统162和应用程序164用于由计算机系统144的CPU146执行。存储器148优选还存储在GPS应用程序164执行期间使用的数据结构，包括GPS伪距和载波相位测量168、从所述中心接收的GPS校正172及在本文档中论述的其它数据结构。

输入端口154用于接收来自GPS接收机142的数据，至少一个辅助接收机155用于经无线电链路124接收来自本地定位系统174或本地RTK系统150中的参考站或地标120的信息，以及用于经卫星链路107接收来自中心105的GPS校正和其它信息。输出端口156用于经无线电链路124或声波或激光装置(未示出)将数据输出到参考站或地标120。在本发明的一个实施例中，计算机系统144的CPU146和存储器148与GPS接收机142集成位于单个外壳内的单个装置，如图2所示。然而，此集成不是执行本发明的方法所必需的。

因此，用户或物体140可同时或在不同时间参与三种不同的操作模式。用户或物体140可在WADGPS模式操作，在该模式中，用户或物体140使用WADGPS系统100为自身定位或导航，可在RTK模式操作，在该模式中，用户或物体140使用本地RTK系统150为自身定位或导航，和/或可在本地定位模式操作，在该模式中，用户或物体140使用本地定位系统174为自身定位或导航。当用户或物体140接近与其相关联的一个或多个地标176并且在所述用户或物体140与所述一个或多个地标176之间的无线电链路可维持时，所述用户可使用所述一个或多个地标176确定自身相对所述一个或多个地标176的位置。当所述用户或物体140接近与其相关联的参考站120并且在所述用户或物体140与所述参考站120之间的无线电链路可维持时，所述用户可使用本地RTK系统150确定自身相对于所述参考站120的位置。所述本地定位系统174和本地RTK系统150比WADGPS系统100更有利，这表现在它更准确并且所述整周整数模糊度可被迅速求解，如下所述。

使用本地RTK系统150，当相对于在参考GPS接收机122和相关联用户GPS接收机142的视野中的n个卫星110进行测量时，测量可用于根据以下阵列格式的方程式求出用户或物体140的位置：

$(\&dtri;\mathrm{\&Phi;}+N)\mathrm{\&lambda;}=\mathrm{Hx}+n_{\mathrm{\&phi;}}---\left(1\right)$

其中， $\&dtri;\mathrm{\&Phi;}={\left[\begin{array}{cccc}{\&dtri;\mathrm{\&phi;}}_1& {\&dtri;\mathrm{\&phi;}}_2& \&CenterDot;\&CenterDot;\&CenterDot;& {\&dtri;\mathrm{\&phi;}}_n\end{array}\right]}^T$ 是由相对于n个卫星110中每个卫星的差分载波相位测量形成的载波相位测量向量；N＝[N₁N₂…N_n]^T是由与载波相位测量向量中每个差分载波相位测量相关联的差分整数模糊度形成的整数模糊度向量；H＝[h₁h₂…h_n]^T是由从用户或物体140到n个卫星110的单位向量形成的测量灵敏度矩阵；x是实未知状态向量(或实向量)，包括本地RTK系统150中从参考站120到用户或物体140的位置向量；以及 $n_{\mathrm{\&phi;}}={\left[\begin{array}{cccc}{n}_{{\mathrm{\&phi;}}_1}& n_{{\mathrm{\&phi;}}_2}& \&CenterDot;\&CenterDot;\&CenterDot;& n_{{\mathrm{\&phi;}}_n}\end{array}\right]}^T$ 是由相对于n个卫星110中每个卫星的差分载波相位噪声形成的测量噪声向量(或相位距离残差向量)。

要使用方程式(1)求出实向量x，需要求解整数模糊度向量N。已开发许多不同的方法来求解在整数模糊度向量N中包括的整数模糊度值，并且这些方法一般使用搜索进程来查找满足某些标准的整数模糊度值组合，如测量残差向量ΔΦ的最小范数，

${\mathrm{\&Delta;}}_{\mathrm{\&Phi;}}=(\&dtri;\mathrm{\&Phi;}+\tilde{N})\mathrm{\&lambda;}-H\hat{x}---\left(2\right)$

其中，Δ_Φ是对应于包括整数模糊度值组合的候选整数模糊度向量

的相位距离残差向量，并且交是方程式(1)的最小二乘解，

$\hat{x}={\left[H^{T}H\right]}^{-1}{H}^T(\&dtri;\mathrm{\&Phi;}+N)\mathrm{\&lambda;}---\left(3\right)$

或者，

$\hat{x}={\left[H^T\mathrm{RH}\right]}^{-1}{H}^T{R}^{-1}(\&dtri;\mathrm{\&Phi;}+N)\mathrm{\&lambda;}---\left(4\right)$

其中

是由σ_i形成的测量协方差矩阵，σ_i是使用常规方法计算得出的差分载波相位噪声

的标准偏差。在“GPS相位和代码观察的精度、互相关和时间相关”(by Peter Bona，GPS Solutions，Vol.4，No.2，Fall2000，p.3-13)或“用于低成本惯性导航的紧密集成姿态确定方法：双天线GPS和GPS磁力计”(by Yang，Y.，Ph.D.Dissertation，Dept.of ElectricalEngineering，University of CaLlfornia，Riverside，CA June 2001)中可找到用于计算σ_i的方法的例子，这两篇文章均通过引用结合于本文中。

在通过引用结合于本文的“瞬时模糊度求解”(＂Instantaneous Ambiguity Resolution＂，by Hatch，R.，in the Proceedings of the KIS Symposium 1990，Banff，Canada)中和在也通过引用结合于本文的专利申请序号为10/338,264的共同拥有的专利申请“用于实时动态观察和导航的快速模糊度求解”(＂Fast Ambiguity Resolution for Real Time KinematicSurvey and Navigation＂)中，可找到搜索方法的其它例子。

利用求解的整数模糊度，用户接收机142的位置、速度和时间(PVT)可作为本地RTK系统150的解来准确地计算。

使用本地定位系统174，用户接收机142的位置、速度、时间(PVT)可作为所述本地定位系统174的解准确地计算。例如，利用信号的飞行时间和/或多普勒频移可以确定相对一个或多个地标的距离和角度信息。在2005年4月11日提交的公布号2005-0270228、名称为“用于本地定位的改进的雷达系统”的美国专利申请中提供进一步说明了确定距离和角度信息，所述专利申请的内容通过引用结合于本文中。此外，使用本地RTK系统150，当相对于在所述一个或多个地标176的一个或多个中的GPS接收机120和相关联用户GPS接收机142的视野中的n个卫星110进行测量时，所述测量可用于根据前面的方程式求出用户或物体140的位置。

尽管它们有许多优势，但由于用户可能移到离所述一个或多个地标176和/或所述参考站120太远的位置或者在所述一个或多个地标176和/或所述参考站120的范围外，使得在用户或物体140与所述地标和/或所述参考站之间的通信链路和/或无线电链路124无法维持，因此，本地定位系统174和/或本地RTK系统150可能对用户或物体140不可用。在这些情况下，无法通过将在用户或物体140的测量与在所述地标176和/或所述参考站120的测量之间的差异考虑在内而满意地消除电离层引起的误差。此误差影响上述搜索整数模糊度向量的进程，这是因为它会导致在测量残差向量Δ_Φ中包括的测量残差增大。

因此，在本地定位系统174和本地RTK系统150由于在用户GPS接收机与所述地标和参考站之间的大间隔而不可用或失去其准确度的情况下，用户可能需要在WADGPS模式操作，在该模式使用不同的方案来求解整数模糊度。使用WADGPS系统100，每个整周模糊度被估计为实值(非整数)变量。此实践经常称为确定“浮点模糊度”值。一个用于确定“浮点模糊度”值的方法涉及基于在用户或物体140进行的原始GPS测量形成折射校正的代码和载波相位测量、将载波相位测量换算成与代码测量相同的单位以及从对应的代码测量减去每个换算的载波相位测量以获得偏移值。在本发明的一个实施例中，标明为P_RC的折射校正的代码测量形成如下：

其中，P₁和P₂是在特殊测量时期分别在L1和L2频率f₁和f₂上的原始伪距代码测量。标明为L_RC的折射校正的载波相位测量以类似的方式形成如下：

其中，L₁和L₂是分别按L1和L2信号的波长换算的载波相位测量，并且每个包括已添加的近似整周模糊度值，以使换算后的载波相位测量接近与对应的代码测量相同的值。因此，

其中，

和

是在同一测量时期分别在L1和L2频率上的原始载波相位测量，并且N₁和N₂的整周值已在载波相位跟踪开始时由用户或物体140初始化，以提供在对应代码测量的一个载波波长内的值，以使在换算后的载波相位测量和对应的代码测量之间保持小的差异。从方程式(7)的形式可注意到折射校正的载波相位测量包括具有由f₁与f₂之和(大约为2.803GHz)确定的波长λ的整周模糊度，因此，λ大约为0.1070米(即，c/(f₁+f₂)。

由于根据方程式(6)-(9)，已从代码和载波相位测量两者中消除电离层效应，并且有关伪距和载波相位测量的卫星时钟和轨道误差的效应是相同的，因此，在步骤310中获得的P_RC和L_RC的值应几乎相同，除与载波相位测量相关联的可能的整周模糊度和在代码测量中更高的多径噪声以外。这通过平滑一系列测量时期在折射校正的代码测量与折射校正的载波相位测量之间的偏移(O＝P_RC-L_RC)，使偏移变为“浮点模糊度”的更加准确的估计，从而允许在L_RC中整周模糊度的求解。通过使用后置(post-fix)测量残差，还可进一步调整平滑的偏移值，以提供附加的载波相位测量调整，使得调整后的测量残差接近零。

在本发明的一个实施例中，通过取如下的偏移的扩展平均值而平滑所述偏移：

$O_i=O_{i-1}+(P_{\mathrm{RC}}^i-L_{\mathrm{RC}}^i-O_{i-1})/\mathrm{\&eta;},---\left(10\right)$

其中，i＝1，2，3，...，用于标明测量时期，并且η值是置信值，在O_i变为浮点模糊度值的更准确估计时增大。在本发明的一个实施例中，η等于i，直至达到最大平均值。例如，如果载波相位测量假设为只具有代码测量噪声的1/100，则“η”的值会限于小于100的平方或10000。方程式(9因此可以递归计算，直至达到浮点模糊度值的预定准确度。

借助于平滑的偏移O_i，通过将当前测量时期的折射校正的载波相位测量加上平滑的偏移，可获得平滑的折射校正的代码测量S，以使

S_i＝O_i+L_i                    (11)

其具有载波相位测量的准确度，而无相关联的模糊度。

对在用户GPS接收机142的视野中的多个卫星中的每一卫星执行如结合方程式(6)-(11)描述的以上过程。利用可用于在用户GPS接收机142的视野中的多个卫星中的每一卫星的平滑的折射校正的代码测量，可获得到这些卫星的伪距。这些伪距通过从中心105接收的WADGPS校正进行调整，并且在加权最小二乘定位中用于计算状态向量x。这样，用户GPS接收机142的位置、速度和时间(PVT)可作为用户GPS接收机142的PVT的WADGPS解来计算。

在通过引用结合于本文的“代码和载波测量的协同作用”(by Hatch，R.in theProceedings of the Third International Geodetic Symposium on Satellite Doppler Positioning，DMA，NOS，Las Cruces，N.M.，New Mexico State University，Vol.II，pp.1213-1232)中以及在也通过引用结合于本文中的美国公布号US-2005-0024263的共同拥有的专利申请“用于生成广域或全球差分GPS系统的时钟校正的方法”中，可找到获得平滑的、折射校正的偏移的方法的其它例子。

也可能将“浮点模糊度”值作为在最小二乘或Kalman滤波解中单独的状态来求出。当模糊度被作为状态包括时，每个浮点模糊度值的估计值根据方差进行调整，以使它在系统几何由于卫星运动而更改时变得更加准确。因此，此技术还随时间进展而产生更加准确的估计。参阅Patrick H.C.Hwang的题为“用于差分定位的动态GPS：求解变化的整数模糊度”(in Navigation Vol.38，No.1，Spring 1991)的论文，该论文通过引用结合于本文。

上述技术有许多种组合和变化可用于佑计“浮点模糊度”值。然而，它们全部涉及在很大的时间间隔内处理数据。在可以确信“浮点模糊度”足够准确、可在用户140的导航位置中产生小于10厘米准确度前，此时间间隔经常可长达一或二小时。要缩短用于获得“浮点模糊度”值的时间间隔，可使用用户GPS接收机142的已知位置按如下所述初始化WADGPS系统。

图3A表示用于初始化WADGPS系统100的方法300。如图3所示，方法300包括步骤310，在该步骤中确定用户是否固定在一已知位置。这可根据用户输入或经允许计算机114确定用户接收机142是否已固定的某一常规机制而实现。如果用户接收机142已固定并且用户接收机142的位置已被准确获悉，则该位置可用于计算浮点模糊度值而无需把本地定位系统174和/或本地RTK系统150的帮助，如下文的详细描述。用户GPS接收机142的测量位置可用作已知位置，或者在一些环境中，可能只是由于用户接收机142已是固定的并且用户位置在先前操作期间已经确定，因而位置已知。

响应于用户固定在已知位置的确定，方法300进行到步骤320，在该步骤中将用户接收机位置设为已知位置。否则，方法300进行到步骤330，在该步骤中启用本地定位系统174和/或本地RTK系统150以使用上述方法来自动更新用户位置。

方法300还包括步骤340，在该步骤中，用户接收机位置无论是在步骤320或步骤330中确定，均用于计算到卫星110的一组理论距离。这可涉及到基于来自WADGPS系统100的广播星历表计算卫星110的位置和根据由WADGPS系统100广播的轨道校正调整那些位置。假设在笛卡尔(Cartesian)坐标中的用户接收机位置和卫星位置，从用户140到每个卫星110的理论距离可计算如下：

$r=\sqrt{{(x_s-x_u)}^2+{(y_s-y_u)}^2+{(z_s-z_u)}^2}---\left(11\right)$

其中，下标s标明卫星坐标，并且下标u标明用户或物体接收机坐标。

方法300还包括步骤350，在该步骤中，通过从计算的理论距离减去相对于同一卫星从折射校正的载波相位测量获得的距离而计算得出对应于每个卫星的初始浮点模糊度值a，因此，

$a=r-L_{\mathrm{RC}}^0---\left(12\right)$

其中，

表示在开始测量时期根据方程式(7)计算的折射校正的载波相位测量。

方法300还包括步骤360，在该步骤中，通过将初始浮点模糊度值加上在随后测量时期的对应的折射校正的载波相位测量，即

${L^i}_{\mathrm{RC}}={L^i}_{\mathrm{RC}}+a,$

并且通过将浮点模糊度值视为众所周知，使得置信度设为高(或方差设为低)，从而求解浮点模糊度值。实际上步骤360通过在用于确定浮点模糊度值的进程中使用小的增益值调整浮点模糊度值来完成。例如，如果通过根据方程式(9)平滑在折射校正的代码测量与折射校正的载波相位测量之间的偏移来确定浮点模糊度值，则小的增益意味着如同在计算浮点模糊度值时已使用大量偏移值处理浮点模糊度值，因此，η＝i+(大的数字)。如果模糊度值在Kalman滤波进程中确定，则通过将模糊度状态的方差设为小的值而实现小的增益。

因此，通过使用固定用户接收机142的已知位置，或者通过使用本地定位系统174和/或本地RTK系统150来初始化浮点模糊度值，避免了在用户接收机位置未知时求出浮点模糊度值所需的一般十五分钟到二小时的“捕捉”时间。这可大大加快WADGPS系统100中用于求解载波相位模糊度的进程，使WADGPS系统100更适用于实时定位和/或导航用途。

为了在方法300中使用本地定位系统174和/或本地RTK系统150来更新用户接收机位置，则本地定位系统174中的一个或多个地标176的位置和/或本地RTK系统150中的参考站120的位置可在WADGPS系统100中准确确定。可以在相对意义上使用常规本地定位系统或本地RTK系统，这表示用户接收机142的位置可相对于所述一个或多个地标和/或参考接收机确定。这样，即使所述一个或多个地标和/或参考接收机的坐标可能或可能不特别准确，并且除普通GPS数据外的坐标数据用于定位所述地标和/或所述参考站，也可获得用户GPS接收机142的准确相对位置。然而，对于本地定位系统174、本地RTK系统150和/或WADGPS系统100的组合使用，需要确定本地定位系统174中的一个或多个地标176和RTK系统150中的参考接收机120的绝对位置。如果将不正确的位置用于本地定位系统174中的一个或多个地标176或本地RTK系统150中的参考站120，则它将导致如上所述计算的浮点模糊度值不正确。当在随后的WADGPS处理期间浮点模糊度值缓慢调整到正确值时，这将导致计算的用户接收机142的位置缓慢漂移。

在本发明的一个实施例中，为了增加可靠性，本地定位系统174中的一个或多个地标176的平均位置和/或RTK系统150中的参考站120的平均位置基于来自WADGPS系统100的数小时的定位数据来确定。在一个替代的实施例中，在所述一个或多个地标176和/或参考站120的计算机系统接受其位置的操作员输入值，并将该位置提供给用户140。这允许使用参考站的该位置立即开始相对本地定位和/或RTK定位。同时，所述一个或多个地标176和/或参考站120的更准确位置由WADGSP系统100确定，并传送给所述一个或多个地标176和/或参考站120。此更准确的位置或操作员输入位置与由WADGPS系统100确定的所述一个或多个地标176和/或参考站120的更准确位置之间的偏移可以随后以较低的速率传送到用户140。

图3B更详细地表示方法300中的步骤330，在该步骤中，使用本地定位系统174和/或本地RTK系统150来更新用户位置。如图3B所示，步骤330包括子步骤331和子步骤333；在子步骤331中，用户或物体140接收本地定位系统174和/或本地RTK系统150中的参考站120的操作员输入位置；在子步骤333中，用户或物体140执行本地定位和/或本地RTK操作以确定相对于所述一个或多个地标176和/或所述参考站120位置的其自身位置。步骤330还包括子步骤335，在该步骤中，用户或物体140接收由WADGPS系统100确定的所述参考站120的更准确位置或在参考站120的操作员输入位置与由WADGPS系统100确定的参考站120的更准确位置之间的偏移。步骤330还包括子步骤337，在该步骤中，用户或物体140使用所述地标和/或所述参考站的用户输入位置或由WADGPS系统100确定的一个或多个地标176和/或参考站120的位置(如果可用)，计算用户GPS接收机142在笛卡尔坐标中的绝对位置。

通过使用方法300可获益的一个例子是定位火车。当火车通过隧道时，本地定位系统链路、本地RTK链路和全球WADGPS链路均会丢失。这种情况下，当火车穿出隧道时，可建立本地定位系统数据链路和/或RTK数据链路以初始化WADGPS浮点模糊度值。这会避免另外确定正确浮点模糊度值所需的长数据间隔。

通过使用方法300可获益的另一例子是定位起飞后的飞机。这种情况下，在飞机准备起飞的机场的本地定位系统和/或本地RTK系统可用于在起飞之前或期间初始化WADGPS模糊度。

因此，包括用户GPS接收机142和藕合至所述用户GPS接收机142的计算机系统144的用户或物体140可在本地定位模式、RTK模式和/或WADGPS模式两者中操作。由于如上所述用于本地定位系统174和本地RTK系统150的搜索进程比WADGPS系统100中用于求解整数模糊度值的平滑方法花费的时间少得多，因此，本地定位系统174和本地RTK系统150比WADGPS系统更有利。在搜索进程中，不需要代码测量的平滑或者执行持续时间短得多的代码测量的平滑，不是直接确定整周模糊度，而是在整数模糊度值的初始集合中提供降低的不确定性，使得随后的搜索进程可受到更紧密的约束。为此，只要几秒的数据便足以获得模糊度值的初始集合。

然而，只在用户GPS接收机142与本地定位系统174中的一个或多个地标176和/或本地RTK系统150中的参考站120之间的通信链路可维持，并且用户或物体140未偏离本地定位系统174中的一个或多个地标176和/或本地RTK系统150中的参考站120太远的情况下，本地定位系统174和/或本地RTK系统150才可用。当这些条件未满足时，即，当本地定位系统174和/或本地RTK系统150不可用或不准确时，通过使用由本地定位系统174和/或本地RTK系统150中最后确定的用户接收机位置来初始化WADGPS系统，用户可求助于WADGPS系统100进行导航，从而避免获得“浮点模糊度”值的长“捕捉(pull-in)”时间。

图4A表示由用户计算机系统144执行的组合本地定位、RTK和WADGPS操作的进程流400。该进程流包括步骤440、450和460。如图4A中所示，在本地定位校正可用时，用户140在本地定位模式操作，而在RTK校正可用时，用户140在RTK模式操作。用户140接收本地定位系统174和/或本地RTK系统150中的地标/参考站120的位置401，并执行步骤440，在该步骤440中，用户接收机的PVT使用从本地定位系统174中的地标176和/或本地RTK系统150中的参考接收机120接收的本地定位/RTK校正410来确定。在执行步骤440期间，用户140可继续从中心105接收WADGPS校正420，以便WADGPS解可在后台生成。用户140还可从中心105以较低的速率接收本地定位系统174中的地标176和/或本地RTK系统150中的参考站120的更新位置430。通过使用所述地标176和/或参考站120的更新位置和用户接收机位置的本地定位/RTK解，可根据上述的方法300在后台持续初始化WADGPS解以与所述本地定位/RTK解一致。

当所述本地定位和RTK校正丢失时，用户140切换到WADGPS操作模式并执行步骤450，在该步骤450中，用户140在本地定位/RTK校正变得不可用之前立即使用在所述本地定位和/或RTK操作模式中确定的用户接收机位置，以根据上述方法300初始化用于WADGPS操作模式的浮点模糊度值。这样，可确定“浮点模糊度”值而没有长的“捕捉”时间。在执行步骤450期间，用户140继续从中心105接收WADGPS校正420。用户140还可从中心105以较低的速率接收本地定位系统174中的一个或多个地标176和/或本地RTK系统150中的参考站120的更新位置430。参考站坐标用于将WADGPS模式中生成的用户接收机位置变换成相对于所述一个或多个地标176和/或本地参考接收机120的位置。这样，由用户计算机系统144生成的PVT结果将在这两种不同的操作模式之间无缝过渡。

当所述本地定位和/或RTK校正再次可用时，用户在步骤460恢复本地定位和/或RTK操作，这类似于步骤440中的所述本地定位和/或RTK操作。

图4B是表示使用本地定位系统、本地RTK系统和/或WADGPS系统的进程流470的流程图。如果可用，则用户的位置可根据来自本地定位系统(480)的信息确定。如果可用，则用户的位置可根据来自RTK系统(482)的信息确定。如果可用，则用户的位置可根据来自WADGPS系统(484)的信息确定。载波相位测量中的浮点模糊度值可被初始化(486)。所述进程流470可以包括更少或附加的操作。两个或更多个操作可以合并，并且至少一个操作的位置可以改变。

虽然基于接收自本地定位或本地RTK系统(例如来自本地参考站)的校正信号的导航典型地比基于接收自WADGPS系统的信息的导航更准确，但是多径信号可能不利地影响本地定位或本地RTK系统中参考站产生和传输的校正信息的准确度。因此，如果本地定位或本地RTK参考站产生的校正信息已被降级(通过多径信号或其它方式)超过一阈值量，则它将需要检测状态并基于来自WADGPS系统的信息导航物体或用户。下面参考图4C描述的进程通过在接收自本地定位或本地RTK参考站的校正信号不准确时进行检测，增加了多模式导航系统的可靠性。当状态被检测时，所述物体或用户的导航是基于接收自WADGPS系统的信息。如果接收自本地定位或本地RTK参考站的所述校正信号的准确度此后得到改善，则所述进程可以过渡到使用接收自本地定位或本地RTK系统的校正信息导航。

图4C是表示用于在从WADGPS系统过渡到本地定位/本地RTK系统时增加定位信息的可靠性的进程流490的流程图。如果可用，则可根据接收自WADGPS系统(492)的信息确定物体或用户的一位置。使用所述WADGPS系统的所述物体或用户的导航可以继续，直到与所述本地定位或本地RTK系统的参考站的通信被建立(或重新建立)。如果可用，则用户的位置校正信息可以接收自本地定位或本地RTK系统(494)，并且可以根据接收自本地定位或本地RTK系统的信息确定所述物体或用户的一位置(有时也称本地RTK位置)(496)。果两个位置(同一物体或用户的)的当前值可用，则比较使用来自所述WADGPS系统的信息确定的所述物体或用户的所述位置与使用来自所述本地定位/本地RTK系统的信息确定的所述物体或用户的所述位置(497)。如果基于来自所述WADGPS系统的信息确定的所述物体或用户的所述位置与基于来自所述本地定位/本地RTK系统的信息确定的所述物体或用户的所述位置相差大于一预定的阈值，则通过WADGPS系统确定的用户的位置被用于导航(498)。如果基于来自所述WADGPS系统的信息确定的所述物体或用户的所述位置与基于来自所述本地定位/本地RTK系统的信息确定的所述物体或用户的所述位置相差小于一预定的阈值，则基于来自所述本地定位/本地RTK系统的信息确定的所述物体或用户的所述位置被用于导航(499)。

在一些实施例中，所述预定的阈值是一个活多个固定值。例如，在一个实施例中，所述预定的阈值是一用于垂直位置差的第一值(如20cm)，及一用于水平位置差的第二值(如15cm)。如果所述WADGPS和本地RTK位置垂直距离相差大于所述第一值，或者水平距离相差大于所述第二值，则来自所述本地定位/RTK系统的所述位置校正信息被拒绝，而所述WADGPS位置被用于导航。在一些其它实施例中，利用物体的所述第一和第二位置(所述WADGPS位置和本地RTK位置)的协方差矩阵的卡方(chi-squared)比率检验法确定所述预定的阈值。协方差矩阵的卡方比率检验法的利用在以下文章中有表述：Satirapod C，Wang J和RizosC.的“用于GPS可观察量的方差-协方差的估算的简化的MINQUE方法(ASimplified MINQUE Procedure for the Estimation of Variance-Covariance Componentsof GPS Observables)”，Survey Review，36卷286期第582-590页，2002；H-K-Lee，C.Rizos和G.I.Jee的“具有一致误差协方差信息的动态DGPS滤波器的设计和分析(Design and Analysis of DGPS Filter with Consistent Error CovarianceInformation)，”卫星导航(SatNav)2003，第6届包括移动定位服务的卫星导航技术国际研讨会(6th International Symposium on Satellite Navigation TechnologyIncluding Mobile Positioning & Loc ation Services)，澳大利亚墨尔本，2003年7月22日；及H.K.Lee，C.Rizos和G.I.Jee的“具有一致误差协方差信息的动态DGPS滤波器的设计(Design of Kinematic DGPS Filters with Consistent Error CovarianceInformation)，”雷达、声纳和导航，IEEE会议文件，151卷6期第382-388页，2004年12月10日，各作为背景信息在此引入作为参考。

在一典型的实施例中，用户GPS接收机142可以在第一模式下操作，该模式下，当与本地定位系统176的通信可用时，使用本地定位系统176确定用户140的第一位置。所述用户140的第二位置可以根据使用在第二操作模式中的WADGPS系统100执行的载波相位测量确定。诸如所述第一位置的所述用户140的一已知位置可以被用于初始化所述载波相位测量中的浮点模糊度值。在一些实施例中，所述用户140的已知位置可以由所述用户提供和/或输入。

在一些实施例中，所述第一操作模式在可用时被用来确定用户140的位置。然而，如果与本地定位系统174的通信丢失，则可以使用所述第二操作模式。如果至本地定位系统174的距离超过一值，诸如100米、500米、1000米、10000米或更多，则与本地定位系统174的通信可能丢失。

在一些实施例中，所述第一操作模式和所述第二操作模式可以基本上同时执行，所述第一位置与所述第二位置之间的差被用于初始化所述载波相位测量中的浮点模糊度值。在一些实施例中，所述第一操作模式和所述第二操作模式可以基本上同时执行，所述第一位置与所述第二位置之间的差被用于确定所述用户140的第三位置。可以根据接收的来自第三操作模式中的实时动态(RTK)系统150中的本地参考接收机122的信息确定所述第三位置。

在一些实施例中，当与本地定位系统174的通信丢失时可使用第三操作模式，而当与所述本地定位系统174的通信再次可用时可以使用所述第一操作模式。

在一些实施例中，当与本地参考接收机122及本地定位系统174的通信丢失时可以使用所述第二操作模式，而当与所述本地定位系统174的通信可用时可以使用所述第一操作模式，其中当与本地参考接收机122的通信可用而与本地定位系统174的通信丢失时使用第三操作模式。

在一些实施例中，如果从本地定位系统174至用户140的距离大于第一值(如10000米)，则使用所述第二操作模式；而如果从本地定位系统174至用户140的距离小于第二值(如1000米)，则使用所述第一操作模式；如果从本地定位系统174至用户140的距离在所述第一预定值与所述第二预定值之间，则使用所述第三操作模式。

在一些实施例中，当所述第二操作模式确定的位置与所述第一操作模式确定的位置进行比较并且所述第二操作模式确定的所述位置与所述第一操作模式确定的所述位置相差大于一预定的阈值时，使用所述第二操作模式。

在一些实施例中，当所述第二操作模式确定的位置与所述第三操作模式确定的位置进行比较并且所述第二操作模式确定的所述位置与所述第三操作模式确定的所述位置相差大于一预定的阈值时，使用所述第二操作模式。

进程400可在许多应用中使用。一个应用涉及将本地定位系统和/或RTK操作扩展到本地定位系统和/或RTK无线电链路无法维持而WADGPS通信链路至少通常可用的区域中。例如，如图5中所示，用户或物体140可以是在起伏的坡面区域501中按行520移动的农用车辆510，而用户接收机142接到农用车辆或接到连接到农用车辆的农用装置。区域501包括从本地RTK系统150中的参考站120可见的区域503和从参考站120看不见的(阴影)区域505和507。由于RTK通信链路通常是直线对传，因此，只要将用户GPS接收机142从区域503移到区域505或507，RTK数据便会丢失。但由于经常得到卫星的帮助，因此，在用户接收机142与WADGPS系统100之间的数据链路通常是可用的。通过只要RTK无线电链路可用并且RTK系统150在操作便在WADGPS系统100中初始化浮点模糊度，实际上可在RTK链路丢失的那些间隔期间保持RTK操作的准确度。

虽然图1中的WADGPS系统100已在上述说明中使用，但将理解，利用来自卫星的载波相位测量实现定位和/或导航用途并因此需要确定与相位测量相关联的模糊度值的任何区域、广域或全球系统也可得益于上述方法300和进程400。这些系统的例子包括由John Deere公司开发的Starfire^TM system和由几家美国政府机构在开发的区域高准确度国家差分(HA-ND)GPS系统。

Claims (22)

1、一种利用广域差分卫星定位系统与本地定位系统二者定位或导航一物体的方法，所述方法包括：

根据来自所述广域差分卫星定位系统的信息确定所述物体的第一位置；

接收来自所述本地定位系统的一基台的位置校正信息，并根据接收自所述本地定位系统的所述基台的所述校正信息确定所述物体的第二位置；

比较所述第一位置与所述第二位置；及

当所述第一位置与所述第二位置相差大于一预定阈值时，利用所述第一位置以导航所述物体，而当所述第一位置与所述第二位置相差小于所述预定的阈值时，利用所述第二位置以导航所述物体。

2、根据权利要求1所述的方法，其特征在于：所述预定的阈值是利用一物体的所述第一和第二位置的协方差矩阵的卡方比率检验法确定的。

3、根据权利要求1所述的方法，其特征在于：所述预定的阈值是一相对于一预定的方向的15至20厘米之间的值。

4、根据权利要求1所述的方法，其特征在于：所述本地定位系统是一本地实时动态(RTK)定位系统。

5、根据权利要求1所述的方法，包括在检测到与所述本地定位系统的所述基台的通信丢失时，从利用所述本地定位系统导航所述物体过渡到利用所述广域差分卫星定位系统导航所述物体。

6、根据权利要求5所述的方法，包括在检测到与所述本地定位系统的所述基台的通信重新取得并确定所述第一位置与所述第二位置相差小于所述预定的阈值时，从利用所述广域差分卫星定位系统导航所述物体过渡到利用所述本地定位系统导航所述物体。

7、根据权利要求1所述的方法，其特征在于：确定所述第一位置包括接收来自所述广域差分卫星定位系统的位置校正信息。

8、一种计算机可读媒体，其存储有计算机可读程序指令，在由处理器执行所述指令时，所述指令使处理器执行一种利用广域差分卫星定位系统与本地定位系统二者定位或导航物体的方法，所述程序指令包括：

用于根据来自所述广域差分卫星定位系统的信息确定所述物体的第一位置的指令；

用于接收来自所述本地定位系统的一基台的位置校正信息并根据接收自所述本地定位系统的所述基台的所述校正信息确定所述物体的第二位置的指令；

用于比较所述第一位置与所述第二位置的指令；及

用于当所述第一位置与所述第二位置相差大于一预定阈值时利用所述第一位置以导航所述物体，而当所述第一位置与所述第二位置相差小于所述预定的阈值时利用所述第二位置以导航所述物体的指令。

9、根据权利要求8所述的计算机可读媒体，其特征在于：所述预定的阈值是利用一物体的所述第一和第二位置的协方差矩阵的卡方比率检验法确定的。

10、根据权利要求8所述的计算机可读媒体，其特征在于：所述预定的阈值是一相对于一预定的方向的15至20厘米之间的值。

11、根据权利要求8所述的计算机可读媒体，其特征在于：所述本地定位系统是一本地实时动态(RTK)定位系统。

12、根据权利要求8所述的计算机可读媒体，其特征在于：所述程序指令还包括用于在检测到与所述本地定位系统的所述基台的通信丢失时从利用所述本地定位系统导航所述物体过渡到利用所述广域差分卫星定位系统导航所述物体的指令。

13、根据权利要求12所述的计算机可读媒体，其特征在于：所述程序指令还包括用于在检测到与所述本地定位系统的所述基台的通信重新取得并确定所述第一位置与所述第二位置相差小于所述预定的阈值时从利用所述广域差分卫星定位系统导航所述物体过渡到利用所述本地定位系统导航所述物体的指令。

14、根据权利要求8所述的计算机可读媒体，包括用于从所述广域差分卫星定位系统接收位置校正信息的指令。

15、一种利用广域差分卫星定位系统与本地定位系统二者定位或导航物体的卫星导航接收机，包括：

一存储器；

一卫星信号接收机；

至少一个辅助接收机，用于接收来自广域差分卫星定位系统和所述本地定位系统的位置校正信息；

一处理器；及

至少一个程序，其存储在所述存储器中并由所述处理器执行，所述程序包括：

用于根据来自所述广域差分卫星定位系统的信息确定所述物体的第一位置的指令；

用于接收来自所述本地定位系统的一基台的位置校正信息并根据接收自所述本地定位系统的所述基台的所述校正信息确定所述物体的第二位置的指令；

用于比较所述第一位置与所述第二位置的指令；及

用于当所述第一位置与所述第二位置相差大于一预定阈值时利用所述第一位置以导航所述物体，而当所述第一位置与所述第二位置相差小于所述预定的阈值时利用所述第二位置以导航所述物体的指令。

16、根据权利要求15所述的卫星导航接收机，其特征在于：所述预定的阈值是利用一物体的所述第一和第二位置的协方差矩阵的卡方比率检验法确定的。

17、根据权利要求15所述的卫星导航接收机，其特征在于：所述预定的阈值是一相对于一预定的方向的15至20厘米之间的值。

18、根据权利要求15所述的卫星导航接收机，其特征在于：所述本地定位系统是一本地实时动态(RTK)定位系统。

19、根据权利要求15所述的卫星导航接收机，其特征在于：所述接收机在检测到与所述本地定位系统的所述基台的通信丢失时从利用所述本地定位系统导航所述物体过渡到利用所述广域差分卫星定位系统导航所述物体。

20、根据权利要求19所述的卫星导航接收机，其特征在于：所述接收机在检测到与所述本地定位系统的所述基台的通信重新取得并确定所述第一位置与所述第二位置相差小于所述预定的阈值时从利用所述广域差分卫星定位系统导航所述物体过渡到利用所述本地定位系统导航所述物体。

21、根据权利要求15所述的卫星导航接收机，包括用于从所述广域差分卫星定位系统接收位置校正信息的指令。

22、一种利用广域差分卫星定位系统与本地定位系统二者定位或导航物体的卫星导航接收机，其特征在于，所述卫星导航接收机：

根据来自所述广域差分卫星定位系统的信息确定所述物体的第一位置；

接收来自所述本地定位系统的一基台的位置校正信息，并根据接收自所述本地定位系统的所述基台的所述校正信息确定所述物体的第二位置；

比较所述第一位置与所述第二位置；及

当所述第一位置与所述第二位置相差大于一预定阈值时利用所述第一位置以导航所述物体，而当所述第一位置与所述第二位置相差小于所述预定的阈值时利用所述第二位置以导航所述物体。

Images