第五章 坐标系统转换与高程系统转换程序设计

第五章坐标系统转换和高程系统转换程序设计

本章重点：一.测量中常用坐标系及这些坐标系的特点；测量中高程系统及其关系、GPS高程二.程序实现时的技巧一、地理空间坐标系

为了确定物体在空间的位置，常用坐标系来描述空间位置。测量中常用的坐标系有以下几种：（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）在相同的基准下,不同坐标系之间的转换，其中由大地坐标转换成空间直角坐标的公式为由空间直角坐标转换成大地坐标的公式为:当Hi-Hi-1小于0.001m且Bi-Bi-1小于0.00001秒，迭代结束。2)在不同的基准下，相同的坐标系之间的转换，实际上是基准间的转换。基准间的转换方法很多，最常见的是布尔沙模型，又称为七参数转换模型。设两个空间直角坐标系间的七个转换参数分别是3个平移参数2椭球面与平面之间的坐标转换(1)在相同的基准下进行的，如在前面所讲的高斯投影坐标正反算(2)在不同的基准下的平面坐标有时也可借助空间直角坐标系作为过渡坐标系完成不同系统间的转换。例如，1954年北京坐标系和1980年国家大地坐标系内坐标间的转换。

同一点在不同坐标系中的高斯平面坐标可通过以下过程来实现：式中，下标1和2分别代表两种坐标系统。

当然,也可以不借助空间直角坐标系,而直接通过大地坐标系进行换算,其过程为：上述两过程，无论采用哪一种，除知道连个参心坐标系所属的地球椭球外，还必须知道(或求得)每一步转换的转换参数。3平面与平面之间的坐标转换(1)在相同的基准，相同的投影方式下进行的平面坐标间的转换发生在邻带之间的坐标换算。例如，高斯投影是按一定的经差宽度分带的，由于各个子午线不同而形成的坐标系也不同，常常需要6度间转换，6度到3度以及3度之间进行邻带坐标换算。其过程是按高斯投影反算求得某点的大地坐标然后,按高斯投影正算公式求得该点在新的中央子午线为投影轴的邻带内的高斯平面坐标式中,x,y的下标1和2表示相邻的两个分带.(2)在相同的基准下,不同的投影方式下也可以产生不同的坐标系.例如，高斯投影平面坐标与墨卡托这类投影之间的坐标换算(3)不同基准下,平面与平面之间的转换除了按上面讲的式子进行外,还可以按如下过程来实现此外，常有这种情况,在一个测区内虽然是在相同的基准下，但各单位所做的控制要求不同，有时也需将地方独立控制网转换到国家网或其他别的新的控制网中，做到成果的相互利用和统一测区的坐标系统。四、高程系统转换

（一）测量中的高程系统及关系测量中常用的高程系统有大地高系统、正高系统以及正常高系统。1）大地高系统

大地高系统：是以参考椭球为基准面的高程系统

大地高：某点的大地高是该点到通过该点的参考椭球的法线与椭球面交点的距离；可用H来表示。

2）正高系统

正高系统

：是以大地水准面为基准面的高程系统。

正高：某点的正高是该点到通过该点的铅垂线与大地水准面的交点之间的距离；可用Hg表示。3）正常高系统

正常高系统：是以似大地水准面为基准面的高程系统

正常高：某点的正常高是该点到通过该点的铅垂线与似大地水准面的交点之间的距离；可用Hr表示。地球表面大地水准面似大地水准面参考椭球面hgHgHHrξ图2高程系统之间的关系高程系统的转换关系

大地水准面与参考椭球面之间的距离，称为大地水准面差异，记为hg。大地高和正高之间的关系：H=Hg+hg………（1）

似大地水准面与参考椭球面之间的距离，称为高程异常，记为ξ。大地高系统和正常高系统之间的关系：H=Hr+ξ………（2）（二）GPS高程在一段时间内利用GPS来建立各类控制网时，绝大多数仅仅局限于解决平面坐标，高程仍沿用常规水准测量方法来测定，如何利用GPS观测中所提供的高程信息来直接为测绘服务就变成了一项很有意义的工作。所谓高程拟合法就是利用范围不大的区域中，高程异常具有一定的几何相关性的原理，利用数学的方法，求解正高、正常高和高程异常。

试想在一局部GPS网中，由若干个点的ξ作为已知值，用数值拟合方法内插出其他GPS测点的高程异常，按式(2)可求得各点的正常高。1）高程拟合算法高程拟合常有六种模型：多项式曲线拟合、三次样多条曲线拟合、Akima曲线拟合、多项式曲面拟合、多面函数法曲面拟合和移动法曲面拟合。前三种属曲线拟合，仅当GPS点布设成测线时采用；后三种属于曲面拟合，当GPS测点分布设成网状时采用。多项式曲面拟合

当GPS测点布成网状时，应用曲面拟合。设测点的ξi和xi、yi存在如下函数关系：ξi=f(xi，yi)+εi

(3)式中，f(x,y)为趋势值；εi为误差。可选用以下空间曲面表达式：