科傻平差软件说明书

前前 言言 地面测量工程控制与施工测量内外业一体化和数据处理自动化系 统 简称科傻系统 将测量基本原理和现代科技相结合 对电子全站 仪 电子水准仪以及常规地面测量仪器进行系统的开发 以地面控制测 量 施工测量和碎部测量等测量工程为对象 实现从外业数据采集 质 量检核 预处理到内业数据处理 成果报表输出的一体化和自动化作业 流程 该系统由两个子系统组成 基于掌上型电脑的测量数据采集和处 理系统 简称 COSA HC 在掌上型电脑 RD EB2 上运行 能自动控 制和引导整个作业过程并进行质量检测 一体化程度高 操作方便 该 子系统具有水准测量 二 三维控制 碎部测量 道路测设 工程放样 等测量作业模块 具有小规模水准网 二 三维工程网的平差功能 具 有文件管理和数据通信功能 该系统灵活方便 适合外业环境 地面测量工程控制测量数据处理通用软件包 简称 CODAPS 或 COSAWIN 在微机 WINDOWS 环境下运行即可独立使用 也可与 COSA HC 联合使用 对 RD EB2 传输过来的原始观测数据进行转换 完成从 概算到平差的数据自动化处理 同时具有粗差探测与剔除 方差分量估 计 闭合差计算 贯通误差影响值估算 报表打印 网图显绘 坐标转 换与换带计算 控制网优化设计以及叠置分析等功能 本手册是为 COSAWIN 用户专门编写的 若有疏漏和不当之处 敬 请读者提出宝贵意见和批评指正 武汉测绘科技大学 武地课题组 2000 5 2 第一章第一章 概述概述 1 11 1 系统简介系统简介 科傻系统 COSA 是 地面测量工程控制与施工测量内 外业一体化和数据处理自动化系统 的简称 包括 COSAWIN 和 COSA HC 两个子系统 COSAWIN 在 IBM 兼容机上运行 COSAWIN 系统除具有概算 平差 精度评定及成果输出 等功能外 还提供了许多实用的功能 如网图显绘 粗差剔 除 方差分量估计 贯通误差影响值计算及闭合差计算等 该系统不同于其它现有控制网平差系统的最大特点是自 动化程度高 通用性强 处理速度快 解算容量大 其自动 化表现在通过和 COSA 子系统 COSA HC 相配合 可以做到由 外业数据采集 检查到内业概算 平差和成果报表输出的自 动化数据处理流程 其通用性表现在对控制网的网形 等级和 网点编号没有任何限制 可以处理任意结构的水准网和平面网 无 须给出冗余的附加信息 其解算速度快 解算容量大表现在 采用稀疏矩阵压缩存储 网点优化排序和虚拟内存等技术 在主频 166MHZ 的 586 微机上 解算 500 个点的平面和水准 控制网不到 1 分钟 在具有 20MB 剩余硬盘空间的微机上 3 可以解算多达 5000 个点的平面控制网 1 21 2 安装及运行安装及运行 1 2 11 2 1 安装安装 运行 COSAWIN 第一张系统盘上的安装程序 SETUP EXE 后 按提示操作即可 安装程序自动在程序组增加 科傻系 统 程序项 1 2 21 2 2 运行运行 在并口上插入加密狗后 点击 开始 菜单栏 程序 组项 会出现 科傻系统 程序组 单击其中程序项 科傻 系统 中的 地面控制测量数据处理系统 进入如图 1 1 所 示的系统界面 所有操作均符合 WINDOWS 操作规范 4 1 31 3 快速入门快速入门 1 3 1 系统菜单系统菜单 1 1 文件 文件 文件菜单的主要功能如 图 1 2 所示 新建新建 新建文本文件 如平面观测文件等 打开打开 打开任意文件 图 1 1 图 1 2 5 打印设置打印设置 打印机设置 单击将打开 WINDOWS 打印机设 置对话框 2 2 平差 平差 平差菜单的主要功能如图 1 3 所 示 平面网平面网 对平面网进行平差 单击 将打开 输入平面观测值文件 对话框 选择平面观测值文件进行平面网平差 高程网高程网 对水准 高程 网进行平 差 单击将打开 打开 对话框 选择 水准 高程 观测值文件进行高程平差 粗差探测粗差探测 自动探测平面控制网观测值中的粗差 若发 现粗差则自动剔除之 方差分量估计方差分量估计 对于平面网中一组或有多组不同种类或 和 精度观测值的情况 通过方差分量估计 可以使各组 观测值的精度获得最佳估计 保证平差随机模型和成果的正 确性 设置与选项 设置与选项 概算 平差 粗差探测以及坐标转换前作 相应的设置和选项 生成概算文件 生成概算文件 作概算时需要调用此项 然后在进行平 差 图 1 3 6 3 3 报表 报表 报表菜单的主要功能如图 1 4 所示 平差结果平差结果 根据平面网或高程网平差结果文件自动生成 平面或高程平差结果报表 原始观测值原始观测值 将掌上型电脑经数据通信所得到的原始观 测值文件自动生成平面高程网或高程网的原始观测值报表 4 4 查看 查看 打开或关闭工具栏和状态栏 5 5 工具 工具 工具菜单的主要功能如图 1 5 所示 其详细说明见第三 章 图 1 4 7 1 3 2 工具条 工具条 快捷按钮快捷按钮 工具条快捷按钮 如图 1 6 所示 共 有 6 类共 18 个快捷 按钮 即 文件类 新建 打开 保存文件 编辑类 剪切 复制 粘贴功能 对 活动文档有效 打印类 打印当前活动文档 包括图形 平差类 平面网平差 高程网平差 平面网网图显绘 平面网网图显绘类 图形放大 缩小 开窗放大 返回 前级 切换误差椭圆 点名和比例显示 关于类 显示系统关于信息 图 1 5 图 1 6 8 第二章第二章 平差平差 使用 COSAWIN 最常用的操作就是进行控制网平差处理 该菜单下包括平面网 高程网平差 粗差探测 方差分量估 计 设置与选项以及生成概算用文件等子菜单 本章主要介 绍平面 高程控制网观测值文件的结构及生成 网平差 设 置与选项以及生成观测值概算文件等内容 粗差探测和方差 分量估计详见第四章 2 12 1 控制网观测文件控制网观测文件 在进行平差之前 必须要准备好控制网观测文件 即平 面观测文件 取名规则为 网名 in2 和高程观测文件 取 名规则为 网名 in1 观测文件采用网点数据结构 除包 含控制网的所有已知点 未知点和观测值信息外 还隐含了 控制网的拓扑信息 可以使用系统菜单中 文件 栏下拉 新建 子菜单项 或单击工具栏左边第一个快捷键建立平面或高程观测文件 2 1 12 1 1 平面观测文件平面观测文件 平面观测文件为标准的 ASC 码文件 可以使用任何文 9 本编辑器建立编辑和修改 其结构如下 方向中误差1 测边固定误差1 比例误差1 精度号1 方向中误差2 测边固定误差2 比例误差2 精度号2 方向中误差n 测边固定误差n 比例误差n 精度号n 已知点点号 X 坐标 Y 坐标 测站点点号 照准点点号 观测值类型 观测值 观测值精度 该文件分为两部分 第一部分为控制网的已知数据 包 括先验的方向观测精度 先验测边精度和已知点坐标 见文件 的 部分 第二部分为控制网的测站观测数据 见文件的 部分 包括方向 边长 方位角观测值 为了文件的简洁和 统一 我们将已知边和已知方位角也放到测站观测数据中 它们和相应的观测边和观测方位角有相同的 观测值类型 但其精度值赋 0 即权为无穷大 第一部分的排列顺序为 第一行为方向中误差 测边固 定误差 测边比例误差 若为纯测角网 则测边固定误差和 比例误差不起作用 若为纯测边网 方向误差也不起作用 这时可输一个默认值 1 程序始终将第一行的方向中误差 值作为单位权中误差 若只有一种 或称为一组 测角 测边 10 精度 则可不输入精度号 这时 从第二行开始为已知点点 号及其坐标值 每一个已知点数据占一行 若有几种测角测 边精度 则需按精度分组 组数为测角 测边中最多的精度 种类数 每一组占一行 精度号输 1 2 参见表 2 2 如两种测角精度 三种测边精度 则应分成三组 方向中误差单位为秒 测边固定误差单位为毫米 测边 比例误差单位为 ppm 第一行的三个值都必须赋值 对于纯测 角网 测边的固定误差和比例误差可输任意两个数值 如 5 3 对于纯测边网 方向中误差赋为 1 0 已知点点号 或 点名 下同 为字符型数据 可以是数字 英文字母 大小 写均可 汉字或它们的组合 测站点 照准点亦然 X Y 坐标以米为单位 第二部分的排列顺序为 第一行为测站点点号 从第二 行开始为照准点点号 观测值类型 观测值和观测值精度 每一个有观测值的测站在文件中只能出现一次 没有设站的 已知点 如附和导线的定向点 和未知点 如前方交会点 在第二部分不必也不能给出任何虚拟测站信息 观测值分三 种 分别用一个字符 大小写均可 表示 L 表示方向 以 度分秒为单位 S 表示边长 以米为单位 A 表示方位角 以度分秒为单位 观测值精度与第一部分中的精度号相对应 若只有一组观测精度 则可省略 否则在观测值精度一栏中 须输入与该观测值对应的精度号 参见表 2 2 已知边长和 11 已知方位角的精度值一定要输 0 在同测站上的方向和边 长观测值按顺时针顺序排列 边角同测时 边长观测值最好 紧放在方向观测值的后面 如果边长是单向观测 则只需在一个测站上给出其边长 观测值 若是对向观测的边 则按实际观测情况在每一测站 上输入相应的边长观测值 程序将自动对往返边长取平均值 并作限差检验和超限提示 如果用户已将对向边长取平均值 则可对往返边长均输入其均值 或第一个边长 如往测 输 均值 第二个边长输一个负数如 1 对向观测边的精度高 于单向观测边的精度 但不增加观测值个数 平面观测文件中的测站顺序可以任意排列 一般来说不 会影响平差效率和结果 但本软件包还特意提供了观测值文 件排序 网点优化排序 的功能 详见 2 3 1 通过优化排序 既有利于网点近似坐标的推算 也可提高解算容量和速度 但一般对于 200 个点以上的大网或一些特殊网才有较明显的 效果 图 2 1 为某一测角网的网图 其相应的平面观测文件 IN2 的数据格式见表 2 1 对于有多组测角 测边精度的 网 其平面观测文件如表 2 2 所示 12 表 2 1 IN2 文件示例 仅一组精度的情况 0 7 3 3 1 3730958 610 264342 591 2 3714636 8876 276866 0832 1 2 L 0 3 L 27 362557 6 L 83 435791 2 4 L 0 3 L 74 593577 1 L 105 481560 4 5 L 0 3 L 41 334905 2 L 77 283653 5 6 L 0 1 6 3 2 5 4 图 2 1 13 3 L 58 405347 4 L 155 514999 6 1 L 0 3 L 57 240198 5 L 117 072390 3 1 L 0 2 L 121 345421 4 L 190 403024 5 L 231 554475 6 L 293 313088 表 2 2 1 800 3 000 2 000 1 3 000 5 000 3 000 2 5 000 5 000 5 000 3 k1 2800 000000 2400 000000 k4 2400 000000 3200 000000 k1 k2 L 0 0000 1 k5 L 44 595993 1 k6 L 89 595993 1 k7 L 135 000120 1 k4 p5 L 0 0000 2 p5 S 200 004728 2 p3 L 90 000031 2 14 2 1 22 1 2 高程观测文件高程观测文件 高程观测文件也是标准的 ASC 码文件 它的结构如下 已知点点号 已知点高程值 测段起点 终点 高差 距离 测段测站数 精度号 该文件的内容也分为两部分 第一部分为高程控制网的 已知数据 即已知高程点点号及其高程值 见文件的第 部分 第二部分为高程控制网的观测数据 它包括测段的起点点号 终点点号 测段高差 测段距离 测段测站数和精度号 见文 件的第 部分 第一部分中每一个已知高程点占一行 已知高程以米为 单位 其顺序可以任意排列 第二部分中每一个测段占一行 对于水准测量 两高程点间的水准线路为一测段 测段高差 以米为单位 测段距离以公里为单位 对于光电测距三角高 程网 测段表示每条光电测距边 测段距离为该边的平距 单位公里 如果平差时每一测段观测按距离定权 则 测 段测站数 这一项不要输入或输入一个负整数如 1 若输了 测站测段数 则平差时自动按测段测站数定权 该文件中测 段的顺序可以任意排列 当只有一种精度时 精度号可以不 15 输 对于多种精度 多等级 的水准网 第一部分的前面还 要增加几行 每行表示一种精度 有三个数据 即 水准等级 每公里精度值 单位 mm km 精度号 A2 A3 表 2 3 TP1 100 Z5 TP1 0 0585 1 000 Z5 Z6 0 0683 1 000 Z6 Z5 0 0634 1 000 Z6 A4 0 0683 1 000 Z6 Z7 0 0489 1 000 TP1 A2 0 0320 1 000 TP1 A3 40 1607 1 415 TP1 A4 0 0562 1 000 A3 TP1 39 8801 1 415 A2 TP1 0 0732 1 000 A4 Z6 0 0780 1 000 TP1 A4 Z5 Z6 Z7 图 2 2 16 A4 TP1 0 0683 1 000 图 2 2 为某一水准网的水准路线图 其相应的高程观测 文件见表 2 3 下面给出按距离定权 按测站数定权和多等级水准网的 高程观测文件的数据结构 1 按距离定权 S0 219 9592 N2 212 5328 N1 S246 24 8433 0 612 N1 S0 62 8298 0 858 N1 N0 50 7066 0 525 N0 S2 34 7798 0 690 N0 N2 4 6745 0 183 2 按测站数定权 9568 30 9584 9568 1 96985 0 10670 4 9568 9567 3 02405 0 07920 2 9567 9566 0 29515 0 05200 2 9568 9584 1 97090 0 10480 4 9584 9585 1 63340 0 10280 2 17 3 多等级水准网 1 表示不按测站数定权 1 1 000 1 2 2 000 2 3 3 000 3 BM1 120 000000 BM3 140 000000 BM1 BM2 19 9942 20 000 1 1 BM2 BM3 40 0073 24 000 1 1 BM3 BM4 10 0314 30 000 1 2 BM2 BM5 30 0088 23 000 1 2 BM3 BM5 10 0314 27 000 1 3 2 1 32 1 3 控制网观测文件的生成控制网观测文件的生成 在 COSAWIN 系统中 控制网观测文件的生成有四种方法 一种是利用通用文本编辑软件人工建立 这适合手工记录野 外观测数据的情况 另一种是在野外利用 COSA 系统的子系统 COSA HC 自动采集数据 传输到微机 通过 工具 栏的 格式转换 功能 全自动化地形成相应的平面和高程观测 文件 文件的第 部分需要人工编辑 第三种方法是人工建 立 SV 文件 调用 工具 栏中的斜距化平功能 自动 生成平面观测文件 参见 3 5 第四种方法是模拟生成平面 18 和高程观测文件 参见第五章 第一种方法一般需要用户对 方向 边长观测值作预处理 如方向改化 倾斜改正 归化 投影改正等 文件中的观测值必须和已知点坐标协调一致 对应于某一高程面和坐标系 第二种方法中 COSA HC 已对 边长作气象改正 仪器常数改正和倾斜改正 对高差已作曲 率和折光改正 但未作方向改化 边长归化投影改正 如需 要进行上述改正计算 可用本程序所提供的概算功能 参见 2 3 1 1 COSAWIN 系统中文件名的约定规则系统中文件名的约定规则 COSAWIN 系统处理对象是控制网 每一个控制网都用一 个字符串来作为控制网名 与某一控制网相关的所有文件 其文件名约定为以控制网名作为主文件名 用不同的后缀来 表示该网的不同类型文件 如对平面控制网的观测文件 其 文件名为 网名 IN2 对高程网观测文件 其文件名为 网 名 IN1 对平面网的平差结果文件 其文件名为 网名 OU2 对高程控制网的平差结果文件 其文件名为 网名 OU1 平面网及高程网的模拟方案文件报表输入输出文件以 及坐标转换文件的扩展名详见第五 七章 19 2 手工建立控制网观测文件 手工建立控制网观测文件 在 WINDOWS95 98 下 启动 COSAWIN 屏幕窗口弹出 如图 1 1 所示的主窗口 用鼠标单击 文件 主菜单中弹 出如图 1 2 所示的下拉菜单 选择其中的 新建 主菜单窗 口便会弹出一个文本编辑窗口 您可以按照 2 1 所述的文件 格式 输入控制网的起算数据和观测数据 输完各种数据后 存盘关闭文本编辑窗口 即完成了控制网观测文件的建立 3 自动建立控制网观测文件 自动建立控制网观测文件 自动建立控制网观测文件是建立在外业使用科傻系统的 子系统 COSA HC 进行外业观测数据自动采集的基础上 使用 COSAWIN 的 手簿通讯 和 数据转换 功能 直接形成控 制网观测文件 参见 3 6 和 3 7 2 22 2 控制网平差控制网平差 准备好控制网观测文件以后 即可进行平差处理 但在 平差前 一般还需要对平差过程中的某些参数进行设置 如 平差迭代限值 边长定权公式 精度评定时是使用先验单位 权中误差还是后验中误差 是否作网点优化排序 是否作观 测值概算 是否设置用边长交会推算网点近似坐标等 设置 20 是通过 COSAWIN 的 平差 主菜单下的 设置 来完成 具 体操作过程详见 2 3 如果控制网的范围较小 高程变化也 较小 且为独立的工程坐标系 已知点的 Y 坐标值较小 如 在 10 公里以下 或者平面观测文件中的观测值已经经过了 各种归化改正 则可直接进行平差处理 如果控制网的已知 点坐标是国家 54 或 80 坐标系下的坐标 且 Y 坐标值较大 即测区离中央子午线较远 平面观测文件中的边长 方向 值也没有经过概算 则需要利用 COSAWIN 的概算功能对方向 和边长观测值进行三差改正以及归化和投影改正 然后才能 进行平差 参见 2 2 1 2 4 和 2 3 1 这里需要说明的是 COSA HC 子系统已对边长观测值作了加乘常数改正 气象改 正以及斜距化平计算 平差时 只需在主菜单中用鼠标单击 平差 则会弹出 下拉菜单 如图 1 3 所示 下面我们对平面网和高程网平差 进一步予以说明 2 2 12 2 1 平面网平差平面网平差 如果观测文件中的边长 方向观测值需要进行改化计算 则须先在 平差 栏的 设置与选项 中进行相应选择 并 在 平差 栏中激活 生成概算文件 参见图 1 3 和 2 4 21 形成概算用文件后 用鼠标单击 平差 栏中的 平面 网 或单击工具条中平差快捷键 主菜单窗口弹出如图 2 3 所示的对话框 在该对话框中选择并打开要进行平差的平面 观测值文件 将自动进行概算 组成并解算法方程 法方程 求逆和精度评定及成果输出等工作 平差结果存于平面平差 结果文件 网名 OU2 并自动打开以供查看 图 2 3 22 在平差过程中若出现迭代次数多且不收敛的情况 或出 现其它提示 平差不能继续进行 首先应检查平面观测文件 是否有大的错误 若平差结果文件的后验单位权中误差显著 偏大 例如是先验单位权中误差的 1 5 倍以上 则应怀疑观 测值可能含有粗差 对于观测值粗差 可以查看观测值改正 数的大小并调用 工具 栏中 闭合差计算 菜单项 检查 闭合差是否超限 参见 3 1 对于图形结构较好 多余观测 数较多的网 还可调用 粗差探测 功能 探测和剔除粗差 参见 4 1 2 2 2 高程网平差高程网平差 用鼠标单击 平差 栏中的 高程网 或单击工具条 中的快捷键 主菜单窗口将弹出如图 2 4 所示的对话框 在 该对话框中选择并打开要进行平差的高程观测值文件 将自 动进行高程网平差 精度评定及成果输出等工作 平差结果 存于高程平差结果文件 网名 OU1 中 并自动打开以供查 看 通过查看和分析后验单位权中误差值以及高差观测值的 图 2 4 23 改正数 可以判断观测值和平差结果的质量 同样也可以调 用 工具 栏中 闭合差计算 功能菜单 检查各水准环线 的闭合差是否超限 2 32 3 设置与选项设置与选项 在 平差 菜单栏中选择 设置与选项 则会弹出如图 3 3 所示的设置对话框 该设置包括 平差设置 坐标常数 和改正数 坐标系统设置以及粗差剔除设置 下面分别予以 介绍 2 3 1 平差设置平差设置 平差设置界面如图 3 3 所示 包括了三个开关选择框 两组单选按钮设置框和一个编辑框 开关选择框用来确定某 项功能的开或关 用鼠标单击左边的方框可以设置开关选择 框的开关状态 当方框中有 标识符时 则表示该选择 框处于 开 状态 否则为 关 状态 对于一组单选按钮 设置框 一次只能选中其中的某一项 选中项的左侧圆圈中 会出现一个黑点 1 1 观测值文件排序 观测值文件排序 当该选项处于选中 开 状态时 则表示平差前先要对 24 原始观测文件进行优化排序 否则表示平差前不排序 这项 选择一般适合于大网 点数 500 或特殊网 对于大型网 观 测值文件优化排序后 可以提高平差计算速度 此外 通过 该项选择 对于较复杂的网 点的近似坐标计算会有影响 如增减迭代计算次数 迭代收敛或不收敛等 因此是否选择 此项 可通过试算确定 图 3 3 25 2 2 观测值概算 观测值概算 当该选项处于 开 状态时 则表示在平差前首先要对 原始观测值进行概算 若要进行概算 需要首先在 平差 栏中点激 生成概 算文件 并对该文件作必要编辑 若不进行概算 则关闭该 项 3 3 近似坐标用边长交会 近似坐标用边长交会 当该选项处于 开 状态时 表示推算近似坐标用边长 前方交会 否则在推算近似坐标时不使用边长交会 这项选 择适用于只有少量方向的边角网或混合网 对于单纯的测边 网 必须打开该项 否则网点近似坐标推算将不能进行 这里需要说明的是 由于边长交会的二义性 当交会某 一点的边只有两条时 交会出的点可能是错误的 这时可以 采用两种方法加以解决 一是建立一个网形信息文件 文件 名为 网名 NET 该文件为标准 ASCII 文件 可以使用任 意文本编辑器形成 其格式为 点名 1 点名 2 点名 3 点名 1 点名 2 点名 3 为边长交会三角形的三个顶点 按逆时针方向排列 每一个三角形组合占一行 二是建立一 个交会点的概略坐标文件 文件名为 网名 XYO 其格式为 26 点名 概略坐标 X0 概略坐标 Y0 概略坐标可以很粗糙 且只需要有二义性的交会点 为 了避免上述问题 布设纯测边网时 最好不要采用单三角形 应多增加跨三角形的长边 每个网点至少有三条边通过 这 样可减少边长交会的二义性 4 4 单位权选择 单位权选择 该选项是用来设置系统在进行精度评定时是使用先验单 位权中误差还是使用后验单位权中误差 用鼠标单击 先验 单位权 按钮 则设置使用先验单位权中误差 用鼠标单击 后验单位权 按钮 则设置使用后验单位权中误差 当多 余观测数较多时 使用后验单位权中误差较好 当多余观测 数很少 例如小于 8 则用先验单位权中误差为宜 在平差 结果文件 网名 OU2 中的最未部分 有先验和后验单位权 中误差信息 若两者相差较大 对于边角网或有多组精度的 网 已知坐标或观测值中可能含有粗差 或边角精度不匹配 若后验单位权特别大 则首先应怀疑观测值文件有错误 或 者近似坐标推算出错 5 5 边长定权公式 边长定权公式 该选项是用来设置系统在平差时采用什么公式来确定边 27 长观测值的中误差 本系统中提供了两种边长定权公式 一 种是按下式计算边长的中误差 另一种计算边长中误差的公式为 式中 A B 分别为测距仪的固定误差和比例误差 取自 网名 in2 文件 S 为边长值 单位为公里 由于边长定 权公式不同 平差结果有一定差别 可以用 工具 中的 叠置分析 进行比较 系统的缺省设置是后一种定权公式 6 6 平差迭代限值 平差迭代限值 平差迭代限值是平差迭代计算中最大的坐标改正数限值 COSAWIN 系统的缺省值为 10 厘米 若需要改变此项设置 可以直接在在编辑框中输入所要设定的值 当最大坐标改正 数小于限值时 停止迭代 进行平差精度评定 对于精度要 求很高的网 可设置小一些 如 1 厘米 如果平差迭代计算 中最大坐标改正数很大且不收敛 则应考虑观测文件的数据 有错 或推算近似坐标出错 设置完上述相应的选项后 用鼠标单击 确认 或 应 用 A 按钮 则接受更改的设置 否则单击 取消 则放弃 更改的设置而保持以前的设置 以后各项设置的确认和取消 SBA 222 SBA 28 操作和此处相同 2 3 2 坐标参数和改正数坐标参数和改正数 在如图 3 3 所示的对话框中 用鼠标击 坐标参数和改 正数 标签 则弹出坐标参数和改正数设置属性表 界面如 图 3 4 所示 29 1 1 坐标加常数 坐标加常数 用来设置 X Y 坐标的百公里加常数 以百公里为单位 当已知点在国家坐标系中且 X 坐标很大 一般为数千公里 可以将整百公里以上的值作为 X 坐标的加常数 以使计算简 便 Y 坐标若加 500 公里 则 Y 坐标的加常数为 500 公里 对于局部坐标系下的工程网 X Y 坐标均比较小 它们的加 常数为零 本系统的缺省值 图 3 4 30 2 2 水准尺每米真长改正数 水准尺每米真长改正数 用来设置水准尺的实际检测长度相对于每米标称尺长的 改正数 在水准网平差中用以对水准测段高差进行改正 缺 省设置值为零 2 3 3 坐标系统设置坐标系统设置 在如图 3 3 所示的对话框中 用鼠标击 坐标系统设置 标签 则弹出坐标系统设置属性表 其界面如图 3 5 所示 该功能在概算和坐标转换使用到 1 1 坐标系统 坐标系统 组合框 坐标系统 中的单选按钮用来设置投影改正时 采用的椭球系统 用户可以从系统所提供的北京 54 坐标椭球 国家 80 坐标椭球 WGS 84 坐标椭球以及用户自定义椭球系 统四种选项中任选一种 当选择用户自定义椭球时 需要用 户自己定义椭球的长半轴和椭球的扁率分母 长半轴输入值 以米为单位 2 2 中央子午线 中央子午线 中央子午线编辑框用来输入测区所处投影带中央子午线 31 的经度值 以度为单位 3 3 投影面高程 投影面高程 在该编辑框中 用户可以输入控展网所要投影到的高斯 投影面的高程值 以米为单位 在 COSAWIN 中的缺省设置是 坐标系统为北京 54 坐标 系统 中央子午线为 111 度 投影面高程为零 用户应特别 注意平面控制网中的已知点坐标是属于哪个坐标系统 哪个 图 3 5 32 投影面 对于独立坐标系下的工程控制网 关键是投影高程 面的正确确定 2 3 4 角度与距离单位设置角度与距离单位设置 角度可设置为 360 度制或 400 度制 欧洲国家采用 缺 省值为 360 度制 2 3 5 粗差剔除设置粗差剔除设置 用于设定方差比及粗差倍值 改正数倍值两个阈值 参 见图 3 6 方差比取值一般为 1 05 1 20 之间 愈小 则 所能探测到的粗差愈小 如 3 10 倍的观测误差 但方差比 不要小于 1 粗差倍值取值应大于 3 改正数倍值可取 1 5 2 5 之间的值 该值愈小 计算工作量愈大 2 42 4 生成生成概算用文件概算用文件 当平面观测文件中的方向 边长观测值需要概算时 调 用此功能可自动生成概算用文件中的近似坐标部分 若有该 网的高差观测值文件 则同时生成近似高程 自动生成概算所需要文件 网名 XYH 后 还须人工编 辑该文件 如添加每个网点的近似高程 对于精度要求较高 的平面网 还要输入每个测站上的觇标高 大地水准面差距 33 以及垂线偏差的子午分量和卯酉分量等附加量 以便对观测 值作三差改正 对于一般网 可不考虑三差改正 不必输后 四项 即上述附加量均自动作 0 处理 概算文件是一个标准 ASCII 格式文件 其结构如下 网点点名 X Y H BH N 该文件中每一网点的概算信息占一行 其中 X Y H 为 该点的近似三维坐标 BH 为觇标高 以米为单位 N 为大地 水准面差距 以米为单位 和 为该点垂线偏差的子午和 图 3 6 34 卯酉分量 以秒为单位 只有在网的精度要求较高时才需要 输入后面的四项 但前面三项是必须有的 当只有前三项时 后四项自动作 0 处理 实例 参见 概算算例网 XYH 文件 1 3730958 6100 264342 5910 535 7 8 0 1 3 3 4 0 7 2 3714636 8876 276866 0832 203 9 8 0 1 1 4 0 2 5 4 3700347 1407 266213 9081 376 0 8 1 1 6 0 6 3 7 5 3709621 8715 258215 6696 182 0 7 9 1 1 1 0 1 4 6 3721646 7827 254621 4564 166 2 7 8 1 0 0 8 3 4 3 3718773 3604 266467 5123 179 6 7 9 1 5 2 8 1 4 2 52 5 附加信息文件附加信息文件 对平面网 若人工建立 网名 IFI 文件 则可生成 网名 IFO 附加信息文件 网名 IFI 的结构为 R PRECISION 点对 1 QXX MATRIX 点号 1 点号 K L CONSTANT 实例 参见 隧道网 IFI R PREcision 2 3 13 31 35 13 21 13 171 21 171 qxx matrix 2 3 13 31 21 171 L constant 附加信息包括给定点对的相对误差信息 给定点的协因 数阵 每个观测值误差方程式的常数项 用于检查粗差 实 例参见 隧道网 IFO 36 第三章第三章 工具工具 COSAWIN 是一套测量控制网通用数据处理软件包 它不 仅能完成任意测量控制网常规的平差解算和精度评定等工作 还提供了一些非常有用的辅助功能 如平面 高程网闭合差 计算 贯通误差影响值计算 网图显绘 叠置分析 手簿通 讯和格式转换等功能 下面分别予以介绍 在主菜单中用鼠标单击 工具 则会弹出下拉菜单 如 图 1 5 所示 3 13 1 平面闭合差计算平面闭合差计算 启动该功能 首先由平面观测文件自动寻找出导线线路 DX 和最小独立闭合环线路 多边形 DBX 生成计算闭合 差线路文件 网名 CLI CLI 文件的结构如下 1 DX F3 D4 P1 P2 P3 F2 F1 2 DX 37 1 DBX D6 S3 S4 2 DBX D1 S1 S5 3 DBX 根据闭合差线路文件 自动计算出导线和多边形的角度 闭合差 坐标闭合差和全长闭合差 并进一步根据多边形角 度闭合差计算方向观测值精度 对闭合差进行评价和超限提 示 计算结果存放于闭合差结果文件 网名 CLO 中 操作 步骤如下 在如图 1 5 所示的下拉菜单中用鼠标单击 闭合差计算 中的 平面网 系统弹出打开文件对话框 提示您选取平面 观测值文件 在该对话框中选择并打开要进行闭合差计算的 平面观测值文件 则自动进行附合导线和最小闭合环的搜索 及闭合差的计算 用户也可按 CLI 文件格式手工建立闭合差 计算的线路文件 例如无定向导线 个别附合导线 若已存 在该网的 CLI 文件 则可直接利用 由于自动寻找附合导线 和最小独立闭合环线路是一个较为耗时的操作 对于点数较 多或较复杂的网 需要等待一段时间 当计算完成后 结果 38 存储在 网名 CLO 中 并自动打开以供查看 3 23 2 高程闭合差计算高程闭合差计算 根据高程观测文件自动寻找出水准 高程 附合线路和 最小独立闭合环线路 存放于闭合差线路文件 网名 GCI 中 根据闭合差线路文件 自动计算附合线路和多边形闭合 环的高程闭合差并作超限提示 同时根据闭合环的闭合差计 算每公里水准 高程 观测值的全中误差 计算结果存放于 闭合差结果文件 网名 GCO 中 其计算步骤同平面网闭合 差计算 3 33 3 贯通误差影响值贯通误差影响值计算计算 该功能主要为隧道施工控制网而设计的 其实质是根据 控制网的洞口点和定向点精度 贯通点的位置以及贯通面的 方向 在完成网平差之后 直接估算隧道贯通误差影响值 为此首先人工建立一个贯通误差引导文件 该文件也是一个 标准的 ASC 文件 命名规则为 网名 GTI 其格式为 进口点号 进口定向点号 出口点号 出口定向点号 贯通 点号 X 坐标 Y 坐标 贯通面方位角 39 A A F F J J G G E E B B C C 如上图所示的贯通方案可建立如下贯通文件 J A C E G XG YG J A C F G XG YG J B C E G XG YG J B C F G XG YG 表 3 1 为某工程网计算贯通误差引导文件实例 其中 171 21 和 100 分别是进口点 出口点和贯通点的点号 参 见 隧道网 GTI 表 3 1 171 8 21 12 100 69140 439260 90 00 171 8 21 11 100 69140 439260 90 00 171 5 21 12 100 69140 439260 90 00 171 5 21 11 100 69140 439260 90 00 171 9 21 12 100 69140 439260 90 00 171 9 21 11 100 69140 439260 90 00 为了选取最优的定向点方案 在一次计算中 可准备多 种不同的进出口点与定向点的组合 每一种组合占一行 准 备好引导文件 网名 GTI 后 用鼠标单击 贯通误差影响 值计算 将自动计算贯通误差影响值 并将结果存放在文件 40 网名 GTO 中 参见 隧道网 GTO 注意 在计算贯通误差影响值之前 除准备好引导文件 外 还须对该控制网进行平差计算 否则贯通误差计算失败 未建立贯通误差引导文件时 将弹出提示建立贯通误差引导 文件的对话框 3 43 4 图形显绘图形显绘 显绘平面网网图 单击 工具 栏中的 网图显绘 或单击工具条中的快捷键 主菜单窗口弹出选择网图信息文件对话框 在该对话框中选择并打开所需要的网图显 绘文件 网名 MAP 该文件是在对控制 网平差时自动形成的 则会自动在窗口显 绘该控制网的网图 用鼠标单击主菜单窗 口 查看 弹出如图 3 1 所示的下拉菜 单 同时在工具条中 对网图操作的一些 工具按钮也被激活 参见图 1 6 可对网 图进行包括放大 缩小 窗口放大 恢复 前级以及误差椭圆的显绘 控制网点的显绘和按等比例尺还 是变比例尺显绘等功能的操作 其中变比例尺显绘功能主要 用来放大隧道网的横向显示范围 单击工具条中的 打印 快捷键 可从打印机输出网图 应预先设置好打印机 图 3 1 41 3 53 5 斜距化平斜距化平 当采用人工记录且希望用该软件作观测值改正计算和概 算时 可调用此功能 该功能可将 网名 SV 文件自动生成 网名 IN2 文件 网名 SV 文件的结构如下 方向中误差 秒 测距仪固定误差 mm 比例误差 ppm 测距仪加常数 MM 乘常数 如 1 0 测距仪光轴和经 纬仪视准轴垂直偏距 mm 测距仪气象改正系数 K1 K2 大气折光系数 K 如 0 13 地球平均曲率半径 m 气压单位 可为 mmHg 或 mb 已知点点号 X 坐标 Y 坐标 海拔高程 H 测站点名 仪器高 m 温度 气压 照准点名 观测值类型 观测值 该文件可分为三部分 第一部分为精度和气象数据 占 文件的前四行 第二部分为已知点信息 点号 X 坐标 Y 坐 标 高程 H 每点占一行 如果该点无相应数据 则用 X Y H 代替 第三部分为测站观测数据 观测值类型为 L 时 表示为方向观测值 观测值类型为 S 时 表示斜距 天 顶距等观测值 且依次为 斜距 m 天顶距 度分秒 目 标高 m 对于整体型全站仪来说 测距仪光轴和经纬仪视准 轴垂直偏距一般为 0 只有当使用组合式全站仪即经纬仪加测 42 距头的情况才给出该垂直偏距值 表 3 2 为某一网的 SV 文件的示例 参见 斜距改化文件 SV 表 3 2 3 53 5 3 0 0 5 0 0 282 3868 0 13 6371 mmhg P1 5000 5000 H P2 5020 5300 180 P3 5046 6 5699 13 179 919 P4 5172 05 5970 5 H Z1 X Y 184 575 P2 1 5 20 0 760 P1 L 0 S1 L 180 0003 S1 S 100 003 90 00004 1 52 S1 1 4 22 760 P2 L 0 S2 L 180 0003 3 63 6 手簿通讯手簿通讯 激活该菜单项 将弹出如图 3 2 所示的对话框 将微机 与 COSA HC 的电子手簿相连接 可对原始观测数据文件 中 间文件文件以及平差结果文件进行双向转输 具体操作步骤 为 43 a 在电子手簿上选择 数据通讯 b 运行 COSAWIN 中的 手簿通讯 c 用 选择串行口 菜单设置通讯端口 d 使用 接收数据 菜单从手簿接收数据或使用 发 送数据 菜单从微机向手簿发送数据 接受或发送数据 时要根据弹出的对话框打开文件和输入存盘文件名 数 据通讯过程将在屏幕上显示 3 73 7 格式转换格式转换 可将由电子手簿通讯到微机的平面网 高程网或三维网 的原始观测值文件转换为平差所需要的平面观测和高程观测 文件 网名 IN2 和 网名 IN1 对于转换后的文件还需 图 3 2 44 作少量的人工编辑 如键入已知点坐标 已知点高程等 3 83 8 叠置分析叠置分析 激活该菜单项 将出现图 3 3 所示的对话框 选择所要 比较的两个文件 将得到同名点平差坐标的差值及其限差 对同一个平面网的不同观测方案 增减观测值 改变观 测精度 或不同设置 如采用不同的边长定权公式 进行平 差计算 调用该功能可得到两套成果平差后坐标之差别 例 如可以对一个网在粗差剔除前后或方差分量估计前后进行平 差结果比较 叠置分析的两种比较方案 其网名要有所区别 如比较同一网的 网名 1 OU2 和 网名 2 OU2 两个文件 图 3 3 是对某网概算前后进行叠置分析 其结果如下 点号 X m Y m x m y m p m 3 0 504 0 075 0 122 0 110 0 165 4 0 913 1 470 0 254 0 218 0 334 5 1 408 0 073 0 255 0 204 0 327 6 0 633 0 864 0 189 0 155 0 244 1 0 000 0 000 0 000 0 000 0 000 2 0 000 0 000 0 000 0 000 0 000 45 图 3 3 46 第四章第四章 粗差剔除和方差分量估计粗差剔除和方差分量估计 4 14 1 粗差剔除粗差剔除 当网平差结果的后验单位权中误差显著偏大时 若不存 在近似坐标推算错误和已知点坐标有问题 应考虑观测值中 可能存在粗差 该功能在网平差前后都可进行 在 平差 菜单栏中点击 粗差剔除 菜单项即可激活粗差剔除操作 根据提示打开所选的 网名 IN2 文件 则开始进行粗差计 算 由于粗差剔除计算工作量较大 比平差花的时间要多一 些 运行 粗差剔除 前 需要在 平差 栏的 设置与选 项 中查看粗差探测的相关设置参数 若不存在粗差 则 网名 IN2 文件不会改动 若存在粗差 将出文件 网名 OUC 该文件不能与 网名 OU2 文件作严格对比 主要查 看相应的粗差信息 见表 4 1 及粗差探测过程 粗差的大小 定值 和所在的观测值 定位 程序将对粗差作自动改正 重新生成剔除粗差后的观测值文件 网名 C IN2 对该文件 平差 相应的结果存放在 网名 C OU2 文件中 47 表 4 1 粗差探测过程 Step From To Error M0 Ratio 0 2 8926 1 1 8075 2 5611 2 3 27 0081 2 1 2284 2 1652 2 3 29 3407 1 2 13 8693 3 1 0275 1 4293 2 3 29 3370 1 2 13 8684 11 10 5 7799 4 0 9741 1 1127 2 3 29 3233 1 2 13 8852 11 10 5 7700 8 171 3 5714 经多维粗差探测 发现有下列粗差 粗差探测结果 FROM TO TYPE VALUE m M sec cm V sec cm G Error sec cm 2 3 L 108 213500 1 80 10 39 29 3370 1 2 L 34 400010 1 80 5 23 13 8684 11 10 S 2650 6625 0 43 4 32 5 7799 48 4 24 2 方差分量估计方差分量估计 当存在两种 如边角网 导线网 或两种以上观测值 或存在两种以上观测值精度时 应作方差分量估计 方差分 量估计应在粗差剔除之后进行 在 平差 栏中点击 方差 分量估计 按提示打开所选的 网名 IN2 文件 即进入方 差分量估计 运行后生成 网名 OUF 的输出文件 与方差 分量有关的输出如表 4 2 所示 此外还生成了方差分量估计 后的平面观测文件 网名 F IN2 该文件中 观测值的精度 应更接近实际观测精度 平差随机模型得到了改善 对该文 件作平差 即为方差分量估计后的平差结果 存于 网名 F OU2 中 该文件与 网点 OUF 是一致的 表 4 2 方差分量估计结果之一 Name VPV sec R Sum 1 6 6622 6 6605 17 2 3 0584 3 0613 9 3 6 6516 6 6551 18 4 5 6686 5 6643 11 5 9 9353 9 9382 28 6 7 0238 7 0206 13 49 方差分量估计结果之二 Sort Proior Error stimated Error Distance m Dir Dis Dir Sec Dis mm Dir Sec Dis cm r1 1 8 1 31 r2 3 0 2 81 r3 5 0 5 34 s1 3 2 0 19 565 6841 s2 5 3 0 38 200 0036 s3 5 5 0 37 99 9946 50 第五章第五章 网的优化设计网的优化设计 5 15 1 生成正态标准随机数生成正态标准随机数 主菜单 设计 栏的下拉菜 单如图 5 1 所示 有三项子菜单 项 现说明于后 单击 生成随机数 将弹出 一对话框 要求您输入生成随机 数的相关参数 第一个参数用于控制生成不相同的随机数序 列 其取值可取 1 10 的任意整数 第二个参数即 随机数个 数 只能选 200 400 或 500 三个值 即最