CN106595583A - 一种rtk测量接收机倾斜测量方法 - Google Patents

Abstract

本发明提供了一种RTK测量接收机倾斜测量方法，包括以下步骤：步骤S1：将对中杆底部固定并进行倾斜摇晃然后进行测量；步骤S2：基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度；步骤S3：基于测量得到的数值得到定位质量标的和待测点大地坐标；步骤S4：判断定位质量标的是否满足要求，决定是否结束测量，本发明根据一定数量倾斜测量时接收机天线相位中心的位置、倾角以及对中杆长度等计算待测点位置，对于厘米级的RTK，当观测值数量足够时，可以获得厘米级的位置计算精度。

Description

一种RTK测量接收机倾斜测量方法

技术领域

本发明涉及RTK测量技术领域，具体涉及一种RTK测量接收机倾斜测量方法。

背景技术

在RTK测量中，流动站对中杆需要保证竖直安置，然后根据GNSS接收机天线相位中心的位置及相位中心至地面的距离归算得到地面上待测点的位置坐标。但是在实际测量工作中，由于受到地形条件的影响和仪器本身的限制，许多测量点无法竖直安置，需要倾斜。导致了测量过程繁琐，测量误差大。

发明内容

为了解决上述不足的缺陷，本发明提供了一种RTK测量接收机倾斜测量方法，根据一定数量倾斜测量时接收机天线相位中心的位置、倾角以及对中杆长度等计算待测点位置，对于厘米级的RTK，当观测值数量足够时，可以获得厘米级的位置计算精度。

本发明提供了一种RTK测量接收机倾斜测量方法，包括以下步骤：

步骤S1：对中杆底部固定并进行倾斜摇晃然后进行测量；

步骤S2：基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度；

步骤S3：基于测量得到的数值得到定位质量标的和待测点大地坐标；

步骤S4：判断定位质量标的是否满足要求，决定是否结束测量。

上述的方法，其中，所述对中杆底部固定并进行倾斜摇晃然后进行测量的步骤包括：

测量中对中杆在使用过程中的长度保持不变，以及对中杆底部保持固定。

上述的方法，其中，所述基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度的步骤包括：

观测得到观测点的大地坐标(B_i,L_i,H_i)与测量倾角θ_i；

将观测点大地坐标(B_i,L_i,H_i)转换成平面坐标(N_i,E_i,H_i)；

基于上述的平面坐标(N_i,E_i,H_i)得到地面待测点平面坐标并转换成大地坐标

上述的方法，其中，所述基于测量得到的数值得到定位质量标的和待测点大地坐标的步骤包括：

通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程；

根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程。

上述的方法，其中，所述判断定位质量标的是否满足要求，决定是否结束测量的步骤包括：

判断定位质量标的是否满足要求，若不满足要求，则继续步骤S1到步骤S3；若满足要求，则结束测量。

上述的方法，其中，所述基于上述的平面坐标(N_i,E_i,H_i)得到地面待测点平面坐标并转换成大地坐标的步骤包括：

建立观测方程，间接平差；

根据平差后的观测值残差进行抗差估计，求解地面待测点平面坐标

上述的方法，其中，所述通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程的步骤中，设GPS接收机天线相位中心平面坐标(N_i,E_i)及大地高H_i(其中i＝1,2,...,n，n为观测点总数，下同)，待测点O的平面坐标及大地高分别为和待测点至观测点长度为则观测方程为

线性化后可得

其中为长度改正数；(v_N,v_E,v_H)为待测点坐标改正数；(N₀,E₀,H₀)为待测点近似坐标；L_i为对中杆长度与天线高之和；为观测点坐标与待测点近似坐标计算的距离。

上述的方法，其中，所述根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程的步骤中，设观测点观测倾角为则观测方程为

其中整理后可得

v_L'＝v_H-(H_i-L_icosθ_i-H₀) (4)

联立式(2)、(4)，可建立误差方程组

V＝BX-W,P (5)

式中，P为权矩阵，根据观测值的精度设定，若没有此信息，可设置为单位阵。

上述的方法，其中，包括以下步骤：根据加权最小二乘准则，可得

最后通过循环迭代计算出地面待测点的平面位置与大地高令

A为系数矩阵，计算ATA得到维度为3行3列的单位矩阵，将单位阵的对角线元素相加后取平方根得到质量标的。

本发明具有以下优点：根据一定数量倾斜测量时接收机天线相位中心的位置、倾角以及对中杆长度等计算待测点位置，对于厘米级的RTK，当观测值数量足够时，可以获得厘米级的位置计算精度。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明及其特征、外形和优点将会变得更明显。在全部附图中相同的标记指示相同的部分。并未刻意按照比例绘制附图，重点在于示出本发明的主旨。

图1为本发明的一种RTK测量接收机倾斜测量方法的流程示意图。

图2为本发明的一种RTK测量接收机倾斜测量方法的原理图。

具体实施方式

在下文的描述中，给出了大量具体的细节以便提供对本发明更为彻底的理解。然而，对于本领域技术人员而言显而易见的是，本发明可以无需一个或多个这些细节而得以实施。在其他的例子中，为了避免与本发明发生混淆，对于本领域公知的一些技术特征未进行描述。

为了彻底理解本发明，将在下列的描述中提出详细的步骤以及详细的结构，以便阐释本发明的技术方案。本发明的较佳实施例详细描述如下，然而除了这些详细描述外，本发明还可以具有其他实施方式。

参照图1-图2所示，本发明提供了一种RTK测量接收机倾斜测量方法，包括以下步骤：

步骤S1：对中杆底部固定并进行倾斜摇晃然后进行测量，测量中对中杆在使用过程中的长度保持不变，以及对中杆底部保持固定，也就是说对测量时接收机的稳定性没有限制，在底部固定的前提下可以随意摇晃。

步骤S2：基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度，具体为如图1中所示的A、B、C，其中A为测量点序列，根据时间顺序记录从开始测量到结束测量期间所有的观测点。其中每一个点都具有其定位信息，信息为实测点(天线相位中心)大地坐标，按照顺序包括大地维度B、大地经度L和大地高H。B、L的单位弧度，H的单位为米。B为测量点倾角序列，根据时间顺序记录了从开始测量到结束测量期间所有的倾角，倾角和测量点是一一对应的关系，表示了得到测量点坐标的同时一起的倾斜角度，单位为弧度。C为测量杆长度和天线相位中心高度，单位为米。D为定位质量标的，称之为几何因子。几何因子是一个正整数，最小值为1，数值越小表示计算的结果精度越高。E为待测点(地面点)大地坐标，为程序计算结果，包括大地纬度B、大地经度L和大地高度H。B、L的单位为弧度，H的单位为米。

具体包括步骤S2a：观测得到观测点的大地坐标(B_i,L_i,H_i)与测量倾角θ_i，步骤S2b：将观测点大地坐标(B_i,L_i,H_i)转换成平面坐标(N_i,E_i,H_i)，步骤S2c：基于上述的平面坐标(N_i,E_i,H_i)得到地面待测点平面坐标并转换成大地坐标具体为建立观测方程，间接平差；根据平差后的观测值残差进行抗差估计，求解地面待测点平面坐标

步骤S3：基于测量得到的数值得到定位质量标的和待测点大地坐标。

具体包括步骤S3a：通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程，具体为设GPS接收机天线相位中心平面坐标(N_i,E_i)及大地高H_i(其中i＝1,2,...,n，n为观测点总数，下同)，待测点O的平面坐标及大地高分别为和待测点至观测点长度为则观测方程为

线性化后可得

其中为长度改正数；(v_N,v_E,v_H)为待测点坐标改正数；(N₀,E₀,H₀)为待测点近似坐标；L_i为对中杆长度与天线高之和；为观测点坐标与待测点近似坐标计算的距离。

步骤S3b：根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程：

具体为设观测点观测倾角为则观测方程为

其中整理后可得

v_L'＝v_H-(H_i-L_icosθ_i-H₀) (10)

联立式(2)、(4)，可建立误差方程组

V＝BX-W,P (11)

式中，P为权矩阵，根据观测值的精度设定，若没有此信息，可设置为单位阵，根据加权最小二乘准则，可得

最后通过循环迭代计算出地面待测点的平面位置与大地高令

A为系数矩阵，计算ATA得到维度为3行3列的单位矩阵，将单位阵的对角线元素相加后取平方根得到质量标的。

步骤S4：判断定位质量标的是否满足要求，决定是否结束测量，具体为判断定位质量标的是否满足要求，若不满足要求，则继续步骤S1到步骤S3；若满足要求，则结束测量。

本发明的工作过程为：参照图1和图2所示，解算步骤如下：

(1)观测得到所有观测点的大地坐标(B_i,L_i,H_i)与测量倾角θ_i；

(2)将观测点大地坐标(B_i,L_i,H_i)转换成平面坐标(N_i,E_i,H_i)；

(3)建立观测方程，间接平差；

(4)根据平差后的观测值残差进行抗差估计，求解地面待测点平面坐标

(5)将待测点平面坐标转换成大地坐标

(6)输出转换后的大地坐标和质量标的(称之为“几何因子”)。

(7)根据质量标的的数值判断精度是否达标，如果达标结束测量，输出转换后的大地坐标。如果未达标则返回过程(1)继续测量，此次观测数据增加到(1)所述的所有观测点中。

步骤(6)中的质量标的(称之为“几何因子”)，是一个最小值为1的正整数，数值越小精度越高。

测量中对中杆在使用过程中的长度为固定不能更改，对中杆底部不能移动。测量时需要摇晃接收机，摇晃的幅度越大效果越好。

参照图2，其解算原理为：，其中O为地面待测点；A为GPS接收机天线相位中心，是倾斜测量的观测点；AB为天线高；OB为对中杆长度；OD为待测点处的参考椭球面法线；AC为观测点A至参考椭球面的法线，在投影后的平面坐标系中AC与OD平行；θ为仪器倾角。

倾斜测量主要包括观测点的位置和倾角两类观测值，分别对应两类观测方程。

通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程。设GPS接收机天线相位中心平面坐标(N_i,E_i)及大地高H_i(其中i＝1,2,...,n，n为观测点总数，下同)，待测点O的平面坐标及大地高分别为和待测点至观测点长度为则观测方程为

线性化后可得

其中为长度改正数；(v_N,v_E,v_H)为待测点坐标改正数；(N₀,E₀,H₀)为待测点近似坐标；L_i为对中杆长度与天线高之和；为观测点坐标与待测点近似坐标计算的距离。

根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程。设观测点观测倾角为则观测方程为

其中整理后可得

v_L'＝v_H-(H_i-L_icosθ_i-H₀) (16)

联立式(2)、(4)，可建立误差方程组

V＝BX-W,P (17)

式中，P为权矩阵，根据观测值的精度设定，若没有此信息，可设置为单位阵。

根据加权最小二乘准则，可得

最后通过循环迭代计算出地面待测点的平面位置与大地高令

A为系数矩阵，计算ATA得到维度为3行3列的单位矩阵，将单位阵的对角线元素相加后取平方根得到质量标的，称之为“几何因子”。几何因子是用来作为精度指标。几何因子越小精度越高。

通过本发明的方法可以实现：

1、只需要带倾角信息的RTK测量接收机即可，不需要额外的硬件成本。

3、与传统的通过测量接收机姿态角信息的倾斜测量模块相比，不受任何外部环境(磁干扰、温度变化等)影响。

4、与传统的倾斜测量方法相比，不需要进行磁场校准。

5、本发明所述的倾斜测量方法与传统的倾斜测量方法相比，对中杆底部位置的测量精度取决于观测点的几何构型健壮性，观测点数量和精度，倾角测量精度，对中杆长度误差和天线相位中心长度误差。测量结果更稳定、测量精度更高。

6、与传统的倾斜测量方法相比，对测量时接收机的稳定性没有限制，在底部固定的前提下可以随意摇晃。

7、与传统倾斜测量方法相比，可以对测量精度阈值进行设置，在精度满足要求的情况下，提高几何因子阈值可以缩短测量时间。

以上对本发明的较佳实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，其中未尽详细描述的设备和结构应该理解为用本领域中的普通方式予以实施；任何熟悉本领域的技术人员，在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例，这并不影响本发明的实质内容。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。

Claims (9)

1.一种RTK测量接收机倾斜测量方法，其特征在于，包括以下步骤：

步骤S1：将对中杆底部固定并进行倾斜摇晃然后进行测量；

步骤S2：基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度；

步骤S3：基于测量得到的数值得到定位质量标的和待测点大地坐标；

步骤S4：判断定位质量标的是否满足要求，决定是否结束测量。

2.如权利要求1所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述将对中杆底部固定并进行倾斜摇晃然后进行测量的步骤包括：

测量中对中杆在使用过程中的长度保持不变，以及对中杆底部保持固定。

3.如权利要求1所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述基于上述测量得到测量点序列、测量点倾斜序列、测量杆长度和天线相位中心高度的步骤包括：

观测得到观测点的大地坐标(B_i,L_i,H_i)与测量倾角θ_i；

将观测点大地坐标(B_i,L_i,H_i)转换成平面坐标(N_i,E_i,H_i)；

基于上述的平面坐标(N_i,E_i,H_i)得到地面待测点平面坐标并转换成大地坐标

4.如权利要求1-3任一项所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述基于测量得到的数值得到定位质量标的和待测点大地坐标的步骤包括：

通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程；

根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程。

5.如权利要求4所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述判断定位质量标的是否满足要求，决定是否结束测量的步骤包括：

判断定位质量标的是否满足要求，若不满足要求，则继续步骤S1到步骤S3；若满足要求，则结束测量。

6.如权利要求3所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述基于上述的平面坐标(N_i,E_i,H_i)得到地面待测点平面坐标并转换成大地坐标的步骤包括：

建立观测方程，间接平差；

根据平差后的观测值残差进行抗差估计，求解地面待测点平面坐标

7.如权利要求4所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述通过观测点位置坐标与观测点至待测点间的距离，根据后方交会原理建立第一类观测方程的步骤中，设GPS接收机天线相位中心平面坐标(N_i,E_i)及大地高H_i(其中i＝1,2,...,n，n为观测点总数，下同)，待测点O的平面坐标及大地高分别为和待测点至观测点长度为则观测方程为

${\tilde{L}}_i=\sqrt{{(N_i-{\tilde{N}}_0)}^2+{(E_i-{\tilde{E}}_0)}^2+{(H_i-{\tilde{H}}_0)}^2}---\left(1\right)$

线性化后可得

$v_{L_i}=\frac{(N_0-N_i)}{L_i^0}{v}_N+\frac{(E_0-E_i)}{L_i^0}{v}_E+\frac{(H_0-H_i)}{L_i^0}{v}_H-(L_i-L_i^0)---\left(2\right)$

其中为长度改正数；(v_N,v_E,v_H)为待测点坐标改正数；(N₀,E₀,H₀)为待测点近似坐标；L_i为对中杆长度与天线高之和；为观测点坐标与待测点近似坐标计算的距离。

8.如权利要求4所述的一种RTK测量接收机倾斜测量方法，其特征在于，所述根据观测点与待测点间高差与距离的倾角余弦关系，建立第二类观测方程的步骤中，设观测点观测倾角为则观测方程为

${\tilde{L}}_i^{\′}={\tilde{H}}_0-H_i---\left(3\right)$

其中整理后可得

v_L'＝v_H-(H_i-L_icosθ_i-H₀)(4)

联立式(2)、(4)，可建立误差方程组

V＝BX-W,P(5)

式中，P为权矩阵，根据观测值的精度设定，若没有此信息，可设置为单位阵。

9.如权利要求8所述的一种RTK测量接收机倾斜测量方法，其特征在于，包括以下步骤：根据加权最小二乘准则，可得

$\hat{X}={\left(B^{T}PB\right)}^{-1}{B}^{T}PW---\left(6\right)$

最后通过循环迭代计算出地面待测点的平面位置与大地高令

$A=\left[\begin{array}{ccc}\frac{N_0-N_1}{L_1^0}& \frac{E_0-E_1}{L_1^0}& \frac{H_0-H_1}{L_1^0}\\ \frac{N_0-N_2}{L_2^0}& \frac{E_0-E_2}{L_2^0}& \frac{H_0-H_2}{L_2^0}\\ \mathrm{...}& & \\ \frac{N_0-N_n}{L_n^0}& \frac{E_0-E_n}{L_n^0}& \frac{H_0-H_n}{L_n^0}\end{array}\right]$

A为系数矩阵，计算ATA得到维度为3行3列的单位矩阵，将单位阵的对角线元素相加后取平方根得到质量标的。