坐标系统转换和高程系统转换程序设计ppt课件

1、第五章第五章 坐标系统转换和高程系坐标系统转换和高程系统转换程序设计统转换程序设计 本章重点：本章重点：一一. .丈量中常用坐标系及这些坐标系的特丈量中常用坐标系及这些坐标系的特点；丈量中高程系统及其关系、点；丈量中高程系统及其关系、GPSGPS高程高程二二. .程序实现时的技巧程序实现时的技巧一、地理空间坐标系一、地理空间坐标系 为了确定物体在空间的位置，常用坐为了确定物体在空间的位置，常用坐标系来描画空间位置。丈量中常用的坐标标系来描画空间位置。丈量中常用的坐标系有以下几种：系有以下几种：1空间直角坐标系空间直角坐标系 空间直角坐标系原点位于参考椭球的中空间直角坐标系原点位于参考椭球的中心

2、心o点，点，Z轴指向参考椭球的北极，轴指向参考椭球的北极，X轴指轴指向起始子午面与赤道的交点，向起始子午面与赤道的交点，Y轴位于赤轴位于赤道面上、且按右手系与道面上、且按右手系与X 轴正交于轴正交于o点处。点处。目前，目前，GPS丈量中采用的丈量中采用的WGS-84坐标就坐标就属于空间直角坐标系。属于空间直角坐标系。 某点在空间直角坐标系中的坐标，可用该点在此坐标系的各个坐标轴上的投影来表示。2空间大地坐标系 空间大地坐标是采用大地经纬度和大地高来描画空间位置的。我国目前广泛采用的1954北京坐标系以及1980年西安坐标系都属于空间大地坐标系一类，但采用的椭球参数不同。3平面直角坐标系 平面直

3、角坐标系是利用投影变换，将空间坐标经过某种数学变换映射到平面上。投影变换的方法有很多，在我国采用的是高斯-克吕格投影高斯投影。二、高斯投影二、高斯投影 高斯投影又称为横轴等角切椭圆柱投高斯投影又称为横轴等角切椭圆柱投影。高斯投影是正形投影的一种，它除了影。高斯投影是正形投影的一种，它除了满足正形投影的普通条件长度比和方向满足正形投影的普通条件长度比和方向无关外，还应该满足高斯投影本身的特无关外，还应该满足高斯投影本身的特殊条件。即必需满足以下殊条件。即必需满足以下3 3个条件：个条件： 1 1中央子午线和地球赤道投影后成为中央子午线和地球赤道投影后成为相互垂直的直线，且为投影的对称轴；相互垂直

4、的直线，且为投影的对称轴； 2 2中央子午线投影后长度不变；中央子午线投影后长度不变； 3 3投影具有正形条件，即等角投影。投影具有正形条件，即等角投影。 根据上述条件导出的高斯坐标正算大地坐标(B，L)求平面坐标(x,y)6425644223422)5861(cossin720)495 (cossin24cossin2lttBBNltBBNl BBNXx 5222425532233)5814185(cos120)1 (cos6cosltttBNltBNlBNy 式中，(x，y)为投影后的高斯平面纵、横坐标；X为经度为零时对应的纵坐标值，也就是赤道至纬度B处中央子午线弧长普通采用积分的方法；B

5、为纬度；l为以秒为单位的经差；N为卯酉圈曲率半径； 206265 ;costan为地球第二偏心率，其中；eBeBt2222bbae22222sin1abaeBeaN，编程实现流程描画：编程实现流程描画：数据的组织、界面设计定义各变量存储常数和中间计算结果结果输出高斯坐标投影坐标反算高斯坐标投影坐标反算 由平面坐标由平面坐标(x,y)(x,y)求大求大地坐标地坐标(B,L)(B,L)的公式为：的公式为：64254223332)459061(720)935(242yttNMtyttNMtyNMtBBffffffffffffffff52224253223)8624285(cos1201)21 (co

6、s61cos1ytttBNytBNyBNlfffffffffffff式中,fB为横坐标值等于零时对应的纬度，也就是将x看做X时由子午线弧长公式反求出的纬度；Mf为横轴坐标值等于零时所对应的子午圈曲率半径 ；其他下标f 的各量也都是类似上述的各自的相应的意义。计算出经差后，即可根据中央线的经度，计算出经度。三、空间坐标转换三、空间坐标转换1 椭球面之间坐标转换椭球面之间坐标转换1 在一样的基准下，不同坐标系之间的在一样的基准下，不同坐标系之间的转换，其中由大地坐标转换成空间直角转换，其中由大地坐标转换成空间直角坐标的公式为坐标的公式为 BeaNabaebaBHbaNBHeNZLBHNYLBHNX

7、22222222sin1,sinsin)1 (coscos)(coscos)(曲率半径分别为；第一偏心率和卯酉圈为地球椭球长、短半轴其中，由空间直角坐标转换成大地坐标的公式为：222222)1 (sin)1 ()()(arctanarctanYXZEBeNBZHHeNYXHNZBXYL的初值可用下式求出，须用迭代的方法，在采用上式进行转换时的值。的初值再次求、求定的初值然后，利用纬度BNHEB当Hi-Hi-1小于0.001m且Bi-Bi-1小于0.00001秒，迭代终了。2) 在不同的基准下，一样的坐标系之间的转换，实践上是基准间的转换。基准间的转换方法很多，最常见的是布尔沙模型，又称为七参数

8、转换模型。设两个空间直角坐标系间的七个转换参数分别是3个平移参数TBBBTAAAzyxzyxAzyxmzy),(),(.1),(3)x(000中的坐标，系为某点在空间直角坐标令个尺度参数，个旋转参数，AAABBBzyxRmzyxzyx)()1 (0002 2 椭球面与平面之间的坐标转换椭球面与平面之间的坐标转换(1)(1)在一样的基准下进展的，如在前面所讲在一样的基准下进展的，如在前面所讲的高斯投影坐标正反算的高斯投影坐标正反算(2)(2)在不同的基准下的平面坐标有时也可借在不同的基准下的平面坐标有时也可借助空间直角坐标系作为过渡坐标系完成助空间直角坐标系作为过渡坐标系完成不同系统间的转换。例

9、如，不同系统间的转换。例如，19541954年北京年北京坐标系和坐标系和19801980年国家大地坐标系内坐标年国家大地坐标系内坐标间的转换。间的转换。 同一点在不同坐标系中的高斯平面坐标同一点在不同坐标系中的高斯平面坐标可经过以下过程来实现：可经过以下过程来实现：),(),(),(),(),(),(22222221111111yxLBZYXZYXLByx式中，下标1和2分别代表两种坐标系统。 当然,也可以不借助空间直角坐标系,而直接经过大地坐标系进展换算,其过程为：),(),(),(),(22221111yxLBLByx上述两过程，无论采用哪一种，除知道连个参心坐标系所属的地球椭球外，还必需

10、知道(或求得)每一步转换的转换参数。3 3 平面与平面之间的坐标转换平面与平面之间的坐标转换(1) 在一样的基准，一样的投影方式下进展的平面坐标间的转换发生在邻带之间的坐标换算。例如，高斯投影是按一定的经差宽度分带的，由于各个子午线不同而构成的坐标系也不同，经常需求6度间转换，6度到3度以及3度之间进展邻带坐标换算。其过程是按高斯投影反算求得某点的大地坐标),(),(1111yxLLyxBB然后,按高斯投影正算公式求得该点在新的中央子午线为投影轴的邻带内的高斯平面坐标),(),(),(),(111122111112yxLyxBFyyxLyxBFx式中,x,y的下标1和2表示相邻的两个分带.(2

11、)在一样的基准下,不同的投影方式下也可以产生不同的坐标系.例如，高斯投影平面坐标与墨卡托这类投影之间的坐标换算(3) 不同基准下,平面与平面之间的转换除了按上面讲的式子进展外,还可以按如下过程来实现),(),(),(),(221111111yxyxLBZYX此外，常有这种情况,在一个测区内虽然是在一样的基准下，但各单位所做的控制要求不同，有时也需将地方独立控制网转换到国家网或其他别的新的控制网中，做到成果的相互利用和一致测区的坐标系统。),(),(2211yxyx例：例：WGS-84WGS-84坐标系到地方平面坐标系的转换流程坐标系到地方平面坐标系的转换流程 四、高程系统转换四、高程系统转换

12、一丈量中的高程系统及关系 丈量中常用的高程系统有大地高系统、正高系统以及正常高系统。1大地高系统 大地高系统：是以参考椭球为基准面的高程系统 大地高：某点的大地高是该点到经过该点的参考椭球的法线与椭球面交点的间隔；可用H来表示。 2正高系统 正高系统 ：是以大地水准面为基准面的高程系统。 正高：某点的正高是该点到经过该点的铅垂线与大地水准面的交点之间的间隔；可用Hg表示。3正常高系统 正常高系统：是以似大地水准面为基准面的高程系统 正常高：某点的正常高是该点到经过该点的铅垂线与似大地水准面的交点之间的间隔；可用Hr表示。地球外表大地水准面似大地水准面参考椭球面hgHgHHr图图2 2 高程系统

13、之间的关系高程系统之间的关系高程系统的转换关系高程系统的转换关系 大地水准面与参考椭球面之间的间隔，称为大地水准面与参考椭球面之间的间隔，称为大地水准面差别，记为大地水准面差别，记为hghg。大地高和正高之间的关系：大地高和正高之间的关系： H=Hg+ hg H=Hg+ hg 1 1 似大地水准面与参考椭球面之间的间隔，称似大地水准面与参考椭球面之间的间隔，称为高程异常，记为为高程异常，记为 。大地高系统和正常高系统之间的关系：大地高系统和正常高系统之间的关系：H=Hr+ H=Hr+ 2 2二二GPSGPS高程高程 在一段时间内利用GPS来建立各类控制网时，绝大多数仅仅局限于处理平面坐标，高程

14、仍沿用常规水准丈量方法来测定 ，如何利用GPS观测中所提供的高程信息来直接为测绘效力就变成了一项很有意义的任务。 所谓高程拟合法就是利用范围不大的区域中，高程异常具有一定的几何相关性的原理，利用数学的方法，求解正高、正常高和高程异常。 试想在一部分GPS网中，由假设干个点的作为知值，用数值拟合方法内插出其他GPS测点的高程异常，按式(2)可求得各点的正常高。 1.1.高程拟合算法高程拟合算法 高程拟合常有六种模型：多项式曲线拟合、高程拟合常有六种模型：多项式曲线拟合、三次样多条曲线拟合、三次样多条曲线拟合、AkimaAkima曲线拟合、多项式曲线拟合、多项式曲面拟合、多面函数法曲面拟合和挪动法

15、曲面拟曲面拟合、多面函数法曲面拟合和挪动法曲面拟合。前三种属曲线拟合，仅当合。前三种属曲线拟合，仅当GPSGPS点布设成测线点布设成测线时采用；后三种属于曲面拟合，当时采用；后三种属于曲面拟合，当GPSGPS测点分布测点分布设成网状时采用。设成网状时采用。(1)(1)多项式曲面拟合多项式曲面拟合 当当GPSGPS测点布成网状时，运用曲面拟合。设测点布成网状时，运用曲面拟合。设测点的测点的ii和和xixi、yiyi存在如下函数关系：存在如下函数关系：i=f(xii=f(xi，yi)+iyi)+i(3)(3) 式中，式中，f(x,y)f(x,y)为趋势值；为趋势值；ii为误差。为误差。可选用以下空

16、间曲面表达式： (4)对于每一个知点，都可列出以上方程，在2=min条件下，解出bi，再按式(4) 求出待求点的高程异常值。 2 2多面函数拟合模型多面函数拟合模型 多面函数拟合的思想是在个数据点上建立一个曲面，经过将这些曲面按一定比例的叠加来最正确地描画所要求的物体外表，并使叠加后的曲面严厉地经过各数据点。多面函数拟合模型如式 miiiiyxyxQK1,上式中， 为模型参数， 为x，y的核函数， ( )为中心点，中心点为从公共点中选取高程异常显著的点，其个数应小于等于知点数，m为核函数的个数，普通选择下面的正曲面作为函数。 ), 2 , 1(miKiiiyx ,iiyxyxQ,式1假设公共数

17、为nnm，那么可构成误差方程式： 基于最小二乘法求得系数K：21222,iiiiyyxxyxyxQ111nmmnnKQVTTQQQK1KQQQPnPPP21计算参数后，将待计算高程异常点计算参数后，将待计算高程异常点P(Xp,Yp)P(Xp,Yp)的坐标的坐标代入式代入式1 1，即可计算出该点所在位置的高程异常，即可计算出该点所在位置的高程异常，即：即： 式2式3式4式5GPSGPS高程拟合程序设计高程拟合程序设计淮南市部分地域的GPS控制网数据一共有81个GPS点其中有26个水准联测点均匀的分布在测区的周围和中间。 公共点点位分布图GPS高程拟合模型二次曲面二次曲面拟合模型是把高程异常近似地

18、看作关于点坐标的函数，用已测得水准的GPS 点的高程异常拟合该函数,再用拟合好的函数模型按平面坐标的方法求解区域内任一点的高程异常值，进而求得该点的正常高。而计算时我们只取二次项之前的项运算。 yxf,xyayaxayaxaayxf52423210,1.1.原始数据文件组织原始数据文件组织2.2.用来存储数据文件中数据的变量、数据用来存储数据文件中数据的变量、数据的定义的定义变量：存储总点数变量：存储总点数t t主要数组：行号主要数组：行号d(),d(),点号点号dm(),dm(),点的点的X X坐标坐标x() ,x() ,点的点的Y Y坐标坐标y() ,y() ,点的大地高点的大地高gpsg() , gpsg() , 点的正常高点的正常高szg() , szg() , 系数矩阵系数矩阵B() , B() , 常数常数项矩阵项矩阵L() ,L() ,权阵权阵p()p()3.3.精度评定精度评定4.4.数据输出数据输出系数矩阵B()的生成直接生成For i = 0 To List2.ListCount -