globk通常的应用
1、重复性分析（glred）
-单独的观测期
-联合的观测期
2、联合多期数据得到基于实验的平均位置
-联合单独观测的测站
3、联合平均位置估计速度或地震偏移和震后运动
卡尔曼滤波
1、等同于序贯最小二乘估计，但是允许随机过程，通常是一阶马尔科夫过程
2、globk允许坐标，EOP，网络平移和缩放，卫星参数的随机游走过程；变量随时间线性增长
3、因为卡尔曼滤波用到协方差阵，而不是normal matrices，所以所有的参数必须有先验约束（通常是松弛的）
globk structural confusions
1、globk和glred是相同的程序，只是在对待h文件清单的方法有稍许不同：
-globk：所有的h文件联合成一个单独的解算
-glred：每个h文件产生一个单独的解算（除非文件清单里h文件后有+号）
2、两种类型的解算文件：
-自己保存的或从外部得到的h文件
-sol文件是内部的，格式因版本不同
3、globk/glred可以执行glorg，或者单独运行glorg以应用对解算的普遍约束
globk文件
用户要提供的：
-控制文件
-h文件清单：gdl文件
-二进制h文件
-apr文件
-EOP（in_pmu文件，可选的）
-eq_file（可选的，但是要出现在顶部）
globk产生的文件：
srt，com，sol，svs
输出文件：
屏幕，log，prt，输出h文件
globk文件操作
-log，prt，org文件是联系在一起的，所以应该删除或重命名，除非你都想得到（例如，glred）
-com，srt，sol文件是overwritten；com/sol不能重命名，因为原始sol文件名字嵌入在com文件中
-可以得到自动命名的com，srt，sol，org文件，和输出的h文件
估计命令的规则
1、globk中要估计的参数，必须使用apr_xxx，其中xxx代表参数类型（例如，neu，svs，wob，ut1，atm）
2、如果控制文件未涉及某个参数，它将不会在解算中出现，但是如果出现在h文件中，将会在globk解算文件中出现；例如：在gamit中估计了轨道，并且在globk中约束，使用apr_svs。
3、如果先验sigma是零值，将不会估计该参数
4、强制参数为先验值，先验sigma设为F
5、glorg估计的参数在globk中要保证松弛；如果在glorg中不估计旋转和缩放，在globk中要严密约束。
卫星轨道参数
通常apr_svs有两种使用格式：
-包含全球站点时，轨道可以估计的很好：apr_svs all 100 100 100 10 10 10 1R
-只使用当地数据，固定先验轨道：apr_svs all 0.1 0.1 0.1 0.01 0.01 0.01 0.01R
-每行末尾的R代表保留辐射参数
地球定向参数EOP
通常是两种格式：
-全球测站网（glorg中允许旋转）
apr_wob 10 10 1 1
apr_ut1 10 1
-区域网（约束的），这种方法的约束，系统不可以自由旋转
apr_wob .2 .2 .02 .02
apr_ut1 .2 .02
大多数分析中，为了在glorg中允许旋转估计，区域网也会使用全球的格式，但是如果网周围没有运动已经很好确定的测站时要小心使用）
数据编辑
为了说明时间序列的暂时相关性，我们对mar_neu命令使用随机游走（RW）过程噪声（单位m2/yr ）
-应用的典型值是2.5E-8 (一年0.5 mm ) 到4E-6 (一年2 mm )
mar_neu all 2.5E-8 2.5E-8 2.5E8 0 0 0
mar_neu chdu 4E-6 4E-6 4e-6 0 0 0
-为了降低噪声部分的权重或者使联合h文件中的持续的和测量模型的数据等同，可以添加随机噪声（单位是m）
sig_neu all .001 .001 .003
sig_neu ankr .005 .005 .020 2002 10 1 0 0 2002 11 30 24 0
sig_neu EMED0504 .010 .010 .1
-移除异常值，可以通过降低权重，或重命名
sig_neu ankr .1 .1 .1 2002 10 1 0 0 2002 10 1 24 0
rename ankr_gps ankr_xcl 2002 10 1 0 0 2002 10 1 24 0 ( eq_file )
glorg
-在h文件排列和执行松弛解算后通过globk执行以应用概括约束；也可以使用globk产生的com/sol文件单独运行程序。
-允许链接欧拉极参数和估计
-glorg中估计的参数在globk中要保证松弛（坐标，EOP，缩放）
内嵌到globk控制文件中glorg控制命令
org_cmd
org_opt <输出选项>
org_out <输出文件的名称>
在org_out没有给出打印文件的扩展名，将用org来代替
glorg命令
apr_file不需要跟globk中的一样，只需要包含用于固定的测站的值和坐标或速度相同的值
pos_org，rot_org：控制固定过程中的要估计的参数
xtran ytran ztran：允许平移
xrot yrot zrot：允许旋转
scale：允许缩放（如果要使用，在globk也要顾及缩放）
cnd_hgtv：控制高度的相对权重（变量，名义上时10）
stab_ite：移除测站的迭代和截止sigma
stab_site：固定过程中使用的测站清单
控制打印选项
crt_opt，prt_opt，org_opt：具体指定屏幕，打印和日志org文件的输出选项
globk/glorg帮助给出了所有的选项，主要的是：
-EARS：在写入时清除文件（通常是添加的文件）
-NOPR：不写输出（例如，执行glorg时的globk0
-BLEN：基线长度
-BRAT：估计速度时的极限比率
-RNRP：重命名后参数估计中产生不同的报告
-FIXA：当glorg中使用equates时，使先验坐标和速度连续（如果apr_file提供了连续值，在复杂的重命名中有时会失败）
-VSUM：速度的经纬总结（需要绘制速度）
-PSUM：经纬位置总结
-GDLF：包含h文件清单和运行中增加的卡方值
-CMDS：显示globk控制文件到输出文件中
地震或仪器变化时的处理步骤
Level 1（通常是必须的）
重命名测站，自动地：
eq_def WC 31.099 103.279 1000 20 2008 5 12 6 28
eq_rename WC forces rename, e.g. CHDU_GPS --> CHDU_GWC
eq_cosei to specify spatially dependent contraints on position adjustments 具体指定：
rename iisc_gps iisc_1ps 1995 10 1 2 0 0 1999 12 1 5 0
Antenna swap from Trimble SST to AOA choke ring
rename lake_gps lake_xhi 2002 1 12 0 0 2002 1 12 24 0
通常是写入eq_file中，不是globk控制文件，在解算中可以有效地使测站速度和位置独立
Level 2 （期望的）
在glorg中同等看待速度，具体地
equate iisc_gps ndot iisc_1ps ndot
equate iisc_gps edot iisc_1ps edot
equate iisc_gps udot iisc_1ps udot
自动地：
eq_dist 1000 ndot
eq_dist 1000 edot
eq_dist 1000 udot
unequate chdu_gps ndot chdu_gwc ndot
unequate chdu_gps edot chdu_gwc edot
unequate chdu_gps udot chdu_gwc udot
有效地链接到平差（应该在打印选项中使用FIXA以确保等同估计速度）
能够使用constrain创建软链接
Level 3 （有效改善far-field速度）
当在EQ半径内的测站有小错位时equate the positions
equate xian_gps npos xian_gwc npos
equate xian_gps epos xian_gwc epos
equate xian_gps upos xian_gwc upos
或许在连接一个模型时有用，重命名命令中应用偏移
rename xian_gwc xian_gwc 2008 5 12 6 28 -0.003 0.004 0.001 NEU
或者在apr_file中改变坐标
程序流
1、读取所有h文件头确定它们的内容（测站，其它参数，历元范围）
2、在eq_file中要求应用的重命名
3、按时间分类h文件清单（srt_dir）
4、用先验约束初始化卡尔曼滤波（apr_xxx）
5、一次读取一个h文件，计算卡方值增加，坐标调整，新数据显示的旋转；如果在限差范围内（max_chii），更新解算，将卡方值增加写入到日志文件中
6、将解算写入到sol_file和prt文件中（对于一个新h文件是可选的）
7、可选glorg应用概括约束
-应用约束（迭代固定0
-应用参数联合（equate，constrain，force），计算每一个卡方值增加值
-估计板块旋转（plate命令）
-将解算写入org文件（glorg prt文件）
globk不能做到的处理
1、修复原始分析中的错误：周跳和错误的天线相位中心模型
2、解算模糊度：将使文件太大
3、克服非线性影响：gamit中平差必须小于30cm
*globk可以删除测站：有效避免污染解算
globk处理中的apr文件
1、gamit
修复周跳后所有测站的精度为10m
最后线性平差小于30cm
模糊度解算时约束测站5cm
2、globk
如果执行glorg定义参考框架，globk中apr_file是可选的如果不执行glorg，需要约束测站的精确apr_file记录
对于复杂的重命名和equates，globk需要apr_file
3、glorg
apr_file只需要参考站和equates的坐标
控制网平差报告 [控制网概况] 计算软件：南方平差易2005 网名：成都学院1-6教导线网 计算日期：2014-07-02 观测人：赵磊 记录人：薛佳丽 计算者：薛佳丽 检查者： 测量单位：成都学院测绘工程一班 备注：第六小组 平面控制网等级：国家四等，验前单位权中误差：2.50(s) 已知坐标点个数：2 未知坐标点个数：8 未知边数：9 最大点位误差[D] = 0.0118 (m) 最小点位误差[B] = 0.0068 (m) 平均点位误差= 0.0109 (m) 最大点间误差=0.0102(m) 最大边长比例误差= 53 平面网验后单位权中误差=1.88(s) [边长统计]总边长：1211.146(m),平均边长：134.572(m),最小边长：51.705(m),最大边长：262.760(m) [闭合差统计报告]
高程网平差 -------------------------------------------------------------------- APPROXIMATE HEIGHT -------------------------------------------------------------------- No. Name Height(m) -------------------------------------------------------------------- 1 A1 500.0000 2 A2 499.6860 3 A3 499.3690 4 A4 499.2295 5 B1 497.9570 6 B2 497.1505 7 C1 495.7295 8 C2 495.0625 9 D1 495.5515 10 D2 494.9110 11 E1 494.5825
近似坐标计算的函数-calcux0y0函数（126页） function [x0,y0]=calcux0y0(x0,y0,e,d,sid,g,f,dir,s,t,az,pn,xyknow,xyunknow,point,aa,bb,cc) %本函数的作用是计算待定点的近似坐标 format short; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% time=0;prelength=length(xyknow);non_orient=0;point_angle=0; while length(xyunknow)>0 %考虑的计算方法有：1.极坐标;2.前方交会；3.测边交会；4.后方交会； %5.无定向导线的两种情况：（1）已知两个点；（2）分离的已知点与方位角；基本思路：%采用循环的方法逐一对每一个未知点进行以上各种方法条件的搜索，满足后即解算。 aa0=[];bb0=[];cc0=[];%记录搜索到两条观测边但需用户给顺序的点，注意要放在while 里面。 time=time+1; % 用于统计循环次数。 way=0; for i=xyunknow %依次循环向量中的各元素 %============================================================= ===== %方法1.极坐标条件搜索与计算-->way=1,基本思路：找到或求出一个方位角，找出一条边。 temp1=[]; temp2=[]; temp3=[]; temp4=[]; temp5=[]; temp6=[]; temp7=[]; temp8=[]; temp9=[]; temp10=[];A=[];B=[];P=[]; %第一步：寻找观测条件：两种情况：一是有已知方位角；二是由两个已知点及方向观测值推出方位角。 temp7=find(t==i); if length(temp7)>0 temp8=find(xyknow==s(temp7(1))); if length(temp8)>0 temp9=find(e==s(temp7(1))&d==t(temp7(1))|e==t(temp7(1))&d==s(temp7(1))); if length(temp9)>0 %第一种情况：有已知方位角（一般适用于闭合导线） S=mean(sid(temp9)); Alfa0=az(temp7(1)); x0(i)=x0(s(temp7(1)))+S*cos(Alfa0); y0(i)=y0(s(temp7(1)))+S*sin(Alfa0); way=1;method(i)=11;%----->由已知方位角算出 end end else temp1=find(f==i);%第二种情况：由两个已知点及方向观测值推出方位角（适用于附合、支导线） if length(temp1)>0 for j=xyknow
FJ－3 省道S229南坑至源头段 二级公路改建工程 GPS静态控制测量网平差报告 萍乡公路勘察设计院 二○一一年九月 目录 一、GPS控制点成果表 (1) 二、GPS控制点网示意图 (1) 三、GPS控制网平差报告……………………………………1~4
一、G PS控制点成果表 二、GPS控制点网示意图 三、GPS控制网平差报告 1 坐标系统 1.1 坐标系统名称 Beijing54 1.2 基准参数
1.3 投影参数 M0 =1.00000000 投影比率 H = 0.0000 投影高 Bm =0投影面的平均纬度 B0 =0:00:00.00N 原点纬度 L0 =113:50:00.00E 中央子午线 N0 =0.0000 北向加常数 E0 =500000.0000 东向加常数 回到顶部 2 三维无约束平差2.1 平差参数 2.2 基线向量及改正数 2.3 τ(Tau)检验表 2.4 τ(Tau)检验直方图
2.5 自由网平差坐标 回到顶部 3 二维约束平差 3.1 平差参数 3.2 平面距离平差值 3.3 平面坐标 ***** 回到顶部
4 高程拟合 4.1 平差参数 4.2 高程拟合坐标 240.7246 回到顶部 5 基线闭合差 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G2->G3.242A 99.9 0.0062 -2063.4456 -1777.5444 1294.6074 3015.5398 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.76ppm EX = 0.0043 EY = -0.0043 EZ = -0.0026 8706.0493 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G2->G4.242B 99.9 0.0072 -4060.9524 -3093.9755 2049.7944 5501.4248 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.48ppm EX = -0.0041 EY = 0.0051 EZ = 0.0010 13683.0814 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G2->GD1.242X 99.9 0.0065 1554.7134 -896.8104 2732.5118 3269.2543 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.80ppm EX = -0.0048 EY = 0.0042 EZ = 0.0017 8261.4927 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G3->G4.242B 99.9 0.0063 -1997.5067 -1316.4322 755.1870 2508.6519 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.12ppm EX = -0.0003 EY = 0.0004 EZ = 0.0015 13695.9047 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G3->GD1.242X 99.9 0.0071 3618.1569 880.7382 1437.9069 3991.7835 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.42ppm EX = 0.0026 EY = -0.0040 EZ = -0.0015 11989.6182 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G4->GD1.242X 99.9 0.0073 5615.6650 2197.1667 682.7190 6068.7182 G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.16ppm EX = 0.0015 EY = -0.0007 EZ = -0.0022 16557.6999
西南林业大学 《全球卫星定位系统原理》GPS控制网平差实习 （2012级） 题目静态GPS控制网平差总结报告 学院土木工程学院 专业测绘工程 学号20120456023 学生姓名施向文 任课教师朱毅 西南林业大学土木工程学院测绘工程系 2015年07月 12 日
目录 1 实习目的 0 2 实习任务 0 3 数据处理依据 0 4 精度要求 0 5 已有成果数据 0 6 数据处理过程 (1) 6.1创建作业及数据导入 (1) 6.2基线预处理 (1) 6.2.1静态基线处理设置 (1) 6.2.2处理基线 (2) 6.2.3搜索闭合环 (2) 6.3设置坐标系 (2) 6.4网平差 (2) 6.5高程内外符合精度检验 (3) 6.5.1内符合精度 (3) 6.5.2外符合精度 (3) 7 数据处理成果 (3) 7.1二维平面坐标平差 (3) 7.1.1 平差参数 (3) 7.1.2 平面坐标 (4) 7.2高程拟合 (7) 7.2.1 平差参数 (7) 7.2.2 外符合精度 (7) 7.2.3内符合精度 (9) 8 质量简评 (11) 9 总结 (12)
静态GPS网平差总结报告 1 实习目的 通过对静态GPS控制网的数据处理，从实践中加深对理论知识的理解。通过本次实习还可以熟悉GPS数据处理软件，现在的数据处理基本用软件处理，使用软件也是必备的一个技能。 2 实习任务 本次实习的任务： （1）静态GPS外业数据基线预处理，预处理基线的方差比应尽量调整在99.9，处理后搜索闭合环要基本合格。 （2）选择/建立坐标系，建立昆明87坐标系。 （3）输入已知点并进行网平差，检测内外符合精度。 （4）撰写数据处理总结报告。 3 数据处理依据 依据《卫星定位城市测量技术规范CJJ/T 73—2010》备案号J990—2010 4 精度要求 二维平差中误差1cm 高程拟合中误差2cm 高程内符合中误差3cm 高程外符合中误差5cm 5 已有成果数据 （1）静态GPS外业数据成果（RINEX） （2）已知点的三维坐标，坐标成果见下表
c o s w i n说明书平差软 件 HUA system office room 【HUA16H-
前言 “地面测量工程控制与施工测量内外业一体化和数据处理自动化系统”（简称科傻系统）将测量基本原理和现代科技相结合，对电子全站仪、电子水准仪以及常规地面测量仪器进行系统的开发，以地面控制测量、施工测量和碎部测量等测量工程为对象，实现从外业数据采集、质量检核、预处理到内业数据处理、成果报表输出的一体化和自动化作业流程。 该系统由两个子系统组成：“基于掌上型电脑的测量数据采集和处理系统”（简称COSA-HC），在掌上型电脑RD-EB2上运行，能自动控制和引导整个作业过程并进行质量检测，一体化程度高，操作方便。该子系统具有水准测量、二、三维控制、碎部测量、道路测设、工程放样等测量作业模块；具
有小规模水准网、二、三维工程网的平差功能；具有文件管理和数据通信功能；该系统灵活方便，适合外业环境。 “地面测量工程控制测量数据处理通用软件包”（简称CODAPS或COSAWIN）在微机WINDOWS环境下运行即可独立使用，也可与COSA-HC联合使用，对RD-EB2传输过来的原始观测数据进行转换，完成从概算到平差的数据自动化处理，同时具有粗差探测与剔除、方差分量估计、闭合差计算、贯通误差影响值估算、报表打印、网图显绘、坐标转换与换带计算、控制网优化设计以及叠置分析等功能。 本手册是为COSAWIN用户专门编写的，若有疏漏和不当之处，敬请读者提出宝贵意见和批评指正。 武汉测绘科技大学武地课题组
2000.5.
第一章概述 1.1 系统简介 科傻系统（COSA）是“地面测量工程控制与施工测量内外业一体化和数据处理自动化系统”的简称，包括COSAWIN和COSA-HC两个子系统。COSAWIN在IBM兼容机上运行。 COSAWIN系统除具有概算、平差、精度评定及成果输出等功能外，还提供了许多实用的功能，如网图
CPIII数据处理系统 （使用说明） 一．安装说明 1.安装硬件狗驱动。点击狗驱动文件夹下的DogInst.exe文件， 在弹出对话框中选择“安装驱动”，即可完成驱动程序的安装。 2.安装软件。插上硬件狗，然后打开setup文件夹，双击setup.exe 文件，按照提示操作即可完成软件的安装。 3.安装完成后需要手动在C盘根目录下新建一个名为“mytemp” 文件夹，然后在D盘根目录下新建一个名为“项目管理中心” 的文件夹即可，安装过程完成。 二．使用说明 1.插上硬件狗，双击CPIIIMAIN.EXE执行文件，启动主程序。 如图2.1.组侧为工具栏，控制平差流程，右侧为项目管理栏， 对项目进行管理。
图2.1 2.新建项目。点击该工具条即新建一个项目，给该项目命名并确 定后，程序将自动在D:/项目管理中心目录下新建一个以该项目名为文件名的文件夹，而后所有对该项目的操作都在该项目下进行，生成文件都保存在该文件夹下。 3.添加观测值。点击该工具条即显示出查找路径对话框，查找观 测原始数据路径，找到数据后点击确定即可将观测数据添加到项目中。如图2.2 图2.2 4.添加控制点。点击该工具条即显示查找路径对话框，查找控制 点数据，点击确定即可将控制点坐标数据添加到项目中。（注意：已知点数据文件格式为.txt格式，文本内容顺序为：点名，X坐标，Y坐标。如图1.3）
图2.3 5.观测边角计算。点击该工具条对观测文件进行边角关系计算， 生成的边角关系文件保存在右侧：项目管理中心/其他格式下 科傻格式，点右键用记事本打开即可查看。（提示：CPII导线 平面数据也可用此方法，生成科傻格式后可用科傻平差）三．数据分析 1.录入点号检错。点击该工具条对导入数据的点号进行检查，可 显示每个点被观测几次，分别在哪几站被观测，放便查找有点 号输入错误的。如图3.1 图3.1
在图 表9-1 试求： （１）1P 、2P 及3P 点高程之最或然值； （２）1P 、2P 点间平差后高差的中误差。 解:(1)列条件方程式,不符值以“mm ”为单位。 已知3,7==t n ，故437=-=r ，其条件方程式为 ??? ? ? ?? =--+=-+--=-+--=++-01030707742643765521v v v v v v v v v v v v （2）列函数式： 555v h x F +== 故 15=f 0764321======f f f f f f （3）组成法方程式。 1）令每公里观测高差的权为1，按1/i i s p =,将条件方程系数及其与权倒数之乘积填于表9-2中。 2)由表9-2数字计算法方程系数,并组成法方程式:
????????????----------5221251021411013????????????d c b a k k k k +????? ???????---1377=0 表9-2 条件方程系数表 （4）法方程式的解算。 1）解算法方程式在表9-3中进行。 2）[]pvv 计算之检核。 [][]wk pvv -= []467.35=-wk 由表9-3中解得[]47.35-=pvv ，两者完全一致，证明表中解算无误。 （5）计算观测值改正数及平差值见表9-4。 （6）计算321,,P P P 点高程最或然值。 359.3611=+=x H H A P m 012.3722=+=x H H A P m
表9-4 改正数与平差值计算表 （7）精度评定。 1）单位权（每公里观测高差）中误差 2）21,P P 点间平差后高差中误差 mm 0.34 47 .35±=±=μ
FJ －3 工程测量技术交流群18874248 省道S 229南坑至源头段 二级公路改建工程 GPS 静态控制测量 网平差报告 萍 乡 公 路 勘 察 设 计 院 二○一一年九月 目 录 一、 GPS 控制点成果表…………………………………………1 二、 GPS 控制点网示意图………………………………………1 三、 GPS 控制网平差报告……………………………………1~4
一、G PS控制点成果表 二、GPS控制点网示意图 三、GPS控制网平差报告 1 坐标系统 1.1 坐标系统名称 Beijing54 1.2 基准参数
1.3 投影参数 M0 =1.00000000 投影比率 H = 0.0000 投影高 Bm =0投影面的平均纬度 B0 =0:00:00.00N 原点纬度 L0 =113:50:00.00E 中央子午线 N0 =0.0000 北向加常数 E0 =500000.0000 东向加常数 回到顶部 2 三维无约束平差2.1 平差参数 2.2 基线向量及改正数 2.3 τ(Tau)检验表 2.4 τ(Tau)检验直方图
2.5 自由网平差坐标 回到顶部 3 二维约束平差 3.1 平差参数 3.2 平面距离平差值 3.3 平面坐标 ***** 回到顶部
4 高程拟合 4.1 平差参数 4.2 高程拟合坐标 240.7246 回到顶部 5 基线闭合差 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G2->G3.242A 99.9 0.0062 -2063.4456 -1777.5444 1294.6074 3015.5398 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.76ppm EX = 0.0043 EY = -0.0043 EZ = -0.0026 8706.0493 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G2->G4.242B 99.9 0.0072 -4060.9524 -3093.9755 2049.7944 5501.4248 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.48ppm EX = -0.0041 EY = 0.0051 EZ = 0.0010 13683.0814 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 G1->G2.242A 99.9 0.0077 -1046.7333 -648.5635 534.7004 1342.4566 G2->GD1.242X 99.9 0.0065 1554.7134 -896.8104 2732.5118 3269.2543 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.80ppm EX = -0.0048 EY = 0.0042 EZ = 0.0017 8261.4927 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G3->G4.242B 99.9 0.0063 -1997.5067 -1316.4322 755.1870 2508.6519 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.12ppm EX = -0.0003 EY = 0.0004 EZ = 0.0015 13695.9047 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G3->GD1.242X 99.9 0.0071 3618.1569 880.7382 1437.9069 3991.7835 G1->G3.242A 99.9 0.0068 -3110.1745 -2426.1123 1829.3052 4348.0529 G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.42ppm EX = 0.0026 EY = -0.0040 EZ = -0.0015 11989.6182 Baseline Type rms dx dy dz distance ------------------------------------------------------------------------------------------- G4->GD1.242X 99.9 0.0073 5615.6650 2197.1667 682.7190 6068.7182 G1->G4.242B 65.6 0.0072 -5107.6816 -3742.5441 2584.4937 6839.1999 G1->GD1.242X 99.9 0.0087 507.9850 -1545.3781 3267.2106 3649.7818 ------------------------------------------------------------------------------------------- 同步环( 3 baselines) 相对误差= 0.16ppm EX = 0.0015 EY = -0.0007 EZ = -0.0022 16557.6999
华南农业大学11测绘GPS控 制网测量 GPS网平差结果 施工单位：11测绘第二小组 负责人：梁永健 观测时间：2013年06月18日
观测数据： 文件名观测日期开始结束点名天线高天线高机号60181691.STH 2013年06月18日 09时37分 12时13分 6018 1.219 1.143 H1186209758 60181693.STH 2013年06月18日 12时46分 13时38分 6018 1.512 1.501 S5******* 60191691.STH 2013年06月18日 09时42分 10时30分 6019 1.359 1.348 S0******* 60501692.STH 2013年06月18日 11时24分 13时38分 6050 1.472 1.395 H1186208939 CL011691.STH 2013年06月18日 09时37分 10时32分 CL01 1.492 1.415 H1186208939 CL021691.STH 2013年06月18日 09时35分 10时39分 CL02 1.535 1.524 S0******* CL031692.STH 2013年06月18日 11时03分 13时38分 CL03 1.329 1.318 S0******* CL041693.STH 2013年06月18日 12时51分 13时59分 CL04 1.156 1.080 H1186209758 GP041692.STH 2013年06月18日 11时18分 13时40分 GP04 0.236 0.236 S0******* XX021691.STH 2013年06月18日 09时07分 12时12分 XX02 0.106 0.106 S5******* 基线解算： 60181691-60191691 观测量L1 L2 P2同步时长 48分钟历元间隔:30 高度截止角:20.0 三差解 0.000 453.731 87.026 247.558 524.147 1/1 双差浮点解 0.040 453.925 86.229 247.285 524.054 1/13095 双差固定解 28.81 0.027 453.992 86.247 247.310 524.127 1/19328 60181691-60501692 观测量L1 L2 P2同步时长 49分钟历元间隔:30 高度截止角:20.0 三差解 0.009 -355.578 -126.826 -49.205 380.712 1/44764 双差浮点解 0.012 -355.579 -126.900 -49.251 380.743 1/30669 双差固定解 89.39 0.013 -355.650 -126.924 -49.253 380.818 1/29569 60501692-60181693 观测量L1 L2 P2同步时长 52分钟历元间隔:30 高度截止角:20.0 三差解 0.007 355.550 126.914 49.256 380.722 1/58172 双差浮点解 0.010 355.550 126.859 49.249 380.703 1/37173 双差固定解 25.33 0.011 355.669 126.916 49.261 380.834 1/34257 60181691-CL011691 观测量L1 L2 P2同步时长 55分钟历元间隔:30 高度截止角:20.0 三差解 0.006 882.655 235.299 378.540 988.806 1/168531 双差浮点解 0.009 882.616 235.317 378.552 988.780 1/107755 双差固定解 33.68 0.010 882.609 235.365 378.590 988.800 1/96010 CL011691-60191691 观测量L1 L2 P2同步时长 48分钟历元间隔:30 高度截止角:25.0 三差解 0.014 -428.657 -149.084 -131.348 472.467 1/32787 双差浮点解 0.037 -428.725 -149.166 -131.311 472.544 1/12944 双差固定解 26.65 0.024 -428.614 -149.113 -131.277 472.418 1/19988 60181691-CL021691 观测量L1 L2 P2同步时长 61分钟历元间隔:30 高度截止角:20.0 三差解 0.006 930.849 131.704 584.022 1106.755 1/170886 双差浮点解 0.012 930.812 131.713 584.021 1106.725 1/92845 双差固定解 37.77 0.013 930.816 131.747 584.051 1106.748 1/86209 CL021691-60191691 观测量L1 L2 P2同步时长 48分钟历元间隔:30 高度截止角:20.0 三差解 0.016 -476.686 -45.576 -336.630 585.344 1/36483 双差浮点解 0.037 -476.957 -45.540 -336.753 585.632 1/15768 双差固定解 28.96 0.024 -476.822 -45.494 -336.730 585.505 1/23925 CL021691-CL011691 观测量L1 L2 P2同步时长 55分钟历元间隔:30 高度截止角:20.0
10.5.1基线解算 1．观测值的处理 GPS基线向量表示了各测站间的一种位置关系，即测站与测站间的坐标增量。GPS基线向量与常规测量中的基线是有区别的，常规测量中的基线只有长度属性，而GPS基线向量则具有长度、水平方位和垂直方位等三项属性。GPS基线向量是GPS同步观测的直接结果，也是进行GPS网平差，获取最终点位的观测值。 若在某一历元中，对k颗卫星数进行了同步观测，则可以得到k-1个双差观测值；若在整个同步观测时段内同步观测卫星的总数为l则整周未知数的数量为l-1。 在进行基线解算时，电离层延迟和对流层延迟一般并不作为未知参数，而是通过模型改正或差分处理等方法将它们消除。因此，基线解算时一般只有两类参数，一类是测站的坐标参数，数量为3；另一类是整周未知数参数（m为同步观测的卫星数），数量为。 2.基线解算 基线解算的过程实际上主要是一个平差的过程，平差所采用的观测值主要是双差观测值。在基线解算时，平差要分三个阶段进行，第一阶段进行初始平差，解算出整周未知数参数的和基线向量的实数解（浮动解）；在第二阶段，将整周未知数固定成整数；在第三阶段，将确定了的整周未知数作为已知值，仅将待定的测站坐标作为未知参数，再次进行平差解算，解求出基线向量的最终解-整数解（固定解）。 （1）初始平差 根据双差观测值的观测方程（需要进行线性化），组成误差方程后，然后组成法方程后，求解待定的未知参数其精度信息，其结果为： 待定参数： 待定参数的协因数阵：， 单位权中误差：。 通过初始平差，所解算出的整周未知数参数本应为整数，但由于观测值误差、随机模型和函数模型不完善等原因，使得其结果为实数，因此，此时与实数的整周未知数参数对应的基线解被称作基线向量的实数解或浮动解。 为了获得较好的基线解算结果，必须准确地确定出整周未知数的整数值。 （2）整周未知数的确定 第二节已提及，此处不再详述。 （3）确定基线向量的固定解 当确定了整周未知数的整数值后，与之相对应的基线向量就是基线向量的整数解。 10.5.2 基线解算的分类 1．单基线解算 （1）定义 当有台GPS接收机进行了一个时段的同步观测后，每两台接收机之间就可以形成一条基线向量，共有条同步观测基线，其中可以选出相互独立的条同步观测基线，至于这条独立基线如何选取，只要保证所选的条独立基线不构成闭合环就可以了。这也是说，凡是构成了闭合环的同步基线是函数相关的，同步观测所获得的独立基线虽然不具有函数相关的特性，但它们却是误差相关的，实际上所有的同步观测基线间都是误差相关的。所谓单基线解算，就是在基线解算时不顾及同步观测基线间的误差相关性，对每条基线单独进行解算。（2）特点 单基线解算的算法简单，但由于其解算结果无法反映同步基线间的误差相关的特性，不利于后面的网平差处理，一般只用在较低级别GPS网的测量中。 2．多基线解算
关于平差易（PA2005） 平差易（Power Adjust 2005，简称PA2005），它是在Windows系统下用VC 开发的控制测量数据处理软件，也是南方测绘PA2002的升级产品。它一改过去单一的表格输入，采用了Ｗindows风格的数据输入技术和多种数据接口（南方系列产品接口、其他软件文件接口），同时辅以网图动态显示，实现了从数据采集、数据处理和成果打印的一体化。成果输出丰富强大、多种多样，平差报告完整详细，报告内容也可根据用户需要自行定制，另有详细的精度统计和网形分析信息等。其界面友好，功能强大，操作简便，是控制测量理想的数据处理工具。 如何安装平差易 平差易可在Windows95、98、2000和XP下安装运行。 安装步骤 平差易（PA2005）的安装光盘中有PA2005文件夹，打开此文件夹并找到setup.exe 文件，双击setup.exe后屏幕上将出现下图的界面。 平差易的安装准备 等待2秒钟，平差易的安装准备完成后即进入平差易的安装。
点击“Next”进入软件安装“用户须知”界面 如果同意安装许可，请点击“yes”。
安装路径 在上图“安装路径”中设置PA2005的安装目录。安装软件给出了默认的安装位置c:\Program Files\South Survey Office\Power Adjust，用户也可以通过单击Browse 按钮从弹出的对话框中修改软件的文件夹。如果已选择好了文件夹，则可以单击Next 按钮开始进行安装。此时平差易的主体程序已安装完毕。 安装完成后屏幕弹出以下界面，确定是否要安装软件狗的驱动程序。
作用： 管网平差是在上面已知条件节点流量、管道直径和某些点压力求得其它节点的水压，各个管段的流量以及由此推算得出的各节点自由水头、各管段的流速、水头损失等。 它的作用有两个：1、根据最不利点压力反算水源点压力，用于确定水源点的水泵供水压力，进而选择水泵。2、根据水源点压力计算最不利点压力，用于校核水源压力是否能满足最不利点的压力要求。 鸿业操作步骤： 1．设置－出图比例20000，设置－文字大小－2.5 2．设置一个tempt（图名，可随便定）图层，颜色为紫色，线宽为0.3（这些可以随便定，主要为了确定给水管道位置，但这个步骤最好有） 3．在这个图层上画出图1中的框框， 1和10节点间有两条管线所以另外加一条折线增加一个节点11。 4．给水－交互布管－新建工程名（名称随意）－确定－管代号（输J），选球墨铸铁管－确定－用鼠标点击节点布置管线，顺序可随意，布完一个环用回车键布其余管段。 5．工具－图层－点击新建图层前的小灯泡（这一步是为了关闭新建图层） 6．给水－节点编号－给水节点自动编号－选择起始节点，在图上点击－输入节点其始编号1，然后回车，全部按主结点编号，回车，全部自动编号，回车，如果出现有些节点未标注，是否自动标注，输Y，回车。（节点编号与所要求的不一致没关系，在下一步中更改） 7．给水－编辑标注－查询－用鼠标点选节点，修改节点编号，输入地面标高。或者采用工具－编辑查询。 8．平差－定义管长－点击管段J10－J11，J11－J1，输600（这是为了使输水管总长为1200）9．平差－定节点流量－定管供水类型－按照提示定义1－4，4－7为单侧供水，10－1，10－11，11－1为不供水 10．平差－定节点流量－定集中流量－输入求得的水源流量，前面加负号 11．平差－定节点流量－按管长分配流量－输H 12．平差－定义管径－定义管线直径。 13．平差－选择反算水源压力，自由水头28，计算公式：海曾威廉，计算温度：20，管道局部损失系数：1.05,球墨铸铁管。点击图面提取，全选图形，回车，输入h，回车，点击节点参数，在控制水压下面，第10个节点处输入根据泵的特性曲线方程求得的水
控制网平差软件介绍 NASEW概况 NASE始于1992年8月的1.0版, 到1992年12月1.6版通过省级鉴定. 1993年5月的2.0版经测绘出版社正式出版, 版本号: ISBN 7-5030-0708-7/P.249 NASEW 95 (v1.0) 是在原NASE v2.40A 的基础上, 结合测绘行业和计算机技术的新发展, 于1995年推出的WINDOWS 版本. 自1992年NASE在许多测绘单位使用以来收到了大量的反馈信息和改进意见, 作者在NASEW 的研制中充分考虑了这些宝贵的意见, 谨借此表示诚挚的谢意! NASEW与NASE数据文件兼容, 保留了原NASE中诸优点, 同时在实用性和易用性上又下了大量的功夫, 且 NASEW与支持各种全站仪进行等级测量的<ELER 电子测量记簿> 保持了程序级的连接, 真正实现了从数据采集, 记簿整理, 平差和成果打印的一体化. NASEW具有如下特点: 1. 适用于任意形式任意规模的平面和高程控制网的概算, 平差和设计, 无需编码. 2. 自动求解控制网各种路线闭合差, 提供了多种粗差定位和自动剔除功能. 具有如验后定权法等多种平差方法可选. 3. 具有现代电子表格效能的数据编辑和操作环境, 输入中坐标, 高程, 差值等自动计算, 并辅以网图动态显示. 观测输入具有标准格式和多种常用网格式可选, 还可自定格式. 具有内联的文本编辑方式可选. 4. 提供了与打印机和纸张自适应的网图打印, 成果打印, 格式和有效位数等可控易控, 并具有打印前的预显功能. 5. 操作上简单易学, 图, 文, 数, 控四者融为一体, 具有丰富的说明和联机帮助. 具有工作现场保护功能. 6. 互用性强, 支持网内多区域数据合并. 可直接读入ELER电子测量记簿的整理数据, 可直接更新EPSW电子平板的KNW文件, 可直接读入武测平差系统的数据文件. NASEW 是商品化的测量专业软件, 保证合理的后续服务, 欢迎用户提出宝贵的意见和建议. NASEW入门 1. NASEW 的运行环境 NASEW 要求在WINDOWS 3.1 上运行, 所有提示, 说明, 帮助均为中文, 所以要求WINDOWS 系统具有中文环境. NASEW 的所有操作和命令既支持键盘又支持鼠标, 鼠标操作更加方便灵活. 2. 控制网处理过程的支持 数据录入, 概算, 闭合差计算, 平差和精度评定是控制网处理的基本过程: 1). 数据录入: 分文件数据读入和直接键入两种. 凡符合NASEW文件系统的数据均可读入, 无需手工键入, 特别是由ELER 电子记簿采集的数据