GPS工程测量网通用平差软件包CosaGPS
科傻系统(COSA)系列软件
GPS工程测量网
通用平差软件包
(CosaGPS V5.1)
使用说明书
2007年11月
CosaGPS
GPS工程测量网通用平差软件包CosaGPS
版权所有 不得翻录
Tel: 13100727893 13628666950
13317123688 13317198509 Email: jmguo@sgg.whu.edu.cn
nxluo@sgg.whu.edu.cn
CosaGPS
GPS工程测量网通用平差软件包CosaGPS
目录
目录 ......................................................................................................1 1.简介 .................................................................................................3
1.1 功能全面 ...............................................................................3 1.2 整体性好 ...............................................................................3 1.3 解算容量大,运算速度快 ...................................................3 1.4 操作简明,使用方便 ...........................................................4 2.“文件”下拉菜单 .............................................................................6
2.1 工程与文件 ...........................................................................6 2.2 “文件”菜单项 ...................................................................8
2.2.1新建 .............................................................................8 2.2.2打开 .............................................................................9 2.2.3关闭 .............................................................................9 2.2.4保存 .............................................................................9 2.2.5另存为 .........................................................................9 2.2.6新建工程 .....................................................................9 2.2.7 打开工程 ..................................................................15 2.2.8 打印 ..........................................................................16 2.2.9 打印预览 ................................................................16 2.2.10 打印设置 ................................................................16 2.2.11 退出 ........................................................................16
3.“GPS数据处理”下拉菜单 .............................................................17
3.1 已知数据 .............................................................................17
3.1.1 三维已知坐标 ..........................................................18 3.1.2 二维已知坐标 ..........................................................19 3.1.3 一维高程点 ..............................................................19 3.1.4 输入地面边长 ..........................................................19 3.1.5 输入地面方位 ..........................................................20 3.2 基线数据 .............................................................................20 3.3 GPS三维向量网平差(无约束平差或约束平差) ..........21 3.4 二维网联合/约束平差 ........................................................22
CosaGPS
1
GPS工程测量网通用平差软件包CosaGPS
3.4.1 联合/约束平差 .........................................................22 3.4.2 输出用户自定义任意两点相对精度 .......................23 3.5 椭球面上三维平差 .............................................................23 3.6 工程网(一点一方向)平差 ...................................................24 3.7 GPS高程拟合 ......................................................................25 3.8 GPS三维秩亏自由网平差 ..................................................26 3.9 稳定性分析 .........................................................................27 3.10 设置 ...................................................................................28 4.“查看”下拉菜单 ...........................................................................28 5.“工具”下拉菜单 ...........................................................................29
5.1 闭合差计算 .........................................................................30 5.2 重复基线差 .........................................................................30 5.3 网图显绘 .............................................................................31 5.4 贯通误差影响值计算 .........................................................31 5.5 GPS网设计 ..........................................................................32 5.6 输出AutoCAD格式的GPS网图 ......................................33 6.“坐标转换”下拉菜单....................................................................34
6.1 XYZ-〉BLH .........................................................................34 6.2 BLH->XYZ ..........................................................................35 6.3 BL->XY ................................................................................36 6.4 XY->BL ................................................................................36 6.5 XY1->XY2 ...........................................................................37 6.6 XYZ1->XYZ2 ......................................................................39 6.7 高程面坐标变换 .................................................................41 7.“帮助”下拉菜单 .........................................................................42 附录1. 功能菜单框图 ......................................................................43 附录2. 算例及说明 ..........................................................................44 附录3. 基线解文件格式说明 ..........................................................46 附录4. 方向及经纬度的角度格式说明 ..........................................54 附录5. 简要操作步骤 ......................................................................55
CosaGPS
2
GPS工程测量网通用平差软件包CosaGPS
1.简介
基于全球卫星定位系统(GPS)的现代测量理论和技术改变了传统的测量模式,使工程测量行业发生了革命性变化,测量外业工作自动化程度大大提高,测量内业软件的作用更加重要。为了满足工程测量单位对GPS数据处理的要求,在分析研究GPS数据处理理论的基础上,我们研制了自主版权的CosaGPS软件系统,该软件具有如下特点:
1.1 功能全面
软件具有在世界空间直角坐标系(WGS-84)进行三维向量网平差(无约束平差和约束平差)、在椭球面上进行卫星网与地面网三维平差、在高斯平面坐标系进行二维联合平差、针对工程独立网的固定一点一方向的平差、高程拟合等功能,并带有常用的工程测量计算工具,可以实现各种坐标转换。
1.2 整体性好
全部软件集成在统一的环境下,编辑器、文档、图形、数据处理模块均自主编写;采用多文档,可同时处理多项任务;采用工程管理模式,可方便进行各类数据的操作。
1.3 解算容量大,运算速度快
软件设计采用节省内存的快速算法,在现有的大部分微机上(Winows98,Windows2000,WindowsXP),可整体解算数千个控制点的GPS控制网,内存不够时则采用外存作缓冲,因而还可解算更
CosaGPS
3
GPS工程测量网通用平差软件包CosaGPS
大规模的GPS工程网。
1.4 操作简明,使用方便
在WIN95/98/2000/XP系统环境下运行,可采用表格方式或文本方式进行数据录入,大部分操作采用“傻瓜“式选项。对于输入量较少的已知数据和参数,采用表格方式输入;对于大批量的数据,则采用文件方式输入。表格方式输入时,屏幕上显示格式如图1.1。
图1.1 输入数据表格
表格中各列的宽度可以改变,将鼠标移到表格中各列标题的结合处,按下左键拖动,调节到合适的宽度即可。表格中的行数是不受限制的,输满后将向上滚动,底下弹出新的空白行。表格式输入的数据被保存到文本文件(文件名参见2.1),用户也可直接对相应的文本文件(例如:工程名. GPS3dKnownXYZ)进行修改,重新进入表格后,表格中的数据将自动进行更新。
表格中各列的宽度可以改变,将鼠标移到表格中各列标题的结合处,按下左键拖动,调节到合适的宽度即可。表格中的行数是不受限制的,输满后将向上滚动,底下弹出新的空白行。表格
CosaGPS
4
GPS工程测量网通用平差软件包CosaGPS
式输入的数据被保存到文本文件(文件名参见2.1),用户也可直接对相应的文本文件(例如:工程名. GPS3dKnownXYZ)进行修改,重新进入表格后,表格中的数据将自动进行更新。
文件方式输入数据时,可以使用本系统的编辑器,也可以使用其它的文本编辑器(比如记事本、书写器、Word等,应选择纯文本方式),其操作基本相同,可以使用快捷键(复制:Ctrl+C,粘贴:Ctrl+V), 或者点击鼠标右键后,根据弹出的快捷菜单进行操作,也可以利用屏幕顶部的“编辑”菜单进行文本处理。对文件编辑完毕后,用屏幕顶部的“文件”-〉“保存”菜单进行同名保存(与鼠标点击相同)或用“另存为”保存为另一文件名。 系统主菜单参见图1.2。
图1.2系统主菜单
CosaGPS
5
GPS工程测量网通用平差软件包CosaGPS
2.“文件”下拉菜单
2.1 工程与文件
工程是指某项确定的任务,它是所有与之相关文档的集合,其中单个的文档称为文件,相关的文件通过工程而联系在一起。该系统是按工程进行管理和处理的,大部分操作是对所选定的工程进行的,这样做的优点是方便用户使用,便于存档和调阅,同时,按工程进行管理也是大部分优秀软件所采用的方法,为软件
图2.1 文件菜单
CosaGPS
6
GPS工程测量网通用平差软件包CosaGPS
使用人员所熟悉。观测数据文件和平差结果文件等都是与工程有关的文档,一个工程会涉及到许多的文档,根据一定的命名规则,系统会调用相应的文档进行处理。“文件”下拉菜单参见图2.1。
工程名一般采用地区或测区名称,这样易于记忆,其构成形式为*.prj。*是用户自己定义的,可由汉字、英文字母、数字、符号等组成,后缀prj是系统指定的,
系统把以prj为后缀的文件看作是工程文件。另外,还有许多的数据文件和结果文件,其命名规则及含义为: ? 与工程有关的GPS文件: *
工程名.GPS1dKnownH 已知高程文件 *
工程名.GPS2dKnownXY 已知平面坐标文件 *
工程名.GPS3dKnownXYZ 已知三维坐标文件 *
工程名.GPS2dAzimuth 地面方位角 *
工程名.GPS2dDistance 地面边长
工程名.GPS3dVector GPS三维基线向量 工程名.GPS2dVector GPS二维坐标差向量 工程名.GPS3dBLHVector GPS三维大地坐标差向量 工程名.GPS1dResult GPS高程拟合结果 工程名.GPS2dResult GPS二维联合平差结果 工程名.GPS3dResult GPS三维向量网平差结果
工程名.GPS3dBLHResult GPS三维网椭球面上联合平差结果 工程名.GPS3dBLH GPS三维大地坐标文件
工程名.GPS3dXYH GPS平面坐标和大地高文件
工程名.GPS3dXYHEFT GPS平面坐标、大地高、误差椭圆元素文件
工程名.GPS2dXYEFT GPS二维联合平差高斯平面坐标及误差椭圆元素文件
工程名.dxf AutoCAD的DXF格式的网图文件 *
表格方式输入的数据文件,也可以用文本编辑器进行编辑 固定一点一方向的工程网有关文件
工程名.OneFix 输入的已知数据文件,与对话框对应
CosaGPS
7
GPS工程测量网通用平差软件包CosaGPS
工程名.GPS2dResult1 GPS二维平差结果
工程名.GPS3dResult1 GPS三维向量网平差结果 闭合差计算文件 工程名.GPS3dLoop
工程名.GPS3dMisclosure
贯通误差影响值计算输入输出文件(参见5.4) 工程名.gti 输入文件 工程名.gto 输出文件 ? 转换参数文件:
Parameter.1d 高程拟合模型系数
Parameter.2d 二维转换旋转角及尺度因子 ? 坐标转换算例文件:
demo.xy 高斯平面直角坐标 demo.BL 大地经纬度
demo.XYZ 三维空间直角坐标 demo.BLH 三维大地坐标
demo.XYXY 不同平面坐标系坐标转换
demo.XYXY_O 不同平面坐标系坐标转换结果 demo.XYZXYZ 不同空间直角坐标系坐标转换
demo.XYZXYZ_O 不同空间直角坐标系坐标转换结果 ? 用户自定义文件 demo.GPS2dRel 用户自定义需要输出相对精度的点对文件
2.2 “文件”菜单项
在主菜单下用鼠标单点“文件”,弹出如图2.2所示的下拉菜单,其中各项含义是:
2.2.1新建
用该系统的编辑器建立新文本文件
CosaGPS
8
GPS工程测量网通用平差软件包CosaGPS
2.2.2打开
用该系统的编辑器打开已有的文本文件
2.2.3关闭
关闭当前活动窗口
2.2.4保存
保存当前活动窗口的文件
2.2.5另存为
换名保存当前活动窗口的文件
2.2.6新建工程
数据处理是按工程进行的,必须首先建立工程。选择此项,弹出如图2.2所示窗口。在该窗口中输入有关的工程参数,其中有:工程、控制网、接收机/基线解类型、投影面大地高、坐标加常数五个组框和中央子午线、测区平均纬度两个编辑框。
2.2.6.1 工程组框
在工程组框中,输入工程名,工程所在路径二项,工程名是工程的标识,路径是工程所在的文件夹或目录。对于工程所在路径也可点取按钮 进行浏览选择,此时会出现如图2.3的浏览文件夹窗口,在此窗口中选择所需文件夹。在“新建工程”时,可立即进行参数设置,系统将记忆有关选项,以后可在“GPS数据处理”->“设置”项中查看和修改。
CosaGPS
9
GPS工程测量网通用平差软件包CosaGPS
图2.3 选择工程路径
CosaGPS
10
GPS工程测量网通用平差软件包CosaGPS
2.2.6.2 控制网组框
在控制网组框中,选定或者新增坐标系统、设定控制网等级。坐标系统是点位坐标的参考系,我国常用的测量坐标系统有: BJ54(北京54坐标), GDZ80(国家80坐标),WGS84坐标,城市坐标系,工程坐标系,独立坐标系。国家坐标系统参照于某个参考椭球,在同一参考椭球下,又有空间直角坐标、大地坐标、平面直角坐标。进行坐标转换需选择相应的椭球参数,椭球的几何参数可由长半轴和扁率分母确定。点压按钮出现如图2.4窗口。在该窗口中输入坐标系统的椭球长半轴和椭球扁率分母,然后可在右边对应的下拉框中选定所需的坐标系统,输入无误后按“确认”按钮。
图2.4 定义坐标系统
其中,“国家80坐标”、“WGS-84坐标”、“北京54坐标”是固定的,不能改变,“工程系1”是用户自定义的。在控制网组框右下角的下拉框中选择要求的坐标系,如图2.5所示。
CosaGPS
11
GPS工程测量网通用平差软件包CosaGPS
图2.5 选择坐标系
控制网等级是按下述系列划分的:
国标A级 国标B级 国标C级 国标D级 国标E级 城市二等 城市三等 城市四等 城市一级 城市二级 铁道B级 铁道C级 铁道D级 铁道E级
公路一级(线路) 公路一级(特殊) 公路二级(线路) 公路二级(特殊) 公路三级(线路) 公路三级(特殊)
CosaGPS
12
GPS工程测量网通用平差软件包CosaGPS
公路四级(线路)
自定义(仪器固定误差,比例误差)
在全球定位系统(GPS)测量规范(GB/T 18314-2001)、全球定位系统城市测量技术规程(CJJ 73-97)、无碴轨道铁路工程测量技术暂行规定、全球定位系统公路测量规范中规定了各个等级GPS控制网的固定误差和比例误差,按照控制网的实际等级进行选择。对于特殊网,如果没有包含在上述等级中,则可选择自定义,即采用仪器框中输入的仪器固定误差和比例误差,用户可以进行编辑和修改。
2.2.6.3 接收机/基线解类型组框
各个GPS接收机生产厂家提供了相应的基线解算软件,比如
Trimble的TGO、Leica的SKI、Topcon(Javad)的Pinnacle、Ashtech的Solution等,不同基线解算软件求得的基线向量的输出格式是不同的,CosaGPS支持的软件格式有:
Trimble(GPSurvey/TGO) Ashtech(GPPS/Solution) Leica(SKI/LGO) Sokkia Rouge Lip
CosaGPS
Topcon/Javad(Pinnacle) Gamit Novatel
Zhonghaida(中海达)
当采用了两种以上软件解算得到网中的基线向量时,首先查看不同软件的基线向量的方差之比是否存在系统性偏差,若其比值为1:m1:m2,则进行匹配处理,对第1种软件的基线输入1作为基线方差因子,生成CosaGPS的基线输入文件(工程名.GPS3dVector),将其名称改为V1.GPS3dVector,对第2种软件的基线输入m1作为基线方差因子,生成CosaGPS的基线输入文
CosaGPS
13
GPS工程测量网通用平差软件包CosaGPS
件(工程名.GPS3dVector),将其名称改为V2.GPS3dVector,对第3种软件的基线输入m2作为基线方差因子,生成CosaGPS的基线输入文件(工程名.GPS3dVector),将其名称改为V3.GPS3dVector,最后将V1、V2、V3三个文件合并在一起并命名为工程名.GPS3dVector,再进行后续平差处理。
接收机框中的固定误差(mm)、比例误差(ppm)、改造基线方差阵是根据GPS接收机的精度指标对基线的方差阵进行修正。一般情况下,不应在检查框中打勾(即不启用修正功能);只有当验后单位权中误差很大时(说明基线向量的方差阵不准确),将该项选中,软件将只利用基线解方差阵的相关性,同时利用仪器的标称精度(接收机的固定误差、比例误差)重新构造方差阵进行网平差。
用验前单位权中误差检查框决定平差结果的精度指标是基于验前值还是验后值,当网中多余观测量较少时,例如当闭合环的个数少于4时,验后单位权中误差是不够准确的,可以采用验前单位权中误差(1cm)。
独立基线条数:省缺值为-1,即认为选定的基线全部为独立基线;若选择了全部基线进行平差(含有同步基线),则平差后的精度指标比实际值偏高,但坐标、边长、方位角仅有微小变化,在此输入独立基线的实际条数,软件将对平差后的精度指标进行修正,从而与独立基线平差结果的精度指标基本一致。 2.2.6.4 坐标加常数组框
坐标加常数是指坐标系常数,例如我国60带高斯坐标在y坐标上加500公里的常数,目的是为了避免出现负值。某些城市坐标系是以过城市中心或某特定点的子午线为中央子午线,往往在高斯坐标上加减一个平移常数。此处的坐标加常数起类似的作用,对GPS三维向量网平差结果中转化的高斯平面坐标起作用,对6.3(BL->XY)和6.4(XY->BL)起作用,对二维联合平差不起作用。
在该组框中输入平面坐标的加常数,以公里为单位。 2.2.6.5 中央子午线、投影类型
在该编辑框中输入中央子午线的经度,格式为:DDD.MMSS,分和秒必须占满两位,该软件所有的角度值(方位角、纬度、经
CosaGPS
14
GPS工程测量网通用平差软件包CosaGPS
度)的输入均采用此格式。例如:114.300751表示114度30分7.51秒,详细说明参见附录4;目前该软件提供的投影类型为高斯投影、UTM投影两类,根据测量项目的需要进行选择,我国测量工程一般采用高斯投影。
2.2.6.6平均纬度、投影面大地高
这两项参数用于“坐标转换”/“高程面坐标变换”,对网平差的其它项目不起作用。计算控制点以不同高程面为参照面的坐标时,在中输入与当前坐标对应的参照面的大地高(正常高+高程异常),在中输入要转换到的参照面对应的
大地高。平均纬度可采用近似值,也可从地图上查取。
2.2.7 打开工程
对于已建立的工程,应选择“打开工程”项,此时弹出如图2.6的选择工程窗口:
图2.6 打开工程
在该窗口中直接输入或选定工程名(以PRJ为后缀的文件)后,用鼠标点
按钮,工程名将显示在主菜单的顶部标题栏中,
CosaGPS
15
GPS工程测量网通用平差软件包CosaGPS
以后的操作都是面向该工程(坐标转换工具除外)。
2.2.8 打印
打印活动窗口的文件
2.2.9 打印预览
标准Windows打印预览窗口
2.2.10 打印设置
打印格式及打印机设置
2.2.11 退出
退出系统
CosaGPS
16
GPS工程测量网通用平差软件包CosaGPS
3.“GPS数据处理”下拉菜单
菜单形式参见图3.1。
图3.1 GPS数据处理
3.1 已知数据
已知数据又分为:三维已知坐标、二维已知坐标、一维高程点、地面边长、地面方位。该项数据有两种用途,一是用于控制网平差处理的解算基准,二是用于解求平面转换参数和高程拟合系数。
CosaGPS
17
GPS工程测量网通用平差软件包CosaGPS
3.1.1 三维已知坐标
作用是为三维平差输入固定点坐标。用鼠标单点该项,弹出如图3.1的窗口,必须至少输入一个点的三维坐标,可以是三维空间直角坐标(X,Y,Z),也可以是大地坐标(纬度B,经度L,大地高H), B、L的格式为:DDD.MMSS,X、Y、Z、H的单位是米。不能将(X,Y,Z)与(B, L, H)混合输入,并注意不要将三维空间直角坐标(X,Y,Z)中的(X,Y)与平面坐标(x,y)弄混。删去点名,该点即被删除,双击格网中的数据单元,其底色变白后,可修改数据,当输至底行时,会自动弹出新的空白行,所有数据向上翻动一行,列宽可用鼠标拖动来变宽或变窄。应特别注意点名必须与基线向量中的点名(起点、终点)完全一致。
图3.2 输入三维已知坐标
CosaGPS
18
GPS工程测量网通用平差软件包CosaGPS
3.1.2 二维已知坐标
操作与3.1.1相似,起作用是为二维联合平差输入地面公共点坐标,一般至少需要两个公共点,若仅有一个公共点,则应采用“固定一点一方位”的平差模式。应特别注意点名必须与基线向量中的点名(起点、终点)完全一致。
3.1.3 一维高程点
操作与3.1.1相似,作用是为高程拟合输入地面公共点的正常高。常数拟合模型至少需要一个公共高程点,平面拟合模型至少需要三个公共高程点,曲面拟合模型至少需要六个公共高程点。应特别注意点名必须与基线向量中的点名(起点、终点)完全一致。
3.1.4 输入地面边长
作用是为二维联合平差输入地面边长。鼠标单点该项,弹出图3.3的窗口,在该窗口中输入地面边长的起点、终点、边长值(m)、中误差(cm),对于已知边长,中误差输入0来表示。应特别注意点名必须与基线向量中的点名(起点、终点)完全一致。
图3.3 输入地面边长
CosaGPS
图3.3
19
GPS工程测量网通用平差软件包CosaGPS
3.1.5 输入地面方位
作用是为二维联合平差输入地面方位角,操作与3.1.4相似。在该窗口中输入地面方位角的起点、终点、方位角值(DDD.MMSS)、中误差(秒),对于固定方位角,中误差输入0来表示。应特别注意点名必须与基线向量中的点名(起点、终点)完全一致。
3.2 基线数据
作用是选择所需的基线解文件,用鼠标单点该项,弹出如图3.4的窗口。在该窗口中选取所需基线解的文件,首先点压“浏览”按钮,弹出选择文件夹对话框,选择需要的文件夹。然后在“类型”下拉组合框指定文件类型,再到“待选基线文件”窗单点鼠标,Shift键+鼠标左键顺序选择文件,Ctrl键+鼠标左键任意选择文件。找到所需的基线文件并标记好之后,单击“选定-〉”按钮,所需的基线文件将出现在已选基线文件窗中,点压“确定”按钮后,形成的基线向量文件(*.GPS3dVector)显示在屏幕上,这也就是三维向量网平差所需要的基线向量输入文件。应特别注意设置对话框中的“接收机/基线解类型”与选择的基线向量文件相对应。
CosaGPS
20
GPS工程测量网通用平差软件包CosaGPS
图3.4 选择基线解算文件
3.3 GPS三维向量网平差(无约束平差或约束平差)
作用是在WGS84空间直角坐标系中进行三维向量网平差,首先需要至少输入一个点的三维坐标(参见3.1.1中的三维已知坐标
CosaGPS
21
GPS工程测量网通用平差软件包CosaGPS
项)并生成基线向量文件(*.GPS3dVector,参见2)。对于独立的GPS网,可取一个点的单点定位解(从基线解文件查取)作为固定坐标,进行无约束平差;若网中联测了多个国家GPS点(比如A级点、B级点),可全部作为固定点输入,进行约束平差。可以用(X,Y,Z)或(B,L,H)的格式输入。
以表格形式显示在窗口中的坐标数据与名称为*.GPS3dKnownXYZ的文件内容互相对应,也可用文本编辑器编辑生成*.GPS3dKnownXYZ文件,表格中的数据随之改变,格式为:
$$$$$$$$ ********.****** ********.****** ********.****** (点名) (X / B) (Y / L) (Z / H) 对于同一控制网,如果采用不同生产厂商的多种类型接收机观测,并用各自的配套软件进行解算,得到的基线向量,有时存在方差阵不匹配问题,可采用不同的方差因子(参见2.2.6.3)对基线的方差阵进行处理。操作方式是先用不同的工程名对每类基线给定相应的方差因子(参见2.2.6.3)分别生成相应的基线向量文件(*.GPS3dVector),然后用编辑器将其合并为一个作为最终的输入文件。
若想改变椭球参数,可到“设置”项选择。
完成坐标输入并生成基线向量文件后,单点该菜单项进行平差计算,结果文件(*.GPS3dResult)将显示在屏幕上。
3.4 二维网联合/约束平差
3.4.1 联合/约束平差
作用是进行二维联合平差,首先需要完成三维向量网平差并至少输入一个公共点的二维平面坐标(参见3.1.2的二维已知坐标项),若只有一个公共点,则还需要输入至少一条地面边长(归算到高斯平面上)和一个地面方位角,当然也可以输入任意多个地面边长和方位角。
CosaGPS
22
GPS工程测量网通用平差软件包CosaGPS
地面边长和方位角可作为观测值进行联合平差,也可作为固定值进行约束平差(参见3.1.4和3.1.5)。
以表格形式显示在窗口中的坐标数据与名称为*.GPS2dKnownXY的文件内容互相对应,也可用文本编辑器编辑生成*.GPS2dKnownXY文件,表格中的数据随之改变,格式为:
$$$$$$$$ ********.****** ********.****** (点名) (x) (y)
若想改变椭球参数,可到“设置”项选择。
完成坐标输入并生成基线向量文件后,点击该菜单项进行平差计算,结果文件(*.GPS2dResult)将显示在屏幕上。
3.4.2 输出用户自定义任意两点相对精度
CosaGPS提供了用户自定义输出任意两点间相对精度的功能。具体方法为,首选需要形成用户要求的“点对”文件,其文件名为:“工程名.GPS2dRel”,其格式为一文本文件,每一行即为一个“点对”(起点点名,终点点名),“点对”间用逗号或空格分隔,如:
A01, A02 A03, A08
系统在平差时,自动判断该文件是否存在,若存在,则读取文件中的点对,并计算其相对精度,输出到二维平差结果文件中的“平差后方位角、边长及精度”信息栏中,为和存在直接观测值的“点对”相对精度有所区别,其序号为“****”。
3.5 椭球面上三维平差
在某一确定的椭球面上进行三维平差,把WGS84椭球到地方
CosaGPS
23
GPS工程测量网通用平差软件包CosaGPS
参考椭球的转换参数作为附加参数,在平差时一并求得。解算是在椭球面上进行的,不受投影变形的限制,可以进行覆盖全国乃至全球的大范围GPS控制网的数据处理。
首先需要完成三维向量网平差并至少输入三个公共点的三维坐标(参见3.1.1的三维已知坐标项),可以输入(X,Y,Z)或(B,L,H)或(B,L)或H。即只知道某点的经纬度时,输入(B,L),H为空(不要输0);只知道某点的大地高时,输入H,(B,L)为空(不要输0)。
结果文件名为*.GPS3dBLHResult. 其它与3.3相似。
3.6 工程网(一点一方向)平差
对于某些工程项目,比如大桥、大坝等,如果采用固定一个点的坐标、指定一个方向的方位角,并且选择相应的工程投影面,从而建立相对独立的坐标系,则可选用该平差项。点击“GPS数据处理/工程网(一点一方向)平差”,屏幕显示图3.5的对话框。
图3.5工程网(一点一方向)平差的输入信息
CosaGPS
24
GPS工程测量网通用平差软件包CosaGPS
在图3.5的对话框中,固定点信息有点名、平面坐标和正常高、大地坐标和投影面正常高,平面坐标可以是工程坐标系的独立坐标或高斯平面坐标;固定方位角一般是工程网中某一特定方向的方位角,例如大桥控制网的桥轴线方向、大坝控制网的坝轴线方向等。在平差前,应在“数据处理/设置”对话框中选择相应的椭球参数和中央子午线,一般是选择工程网对应的地方椭球的参数。相应的数据文件为: “工程名.OneFix”(输入的已知数据文件,与对话框对应),“工程名.GPS2dResult1”(GPS二维平差结果),“工程名.GPS3dResult1”(GPS三维向量网平差结果)
3.7 GPS高程拟合
作用是进行高程拟合,首先需要完成三维向量网平差并至少输入一个公共点的高程(参见1的已知高程点项),点取该菜单项后,弹出如图3.5的窗口,在该窗口中选择拟合模型,其中常数模型需要一个以上的公共高程点,平面模型需要三个以上的公共高程点,曲面模型需要六个以上的公共高程点。高程结果文件名为*.GPS1dResult.
图3.6 选择高程拟合模型
CosaGPS
25
GPS工程测量网通用平差软件包CosaGPS
3.8 GPS三维秩亏自由网平差
采用秩亏自由网平差时,选择该菜单项。首先进行“三维向量网无约束平差”(参见3.3),得到近似坐标文件“工程名.GPS3dApproximateXYZ”,该文件是进行秩亏自由网平差平差的输入文件;然后进行“GPS三维秩亏自由网平差”,由程序自动生成的用于秩亏自由网平差的输入文件有:
工程名.GPS3dApproximateXYZ 工程名.GPS3dFreeXXInput 工程名.GPS3dFreeGxxInput 工程名.GPS3dFreeQxxInput 其中,“工程名.GPS3dFreeXXInput”的内容为:
总点数 多余观测数 单位权中误差(mm)
点名 1.0 dX(mm) dY(mm) dZ(mm) X(m) Y(m) Z(m) ……
点名之后的“1.0”是一个标识数字,含义为该点属于稳定点组,可以根据实际稳定情况进行手工修改为“0.0”,则认为该点是一个不稳定点,不属于稳定点组。手工修改该文件之后,在屏幕提示“工程名.GPS3dFreeXXInput已存在,重新产生?”,应选择“否”。
用于稳定性分析的输入文件有:
工程名.GPS3dFreeXXQxx 工程名.GPS3dFreeENUQenu
文件“工程名.GPS3dFreeXXQxx”的内容为:
总点数 多余观测数 单位权中误差(mm) 点名 dX(mm) dY(mm) dZ(mm) …… 协因数阵
文件“工程名.GPS3dFreeENUQenu”的内容为:
总点数 多余观测数 单位权中误差(mm) 点名 dE(mm) dN(mm) dU(mm)
CosaGPS
26
GPS工程测量网通用平差软件包CosaGPS
…… 协因数阵
秩亏自由网平差的结果文件为:
工程名.GPS3dFreeXXOut 工程名.GPS3dFreeQxxOut 工程名.GPS3dFreeMxxOut
工程名.GPS3dFreeMxxOut的内容为:前半部分为经典自由网平差(网中含有一个固定点)的空间直角坐标及其中误差,后半部分为秩亏自由网平差(网中含有一个固定点)的空间直角坐标及其中误差。
3.9 稳定性分析
对于两期观测的变形监测网,各自先进行三维向量网无约束平差,再进行三维秩亏自由网平差(参见3.8),最后再进行稳定性分析。
特别需要注意的是两期控制网的控制点应相同,进行三维向量网无约束平差时的未知数序号应相同。为了满足未知数序号相同的要求,两期网中的固定点应是同一个点,并且基线向量文件中前面的基线向量的顺序在两期中应相互对应,保证推算近似坐标的顺序相同,从而保证了未知数的序号相一致。进行两期自由网平差时,应采用相同的近似坐标。
输入文件为:
第一期的.GPS3dFreeXXQxx或者.GPS3dFreeENUQenu 第二期的.GPS3dFreeXXQxx或者.GPS3dFreeENUQenu
结果文件为:
工程名.GPS3dFreeQdd 工程名.GPS3dFreeD
CosaGPS
27
GPS工程测量网通用平差软件包CosaGPS
文件“工程名.GPS3dFreeD”中的结果与输入文件有关,当输入文件是“.GPS3dFreeXXQxx”时,其中各行的数据为:
点名 dX(mm) dY(mm) dZ(mm) 自由度 标准化统计量XYZ ……
当输入文件是“.GPS3dFreeENUQenu”时,其中各行的数据为:
点名 dE(mm) dN(mm) dU(mm) 自由度 标准化统计量ENU ……
3.10 设置
作用是改变所需的参数设置,操作与第一部分的新建工程相同,参见图2.2。
当工程已建好之后,若想再改变其中的参数,就必须选择该项。
注意:应经常查看设置项,确信有关参数正确。
4.“查看”下拉菜单
作用是查看主要的数据文件内容(平差结果、观测值)、设置屏幕显示项(工具栏、状态栏),参见图4.1。
图4.1 查看数据文件
CosaGPS
28
GPS工程测量网通用平差软件包CosaGPS
5.“工具”下拉菜单
该工具菜单参见图5.1,包含:闭合差计算,重复基线差,网图显绘,贯通误差影响值计算
图5.1 工具菜单
图5.2 闭合环
CosaGPS
29
GPS工程测量网通用平差软件包CosaGPS
5.1 闭合差计算
闭合差计算是自动搜索和计算所有最小环路的闭合差。在执行3.2所述“选取基线”生成“工程名.GPS3dVector”后,利用该菜单项,将生成2个文件:“工程名.GPS3dLoop”和“工程名.GPS3dMisclosure”。“工程名.GPS3dLoop”是各个环路的端点点名,每个环路占一行,各点按连接顺序排列,如图5.2的环路可表示为:A01 A02 A03 A04。若需计算某一指定环路的闭合差,可以在“工程名.GPS3dLoop”中输入与该环路对应的点名构成的一行,保存文件后,再进行闭合差计算,屏幕弹出如图5.3所示的对话框,应选“否”(不重新生成闭合环路文件),计算完毕后,结果存到“工程名.GPS3dMisclosure”中。
图5.3 闭合环提示窗口
5.2 重复基线差
为了计算重复基线向量各分量的差值,可在选取基线时(参见3.2),选取全部基线(不是为了平差目的,只是为了计算闭合差、重复基线差),在执行3.2所述“选取基线”生成“工程名.GPS3dVector”后,利用该菜单项,将生成文件“工程名.GPS3dRepeatBaseline”,其内容是重复基线有关的数据。
CosaGPS
30
GPS工程测量网通用平差软件包CosaGPS
5.3 网图显绘
网图显绘[WGS84]显示三维平差时的网图。由于二维联合平差时,增加了地面控制点,已知点将与三维平差时[WGS84]有所不同,因而应采用“网图显绘[地方坐标]”进行二维联合平差的网图绘制。
5.4 贯通误差影响值计算
该功能主要为隧道施工控制网而设计的,其实质是根据控制网的洞口点和定向点精度、贯通点的位置以及贯通面的方向,在完成网平差之后,直接估算隧道贯通误差影响值。为此首先人工建立一个贯通误差引导文件,该文件也是一个标准的ASCⅡ文件,命名规则为“工程名.gti”,其格式为:
进口点号,进口定向点号,出口点号,出口定向点号,贯通点号,X坐标,Y坐标,贯通面方位角.
如图5.4所示的贯通方案(贯通点为MMM)可建立如下贯通文件 A12,A02,A03,A14,MMM,X,Y,α 算例数据(demo.gti)如下所示:
A12,A02,A03,A14,MMM,3320821.76615,2337.9802,85.3045 为了选取最优的定向点方案,在一次计算中,可准备多种不同的进出口点与定向点的组合,每一种组合占一行。准备好引导文件“工程名.gti”后,用鼠标单击“贯通误差影响值计算”,将自动计算贯通误差影响值,并将结果存放在文件“工程名.gto”中(参见“demo.gto”)。
CosaGPS
31
GPS工程测量网通用平差软件包CosaGPS
图5.4 贯通误差影响值计算
注意:在计算贯通误差影响值之前,除准备好引导文件外,还须对该控制网进行平差计算,否则贯通误差计算失败。未建立贯通误差引导文件时,将弹出提示建立贯通误差引导文件的对话框。
5.5 GPS网设计
选择“工具/GPS网设计”(参见图5.5)进行方案设计。在进行GPS网设计时,需要准备两个文件: 1)各测站的设计坐标文件(*.xyz) 格式为:
基线固定误差(mm) 基线比例误差(ppm) 基线水平方位误差(“) 点名 X(m) Y(m) Z(m) ……
各项之间用空格分隔,(X Y Z)为地心为原点的空间直角坐标。如
CosaGPS
32
GPS工程测量网通用平差软件包CosaGPS
果只有各测站的高斯平面直角坐标的设计值,则可以利用“坐标转换”功能得到经纬度(XY->BL),再把大地高加入形成BLH文件,进一步转换为XYZ文件(BLH->XYZ)。 2)GPS独立基线定义文件(*.line) 格式为: 起点 终点 ……
图5.5 GPS网设计
5.6 输出AutoCAD格式的GPS网图
选择“工具/输出ACAD_DXF网图”(参见图5.5)则输出AutoCAD的DXF格式的控制网网图。
输入文件为:
工程名.GPS3dXYHEFT 工程名.GPS3dVector
这两个文件是在“三维向量网平差”(参见3.3)时生成的。 结果文件为:*.dxf
CosaGPS
33
GPS工程测量网通用平差软件包CosaGPS
6.“坐标转换”下拉菜单
作用是进行各种坐标变换计算,在这里提供了多项非常有用的工具,参见图6.1。
该菜单项中的坐标变换与网平差的工程名无关,是独立的测量计算辅助工具,执行完本项功能后,若再执行其它平差功能时,出现文件找不到或其他异常情况,是由于路径已被改变,执行“文件”-〉“打开工程”即可恢复原来的工程路经和其他设置好的参数。
图6.1 坐标转换
6.1 XYZ-〉BLH
功能:三维空间直角坐标转换为大地坐标。 输入文件 *.XYZ 输出文件 *.BLH
CosaGPS
34
GPS工程测量网通用平差软件包CosaGPS
上述两文件名称相同,后缀不同。操作步骤:
1) 在进行转换前选择“文件”菜单的“新建”项,编辑生成输入
文件 *.XYZ,其格式为:
点名 X Y Z
也可采用任何其它文本编辑器生产该文件。各项以空格分隔,X、Y、Z以米为单位。
2) 设置坐标系参数:椭球长轴,椭球扁率分母(参见3.7)。 3) 在“坐标转换”菜单下选“XYZ-〉BLH”子项进行转换,此时屏幕出现打开文件窗口,选取输入文件(*.XYZ)
4) 转换结果将保存到*.BLH文件中,若已存在同名点文件,则在其尾部追加数据,同时屏幕显示输入文件(*.XYZ)和输出文件(*.BLH)的窗口,供用户查看。 5) 结果文件格式为: 点名 B L H
经纬度B、L以度分秒(DDD.MMSS)为单位,大地高H以米为单位。
6) 若想查看其它文件,可选取“文件”菜单的“打开”子项,在屏幕出现“打开文件”窗口后点取打算查看的文件名,该文件的内容将出现在屏幕窗口上。
6.2 BLH->XYZ
功能:大地坐标转换为三维空间直角坐标。 输入文件 *. BLH 输出文件 *. XYZ
上述两文件名称相同,后缀不同。操作步骤:
1) 在进行转换前需编辑生成输入文件 *.BLH,其格式为:
点名 B L H
各项以空格分隔,经纬度B、L以度分秒为单位,大地高H以米为单位。
2) 设置坐标系参数:椭球长轴,椭球扁率分母(参见3.7)。
CosaGPS
35
GPS工程测量网通用平差软件包CosaGPS
3) 执行程序进行转换,若已存在同名点文件,则在其尾部追加数据。
4) 结果文件格式为:
点名 X Y Z
6.3 BL->XY
将大地经纬度转换成高斯平面坐标,与6.4联合使用可进行换
带计算。
其操作与6.4相似。
6.4 XY->BL
功能:高斯平面坐标转换为大地经纬度。 输入文件 *. XY 输出文件 *. BL
上述两文件名称相同,后缀不同。操作步骤:
1) 在进行转换前选择“文件”菜单的“新建”项,编辑生成输入文件 *.XY,其格式为:
点名 X Y
各项以空格分隔,X、Y以米为单位。若X、Y含有固定的加常数,则需在“参数设置”项输入X加常数和Y加常数(参见3.7)。 2) 设置坐标系参数:椭球长轴,椭球扁率分母,中央子午线(以度分秒为单位DDD.MMSS)(参见3.7)。
3) 在“坐标转换”菜单下选“XY-〉BL”子项进行转换,此时屏幕出现打开文件窗口,选取输入文件(*.XY)
4) 转换结果将保存到*.BL文件中,若已存在同名点文件,则在其尾部追加数据,同时屏幕显示输入文件(*.XY)和输出文件(*.BL)的窗口,供用户查看。
5) 结果文件*.BL格式为:
CosaGPS
36
GPS工程测量网通用平差软件包CosaGPS
点名 B L
各项以空格分隔,经纬度B、L以度分秒(DDD.MMSS)为单位。
6) 同6.1.6)。
6.5 XY1->XY2
功能:二维直角坐标变换,转换模型为:
?x??x0??cos? sin???x??????=?(1+k/1000000)??y??y0??-sin? cos?????y?? ??新??????旧x0,y0:平移参数 k:尺度参数(ppm) α:旋转角(弧度)
注意:CosaGPS输出的旋转角是以度分秒为单位(ddd.mmss) 输入文件 *. XYXY 输出文件 *. XYXY_0 坐标变换方法:
? 三参数法(平移、旋转)
转换模型为:
?x??x0??cos? sin???x????=?y??y0?????-sin? cos?????y?? ??新??????旧x0,y0:平移参数 α:旋转角(弧度)
注意:CosaGPS输出的旋转角是以度分秒为单位(ddd.mmss) ? 四参数法(赫尔默特法:平移、旋转、缩放)
转换模型为:
?x??x0??cos? sin???x??????=?(1+k/1000000)??y??y0??-sin? cos?????y?? ??新??????旧x0,y0:平移参数
CosaGPS
37
GPS工程测量网通用平差软件包CosaGPS
k:尺度参数(ppm) α:旋转角(弧度)
注意:CosaGPS输出的旋转角是以度分秒为单位(ddd.mmss) ? 六参数法(仿射变换法:平移、旋转、缩放)
转换模型为:
?ax??bx?x?????y??=??a???b??新?y??yax?x0cx??x????y?? cy????旧bx?(1?X尺度/1000000)?cos(X旋转)cx?(1?Y尺度/1000000)?sin(Y旋转)ay?y0cy?(1?Y尺度/1000000)?cos(Y旋转)
x0,y0:平移参数
X尺度:X尺度参数(ppm)
by??(1?X尺度/1000000)?sin(X旋转)Y尺度:Y尺度参数(ppm) X旋转:X旋转角(弧度) Y旋转:Y旋转角(弧度)
注意:CosaGPS输出的旋转角是以度分秒为单位(ddd.mmss)
操作步骤:
1) 编辑生成输入文件(*.XYXY),格式为(参见“demo.xyxy”):
点名 X旧 Y旧 X新 X新
公共点 ·······
·······
CosaGPS
38
GPS工程测量网通用平差软件包CosaGPS
点名 X旧 Y旧
非公共点 ········
········
下面为某一“*. XYXY”的实例:
A001 2886585.7645 531481.3787 29459.2790 31382.6500 A002 2891839.8474 479813.8408 34711.8360 -20285.9390 A003 2861398.5526 509387.9347 4271.0320 9289.5760 A004 2819325.6136 524901.0211 -37801.9780 24804.2570 B001 2860046.5165 524597.2659 2919.4310 24499.2250 P001 2865890.2564 519087.7625
2) 选择转换方法(3参数、4参数、6参数)。 3) 屏幕出现打开文件窗口后,选择输入文件。 4) 结果保存到文件*.XYXY_0中。
6.6 XYZ1->XYZ2
功能:三维直角坐标变换 转换模型为:
?X??X0?????Y=Y0(1+k/1000000)??????Z?????新?Z0?
?1?????z???y?z1??x??y??X?????x??Y??Z?1????旧X0,Y0,Z0:平移参数
k:尺度参数(ppm)
?x,?y,?z: X,Y,Z的旋转角(弧度)
注意:CosaGPS输出的旋转角是以度分秒为单位(ddd.mmss) 输入文件 *. XYZXYZ 输出文件 *. XYZXYZ_0
CosaGPS
39
GPS工程测量网通用平差软件包CosaGPS
操作步骤:
1) 编辑生成输入文件(*.XYZXYZ),格式为(参见
“demo.XYZXYZ”):
点名X旧 Y旧 Z旧 X新 Y新 Z新
公共点 ········
········ 点名 X旧 Y旧 Z旧
非公共点 ········
········
2) 选取“坐标转换”菜单的子项“XYZ->XYZ”,屏幕将出现“打
开文件”窗口,选择输入文件(*.XYZXYZ)。 3) 转换结果保存到文件*.XYZXYZ_0中。
CosaGPS
40
GPS工程测量网通用平差软件包CosaGPS
6.7 高程面坐标变换
对于工程控制网,例如桥梁、大坝等的施工控制网,为了减小投影变形,从而使用坐标反算的边长值与实测值之间的系统性差异较小,一般是以最重要的某一施工高程面作为坐标系的投影面,“高程面坐标变换”是将当前高程投影面对应的平面坐标变换到另一个高程面,须在工程的设置对话框中输入这两个高程面的大地高,并输入测区平均纬度,同时还应检查坐标系统对应的椭球参数是否正确。点击该项菜单后,屏幕提示输入数据文件,其格式为:点名 X Y。计算完毕之后,结果文件内容将显示在屏幕上,结构同输入文件,其名称是在输入文件名后加“1”。
图6.2 高程面坐标变换
CosaGPS
41
GPS工程测量网通用平差软件包CosaGPS
7.“帮助”下拉菜单
点击“关于CosaGPS2000”,屏幕显示程序标志,双击标志区或按Esc键返回程序状态, 参见图7.1。
图7.1 程序图标
CosaGPS
42
GPS工程测量网通用平差软件包CosaGPS
附录1. 功能菜单框图
文件(F) GPS数据处理(G) 查看(V) 工具(T) 坐标转换(C) 帮助(H) 新建(N) 打开(O) 关闭(C) 保存(S) 另存为(A) 新建工程( 打开工程 打印(P) 打印预览(V) 打印设置(R) 最近文档 退出(X) 已知数据 基线数据 三维向量网平差 二维联合平差 椭球面三维平差 工程网(一点一方向) GPS高程拟合 秩亏自由网平差 稳定性分析 设置 三维已知坐标 二维已知坐标 一维高程点 输入地面边长 输入地面方位 平差结果 观测值 闭合差计算 重复基线差 网图显绘 贯通误差 GPS网设计 输出CAD网图 三维基线向量 二维基线向量 XYZ->BLH BLH->XYZ BL->XY XY->BL XY1->XY2 XYZ->XYZ 二维转换参数 一维转换参数 帮助主题 关于CosaGPS 三维向量网 二维联合平差 一维高程拟合 三参数相似变换 四参数相似变换 六参数仿射变换 CosaGPS
43
GPS工程测量网通用平差软件包CosaGPS
附录2. 算例及说明
参见Example目录下的算例文件。 1) 工程文件 demo.prj
2) 已知坐标和观测值 demo.asc(基线解文件)
demo.GPS3dKnownXYZ(已知三维空间直角坐标) demo.GPS2dKnownXY(已知二维平面直角坐标) demo.GPS2dDistance(地面边长值) demo.GPS2dAzimuth(地面方位角值) demo.GPS1dKnownH(已知正常高)
demo.GPS3dVector(GPS三维向量观测值) 3) 平差结果
demo.GPS1dResult(高程拟合结果) demo.GPS2dResult(二维联合平差结果) demo.GPS3dResult(三维向量网平差结果)
demo.GPS3dBLHResult(椭球面上三维平差结果) demo.GPS3dBLH(大地经纬度、大地高)
demo.GPS3dXYH(高斯平面直角坐标和大地高) demo.GPS3dXYHEFT(三维向量网平差后转换得到的高斯平面直角坐标、大地高及误差椭圆元素)
demo.GPS2dXYEFT(二维联合平差的平面直角坐标和点位误差椭圆元素) 4) 坐标转换 demo.BL demo.XY
CosaGPS
44
GPS工程测量网通用平差软件包CosaGPS
demo.XYZ demo.BLH demo.XYXY demo.XYXY_O demo.XYZXYZ demo.XYZXYZ_O 5)闭合差计算文件 demo.GPS3dLoop
demo.GPS3dMisclosure
CosaGPS
45
GPS工程测量网通用平差软件包CosaGPS
附录3.基线解文件格式说明
1 CosaGPS基线文件格式
起点终点?X?Y2?Z??X?X2??Y?Y2??Z?Z??……
每条基线向量占一行,基线向量的各分量的单位为米,其方差/协方差的单位是平方厘米。 2 TGO基线文件格式
附图 1 TGO导出文件
CosaGPS
46
GPS工程测量网通用平差软件包CosaGPS
TGO软件可以根据用户指定的格式输出基线数据,输出CosaGPS相兼容的数据格式的操作如下(详情参见TGO使用说明书):
1). 选主菜单 File / Export(导出) 参见附图 1
2). 选 Custom / New format (自定义/新格式)创建新格式 参见附图 2
附图 2 TGO自定义新格式
CosaGPS
47
GPS工程测量网通用平差软件包CosaGPS
3) 按附图3所示定义格式文件 Name: CosaGPS
Export from: GPS Vector data (中文版:GPS基线向量)
中间大方格(Format body)的内容应在一行内输完,如下: [From Point Name:10] [To Point Name:10] [DeltaX:15.3] [DeltaY:15.3] [DeltaZ:15.3] [CovarianceXX:20.14] [CovarianceYY:20.14] [CovarianceZZ:20.14] [CovarianceXY:20.14] [CovarianceXZ:20.14] [CovarianceYZ:20.14]
附图3 文件格式
CosaGPS
48