当前位置:文档之家› 地理坐标到投影坐标转化方法理论1概要

地理坐标到投影坐标转化方法理论1概要

地理坐标系统和投影变换基础知识一、理论知识和背景介绍GIS处理的是空间信息,而所有对空间信息的量算都是基于某个坐标系统的,因此GIS中坐标系统的定义是GIS系统的基础,正确理解GIS中的坐标系统就变得尤为重要。

坐标系统又可分为两大类:地理坐标系统、投影坐标系统。

本文就对坐标系和投影及其在ArcGIS桌面产品中的应用做一些简单的论述。

GIS中的坐标系定义由基准面和地图投影两组参数确定,而基准面的定义则由特定椭球体及其对应的转换参数确定,因此欲正确定义GIS系统坐标系,首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。

1、地球椭球体(Ellipsoid)众所周知我们的地球表面是一个凸凹不平的表面,而对于地球测量而言,地表是一个无法用数学公式表达的曲面,这样的曲面不能作为测量和制图的基准面。

假想一个扁率极小的椭圆,绕大地球体短轴旋转所形成的规则椭球体称之为地球椭球体。

地球椭球体表面是一个规则的数学表面,可以用数学公式表达,所以在测量和制图中就用它替代地球的自然表面。

因此就有了地球椭球体的概念。

地球椭球体有长半径和短半径之分,长半径(a)即赤道半径,短半径(b)即极半径。

f =(a-b)/a为椭球体的扁率,表示椭球体的扁平程度。

由此可见,地球椭球体的形状和大小取决于a、b、f 。

因此,a、b、f被称为地球椭球体的三要素。

ArcGIS(ArcInfo)桌面软件中提供了30种地球椭球体模型;常见的地球椭球体数据见下表:椭球名称年代长半径扁率附注德兰勃1800 6 375 653 1:334.0 法国瓦尔别克1819 6 376 896 1:302.8 俄国埃弗瑞斯特1830 6 377 2761:300.801英国艾黎1830 6 376 5421:299.3 英国贝塞尔1841 6 377 3971:299.152 德国克拉克1856 6 377 8621:298.1英国克拉克1863 6 378 2881:294.4 英国克拉克1866 6 378 2061:294.978英国克拉克1880 6 378 2491:293.459 英国日丹诺夫1893 6 377 717 1:299.7俄国赫尔默特1906 6 378 1401:298.3 德国海福特1906 6 378 283 1:297.8美国赫尔默特1907 6 378 2001:298.3德国海福特1910 6 378 3881:297.01942年国际第一个推荐值热海景良1933 6 376 918 1:310.6日本川烟辛夫1935 6 377 0871:304.0日本克拉索夫斯基1940 6 378 2451:298.3苏联柯洛柯夫1955 6 378 2031:298.3苏联霍夫1956 6 378 2701:297.0美国WGS1960 6 378 1561:298.3美国国防部1960年世界大地坐标系弗希尔1960 6 378 1601:298.329美国凡氏(C一5)1965 6 378 1691:298.25美国施密森天文台美国施密森凡氏(C一5)1966 6 378 165 1:298.25天文台对地球椭球体而言,其围绕旋转的轴叫地轴。

地轴的北端称为地球的北极,南端称为南极;过地心与地轴垂直的平面与椭球面的交线是一个圆,这就是地球的赤道;过英国格林威治天文台旧址和地轴的平面与椭球面的交线称为本初子午线。

以地球的北极、南极、赤道和本初子午线等作为基本要素,即可构成地球椭球面的地理坐标系统(A geo graphic coordinate system (GCS) uses a three dimensional 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表达:图1 地理坐标系统2、大地基准面(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椭球体的参数可以参考常见的地球椭球体数据表。

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

地球椭球体和基准面之间的关系以及基准面是如何结合地球椭球体从而实现来逼近地球表面的可以通过图2一目了然。

图2 基准面定义椭球体拟合地表某一区域表面3、投影坐标系统(Projected Coordinate Systems )地球椭球体表面也是个曲面,而我们日常生活中的地图及量测空间通常是二维平面,因此在地图制图和线性量测时首先要考虑把曲面转化成平面。

由于球面上任何一点的位置是用地理坐标(λ,φ)表示的,而平面上的点的位置是用直角坐标(χ,у)或极坐标(r,)表示的,所以要想将地球表面上的点转移到平面上,必须采用一定的方法来确定地理坐标与平面直角坐标或极坐标之间的关系。

这种在球面和平面之间建立点与点之间函数关系的数学方法,就是地图投影方法。

接下来首先让我们来看看ArcGIS产品中对于北京54投影坐标系统的定义参数:Projection: Gauss_KrugerParameters:False_Easting: 500000.000000False_Northing: 0.000000Central_Meridian: 117.000000Scale_Factor: 1.000000Latitude_Of_Origin: 0.000000Linear Unit: Meter (1.000000)Geographic Coordinate System:Name: GCS_Beijing_1954Alias:Abbreviation:Remarks:Angular Unit: Degree (0.017453292519943299)Prime Meridian: Greenwich (0.000000000000000000)Datum: D_Beijing_1954Spheroid: Krasovsky_1940Semimajor Axis: 6378245.000000000000000000Semiminor Axis: 6356863.018773047300000000Inverse Flattening: 298.300000000000010000从参数中可以看出,每一个投影坐标系统都必定会有Geographic Coordinate Sys tem(地理坐标系统)。

那么我们从这一角度上解释一下投影和投影所需要的必要条件:将球面坐标转化为平面坐标的过程便是投影过程;投影所需要的必要条件是:第一、任何一种投影都必须基于一个椭球(地球椭球体),第二、将球面坐标转换为平面坐标的过程(投影算法)。

简单的说投影坐标系是地理坐标系+投影过程。

让我们从透视法(地图投影方法的一种)角度来直观的理解投影,图3。

几何透视法是利用透视的关系,将地球体面上的点投影到投影面(借助的几何面)上的一种投影方法。

如假设地球按比例缩小成一个透明的地球仪般的球体,在其球心或球面、球外安置一个光源,将球面上的经纬线投影到球外的一个投影平面上。

图3 透视法投影示意图投影既然是一种数学变换方法,那么任何一种投影都存在一定的变形,因此可以按照变形性质将投影方法如下分类:等角投影(Conformal Projection)、等积投影(E qual Area Projection)、等距投影(Equidistant Projection)、等方位投影(True-dir ection Projection)四种。

每种投影根据其名称就可以知道其方法保证了数据的那些几何属性,在实际应用过程中应根据需求来选取某种投影。

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

相关主题