第五章第五章 坐标系统转换和高程系坐标系统转换和高程系统转换程序设计统转换程序设计 本章重点：本章重点：一一. .测量中常用坐标系及这些坐标系的特测量中常用坐标系及这些坐标系的特点；测量中高程系统及其关系、点；测量中高程系统及其关系、GPSGPS高程高程二二. .程序实现时的技巧程序实现时的技巧 一、地理空间坐标系一、地理空间坐标系 为了确定物体在空间的位置，常用坐标系来描述空间位置测量中常用的坐标系有以下几种：（1）空间直角坐标系 空间直角坐标系原点位于参考椭球的中心o点，Z轴指向参考椭球的北极，X轴指向起始子午面与赤道的交点，Y轴位于赤道面上、且按右手系与X 轴正交于o点处目前，GPS测量中采用的WGS-84坐标就属于空间直角坐标系 某点在空间直角坐标系中的坐标，可用该点在此坐标系的各个坐标轴上的投影来表示 （2）空间大地坐标系 空间大地坐标是采用大地经纬度和大地高来描述空间位置的我国目前广泛采用的1954北京坐标系以及1980年西安坐标系都属于空间大地坐标系一类，但采用的椭球参数不同3）平面直角坐标系 平面直角坐标系是利用投影变换，将空间坐标通过某种数学变换映射到平面上。
投影变换的方法有很多，在我国采用的是高斯-克吕格投影（高斯投影） 二、高斯投影二、高斯投影 高斯投影又称为横轴等角切椭圆柱投影横轴等角切椭圆柱投影高斯投影是正形投影的一种，它除了满足正形投影的一般条件（长度比和方向无关）外，还应该满足高斯投影本身的特殊条件即必须满足以下3个条件： （1）中央子午线和地球赤道投影后成为相互垂直的直线，且为投影的对称轴； （2）中央子午线投影后长度不变； （3）投影具有正形条件，即等角投影 根据上述条件导出的高斯坐标正算高斯坐标正算[大地坐标(B，L)求平面坐标(x,y)] 式中，(x，y)为投影后的高斯平面纵、横坐标；X为经度为零时对应的纵坐标值，也就是赤道至纬度B处中央子午线弧长（一般采用积分的方法）；B为纬度；l’’为以秒为单位的经差；N为卯酉圈曲率半径； 编程实现流程描述：编程实现流程描述：数据的组织、界面设计定义各变量存储常数和中间计算结果结果输出 高斯坐标投影坐标反算高斯坐标投影坐标反算[由平面坐标(x,y)求大地坐标(B,L)]的公式为：式中,为横坐标值等于零时对应的纬度，也就是将x看做X时由子午线弧长公式反求出的纬度；Mf为横轴坐标值 等于零时所对应的子午圈曲率半径 ；其余下标f 的各量也都是类似上述的各自的相应的意义。
计算出经差后，即可根据中央线的经度，计算出经度Ｌ三、空间坐标转换三、空间坐标转换1 1 椭球面之间坐标转换椭球面之间坐标转换（1） 在相同的基准下，不同坐标系之间的转换，其中由大地坐标转换成空间直角坐标的公式为 由空间直角坐标转换成大地坐标的公式为：当Hi-Hi-1小于且Bi-Bi-1小于秒，迭代结束 2) 在不同的基准下，相同的坐标系之间的转换，实际上是基准间的转换基准间的转换方法很多，最常见的是布尔沙模型，又称为七参数转换模型七参数转换模型设两个空间直角坐标系间的七个转换参数分别是3个平移参数 2 2 椭球面与平面之间的坐标转换椭球面与平面之间的坐标转换(1)在相同的基准下进行的，如在前面所讲的高斯投影坐标正反算(2)在不同的基准下的平面坐标有时也可借助空间直角坐标系作为过渡坐标系完成不同系统间的转换例如，1954年北京坐标系和1980年国家大地坐标系内坐标间的转换 同一点在不同坐标系中的高斯平面坐标可通过以下过程来实现：式中，下标1和2分别代表两种坐标系统 当然,也可以不借助空间直角坐标系,而直接通过大地坐标系进行换算,其过程为：上述两过程，无论采用哪一种，除知道连个参心坐标系所属的地球椭球外，还必须知道(或求得)每一步转换的转换参数。
3 3 平面与平面之间的坐标转换平面与平面之间的坐标转换(1) 在相同的基准，相同的投影方式下进行的平面坐标间的转换发生在邻带之间的坐标换算例如，高斯投影是按一定的经差宽度分带的，由于各个子午线不同而形成的坐标系也不同，常常需要6度间转换，6度到3度以及3度之间进行邻带坐标换算其过程是按高斯投影反算 求得某点的大地坐标然后,按高斯投影正算公式求得该点在新的中央子午线为投影轴的邻带内的高斯平面坐标式中,x,y的下标1和2表示相邻的两个分带.(2)在相同的基准下,不同的投影方式下也可以产生不同的坐标系.例如，高斯投影平面坐标与墨卡托这类投影之间的坐标换算 (3) 不同基准下,平面与平面之间的转换除了按上面讲的式子进行外,还可以按如下过程来实现此外，常有这种情况,在一个测区内虽然是在相同的基准下，但各单位所做的控制要求不同，有时也需将地方独立控制网转换到国家网或其他别的新的控制网中，做到成果的相互利用和统一测区的坐标系统 例：例：WGS-84WGS-84坐标系到地方平面坐标系的转换流程坐标系到地方平面坐标系的转换流程 四、高程系统转换四、高程系统转换 （一）测量中的高程系统及关系（一）测量中的高程系统及关系 测量中常用的高程系统有大地高系统、正高系统以及正常高系统。
1）大地高系统 大地高系统：大地高系统：是以参考椭球为基准面的高程系统 大地高：大地高：某点的大地高是该点到通过该点的参考椭球的法线与椭球面交点的距离；可用H来表示 2）正高系统 正高系统正高系统 ：：是以大地水准面为基准面的高程系统 正高：正高：某点的正高是该点到通过该点的铅垂线与大地水准面的交点之间的距离；可用Hg表示3 3）正常高系统）正常高系统 正常高系统：正常高系统：是以似大地水准面为基准面的高程系统 正常高：正常高：某点的正常高是该点到通过该点的铅垂线与似大地水准面的交点之间的距离；可用Hr表示 地球表面大地水准面似大地水准面参考椭球面hgHgHHrξ图图2 2 高程系统之间的关系高程系统之间的关系 高程系统的转换关系高程系统的转换关系 大地水准面与参考椭球面之间的距离，称为大地水准面差异，记为hg大地高和正高之间的关系： H=Hg+ hg ………（1） 似大地水准面与参考椭球面之间的距离，称为高程异常，记为ξ 大地高系统和正常高系统之间的关系：H=Hr+ξ ………（2） （二）（二）GPSGPS高程高程 在一段时间内利用GPS来建立各类控制网时，绝大多数仅仅局限于解决平面坐标，高程仍沿用常规水准测量方法来测定 ，如何利用GPS观测中所提供的高程信息来直接为测绘服务就变成了一项很有意义的工作。
所谓高程拟合法就是利用范围不大的区域中，高程异常具有一定的几何相关性的原理，利用数学的方法，求解正高、正常高和高程异常 试想在一局部GPS网中，由若干个点的ξ作为已知值，用数值拟合方法内插出其他GPS测点的高程异常，按式(2)可求得各点的正常高 1.1.高程拟合算法高程拟合算法 高程拟合常有六种模型：多项式曲线拟合、三次样多条曲线拟合、Akima曲线拟合、多项式曲面拟合、多面函数法曲面拟合和移动法曲面拟合前三种属曲线拟合，仅当GPS点布设成测线时采用；后三种属于曲面拟合，当GPS测点分布设成网状时采用1)(1)多项式曲面拟合多项式曲面拟合 当GPS测点布成网状时，应用曲面拟合设测点的ξi和xi、yi存在如下函数关系：ξi=f(xi，yi)+εi (3) 式中，f(x,y)为趋势值；εi为误差 可选用以下空间曲面表达式：…… (4) 对于每一个已知点，都可列出以上方程，在Σε2=min条件下，解出bi，再按式(4) 求出待求点的高程异常ξ值 （（2 2）多面函数拟合模型）多面函数拟合模型 多面函数拟合的思想是在个数据点上建立一个曲面，通过将这些曲面按一定比例的叠加来最佳地描述所要求的物体表面，并使叠加后的曲面严格地通过各数据点。
多面函数拟合模型如式 上式中， 为模型参数， 为x，y的核函数， ( )为中心点，中心点为从公共点中选取高程异常显著的点，其个数应小于等于已知点数，m为核函数的个数，一般选择下面的正曲面作为函数 式（1） 若公共数为n（n≥m），则可构成误差方程式： 基于最小二乘法求得系数K：计算参数后，将待计算高程异常点计算参数后，将待计算高程异常点P(Xp,Yp)P(Xp,Yp)的坐标的坐标代入式（代入式（1 1），即可计算出该点所在位置的高程异常，），即可计算出该点所在位置的高程异常，即：即： 式（2）式（3）式（4）式（5） GPSGPS高程拟合程序设计高程拟合程序设计淮南市部分地区的GPS控制网数据一共有81个GPS点其中有26个水准联测点均匀的分布在测区的周围和中间 公共点点位分布图 GPS高程拟合模型（二次曲面）二次曲面拟合模型是把高程异常ζ近似地看作关于点坐标的函数，用已测得水准的GPS 点的高程异常拟合该函数,再用拟合好的函数模型按平面坐标的方法求解区域内任一点的高程异常值，进而求得该点的正常高。
而计算时我们只取二次项之前的项运算 1.1.原始数据文件组织原始数据文件组织2.2.用来存储数据文件中数据的变量、数据用来存储数据文件中数据的变量、数据的定义的定义u变量：变量：存储总点数tu主要数组：主要数组：行号d(),点号dm(),点的X坐标x() ,点的Y坐标y() ,点的大地高gpsg() , 点的正常高szg() , 系数矩阵B() , 常数项矩阵L() ,权阵p()3.精度评定4.数据输出 系数矩阵B()的生成（直接生成）For i = 0 To List2.ListCount - 1B(i + 1, 1) = 1B(i + 1, 2) = x1(i + 1) - vB(i + 1, 3) = y1(i + 1) - mB(i + 1, 4) = (x1(i + 1) - v) ^ 2B(i + 1, 5) = (y1(i + 1) - m) ^ 2B(i + 1, 6) = (x1(i + 1) - v) * (y1(i + 1) - m)Next L(i, 1) = szg (i) - gpsg (i)Next i系数矩阵L()的生成（直接生成） If i = n Then p(i, n) = 1 / Sqr((x1(i) - v) ^ 2 + (y1(i) - m) ^ 2) Else p(i, n) = 0 End If Next nNext i权阵P()的生成（直接生成） 精度评定精度评定------内符合精度和外符合精度内符合精度和外符合精度n内符合精度内符合精度根据参与计算的公共点来进行计算根据参与计算的公共点来进行计算n外符合精度外符合精度根据没有参与计算的公共点（检查点）来进行根据没有参与计算的公共点（检查点）来进行计算计算 上机内容（四）上机内容（四）水准网平差程序设计水准网平差程序设计1）输出高程平差值和高程中误差；2）输出观测值及其改正数与平差值；3）提交成果。