1、测绘程序设计,控制网程序设计,共8学时 上机实验4学时,4-3 平面控制网平差程序设计控制网的平差方法基本上分为两类：条件平差法和间接平差法。 在手工计算时，主要考虑的是选择哪种法方程使得阶数较少，又便于解算的方法，因此平差方法的选择依控制网的具体情况而定。 然而对于电算来说，主要考虑的是程序编制的难易度和通用性。于是用条件平差方法编制通用程序就必须顾及所有可能的条件式，而要做到这一点，即使只要求程序适合于测角网，程序编起来也是相当困难。但是对于间接平差来说，情况就完全不同了。因为无论任何一个平面控制网，通常只包含两类观测值，即：方向观测值和边长观测值，,间接平差的基本公式设：观测值为L，其权
2、为P，相应的改正数为V，必要观测值的个数为t，选定未知数X，则有误差方程： V=BX+L 法方程： BTPBX+BTPL=0 法方程的解： X=-(BTPB)-1 BTPL 单位权中误差：,一、误差方程式的列立间接平差以未知数(待定点的坐标)为平差对象，一个观测值对应着一个误差方程式，误差方程式包括两类，一类为方向观测值误差方程式，另一类为边长观测值误差方程式。,ajk、bjk称j、k方向的方向系数,1方向观测值误差方程式1）方向观测值误差方程式的一般形式Vjk-dzj+ajkdxj+bjkdyj-ajkdxk-bjkdyk+ljk 式中，dZi为测站i的定向角改正数，dX，dY为坐标改正数，
3、其权为Pij。,Zj定向角,对于某i测站其误差方程式为,其和方程式为,(当Mij 时),(当Mij 时),2）近似定向角的计算定向角的近似值可任选，一般采用下列方法计算。(1)以零方向的方位角 为定向角近似值 ，式中j为本测站的零方向。(2)以一个测站上各方向的近似定向角 的平均值作为测站的近似定向角,Mij为从i到j的观测值,3）方向系数aij、bij的计算,当改正数以厘米为单位， 以米为单位时,4）常数项Lij的计算,5) 方向误差方程式的列立：以点的编号为序，从各测站的零方向开始，依次列出各方向观测值的误差方程式，并且列出一个误差方程式后直接组成法方程，即误差方程式不予存贮。用mi表示第
4、i测站在DD IODUO数组中方向观测值的起始地址，用ni表示终了地址，用dv表示误差方程式的常数项，用p表示观测值的权。程序如下：,6)法方程式的组成,2边长观测值误差方程式1）边长观测值误差方程式的一般形式,或,其中 ，Sij0 由近似坐标计算而得，Sij为实测边长。,2）边长观测值误差方程式的列立以观测边长在DLIOS数组中的存放顺序，逐个列出边长观测值的误差方程式，并组成法方程式。用dv表示误差方程式的常数项，P表示观测值的权。程序如下：,二、方位角条件与边长条件的列立1方位角条件,2边长条件,3条件式的列立平面控制网按坐标间接平差中，当存在条件式时未知数之间是不独立的。解决的方法一般
5、有以下两种。 (1)利用未知数之间的条件式，消去相关未知数，按一般的间接平差法解算。 (2)列出未知数之间的条件式，按附有条件的间接平差法解算。 由于消去未知数的方法不便于程序化，因此一般采用附有条件的间接平差来解算。其中a、b系数的计算与误差方程式中的计算完全一样，而且条件式不需要法化，直接附在法方程后即可。,三、法方程式的组成与解算1法方程式的组成2法方程式的解算高斯消去法,下图中A为已知水准点，其高程HA237.483m，为了确定B、C、D点高程，共观测了5个高差，高差观测值及相应水准路线的路线长度为,试按间接平差法求B、C、D点高程的平差值。,实例：,4-4 高程控制网平差程序设计这里
6、主要介绍以测段的高差为平差元素、待定点的高程为未知数的间接平差法的解算过程。 一、高程控制网间接平差的步骤及总框图1高程控制网间接平差计算步骤(1)计算待定点的近似高程X(x1，x2 )T(2)列出误差方程(3)组成法方程(4)解算法方程，求得dX(dx1，dx2)T(5)求得平差后的高程XX十dX(6)如果dX中的绝对值最大的一个小于给定的值，则平差结束，否则由新的近似值再一次平差， 即重复第(2)、(3)、(4)、(5)步。(7)精度评定。,2高程控制网间接平差流程图高程控制网间接平差是以测段高差为平差元素。一般为三角高程网和水准高程网的平差计算。要能够评定网中任意点高程和任意两点间高差的
7、精度。主要包括以下三个主模块。(1)数据输入模块(2)平差计算模块(3)精度评定及成果输出模块,3数据结构(1)网形的编号为了编程方便，首先约定各高程点编号由小到大按自然数的顺序编 码，不可重复也不可缺少。已知点优先于待定点，靠近已知点的待定点编号要尽量小。,(2)简变约定(WORD为无符号整型)WORD IKP：已知点点数 WORD IUP：未知点点数 WORD IQ：总点数。IQ：IKP+IUP WORD IGL：已知概略高程点数 WORD IOH：高差观测个数 WORD IZ：未知数总个数。IZ2*IUP WORD UD：评定精度的个数 Double EM：验前单位权中误差(mm),(3
8、)数组(为了叙述方便，忽略掉各一维数组的零单元，因此说明数组时下标加1) Double DH10H+1：高差观测值数组 Double PDH10H+1)：定权元素数组 double HIQ：高程数组 WORD ZDHAI0H+1，ZDHBI0H+1：高差的两个端点点号数组。高差观测值是有方向的，ZDHAI0H+1为起点，ZDHB IOH+1为终点。 DWORDZJD2*IQ+1评定精度元素数组，以XXXXXXXX方式存放，前四位为一个点号，后四位为另一个点号。评定高程精度时只有一个点号。 Char NameIQ+118：以点号为序存放各点的点名。不得大于18个字符。 Double RMIZ*(
9、IZ+1)2+1，RWIQ+1：法方程系数和法方程常数项数组。,试按间接平差法求B、C、D点高程的平差值。,实例：,下图中A为已知水准点，其高程HA237.483m，为了确定B、C、D点高程，共观测了5个高差，高差观测值及相应水准路线的路线长度为,例42 如图所示的水准高程控制网，已知点2个， 未知点3个，总点数5，已知概略高程点数0，高差观测值个数7，评定精度的个数3。网形编号如图42所示。,二、待定点高程计算1原理待定点的高程计算是为了得到未知数的高程近似值。在端点号数组ZDHAIOH+1和ZDHBI0H+1中，如果该两点有一点是已知或已算出的，而另一点未知，则可由其中已知高程的点与未知点
10、间的观测高差运算得到未知点的高程。为了区别哪些点已知或已算出，哪些点待计算，需要设立标志数组NameBzIQ，当NameBzI1时，说明第1点高程已知或已算出，NameBzI=0时，说明第1点待计算。为了防止由于编号不当或观测值排列顺序的限制而一次不能完全推算出所有的待定点，采用多次推算的方法来解决。2程序,三、高程控制网的误差方程和法方程1误差方程式的列立(1)设有两点i和j，则误差方程式为,(2)权的确定 以测站数ni定权， Pi=1ni 以中误差定权， Pi=lmi2； 以距离定权： Pi=1Si 以距离的平方定权： Pi=1Si2 (3)未知数的排列 以点的编号为序，从第一个待定点开始
11、，依次排为第一未知数、第二未知数。即第一未知数为IKP+1点，第二未知数为IKP+2点等依次类推。,2法方程的组成 (1)法方程的存放 只存放其包括主对角线的下三角元素，并存放在一维数组RM中，RM1至RM IZ*(1Z+1)2的内容依次为paa、pba、pbb、pca、pcb、pcc 。 (2)法方程式的组成 高差观测值的误差方程式的一般形式为,列出一个误差方程式后立即法化，并且要对i和j是否为已知点作处理。每一个误差方程最多只能 有三个法方程式的系数和二个常数项，其组成分别为：paa、pab、(或pba)、pbb、pal、pbl，将其分别 累加于RM敷组和RW敷组的相应单元，当所有的误差方
12、程式都法化后组成最终的法方程系数 和其对应的常数项误差方程式的常数项以毫米为单位。 3程序 字符串ADJC用于控制定权方式，如果ADJC：“水准网”，则用定权元素的倒数定权；否则用定权元素倒数的平方定权。,4-5 坐标换带与新旧坐标转换程序设计一、高斯投影坐标换带高斯投影坐标换带计算常用的两种方法是，查表和使用高斯投影正、反算公式间接进行换带计算。后一种方法其实质是把椭球面上的大地坐标作为过渡坐标。首先把I投影带内的平面坐标(x，y)1，利用高斯投影反算公式换算成椭球面上的大地坐标(B，L)，进而得到 （l为经差)；然后再由此大地坐标(B，L)，利用投影正算公式换算成相邻带的平面坐标(x，y)
13、，下面是高斯投影坐标换带计算常用的数学模型。1.高斯投影正算公式,2.高斯投影反算公式,利用高斯投影正、反算公式间接进行换带计算，精度可靠，通用性强，它不仅适用于66带，33带以及63带互相之间的邻带坐标换带，而且也适用于任意带之间的坐标换带。,高斯投影坐标变换计算器,输入区 输出区纬度B X=经度L Y=,坐标正算,0,0,0,0,输入区 输出区X= 纬度B=Y= 经度L=,坐标反算,0,0,0,0,二、新旧坐标换算旧坐标系指先前建立和计算了的控制点坐标系，而新坐标系指后来建立或重新计算的控制点坐标系，这就必然存在着坐标系的换算问题。换算方法不外乎有两种严密方法和近似方法。所谓严密方法，就是
14、将旧网的全部观测资料，重新归算到新坐标系中，重新平差并算出各点的新坐标值。而近似方法是在旧网原始观测资料不足或其它工程急需的情况下，常采用的一种方法。采用近似方法对新旧坐标换算，必须有足够的新旧网重合点，根据重合点的差值，按一定的规律修正旧网的各点坐标值，使旧网与新网达到最佳吻合。下面是常见几种方法的数学模型。,1简单变换方法(赫尔默特法)这种变换法的实质是使旧网坐标系平移、旋转和进行尺度因子改正。其敷学模型如下,2正形变换法,式中，p0，pl，p2，p3，q0，q1，q2，q3，分别为坐标交换参数的最或是值。,实验三：控制测量常用程序设计 一、实验目的与要求初步掌握控制测量数据处理过程中的数据结构类型；初步学会控制测量常用软件设计方法、过程、程序代码编写方法及程序调试过程。 二、实验安排本实验需要4学时。实验每个小组1人，用C语言在计算机上完成程序代码的编写，并调试通过，然后用实际算例印证程序的正确性。2. 时间安排：周二3-4节。周三7-8节。 三、实验步骤及要点 课余时间完成预先布置的控制测量数据处理问题，完成程序代码的纸上编写和设计工作； 在计算机上写入程序代码，进行调试； 用预先准备好的算例进行印证，要求每人至少准备2组典型算例，并打印测试结果，附在实验报告中。 四、要求1.代码编写要独立完成，算法可以小组为单位进行研讨。2.每次完成实验后应该撰写实验报告。,
