YG测绘数据处理系统操作手册

YG平差工具箱2009版

本软件使用Visual Basic 6.0 软件编写，是集合在测绘工作中常用的平差计算、投影换代、坐标转换计算等一体的计算软件，内置程序共22个，程序经过调试运行，计算成果准确可靠。程序中的有关各项限差均以《城市测量规范》CJJ8-99中的规定计算。

为方便使用，本软件中的部分计算程序均设计有数据文件，数据文件包括：原始数据文件、平差或计算成果文件，原始数据文件均为.txt格式，平差或计算成果文件均为.doc格式（需先安装Word）。

一 高程控制平差

本菜单包括三个计算程序：四等水准近似平差、三角高程近似平差、水准网严密平差 1. 四等水准近似平差

（1）程序功能：此平差程序的平差方式为按每测段的距离进行配赋平差，可对闭合水准、附

合水准、支水准路线平差。对于闭合、附合水准路线，平差后计算其闭合差和闭合差限差，闭合差超限的数据不予以平差，且有提示对话框。

（2）数据文件：此程序设有原始数据、平差成果文件。可先编辑好数据文件，然后直接读取

即可，平差后可将平差成果保存为.doc格式文档，方便打印。

原始数据格式：文件原始数据文件格式如下所示（附合水准）：

A(已知),21.453,B(已知),25.006, A,1,3,3.748, 1,2,2,-12.365, 2,3,4,11.430, 3,B,5,0.715,

上述数据文件格式为.txt格式(记事本格式)，且为附合水准路线原始数据文件，原始数据所代表的意义如下，表格如图1所示：

第一行：已知点A，已知点A的高程21.453，已知点B，已知点B的高程25.006， 第二行：起点A，待定点1，线路长度3公里，高差3.748米， 第三行：待定点1，待定点2，线路长度2公里，高差-12.365米， ??

第五行：待定点3，已知点B，线路长度5公里，高差0.715米，

对于闭合水准路线计算，只需把第一行改为：A(已知),21.453,A(已知),21.453, （3）操作步骤：

第一步：点击四等水准近似平差进入计算界面，点击“计算设置”，按要求填写，填写完毕后

单击“确定”。

第二步：单击“调入表格”，可以手工输入观测数据，也可调用编辑好的原始数据，数据输入完毕后，单击“确定”。

第三步：单击“平差计算”，如下图：

第四步：单击“保存平差成果”，成果保存为Word文档，如下：

闭合水准近似平差报告

计算日期: 计算者:

---------------------------------------------------------------------------

[闭合差计算]

水准路线闭合差允许值fh允 =0.0849 (m) 水准路线闭合差fh =-0.0013 (m) 水准路线长度∑D =18.0000 (km)

---------------------------------------------------------------------------

[配赋平差表格]

点号 距离D(km) 高差 h(m) 改正数改正后高差高程 H(m)

v(m) h(m)

已知点A 1000.000

1 5 1.1234 0.0004 1.1238 1001.1238 2 6 1.2321 0.0004 1.2325 1002.3563 3 3 2.1251 0.0002 2.1253 1004.4816 已知点A 4 -4.4819 0.0003 -4.4816 1000.0000

（4）下图为附合水准和支水准计算界面，供参考：

2. 三角高程近似平差

（1）程序功能：此程序专为全站仪测量编写，可对闭合、附合、支三角高程导线进行平差计

算，由于全站仪可直接测量初算高差，初算高差为没有经过加入仪器高和目标高计算的高差，具体各个符号代表的意义在三角高程测量设置中均给予说明。

（2）操作步骤：三角高程观测必须是对向观测，程序应用具体步骤如下：

第一步：按软件要求输入观测数据后，如图4，单击“计算

正数、往返高差较差允许值、往返高差平均值，如图5；

图 4

图 5

第二步：单击“闭合差调整”按钮，出现调整闭合差界面，然后单击“读取数据”按钮，自

动导入第一步计算出的数据，再单击“计算”按钮，进行平差计算及闭合差的计算，如图6。

图 6

3. 水准网严密平差

（1）程序功能：本平差程序适用国家各等级水准网平差，可对各种水准网进行平差计算，采

用间接平差方法，利用线路长度定权，可以对观测高差加入正常位水准面不平行改正数、水准尺尺长改正数。由于本程序在编写过程中，程序中的各种过程性数据，在计算时没有取位和四舍五入的问题，并且在进行NBB求逆矩阵的时候，采用了全选主元的高斯—约当法，虽然计算过程复杂，但避免了求逆过程中舍入误差对结果的影响，保证了程序的计算精度。此程序设有原始数据文件和平差成果文件。

（2）计算方法：此程序的计算方法，依据武汉测绘科技大学出版《误差理论与测量平差基础》

书中的水准网间接平差计算。

（3）原始数据格式：

1、不加正常位水准面不平行改正的水准网原始数据文件，已知点为BM1，待定点为BM2、

BM3、BM4，有5段观测路线，有几段观测路线，就有几行数据，数据文件如下：

1,237.483,2,,5.835,3.5, 2,,3,,3.782,2.7, 1,237.483,3,,9.640,4.0, 4,,3,,7.384,3.0,

1,237.483,4,,2.270,2.5,

第一行：线路起点1，已知点高程237.483，线路终点2（待定点高程不写)，线路观测高差5.835

米，线路长度3.5公里，

第二行：线路起点2，起点高程不写，线路终点3（待定点高程不写），线路观测高差3.782

米，线路长度2.7公里，

第三行：线路起点1，已知点高程237.483，线路终点3（待定点高程不写），线路观测高差

9.640米，线路长度4.0公里，

??

第五行：线路起点1，已知点高程237.483，线路终点4（待定点高程不写），线路观测高差

2.270米，线路长度2.5公里，

2、加正常位水准面不平行改正的水准网原始数据文件，已知点为1，待定点为1、2、3，有

5段观测路线，有几段观测路线，就有几行数据，数据文件如下：

1,237.483,41.21,2,,41.22,5.835,3.5, 2,, 41.22,3,,41.30,3.782,2.7, 1,237.483,41.21,3,,41.30,9.640,4.0, 4,,41.28,3,,41.30,7.384,3.0, 1,237.483,41.21,4,,41.28,2.270,2.5,

第一行：线路起点1，已知点高程237.483，起点纬度41.21，线路终点2（待定点高程不写），

终点纬度41.22，线路观测高差5.835米，线路长度3.5公里，

第二行：线路起点2，起点高程不写，起点纬度41.21，线路终点3（待定点高程不写），终点

纬度41.30，线路观测高差3.782米，线路长度2.7公里，

第三行：线路起点1，已知点高程237.483，起点纬度41.21，线路终点3（待定点高程不写），

终点纬度41.30，线路观测高差9.640米，线路长度4.0公里，

??

第五行：线路起点1，已知点高程237.483，起点纬度41.21，线路终点4（待定点高程不写），

终点纬度41.28，线路观测高差2.270米，线路长度2.5公里，

（4）操作步骤：

第一步：平差设置，输入已知点数、总点数、总路线数、单位权观测高差的公里数，根据需

要选择是否加入正常位水准面不平行改正数、水准尺尺长改正数，当需要加入正常位水准面不平行改正时，须填写每个点所在位置的纬度，设置完毕以后，单击“确定”按钮，如图：

第二步：单击“调入表格”，出现水准网观测数据输入表格，数据的输入可以直接输入表格，

根据表头显示的内容，依次填写观测数据，特别注意：表格中的起点点号和终点点号必须是数字，且已知点点号在前，如有三个已知点，已知点点号应为1、2、3，四个待定点，待定点点号接着已知点点号，为4、5、6、7。也可以单击“读取数据”，读取已编好的原始数据文件，然后单击“确定”，如图：

第三步：当按要求填写完毕观测数据以后，单击“平差计算及精度评定”按钮，程序将计算

每测段的高差改正数、高差平差值，待定点近似高程、平差后高程，在精度评定方面，可计算待定点高程中误差，计算单位权中误差，如图：

第四步：如果是手工输入的观测数据，此时最好先单击“原始数据保存”，将原始观测数据保

存，方便下次使用。如果输入观测数据时使用的读取数据，就不需要再保存原始数据了，直接单击“平差成果保存”。

（5）成果文件：

水准网严密平差报告

---------------------------------------------------------------------------

[观测记录手薄]

项目名称: 计算者: 观 测 者: 仪器号: 仪器名称: 计算者: 高程系统: 计算日期:

---------------------------------------------------------------------------

[网平差设置]

已知点数:1 总点数:4 总路线数:5 平差前对观测高差加水准尺长改正数,改正数(mm)=0.2 平差前对观测高差加正常位水准面不平行改正数。 水准网线路总长∑D= 206.4(km) 1公里观测高差中误差=± 7.7(mm)

---------------------------------------------------------------------------

[观测高差平差结果]

起终尺长改不平行改正观测高差(m) 高差改正(mm) 高差平差值点 点 正(mm) (mm) (m) 1 2 1.2 -0.4 5.835 11.90 5.8469 2 3 0.8 -3.0 3.782 9.14 3.7911 1 3 1.9 -3.3 9.640 -2.00 9.6380 4 3 1.5 -0.7 7.384 -8.73 7.3753 1 4 0.5 -2.5 2.270 -7.31 2.2627

---------------------------------------------------------------------------

[待定点高程及其精度评定]

点名 高程近似值(m) 高程平差值(m) 高程中误差(m) 2 243.3180 243.3307 11.1 3 247.1230 247.1196 10.0 4 239.7530 239.7436 10.1

二 平面控制平差

本菜单包括两个计算程序：单导线近似平差、单导线严密平差 1. 单导线近似平差

（1）程序说明：此程序专为城市一级导线编写，各项限差按照《城市测量规范》CJJ8-99中

的规定计算。

（2）程序功能：此程序可对单一附合（闭合）导线、无定向导线、支导线近似平差，在第三

个选项卡中，附有各中导线的形式。程序计算后，可在表格中显示出计算过程，即闭合差配赋过程。

（3）操作步骤：

第一步：输入计算数据，表格中呈灰色单元格的为需要输入的观测数据和已知数据，已知数

据包括起算边的坐标方位角和连接点的已知坐标，观测数据按表头显示的内容依顺序填写，填写完毕后可单击“平差计算”，进行平差计算。此程序还可同时计算出导线的角度闭合差允许值、角度闭合差、导线全长绝对闭合差、导线全长相对闭合差。

第二步：在计算完毕后，单击“成果保存”，将平差后待定点坐标、闭合差项目保存在WORD

文档中，方便以后使用和打印成果。

2. 单一附合（闭合）导线间接平差

（1）程序功能：此程序专为城市一级导线编写，各项限差按照《城市测量规范》CJJ8-99中

的规定计算。程序平差方法为按观测角度平差，可计算闭合差项目，计算观测值的改正数精度评定计算。此程序界面由四个选项卡组成，使用之前先在第四个选项卡中了解程序规定导线的形式，以免在输入数据时发生错误。

（2）操作步骤：

第一步：单击“数据输入及平差计算”选项卡，填写导线严密平差设置中需要填写的数据（先

验方向中误差，对于城市一级导线为5秒）填写完毕后，单击“调入表格”，显示计算表格，观测数据和已知数据填写在表格中灰色的单元格中，已知数据为各已知点的坐标，程序规定起算边为AB边，且B点的点号为1，待定点为2 、3??。在附合导线中，CD为附合边，且C点的点号为：待定点点数+2。

第二步：待所有必要的数据输入完毕后，单击“平差计算及精度评定”按钮，计算出的数据

有：角度、边长的平差值，导线个待定点的近似坐标（直接用观测值计算而出的坐标），经间接平差后的待定点坐标，各待定点的平面点位误差和误差椭圆参数，闭合差项目，验后的方向中误差。

第三步：单击“观测手薄记录”选项卡，按要求填写各内容，可不填。 第四步：单击“成果保存” 按钮，将成果保存为WORD文档。

三 计算工具

本菜单包括五个程序：坐标正反算、单点坐标推算、度分秒转度、交会计算、i角计算 1. 坐标正反算

程序功能：此程序可进行坐标正算（已知一边的方位角、连接点的坐标、转角和连接点与待

定点之间的边长，求待定点坐标）和坐标反算（通过两个已知坐标，求已知边的方位角和距离）。

2. 单点坐标推算

程序功能：此程序为支点计算程序。 3. 度分秒转度

4. 交会计算

程序功能：此程序可对前方交会进行计算。 5. i角计算

程序功能：此程序可进行水准仪检较计算，计算方法为I1ABI2。程序计算方法参考《国家三、

四等水准测量规范》 GB12898-91 。

四 高斯投影、改化和坐标系转换

本菜单下有五个程序：高斯投影正反算、坐标系转换、WGS-84大地坐标转国家直角坐标、投影变形计算、椭球面到高斯投影面的改化。 1. 高斯投影正反算

（1）程序功能：此程序可进行高斯投影正算、高斯投影反算、换带计算。计算的坐标系统目

前有北京54坐标系和80西安坐标系两种，程序可计算大地坐标系以任意中央子午线在两种系统中的坐标，计算出的Y坐标可以选择是否加500公里，另外还可以计算出子午线收敛角r。

（2）计算方法：本程序应用的计算方法参考《控制测量学》下册，武汉大学出版社2006年

11月出版。

测绘通报1997年3月《一组高精度椭球面计算实用公式》 顾旦生 （3）原始数据格式： ①、高斯投影正算数据格式

III呼市,31.04416832,111.4724897,

III大台,30.45254425,111.1758359,

②、高斯投影反算数据格式

III呼市,3439978.970363,575412.861661,

III大台,3404139.839337,528680.555054,

③、换带计算数据格式

II大台,1945024.114,739233.054, （4）操作步骤：

第一步：点击高斯投影正反算进入计算界面，点击“计算设置”，按要求填写，如图11，填

写完毕后单击“确定”， 本例为高斯投影正算。

图11

第二步：调入表格，出现计算表格，输入已知数据。已知数据的输入方法有两种，一种

是手工输入，另一种将编辑好的数据文件直接调用。调入数据后如图12，如果 已知数据是手工数据，最好先单击“保存原始数据”，然后单击“确定”。

图12

第三步：转换计算，计算后如图13，然后单击“成果保存”，将计算成果保存为word文档。

图13

（4）成果保存

高斯投影正算成果 计算日期: 计算者:

---------------------------------------------------------------------------

[计算设置]

计算点数 = 2

坐标系统 ： 北京54坐标系 中央子午线经度（dms）= 111 横坐标带号： 3度带 36带 Y坐标加常数= 500公里

---------------------------------------------------------------------------

[已知点经纬度]

序 号 点 名 纬 度 B(dms) 经 度 L(dms) 1 III大正 31.04416832 111.4724897 2 III曲沟 30.45254425 111.1758359

---------------------------------------------------------------------------

[待定点坐标值]

序 号 点 名 坐标 X(m) 坐标 Y(m) 收敛角r(dms) 1 III大正 3439978.97041 36575412.86181 0.24286284 2 III曲沟 3404139.83928 36528680.55499 0.09114756

2. 平面坐标转换

（1）程序功能：在某些测区内，由于各单位所作控制要求不同，往往会存在几个不同的坐标

系统。有时为了相互利用成果，有时为了统一测区坐标系统，于是就产生了坐标转换问题。利用本程序可以实现任意两个坐标系相互转换，并求出转换元素和检核元素。

（2）计算说明：本程序应用的计算方法参考《控制测量计算手册》（冶金工业出版社 1972

年12月第一版）中282页“应用最小二乘法进行平面坐标换算”一节。当新旧坐标系重合点为2个时，采用重心坐标公式求解，当新旧坐标系重合点大于2个时，采用转换精度较高的最小二乘法求解。

（3）操作步骤：

第一步，转换设置，填写两坐标系重合点点数、待转换坐标点点数，“确定”，图14。

图14

第二步，调入表格，现在可以直接输入数据或者读取原始数据文件，待数据输入完毕后，最好先点击“保存原始数据”菜单，将数据的数据保存为.txt文本格式，方便重新计算和资料备份，如图15

右图15

第三步：转换计算，即计算出待转换点在新坐标系中的坐标值和转换参数，如图16：

右图16

第四步：保存转换成果，将转换后的成果保存为.doc Word文档格式，方便排版打印。

坐标转换成果报告

计算日期: 计算者:

--------------------------------------------------------------------------- 新旧坐标系重合点点数 =8 待转换坐标点点数 = 2

---------------------------------------------------------------------------

[转换参数及检核计算]

纵坐标X平移量(m) =4059.25816342101 横坐标Y平移量(m) =10125.7892878279 旋转角α(dms)= 6.33240 尺度k = 0.9998409989

重合点在新系统内的x不符值= 0.000 重合点在新系统内的y不符值= 0.000

---------------------------------------------------------------------------

[新旧坐标系重合点坐标值]

序号 点名 旧坐标X(m) 旧坐标Y(m) 新坐标x(m) 新坐标y(m) 1 A 2000.00 2000.00 5817.54 12340.67 2 B 2212.64 2353.61 5988.31 12716.39 3 C 2611.87 2711.97 6344.12 13117.67 4 D 2985.42 2524.12 6736.39 12974.01 5 E 2535.31 1725.17 6380.47 12128.71 6 F 2332.15 1649.85 6187.56 12030.66 7 G 1951.57 1711.72 5802.24 12048.91 8 H 1894.81 1765.73 5739.96 12096.14

---------------------------------------------------------------------------

[新旧坐标系待定点坐标值]

序号 点名 旧坐标X(m) 旧坐标Y(m) 新坐标x(m) 新坐标y(m) 1 Z1 2710.17 2063.91 6515.6423 12485.2877 2 Z2 2861.95 2359.18 6632.6954 12795.9082

3. WGS-84大地坐标转国家直角坐标

（1）程序功能：本程序可以实现WGS-84大地坐标(B,L)转换为国家直角坐标或地方坐标

(X,Y)。由WGS-84大地坐标转换北京54或国家80坐标，在重合点精度较高时，转换精度可以达到毫米位。

（2）计算说明：本程序应用的计算方法参考测绘通报1997年3月《一组高精度椭球面计算

实用公式》、《高精度底点纬度公式》、冶金工业出版社《控制测量计算手册》、武汉科技大学《GPS测量原理及应用》。

（3）操作步骤：

第一步：转换设置，填写两坐标系重合点点数、 待转换坐标点点数、中央子午线经度等， 单击 “确定”，如右图17。

第二步，调入表格，现在可以直接输入数据或者读取原始数据文件，待数据输入完毕后，最好先点击“保存原始数据”菜单，将数据的数据保存为.txt文本格式，如图18

图 18

第三步：转换计算，计算待转换WGS-84大地坐标的国家坐标，如图19：

图 19 （4）成果保存：

WGS-84大地坐标转国家直角坐标成果报告

计算日期: 计算者:

--------------------------------------------------------------------------- 两坐标系重合点点数 =45 待转换坐标点点数 = 5

---------------------------------------------------------------------------

[两坐标系重合点坐标值]

序号 点名 WGS-84大地坐WGS-84大地坐标国家平面坐标国家平面坐标

标纬度B(dms) 经度L(dms) X(m) Y(m)

1 A001 40.502174513 111.443384913 4523039.9842 562597.4353 2 A002 40.532418540 111.500618650 4528737.8294 570329.6085 3 A003 40.451793539 111.500051222 4513737.1139 570339.2895 4 A004 40.405818692 111.414959189 4505624.4657 558887.2656 5 A005 40.445244048 111.344065435 4512777.3742 548767.3427 6 C001 40.562300255 111.523780050 4534288.6302 573823.5715 7 C002 40.532336605 111.541829148 4528771.1745 576231.5497 8 C003 40.551466060 111.471027118 4532107.6233 566180.9279 9 C004 40.545865535 111.423141557 4531558.2037 559660.0726 10 C005 40.531671240 111.440108375 4528430.7937 561784.6638 11 C006 40.513218877 111.485362589 4525267.0261 568663.1261 12 C007 40.501469371 111.525085053 4522930.2420 574243.0694 13 C008 40.464655844 111.544235559 4516536.5452 576922.1462 14 C009 40.481203994 111.493313215 4519101.6467 569646.5157

15 C010 40.471907027 111.440004728 4517398.3393 561852.6363 16 C011 40.513827811 111.432317139 4525386.9789 560922.1768 17 C012 40.532074172 111.393448685 4528505.4760 555542.6451 ---------------------------------------------------------------------------

[两坐标系待定点坐标值]

序号 点名 WGS-84大地坐WGS-84大地坐标国家平面坐标国家平面坐标

标纬度

B(dms) 经度L(dms) X(m) Y(m)

1 A001 40.502174513 111.443384913 4523039.9849 562597.4343 2 A002 40.532418540 111.500618650 4528737.8338 570329.6127 3 A003 40.451793539 111.500051222 4513737.1118 570339.2881 4 A004 40.405818692 111.414959189 4505624.4653 558887.2647 5 A005 40.445244048 111.344065435 4512777.3744 548767.3436

4. 投影长度变形计算

程序功能：实现两种投影变形计算，即由实测边长归算到参考椭球面上的变形计算，由

参考椭球面上的边长投影到高斯平面上的变形计算。如下图20：

4. 椭球面到高斯投影面的改化

程序功能：椭球面到高斯投影面的改化：此程序的参考椭球面有克拉索夫斯基椭球和1975年

国际椭球两种，通过已知AB两点的坐标和纬度，计算出在高斯投影面的方向和距离改正数以及改化后AB两点间的距离，如下图21：

五 独立坐标系建立

本菜单下有三个计算程序：独立坐标系的选择、计算椭球常数、独立坐标系坐标计算。 1. 独立坐标系的选择

（1）程序功能：在《城市测量规范》中有规定：城市平面控制测量坐标系统的选择应以投影

长度变形值不大于2.5cm/KM为原则，并根据城市地理位置和平均高程而定。本程序提供三种方法以供选择：抵偿投影面的3度带高斯正形投影平面直角坐标系、任意带高斯正形投影平面直角坐标系、具有高程抵偿面的任意带高斯正形投影平面直角坐标系。

（2）计算说明：本程序参考测绘通报1997年10月《建立地方独立坐标系的方法》、测绘通

报2008年8月《独立网椭球变换与坐标转换的研究》，武汉科技大学出版《控制测量学》、武汉科技大学出版《GPS测量原理及应用》等。

（3）操作步骤：（具有高程抵偿面的任意带高斯正形投影平面直角坐标系做为算例） 第一步：填写已知数据，有：测区平均纬度、测区中心与国家中央子午线的距离、测区平均

大地高、投影椭球选择，如图22：

图 22

第二步：选择工程测量坐标系，如图23：

图 23

第三步：单击“智能计算”，计算出的坐标系建立条件是在能保证测区控制范围内的最佳条件，

如图24的：

图 24

第四步：如果想改变第三步中程序所选的条件，如高程抵偿面和中央子午线的距离等，可单

击“人工调整”菜单下的“改变条件”，输入改变后的数据，在单击“人工调整”菜单下的“调整后计算”，计算结果如图25：

图 25

2. 计算椭球常数

（1）程序功能：本程序提供两种椭球变换方法，分别为：大地高直接改变新椭球长轴、大地

高改变新椭球的平均曲率半径。通过克拉索夫斯基椭球或1975年国际椭球，利用椭球膨胀法，计算出新椭球的常数15个。计算出的新椭球常数供独立坐标系坐标计算所用。

（2）计算说明：同上。

（3）操作步骤：先选择椭球的变换方法，在输入已知数据，即可计算，因操作较简单，这里

不详述，操作结果如图26：

图 26

3. 独立坐标系坐标计算

（1）程序功能：本程序仅限于具有高程抵偿面的任意带高斯正形投影平面直角坐标系，计算

时需调用上一程序中计算出的椭球常数，实现北京54坐标或国家80坐标与独立坐标系坐标之间的转换。

（2）计算说明：同上。 （3）操作步骤：

第一步：单击“计算设置“，按要求填写已知数据，选择椭球常数文件，如图27：

图 27

第二步：，

，现在可以直接输入数据或者读取原始数据文件，待数据输入完毕 后，最好先点击“保存原始数据”菜单，将数据的数据保存为.txt文本格式，如图28：

图28

第三步：确定后单击“转换计算“，如图29：

图 29

（4）成果保存

独立坐标系坐标计算成果

计算日期: 计算者:

---------------------------------------------------------------------------

[计算设置]

独立坐标系建立方法： 具有高程抵偿面的任意带高斯正形投影平面直角坐标系。 转换方式： 1954北京系→独立坐标系 计算点数 = 2

中央子午线经度（dms）= 111

测区中央子午线经度（dms）= 113.03 测区平均高程面（m）= 1510.000

---------------------------------------------------------------------------

[独立椭球参数]

椭球长半轴a(m) = 6379795.0000000000000 椭球短半轴b(m) =6358407.8226617500000 椭球的扁率α = 1/298.3

椭球第一偏心率的平方e = 0.0066934216229659500 椭球第二偏心率的平方e′=0.0067385254146835000

椭球极点处的子午线曲率半径C(m) =6401254.11537168000

---------------------------------------------------------------------------

[已知点坐标]

序 号 点 名 坐标 X(m) 坐标 Y(m) 1 1 4547070.860 685545.770 2 2 4531491.400 661461.700

---------------------------------------------------------------------------

序 号 1 2

[待定点坐标]

点 名 坐标 X(m) 1 4545846.0123 2 4530837.8204

坐标 Y(m) 513158.1277 488718.2823

六 GPS高程拟合

（1）程序功能：由于GPS测量中地面点的高程H84采用的是WGS—84坐标系统中的大

地高，但在实际应用中，地面点的高程采用正常高系统Hr，所以必须将H84

转换为Hr 。由图可知，只要知道了各GPS点的大地高和高程异常值?，就可以将H84转换为Hr 。其换算关系为：

Hr = H84 - ?

（2）计算说明：本程序参考武汉科技大学出版《控制测量学》、武汉科技大学出版《GPS测

量原理及应用》、测绘通报1999年7月《GPS高程拟合中二次曲面的特征根和高斯曲率》、世界地质2006年6月《小区域GPS高程拟合方法精度研究》等。

（3）操作步骤：

第一步：单击“计算设置“，按要求填写已知数据，如图：

第二步：，调入表格，现在可以直接输入数据或者读取原始数据文件，待数据输入完毕 后，最好先点击“保存原始数据”菜单，将数据的数据保存为.txt文本格式，如图：

第三步：确定后单击“计算“菜单下的“拟合计算”，如图29：

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