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

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

【技术实现步骤摘要】

【技术保护点】
一种基于三角高程网的大气折光系数反演方法，其特征在于：包括以下步骤：1）、通过全圆观测法观测三角高程网的数据，所述三角高程网包括m个测点，其中，m≥4，在每个测点上观测其他测点，其中，观测值包括天顶距和斜距；2）、利用三角高程测量方法计算所述三角高程网中测段的往测高差h11，h1i，…h1n和返测高差h21，h2i，…h2n，n为所述三角高程网中测段的数量，h1i和h2i均为与大气折光系数k相关的函数；3）、计算所述三角高程网能形成的三角形个数2t；4）、计算三角形的高差闭合差在所述三角高程网中的权重P1，Pj，…P2t，所述高差闭合差包括往测高差闭合差和返测高差闭合差，其中，Pj＝1/Cj，其中，Cj为第j个三角形的周长，j＝1,2,…,2t；5）、利用步骤4）得到的权重Pj计算所述三角形的高差闭合差ω1，ω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，计算得到所述三角高程网中测段的往测高差h11，h1i，…h1n和返测高差h21，h2i，…h2n以及所述三角形的高差闭合差ω1，ωj，…ω2t；计算所有测段的往返差不符值Δi＝h1i?h2i，其中，i＝1,2,…n；如果并且Δi≤100Si，其中，ρ为206265″，So为第j个三角形的三条边的长度，o＝1,2,3，Si为测段长度，则当前k值为所述三角高程网施测时的 大气折光系数，否则，重复步骤8），并验证新的大气折光系数k。FDA00003550694900012.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：黄腾，沈月千，黄昱旻，欧乐，王成，
申请(专利权)人：河海大学，
类型：发明
国别省市：