正确区分测量学中的坐标系类型

正确区分测量学中的坐标系类型

黄河水利职业技术学院 地理信息0901 杨鹏飞 [email protected]

1与坐标系相关的几个基本概念

1.1对地球形状的描述

宇宙中地球是一个美丽的湛蓝湛蓝的球状星体，从机舱窗口俯视大地，地球是一个有些微起伏，极其复杂的表面，珠穆朗玛峰海拔8848.43米，马里亚纳海沟-11022米，实时是地球不是一个正球体，而是一个极半径略短，赤道半径略突出，南极略扁平，近似于梨形的椭球体

一级逼近:水准面水准面(大地体)

设想一个与静止的平均海水面重合并延伸到大陆内部的包围整个地球的封闭的重力位水准面。

二级逼近:旋转椭球面 地球椭球

大地水准面的形状接近于一个两级略短的旋转椭球，通常用旋转椭球代表地球，作为描述地球空间位置的基准 长轴a 、短轴b 、椭球扁率f=a-b/a

三级逼近:参考椭球体定位

对地球的a 、b 、f 测定后，还必须测定大地水准面相对关系，即确定一个局部地区与大地水准面符合最好的椭球体，并通过数学方法计算出两者各点间的偏差。

1.2地球椭球体(Ellipsoid)

众所周知我们的地球表面是一个凸凹不平的表面，而对于地球测量而言，地表是一个无法用数学公式表达的曲面，这样的曲面不能作为测量和制图的基准面。假想一个扁率极小的椭圆，绕大地球体短轴旋转所形成的规则椭球体称之为地球椭球体。地球椭球体表面是一个规则的数学表面，可以用数学公式表达，所以在测量和制图中就用它替代地球的自然表面。因此就有了地球椭球体的概念。

地球椭球体有长半径和短半径之分，长半径(a)即赤道半径，短半径(b)即极半径。f=（a-b ）/a为椭球体的扁率，表示椭球体的扁平程度。由此可见，地球椭球体的形状和大小取决于a 、b 、f 。因此，a 、b 、f 被称为地球椭球体的三要素。

对地球椭球体而言，其围绕旋转的轴叫地轴。地轴的北端称为地球的北极，南端称为南极；过地心与地轴垂直的平面与椭球面的交线是一个圆，这就是地球的赤道；过英国格林威治天文台旧址和地轴的平面与椭球面的交线称为本初子午线。以地球的北极、南极、赤道和本初子午线等作为基本要素，即可构成地球椭球面的地理坐标系统(A geographic coordinate system (GCS) uses a threedimensional spherical surface to define locations on the earth.A GCS includes an angular unit of measure, a prime meridian,and a datum (based on a spheroid).) 。可以看出地理坐标系统是球面坐标系统，以经度/维度（通常以十进制度或度分秒(DMS)的形式）来表示地面点位的位置。

地理坐标系统以本初子午线为基准（向东，向西各分了1800）之东为东经其值为正，之西为西经其值为负；以赤道为基准（向南、向北各分了900）之北为北纬其值为正，之南为南纬其值为负。

1.3大地基准面（Geodetic datum）

大地基准面（Geodetic datum），设计用为最密合部份或全部大地水准面的数学模式。它由椭球体本身及椭球体和地表上一点视为原点间之关系来定义。此关系能以 6个量来定

义，通常（但非必然）是大地纬度、大地经度、原点高度、原点垂线偏差之两分量及原点至某点的大地方位角。

让我们先抛开测绘学上这个晦涩难懂的概念，看看GIS 系统中的基准面是如何定义的，GIS 中的基准面通过当地基准面向WGS1984的转换7参数来定义，转换通过相似变换方法实现，具体算法可参考科学出版社1999年出版的《城市地理信息系统标准化指南》第76至86页。假设Xg 、Yg 、Zg 表示WGS84地心坐标系的三坐标轴，Xt 、Yt 、Zt 表示当地坐标系的三坐标轴，那么自定义基准面的7参数分别为：三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值；三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时，分别绕Xt 、Yt 、Zt 的旋转角；最后是比例校正因子，用于调整椭球大小。

那么现在让我们把地球椭球体和基准面结合起来看，在此我们把地球比做是“马铃薯”，表面凸凹不平，而地球椭球体就好比一个“鸭蛋”，那么按照我们前面的定义，基准面就定义了怎样拿这个“鸭蛋”去逼近“马铃薯”某一个区域的表面，X 、Y 、Z 轴进行一定的偏移，并各自旋转一定的角度，大小不适当的时候就缩放一下“鸭蛋”，那么通过如上的处理必定可以达到很好的逼近地球某一区域的表面。因此，从这一点上也可以很好的理解，每个国家或地区均有各自的基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的1975地球椭球体（IAG75）建立了我国新的大地坐标系--西安80坐标系，目前大地测量基本上仍以北京54坐标系作为参照，北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。 WGS1984基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心，目前GPS 测量数据多以WGS1984为基准。

克拉索夫斯基(Krassovsky)、1975地球椭球体（IAG75）、WGS1984椭球体的参数可以参考常见的地球椭球体数据表。

椭球体与基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面。

1.4投影坐标系统（Projected Coordinate Systems ）

地球椭球体表面也是个曲面，而我们日常生活中的地图及量测空间通常是二维平面，因此在地图制图和线性量测时首先要考虑把曲面转化成平面。由于球面上任何一点的位置是用地理坐标（λ，φ）表示的，而平面上的点的位置是用直角坐标（χ，у）或极坐标（r ， ）表示的，所以要想将地球表面上的点转移到平面上，必须采用一定的方法来确定地理坐标与平面直角坐标或极坐标之间的关系。这种在球面和平面之间建立点与点之间函数关系的数学方法，就是地图投影方法。

接下来首先让我们来看看ArcGIS 产品中对于北京54投影坐标系统的定义参数： Projection: Gauss_Kruger

Parameters:

False_Easting: 500000.000000

False_Northing: 0.000000

Central_Meridian: 117.000000

Scale_Factor: 1.000000

Latitude_Of_Origin: 0.000000

Linear Unit: Meter (1.000000)

Geographic Coordinate System:

Name: GCS_Beijing_1954

Alias:

Abbreviation:

Remarks:

Angular Unit: Degree (0.[***********])

Prime Meridian: Greenwich (0.[***********])

Datum: D_Beijing_1954

Spheroid: Krasovsky_1940

Semimajor Axis: 6378245.[***********]

Semiminor Axis: 6356863.[***********]

Inverse Flattening: 298.[***********]

从参数中可以看出，每一个投影坐标系统都必定会有Geographic Coordinate System（地理坐标系统）。那么我们从这一角度上解释一下投影和投影所需要的必要条件：将球面坐标转化为平面坐标的过程便是投影过程；投影所需要的必要条件是：第一、任何一种投影都必须基于一个椭球（地球椭球体），第二、将球面坐标转换为平面坐标的过程（投影算法）。简单的说投影坐标系是地理坐标系+投影过程。

让我们从透视法（地图投影方法的一种）角度来直观的理解投影，图2。

几何透视法是利用透视的关系，将地球体面上的点投影到投影面（借助的几何面）上的一种投影方法。如假设地球按比例缩小成一个透明的地球仪般的球体，在其球心或球面、球外安置一个光源，将球面上的经纬线投影到球外的一个投影平面上。

投影既然是一种数学变换方法，那么任何一种投影都存在一定的变形，因此可以按照变形性质将投影方法如下分类：等角投影（Conformal Projection） 、 等积投影（Equal Area Projection）、等距投影（Equidistant Projection）、等方位投影（True-direction Projection ）四种。每种投影根据其名称就可以知道其方法保证了数据的那些几何属性，在实际应用过程中应根据需求来选取某种投影。

如果按照投影的构成方法分类又可分为方位、圆柱、圆锥投影三种，在上述三种投影中由于几何面与球面的关系位置不同，又分为正轴、横轴和斜轴三种。

接下来我们来看看我们国家通常采用的投影——高斯—克吕格（Gauss-Kruger) 投影，是一种“等角横切圆柱投影”。德国数学家、物理学家、天文学家高斯（Carl Friedrich Gauss ，1777一 1855）于十九世纪二十年代拟定，后经德国大地测量学家克吕格（Johannes Kruger ，1857～1928）于 1912年对投影公式加以补充，故名。设想用一个圆柱横切于球面上投影带的中央经线，按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件，将中央经线两侧一定经差范围内的球面正形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开展平，即获高斯-克吕格投影平面。高斯—克吕格投影后，除中央经线和赤道为直线外，其他经线均为对称于中央经线的曲线。高斯—克吕格投影没有角度变形，在长度和面积上变形也很小，中央经线无变形，自中央经线向投影带边缘，变形逐渐增加，变形最大处在投影带内赤道的两端。按一定经差将地球椭球面划分成若干投影带, 这是高斯投影中限制长度变形的最有效方法。

分带时既要控制长度变形使其不大于测图误差，又要使带数不致过多以减少换带计算工作，据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带, 以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自 0度子午线起每隔经差6度自西向东分带，带号依次编为第 1、2„60带。三度带是在六度带的基础上分成的，它的中央子午线与六度带的中央子午线和分带子午线重合，即自 1.5度子午线起每隔经差3度自西向东分带，带号依次编为三度带第 1、2„120带。我国的经度范围西起 73°东至135°，可分成六度带十一个，各带中央经线依次为75°、81°、87°、„„、117°、123°、129°、135°，

或三度带二十二个。

我国大于等于50万的大中比例尺地形图多采用六度带高斯—克吕格投影，三度带高斯—克吕格投影多用于大比例尺1：1万测图，如城建坐标多采用三度带的高斯—克吕格投影。高斯—克吕格投影按分带方法各自进行投影，故各带坐标成独立系统。以中央经线（L0）投影为纵轴X, 赤道投影为横轴Y ，两轴交点即为各带的坐标原点。为了避免横坐标出现负值，高斯—克吕格投影北半球投影中规定将坐标纵轴西移500公里当作起始轴。由于高斯—克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值，所以各带的坐标完全相同，为了区别某一坐标系统属于哪一带，通常在横轴坐标前加上带号，如(4231898m,21655933m)，其中21即为带号。

2坐标系的分类

2.1按照对数据的存储单位(对同一地面点的不同描述)

地理坐标系、投影平面直角坐标系、空间直角坐标系、用户自定义坐标系

2.2按照基准面不同分类

我国两大基准面：北京54坐标系、西安80坐标系

3 mapgis 中大地坐标系的解释

MAPGIS 中的大地坐标系其实是投影平面直角坐标系高斯克吕格投影类型中的一个情况，比例尺分母为1，单位为米。因为此时的图形坐标和实际测量的大地坐标是一致的，所以成为大地坐标系。

测量学中的大地坐标系并不是上述的含义，它是大地地理坐标系的简称。地球椭球面上任一点的位置，可由该点的纬度(B)和经度(L)确定，即地面点的地理坐标值，由经线和纬线构成两组互相正交的曲线坐标网叫地理坐标网。由经纬度构成的地理坐标系统又叫地理坐标系。地理坐标分为天文地理坐标和大地地理坐标，天文地理坐标是用天文测量方法确定的，大地地理坐标是用大地测量方法确定的。我们在地球椭球面上所用的地理坐标系属于大地地理坐标系，简称大地坐标系。

mapgis 如何知道它是什么坐标系的.1、经纬网的形状（矩形表示是地理坐标，梯形表示平面真角坐标）； 2、公里网的形状（水平垂直表示原图已几何变换，有可能是老的地图（北京54坐标的图常这样））； 3、看左下角的坐标与标出来的理论坐标（经纬度或是大地坐标）的关系。 4、试图寻找出中央垂直经线。 5、试图寻找出图的比例尺。（可以在工程输出版面编排得出图的大概大小）。根据经验，125万、1：20万、1：10万、1：5万、1：

2.5万的图大多是平面直角坐标系高斯投影；1：50万大多是等角割圆锥投影