一种GPSRTK中测量隐蔽点点位坐标的方法

摘要

本发明公开了一种GPS RTK中测量隐蔽点点位坐标的方法：安置辅助测量装置；所述辅助测量装置包括三根可伸缩的测杆，该三根测杆分别通过铰链首尾相接组成三角形；将辅助测量装置的一个顶点安置在待测的隐蔽点，将另外两个顶点分别置于两个测量点，该两个测量点是可通视点；保持该三角形的边都处于水平状态；测量点数据采集；信息记录；内业数据处理：解三角形；计算AB的方位角αAB：计算AO方位角αAO：求隐蔽点O点的坐标。该方法基于在隐蔽点附近建立局部极坐标系，并通过建立隐蔽点与测量点间的计算关系，测量获取局部极坐标系参数，同时测量隐蔽点在局部极坐标间的极坐标，从而快速确定隐蔽点点位坐标。

权利要求

1.一种GPS RTK中测量隐蔽点点位坐标的方法，其特征在于，具体包括如下步骤：步骤1：外业数据采集

1）安置辅助测量装置

所述辅助测量装置包括三根可伸缩的测杆，该三根测杆分别通过铰链首尾相接组成三角形，测杆的长度可调整范围为1.5m至2.5m；

依据现场具体情况调整辅助测量装置的各边长；将辅助测量装置的一个顶点安置在待测的隐蔽点O点，将另外两个顶点分别置于两个测量点A、B，测量点A、B均为可通视点；保持每根测杆都处于水平状态；

2）测量点数据采集

利用GPS RTK流动站获取测量点A、B点的GPS RTK坐标；

3）信息记录

需要记录的数据：隐蔽点点号NoO；测量点点号NoA；A点X坐标XA；A点Y坐标YA；A点高程值hA；测量点点号NoB；B点X坐标XB；B点Y坐标YB；B点高程值hB；辅助测量装置的边长SOA、SAB、SBO；

步骤2：内业数据处理：

1）解三角形

已知△ABO的边长SOA、SAB和SBO，利用余弦定理分别求出∠AOB，∠OAB，∠OBA；

2）计算AB的方位角αAB；

3）计算AO方位角αAO：

αAO＝αAB-∠OAB；

4）求隐蔽点O点的坐标

由A、O两点边长SAO和坐标方位角αAO计算坐标增量ΔxAO和ΔyAO：ΔxAO＝SAO cosαAOΔyAO＝SAO sinαAO求得隐蔽点O点坐标为：

XO=XA+ΔxAO；YO=YA+ΔyAO；hO=（hA+hB）/2。

2.如权利要求1所述的方法，其特征在于，每根测杆上均安装有水准管（3）。

3.如权利要求1所述的方法，其特征在于，每根测杆由两根重叠相接的木质的量尺（1）组成，该两根量尺（1）通过螺栓（2）连接且该两根量尺（1）能够相对运动调整测杆的长度。

4.如权利要求1所述的方法，其特征在于，所述步骤2中的2）所述计算AB的方位角αAB：首先，分别计算X方向和Y方向上的坐标增量ΔxAB和ΔyAB：ΔxAB=XB-XA

ΔyAB=YB-YA

其次，计算坐标增量的反正切值：αAB锐最后，通过ΔxAB、ΔyAB的正负号来判断αAB所在的象限，得到AB的方位角αAB。

说明书

一种GPS RTK中测量隐蔽点点位坐标的方法

技术领域

[0001] 本发明属于工程测量领域，具体涉及一种GPS RTK中测量隐蔽点点位坐标的方法。

背景技术

[0002] 一、GPS-RTK测量背景简介

[0003] 全球定位系统(Global Positioning System，通常简称GPS)，又称全球卫星定位系统，是一个中距离圆型轨道卫星导航系统。它可以为地球表面绝大部分地区(98％)提供准确的定位、测速和高精度的时间标准。系统由美国国防部研制和维护，可满足位于全球任何地方或近地空间的军事用户连续精确的确定三维位置、三维运动和时间的需要。该系统包括太空中的24颗GPS卫星；地面上的1个主控站、3个数据注入站和5个监测站及作为用户端的GPS接收机。通过测量卫星，就能迅速确定用户端在地球上所处的位置及海拔高度；所能收联接到的卫星数越多，解码出来的位置就越精确。

[0004] 实时动态差分测量(Real Time Kkinematic，通常简称RTK)。这是一种常用的基于GPS技术的测量方法，高精度的GPS测量必须采用载波相位观测值，GPS RTK定位技术就是基于载波相位观测值的实时动态定位技术，它能够实时地提供测站点在指定坐标系中的三维定位结果，并达到厘米级精度。在GPS RTK作业模式下，基准站通过数据链将其观测值和测站坐标信息一起传送给流动站。流动站不仅通过数据链接收来自基准站的数据，还要采集GPS观测数据，并在系统内组成差分观测值进行实时处理，同时给出厘米级定位结果，历时不足一秒钟。流动站可处于静止状态，也可处于运动状态；可在固定点上先进行初始化后再进入动态作业，也可在动态条件下直接开机，并在动态环境下完成周模糊度的搜索求解。在整周末知数解固定后，即可进行每个历元的实时处理，只要能保持四颗以上卫星相位观测值的跟踪和必要的几何图形，则流动站可随时给出厘米级定位结果。

[0005] 常规GPS静态、快速静态、动态测量都需要事后进行解算才能获得厘米级的精度，而GPSRTK是能够在野外实时得到厘米级定位精度的测量方法，它采用了载波相位动态实时差分方法，它的出现为工程放样、地形测图，各种控制测量带来了新途径，极大地提高了外业作业效率。

[0006] 二、全球定位系统实时动态差分测量(GPS RTK)的局限性

[0007] GPS RTK依靠的是接收从地面以上约两万公里的卫星发射来的无线电信号，相对而言，这些信号频率高、功率低，不易穿透可能阻挡卫星和GPS接收机之间视线的障碍物。事实上，存在于GPS接收机和卫星之间路径上的任何物体都会对系统的操作产生有害影响。有些物体如房屋，会完全屏蔽卫星信号，有些物体如树木会部分阻挡、反射或折射信号，GPS信号的接收在树林茂密的地区会很差。树林中有时会有足够的信号来计算概略位置，但信号清晰度难以达到厘米水平的精确定位。

[0008] GPS RTK要求能观测到足够多的卫星来精确可靠地实现定位，一般情况下，要求天空中至少有5颗适当分布的卫星，就可作精确可靠的定位。这就要求测量点附近留有足够的开放空间，使GPS RTK系统可观测到至少5颗卫星，RTK测量才有能够实施的条件。事实上，在林间或城市建筑群中测量时，这一条件不易实现，在GPS RTK测量中，将不能达到GPS RTK观测条件的点称之为隐蔽点，隐蔽点一般位于高层建筑物L形拐角处、大树底下等GPS信号不易覆盖处，隐蔽点的点位坐标是难以被测量到的。

[0009] 因此，研究一种全球定位系统实时动态差分测量中测量隐蔽点坐标的方法有着重要的现实意义。

发明内容

[0010] 针对现有技术存在的缺陷或不足，本发明的目的在于，提供一种全球定位系统实时动态差分测量中测量隐蔽点坐标的方法，该方法基于在隐蔽点附近建立局部极坐标系，并通过建立隐蔽点与测量点间的运算关系，测量获取局部极坐标系参数，同时测量隐蔽点在局部极坐标间的极坐标，从而快速确定隐蔽点点位坐标。

[0011] 为了达到上述任务，本发明采用如下技术解决方案：

[0012] 一种GPS RTK中测量隐蔽点点位坐标的方法，其特征在于，具体包括如下步骤：

[0013] 步骤1：外业数据采集

[0014] 1)安置辅助测量装置

[0015] 所述辅助测量装置包括三根可伸缩的测杆，该三根测杆分别通过铰链首尾相接组成三角形，测杆的长度可调整范围为1.5m至2.5m；

[0016] 依据现场具体情况调整辅助测量装置的各边长；将辅助测量装置的一个顶点安置在待测的隐蔽点O点，将另外两个顶点分别置于两个测量点A、B，测量点A、B均为可通视点；保持每根测杆都处于水平状态；

[0017] 2)测量点数据采集

[0018] 利用GPS RTK流动站获取测量点A、B点的GPS RTK坐标；

[0019] 3)信息记录

[0020] 需要记录的数据：隐蔽点点号NoO；测量点点号NoA；A点X坐标XA；A点Y坐标YA；A点高程值hA；测量点点号NoB；B点X坐标XB；B点Y坐标YB；B点高程值hB；辅助测量装置的边长SOA、SAB、SBO；

[0021] 步骤2：内业数据处理：

[0022] 1)解三角形

[0023] 已知△ABO的边长SOA、SAB和SBO，利用余弦定理分别求出∠AOB，∠OAB，∠OBA；

[0024] 2)计算AB的方位角αAB；

[0025] 3)计算AO方位角αAB：

[0026] αAO＝αAB-∠OAB；

[0027] 4)求隐蔽点O点的坐标

[0028] 由A、O两点边长SAO和坐标方位角αAO计算坐标增量ΔxAO和ΔyAO：

[0029] ΔxAO＝SAO cosαAO

[0030] ΔyAO＝SAO sinαAO

[0031] 求得隐蔽点O点坐标为：

[0032] XO＝XA+ΔXAO；YO＝YA+ΔYAO；hO＝(hA+hB)/2。

[0033] 进一步的，每根测杆上均安装有水准管。

[0034] 进一步的，每根测杆由两根重叠相接的木质的量尺组成，该两根量尺通过螺栓连接且两者能够相对运动调整测杆的长度。

[0035] 进一步的，所述步骤2中的2)所述计算AB的方位角αAB：

[0036] 首先，分别计算X方向和Y方向上的坐标增量ΔxAB和ΔyAB：

[0037] ΔxAB＝xB-xA

[0038] ΔyAB＝yB-yA

[0039] 其次，计算坐标增量的反正切值：αAB锐

[0040]

[0041] 最后，通过ΔxAB、ΔyAB的正负号来判断αAB所在的象限，得到AB的方位角αAB。

[0042] 本发明利用GPS RTK以及辅助测量装置来解决隐蔽点测量问题，具有适用性广、成本低、操作简便、理论模型便于计算机编程实现自动化处理等优点。

附图说明

[0043] 图1是本发明采用的辅助测量装置的结构示意图。

[0044] 图2是隐蔽点与测量点量测关系示意图。图中N方向为正北方向。

[0045] 以下结合附图和具体实施方式对本发明进一步解释说明。

具体实施方式

[0046] 本发明设隐蔽点为O点，测量点(通视点)为A、B点，隐蔽点O点坐标为(XO，YO，hO)，A、B点坐标为(XA，YA，hA)、(XB，YB，hB)，测量点A点和B点坐标已知，O点坐标未知。本发明基于可量测几何关系，建立隐蔽点与测量点间的计算关系，在具体作业过程中，通过辅助测量装置，确立隐蔽点与测量点间的具体计算公式，从而快速确定隐蔽点点位坐标。

[0047] 本发明的GPS RTK中测量隐蔽点点位坐标的方法，具体包括如下步骤：

[0048] 1、制作辅助测量装置

[0049] 如图1所示，辅助测量装置包括三根可伸缩的测杆，该三根测杆分别通过铰链首尾相接组成三角形，测杆的长度可调整范围为1.5m至2.5m。每根测杆上均安装有水准管3。本实施例中，每根测杆由两根重叠相接的木质的量尺1组成，该两根量尺1通过螺栓2连接且两者能够相对运动调整测杆的长度，每根量尺量程为1.5m，测杆的长度可调整范围为1.5m至2.5m。

[0050] 2、外业数据采集

[0051] 1)安置辅助测量装置

[0052] 依据现场具体情况调整辅助测量装置的各边长；将辅助测量装置的三角形的一个顶点安置在待测的隐蔽点O点，将另外两个顶点A、B点分别置于两个测量点，该两个测量点是可通视点；观察每根测杆上安装的水准管3，保持每根测杆都处于水平状态。为计算方便，O、A、B三个点的摆放按照固定顺序，本实施例以逆时针顺序摆放。

[0053] 2)测量点数据采集

[0054] 将GPS RTK流动站分别置于辅助测量装置的两个顶点位置，获取测量点A、B的GPS RTK坐标。如果由于卫星信号问题不能获取A、B点坐标值，可适当改变辅助测量装置的安放位置，或者重新调整其三边尺寸，直至获取测量点A、B的坐标。

[0055] 3)信息记录

[0056] 需要记录的数据：隐蔽点点号NoO；测量点点号NoA；A点X坐标XA；A点Y坐标YA；A点高程值hA；测量点点号NoB；B点X坐标XB；B点Y坐标YB；B点高程值hB；辅助测量装置的边长SOA、SAB、SBO。

[0057] 3、参照图2，内业数据处理如下：

[0058] 1)解三角形

[0059] 已知△ABO的边长SOA、SAB和SBO，则△ABO可解，利用余弦定理分别求出∠AOB，∠OAB，∠OBA；

[0060] 2)计算AB的方位角αAB

[0061] 方位角定义：由北方向沿顺时针方向量测到测线上的夹角，称为该直线的方位角，其范围为0°～360°。AB的方位角αAB的具体计算方法如下：

[0062] 首先，分别计算X方向和Y方向上的坐标增量ΔxAB和ΔyAB：

[0063] ΔxAB＝xB-xA

[0064] ΔyAB＝yB-yA

[0065] 其次，计算坐标增量的反正切值：αAB锐

[0066]

[0067] 最后通过ΔxAB、ΔyAB的正负号来判断αAB所在的象限，得到AB的方位角αAB：

[0068] a)ΔxAB＞0且ΔyAB＞0则为一象限，αAB＝αAB锐；

[0069] b)ΔxAB＜0且ΔyAB＞0则为二象限，αAB＝180°-αAB锐；

[0070] c)ΔxAB＜0且ΔyAB＜0则为三象限，αAB＝180°+αAB锐；

[0071] d)ΔxAB＞0且ΔyAB＜0则为四象限，αAB＝360°-αAB锐；

[0072] e)ΔyAB＝0且ΔyAB＞0则αAB＝0°；

[0073] f)ΔxAB＝0且ΔyAB＞0则αAB＝90°；

[0074] g)ΔyAB＝0且ΔxAB＜0则αAB＝180°；

[0075] h)ΔxAB＝0且ΔyAB＜0则αAB＝270°；

[0076] 3)计算AO方位角αAO：

[0077] αAO＝αAB-∠OAB；

[0078] 4)求隐蔽点坐标

[0079] 由A、O两点边长SAO和坐标方位角αAO计算坐标增量ΔxAO和ΔyAO：

[0080] ΔxAO＝SAO cosαAO

[0081] ΔyAO＝SAO sinαAO

[0082] 求得隐蔽点O点坐标为：

[0083] XO＝XA+ΔXAO；YO＝YA+ΔYAO；

[0084] 因测量时三角形保持水平状态，所以O点高程值hO＝hA＝hB，实际操作中一般取hO＝(hA+hB)/2。

[0085] 具体实施例：

[0086] 1、外业数据采集

[0087] 1)安置辅助测量装置

[0088] 设O点为隐蔽点，首先将辅助测量装置的一个顶点放置于O点处，另外两个顶点A、B两点分别置于观测条件较好的位置，将辅助测量装置的三边边长调整为2米；

[0089] 2)测量点数据采集

[0090] 利用GPS RTK流动站分别测量得到A、B两点的点位坐标(XA，YA，hA)、(XB，YB，hB)；

[0091] 3)信息记录

[0092] 记 录 数 据 为：O；A；3766214.962；611671.862；450.125；B；3766213.529；611673.257；450.127；2、2、2。

[0093] 2、内业数据处理

[0094] 1)解三角形

[0095] 已知△ABO的边长均为2m，则△ABO为等边三角形，∠AOB，∠OAB，∠OBA均为60°。

[0096] 2)计算AB方位角：αAB＝135°45′32″；

[0097] 3)计算AO方位角：αAO＝αAB-∠OAB＝75°45′32″；

[0098] 4)求隐蔽点坐标：

[0099] 由A、O两点边长SAO和坐标方位角αAO计算坐标增量ΔxAO和ΔyAO：

[0100] ΔxAO＝SAO cosαAO＝0.492；

[0101] ΔyAO＝SAO sinαAO＝1.939；

[0102] 求得隐蔽点O点坐标为：

[0103] XO＝XA+ΔXAO＝3766214.962+0.492＝3766215.454；

[0104] YO＝YA+ΔYAO＝611671.862+1.939＝611673.801；

[0105] hO＝(hA+hB)/2＝(450.125+450.127)/2＝450.126。