一种测边前方交会测点的方法

摘要

本发明旨在提出一种新的测边前方交会测点的方法，通过最小二乘计算未知点的坐标，其有益效果是：1）当已知点相距较近时，其连线方位角的误差不会传递给未知点；2）对边长观测值进行合理定权，并采用最小二乘进行解算，必要时进行迭代，便于编程实现，能获得未知点坐标的最优解；3）能同时获得未知点的高程。

说明书

技术领域

本发明涉及一种测边前方交会测点的方法，属于测绘地理信息领域，可以用 于测定未知点的坐标。

背景技术

交会测点是测绘领域常用的一种测点方式，根据交会方式可分为前方交会、 后方交会；根据观测值类型可以分为测角交会、测边交会和边角交会。其原理是 通过观测未知点同已知点之间的几何关系，推算未知点的空间坐标。

已有的测边前方交会原理如附图1所示，图中，A、B为已知点，P为未知点，为 求P点坐标，测量了边长D_AP和D_BP，传统的计算公式为：

$\left(\begin{array}{c}{x}_P=x_A+{\mathrm{\Delta x}}_{AP}=x_A+D_{AP}{\mathrm{cos\alpha }}_{AP}\\ y_P=y_A+{\mathrm{\Delta y}}_{AP}=y_A+D_{AP}{\mathrm{sin\alpha }}_{AP}\end{array}\right)$

其中${\alpha }_{AP}={\alpha }_{AB}-\angle BAP={\alpha }_{AB}-arccos\left(\frac{D_{AP}^2+D_{AB}^2-D_{BP}^2}{2D_{AB}{D}_{AP}}\right),$α_AB为根据A、B已知坐 标反算的方位角。此种方法的缺点在于1)由于已知点的坐标不可能绝对精确， 若A、B相距较近，则依据其已知坐标计算的方位角α_AB之误差将会较大，导致P 点的坐标精度较低；2)当有多于2个已知点可以利用时，传统方法要将相邻的 已知点两两组合，分别按照上述公式进行计算，计算步骤较为复杂，不利于编程 实现，最后对各组结果采取简单平均的方法，模型不够严密；3)传统方法只考 虑平面坐标的获取，而没有考虑高程问题。

发明内容

本发明旨在提出一种新的测边前方交会的计算方法，以便当已知点距离较近 时，测边交会仍能获得较高精度；旨在能同时求解出未知点的平面坐标和高程； 旨在能更有效地利用多个已知点，并通过最小二乘原理进一步求得未知点的最优 解。

本发明采用如下技术方案：

一种测边前方交会测点的方法，包括以下步骤：

步骤1，在待定点P附近，选择n个(n≥2)坐标已知点P1,P2,……,Pn, 测量各点到点P的平距，记点Pi(i＝1,2,……,n)到点P的平距观测值为D_i；

步骤2，求取点P的近似坐标，

$\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)=\left(\begin{array}{c}{x}_j\\ y_j\end{array}\right)+\left(\begin{array}{c}{D}_j{\mathrm{cos\alpha }}_{j0}\\ D_j{\mathrm{sin\alpha }}_{j0}\end{array}\right)$

其中，${\alpha }_{j0}={\alpha }_{jm}-\arccos \left(\frac{D_j^2+D_{jm}^2-D_m^2}{2D_{jm}{D}_j}\right),$j、m代表n个已知点中挑选的任 意两个，$\left(\begin{array}{c}{x}_j\\ y_j\end{array}\right)$表示点j的已知坐标，D_jm表示已知点j、m间的平距，α_jm表示已 知点j、m连线的已知方位角；

步骤3，列立如下误差方程，

$\underset{n\times 1}{V}=\underset{n\times 2}{B}\underset{2\times 1}{\omega }-\underset{n\times 1}{l}$

其中，$V=\left(\begin{array}{c}{v}_{D_1}\\ .\\ .\\ .\\ v_{D_n}\end{array}\right)$表示各距离观测值的改正数，$B=\left(\begin{array}{cc}\frac{(x_0-x_1)}{D_1^0}& \frac{(y_0-y_1)}{D_1^0}\\ .& .\\ .& .\\ .& .\\ \frac{(x_0-x_n)}{D_n^0}& \frac{(y_0-y_n)}{D_n^0}\end{array}\right),$$\omega =\left(\begin{array}{c}{\Delta }_x\\ {\Delta }_y\end{array}\right)$表示点P的坐标改正数，$l=\left(\begin{array}{c}{D}_1-D_1^0\\ .\\ .\\ .\\ D_n-D_n^0\end{array}\right),D_i^0=\sqrt{{(x_i-x_0)}^2+{(y_i-y_0)}^2};$

步骤4，计算点P的坐标改正数，

ω＝(B^TWB)^-1B^TWl

其中，$W=diag(\frac{1}{{\delta }_1^2},\frac{1}{{\delta }_2^2},...,\frac{1}{{\delta }_n^2}),$δ_i表示D_i的测距中误差；

步骤5，计算点P的坐标平差值，

$\left(\begin{array}{c}{\hat{x}}_0\\ {\hat{y}}_0\end{array}\right)=\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)+\omega =\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)+\left(\begin{array}{c}{\Delta }_x\\ {\Delta }_y\end{array}\right)$

上述步骤2中所述P点的近似坐标可取$\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)=\frac{1}{n}\left(\begin{array}{c}\underset{i=1}{\overset{n}{\Sigma }}{x}_i\\ \underset{i=1}{\overset{n}{\Sigma }}{y}_i\end{array}\right),$以简化计算，后续 可通过迭代保证精度不会损失。

步骤3中权矩阵W可取为W＝diag(1,1,…,1)，以简化计算；或 $W=diag(\frac{1}{{\delta }_1^2+{\delta }_{P1}^2},\frac{1}{{\delta }_2^2+{\delta }_{P2}^2},...,\frac{1}{{\delta }_n^2+{\delta }_{Pn}^2}),$其中δ_Pi表示点Pi的平面点位中误差， 后者能更合理地给各观测值定权，进一步提高平差结果的精度。

此外，还可采用迭代方式进行m(m≥1)次平差，其中1次平差是指步骤3 至5的一次顺序执行，即将第i次平差所得结果视为第i+1次平差的输 入x₀、y₀，反复计算，直至Δ_x、Δ_y均小于某一阈值，使平差结果更趋近真值。

此外，还可在步骤1中增加以下操作：量测点Pi与点P间的高差，记为h_i；

步骤5中增加如下过程：

1)计算点P的k(1≤k≤n)个初始高程值H_P(i)，

H_P(i)＝H_i+h_i

其中H_i表示点Pi的已知高程；

2)计算点P的高程平差值，

${\hat{H}}_0=\frac{\underset{i=1}{\overset{k}{\Sigma }}{H}_P\left(i\right)}{k}$

因此，除了平面坐标，高程也可获得。

此外，还可以采用下式计算点P的高程平差值，

其中δ_Hi表示已知点Pi的高程中误差，δ_hi表示点Pi与点P间高差观测值的 中误差，由此通过加权平均的方式，可获得比简单取平均更为精确的高程平差值。

与现有技术相比，本发明的有益效果是：1)当已知点相距较近时，其连线 方位角的误差不会传递给未知点；2)对边长观测值进行合理定权，并采用最小 二乘进行解算，必要时进行迭代，便于编程实现，能获得未知点坐标的最优解； 3)能同时获得未知点的高程。

附图说明

图1为两个已知点的测边前方交会原理示意图；

图2为n个已知点的测边前方交会原理示意图。

具体实施方式

下面结合附图和实施例对本发明作进一步说明。

如附图2所示，P1、P2、P3、P4均为已知点，将点Pi的北坐标、东坐标、 高程分别记为x_i、y_i、H_i；点Pi的平面点位中误差记为δ_Pi，高程中误差记为δ_Hi。 点P的坐标按以下步骤进行求算：

步骤1，量测各点与点P的平距和高差，其中平距可采用钢尺、激光测距仪、 全站仪等仪器测定；高差可采用水准仪、全站仪、卷尺等仪器测定，点Pi点P 的平距记为D_i，测距中误差记为δ_i；点Pi与点P的高差观测值记为h_i，高差中 误差记为δ_hi；

步骤2，在已知点中，任意挑选出两个点，如P1和P3，采用如下公式计算 点P的坐标初始值：

$\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)=\left(\begin{array}{c}{x}_1\\ y_1\end{array}\right)+\left(\begin{array}{c}{D}_1{\mathrm{cos\alpha }}_{10}\\ D_1{\mathrm{sin\alpha }}_{10}\end{array}\right)$

其中，D₁₃表示已知点P₁和P₃间的平距，α₁₃表 示P₁和P₃连线的已知方位角；

步骤3，列立误差方程式：

$\underset{4\times 1}{V}=\underset{4\times 2}{B}\underset{2\times 1}{\omega }-\underset{4\times 1}{l}$

其中，$V=\left(\begin{array}{c}v1\\ v2\\ v3\\ v4\end{array}\right),B=\left(\begin{array}{cc}\frac{(x_0-x_1)}{D_1^0}& \frac{(y_0-y_1)}{D_1^0}\\ \frac{(x_0-x_2)}{D_2^0}& \frac{(y_0-y_2)}{D_2^0}\\ \frac{(x_0-x_3)}{D_3^0}& \frac{(y_0-y_3)}{D_3^0}\\ \frac{(x_0-x_4)}{D_4^0}& \frac{(y_0-y_4)}{D_4^0}\end{array}\right),\omega =\left(\begin{array}{c}{\Delta }_x\\ {\Delta }_y\end{array}\right),l=\left(\begin{array}{c}{D}_1-D_1^0\\ D_2-D_2^0\\ D_3-D_3^0\\ D_4-D_4^0\end{array}\right),$

$D_i^0=\sqrt{{(x_i-x_0)}^2+{(y_i-y_0)}^2};$

步骤4，计算点P的坐标改正数，

ω＝(B^TWB)^-1B^TWl

其中，$W=diag(\frac{1}{{\delta }_1^2+{\delta }_{P1}^2},\frac{1}{{\delta }_2^2+{\delta }_{P2}^2},\frac{1}{{\delta }_3^2+{\delta }_{P3}^2},\frac{1}{{\delta }_4^2+{\delta }_{P4}^2});$

步骤5，计算点P的坐标平差值：

$\left(\begin{array}{c}{\hat{x}}_0\\ {\hat{y}}_0\end{array}\right)=\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)+\omega =\left(\begin{array}{c}{x}_0\\ y_0\end{array}\right)+\left(\begin{array}{c}{\Delta }_x\\ {\Delta }_y\end{array}\right)$

并判断Δ_x、Δ_y的绝对值是否小于1mm，若否则进行迭代，重复步骤3至5，直 至满足前述条件。

步骤6，计算点P的4个高程初始值：

$\left(\begin{array}{c}{H}_P\left(1\right)=H_1+h_1\\ H_P\left(2\right)=H_2+h_2\\ H_P\left(3\right)=H_3+h_3\\ H_P\left(4\right)=H_4+h_4\end{array}\right)$

计算点P的高程平差值：

${\hat{H}}_0=\frac{\underset{i=1}{\overset{4}{\Sigma }}\frac{H_P\left(i\right)}{{\delta }_{Hi}^2+{\delta }_{hi}^2}}{\underset{i=1}{\overset{4}{\Sigma }}\frac{1}{{\delta }_{Hi}^2+{\delta }_{hi}^2}}$

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技 术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用 类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的 范围。