GIS坐标系的分类.docVIP

GIS坐标系的分类 对于坐标系的概念，相信从事GIS行业的人都不陌生，曾经有段时间我对它是一知半解的，头痛了很久，现在有空就做个总结，希望对以后的学习会有所帮助。文章主要参阅了青岛海洋地质研究所，戴勤奋老师的《坐标转换公式》译文，另外还参考了同济大学，王解元、王军、陆彩萍的《WGS84与北京54的转换问题》，以及网络上发布的测量学的基本知识等文章。 2 基础知识 2.1大地椭球体 同大地体相近的可以用数学方法来表达的旋转椭球体。（大地体是由大地水准面包围而成。大地水准面是假定在重力作用下海水面静止时的平均水面，并设想此面穿过大陆与岛屿，连续扩展形成处处与铅垂线成正交的闭合曲面。） 2.2基准面 大地基准面是利用特定椭球体对特定地区地球表面的逼近，用以对椭球体进行定位的。 2.3投影 指按照一定的数学法则将地球面上的经纬网转换到平面上，使地面的地理坐标与平面直角坐标建立起函数关系。实现由曲面向平面的转化。常用的投影大概有二三十种，投影的选取要考虑地图的用途，投影的形变大小等众多因素。 2.4大地椭球体和基准面的关系 基准面是在椭球体的基础上建立的，两者不能等同，同一个椭球体可以定义不同的基准面，两者是一对多的关系。 3.坐标系的分类 3.1地理坐标系 大地坐标系：大地坐标系是大地测量中以参考椭球面为基准面建立起来的坐标系。也叫经纬坐标系，地面点的位置用大地经度（L）、大地纬度(B)和大地高度(H)表示。 地心坐标系：用X、Y、Z表示的直角坐标系，其中，直角坐标系原点位于地心；Z轴为极轴，向北为正；X 轴穿过本初子午线与赤道的交点；Y轴穿过赤道与东经90°的交点。 3.2投影坐标系 一种平面坐标系，依赖于地理坐标系和投影方式两个条件。 4.坐标系的转换 4.1大地坐标与地心坐标的互换 其实这两种坐标系都是描述地理坐标的，只是描述的方式不一样，其中大地坐标系的表述应该更常见些，依据不同的表示格式要求，可以对这两种表示方式进行相互转换。 在进行转换之前先引入一些变量： L 经度 B 纬度 H 相对于椭球体的高度，非重力高度 X 地心坐标X Y 地心坐标Y Z 地心坐标Z a 椭球体长半轴 b 椭球体短半轴 f 椭球体扁率 e 椭球体第一偏心律 e′椭球体第二偏心律 N 纬度B处的卯酉圈曲率半径。 约定：大地坐标系中的零度经线为格林威治子午线。如果定义不一致，需将本文中的各个公式中的零度经纬线转换到格林威治子午线，公式中的角度单位是弧度。 大地坐标转地心坐标公式： X=(N+H)cosBcosL Y=(N+H)cosBsinL Z=[(1*e2)N+H]sinB 地心坐标转大地坐标公式： B =tg-1[(Z+e2NsinB)-(X2+Y2)0.5] L =tg-1（Y/X） H =XsecLsecB-N 4.2不同大地坐标系的互换 这种变换俗称基准面变换，可以分为如下三步：A大地坐标系—A地心坐标系—B地心坐标系—B大地坐标系。其中第二步转换是关键。对于这类转换，前人已有深入的研究，现列出常用的三参数和七参数的转换方法。 三参数法 Xt=Xf+dX Yt=Yf+dY Zt=Zf+dZ 其中Xt、Yt、Zt表示新坐标坐标值, Xf、Yf、Zf表示原来的坐标系坐标值，dX、dY、dZ则表示新坐标系相对于就坐标系的偏移量。该方法假定两个大地坐标系相对应的地心坐标系的坐标轴是各自相互平行的，虽然这种假设通常是不成立的，但是在小范围测量，而且精度要求不是很高的情况下，该假设引起的误差可以忽略。 通常情况下，dX、dY、dZ 是每个椭球体相对WGS84椭球体的偏移量，比如由A坐标系转换成B坐标系,那么它们的偏移量是B相对于WGS84的偏移量减去A相对于WGS84的偏移量。 七参数法：（Bursa—Wolf） （dX,dY,dZ）：两坐标系的原点平移矢量（平移参数），原坐标系中的点位置矢量加上原点平移矢量即得到该点在新坐标系中的位置矢量。平移参数也就是原坐标系的原点在新坐标系中的坐标值。 （Rx, Ry, Rz）：位置矢量的旋转角（旋转参数）。参数符号约定如下：从直角坐标系原点，沿轴正向看，位置矢量绕轴顺时针旋转为正。从原坐标系转换到新坐标系，如果绕 Z 轴的旋转角度为正，那么转换后坐标点的经度将增大。。 M：位置矢量的比例因子（尺度比参数），位置矢量从原坐标系转换到新坐标系的尺度伸缩量。M = (1 + dS*10-6)，其中 dS 为尺度校正量，以百万分之一计（ppm）。 4.3大