测绘程序设计课件来自测绘魏志刚.ppt

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