COSA（科傻）-CODAPS软件说明包教学版

科傻系统系列软件之二 （CODAPS 2003）

地 面 测 量 工 程

控制测量数据处理通用软件包

Version 5.0

武 地 课 题 组

2003.1．武 汉 前 言 .................................. 1 第一章 概述 .............................. 3

1.1 系统简介 .......................... 3 1.2 安装及运行 ........................ 3 1.3 快速入门 .......................... 5

第二章 平差 ............................. 10

2.1 控制网观测值文件 .................. 10 2.2 控制网平差 ....................... 20 2.3 设置与选项 ....................... 23

2.4 生成概算用文件 .................. 32

2.5 附加信息文件 ..................... 33

第三章 工具 ............................. 34

3.1 平面闭合差计算 .................. 34

3.2 高程闭合差计算 .................... 36 3.3贯通误差影响值计算 ................. 36

3.4图形显绘 ........................ 38

3.5斜距化平 ........................ 38 3.6手簿通讯 ........................ 40

3.7格式转换 ......................... 41 3.8叠置分析 ......................... 41

第四章 粗差探测、剔除和方差分量估计 ......... 43

4.1粗差探测与剔除 .................... 43

4.2方差分量估计 ...................... 48

第五章 网的模拟计算和优化设计 .............. 50

5.1生成正态标准随机数 ................. 50 5.2网的模拟计算 ...................... 50 5.3平面网优化设计 .................... 58

第六章 报表输出 ......................... 60

6.1原始数据报表 ...................... 60 6.2 平差结果报表 ..................... 67

第七章 坐标转换 ......................... 70

7.1 XYZ-〉BLH....................... 70 7.2 BLH->XYZ ....................... 71 7.3 XY->BL ......................... 72 7.4 BL->XY ......................... 73 7.5 XY1->XY2 ....................... 73 7.6 XY1->XY2 ....................... 74

附录2 CODAPS的文件组织 ................. 76

1平面控制网 ......................... 76 2 水准（高程）网 ...................... 77 3 GPS网 ............................ 78

附录3 所附实例文件目录 .................. 80 附录4 有关参考文献 ..................... 83 附录5 有关获奖情况 ..................... 87

3

前 言

“测量工程控制与施工测量内外业一体化和数据处理自动化系统”简称科傻系统(COSA),意即用高科技集成的傻瓜式测量数据采集和处理系统。它将测量学基本原理和现代科技相结合，对电子全站仪（TPS）、电子水准仪(ELI)、全球定位系统(GPS)以及常规地面测量仪器进行了系统的开发，以工程控制测量、施工测量和碎部测量等测量工程为对象，实现了从外业数据采集、质量检核、预处理到网平差等内业数据处理、成果报表输出的一体化和自动化作业流程。

科傻系统由以下三个子系统组成：

1） “基于掌上型电脑的测量数据采集和处理系统”（简称EREPS或COSA1）,在商品化掌上型电脑PDA上运行,该子系统灵活方便,适合外业环境。它能自动控制和引导整个作业过程并进行质量检测，一体化程度高,操作方便。该子系统具有水准测量、二三维控制测量、碎部测量、断面测量、道路和桥梁测设、工程放样等外业测量作业模块；具有小规模水准网、二三维工程网的平差功能；具有文件管理和数据通信功能。可实现从数据采集、质量检核、预处理到网平差等作业流程的一体化和自动化。

2) “地面测量工程控制测量数据处理通用软件包”（简称CODAPS或COSA2）,在便携式或台式微机的WINDOWS环境下运行,既可独立使用,也可与COSA1联合使用。可对掌上型电脑传输过来的水准测量、二三维控制测量原始观测数据进行转换,完成从概算到高程（水准高程、三角高程）网、平面网平差的自动化数据处理流程,同时具有控制网优化设计、粗差探测与剔除、方差分量估计、闭合差计算、

1

贯通误差影响值估算、坐标转换、换带计算、网图显绘、报表打印以及叠置分析等功能。该子系统经历了从DOS 版到WINDOWS版的改版和WINDOWS版的多次改进、升级,目前已达到炉火纯青的程度。

3) “GPS工程测量控制网通用平差软件包”（简称COSAGPS）,在便携式或台式微机的WINDOWS环境下运行,能读取各种GPS接收机的基线向量解算文件，进行网的三维无约束平差和二维约束平差，并具有坐标转换、网图显绘、报表打印以及与地面观测边长联合平差等功能。

本手册是在对WINDOWS版作多次改进、升级的基础上专门为CODAPS用户编写的，若仍有疏漏和不当之处,敬请读者提出宝贵意见和批评指正。

武 地 课 题 组

2003.1．武 汉 2

第一章 概述

1.1 系统简介

科傻系统（COSA）是“测量工程控制与施工测量内外业一体化和数据处理自动化系统”的简称，包括CODAPS、EREPS和COSAGPS三个子系统。CODAPS在IBM兼容机上运行。

CODAPS系统除具有概算、平差、精度评定及成果输出等功能外，还提供了许多实用的功能，如网的模拟设计、网图显绘、粗差剔除、方差分量估计、贯通误差影响值计算及闭合差计算等。

该系统不同于其它现有控制网平差系统的最大特点是自动化程度高，通用性强，处理速度快，解算容量大。其自动化表现在通过和COSA子系统EREPS相配合,可以做到由外业数据采集、检查到内业概算、平差和成果报表输出的自动化数据处理流程;其通用性表现在对控制网的网形、等级和网点编号没有任何限制, 无须给出冗余的附加信息，就可以处理任意结构的水准网和平面网；其解算速度快，解算容量大表现在采用稀疏矩阵压缩存储、网点优化排序和虚拟内存等技术，在主频166MHZ的586微机上，解算800个点的平面和水准控制网不到1分钟；在具有30MB剩余硬盘空间的微机上，可以解算多达6000个点的平面控制网。

1.2 安装及运行

1.2.1 安装

运行CODAPS系统盘上的安装程序SETUP.EXE后，按提示操作即

3

可。安装程序自动在程序组增加“科傻系统”程序项。

1.2.2 运行

在并口上插入加密狗后，点击“开始”菜单栏“程序”组项，会出现“科傻系统”程序组，单击其中程序项“科傻系统”中的“地面控制测量数据处理系统”，进入如图1-1所示的系统界面，所有操作符合WINDOWS 操作规范。系统的功能参见附录一。

4

图1-1

1.3 快速入门

1.3.1 系统菜单

（1）文件

文件菜单的主要功能如图1－2所示：

新建：新建文本文件，如平面观测值文件等。

打开：打开任意文件。

打印设置：打印机设置，单击将打开打印机设置对话框。

图1-2

5

（2）平差

平差菜单的主要功能如图1－3所示：

平面网：对平面网进行平差。单击将打开“输入平面观测值文件”对话框，选择平面观测值文件进行平面网平差。

高程网：对水准（高程）网进行平差。单击将打开“打开”对话框，选择水准（高程）观测值文件进行高程平差。

粗差探测：自动探测平面网观测值中的粗差，若发现粗差则自动剔除之。

方差分量估计：对于平面网中一组或多组不同种类或/和精度观测值的情况，通过方差分量估计，可以使各组观测值的精度获得最佳估计，保证平差随机模型的正确性和结果的可靠性。

设置与选项：在概算、平差、粗差探测以及坐标转换前作相应的设置和选项。

生成概算文件：作概算时需要调用此项，然后再进行平差。

图1-3

6

图1-4

（3）报表

报表菜单的主要功能如图1－4所示：

平差结果：根据平面网或高程网平差结果文件自动生成平面或高程平差结果报表。

原始观测值：将掌上型电脑经数据通信所得到的原始观测值文件自动生成平面高程网或高程网的原始观测值报表。其详细说明见第六章。

（4）查看

打开或关闭工具栏和状态栏。 （5）工具

工具菜单的主要功能如图1－5所示，其详细说明见第三章。

图1-5

7

（6）设计

用于平面网、水准网和GPS 网的模拟计算和优化设计。 功能菜单如图1－6所示， 其详细说明见第五章。

图1-6 （7）坐标转换

主要用于高斯坐标的 正反算、换带计算和 坐标转换。主要功能 菜单如图1－7所示， 其详细说明见第七章。

图1-7

1.3.2．工具条

快捷按钮

工具条快捷按钮如图1－8所示，共有6类共18个快捷按钮,即：

文件类：新建，打开，保存文件。

编辑类：剪切，复制，粘贴功能（对活动文档有效）。

8

打印类：打印当前活动文档（包括图形）。 平差类：平面网平差，高程网平差。

平面网网图显绘类：图形放大、缩小，开窗放大，返回前级，切换误差椭圆、点名和比例显示。

关于类：显示系统关于信息。

图1-8

1.3.3．文件命名规则及结构

本系统的文件命名规则及结构，请参见附录2“CODAPS的文件组织”和附录3“所附实例文件目录”。凡本说明书中所涉及的须人工生成的文件的详细结构，均可从所附实例文件找到；自动 生成的文件的可方便地从程序运行中得到。

1.3.4．算法说明和深层次研究

本说明书在附录4中列出了“有关参考文献”, 供深层次研究之用。

9

第二章 平差

使用CODAPS最常用的操作就是进行控制网平差，该菜单下包括平面网、高程网平差，粗差探测，方差分量估计，设置与选项以及生成概算用文件等子菜单。本章主要介绍平面、高程控制网观测值文件的结构及生成，平面、高程网平差，设置与选项以及生成观测值概算文件等内容。粗差探测和方差分量估计详见第四章。

2.1 控制网观测值文件

在进行平差之前，必须要准备好控制网观测值文件，即平面观测值文件（取名规则为“网名.in2”）和高程观测值文件（取名规则为“网名.in1”）。观测值文件采用网点数据结构，除包含控制网的所有已知点、未知点和观测值信息外，还隐含了控制网的拓扑信息。

可以使用系统菜单中“文件”栏下拉“新建”子菜单项或单击工具栏左边第一个快捷键建立平面或高程观测值文件。

2.1.1 平面观测值文件

平面观测值文件为标准的ASCⅡ码文件，可以使用任何文本编辑器建立编辑和修改。其结构如下：

10

方向中误差1，测边固定误差1，比例误差1[，精度号1] 方向中误差2，测边固定误差2，比例误差2，精度号2

Ⅰ

?，?，?，?

方向中误差n，测边固定误差n，比例误差n，精度号n 已知点点号，X坐标，Y坐标 ?，?，?

测站点点号

Ⅱ 照准点点号，观测值类型，观测值[，观测值精度] ?，?，?[，?]

该文件分为两部分：第一部分为控制网的已知数据，包括先验

的方向观测精度，先验测边精度和已知点坐标(见文件的Ⅰ部分)；第二部分为控制网的测站观测数据(见文件的Ⅱ部分)，包括方向、边长、方位角观测值。为了文件的简洁和统一，我们将已知边和已知方位角也放到测站观测数据中，它们和相应的观测边和观测方位角有相同的“观测值类型”，但其精度值赋“0”，即权为无穷大。

第一部分的排列顺序为：第一行为方向中误差，测边固定误差，测边比例误差。若为纯测角网，则测边固定误差和比例误差不起作用；若为纯测边网，方向误差也不起作用，这时可输一个默认值“1”。程序始终将第一行的方向中误差值作为单位权中误差。若只有一种(或称为一组)测角、测边精度，则可不输入精度号。这时，从第二行开始为已知点点号及其坐标值，每一个已知点数据占一行。若有几种测角测边精度，则需按精度分组，组数为测角、测边中最多的精度种类数，每一组占一行，精度号输1、2、...（参见表2-2）。如两种测角精度，三种测边精度，则应分成三组。

11

方向中误差单位为秒，测边固定误差单位为毫米，测边比例误差单位为ppm。第一行的三个值都必须赋值，对于纯测角网，测边的固定误差和比例误差可输任意两个数值，如5,3；对于纯测边网，方向中误差赋为1.0。已知点点号（或点名，下同）为字符型数据，可以是数字、英文字母（大小写均可）、汉字或它们的组合（测站点，照准点亦然），X、Y坐标以米为单位。

第二部分的排列顺序为：第一行为测站点点号，从第二行开始为照准点点号，观测值类型，观测值和观测值精度。在同测站上的方向和边长观测值必须按顺时针顺序排列。边角同测时，边长观测值最好紧放在方向观测值的后面。每一个有观测值的测站在文件中只能出现一次。没有设站的已知点（如附和导线的定向点）和未知点（如前方交会点）在第二部分不必也不能给出任何虚拟测站信息。观测值分三种，分别用一个字符（大小写均可）表示：L—表示方向，以度分秒为单位。S—表示边长，以米为单位。A—表示方位角，以度分秒为单位。观测值精度与第一部分中的精度号相对应，若只有一组观测精度，则可省略；否则在观测值精度一栏中须输入与该观测值对应的精度号（参见表2-2）。已知边长和已知方位角的精度值一定要输“0”。

如果边长是单向观测，则只需在一个测站上给出其边长观测值。若是对向观测的边，则按实际观测情况在每一测站上输入相应的边长观测值，程序将自动对往返边长取平均值并作限差检验和超限提示；如果用户已将对向边长取平均值，则可对往返边长均输入其均值，或第一个边长（如往测）输均值，第二个边长输一个负数如“-1”。对向观测边的精度高于单向观测边的精度，但不增加观测值个数。

12

平面观测值文件中的测站顺序可以任意排列。一般来说不会影响平差效率和结果，但本软件包还特意提供了观测值文件排序(网点优化排序)的功能（详见2.3.1）。通过优化排序，既有利于网点近似坐标的推算，也可提高解算容量和速度，但一般对于200个点以上的大网或一些特殊网才有较明显的效果。

图2-1为某一测角网的网图，其相应的平面观测值文件*.IN2”的数据格式见表2-1。对于有多组测角、测边精度的网，其平面观测值文件如表2-2所示。

1 6

3 2

5 4 图 2-1

表 2-1 IN2文件示例（仅一组精度的情况）

0.7,3,3

1,3730958.610,264342.591 2,3714636.8876,276866.0832 1

13

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 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

14

表2-2 多组精度情况的IN2文件示例

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 ??

2.1.2 高程观测值文件

高程观测值文件也是标准的ASCⅡ码文件，它的结构如下：

Ⅰ 已知点点号，已知点高程值 ?，?

测段起点，终点，高差，距离，测段测站数，精度号

Ⅱ ????[，?]

该文件的内容也分为两部分，第一部分为高程控制网的已知数

15

据，即已知高程点点号及其高程值(见文件的第Ⅰ部分)。第二部分为高程控制网的观测数据，它包括测段的起点点号，终点点号，测段高差，测段距离、测段测站数和精度号(见文件的第Ⅱ部分)。

第一部分中每一个已知高程点占一行，已知高程以米为单位，其顺序可以任意排列。第二部分中每一个测段占一行，对于水准测量，两高程点间的水准线路为一测段，测段高差以米为单位，测段距离以公里为单位。对于光电测距三角高程网，测段表示每条光电测距边，测段距离为该边的平距（单位公里）。如果平差时每一测段观测按距离定权，则“测段测站数”这一项不要输入或输入一个负整数如-1。若输了测站测段数，则平差时自动按测段测站数定权。该文件中测段的顺序可以任意排列。当只有一种精度时，精度号可以不输。对于多种精度（多等级）的水准网，第一部分的前面还要增加几行，每行表示一种精度，有三个数据。

水准等级，每公里精度值（单位mm/km）,精度号。 A2 A3 TP1 A4 Z5 Z6 Z7 图 2-2

图2-2为某一水准网的水准路线图，其相应的高程观测值文件

16

见表2-3。

表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 A4,TP1,0.0683,1.000

下面给出按距离定权、按测站数定权和多等级水准网的高程观测值文件的数据结构

（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

17

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 ......

（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 ......

18

2.1.3 控制网观测值文件的生成

在CODAPS系统中，控制网观测值文件的生成有四种方法：一种是利用通用文本编辑软件人工建立，这适合手工记录野外观测数据的情况；另一种是在野外利用COSA系统的子系统EREPS自动采集数据，传输到微机，通过“工具”栏的“格式转换”功能，全自动化地形成相应的平面和高程观测值文件（文件的第Ⅰ部分需要人工编辑）第三种方法是人工建立“*.SV”文件，调用“工具”栏中的斜距化平功能，自动生成平面观测值文件（参见3.5）。第四种方法是模拟生成平面和高程观测值文件(参见第五章)。第一种方法一般需要用户对方向、边长观测值作预处理，如方向改化、倾斜改正、归化投影改正等，文件中的观测值必须和已知点坐标协调一致，对应于某一坐标系和投影高程面。第二种方法中，EREPS已对边长作气象改正、仪器常数改正和倾斜改正，对高差已作曲率和折光改正，但未作方向改化、边长归化及投影改正。如需要进行上述改正计算，可使用本程序所提供的概算功能（参见2.3.1）。 （1）CODAPS系统中文件名的约定规则

CODAPS系统处理对象是控制网，每一个控制网都用一个字符串来作为控制网名。与某一控制网相关的所有文件，其文件名约定为以控制网名作为主文件名，用不同的后缀来表示该网的不同类型文件。如对平面控制网的观测值文件，其文件名为“网名.IN2”，对高程网观测值文件，其文件名为“网名.IN1”；对平面网的平差结果文件，其文件名为“网名.OU2”，对高程控制网的平差结果文件，其文件名为“网名.OU1”。平面网及高程网的模拟观测方案文件、报表输入输出文件以及坐标转换文件的扩展名详见第五、七章。

19

（2）手工建立控制网观测值文件

在WINDOWS95/98下，启动CODAPS，屏幕窗口弹出如图1-1所示的主窗口。用鼠标单击“文件”， 主菜单中弹出如图1-2所示的下拉菜单。选择其中的“新建”，主菜单窗口便会弹出一个文本编辑窗口，您可以按照2.1所述的文件格式，输入控制网的起算数据和观测数据。输完各种数据后，存盘关闭文本编辑窗口，即完成了控制网观测值文件的建立。 （3）自动建立控制网观测值文件

自动建立控制网观测值文件是建立在外业使用科傻系统的子系统EREPS进行外业观测数据自动采集的基础上，使用CODAPS的“手簿通讯”和“数据转换”功能，直接形成控制网观测值文件，参见3.6和3.7。

2.2 控制网平差

准备好控制网观测值文件以后，即可进行平差处理。但在平差前，一般还需要对平差过程中的某些参数进行设置，如平差迭代限值，边长定权公式；精度评定时是使用先验单位权中误差还是后验中误差；是否作网点优化排序，是否作观测值概算，是否设置用边长交会推算网点近似坐标等。设置是通过CODAPS的“平差”主菜单下的“设置”来完成，具体操作过程详见2.3。如果控制网的范围较小，高程变化也较小，且为独立的工程坐标系，已知点的Y坐标值较小，或者平面观测值文件中的观测值已经经过了各种归化改正，则可直接进行平差处理。如果控制网的已知点坐标是国家54或80坐标系下的坐标，且Y坐标值较大（即测区离中央子午线较远），平

20

面观测值文件中的边长、方向值也没有经过概算，则需要利用CODAPS的概算功能对方向和边长观测值进行三差改正以及归化和投影改正计算，然后才能进行平差（参见2.2.1,2.4和2.3.1）。这里需要说明的是，EREPS子系统已对边长观测值作了加乘常数改正、气象改正以及斜距化平计算。

平差时，只需在主菜单中用鼠标单击“平差”，则会弹出下拉菜单，如图1-3所示。下面我们对平面网和高程网平差进一步予以说明。

2.2.1 平面网平差

如果观测值文件中的边长、方向观测值需要进行改化计算，则须先在“平差”栏的“设置与选项”中进行相应选择，并在“平差”栏中激活“生成概算文件”(参见图1-3和2.4)。

形成概算用文件后，用鼠标单击“平差”栏中的“平面网”或单击工具条中平差快捷键，主菜单窗口弹出如图2-3所示的对话框。在该对话框中选择并打开要进行平差的平面观测值文件，将自动进行概算、组成并解算法方程、法方程求逆和精度评定及成果输出等工作，平差结果存于平面平差结果文件“网名.OU2”，并自动打开以供查看。

21

图 2-3

在平差过程中若出现迭代次数多且不收敛的情况，或出现其它提示，平差不能继续进行，首先应检查生成的平面观测值文件是否有错误。若平差结果文件的后验单位权中误差显著偏大（例如是先验单位权中误差的2倍以上），则应怀疑观测值可能含有粗差。对于观测值粗差，可以查看观测值改正数的大小并调用“工具”栏中“闭合差计算”菜单项，检查闭合差是否超限（参见3.1）。对于图形结构较好，多余观测数较多的网，还可调用“粗差探测”功能，探测和剔除粗差（参见4.1）。

22

2.2.2 高程网平差

用鼠标单击“平差”栏中的“高程网”， 或单击工具条中的快捷键，主菜单窗口将弹出如图2-4所示的对话框。在该对话框中选择并打开要进行平差的高程观测值文件，将自动进行高程网平差、精度评定及成果输出等工作。平差结果存于高程平差结果文件“网名.OU1”中，并自动打开以供查看。通过查看和分析后验单位权中误差值以及高差观测值的改正数，可以判断观测值和平差结果的质量；同样也可以调用“工具”栏中“闭合差计算”功能菜单，检查各水准环线的闭合差是否超限

2.3 设置与选项

图2-4

在“平差”菜单栏中选择“设置与选项”，则会弹出如图3-3所示的设置对话框，该设置包括：平差设置、坐标常数和改正数、坐标系统设置以及粗差剔除设置，下面分别予以介绍。

23

2.3.1 平差设置

平差设置界面如图2-5所示，包括了三个开关选择框、两组单选按钮设置框和一个编辑框。开关选择框用来确定某项功能的开或关，用鼠标单击左边的方框可以设置开关选择框的开关状态，当方框中有“√”标识符时，则表示该选择框处于“开”状态，否则为“关”状态。对于一组单选按钮设置框，一次只能选中其中的某一项，选中项的左侧圆圈中会出现一个黑点。

图2-5

24

（1）观测值文件排序

当该选项处于选中（开）状态时，则表示平差前先要对原始观测值文件进行优化排序，否则表示平差前不排序。这项选择一般适合于大网(点数>500)或特殊网。对于大型网，观测值文件优化排序后，可以提高平差计算速度。此外，通过该项选择，对于较复杂的网，点的近似坐标计算会有影响，如增减迭代计算次数，迭代收敛或不收敛等。因此是否选择此项，可通过试算确定。 （2）观测值概算

当该选项处于“开”状态时，则表示在平差前首先要对原始观测值进行概算。

若要进行概算，需要首先在“平差”栏中点激“生成概算文件”，并对该文件作必要编辑。若不进行概算，则关闭该项。 （3）近似坐标用边长交会

当该选项处于“开”状态时，表示推算近似坐标用边长前方交会，否则在推算近似坐标时不使用边长交会。这项选择适用于只有少量方向的边角网或混合网，对于单纯的测边网，必须打开该项，否则网点近似坐标推算将不能进行。

这里需要说明的是：由于边长交会的二义性，当交会某一点的边只有两条时，交会出的点可能是错误的，这时可以采用两种方法加以解决，一是建立一个网形信息文件，文件名为“网名.NET”，该文件为标准ASCII文件，可以使用任意文本编辑器形成，其格式为：

点名1，点名2，点名3

点名1、点名2、点名3为边长交会三角形的三个顶点，按逆时

25

针方向排列，每一个三角形组合占一行。二是建立一个交会点的概略坐标文件，文件名为“网名.XYO”，其格式为：

点名 概略坐标X0 概略坐标Y0

概略坐标可以很粗糙，且只需要有二义性的交会点。为了避免上述问题，布设纯测边网时，最好不要采用单三角形，应增加跨三角形的边，每个网点至少有三条边通过，这样可减少边长交会的二义性。

（4）单位权选择

该选项是用来设置系统在进行精度评定时是使用先验单位权中误差还是使用后验单位权中误差，用鼠标单击“先验单位权”按钮，则设置使用先验单位权中误差；用鼠标单击“后验单位权”按钮，则设置使用后验单位权中误差。当多余观测数较多时，使用后验单位权中误差较好。当多余观测数很少（例如小于10），则用先验单位权中误差为宜。在平差结果文件“网名.OU2”中的最未部分，有先验和后验单位权中误差信息，若两者相差较大，对于边角网或有多组精度的网，已知坐标或观测值中可能含有粗差，或边角精度不匹配。若后验单位权特别大，则首先应怀疑观测值文件有错误，或者近似坐标推算出错。 （5）边长定权公式

该选项是用来设置系统在平差时采用什么公式来确定边长观测值的中误差，本系统中提供了两种边长定权公式，一种是按下式计算边长的中误差：

A?B?S另一种计算边长中误差的公式为：

A2?B2?S2 26

式中，A、B分别为测距仪的固定误差和比例误差，取自“网名.in2”文件，S为边长值，单位为公里。由于边长定权公式不同，平差结果有一定差别，可以用“工具”中的“叠置分析”进行比较。系统的缺省设置是后一种定权公式。 （6）平差迭代限值

平差迭代限值是平差迭代计算中最大的坐标改正数限值，CODAPS系统的缺省值为10厘米。若需要改变此项设置，可以直接在在编辑框中输入所要设定的值。当最大坐标改正数小于限值时，停止迭代，进行平差精度评定。对于精度要求很高的网，可设置小一些（如1厘米）。如果平差迭代计算中最大坐标改正数很大且不收敛，则应考虑观测值文件的数据有错，或推算近似坐标出错。

设置完上述相应的选项后，用鼠标单击〖确认〗或〖应用(A)〗按钮，则接受更改的设置，否则单击〖取消〗则放弃更改的设置而保持以前的设置（以后各项设置的确认和取消操作和此处相同）。

2.3.2 坐标参数和改正数

在如图2-5所示的对话框中，用鼠标击“坐标参数和改正数”标签，则弹出坐标参数和改正数设置表，界面如图2-6所示。

27

图2-6

（1）坐标加常数

用来设置X、Y坐标的百公里加常数，以百公里为单位。当已知点在国家坐标系中且X坐标很大（一般为数千公里），可以将整百公里以上的值作为X坐标的加常数，以使计算简便。若设定Y坐标的加常数为500公里，则Y坐标将加上500公里（Y坐标前未添加带号）。对于局部坐标系下的工程网，X、Y坐标均比较小，它们的加常数为零（本系统的缺省值）。 （2）水准尺每米真长改正数

用来设置水准尺的实际检测长度相对于每米标称尺长的改正

28

数，在水准网平差中用以对水准测段高差进行改正。缺省设置值为零。

2.3.3 坐标系统设置

在如图2-5所示的对话框中，用鼠标击“坐标系统设置”标签，则弹出坐标系统设置表，其界面如图2-7所示。该功能在概算和坐标转换使用到。 （1）坐标系统

组合框“坐标系统”中的单选按钮用来设置投影改正时采用的椭球系统，用户可以从系统所提供的北京54坐标椭球、国家80坐标椭球、WGS-84坐标椭球以及用户自定义椭球系统四种选项中任选一种。当选择用户自定义椭球时，需要用户自己定义椭球的长半轴和椭球的扁率分母，长半轴输入值以米为单位。 （2）中央子午线

中央子午线编辑框用来输入测区所处投影带中央子午线的经度 值，以度为单位。

29

图2-7

（3）投影面高程

在该编辑框中，用户可以输入控展网所要投影到的高斯投影面的高程值，以米为单位。

在CODAPS中的缺省设置是：坐标系统为北京54坐标系统，中央子午线为111度，投影面高程为零。用户应特别注意平面控制网中的已知点坐标是属于哪个坐标系统，哪条中央子午线，哪个投影高程面。对于独立坐标系下的工程控制网，关键是投影高程面的正确确定。

30

2.3.4 角度与距离单位设置

角度可设置为360度制或400度制（欧洲国家采用）缺省值为360度制。

2.3.5粗差剔除设置

用于设定方差比及粗差倍值、改正数倍值两个阈值（参见图2-8。方差比取值一般为1.05 - 1.20之间，愈小，则所能探测到的粗差愈小（如3 - 4倍的观测误差）。但方差比不要小于1，粗差倍值取值应大于3。改正数倍值可取1.5 - 2.5之间的值，该值愈小，计算工作量愈大。

图2-8

31

2.4 生成概算用文件

当平面观测值文件中的方向、边长观测值需要概算时，调用此功能可自动生成概算用文件中的近似坐标部分，若有该网的高差观测值文件，则同时生成近似高程。

自动生成概算所需要文件“网名.XYH”后，还须人工编辑该文件，如添加每个网点的近似高程，对于精度要求较高的平面网，还要输入每个测站上的觇标高、大地水准面差距以及垂线偏差的子午分量和卯酉分量等附加量，以便对观测值作三差改正。对于一般网，可不考虑三差改正，不必输后四项，即上述附加量均自动作“0”处理。

概算文件是一个标准ASCII格式文件，其结构如下：

网点点名，X，Y，H，BH，N，ζ，η 该文件中每一网点的概算信息占一行，其中X、Y、H为该点的近似三维坐标；BH为觇标高，以米为单位；N为大地水准面差距，以米为单位；ζ和η为该点垂线偏差的子午和卯酉分量，以秒为单位。只有在网的精度要求较高时才需要输入后面的四项。但前面三项是必须有的，当只有前三项时，后四项自动作“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

32

2.5 附加信息文件

对平面网，若人工建立“网名.IFI”文件，则可生成“网名.IFO”附加信息文件。“网名.IFI”的结构为：

R_PRECISION 点对1 ...... QXX_MATRIX

点号1，......，点号K L_CONSTANT

实例（参见“隧道网.IFI”文件） R_PREcision 2,3 13 31 13 21 13 171 21,171 qxx_matrix 2,3,13 31 21 171 L_constant

附加信息包括给定点对的相对误差信息，给定点的协因数阵，每个观测值误差方程式的常数项（用于检查粗差）。

33

第三章 工具

CODAPS是一套完整的测量控制网数据处理通用软件包，它不仅能完成任意测量控制网常规的平差计算和精度评定等工作，还提供了一些非常有用的辅助功能。如平面网、高程网的闭合差计算，隧道平面网的贯通误差影响值计算，网图显绘，叠置分析，手簿通讯和格式转换等功能。下面分别予以介绍。

在主菜单中用鼠标单击“工具”，则会弹出下拉菜单，如图1-5所示。

3.1 平面闭合差计算

启动该功能，首先由平面观测值文件自动寻找出线路（导线DX）和最小独立闭合环线路（多边形DBX），并自动生成计算闭合差的输入文件“网名.CLI”（该文件亦可人工编辑生成）。

CLI文件的结构如下： 1,DX F3 D4 P1 P2 P3 F2 F1 2,DX

34

...... 1,DBX D6 S3 S4 2,DBX D1 S1 S5 3,DBX ......

根据闭合差线路文件，自动计算出导线和多边形的角度闭合差、坐标闭合差和全长闭合差，并进一步根据多边形角度闭合差计算方向观测值精度，对闭合差进行评价和超限提示。计算结果存放于闭合差结果文件“网名.CLO”中。操作步骤如下：

在如图1-5所示的下拉菜单中用鼠标单击“闭合差计算”中的“平面网”，系统弹出打开文件对话框，提示您选取平面观测值文件。在该对话框中选择并打开要进行闭合差计算的平面观测值文件，则自动进行附合导线和最小闭合环的搜索及闭合差的计算。用户也可按CLI文件格式手工建立闭合差计算的线路文件（例如无定向导线，个别附合导线）。若已存在该网的CLI文件，则可直接利用。由于自动寻找附合导线和最小独立闭合环线路是一个较为耗时的操作，对于点数较多或较复杂的网，需要等待一段时间。当计算完成后，结

35

果存储在“网名.CLO”中，并自动打开以供查看。

3.2 高程闭合差计算

根据高程观测值文件自动寻找出水准（高程）附合线路和最小独立闭合环线路，存放于闭合差线路文件“网名.GCI”中，根据闭合差线路文件，自动计算附合线路和多边形闭合环的高程闭合差并作超限提示，同时根据闭合环的闭合差计算每公里水准（高程）观测值的全中误差，计算结果存放于闭合差结果文件“网名.GCO”中。其计算步骤同平面网闭合差计算。

3.3贯通误差影响值计算

该功能主要为隧道施工控制网而设计的，其实质是根据控制网的洞口点和定向点精度、贯通点的位置以及贯通面的方向，在完成网平差之后，直接估算隧道贯通误差影响值。为此首先人工建立一个贯通误差引导文件，该文件也是一个标准的ASCⅡ文件，命名规则为“网名.GTI”，其格式为：

进口点号，进口定向点号，出口点号，出口定向点号，贯通点号，X坐标，Y坐标，贯通面方位角 A J

如上图所示的贯通方案可建立如下贯通文件 J，A，C，E，G，XG，YG，α

B F E C G 36

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”后，用鼠标单击“贯通误差影响值计算”，将自动计算贯通误差影响值，并将结果存放在文件“网名.GTO”中（参见“隧道网.GTO”）。

注意：在计算贯通误差影响值之前，除准备好引导文件外，还须对该控制网进行平差计算，否则贯通误差计算失败。未建立贯通误差引导文件时，将弹出提示建立贯通误差引导文件的对话框。

37

3.4图形显绘

显绘平面网网图。单击“工具”栏中的“网图显绘”或单击工具条中的快捷键，主菜单窗口弹出选择网图信息文件对话框。在该对话框中选择并打开所需要的网图显绘文件“网名.MAP”（该文件是在对控制网平差时自动形成的），则会自动在窗口显绘该控制网的网图。用鼠标单击主菜单窗口“查看”，弹出如图3-1所示的下拉菜单，同时在工具条中，对网图操作的一些工具按钮也被激活（参见图1-6）。可对网图进行包括放大、缩小、窗口放大、误差椭圆显绘和控制网点显绘的关闭，和按等比例尺还是变比例尺显绘等功能的操作。其中变比例尺显绘功能主要

用来放大隧道网的横向显示范围。单击工具条中的“打印”快捷键，可从打印机输出网图（应预先设置好打印机）。

图3-1

3.5斜距化平

当采用人工记录且希望用该软件作观测值改正计算和概算时，可调用此功能。该功能可将“网名.SV”文件自动生成“网名.IN2”文件。“网名.SV”文件的结构如下：

方向中误差(秒)，测距仪固定误差(mm),比例误差(ppm) 测距仪加常数(MM)，乘常数(如1.0)，测距仪光轴和经纬仪视准轴垂直偏距(mm)

测距仪气象改正系数K1，K2，大气折光系数K(如0.13)

38

地球平均曲率半径(m)，气压单位(可为mmHg 或 mb) 已知点点号，X坐标，Y坐标，海拔高程H 测站点名

仪器高(m)，温度(℃)，气压 照准点名，观测值类型，观测值

该文件可分为三部分，第一部分为精度和气象数据，占文件的前四行。第二部分为已知点信息（点号，X坐标，Y坐标，高程H），每点占一行，如果该点无相应数据，则用X，Y，H代替。第三部分为测站观测数据，观测值类型为L时，表示为方向观测值；观测值类型为S时，表示斜距、天顶距等观测值，且依次为：斜距(m)，天顶距(度分秒)，目标高(m)。对于整体型全站仪来说，测距仪光轴和经纬仪视准轴垂直偏距一般为0；只有当使用组合式全站仪即经纬仪加测距头的情况才给出该垂直偏距值。表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

39

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.6手簿通讯

激活该菜单项，将弹出如图3-2所示的对话框，将微机与EREPS

图3-2

40

的电子手簿相连接，可对原始观测数据文件，中间文件文件以及平差结果文件进行双向转输。具体操作步骤为：

(a) 在电子手簿上选择“数据通讯”； (b) 运行CODAPS中的“手簿通讯”； (c) 用“选择串行口”菜单设置通讯端口；

(d) 使用“接收数据”菜单从手簿接收数据或使用“发送数据”菜单从微机向手簿发送数据，接受或发送数据时要根据弹出的对话框打开文件和输入存盘文件名，数据通讯过程将在屏幕上显示。

3.7格式转换

可将由电子手簿通讯到微机的平面网、高程网或三维网的原始观测值文件转换为平差所需要的平面观测和高程观测值文件“网名.IN2”和“网名.IN1”。对于转换后的文件还需作少量的人工编辑，如键入已知点坐标、已知点高程等。

3.8叠置分析

激活该菜单项，将出现图3-3所示的对话框。选择所要比较的两个文件，将得到同名点平差坐标的差值及其限差。

对同一个平面网的不同观测方案（增减观测值，改变观测精度）或不同设置（如采用不同的边长定权公式）进行平差计算，调用该功能可得到两套成果平差后坐标之差别。例如可以对一个网在粗差剔除前后或方差分量估计前后进行平差结果比较。叠置分析的两种比较方案，其网名要有所区别，如比较同一网的“网名1.OU2”和“网名2.OU2”两个文件。

41

图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

图3-3

42

第四章 粗差探测、剔除和方差分量估计

4.1粗差探测与剔除

4.1.1．平面网的粗差探测与剔除

当网平差结果的后验单位权中误差较先验单位权中误差显著偏大（如5-10倍）时，若不存在输入文件错误、近似坐标推算错误（若存在这类错误，则网的平差结果文件和后验单位权中误差将出现明显出错显示）和已知点坐标有问题，应考虑观测值中可能存在粗差。

该功能只有在网的图形强度好（如边角全测的由三角形或/和大地四边形组成的施工控制网和变形监测网）、有足够多余观测（如多余观测数大于20且平均多余观测分量大于0.4）时才有效和可靠。

对其它的图形强度一般的网，粗差探测的结果可作参考，还应与平面网的闭合差计算结果相结合进行分析，以确定粗差观测值。请参见本说明书附录4的“有关参考文献”。

运行“粗差剔除”前，需要在“平差”栏的“设置与选项”中查看粗差探测的相关设置参数。在 “平差”菜单栏中点击“粗差探测”下的“平面网”菜单项，根据提示打开所选的“网名.IN2”文件，则开始进行粗差计算。由于粗差剔除计算工作量较大，比平差花的时间要多一些。若不存在粗差，则“网名.IN2”文件不会改动；若存在粗差，将出文件“网名.OUC”。请注意：该文件的其它结果不能与“网名.OU2”文件作严格对比,主要是查看相应的粗差信息(见表4-1)及粗差探测过程、粗差的大小（定值）和所在的观

43

测值（定位），参见表4-1。程序将对粗差作自动改正（即粗差剔除），重新生成剔除粗差后的观测值文件——“网名C.IN2”。再对文件“网名C.IN2”平差，相应的结果将存放在“网名C.OU2”文件中，为粗差剔除后的网平差结果，其后验单位权中误差将显著减小。

表 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

44

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

4.1.2．附合导线的粗差探测与剔除

本系统专门提供了附合导线的粗差探测与剔除功能，对于附合导线，当且仅当只存在一个方向含粗差或只存在一条边含粗差的情况，则可应采用本模块，根据导线平面观测值文件（*.IN2），自动确定出粗差的位置和大小。在附合导线平差和闭合差计算中，若发现计算结果不满足精度要求，则导线观测值中可能含有粗差。在进行附合导线粗差探测时，导线的方位角闭合差限值取为10?n，n为导线观测测站数，边长闭合差限值取为5厘米。

在进行“附合导线粗差探测”时，对导线平面观测值文件（*.IN2）

45

的要求是，文件中的第一个测站应是导线的起点，后面各站的顺序不作要求。文件中的已知点数为附合导线的四个已知点，点的先后顺序不作要求。

在 “平差”菜单栏中点击“粗差探测”下的“附合导线”菜单项，根据提示打开所选的“网名.IN2”文件，则开始进行粗差计算。

若导线的方位角闭合差和边长闭合差小于给定的限值，则认为没有粗差，报告如图4-1和图4-2所示。

图4-1 没有粗差报告

图4-1

46

图4-2

如果方向观测值含有粗差，报告如图4-3所示。

图4-3 发现有方向粗差的报告

探测结果表明导线中测站ZF3的方向观测值含有粗差，该站的角度观测值大了约122秒。应对ZF3站进行重测。

如果边长观测值含有粗差，报告如图4-4所示。

47

本文来源：https://www.bwwdw.com/article/qo9o.html