[发明专利]一种基于三角高程网的大气折光系数反演方法

摘要

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

主权项

1.一种基于三角高程网的大气折光系数反演方法，其特征在于：包括以下步骤：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=Σ12tωj2/Cj2Σ12tPj;]]>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。