1、 南阳师范学院 20XX 届毕业生毕业论文（设计）题 目：测量控制网平差程序设计与实现 完 成 人： 班 级： 学 制： 专 业： 测绘工程 指导教师： 完成日期： 目 录摘要 (1)0 引言 (1)1 整体框架设计 (2)2 特殊情况控制网近似坐标推算方法 (4)2.1 测边网近似坐标推算 (5)2.2 已知点间没有直接观测量的近似坐标推算 (6)3 不同类型观测值的定权方法讨论 (7)3.1 按常规方法定权 (7)3.2 根据三角形闭合差定权 (8)3.3 用平差得到的验后中误差来定权 (8)4 程序应用算例 (9)4.1 程序验证 (10)42 不同定权方法平差结果精度分析 (11)5
2、结束语 (13)参考文献 (14)Abstract (15)第 1 页 共 16 页测量控制网平差程序设计与实现摘要：根据控制网间接平差的基本原理，基于 Visual C+6.0 设计和实现了一款控制网平差计算程序 GNetAdj，该程序可实现水准网、测角网、测边网、边角网等常规控制网的平差计算与精度评定重点探讨了特殊情况控制网的近似坐标推算方法，以及不同类型观测值较为合理的定权方法最后通过一算例验证了程序的正确性关键词：控制网；方差估计；平差；定权方法；近似坐标推算0 引言控制网平差计算相当复杂，计算量很大，对于规模较大的控制网必需借助计算机及计算软件才能计算出来虽然现在大多数的控制网平差软
3、件功能很齐全，使用也比较方便，但有的地方还不够理想，如缺乏对角度观测量的平差，对于已知点之间没有直接观测量的网不能推算坐标；对于边角同测的网，两类观测值仅按常规方法来定权，往往导致先验和验后中误差相差较大，对平差结果产生一定影响鉴于以上问题，本文将重点讨论特殊情况下控制网的近似坐标推算方法，如单纯测边网、已知点之间没有直接观测值的情况，另外还对不同类观测值的定权方法进行探讨1 整体框架设计图 l 是 GNetAdj 的运行界面，程序主要包含的类及其功能如图2 所示，其中 WangView 类显示网的图脉形需要将实际网的坐标转换第 2 页 共 16 页到程序视图的坐标，在图上通过点击鼠标还能显示
4、待定点的误差椭圆及其参数图 1 程序运行界面图 2 程序涉及的几个重要的类GNetAdj 是基于 Visual C+6.0 设计的多文档应用程序，其菜单栏包含程序的所有功能，除了常规控制网平差计算外，还有坐标第 3 页 共 16 页转换、交会法计算坐标等功能，工具栏的各个按钮对应各项功能由于控制网平差涉及大量的计算，而且不同类型的网的属性和计算大致相同，因此构造一个间接平差类，把控制网的大部分属性和计算封装起来对于控制网平差，通过新建一个文本，编辑已知点坐标和观测数据，或打开已编辑好的数据文本，在菜单栏或工具栏上选择所要平差的类型，在弹出的对话框(见图 3)中设置好相关参数，点击“确定”按钮后
5、程序就进行平差计算并显示结果对话框设置的参数有已知点点数、待定点点数、各类误差及观测值个数等，根据这些参数，程序先读取数据文本中的已知点，接着读取方向或角度观测值，最后读取边长观测值，并存储到相应的变量中程序读完数据后便开始计算待定点的近似坐标，进而计算权阵、误差方程的系数阵及常数项，计算过程涉及到的矩阵均采用一维数组来存储，通过对矩阵进行转置、相乘、求逆等运算，最后求得坐标改正数为提高精度，平差计算采用迭代方式，直到改正数收敛为止最后进行精度评定并输出平差结果如果读取数据或计算过程出错，程序会中止平差计算，并返回错误信息，好让用户检查哪第 4 页 共 16 页里有问题结果有文本和图像，前者是
6、待定点的坐标平差值、观测值的平差值、各类中误差等数据，后者是控制网的图形，利用消息映射机制，当光标在待定点附近时，通过点击鼠标左键会显示该点的误差椭圆及其参数，点击右键则取消显示另外，程序的其它功能通过对话框的形式来操作，在相应的编辑框中输人数据，点击相关按钮，计算结果就显示出来，非常直观明了第 5 页 共 16 页2 特殊情况控制网近似坐标推算方法2.1 测边网近似坐标推算虽然 GNetAdj 采用多种方法(极坐标法、前方交会法、测边交会法)来推算待定点的近似坐标，但如果只根据边长观测数据是不能唯一确定待定点的坐标，因为单靠边长程序判读不了点的位置关系，而能否平差测边网关键就在于此若程序判断
7、不了，则由用户来判断，如遇到不能确定点的位置时，程序会弹出图 4 的对话框让用户判断用户根据实地点的分布或草图就可作出正确判断，这样就能推算坐标并进行平差为尽量减小用户的判断，对有 3 个或以上已知点的测边网(如图5 所示)，GNetAdj 先用以下方法来推算坐标：取其中两个已知点M1、M2，三角形 M1M2 P 中 3 条边长均知道，可用余弦定理求出三角第 6 页 共 16 页形中各个内角，再用极坐标法或前方交会法分别按顺时针和逆时针各计算出一个 P 点可能的坐标，再由这两个坐标分别与第 3 个已知点 M3 求距离，哪一个距离与 P 到 M3 的距离观测值相接近，就取相应的坐标作为 P 的近
8、似坐标心 J程序会先检验 3 个已知点是否在同一条直线上，不在才采用此法若经此方法推算后还有待定点算不出坐标，将会弹出图 4 的对话框要求用户作判断22 已知点间没有直接观测量的近似坐标推算有些控制网已知点比较分散，已知点之间没有进行观测，如图6 所示的网，只有 A、B 是已知点，这样就不能直接计算待定点的近似坐标GNetAdj 采用平面相似变换法(四参数法)来解决这个问题若是测角网，搜索某一观测值得到两个端点，对这两个点模拟一套坐标并设为已知点，其它点均为待定点，计算所有点在模拟坐标系中的坐标，然后根据原来已知点的坐标及其模拟坐标来计算两第 7 页 共 16 页个坐标系的转换参数(两个平移参
9、数、一个旋转参数和一个缩放参数)，利用转换参数将模拟坐标转为实际网的近似坐标若是边角网，由于推算坐标时可能用到边长观测值，因此程序选择某一观测边的一个端点并赋予模拟坐标，然后假设该边的方位角为 ，用坐标正算可求得另一个端点的坐标，以这两个点作为已03知点来推算其它待定点的模拟坐标，然后计算转换参数，最后将模拟坐标转为实际网的近似坐标3 不同类型观测值的定权方法讨论角度和边长的权是一种比例关系，为定权故必须精确知道观测值的中误差，但在平差前这个值是不知道的，只能设法找出它们尽量好的估值来GNetAdj 使用 3 种方法来确定观测值的先验中误差，并用其来对观测值定权第 8 页 共 16 页3.1
10、按常规方法定权测角中误差 根据测角仪器的类型和观测测回数，参照相应m等级的三角测量精度来决定，测边中误差 按测距仪器出厂时的标ism称精度公式 =amm+b D km 来确定，然后按is610（31）2mP（32）20iis来定权3.2 根据三角形闭合差定权对于测角网，GNetAdj 先找出网中三角形的个数 k，并计算各三角形的闭合差，若闭合差大于 3 倍的测角中误差，则终止平差并显示超限信息，供用户检查如不超限则继续往下计算，若 k20 则按三角形闭合差并依菲列罗公式3wmn计算测角 中误差。3.3 用平差得到的验后中误差来定权当控制网采用边角同测的观测方案时，将该网分成角网和边网，然后各自
11、按角网和边网单独平差，由改正数分别计算测角中误差及测边中误差：1vmr（33）第 9 页 共 16 页2ssvmr（34）（式中 1r、 2为测角网和测边网的多余观测数）2200,i ii is ss mmPp:再 按 计 算 侧 边 中 误 差 ， 然 后 按 P定 权 。若能构成独立完整的角网，但不能构成边网，则按下面的方法：首先对角网进行平差，求其中误差 1，令其自由度(多余观测数)为1f；按 =P，2iSsm定权，对含有边长观测量的角网进行整体平差，求得平差后观测中误差的估值 2，令其自由度为 2f；2211fm，求边长观测值的中误差估值。可以证明中误差估值21与是互相独立的；利用先后
12、求得的中误差估值 12m与 重新定权，即有21,sP，然后对整体网平差，求出新的中误差估值 21，如果 21m满足 21210.m， ，则计算可终止否则，取用 代替重新定权平差，迭代计算直到满足 21210.m4为 止。4 程序应用算例图 7 是用 Auto CAD 2006 模拟出来的控制网，通过标注可以得到边长、方向、坐标的真值，网中 M1、M2、M3 为已知点，其它 6 个为待定点，共获得 34 个方向观测值和 16 个边长观测值，将其作为平差计算的基础数据第 10 页 共 16 页4.1 程序验证为验证程序的正确性，选择“平差易”(南方测绘仪器公司开发的控制测量数据处理软件)与 GNe
13、tAdj 作对比，用不加任何误差的34 个方向观测值和 16 个边长观测值进行平差计算，结果见表 1表 1 观测值无误差情况下两个程序的平差结果与真值对比 m真值 平差易结果 GNetAdj 结果点名X Y X Y X YP101 30654.5073 41268.2041 30654.5072 41268.2042 30654.5072 41268.2043P102 30886.0734 41367.9236 30886.0735 41367.9236 30886.0735 41367.9237P103 30702.8724 41590.8252 30702.8723 41590.8253
14、30702.8723 41590.8254P104 30557.7772 41884.5226 30557.7772 41844.5226 30557.7772 41844.5227P105 30657.4386 42171.5431 30657.4388 42171.5431 30657.4389 42171.5432P106 30953.4914 41618.6879 30953.4913 41618.6880 30953.4913 41618.6881第 11 页 共 16 页从表 1 可看出，平差易和 GNetAdj 计算的结果与真值相差很小，均不超过 03 mm，平差易和 GNetA
15、dj 计算的结果相比不超过01 mm，从统计意义上来讲，两者计算的结果是一致的，这就证明了 GNetAdj 程序平差计算是正确的4.2 不同定权方法平差结果精度分析为了分析不同类型观测值的定权方法对平差结果的影响，利用Matlab 中的 random 对方向和边长观测值分别模拟5函 数“632m+10kD和随机误差，然后分别用南方平差易和GNetAdj 进行平差。图 8 是根据平差结果利用公式22sxy:（式中 xy为坐标平差值， xy:为坐标真值）计算坐标平差值相对于真值的偏离量从以上可看出，在先验中误差能够准确得到的情况下，先验和验后中误差趋于一致，平差结果最接近真值若先验中误差偏离实际值
16、越大，平差结果偏离真值也越大而 GNetAdj 用前述 3 中 3)的方法来定权，平差结果和验后中误差并不随着先验中误差的变化而变化图 8 比较 4 套结果与真值的偏离量，显然，利用平差得到的验后中误差来定权得到的结果更接近真值，证明该方法是合理可行的第 12 页 共 16 页5 结束语GNetAdj 可以对水准网、测角网、测边网、边角网等常规控制网进行平差计算与精度评定，在此基础上加入了 3 个特色功能：通过用户的判断信息能够对纯粹的测边网进行平差；在已知点较少且分散或障碍物较多的地区布网，即使已知点间没有直接观测量，利用 GNetAdj 仍能进行平差，这就拓宽了网型的选择和布网的形式，充分
17、利用已有控制点从而减少布网时间和建网费用；通过对不同类型的观测量采取相应的定权方法使得先验中误差和验后中误差趋于一致，平差结果更为精确通过与平差易的比较，验证了 GNetAdj的正确性，同时也表明 GNetAdj 具有很好的实用价值另外，程序还集成几个常用的计算工具，如坐标转换的功能，这样可以方便地第 13 页 共 16 页将平差出来的结果转换到其它坐标系中参考文献1成虎工程项目管理M北京：中国建筑工业出版社2002：58-72.2黄渝祥后评价的内涵及其在项目管理中的应用J石油化工技术经济，1997，12(1)：55-573中国水利经济研究会水利建设项目后评价理论与方法M北京：中国水利水电出版
18、社，2003：33-46.4张道军逻辑框架法在水利工程项目后评价中应用实例J水利经济，1998，6(3)：34-455叶义成，柯丽华，黄德育系统综合评价技术及其应用M北京：冶金工业出版社 2006：136-148.6丁桑岚环境评价概论M北京：化学工业出版社，2001：96-134.第 14 页 共 16 页Survey control net adjustment program design and implementationAbstract：GNetAdj is based on the principles of indirect adjustment，and it is design
19、ed by Visual C+6 0It can realize the adjustment and estimate of precision for level network，triangulation network，trilateration network，triangulateration network and other conventional control networksSome methods of calculating approximation coord；nates for special control network are discussed，and some methods of fixing the weight of different observations are analyzedFinally，the program is verified by a case，and the result shows that it has good practical valueKey words：control network；estimation of variance；adjustment ；methods of fixing weight；calculation of approximate coordinates
