您好,欢迎来到爱站旅游。
搜索
您的当前位置:首页坐标转换的方法及应用

坐标转换的方法及应用

来源:爱站旅游
第37卷第5期 现代测绘 V01.37。No.5 2014年O9月 Modern Surveying and Mapping Sep.2014 坐标转换的方法及应用 鲍建宽 (黑龙江工程学院测绘工程学院,黑龙江哈尔滨150050) 摘要本文介绍了几种主要的坐标转换模型,综合分析各模型的特点和适用范围,针对坐标转换模型的选用、转 换参数的解算、转换计算的方法、转换计算中值得注意的问题加以研究和探讨。结论内容,对于在测量实践中进行 局部坐标系和地心坐标系的转换、各局部坐标系间的转换具有参考价值。 关键词坐标转换相似转换多项式逼近公共点布尔莎转换模型 中图分类号:P226 文献标识码:A 文章编号:1672—4097(2014)05—0OO3—06 为确定地面点或空间目标的位置而建立的大 顾及两个坐标系尺度不尽一致,从而还有一个尺度 地坐标系是通过参考椭球的定位和定向实现的。 变化参数( ),共计有7个参数。当旋转角e较小时, 随着测量技术的发展,在不同时期建立的大地坐标 相应的布尔莎坐标转换公式为: 系采用的参考椭球参数及定位定向的方式是不相 同的,因此,就存在不同大地坐标系统的坐标互相 转换问题。目前,我国的测绘成果资料中就有1954 I 三 ] 北京坐标系、新1954北京坐标系、1980年国家大地 坐标系、2000中国大地坐标系几种坐标系共存。还 1.1.2布尔莎模型的转换方法 有GPS系统采用的WGS一84坐标系、GLONASS 为实现两个坐标系的坐标转换,在实践中,通 系统采用的PZ一9O坐标系,这些卫星定位测量成 过将两坐标系中公共已知点的坐标代入(1.1)式, 果也需要转换为我国的相应坐标系中。此外,在一 解算转换参数,然后再由(1.1)式把原坐标系内的 些地区,由于工程建设的急需(如矿山开发、新城市 坐标(X , , )转换为新坐标(X , , ),最后 建设等)而建立的地方坐标系有时也需要把独 把空间直角坐标(Xi,Yi,Zi)化为大地坐标(B ,L , 立坐标系的坐标转换到相应国家统一的坐标系中 H ),再化为测绘工程中使用的高斯平面坐标。 去。本文就不同的坐标系转换模型的方法问题进 为了可靠地求定两个不同空间直角坐标系的 行探讨。 坐标转换公式中的7个转换参数X。,y0,Z0,ex,e'Y, 1基于三维直角坐标的坐标转换模型 e ,/1,则至少需要3个公共点,当具有3个以上的公 共点时,可按最小二乘法求得7个参数的最或然值。 1.1布尔莎(Bursa)转换模型 若令al=1+ ,a2=以lex,a3=a1ey,a4=:=口1ez,则可 两空间直角坐标系的关系如图1所示,其中坐 将(1)式改写为: 标系o,一X y,Z 是由原坐标系0一XYZ通过绕坐 △ 标轴旋转(3个旋转参数分别为ex,e ,ez)、坐标原点  .平移(3个平移参数分别为△ ,△lyO,△Zo)得到,再 冈一J  △ n1 122 123 口4 (2) 若取 VyY'i, 图1 两个三维直角坐标的关系 值 4 由(1.2)式得误差方程 [-1 0 0 X 0 现代测绘 第37卷 一Zi Yi-I l_ o0 1  o o 1 z o —Xl l 一 x o j 位位移将不可避免。为削弱误差较大甚至存在粗 差的公共点对解算转换参数的影响,可采用稳健估 计的方法进行求解以确保转换参数计算的可靠性。 (4)当利用3个以上的公共点求解转换参数时 存在多余观测,由于公共点误差的影响而使得转换 的公共点的坐标值与已知值不完全相同,而实际工 作中又往往要求所有的已知点的坐标值保持固定 ● ● ● ● ● △Xn . △ al 乙 ● (4) 不变。为了解决这一矛盾,可采用配置法,将公共 K K ¨ a2 点的转换值改正为已知值,对非公共点的转换值进 。 n3 行相应的配置。具体方法是:先计算公共点坐标转 P 换值的改正数 一已知值 一转换值 (i一1、2、…、 (24 ~P 。 P —P 当有 (≥3)个公共点时,可分别列出3 Trt个误 差方程。如果平移参数以米为单位;尺度比参数的 单位取用10—6;旋转角参数以秒为单位时,误差方 程系数矩阵B应表示为如下形式: B 根据最小二乘法VTPV--rain的原则,可由误 差方程列出法方程,进一步解算出转换参数。 1.1.3 实用中应注意的问题 (1)布尔莎模型是舍去高次项后的线性化模 型,当旋转角e≤5 时,对计算转换坐标的影响值最 大不过5 mm;当旋转角£≤10”时,对计算转换坐标 的影响值最大不过18 mm;当旋转角£≤100”时,对 计算转换坐标的影响最大可达1.8 1TI。可见,使用 该模型时,坐标轴旋转角宜小于1O”。 (2)所求得转换参数的精度与公共点坐标的精 度和误差方程式系数阵B有关,而误差方程式的系 数阵B又取决于公共点的个数和几何分布情况,即 转换参数的精度又与公共点的个数与几何分布情 况有关。因此,为了求得较好的转换参数,应选择 一定数量的精度较高且分布较均匀并有较大覆盖 面的公共点。 (3)由于1954北京坐标系统的大地点坐标是 经过局部分区平差得到的,这导致整个1954北京坐 标系统的内符合性非常差。另外,随时间的推移、 自然条件的变化、地球物理方面的变迁等,局部点 ~优);P 。再计算非公共点转换值的改正数。 蜀 P; ~P 。 一P  P墨一P 。 P 一P— 。J ‘i7J  一上∑Pj ∑  f一1 式中irft为公共点的个数; 一1/S 为权;S 为非公共点 与公共点i间的距离。 (5)当旋转角e>10 时,(1)式中的旋转矩阵为 CO SCvCOSEz coseyslr ̄z  Iz co踺x∞sEz+SII ̄XSlrlEySlnEz l sinexC0S ̄2+COSE"xsinevsin ̄z_d 可采用导数的方法、高斯一牛顿法或者无约束 优化中的Rosenbrock方法等等计算转换参数。 1.2莫罗金斯基(Molodensky)转换模型 1.2.1莫罗金斯基模型公式 布尔莎坐标转换公式(1.1)是对坐标点在原坐 标系 一XyZ中的位置向量施加旋转和尺度转换, 旋转中心就在坐标原点o,尺度变换作用于以原点 0为起点的向径OP 。但在大量的实际应用中,往 往只是在局部地区进行坐标系的转换,在已有的局 部地面控制网中,由于所有的控制点坐标都是由网 中某一参考点(例如国家大地控制网中的大地原 点、首级工程控制网的起始点)依次推算得出的。 图2地方坐标系与目标坐标系的关系 第5期 鲍建宽:坐标转换的方法及应用 5 因此就应将旋转和尺度变换的作用范围在控制 网区域以内,宜以该参考点PK作为变换中心(图2), 与此相应的即为莫洛金斯基坐标转换模型。 网参考点PK的站心地平坐标系的三个坐标轴 P z,P ,P 2为旋转轴,这点是与莫洛坚斯基转 换模型所不同的;而(AXo,△y0,△ )为坐标系O— XYZ的原点O在坐标系0 一X y Z 中位置向 量,亦即两坐标系间的三个平移参数; 为参考点 PK到网中任一点P 间向径的尺度变化参数。这构 1 因 ㈣ 成了范士转换模型的7个坐标转换参数。 在莫洛金斯基坐标转换模型公式(1.5)中, (△x0,△y0,△Zo)为坐标系0一XyZ的原点0在坐 标系0 一X y Z 中位置向量,亦即两坐标系间的 三个平移参数;(XK,YK,ZK)则为参考点P 在O— XYZ坐标系中的位置参数;(£x,£y,£z)是以参考点 PK为旋转中心的旋转参数; 为参考点PK到网中 任一点P 间向径的尺度变化参数。 1.2.2 实用中应注意的问题 (1)莫洛金斯基坐标转换模型中转换参数的解 算方法同布尔莎转换模型。 (2)莫洛金斯基坐标转换模型更适用于参心坐 标系与地心坐标系之间的转换,以及地方坐标系与 国家坐标系之间的转换。 (3)在进行参心坐标系与地心坐标系之间的转 换时,由于参心大地坐标系都是用传统测量技术建立 的,并非严格意义上的三维坐标系,难以实现与三维 地心坐标系的严格转换,求得的参数就不够可靠。又 因为参心大地坐标系在参考椭球定位定向时难以满 足理想条件,坐标系本身的偏差,也存在着地面网定 位定向的系统误差,且由于地壳形变等多种原因使地 面网点的坐标值随时间延续而发生变化。因此,地面 网参心坐标系的系统误差,使得数值较小的三个旋转 角难以准确求定,旋转角参数的估值在数值上往往与 其本身的中误差相接近。这样,有时就不考虑旋转和 尺度变换,而仅取用三个平移参数。 1.3武测转换模型 在莫洛金斯基转换模型中,若把旋转中心由参 考点PK变为坐标原点o,亦即旋转参数(ex,ey,ez) 仍为坐标系0一XYZ通过以为0中心绕坐标轴旋 转的角度,则得到武测转换模型的转换公式: I y, I=I△ I+l yK I+(1+ )l Yi—yK l+ z 圄 ㈣ 1.4范士转换模型 如图3所示,三个旋转角(叫 ,∽,叫 )是以控制 据(sx,£y,£z)与( x,∽,叫z)的关系 LE J Lin 一B。K。c OBsL oK ÷inL 一K二 旧  i B J j 把上式代人(5)式既得范士转换模型的转换 公式。 1.5广义大地坐标微分公式 1.5.1基本公式 对于不同大地坐标系的换算,除了包含3个平 移参数、3个旋转参数和1个尺度变化参数外,还包 括2个地球椭球元素变化参数。不同大地坐标系的 换算公式为: 豳一 A …[ ]㈦ 式中: fl ~(N+H)S1I1L cosD ̄。(N+H)COsL cosD ̄ ^  lA。一lI 一 再 。sinBcosL  一 而sir出sinL P M ̄-cosB.o - ll  l cosBc0sL cosBsinL sinB J r anBc。sL nB。inL 一 ] A—ll 札 ∞ o I, Ne2 sinBcosBsinL—T——TNe2 sinBcosBcosL o j 一 si衄c。s l LN(]一 i zB)+Hl_ 6 现代测. 绘 第37卷 嚣0M((+2一HP) (si1n一。B 口) ……。叶  I(1 si㈣sin B j 式(8)通常称为广义大地坐标微分公式或广义 变换椭球微分公式。 1.5.2 实用中应注意的问题 (2)与三维坐标转换一样,由于公共点的坐标存 在误差,求得的转换参数将受其影响,公共点坐标误差 对转换参数的影响与点位的几何分布及点数的多少有 关,因而为了较好地求得转换参数,应选择一定数量的 精度较高且分布较均匀并有较大覆盖面的公共点。同 样,可采用稳健估计的方法进行求解,以此削弱误差较 大甚至存在粗差的公共点对解算转换参数的影响。 (1)根据 (>3)个以上公共点的两套大地坐 标值,由(8)式参照(3)和(4)式可以列出3 个误差 方程,采用最小二乘原理可求出其中的转换参数。 (2)在通常隋况下,两个不同大地坐标系的 、如 是已知的。对于同一国家或地区,尽管大地坐标系因 参考椭球元素和定位定向改变,但坐标轴总设定是平 行的。对于同一大地网,长度变化参数IX总为零,对于 不同大地网, 相差也不会太大。因此,在一般情况下 可以略去科、 诸参数,而只解求XO, , 三个 图4两个平面直角坐标系的关系 椭球定位变动量,即就是坐标系原点的平移量。 2基于平面直角坐标的坐标转换模型 2.1相似变换模型 2.1.1基本公式 (3)相似转换仅适用于投影方式、椭球元素、椭 球定位均相同的情形。其优点是转换参数的意义 较为明确,计算简便,其常数能保证转换过程的正 形性。不足之处在于,它只适用于小范围内、局部 如图4所示,旧坐标与新坐标的关系为: [ 一 ][ [ 坐标系间的坐标转换,且没有考虑两坐标系问的局 部变形和误差积累,转换精度较低。 (4)转换公式的转换元素 、6、0、k中,口、6、0是 常量,而k由于存在边长高斯改正,则有所变动。坐 若取c=kcos ̄,d--ksinO,上式写为 [ 一[ 三 ][ ] 此即为相似转换模型。各参数的几何意义:尼 为尺度参数; 为旋转参数;n、b为平移参数。 2.1.2 实用中应注意的问题 标转换时,在一定范围内,可以将测区投影高程面 及高斯平面各边的长度比k视为一个常量,其原则 是在一个测区内用一个k值进行转换,其长度变形 值不超过25 cm/km。由于长度变形与线段距X轴 的距离(即端点y坐标的平均值 )、线段长度及该 地的平均地球曲率半径有关,据此可计算出相似转 (1)当测区内具有新旧两种坐标系的 (>2)个 重合点时,即可利用最小二乘法求定新旧坐标系的转 换参数。误差方程的列立参照参照(3)和(4)式。 换模型的适用区域,结果列于表1。 单位:km 1OO 11O 12O 13O 表1 测区范围与测区中心所处位置关系 40 50 60 70 80 90 两坐标平行时测区范围Ay 两坐标垂直时测区范围Ay 25.4 2O.3 2O.3 17.3 16.9 15.0 14.5 13.2 12.7 11.8 11.3 1O.O 1O.2 9.7 9.2 8.9 8.5 8.2 7.8 7.6 2.2正形转换模型 (10) 2.2.1基本公式 2.2.2转换方法及应注意的问题 依据正形投影的柯西一黎曼条件,可推导出实 现两坐标系间坐标转换的正形转换模型 X 一X+p。+却 一Yqt+(X2一yz)pz一2XYqz+( 一3X ̄)Ps] (1)在式(2.2)中仅取一次项,则与式(2.1)表 示的相似变换模型一样。与相似变换相比,正形变 换法比相似变换法精度要高,它适用于变换区域相 y— +‘ 一。 + P4一‘ y一 q4+…l 当大的情形。 3_+(  一3 )q、o 4- yqpt -'}-Xq ̄+2 ‘ +(4X3Y一4 )P4+( 一6 。 、-Y3 )p3 J + )吼+…J (2)同样,正形转换也要求两坐标系的投影方 、……………一………………… 第5期 鲍建宽:坐标转换的方法及应用 7 X 一ao+a1X+a2Y+口3 +a4XY+a5y2] Y,+n6一6o +6+n7 x+6 y+a82y+6x +a9。 +64XY+ }y3+… l  1) +b6 +b7 Y+b8x +b9y3+… J 一X+bo+61(x—Xo)+b2(Y—Yo)+ (x—x。) 定系数到a 和b 至少需要6个公共点坐标,当公 共点多于6个时,利用最t] ̄--乘法解算待定系数。 (2)多项式拟合法是一种采用平面坐标位置作 为因变量的二维坐标转换方法。当公共点分布较 均匀且点间距离适当的情况下,所建立数学模型能 够达到较好转换精度。 (3)多项式拟合法的最大特点是以转换点的高 斯平面坐标为因变量,不需知道目标椭球似大地水 准面差距 ,就能够实现不同坐标系问的二维坐标 转换,具有坐标转换精度高、成果可靠、便于应用等 优点,适合任意坐标系间平面坐标转换计算。 (4)如果需要转换的区域较大,公共点较多,可 以选择更高阶数的多项式进行转换。 3结束语 在多种测量坐标系统并存共用的今天,坐标转 换计算工作是测量计算的重要内内容之一。本文 介绍了几种主要的坐标转换模型,并分析了它们间 的联系、各自的意义和适用范围。综合前面的论 述,对坐标转换计算问题做一小结,供参考探讨。 (1)在选用坐标转换模型时,要综合考虑待转 换的两个坐标系的基本信息、待转区域的地理位置 (距离央子午线的距离)、区域的面积、公共点的多 少与分布等情况,来选择适宜的坐标转换模型。 (2)对公共点的选择,精度上要择优选控制网 的起算点及高精度控制点;其个数要多于模型要求 的最少个数,当然要顾及转换区域的面积与区域形 状,一般的,采用二维转换模式至少选取2个以上的 公共点,采用三维转换模式至少选取3个以上的重 合点,大中城市至少选取5个公共点,小城市至少选 取4个重合点;对公共点位置的要求,公共点要要分 布均匀,包围城市区域,大城市内部还要有一公共 点点;在城市内部选定若干个均匀分布的公共点对 坐标转换精度进行检核。 (3)对于不同坐标系之间的坐标转换,目前最 广的是平面四参数转换模型与布尔莎七参数坐标 转换模型两种。平面四参数转换模型原理简单,数 值稳定可靠;对较小区域,它转换的精度较高,但当 范围较大时,由于受投影变形误差的影响,其转换 精度就较差,因而它只适合于较小区域的坐标转 换。布尔莎七参数转换模型为三维模型,在空间直 角坐标系中,两坐标系之间存在严密的转换模型, 由于理论比较严密,不存在模型误差和投影变形误 差,因而它适合于任何区域的坐标转换。 【下转第28页】 28 方便。 参考文献 1 现代测绘 社,2004. 第37卷 4邬伦,刘瑜,张晶,等.地理信息系统原理方法和应用 [M].科学出版社.2001 5冯克忠,姜遵锋,徐杨,等.ArcObjects开发指南[M].北 京:电子工业出版社,2007. Peng Z R,Huang R Design and development of interactive trip planning for web-based transit information systems[J].Transportation Research Part C:Emerging Technologies,2000,8(1):409—425. 6何涛,张世禄.基于ArcGIS的县级林业资源管理信息系 统研究I-J],2009,19(2):183—186. 7 ESRI.ArcGISDeveloper Help[S]. ESRIProduct Document,2004. 、 2张存保,李华,严新平.基于Web GIS的城市公交问路系 统I-J].武汉理工大学学报(交通科学与工程版),2004, 28(1):99—102. 8杨利兵,陈艳红,刘亚立.基于Geodatabase的地理数据 3黄杏元.地理信息系统概论EM].北京:高等教育出版 库的设计[J].山西建筑,2007,33(25):367—368. The Research and Design Bus Route of Points of Interest Query System based on GIS GU Peng—fei,GAO Cheng-fa (School of Transportation,Southeast University,Nanjing Jiangsu 210096,China) Abstract In order to use the public transportation more efficiently and alleviate traffic congestion in the city,this paper presents a new method of travel route query in the field of public transport.It can provide bus route of points of interest query based on the location of bus stop without having to transfer.According to the nearest bus—stop information around the location of the user,it can find out points of interest which meet the certain conditions and give user the bus routes to go to the destination.It set up a corresponding spatial database to store the information such as bus stops.Based on ArcEngine development platform,using C#.net language,it development the Point of Interest Around Bus Stops Query System.By entering the point of interest information and bus stop directly,the user could obtain a path of trave1.In this paper,the system takes Shenzhen as an example for testing and the result can meet the demand,so it can provide a new public transportation route query platfotin for people’s trave1. Key words GIS;routequery;ArcEngine;bus application 【上接第7页】 (4)对于各局部坐标系间的转换,当投影方式 等不同时,可采用多项式转换模型。当投影方式相 大,精度要求较高,或基准点不匹配时,可采用稳健估 计和最小二乘配置法提高转换精度。也可以采用综 合法进行坐标转换,就是在相似变换(布尔莎七参数 转换)的基础上再对空间直角坐标残差进行多项式拟 同时,二维坐标转换模式通常有平面四参数转换模 型、二维七参数转换模型,三维坐标转换模式通常 有布尔莎七参数转换模型、三维七参数转换模型; 如果变换区域范围较大,精度要求较高,可采用正 形转换。 (5)由局部坐标系向地心坐标系的转换时,多项 式转换模型是一种切实可行的方法。当转换区域较 合,系统误差通过多项式系数得到消弱,使统一后的 坐标系框架点坐标具有较好的一致眭,从而提高坐标 转换精度。或者在相似变换(布尔莎七参数转换)的 基础上,再利用正形转换模型在两坐标系间的平面 坐标间进行转换。 Methods of Coordinate Transformation and Its Application BAO Jian.kuan (Heilongjiang Institute of Technology,Harbin Heilongjiang 150050,China) Abstract This paper gives all kinds of coordinate transformation models,synthesized analysis of the characteristics of each model and applicable scope&About basic model selection of coordinate transforrnation,calculation of transformation parameters,calculation method of transforrnation and problems existing in transformation calculation are researched and discussed.Conclusions content will have a reference value for practical transformation from local geodetic coordinate syste ̄n to geocentric geodetic coordiate system annd from loca1 geodetic coordinate system to 1oea1. Key words coordinate transformation;similarity transformation;polynomial approximation;common point;Bursa model 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- azee.cn 版权所有 赣ICP备2024042794号-5

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务