GIS学习笔记——地图投影与GPS参数计算要点_第1页
GIS学习笔记——地图投影与GPS参数计算要点_第2页
GIS学习笔记——地图投影与GPS参数计算要点_第3页
GIS学习笔记——地图投影与GPS参数计算要点_第4页
GIS学习笔记——地图投影与GPS参数计算要点_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

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

文档简介

GIS学习笔记地图投影与GPS参数计算1、椭球体 GIS中的坐标系定义由基准面和地图投影两组参数确定,而基准面的定义则由特定椭球体及其对应的转换参数确定。 基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的基准面。基准面是在椭球体基础上建立的,椭球体可以对应多个基准面,而基准面只一个椭球体。 椭球体的几何定义: O是椭球中心,NS为旋转轴,a为长半轴,b为短半轴。 子午圈:包含旋转轴的平面与椭球面相截所得的椭圆。 纬圈:垂直于旋转轴的平面与椭球面相截所得的圆,也叫平行圈。 赤道:通过椭球中心的平行圈。 基本几何参数: 其中a、b称为长度元素;扁率反映了椭球体的扁平程度。偏心率e和e是子午椭圆的焦点离开中心的距离与椭圆半径之比,它们也反映椭球体的扁平程度,偏心率愈大,椭球愈扁。 套用不同的椭球体,同一个地点会测量到不同的经纬度。下面是几种常见的椭球体及参数列表。 几种常见的椭球体参数值2、地图投影 地球是一个球体,球面上的位置,是以经纬度来表示,我们把它称为“球面坐标系統”或“地理坐标系統”。在球面上计算角度距离十分麻烦,而且地图是印刷在平面纸张上,要将球面上的物体画到紙上,就必须展平,这种将球面转化为平面的过程,称为“投影”。经由投影的过程,把球面坐标换算为平面直角坐标,便于印刷与计算角度与距离。由于球面無法百分之百展为平面而不变形,所以除了地球仪外,所有地图都有某些程度的变形,有些可保持面积不变,有些可保持方位不变,视其用途而定。 目前国际间普遍采用的一种投影,是即横轴墨卡托投影(Transverse Mecator Projection),又称为高斯-克吕格投影(Gauss-Kruger Projection),在小范围内保持形状不变,对于各种应用较为方便。我们可以想象成将一个圆柱体橫躺,套在地球外面,再将地表投影到这个圆柱上,然后将圆柱体展开成平面。圆柱与地球沿南北经线方向相切,我们将这条切线称为“中央经线”。 在中央经线上,投影面与地球完全密合,因此图形没有变形;由中央经线往東西两侧延伸,地表图形会被逐渐放大,变形也会越来越严重。 为了保持投影精度在可接受范围内,每次只能取中央经线两侧附近地区来用,因此必须切割为许多投影带。就像将地球沿南北子午线方向,如切西瓜一般,切割为若干带状,再展成平面。目前世界各国军用地图所采用之 UTM 坐标系統 (Universal Transverse Mecator Projection System),即为横轴投影的一种。是将地球沿子午线方向,每隔 6 度切割为一带,全球共切割为 60 个投影带 地图投影几何分类主要包括: 结合变形性质和几何投影,投影分类包括: 3、GIS中地图投影的定义 我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯克吕格投影(Gauss-Kruger);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(Lambert Conformal Conic);海上小于50万的地形图多用正轴等角园柱投影,又叫墨卡托投影(Mercator),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统。 相应高斯-克吕格投影、兰勃特投影、墨卡托投影需要定义的坐标系参数序列如下: 高斯克吕格:投影代号(Type),基准面(Datum),单位(Unit), 中央经度(OriginLongitude),原点纬度(OriginLatitude), 比例系数(ScaleFactor), 东伪偏移(FalseEasting),北纬偏移(FalseNorthing)兰勃特:投影代号(Type),基准面(Datum),单位(Unit),中央经度(OriginLongitude),原点纬度(OriginLatitude), 标准纬度1(StandardParallelOne),标准纬度2(StandardParallelTwo), 东伪偏移(FalseEasting),北纬偏移(FalseNorthing)墨卡托:投影代号(Type),基准面(Datum),单位(Unit), 原点经度(OriginLongitude),原点纬度(OriginLatitude), 标准纬度(StandardParallelOne) 在城市GIS系统中均采用6度或3度分带的高斯克吕格投影,因为一般城建坐标采用的是6度或3度分带的高斯克吕格投影坐标。高斯克吕格投影以6度或3度分带,每一个分带构成一个独立的平面直角坐标网,投影带中央经线投影后的直线为X轴(纵轴,纬度方向),赤道投影后为Y轴(横轴,经度方向),为了防止经度方向的坐标出现负值,规定每带的中央经线西移500公里,即东伪偏移值为500公里,由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,因此规定在横轴坐标前加上带号,如(,)其中21即为带号,同样所定义的东伪偏移值也需要加上带号,如21带的东伪偏移值为米。假如你的工作区位于21带,即经度在120度至126度范围,该带的中央经度为123度,采用Pulkovo 1942基准面,那么定义6度分带的高斯-克吕格投影坐标系参数为:(8,1001,7,123,0,1,0)。 4、大地坐标系 有了椭球体以及地图投影,坐标系就能确定下来了。北京54和西安80是我们使用最多的坐标系。我们通常称谓的北京54坐标系、西安80坐标系实际上使用的是我国的两个大地基准面北京54基准面和西安80基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系西安80坐标系,目前大地测量基本上仍以北京54坐标系作为参照,北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。 WGS84坐标系采用WGS1984基准面及WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心,目前GPS测量数据多以WGS1984为基准。北京54坐标系北京54坐标系为参心大地坐标系,大地上的一点可用经度L54、纬度M54和大地高H54定位,它是以格拉索夫斯基椭球为基础,经局部平差后产生的坐标系,与苏联1942年建立的以普尔科夫天文台为原点的大地坐标系统相联系,相应的椭球为克拉索夫斯基椭球。到20世纪80年代初,我国已基本完成了天文大地测量,经计算表明,54坐标系统普遍低于我国的大地水准面,平均误差为29米左右。 西安80坐标系 西安80是为了进行全国天文大地网整体平差而建立的。根据椭球定位的基本原理,在建立西安80坐标系时有以下先决条件: (1)大地原点在我国中部,具体地点是陕西省径阳县永乐镇; (2)西安80坐标系是参心坐标系,椭球短轴Z轴平行于地球质心指向地极原点方向,大地起始子午面平行于格林尼治平均天文台子午面;X轴在大地起始子午面内与 Z轴垂直指向经度 0方向;Y轴与 Z、X轴成右手坐标系; (3)椭球参数采用IUG 1975年大会推荐的参数,因而可得西安80椭球两个最常用的几何参数为: 长轴:5(m); 扁率:1:298.257椭球定位时按我国范围内高程异常值平方和最小为原则求解参数。(4)多点定位;(5)大地高程以1956年青岛验潮站求出的黄海平均水面为基准。 WGS84坐标系WGS84(World Geodetic System,1984年)是美国国防部研制确定的大地坐标系,其坐标系的几何定义是:原点在地球质心,z轴指向 BIH 1984.0定义的协议地球极(CTP)方向,X轴指向 BIH 1984.0 的零子午面和 CTP赤道的交点。Y轴与 Z、X轴构成右手坐标系(如图所示)。 WGs-84椭球及有关常数: 对应于WGS-8大地坐标系有一个WGS-84椭球,其常数采用IUGG第17届大会大地测量常数的推荐值。 WGS-84椭球的几何常数: 长半轴: 2(m) 扁率:1 / 298. 地球引力常数(含大气层)GM 正常化二阶带谐系数C2.0 -484.1668510-6 地球自转角速度 10-11 rads -1 主要几何和物理常数 短半径 .3142 扁率 1/298. 第一偏心率平方 e2 0.013 第二偏心率平方 e2 0.2227 橢球正常重力位 U0 .8497m2s-2赤道正常重力 0 9.ms-25.几个名词解释将球面坐标转化为平面坐标的过程便称为投影 大地坐标(Geodetic Coordinate):大地测量中以参考椭球面为基准面的坐标。地面点P的位置用大地经度L、大地纬度B和大地高H表示。当点在参考椭球面上时,仅用大地经度和大地纬度表示。大地经度是通过该点的大地子午面与起始大地子午面之间的夹角,大地纬度是通过该点的法线与赤道面的夹角,大地高是地面点沿法线到参考椭球面的距离。 方里网:是由平行于投影坐标轴的两组平行线所构成的方格网。因为是每隔整公里绘出坐标纵线和坐标横线,所以称之为方里网,由于方里线同时 又是平行于直角坐标轴的坐标网线,故又称直角坐标网。 在1:1万1:20万比例尺的地形图上,经纬线只以图廓线的形式直接表现出来,并在图角处注出相应度数。为了在用图时加密成 网,在内外图廓间还绘有加密经纬网的加密分划短线(图式中称“分度带”),必要时对应短线相连就可以构成加密的经纬线网。1:2 5万地形图上,除内图廓上绘有经纬网的加密分划外,图内还有加密用的十字线。 我国的1:50万1:100万地形图,在图面上直接绘出经纬线网,内图廓上也有供加密经纬线网的加密分划短线。 直角坐标网的坐标系以中央经线投影后的直线为X轴,以赤道投影后的直线为Y轴,它们的交点为坐标原点。这样,坐标系中就出现了四 个象限。纵坐标从赤道算起向北为正、向南为负;横坐标从中央经线算起,向东为正、向西为负。 虽然我们可以认为方里网是直角坐标,大地坐标就是球面坐标。但是我们在一副地形图上经常见到方里网和经纬度网,我们很习惯的称经纬度网为大地坐标,这个时候的大地坐标不是球面坐标,她与方里网的投影是一样的(一般为高斯投影),也是平面坐标。Albers:等面积割圆锥投影,一般涉及到面积的图用它,比如土地利用图等。Lambert(Lambert Conformal Conic):等角割圆锥投影,小比例尺地形图、政区图等常用。在需要计算面积时,可以转到albers等面积投影坐标系下计算。中国全境Xian_80 投影参数如下:中国全境Beijing54 投影参数如下:投影参数如下:Projection: Albers False_Easting: 0.False_Northing: 0.Central_Meridian: 105.Standard_Parallel_1: 25.Standard_Parallel_2: 47.Latitude_Of_Origin: 0.Linear Unit: Meter (1.)(GCS_Xian_1980 Datum: D_Xian_1980)地理坐标参数如下:Geographic Coordinate System: GCS_Xian_1980Angular Unit: Degree (0.)Prime Meridian: Greenwich (0.)Datum: D_Xian_1980Spheroid: Xian_1980 Semimajor Axis: . Semiminor Axis: . Inverse Flattening: 298.投影方式东偏移南北偏移中央经线第一纬线第二纬线起点纬度线性单位坐标名称度转弧度本初子午线大地原点椭球名称A长半轴B短半轴偏心率(扁率(a-b)/a)投影参数如下:Projection: AlbersFalse_Easting: 0.False_Northing: 0.Central_Meridian: 105.Standard_Parallel_1: 25.Standard_Parallel_2: 47.Latitude_Of_Origin: 0.Linear Unit: Meter (1.)(GCS_Beijing_1954,DATUM:D_Beijing_1954)地理坐标参数如下:Name: GCS_Beijing_1954Angular Unit: Degree (0.)Prime Meridian: Greenwich (0.)Datum: D_Beijing_1954 Spheroid: Krasovsky_1940 Semimajor Axis: . Semiminor Axis: . Inverse Flattening: 298.6.我国分省地图投影标准纬线正轴圆锥投影和圆柱投影最适宜于沿纬线伸展的地区,特别是正轴圆锥投影适宜于中纬度地区,正轴圆柱投影最适宜于低纬度和赤道地区。对于沿经线伸展的地区,宜采用横轴圆柱投影。6.1中国分省(区)地图投影的选择 (1)从制图区域的形状和位置来看:我国绝大多数省(区)处于中纬度地区,因此最适宜采用圆锥投影;对于个别省区,如广东省包括南海诸岛及南中国海域,它位于赤道附近地区,可采用正轴圆柱投影;对于经差较小的地区,亦可采用高斯克吕格投影。即正轴等角圆锥投影;正轴等角割圆柱投影;宽带高斯克吕格投影。我国目前各省(区)按制图区域单幅地图选择投影时,所采用的两条标准纬线如下:注:北京市、天津市标准纬线同河北省,上海市标准纬线同江苏省。 南海诸岛采用正圆柱投影。另一种情况,是采用分带投影的方法,即把相近的同纬度省(区)合用一个投影,把全国各省(区)分别采用若干个正轴等角圆锥投影,下表是将全国各省(区)分为10个投影带,计算得采用正轴等角圆锥投影时长度变形小于0.5%, 中国全图(南海诸岛作插图)正轴等面积割圆锥投影目前常采用1=25 00,2=47 00在GIS中,地理数据的显示往往可以根据用户的需要,指定各种投影。但当所显示的地图与国家基本地图系列的比例尺一致时,往往采用与国家基本系列地图所用的投影。我国常用的地图投影的情况为: (1)、我国基本比例尺地形图(1:100万、1:50万、1:25万、1:10万、1:5万、1:2.5、1:1万、1:5000)除1:100万外均采用高斯克吕格投影为地理基础; (2)、我国1:100万地形图采用了Lambert投影,其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影保持一致。 (3)、我国大部分省区图以及大多数这一比例尺的地图也多采用Lambert投影和属于同一投影系统的Albers投影(正轴等面积割圆锥投影);(4)、Lambert投影中,地球表面上两点间的最短距离(即大圆航线)表现为近于直线,这有利于地理信息系统中和空间分析量度的正确实施。Lambert投影一种等角圆锥投影,我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),又叫横轴墨卡托投影(Transverse Mercator);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(Lambert Conformal Conic),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统。我国的全国地图及分省地图和小于1比50万的地图一般使用lambert投影,全国地图的标准纬线现在是使用 25度和47度(之前使用过 25,45)。而我国1:100万地图的投影是按百万分之一地图的纬度划分原则分带投影的。即从0开始,每隔纬差4为一个投影带,每个投影带单独计算坐标,建立数学基础。同一投影带内再按经差6分幅,各图幅的大小完全相同,故只需计算经差6、纬差4的一幅图的投影坐标即可。每幅图的直角坐标,是以图幅的中央经线作为轴,中央经线与图幅南纬线交点为原点,过原点切线为轴,组成直角坐标系。每个投影带设置两条标准纬线,其位置是:1=S+30´2=N-30´该投影的变形分布规律:没有角度变形;两条标准纬线上没有任何变形;由于采用了分带投影,每带纬差较小,因此我国范围内的变形几乎相等,最大长度变形不超过0.03%(南北图廓和中间纬线),最大面积变形不大于0.06。7.GPS参数DXDYDZ三个参数代表的是坐标系上三个坐标轴的平移参量,da代表的是需要转换的两个坐标系之间椭球半长轴的差,df代表的是需要转换的两个坐标系之间椭球扁率倒数的差。广东省转换参数地区名称地区代码中心经度地区名称地区代码中心经度广州114梅州117韶关114汕尾114深圳114河源114珠海114阳江111汕头117清远114佛山114东莞114江门114中山114湛江111潮州117茂名111揭阳117肇庆111云浮111惠州114WGS-84转北京54坐标参数:序号DXDYDZDADF适 宜 地 区1-14-110-52-1080.广州、佛山、江门、中山、珠海、肇庆、云浮、顺德2-8-103-61-1080.韶关、清远3-9-98-71-1080.茂名、阳江、湛江、4-19-112-55-1080.惠州、深圳、东莞、河源5-30-119-58-1080.揭阳、潮州、汕头、汕尾6-16-114-53-1080.梅州WGS-84转西安80坐标参数:DXDYDZDADF适宜地区-96-5112-30广东省8.手持GPS坐标系统转换参数的求解方法目前,市面上出售的手持GPS所使用的坐标系统基本都是WGS-84坐标系统,而我们使用的地图资料大部分都属于1954年北京坐标系或1980年西安国家大地坐标系。不同的坐标系统给我们的使用带来了困难,于是就出现了如何把WGS-84坐标转换到1954北京坐标系或1980西安国家大地坐标系上来的问题。大家知道,不同坐标系之间存在着平移和旋转的关系,要使手持GPS所测量的数据转换为自己需要的坐标,必须求出两个坐标系(WGS-84和北京54坐标系或西安80坐标系)之间的转换参数。两坐标系之间的转换有七参数法、五参数法、和三参数法。七参数法一般用于转换精度要求较高的计算,而手持GPS接收机内部设置的是五参数法,因此只要用户计算出五个参数(DX、DY、DZ、DA、DF)并按提示输入即可在仪器上进行坐标转换。下面以1954北京坐标系为例,求手持GPS接收机坐标转换五个参数的方法。计算流程见下图:8.1收集测区高等级控制点资料在应用手持GPS接收机观测的区域内找出三个以上分布均匀的等级点(精度越高越好)或GPS“B”级网网点,点位最好是周围无电磁波干扰,视野开阔,卫星信号强。并到测绘管理部门抄取这些点的54北京坐标系的高斯平面直角坐标(x、y),大地经纬度(B、L),高程h ,高程异常值和WGS-84坐标系的大地经纬度(B、L),大地高H。8.2直接用手持GPS测定已知点B、L、H值若同时收集到北京坐标系x、y、B、L、h、值和WGS-84坐标系B、L、H值,则不需此步骤。用户如果收集到的只是54北京坐标,必须进行此步工作。由于WGS-84坐标系与我国坐标系之间的平面差异较大,要消除这个误差,应借助收集到的控制点坐标进行转换参数的计算,此时应在收集到的高等级控制点上分别测量B、L、H值(即WGS-84坐标),供计算转换参数时使用。 8.3计算三维直角坐标X、Y、Z 大地坐标系与空间直角坐标系之间的关系如下图。对于同一空间点,大地坐标系与空间直角坐标系有下列转换关系式: (3-1)其中,N=A/(1-E2sin2B)1/2,1954北京坐标系的大地高H=h+,X、Y、Z为大地坐标系中的三维直角坐标,A为大地坐标系对应椭球之长半轴,E为大地坐标系对应第一偏心率,F为大地坐标系对应之扁率,N为该点的卯酉圈曲率半径。根据测量到的大地坐标值BWGS84、L WGS84、H WGS84和收集到的BBJ54、LBJ54、HBJ54分别代入公式(3-1)求得三维直角坐标X1、Y1、Z1和X2、Y2、Z2。如果收集到的只是高斯平面直角坐标(x,y),则应把平面直角坐标(x,y)代入高斯投影反算公式(3-2)求出大地坐标值(B,L)再代入式(3-1)求X2、Y2、Z2,此时(3-2)式中不同坐标系对应的椭球参数见下表:项目WGS-84北京54西安80长半轴A M M M第一偏心率平方E20.0130.0.扁率F1/298.1/298.31/298.2578.4求五个转换参数DX、DY、DZ、DA、DF(4-1)上述函数模型是把WGS-84坐标系的空间直角坐标原点平移到1954北京坐标系的原点,用1954北京坐标系的椭球参数重算以达到两坐标系转换的方法。(注:手持GPS存储的是WGS-84大地坐标,转换到用户坐标系后显示的是经过转换后的坐标值)把第三步求得的两组三维直角坐标值分别代入(4-1)式,求出DX、DY、DZ、DA、DF值。一旦需转换的两个坐标系统确定以后,DA、DF是常值,但是DX、DY、DZ对于不同地区有不同的值。8.5参数检验DX、DY、DZ、DA、DF五个转换参数求出后,必须按提示分别输入手持GPS,同时输入测区中央子午线经度。E代表东经,投影比例参数为1.00000,东西偏差为m,南北偏差为0,并设单位为米。输入这些参数后,应拿到实地检测,检验这五个参数是否正确。方法是,在野外选定视野开阔、GPS接收信号强的特征点(如线状地物交叉点、独立地物等),最好是埋石控制点进行测量,然后找出这些点的理论坐标与之比较。如比较结果超过仪器标称精度,则应重新测算转换参数。下面是我们在河池市国际疫苗研究所“DOMI伤寒Vi疫苗国际合作项目研究” 项目中,测量GIS基础数据图时,对手持GPS进行坐标转换后,进行实地检测的结果。 检测坐标比较表:特征点实测Y实测X已知Y已知XYX130511948305219461.0-2.02.2230441956304419540.0-2.02.033039195230381949-1.0-3.03.2430451942304719422.00.02.0531962058320020534.0-5.06.4630012232310322292.0-3.03.6731392196314221933.0-3.04.2从上表可以看出,最大点位误差6.4m,最小点位误差2.0m,点位误差并不大,以目前市面上出售的手持GPS(标称的单机定位精度一般为10-15米)来说,在没有进行差分的情况下,达到这种精度,应该说还是比较高的。而我们上述转换参数的求解结果仅是在图根点上进行,如果在GPS B级网点上进行求解,相信精度会更高。使用ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标的转换1 坐标转换简介坐标系统之间的坐标转换既包括不同的参心坐标之间的转换,或者不同的地心坐标系之间的转换,也包括参心坐标系与地心坐标系之间的转换以及相同坐标系的直角坐标与大地坐标之间的坐标转换,还有大地坐标与高斯平面坐标之间的转换。 在两个空间角直坐标系中,假设其分别为O-XYZ和O-XYZ,如果两个坐标系的原来相同,通过三次旋转,就可以两个坐标系重合;如果两个直角坐标系的原点不在同一个位置,通过坐标轴的平移和旋转可以取得一致;如果两个坐标系的尺度也不尽一致,就需要再增加一个尺度变化参数;而对于大地坐标和高斯投影平面坐标之间的转换,则需要通过高斯投影正算和高斯投影反算,通过使用中央子午线的经度和不同的参考椭球以及不同的投影面的选择来实现坐标的转换。如何使用ArcGIS实现WGS84经纬度坐标到BJ54高斯投影坐标的转换?这是很多从事GIS工作或者测绘工作者普遍遇到的问题。本文目的在于帮助用户解决这个问题。我们通常说的WGS-84坐标是指经纬度这种坐标表示方法,北京54坐标通常是指经过高斯投影的平面直角坐标这种坐标表示方法。为什么要进行坐标转换?我们先来看两组参数,如表1所示:表1 BJ54与WGS84基准参数参考椭球体长半轴短半轴扁率BJ54基准参数Krasovsky_1940.0188298.3WGS84基准参数WGS 84.3142298.很显然,WGS84与BJ54是两种不同的大地基准面,不同的参考椭球体,因而两种地图下,同一个点的坐标是不同的,无论是三度带六度带坐标还是经纬度坐标都是不同的。当要把GPS接收到的点(WGS84坐标系统的)叠加到BJ54坐标系统的底图上,那就会发现这些GPS点不能准确的在它该在的地方,即“与实际地点发生了偏移”。这就要求把这些GPS点从WGS84的坐标系统转换成BJ54的坐标系统了。有关WGS84与BJ54的坐标转换问题,实质是WGS-84椭球体到BJ54椭球体的转换问题。如果我们是需要把WGS84的经纬度坐标转换成BJ54的高斯投影坐标,那就还会涉及到投影变换问题。因此,这个转换过程,一般的GPS数据处理软件都是采用下述步骤进行的:1)(B,L)84(X,Y,Z)84,空间大地坐标到空间直角坐标的转换。2)(X,Y,Z)84(X,Y,Z)54,坐标基准的转换,即Datum转换。通常有三种转换方法:七参数、简化三参数、Molodensky。3)(X,Y,Z)54(B,L)54,空间直角坐标到空间大地坐标的转换。4)(B,L)54(x,y)54, 高斯投影正算。从以上步骤不难看出,转换的关键是第二步,转换的参数。鉴于我国曾使用不同的坐标基准(BJ54、State80、Correct54),各地的重力值又有很大差异,所以很难确定一套适合全国且精度较好的转换参数。在WGS-84坐标和北京54坐标之间是不存在一套转换参数可以全国通用的,在每个地方会不一样。必须了解,在不同的椭球之间的转换是不严密的。那么,两个椭球间的坐标转换应该是怎样的呢?一般而言比较严密的是用七参数法,即3个平移因子(X平移,Y平移,Z平移),3个旋转因子(X旋转,Y旋转,Z旋转),一个比例因子(也叫尺度变化K)。国内参数来源的途径不多,一般当地测绘部门会有。通行的做法是:在工作区内找三个以上的已知点,利用已知点的BJ54坐标和所测WGS84坐标,通过一定的数学模型,求解七参数。若多选几个已知点,通过平差的方法可以获得较好的精度。如果区域范围不大,最远点间的距离不大于30Km(经验值),这可以用三参数,即只考虑3个平移因子(X平移,Y平移,Z平移),而将旋转因子及比例因子(X旋转,Y旋转,Z旋转,尺度变化K)都视为0,所以三参数只是七参数的一种特例。北京54和西安80也是两种不同的大地基准面,不同的参考椭球体,他们之间的转换也是同理。在ArcGIS中提供了三参数、七参数转换法。而在同一个椭球里的转换都是严密的,在同一个椭球的不同坐标系中转换需要用到四参数转换,举个例子,在深圳既有北京54坐标又有深圳坐标,在这两种坐标之间转换就用到四参数,计算四参数需要两个已知点。2 ArcGIS坐标转换例子2.1 应注意问题使用ArcGIS如何实现WGS84经纬度坐标到BJ54高斯投影坐标的转换呢?在ArcGIS中,这个坐标转换步骤简化了,用户只需要两个步骤就能够直接从最初的WGS84经纬度坐标转换到BJ54高斯投影坐标。这就是ArcGIS的强大之处。接下来,我们做一个例子。假设我们已经知道了7参数,应该如何操作呢?在具体的操作前,请大家一定注意以下三点:WGS84的经纬度坐标值是用度来表示,而不能是度分秒表示七参数的平移因子单位是米,旋转因子单位是秒,比例因子单位是百万。在ArcGIS中,7参数法的名字是Coordinate_Frame 方法。有人在用ArcGIS进行不同椭球体间的坐标转换时,转换出来的结果不对,然后就写文章说变形如何如何,很可能是由于他们没有注意上面这三个关键的问题造成的。2.2 转换步骤a、定义7参数的地理转换(Create Custom Geographic Transformation)在Arctool中打开Create Custom Geographic Transformation工具,如图1所示:在弹出的窗口中,输入一个转换的名字,如wgs84ToBJ54。在定义地理转换方法下面,在Method中选择合适的转换方法如 COORDINATE_FRAM

温馨提示

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

评论

0/150

提交评论