一种基于三角高程网的大气折光系数反演方法

摘要

本发明公开了一种基于三角高程网的大气折光系数反演方法，包括以下步骤:三角高程网数据采集；利用三角高程测量方法计算网中各测边往测高差h_1i与返测高差h_2i；计算网中独立三角形个数n；计算各个三角形的三边之和C_j；计算各个三角形在网中的权P_j；计算各个三角形的高差闭合差ω_j；施加约束条件y；根据约束条件y与大气折光系数k的关系，绘制曲线图；利用曲线图与约束条件，取k值；计算往返测高差不符值Δ_i，在满足Δ_i与ω_j不超限的前提下，找出使y值达到最小的k值。本发明对于大气折光系数的选择更具动态性，而不是单纯使用大气折光参考系数；不需求解每条边测量时的k值，而是根据相关条件，求解出一个最适合全网的k值，该值与实际吻合度较高，在一定程度上提高了该方法的适用性。

主权项

一种基于三角高程网的大气折光系数反演方法，包括以下步骤：1)、通过全圆观测法观测三角高程网的数据，所述三角高程网包括m个测点，其中，m≥4，在每个测点上观测其他测点，其中，观测值包括天顶距和斜距；其特征在于：还包括以下步骤：2)、利用三角高程测量方法计算所述三角高程网中测段的往测高差h₁₁，h_1i，…h_1n和返测高差h₂₁，h_2i，…h_2n，n为所述三角高程网中测段的数量，h_1i和h_2i均为与大气折光系数k相关的函数；3)、计算所述三角高程网能形成的三角形个数2t；4)、计算三角形的高差闭合差在所述三角高程网中的权重P₁，P_j，…P_2t，所述高差闭合差包括往测高差闭合差和返测高差闭合差，其中，P_j＝1/C_j，其中，C_j为第j个三角形的周长，j＝1,2,…,2t；5)、利用步骤4)得到的权重P_j计算所述三角形的高差闭合差ω₁，ω_j，…ω_2t，其中，j＝1,2,…,2t，ω_j为与大气折光系数k相关的函数；6)、令约束条件$y=\frac{\sqrt{\underset{1}{\overset{2t}{\Sigma }}{\omega }_j^2/C_j^2}}{\underset{1}{\overset{2t}{\Sigma }}{P}_j};$7)、根据步骤6)绘制y与大气折光系数k的曲线图；8)、根据步骤7)得到的曲线图，得到y的最小值，根据步骤6)计算大气折光系数k；9)、利用步骤8)得到的k，计算得到所述三角高程网中测段的往测高差h₁₁，h_1i，…h_1n和返测高差h₂₁，h_2i，…h_2n以及所述三角形的高差闭合差ω₁，ω_j，…ω_2t；计算所有测段的往返差不符值Δ_i＝h_1i‑h_2i，其中，i＝1,2,…n；如果并且Δ_i≤100S_i，其中，ρ为206265″，S_o为第j个三角形的三条边的长度，o＝1,2,3，S_i为测段长度，则当前k值为所述三角高程网施测时的大气折光系数，否则，重复步骤8)，并验证新的大气折光系数k。