地理坐标到投影坐标转化方法理论_第1页
地理坐标到投影坐标转化方法理论_第2页
地理坐标到投影坐标转化方法理论_第3页
地理坐标到投影坐标转化方法理论_第4页
地理坐标到投影坐标转化方法理论_第5页
已阅读5页,还剩18页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、地理坐标系统和投影变换基础知识一、理论知识和背景介绍GIS处理的是空间信息,而所有对空间信息的量算都是基于某个坐标系统的,因此GIS中坐标系统的定义是 GIS系统的基础,正确理解 GIS中的坐标系统就变得尤为重 要。坐标系统又可分为两大类:地理坐标系统、投影坐标系统。本文就对坐标系和投影及其在ArcGIS桌面产品中的应用做一些简单的论述。GIS中的坐标系定义由基准面和地图投影两组参数确定,而基准面的定义则由特定椭球体及其对应的转换参数确定,因此欲正确定义GIS系统坐标系,首先必须弄清地球 椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它

2、们之 间的关系。1、地球椭球体(Ellipsoid)众所周知我们的地球表面是一个凸凹不平的表面,而对于地球测量而言,地表是一个无法用数学公式表达的曲面,这样的曲面不能作为测量和制图的基准面。假想一个扁率极小的椭圆,绕大地球体短轴旋转所形成的规则椭球体称之为地球椭球体。地球椭球 体表面是一个规则的数学表面,可以用数学公式表达,所以在测量和制图中就用它替代地球的自然表面。因此就有了地球椭球体的概念。地球椭球体有长半径和短半径之分,长半径(a)即赤道半径,短半径(b)即极半径。f=(a-b) /a为椭球体的扁率,表示椭球体的扁平程度。由此可见,地球椭球体的形状和 大小取决于a、b、f。因此,a、b、

3、f被称为地球椭球体的三要素。ArcGIS(ArcInfo)桌面软件中提供了 30种地球椭球体模型;常见的地球椭球体数据见 下表:椭球名称年代长半径扁率附注1德兰勃18006 375 6531:334.0法国瓦尔别克18196 376 8961 : 302.8俄国埃弗瑞斯特18306 377 2761: 300.801英国18306 376 5421: 299 . 3英国贝塞尔18416 377 3971: 299 . 152德国克拉克18566 377 8621: 298 . 1英国克拉克18636 378 2881: 294 . 4英国克拉克18666 378 2061: 294 , 978

4、英国克拉克18806 378 2491: 293 . 459英国日丹诺夫18936 377 7171: 299 . 7俄国赫尔默特19066378 1401: 298 . 3德国海福特19066 378 2831 : 297.8美国赫尔默特19076 378 2001: 298 . 3德国海福特19106 378 3881: 297 . 01942年国际第一个推荐值热海景良19336 376 9181: 310 . 6日本川烟辛夫19356 377 0871: 304 . 0日本克拉索夫斯基19406 378 2451 : 298.3苏联柯洛柯夫19556 378 2031 : 298.3苏联

5、口19566 378 2701 : 297.0美国WGS19606378 1561 : 298.3美国国防部1960 ”界大地坐标系弗希尔19606378 1601: 298.329美国凡氏(C 一 5)19656378 1691 : 298.25美国施密森天文台凡氏(C 一 5)19666378 1651 : 298.25美国施密森天文台对地球椭球体而言,其围绕旋转的轴叫地轴。地轴的北端称为地球的北极,南端称 为南极;过地心与地轴垂直的平面与椭球面的交线是一个圆,这就是地球的赤道;过英 国格林威治天文台旧址和地轴的平面与椭球面的交线称为本初子午线。以地球的北极、 南极、赤道和本初子午线等作为

6、基本要素,即可构成地球椭球面的地理坐标系统(A geographic coordinate system (GCS) uses a three dimensional spherical surfaceto define locations on the earth. A GCS includes an angular unit of measure, a prime meridian, and a datum (based on a spheroid).) 。可以看出地理坐标系 统是球面坐标系统,以经度/维度(通常以十进制度或度分秒(DMS)的形式)来表示地面 点位的位置。地理坐标系统以本初

7、子午线为基准 (向东,向西各分了 1800)之东为东经其值为正, 之西为西经其值为负;以赤道为基准(向南、向北各分了90°)之北为北纬其值为正,之南为南纬其值为负。地表任意位置的坐标值可由图 1表达:Pamllftlft(LfintB of IfllttudeIAerk9taM(LJnM of longhiidv)Ntwork2 、大地基准面(Geodetic datum )大地基准面( Geodetic datum ) , 设计用为最密合部份或全部大地水准面的数学模式。它由椭球体本身及椭球体和地表上一点视为原点间之关系来定义。此关系能以 6个量来定义,通常(但非必然)是大地纬度、大

8、地经度、原点高度、原点垂线偏差之两分量及原点至某点的大地方位角。让我们先抛开测绘学上这个晦涩难懂的概念,看看GIS 系统中的基准面是如何定义的, GIS 中的基准面通过当地基准面向 WGS1984 的转换 7 参数来定义, 转换通过相似变换方法实现,具体算法可参考科学出版社1999 年出版的城市地理信息系统标准化指南第 76 至 86 页。假设 Xg 、 Yg 、 Zg 表示 WGS84 地心坐标系的三坐标轴, Xt、Yt 、 Zt 表示当地坐标系的三坐标轴, 那么自定义基准面的 7 参数分别为: 三个平移参数 X、AY、 Z表小两坐标原点的平移值;三个旋转参数ex e y e z表不当地坐标

9、系旋转至与地心坐标系平行时,分别绕Xt 、 Yt 、 Zt 的旋转角;最后是比例校正因子,用于调整椭球大小。那么现在让我们把地球椭球体和基准面结合起来看,在此我们把地球比做是 “马铃薯 ” ,表面凸凹不平,而地球椭球体就好比一个 “鸭蛋 ”,那么按照我们前面的定义,基准面就定义了怎样拿这个 “鸭蛋 ”去逼近 “马铃薯 ”某一个区域的表面, X 、 Y 、 Z 轴进行一定的偏移,并各自旋转一定的角度,大小不适当的时候就缩放一下 “鸭蛋 ”,那么通过如上的处理必定可以达到很好的逼近地球某一区域的表面。因此,从这一点上也可以很好的理解,每个国家或地区均有各自的基准面,我们通常称谓的北京 54 坐标系

10、、 西安 80 坐标系实际上指的是我国的两个大地基准面。 我国参照前苏联从1953 年起采用克拉索夫斯基(Krassovsky) 椭球体建立了我国的北京 54 坐标系, 1978 年采用国际大地测量协会推荐的 1975 地球椭球体( IAG75 )建立了我国新的大地坐标系-西安80 坐标系,目前大地测量基本上仍以北京54 坐标系作为参照,北京54 与西安 80 坐标之间的转换可查阅国家测绘局公布的对照表。 WGS1984 基准面采用 WGS84 椭球体,它是一地心坐标系,即以地心作为椭球体中心,目前GPS 测量数据多以 WGS1984 为基准。克拉索夫斯基(Krassovsky)、1975地球

11、椭球体(IAG75 )、WGS1984椭球体的参数可以参考常见的地球椭球体数据表。椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面。地球椭球体和基准面之间的关系以及基准面是如何结合地球椭球体从而实现来逼近地球表面的可以通过图2 一目了然。LocbI geographic coordmate systemEarth's surfaceEarth-centered (WGSA4) datumLocal (NA027) datumEarlh-rnitered geographic coordinate syst

12、em图2基准面定义椭球体拟合地表某一区域表面3、投影坐标系统(Projected Coordinate Systems )地球椭球体表面也是个曲面,而我们日常生活中的地图及量测空间通常是二维平面,因此在地图制图和线性量测时首先要考虑把曲面转化成平面。由于球面上任何一点的位置是用地理坐标(入,6)表示的,而平面上的点的位置是用直角坐标( 为y)或极坐标(r,)表示的,所以要想将地球表面上的点转移到平面上,必须采用一定的方法 来确定地理坐标与平面直角坐标或极坐标之间的关系。这种在球面和平面之间建立点与点之间函数关系的数学方法,就是地图投影方法。接下来首先让我们来看看 ArcGIS产品中对于北京54

13、投影坐标系统的定义参数: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: De

14、gree (0.017453292519943299)Prime Meridian: Greenwich (0.000000000000000000)Datum: D_Beijing_1954Spheroid: Krasovsky_1940Semimajor Axis: 6378245.000000000000000000Semiminor Axis: 6356863.018773047300000000Inverse Flattening: 298.300000000000010000从参数中可以看出, 每一个投影坐标系统都必定会有Geographic Coordinate System (

15、地理坐标系统)。那么我们从这一角度上解释一下投影和投影所需要的必要条件: 将球面坐标转化为平面坐标的过程便是投影过程;投影所需要的必要条件是:第一、任 何一种投影都必须基于一个椭球(地球椭球体),第二、将球面坐标转换为平面坐标的 过程(投影算法)。简单的说投影坐标系是地理坐标系+投影过程。让我们从透视法(地图投影方法的一种)角度来直观的理解投影,图 3。几何透视法是利用透视的关系,将地球体面上的点投影到投影面(借助的几何面)上的一种投影方法。如假设地球按比例缩小成一个透明的地球仪般的球体,在其球心或球面、球外安置一个光源,将球面上的经纬线投影到球外的一个投影平面上。图3透视法投影示意图投影既然

16、是一种数学变换方法,那么任何一种投影都存在一定的变形,因此可以按照变形性质将投影方法如下分类:等角投影(Conformal Projection )、等积投影(Equal Area Projection )、等距投影(Equidistant Projection )、等方位投影(True-dir ection Projection )四种。每种投影根据其名称就可以知道其方法保证了数据的那些几 何属性,在实际应用过程中应根据需求来选取某种投影。如果按照投影的构成方法分类又可分为方位、圆柱、圆锥投影三种,在上述三种投影中由于几何面与球面的关系位置不同,又分为正轴、横轴和斜轴三种。图3-5将直观展现

17、上述各种投影。图4正、横、斜圆柱投影示意图图5正、横、斜圆锥投影示意图图6正、横、斜方位投影示意图接下来我们来看看我们国家通常采用的投影 高斯一克吕格(Gauss-Kruger)投 影,是一种 等角横切圆柱投影德国数学家、物理学家、天文学家高斯( Carl Friedri ch Gauss ,1777 1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(Johannes Kruger, 18571928 )于1912年对投影公式加以补充,故名。设想用一个圆柱横切于球面上投影带的中央经线,按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件,将中央经线两侧一定经差范围内的球面正

18、形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开展平,即获高斯-克吕格投影平面。高斯 一克吕格投影后,除中央经线和赤道为直线外,其他经线均为对称于中央经线的曲线。高斯 一克 吕格投影没有角度变形,在长度和面积上变形也很小,中央经线无变形,自中央经线向 投影带边缘,变形逐渐增加,变形最大处在投影带内赤道的两端。按一定经差将地球椭 球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自0度子午线

19、起每隔经差6度自西向东分带,带号依次编为第1、2 -60带。三度带是在六度带的基础上分成的,它的中央子午线与六度带的中央子午线和分带子午线重合,即自1.5度子午线起每隔经差3度自西向东分带,带号依次编为三度带第1、2- -120带。我国的经度范围西起 73。东至135° ,可分成六度带十一个,各带中央经线依次为75°、81°、87°、117°、123°、129°、135°,或三度带二十二个。我国大于等于50万的大中比例尺地形图多采用六度带高斯一克吕格投影,三度带高斯一克吕格投影多用于大比例尺 1: 1万测图,如城建

20、坐标多采用三度带的高斯一克吕 格投影。高斯一克吕格投影按分带方法各自进行投影,故各带坐标成独立系统。以中央 经线(L0)投影为纵轴X,赤道投影为横轴Y,两轴交点即为各带的坐标原点。为了避 免横坐标出现负值,高斯一克吕格投影北半球投影中规定将坐标纵轴西移500公里当作起始轴。由于高斯一克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所 以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,通常在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。高斯一克吕格投影及分带示意图 如下:图7高斯一克吕格投影及分带示意图接下来让我们来看看ArcGIS 中对我们国家经

21、常采用北京 54和西安 80 坐标系统是怎么样描述的,在ArcMap 或是 ArcCatalog 中选择系统预定义的北京54 和西安 80 坐标系统。在$ArcGISHomeCoordinateSystemsCoordinate SystemsProjected CoordinateSystemsGauss KrugerBeijing 1954 目录中,我们可以看到四种不同的命名方式:Beijing 1954 3 Degree GK CM 75E.prjBeijing 1954 3 Degree GK Zone 25.prjBeijing 1954 GK Zone 13.prjBeijing

22、1954 GK Zone 13N.prj对它们的说明分别如下:三度分带法的北京 54 坐标系,中央经线在东75 度的分带坐标,横坐标前不加带号三度分带法的北京 54 坐标系,中央经线在东75 度的分带坐标,横坐标前加带号六度分带法的北京 54 坐标系,分带号为 13 ,横坐标前加带号六度分带法的北京 54 坐标系,分带号为 13 ,横坐标前不加带号在 Coordinate SystemsProjected Coordinate SystemsGauss KrugerXian 1980目录中,文件命名方式又有所变化:Xian19803 DegreeGK CM 75E.prjXian19803 D

23、egreeGK Zone 25.prjXian1980GK CM 75E.prjXian1980GK Zone13.prj西安 80 坐标文件的命名方式、含义和北京 54 前两个坐标相同,但没有出现 “带号 +N” 这种形式,为什么没有采用统一的命名方式?让人看了的确有些费解, 大家在应用过程中需要特别注意一下。ArcGIS 的应用在我们了解了坐标系统和投影的定义和其内在的联系后, 本文着重总结一下坐标系统和投影变换在桌面产品(版本9.2 )中的应用(分ArcMap 、 ArcCatalog 、 ArcToolBox 三大主要应用模块)。1、动态投影(ArcMap)所谓动态投影指:改变 Arc

24、Map中的Data Frame (工作区)的空间参考或是对后 加入到ArcMap工作区中数据的投影变换。ArcMap的Data Frame (工作区)的坐标系统默认为第一个加载到当前Data Frame (工作区)的那个文件的坐标系统,后加入的数据,如果和当前工作区坐标系统不同,则 ArcMap会自动做投影变换,把后加入的 数据投影变换到当前坐标系统下显示,但此时数据文件所存储的实际数据坐标值并没有 改变,只是显示形态上的变化!因此叫动态投影。表现这一点最明显的例子就是在 Exp ort Data时,用户可以选择是按 this layer's source data (数据源的坐标系统

25、导出), 还是按照the Data Frame (当前工作区的坐标系统)导出数据。关于ArcMap的这种动态投影机制,我们可以利用一个北京 54投影坐标系数据(乡镇.shp )和 ArcGIS Installation目录下的world30.shp数据来做一个实验说明乡镇.shp数据的坐标系统为北京 54投影坐标系(Krasovsky_1940_Transverse_Mercator)。在 ArcMap或ArcCatalog中预览形态如图 8所示:(knerfl 门“财前占。£炉工”1h嘉氧工也门”g山M1: 士H Ml:r fed :i! aA川诵3g区£皿加皿 O WO

26、OWLcr射;”露©0WW3 i.OCIMOOLifW«r UM: HSfer (I C£»0(WEastern.叱”占小小9+0.MquLtf LMJ Degr«t (Q,OI 7*532*125 W3ZW| 正出悖七方林:加的曲酬趣酬或的OjOeJfeun: DJ1好加.小Jphsfod: 口必切密9M 知加部Df *1的 63海石.OCgXKOKHKMXXMMKv二 I *KtQpiX中igiwc小3彳皆 fEi,即E事g/*MW 杓例 冲J HY £ 皿 MhwJ das£r rsbefKDcatc- a new 0

27、«丽虚2I Edt g ecmbtIW &审*由寸 4f4 Mlec 睡&M如灯,二J1k4事 I M (Ft犷0刘h S ggwrb图8北京54投影坐标系数据单独显示几何形态world30.shp数据的坐标系统为 WGS84坐标系(GCS_WGS_1984)。在ArcMap或ArcCatalog中预览形态如图9所示:卜 _J StTWl411.*1SlTXT Can-lijii.Lv 5nt- Paa-ldc | InlvriT硒5irfi lt( 沪 r *u|wt +事? I'X |g; >qr=e (0.0l74M292S+:32?9)Pmw M

28、eocfep- eemnehtWtun; D_WG5如WS5_L9MMnEor 命佰:37训#.000000000001»00(»0Select & pEdrfTwdImport 导 coorctaute systMn srd ;小/Z and M tkxr&金 frem -an t:曲 E '3«xiaft55et,底而re 出t 必* f 2b5 P讪4, r闻erj.(rwtbR m rws mminxMSenwnrw 苦ds: 63567521112453 万:KCDOEW liiwse Flattening: 2W.M7223S&

29、amp;3CI0M3IMMC-e赞.Mm图9 WGS84坐标系数据单独显示几何形态而在ArcMap中先加载北京54坐标系数据后再加入 WGS84坐标系数据,让 ArcMap对WGS84坐标系数据进行动态投影后两数据叠加显示效果如图10所示:Ida. "口步占岂M - =1针汇/HAL 迦 油1- a 一5H JT 蜃 A图10 ArcMap对WGS84数据进行动态投影后的显示状态可以非常明显的看到 ArcMap对WGS84数据做完动态投影后的数据几何形态上的改变,并且此时从ArcMap右下角的状态栏上也可以看到当前Data Frame (工作空间)的坐标系统为北京54平面投影坐标系统

30、。反之在ArcMap中先加载 WGS84坐标系数据后再加入北京 54坐标系数据,让ArcMap对北京54坐标系数据进行动态投影后两数据叠加显示效果如图11所示:En IA.I Ei.卬a IErJd 3*口 好百也 =1钞IVi.MI.H ee”8+ ETWEn*Gcwi-d CrN FraTri'-«nwIw" Lim 坤mmiiH三向三1 » - A - r*T=aIW 71r ZW :lMZ图11 ArcMap对北京54坐标系数据进行动态投影后的显示状态将在图10中动态投影后的 WGS84坐标系统数据按系统框架坐标系统导出后,单独加载或预览的数据几何

31、形态如图12 :tih X4ii m/q / A r* A 、i向%?I- rrs n班0单*¥ ®边 O Si!fi-airl 1LMMiU ,'AL'lMr:! ! d*.,d *卡/i_h,_j":洱 EwXniFiEi* Ai hruMihl。空出sijii4liiHi. f图12按北京54坐标系统框架导出 WGS84数据后的数据几何形态显示通过上述实验能够很好的说明ArcMap的动态投影特性。2、坐标系统描述(ArcCatalog)大家都知道在ArcCatalog中可以给数据定义坐标系统描述!即在数据上鼠标右键->Propertie

32、s->XY Coordinate System 选项卡,这里可以通过 New、Modify、Select、I mport方式来为数据定义坐标系统描述。但有许多用户都认为在这里定义了数据的坐标 系统信息后,其数据本身就发生了投影变换。其实不然,这里定义的数据坐标系统信息 都对应到与该数据同名而后缀名为.prj的文件当中!如果把该文件删除,在 ArcCatalog中重新查看(要在该数据的上层节点上Refresh刷新一下)该文件的坐标信息时,一样会显示为Unknown ,并且数据的坐标值并没有发生实质上的投影变换,这里改的仅仅 是对数据坐标系统信息的一个描述而已,这就好比我们每个人的基本信息登

33、记卡,更改了登记信息,但并没有改变你这个人本身。因此数据文件中所存储数据的坐标值并没有 真正的投影变换到你想要更改到的坐标系统下。我们同样拿上述的两个数据做一下实验,在ArcCatalog中更改world30.shp的坐标 系统描述,在 world30.shp 文件上鼠标右键->Properties->XY Coordinate System 选项 卡中,通过 Import 方式导入乡镇.shp 文件的 Krasovsky_1940_Transverse_Mercator 投影坐标系统描述,之后看一下结果图13。13 更改坐标系统描述后的数据几何形态从上述示例我们可以很明显的看到更

34、改数据的坐标系统描述并不能使数据做投影变换,从而使数据投影到平面上来,但该数据的 prj 文件已经记录了更改后的坐标系统描述,PROJCS"Krasovsky_1940_Transverse_Mercator",GEOGCS"GCS_Krasovsky_194 0",DATUM"D_Krasovsky_1940",SPHEROID"Krasovsky_1940",6378245.0,298.3,PRIMEM"Greenwich",0.0,UNIT"Degree",0.0174

35、532925199433,PROJECTION"Transverse_Mercator",PARAMETER"False_Easting",500000.0,PARAMETER"False_Northing",0.0,PARAMETER"Central_Meridian",111.0,PARAMETER"Scale_Factor",1.0,PARAMETER"Latitude_Of_Origin",0.0,UNIT"Meter",1.0但对数据坐标系统的这

36、个描述也是非常重要的, 如果我们拿到一份数据, 从 ArcMa p 下所显示的坐标来看,像是投影坐标系统下的平面坐标,但不知道是基于哪个椭球体的什么投影方法,因此就无法再对数据做进一步的处理,如:投影变换、量测等操作。因为我们无法得知从什么坐标系统下开始变换,以及该坐标系统下的量测单位是什么。因此大家一定要更正对 ArcCatalog 中数据属性中关于坐标系统描述的认识。3 、投影变换(ArcToolBox)上面说了这么多,可能有人要问: “要真正的改变数据的坐标值该怎么办? ”也就是 做真正的投影变换。在ArcToolBox->Data Management Tools->Pro

37、jections and Transformations 中提供了如下工具:在这个工具集下有这么几个工具最为常用:1、 Define Projection2、 Feature->Project3、 Raster->Project Raster4、 Create Custom Geographic Transformation当数据在没有任何空间参考信息时,在 ArcCatalog 的坐标系统描述( XY Coordinate System )选项卡中会显示为Unknown !这时如果要对数据进行投影变换就要先利用 Define Projection 工具来给数据定义一个Coordi

38、nate System , 然后再利用Feature->Project 或 Raster->Project Raster 工具来对数据进行投影变换!由于我们国家经常使用的坐标系统为北京 54 和西安 80 。 这两个坐标系统变换到其他坐标系统下时,通常需要提供一个Geographic Transformation ,因为不同投影所基于的椭球体及Datum 不同!关键是Datum 不同,也就是说当两个投影基于不同的 Datum 时就需要制定参数做Geographic Transformation 。 这里就用到我们前面所说的转换3参数、转换7参数了(三个平移参数 AX、AY、AZ表示

39、两坐标原点的平移值;三个旋转参数ex £y e/示当地坐标系旋转至与地心坐标系平行时,分别绕Xt 、 Yt 、 Zt 的旋转角;最后是比例校正因子,用于调整椭球大小),而我们国家的转换参数是保密的, 因此可以自己计算或在购买数据时向国家测绘部门索要。实际工作中一般都根据工作区内已知的北京54 坐标控制点计算转换参数,如果工作区内有足够多的已知北京54 与 WGS84 坐标控制点,可直接计算坐标转换所需的 7参数或 3 参数;当工作区内有3 个已知北京54 与 WGS84 坐标控制点时,可用下式计算 WGS84 到北京 54 坐标的转换参数(A、 B、 C 、 D 、 E 、 F) :

40、x54 = AX84 + BY84 + Cy54 = DX84 + EY84 + F多余一点用作检验;在只有一个已知控制点的情况下 (往往如此) ,用已知点的北京54与 WGS84 坐标之差作为平移参数,当工作区范围不大时精度也足够了。那么当精度要求较高,实测数据为 WGS1984 坐标数据时,欲转换到北京 54 基准面的高斯 -克吕格投影坐标,如何定义坐标系参数呢?你可选择WGS 1984 作为基准面,当只有一个已知控制点时, 根据平移参数调整东伪偏移、北纬偏移值实现WGS84 到北京54的转换,此时的转换系数 (A、B、C、D、E、F)中A、B、D、E为0,只有X、Y方向的平移值C、 F

41、;当有 3 个已知控制点时,可利用得到的转换系数(A 、 B、 C 、 D、 E 、F)定义AffineTransform 坐标系变换对象,实现坐标系的转换,当然有足够多已知控制点时,直接求定7 参数自定义基准面就行了。Geographic Transformation通常需要指定变换方向的,如果没有指定变换的方向,ArcMap会自动正确应用变换方法如从WGS 1984到NAD1927 ,我们就可以选择 NAD_1927_to_WGS_1984_3 这种变换方法,ArcMap会自动判定转换方向从而正确实现 数据的投影变换。有些变换方法是先把地理坐标系(经纬度)变换到地心坐标系(X YZ),然后

42、再将地心坐标系坐标(X Y Z),变换到地理坐标系统,图 14。地理坐标系统地心坐标系统地理坐标系统,图14从地理坐标系统变换到地心坐标系统再到地理坐标系统Geocentric Transformation© 参数)、Coordinate Frame(7 参数)、Molodensky-Badekas(10参数,是7参数变换方法的一种变体,多了定义XYZ旋转轴的起始位置)都是这样的方法。Molodensky(3 参数)和 Abridged Molodensky(3 参数是 Molodensky 方法的一种变 体)变化方法则是直接在不同的地理坐标系统下进行变换无需借助地心坐标系进行转 换

43、。那么上述的这些方法的精度取决与用于定义转换参数的控制点的质量个数及选择的 变换方法,精度相差可从厘米到米。美国国家大地测量局用一种基于格网的变换方法(NADCON和HARN)在NAD1927和NAD1983及其他地理坐标系统下做变换,彼此相临的州的精度可达0.15米,阿拉斯加州及周边岛屿的精度可达 0.5米,夏威夷地区精度可达 0.2米等等,精度取决于计 算栅格大小时测量数据的质量。随着卫星测距和测量技术的改善,美国在更新了测量控 制网后,推出了这种 High Accuracy Reference Network (HARN)方法或说是栅格格 网,其精度可达0.05米。加拿大采用同 NADC

44、ON 类似的基于格网的方法 (NTV2)在NAD 1927和NAD 1983之间做变换,这种方法采用了双线性内插的方法来计算点的坐标。关于这两种方法 本文不做更多描述,详细信息可参考ArcGIS帮助。知道转换参数后,可以利用Create Custom Geographic Transformation 工具来定 义一个Geographic Transformation方法,变换方法可以根据已知的转换参数个数选择 变换方法,这就完成了对数据的投影变换,数据本身坐标值就发生了变化。当然这种投影变换工作也可以在 ArcMap中通过改变 Data Frame的Coordinate System来实现,

45、 只是要在做完之后按照 Data Frame的坐标系统导出数据即可,即为做动态投影变换”相关的实验可选用 ArcTutor(ArcDesktop 中的练习数据)中的 WGS84和NAD1927 或NAD1983数据进行试验。三、World files 文件1、CAD world files该文件在ArcGIS中对于CAD来说,是定义了 2个点之间的相似变换的坐标对。即源点、目标点坐标值,一般来说源数据点的坐标是CAD图层上一个已知控制点(也可以是任意点)的坐标,而目标点的坐标值为想要该已知点重新定位到的坐标值。CAD的world files文件是一个以wld为后缀的文本文件。因此我们可以用记事

46、本 或写字板任意创建或编辑 world files文件,从而对CAD数据进行正确的配准和矫正。 world file文件的示例如图15所示:M Pdicfrh wld , HoiepddFjlo Ffiniiat 出dw 阴口26315.29965 0,11063.15 3340 519481. 78700, W276. 54 700028854.942990,11110.990860 5195O5.855OO,93502,95800图15 world file文件示例CAD world file里的坐标对最多只有四对,即做一个 2点的变换,world file文件 既可以通过我们自己用文件编

47、辑工具进行创建,同样 ArcGIS也支持world file文件的 创建与保存和加载等操作。在ArcMap中可以借助 Georeferencing 工具条的Add Contorl Point来增加控制点的映射关系(源和目标点),通过 View Link Table打开Link Table对话框,利用Save工具按钮可以将建立好的Link关系保存成 world file文件,或是使用Update Georeferencing 菜单来保存 Link关系从而生成 world file文件。一个CAD只能有一个 world file文件,如果该数据集已经有了 world file文件,则 上述的Save或是Update Georeferencing 操作将重写现有的 world file文件。ArcMap在加载CAD数据时,如果该数据有 world file文件,ArcMap将使用worl d file所定义的转换关系对整个数据集的数据做变换。2、Raster world files影像数据的投影及坐标系统原理与矢量数据一致,只是因为其存储方式是以栅格行列形式存储的,因此需要将栅格数据与现实世界的坐标系建立起一个联系,因此部分影像数据也具有world file文件,

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论