关于MAPGIS下的投影问题

发布于:2021-07-26 15:03:15

关于 MAPGIS 下的投影问题
默认分类 2009-09-06 17:26:31 阅读 204 评论 2 字号:大中小

投影变换的知识 1 地图投影,是 GIS 知识体系中重要的组成部分,每个 GIS 软件都会涉及到这一 部分知识,并不是只有 MAPGIS 软件中才有,MAPGIS 软件中的投影变换相比 国外的软件更具有针对性,更符合我们国家的国情,比如标准框等我这里只是给 大家说说我对投影变换的一个理解,讲很多的知识点串起来,不正确的地方,还 请大家给予批评指正 投影变换,我个人理解,就是对投影进行变换只要把握住了这个核心的思想,其 他的就不在那么难理解了那么下面就要搞清楚两个问题,就是什么是投影?为什 么要进行投影?然后再来理解如何变换 那么什么是投影呢? 我们知道,地球是一个*似于梨型的不规则椭球体,而 GIS 软件所处理的都是 二维*面上的地物要素的信息所以首先要考的一个问题,就是如果如何将地球表 面上的地物展到*面上去 最简单的一个方法,或者说是最容易想到的一个方法就是将地球表面沿着某个经 线剪开,然后展成*面,即采用这种物理的方法来实现可采用物理的方法将地球 表面展开成地图*面必然产生裂隙或褶皱,大家可以想象一下,如果把一个足球 展成*面的,会是什么结果所以这种方法存在着很大的误差和变形,是不行的 那么我们就可以采用地图投影的方法,就是建立地球表面上的点与地图*面上点 之间的一一对应关系,利用数学法则把地球表面上的经纬线网表示到*面上,这 样就可以很好的控制变形和误差凡是地理信息系统就必然要考虑到地图投影,地 图投影的使用保证了空间信息在地域上的联系和完整性,在各类地理信息系统的 建立过程中,选择适当的地图投影系统是首先要考虑的问题 所以一句话,投影:就是建立地球表面上点(Q,)和*面上的点(x,y)之 间的函数关系式的过程 这时候就有一个问题要问了,就是随着地图制图理论及科学技术的不断发展,就 会有不同的国家,不同的人,提出了不同的数学法则这就表示存在着很多的投影 方式有时候我们需要将不同的投影方式变换成同一种投影方式,或者将不同的投 影参数,变换成相同的投影参数,这都需要进行投影变换 所以一句话,投影变换:就是将不同的地图投影函数关系式变换的过程 在 MAPGIS 中的投影变换的定义如下:将当前地图投影坐标转换为另一种投影 坐标,它包括坐标系的转换不同投影系之间的变换以及同一投影系下不同坐标的 变换等多种变换 下面我们就来看看投影和变换过程中所涉及到的知识点 地球椭球体 地图投影是指建立地球表面上点(Q,)和*面上的点(x,y)之间的函数关 系式的过程那我们先来看看,如何在地球表面上表示地物要素的空间信息只有先 将地球表面上的地物要素的空间信息描述好了以后,在将它们通过函数关系式,

投影到地图*面上去,这样才可以进行空间分析或者其它的运算 我们知道:如果要描述地物要素的空间信息,或者不同地物要素之间的相对空间 关系,首先要在地球上建立一个参考系,只有建立了参考系,才能去准确的描述 每个地物的坐标等信息这涉及到很多地球的形状及椭球体方面的知识 1 地球的形状 地球自然表面是一个起伏不*十分不规则的表面,有高山丘陵和*原,又有江河 湖海地球表面约有 71%的面积为海洋所占用,29%的面积是大陆与岛屿陆地上 最高点与海洋中最深处相差* 20 公里这个高低不*的表面无法用数学公式表 达,也无法进行运算所以在量测与制图时,必须找一个规则的曲面来代替地球的 自然表面当海洋静止时,它的自由水面必定与该面上各点的重力方向(铅垂线方 向)成正交,我们把这个面叫做水准面但水准面有无数多个,其中有一个与静止 的*均海水面相重合可以设想这个静止的*均海水面穿过大陆和岛屿形成一个 闭合的曲面,这就是大地水准面 大地水准面所包围的形体,叫大地球体由于地球体内部质量分布的不均匀,引起 重力方向的变化,导致处处和重力方向成正交的大地水准面成为一个不规则的, 仍然是不能用数学表达的曲面大地水准面形状虽然十分复杂,但从整体来看,起 伏是微小的它是一个很接*于绕自转轴(短轴)旋转的椭球体所以在测量和制图 中就用旋转椭球来代替大地球体,这个旋转球体通常称地球椭球 2 地球的大小 关于地球椭球体的大小,由于采用不同的资料推算,椭球体的元素值是不同的现 将世界各国常用的地球椭球体的数据列表如下: 各种地球椭球体模型 椭球体名称 年代 长半轴(米) 短半轴(米) 扁率 白塞尔(Bessel) 1841 6377397 6356079 1:299.15 克拉克(Clarke) 1880 6378249 6356515 1:293.5 克拉克(Clarke) 1866 6378206 6356584 1:295.0 海福特(Hayford) 1910 6378388 6356912 1:297 克拉索夫斯基(北京 54)1940 6378245 6356863 1:298.3 I.U.G.G (西安 80) 1975 6378140 6356755 1:298.25 WGS-84 1979 6378137 6356759 1:298.26 3MAPGIS 中的椭球体 在 MAPGIS 软件中最常用的就是两种椭球体,它们在 MAPGIS 软件中是以选择 北京 54 坐标系或西安 80 坐标系的方式表现出来的比如在做标准框时,系统提 示我们选择椭球体,这时我们要么选择北京 54,要么选择西安 80 或者其他所以 说在 MAPGIS 中,当提到北京 54 坐标系或西安 80 坐标系时,它们所代表的含 义不是大地测量中的大地坐标系,而是指不同的椭球参数,这个一定要搞清楚 下面我们就了解一下我们国家的坐标系 当前我国采用坐标系主要有 1954 年北京坐标系 1980 年西安坐标系新 1954 年 : 北京坐标系 WGS84 坐标系 该坐标系是通过与原苏联 1942 年坐标系联测而建立的解放后,为了建立我国天 文大地网,鉴于当时历史条件,在东北黑龙江边境上同苏联大地网联测,推算出 其坐标作为我国天文大地网的起算数据;随后,通过锁网的大地坐标计算,推算 出北京点的坐标,并定名为 1954 年北京坐标系因此,1954 年北京坐标系是苏 联 1942 年坐标系的延伸,其原点不在北京,而在苏联普尔科沃该坐标系采用克

拉索夫斯基椭球作为参考椭球,高程系统采用正常高,以 1956 年黄海*均海水 面为基准 该坐标系有两个缺陷:因为它是在东北黑龙江边境上同苏联大地网联测,推算出 其坐标作为我国天文大地网的起算数据,所以随着误差的不断累计,到了中国西 部以后,测量的数据必须经过严格修正后,才能达到要求 1954 年北京坐标系采 用克拉索夫斯基椭球作为参考椭球,这一点和其他国家的参考椭球不一致,所以 该坐标系的数据必须经过变换后才可以在国际上得到认可 1980 年西安坐标系 1978 年 4 月召开的全国天文大地网*差会议上决定建立我国新的坐标系 称为 , 1980 年国家大地坐标系其大地原点设在西安西北的永乐镇,简称西安原点椭球 参数选用 1975 年国际大地测量与地球物理联合会第 16 界大会的推荐值简称 IUUG-75 地球椭球参数或 IAG-75 地球椭球 新 1954 年北京坐标系 将全国大地网整体*差的结果整体换算到克拉索夫斯基椭球体上,形成一个新的 坐标系,称为新 1954 年北京坐标系该坐标系与 1980 年国家大地坐标系的轴定 向基准相同,网的点位精度相同 WGS84 坐标系 在 GPS 定位中,定位结果属于 WGS-84 坐标系该坐标系是使用了更高精度的 VLBLSLR 等成果而建立的坐标系原点位于地球质心,Z 轴指向 BIH1984.0 协 议地极(CTP) 不同的投影方式 前面提到,随着地图制图理论及科学技术的不断发展,就会有不同的国家,不同 的人,提出了不同的数学法则这就表示存在着很多的投影方式下面对不同投影方 式做一下归类,详细的资料可以参考有关的书籍 按地图投影的构成方法分: a 几何投影:几何投影源于透视几何学原理,并以几何特征为依据,将地球椭球 面 的经纬网投影到*面上或投影到可以展成*面的圆柱表面和圆锥表面等几何面 上,从而构成:方位投影圆柱投影圆锥投影; 方位投影:以*面作为辅助投影面,使球体与*面相切或相割,将球体表面上的 经纬网投影到*面上构成的一种投影; 圆柱投影:以圆柱表面作为辅助投影面,使球体和圆柱表面相切或相割,将球体 表面上的经纬网投影到圆柱表面上,然后再将圆柱表面展成*面而构成的一种投 影; 圆锥投影:以圆锥表面作为辅助投影面,使球体和圆柱表面相切或相割,将球体 表面上的经纬网投影到圆柱表面上,然后再将圆锥表面展成*面而构成的一种投 影 据球面和投影面的相对部位不同,上述投投影影有可分为:正轴投影横轴投影斜 轴投影; 在圆柱投影中,以正轴和横轴常见;在圆锥投影中以正轴常见; 正横斜轴方位投影 正横斜轴圆柱投影

正横斜轴圆锥投影 正轴投影经纬线形状 b 非几何投影:通过一系列数学解析法,由几何投影演绎产生了非几何投影,它 们并 不借助投影面,而是根据制图的某些特定要求,如考虑制图区域形状等特点,选 用合适的投影条件,用数学解析方法,求出投影公式,确定*面和球面之间点与 点间的函数关系 据经纬线的形状,可将非几何投影分为伪方位投影伪圆柱投影伪圆锥投影多圆锥 投影;(新编地图学 P59) 伪方位投影:在正轴情况下,伪方位投影的纬线仍投影为同心圆,除中央经线 投影成直线外,其余经线均投影成对称于中央经线的曲线,且交 于纬线的共同圆心; 伪圆柱投影:在正轴圆柱投影基础上,规定纬线仍为*行直线,除中央经线投影 成直线外,其余经线均投影成对称于中央经线的曲线; 伪圆锥投影:在圆锥投影基础上,规定纬线仍为同心圆弧,除中央经线仍为直线 外,其余经线则投影成对称于中央经线的曲线; 多圆锥投影:这是一种假想借助多个圆锥表面与球体相切设计而成的投影,纬线 为同轴为同轴圆弧,其圆心位于中央经线上,中央经线为直线,其余经线则投影 成对称与中央经线的曲线; 按地图投影的变形性质分: 等角投影地球表面上无穷小图形投影后仍保持相似或两微分线段所组成的角 度,在投影后仍保持不变,称等角投影,又称正形投影在等角投影中变形椭圆为 不同大小的圆,它满足条件,极值长度比a=b或经纬线夹角=90°和沿经纬 度长度比相圆锥等(m=n) 等积投影地球面上的图形在投影后保持面积不变,称等面积投影投影中变形椭圆 为面积相等而形状不同的椭圆,这满足条件,面积比P=a×b=mnsin=1 任意投影即不具备以上两种投影的,在任意投影中,如果沿某一主方向的长度比 等于1,即a=1或b=1,则这种投影称为等距离投影 前面对投影方式做了大体的分类后,下面讲解一个具体的投影方式:高斯-克吕 格投影 高斯-克吕格投影由德国数学家高斯提出,后经克吕格扩充并推导出计算公式, 故称为高斯-克吕格投影,简称高斯投影为了控制变形,本投影采用分带的思想 6 度带是从 0 度子午线(在英国的格林尼治天文台附*)起,自西向东每隔经 差 6 为一投影带,全球分为 60 带,各带的带号用自然序数 1,2,3,60 表示 即以东经 0-6 为第 1 带,其中央经线为 3E,东经 6-12 为第 2 带,其中央经线 为 9E,其余类推 3 度带,是从东经 1 度 30 分的经线开始,每隔 3 度为一带,全球划分为 120 个 投影带这样的目的就是为了让 6 度分带的中央经线全部和 3 度分带的中央经线 重合,3 度分带的中央经线只有一半和 6 度分带的中央经线重合 在高斯-克吕格投影上,规定以中央经线为 X 轴,赤道为 Y 轴,两轴的交点为 坐标原点 如下图所示:

我国规定 1:1 万 1:2.5 万 1:5 万 1:10 万 1:25 万 1:50 万比例尺地形 图,均采用高斯-克吕格投影 1:2.5 至 1:50 万比例尺地形图采用经差 6 度分 带,1:1 万比例尺地形图采用经差 3 度分带 MAPGIS 中图框的制作 由于图框和投影变换紧密相连,故 MAPGIS 将其放在同一个系统中在 MAPGIS 中生成图框,大家应该用的很多了,这里就不再重复只是将其中用到的一些重要 的知识点做一下归纳和总结 首先搞清楚在 MAPGIS 大小比例尺的分界,如下:它以 1:5000 为界 小于或等于 1:5000 时,小比例尺,图幅为梯形图幅(在后面地图入库的时候, 是选择矩形图幅,还是梯形图幅就看这里),单位为经纬度; 大于 1:5000 时,大比例尺,图幅为矩形图幅,单位为公里值; 根据这个标准,在 MAPGIS 中我将图框分为一下四类: 小比例尺的标准框:在系统标准框菜单下,选择相应的比例命令即可 小比例尺的非标准框:在投影变换菜单下绘制投影经纬网命令 大比例尺的标准框:在系统标准框菜单下,选择相应的比例尺命令,在矩形分幅 方法中选择正方形或者矩形 大比例尺的非标准框:在系统标准框菜单下,选择相应的比例尺命令,在矩形分 幅方法中选择任意矩形分幅 所以总这里可以看出,小比例尺的标准框和小比例尺的非标准框是通过不同的菜 单下不同的命令生成的,而大比例尺的标准框和大比例尺的非标准框则是通过同 一个命令生成的,只是矩形分幅方法不一致而已 北京 54 坐标系转西安 80 坐标系 首先将 MAPGIS *台的工作路径设置为..\北京 54 转西安 80 文件夹下 下面我们来讲解北京 54 坐标系转西安 80 坐标系的转换方法和步骤 一 数据说明 北京 54 坐标系和西安 80 坐标系之间的转换其实是两种不同的椭球参数之间的 转换,一般而言比较严密的是用七参数布尔莎模型,即 X *移,Y *移,Z *移, X 旋转(WX),Y 旋转(WY),Z 旋转(WY),尺度变化(DM)若得七参数 就需要在一个地区提供 3 个以上的公共点坐标对(即北京 54 坐标下 xyz 和西安 80 坐标系下 xyz),可以向地方测绘局获取 二 北京 54 坐标系转西安 80 坐标系的操作步骤 启动投影变换模块,单击文件菜单下打开文件命令,将演示数据演示数据_北京 54.WT 演示数据_北京 54.WL 演示数据_北京 54.WP 打开,如图 1 所示: 图1 1 单击投影转换菜单下 S 坐标系转换命令,系统弹出转换坐标值话框,如图 2 所示: 图2 在输入一栏中,坐标系设置为北京 54 坐标系,单位设置为线类单位-米; 在输出一栏中,坐标系设置为西安 80 坐标系,单位设置为线类单位-米; 在转换方法一栏中,单击公共点操作求系数项; 在输入一栏中,输入北京 54 坐标系下一个公共点的(xyz),如图 2 所示; 在输出一栏中,输入西安 80 坐标系下对应的公共点的(xyz),如图 2 所示;

在窗口右下角,单击输入公共点按钮,右边的数字变为 1,表示输入了一个公共 点对,如图 2 所示; 依照相同的方法,再输入另外的 2 个公共点对; 在转换方法一栏中,单击七参数布尔莎模型项,将右边的转换系数项激活; 单击求转换系数菜单下求转换系数命令 系统根据输入的 3 个公共点对坐标自动 , 计算出 7 个参数,如图 3 所示,将其记录下来; 然后单击确定按钮; 图3 2 单击投影转换菜单下编辑坐标转换参数命令 系统弹出不同地理坐标系转换参 , 数设置对话框,如图 4 所示; 图4 在坐标系选项一栏中,设置各项参数如下: 源坐标系:北京 54 坐标系; 目的坐标系:西安 80 坐标系; 转换方法:七参数布尔莎模型; 长度单位:米; 角度单位:弧度; 然后单击添加项按钮,则在窗口左边的不同椭球间转换列表中将该转换关系列 出; 在窗口下方的参数设置一栏中,将上一步得到的七个参数依次输入到相应的文本 框中,如图 4 所示; 单击修改项按钮,输入转换关系,并单击确定按钮; 接下来就是文件投影的操作过程了 3 单击投影转换菜单下 MAPGIS 投影转换/选转换线文件命令,系统弹出选择文 件对话框,如图 5 所示: 图5 选中待转换的文件演示数据_北京 54.WL,单击确定按钮; 4 设置文件的 Tic 点,在投影变换模块下提供了两种方法:手工设置和文件间拷 贝,这里不作详细的说明; 5 单击投影转换菜单下编辑当前投影参数命令,系统弹出输入投影参数对话框, 如图 6 所示,根据数据的实际情况来设置其地图参数,如下: 坐标系类型:大地坐标系 椭球参数:北京 54 投影类型:高斯-克吕格投影 比例尺分母:1 坐标单位:米 投影中心点经度(DMS):1230000 然后单击确定按钮; 图6 6 单击投影转换菜单下设置转换后参数命令,系统弹出输入投影参数对话框,如 图 7 所示,转换后的参数设置为: 坐标系类型:大地坐标系

椭球参数:西安 80(注意椭球参数的变换) 投影类型:高斯-克吕格投影 比例尺分母:1 坐标单位:米 投影中心点经度(DMS):1230000(注意前后中央经线保持一致) 图7 7 单击投影转换菜单下进行投影变换命令,系统弹出输入转换后位移值对话框, 单击开始转换按钮,系统开始按照设定的参数转换线文件,如图 8 所示: 图8 以同样的操作步骤和参数设置,将演示数据_北京 54.WL 演示数据_北京 54.WP 文件进行投影转换; 8 单击鼠标右键,选择复位命令,系统弹出选择文件名对话框,可以看到系统生 成了三个新的文件:NEWLIN.WLNEWPNT.WTNEWPNT.WP,依次选中这三个 文件,单击确定按钮,如图 7 所示: 这时新生成的三个文件就是西安 80 坐标系下的文件; 补充:通常情况下,转换过来的数据会有一定的误差存在,所以有时为了保证数 据的精度,在转换的过程中通过设置横坐标和纵坐标的偏移量来修正转换后的坐 标值; 跨带投影 我们知道高斯-克吕格投影采用了分带投影的思想,每一个投影带的坐标都是对 本带坐标原点的相对值,所以各带的坐标完全相同,使用时只需变一个带号即 可,这样就存在着如果不考虑带号的情况下,会有重叠的情况出现,如果要想将 重叠的图框回到其实际所在的位置上,这时就会用到跨带投影 跨带投影是 MAPGIS 投影变换中一个很重要的方面 下面来讲解跨带投影的操作方法和步骤,共分为两部分: 一 演示数据的生成和说明: 在投影变换模块下分别生成 3 幅 1:50 万的标准框,并在输入编辑模块中将其 改成不同的颜色(FRAM_50_左.W~表示 FRAM_50_左.WL 和 FRAM_50_ 左.WT): 名称 起始经度(DDMMSS) 起始纬度(DDMMSS) 中央经线(DDMMSS) FRAM_50_左.W~ 1140000 280000 1170000 FRAM_50_中.W~ 1170000 280000 1170000 FRAM_50_右.W~ 1200000 280000 1230000 因为在投影过程中采用的是高斯克吕格投影,且 1:50 万的标准图框的经线跨 度为 3°,所以当同时打开这三幅标准图框时,会发现 FRAM_50_左.W~和 FRAM_50_右.W~二者重叠在了一起,如图 1 所示: 图1 如果想实现三个标准框连续排列,则需要经过跨带投影 二 跨带投影的操作步骤 启动投影变换模块,单击文件菜单下打开文件命令,将 FRAM_50_ 左.W~FRAM_50_中.W~FRAM_50_右.W~三个标准框添加进来,如图 2 所示:

图2 1 单击投影转换菜单下 MAPGIS 文件投影/选转换线文件文件命令,系统弹出选 择文件对话框,选择 FRAM_50_右.WL,单击确定按钮,如图 3 所示: 图3 2 设置文件的 Tic 点,因为在生成标准图框时 MAPGIS 为自动为其添加 4 个 Tic 点,所以这里不再作详细的说明; 3 单击投影转换菜单下编辑当前投影参数命令,系统弹出输入投影参数对话框, 如图 4 所示: 坐标系类型:投影*面直角坐标系 椭球参数:西安 80 投影类型:高斯-克吕格投影 比例尺分母:500000 坐标单位:毫米 投影中心点经度(DMS):1230000 通常情况下,因为是标准框,所以系统会自动的读取其各项参数,所以只需检查 各项参数设置是否有错即可; 图4 4 单击投影转换菜单下设置转换后参数命令,系统弹出输入投影参数对话框,如 图 5 所示: 坐标系类型:投影*面直角坐标系 椭球参数:西安 80 投影类型:高斯-克吕格投影 比例尺分母:500000 坐标单位:毫米 投影中心点经度(DMS):1170000(注意前后中央经线发生了变化) 图5 5 单击投影转换菜单下进行投影变换命令,系统弹出输入转换后位移值对话框, 单击开始转换按钮,系统开始按照设定的参数转换线文件,如图 6 所示: 图6 以同样的操作步骤和参数设置,将 FRAM_50_右.WT 文件进行投影转换; 6 单击鼠标右键,选择复位命令,系统弹出选择文件名对话框,可以看到系统生 成了两个新的文件:NEWLIN.WLNEWPNT.WT,依次选中 FRAM_50_ 左.W~FRAM_50_中.W~及两个新生成的文件,然后单击确定按钮,如图 7 所示: 图7 最终结果如图 8 所示: 图8 补充:中央经线的设置方法 跨带投影的过程中设计到一个很重要的参数就是中央经线,因为高斯-克吕格投 影采用的是分带的思想,所以在每个投影带都会有一个中央经线,中央经线设置 错误,则投影变换的结果就会有问题,尤其是跨带投影的情况下 那如何查阅一个标准框的中央经线呢?

我们国家规定:高斯-克吕格投影,1:2.5 万~1:50 万地形图均采用 6 度分 带;1:1 万及更大比例尺采用 3 度分带,所以上诉 3 幅标准图框都采用的 6 度 分带 由标准框的起始经纬度,如 FRAM_50_左.W~的起始经纬度 1140000,我们可 以查阅出其对应的中央经线 单击投影变换模块帮助菜单下帮助目录命令,在系统弹出的对话框中,选择索引 页,找到 6 度分带表,单击显示按钮,如图 9 所示: 图9 则 6 度分带表如图 10 所示: 6 度分带表 根据标准框的起始经纬度,可以分别查阅到 FRAM_50_左.W~的中央经线为: 1170000 FRAM_50_中.W~的中央经线为:1170000FRAM_50_右.W~的中央经线为: 1230000 地图坐标常识 1 椭球面 地图坐标系由大地基准面和地图投影确定,大地基准面是利用特定椭球体对 特定地区地球表面的*,因此每个国家或地区均有各自的大地基准面,我们通 常称谓的北京 54 坐标系西安 80 坐标系实际上指的是我国的两个大地基准面我 国参照前苏联从 1953 年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的 北京 54 坐标系,1978 年采用国际大地测量协会推荐的 IAG 75 地球椭球体建立 了我国新的大地坐标系--西安 80 坐标系, 目前 GPS 定位所得出的结果都属于 WGS84 坐标系统,WGS84 基准面采用 WGS84 椭球体,它是一地心坐标系, 即以地心作为椭球体中心的坐标系因此相对同一地理位置,不同的大地基准面, 它们的经纬度坐标是有差异的 采用的 3 个椭球体参数如下 (源自全球定位系统测量规范 GB/T 18314-2001) : 理解:椭球面是用来*厍虻,应该是一个立的椭圆旋转而成的 2 大地基准面 椭球体与大地基准面之间的关系是一对多的关系,也就是基准面是在椭球体 基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面, 如前苏联的 Pulkovo 1942 非洲索马里的 Afgooye 基准面都采用了 Krassovsky 椭球体,但它们的大地基准面显然是不同的在目前的 GIS 商用软件中,大地基 准面都通过当地基准面向 WGS84 的转换 7 参数来定义,即三个*移参数 XYZ 表示两坐标原点的*移值;三个旋转参数 xyz 表示当地坐标系旋转至与地心坐标 系*行时,分别绕 XtYtZt 的旋转角;最后是比例校正因子,用于调整椭球大小 北京 54 西安 80 相对 WGS84 的转换参数至今没有公开,实际工作中可利用工 作区内已知的北京 54 或西安 80 坐标控制点进行与 WGS84 坐标值的转换,在 只有一个已知控制点的情况下(往往如此),用已知点的北京 54 与 WGS84 坐标 之差作为*移参数,当工作区范围不大时,如青岛市,精度也足够了 以(32°,121°)的高斯-克吕格投影结果为例,北京 54 及 WGS84 基准面, 两者投影结果在南北方向差距约 63 米(见下表),对于几十或几百万的地图来说, 这一误差无足轻重,但在工程地图中还是应该加以考虑的

输入坐标(度) 北京 54 高斯投影(米) WGS84 高斯投影(米) 纬度值(X) 32 3543664 3543601 经度值(Y) 121 21310994 21310997 理解:椭球面和地球肯定不是完全贴合的,因而,即使用同一个椭球面,不同的 地区由于关心的位置不同,需要最大限度的贴合自己的那一部分,因而大地基准 面就会不同 3 高斯投影 (1)高斯-克吕格投影性质 高斯-克吕格(Gauss-Kruger)投影简称高斯投影,又名"等角横切椭圆柱投影, 地球椭球面和*面间正形投影的一种德国数学家物理学家天文学家高斯(Carl FriedrichGauss,1777 一 1855)于十九世纪二十年代拟定,后经德国大地测 量学家克吕格(Johannes Kruger,1857~1928)于 1912 年对投影公式加以补 充,故名该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线 的条件,确定函数的形式,从而得到高斯一克吕格投影公式投影后,除中央子午 线和赤道为直线外, 其他子午线均为对称于中央子午线的曲线设想用一个椭圆 柱横切于椭球面上投影带的中央子午线,按上述投影条件,将中央子午线两侧一 定经差范围内的椭球面正形投影于椭圆柱面将椭圆柱面沿过南北极的母线剪开 展*,即为高斯投影*面取中央子午线与赤道交点的投影为原点,中央子午线的 投影为纵坐标 x 轴,赤道的投影为横坐标 y 轴,构成高斯克吕格*面直角坐标系 高斯-克吕格投影在长度和面积上变形很小,中央经线无变形,自中央经线向 投影带边缘,变形逐渐增加,变形最大之处在投影带内赤道的两端由于其投影精 度高,变形小,而且计算简便(各投影带坐标一致,只要算出一个带的数据,其 他各带都能应用),因此在大比例尺地形图中应用,可以满足军事上各种需要, 能在图上进行精确的量测计算 (2)高斯-克吕格投影分带 按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形 的最有效方法分带时既要控制长度变形使其不大于测图误差,又要使带数不致过 多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的* 形地带,以便分带投影通常按经差 6 度或 3 度分为六度带或三度带六度带自 0 度 子午线起每隔经差 6 度自西向东分带 带号依次编为第 1260 带三度带是在六度 , 带的基础上分成的,它的中央子午线与六度带的中央子午线和分带子午线重合, 即自 1.5 度子午线起每隔经差 3 度自西向东分带 带号依次编为三度带第 12120 , 带我国的经度范围西起 73°东至 135°,可分成六度带十一个,各带中央经线依 次为 75° 87° 81° 117° 123° 129° 135° , 或三度带二十二个六度带可用于中小比例尺 (如 1:250000)测图,三度带可用于大比例尺(如 1:10000)测图,城建 坐标多采用三度带的高斯投影 (3)高斯-克吕格投影坐标 高斯- 克吕格投影是按分带方法各自进行投影,故各带坐标成独立系统以中 央经线投影为纵轴(x), 赤道投影为横轴(y),两轴交点即为各带的坐标原点纵坐标 以赤道为零起算,赤道以北为正,以南为负我国位于北半球,纵坐标均为正值横 坐标如以中央经线为零起算,中央经线以东为正,以西为负,横坐标出现负值, 使用不便,故规定将坐标纵轴西移 500 公里当作起始轴,凡是带内的横坐标值 均加 500 公里由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点 的相对值,所以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,在横

轴坐标前加上带号,如(4231898m,21655933m),其中 21 即为带号 (4)高斯-克吕格投影与 UTM 投影 某些国外的软件如 ARC/INFO 或国外仪器的配套软件如多波束的数据处理软 件等,往往不支持高斯-克吕格投影,但支持 UTM 投影,因此常有把 UTM 投影 坐标当作高斯-克吕格投影坐标提交的现象 UTM 投影全称为通用横轴墨卡托投影,是等角横轴割圆柱投影 (高斯-克吕格 为等角横轴切圆柱投影),圆柱割地球于南纬 80 度北纬 84 度两条等高圈,该 投影将地球划分为 60 个投影带,每带经差为 6 度,已被许多国家作为地形图的 数学基础 UTM 投影与高斯投影的主要区别在南北格网线的比例系数上,高斯克吕格投影的中央经线投影后保持长度不变,即比例系数为 1,而 UTM 投影的 比例系数为 0.9996UTM 投影沿每一条南北格网线比例系数为常数,在东西方向 则为变数,中心格网线的比例系数为 0.9996,在南北纵行最宽部分的边缘上距 离中心点大约 363 公里,比例系数为 1.00158 高斯-克吕格投影与 UTM 投影可*似采用 Xutm=0.9996 * X 高斯, Yutm=0.9996 * Y 高斯进行坐标转换以下举例说明(基准面为 WGS84): 输入坐标(度) 高斯投影(米) UTM 投影(米) Xutm=0.9996 * X 高斯, Yutm=0.9996 * Y 高斯 纬度值(X) 32 3543600.9 3542183.5 3543600.9*0.9996 3542183.5 经度值 (Y) 121 21310996.8 311072.4 (310996.8-500000)*0.9996+500000 311072.4 注:坐标点(32,121)位于高斯投影的 21 带,高斯投影 Y 值 21310996.8 中前 两位 21 为带号;坐标点(32,121)位于 UTM 投影的 51 带,上表中 UTM 投影 的 Y 值没加带号因坐标纵轴西移了 500000 米,转换时必须将 Y 值减去 500000 乘上比例因子后再加 500000 理解:高斯投影的方法就是保持赤道和中央经线不变形,把球面摊*方法:用一 个椭圆柱套住椭球,把它投影到椭圆柱上,然后打开椭圆柱即可 注:坐标点(32,121)位于高斯投影的 21 带,高斯投影 Y 值 21310996.8 中前 两位 21 为带号;坐标点(32,121)位于 UTM 投影的 51 带,上表中 UTM 投影 的 Y 值没加带号因坐标纵轴西移了 500000 米,转换时必须将 Y 值减去 500000 乘上比例因子后再加 500000 理解:高斯投影的方法就是保持赤道和中央经线不变形,把球面摊*方法:用一 个椭圆柱套住椭球,把它投影到椭圆柱上,然后打开椭圆柱即可 (4)高斯-克吕格投影与 UTM 投影 某些国外的软件如 ARC/INFO 或国外仪器的配套软件如多波束的数据处理软 件等,往往不支持高斯-克吕格投影,但支持 UTM 投影,因此常有把 UTM 投影 坐标当作高斯-克吕格投影坐标提交的现象 UTM 投影全称为通用横轴墨卡托投影,是等角横轴割圆柱投影 (高斯-克吕格 为等角横轴切圆柱投影),圆柱割地球于南纬 80 度北纬 84 度两条等高圈,该 投影将地球划分为 60 个投影带,每带经差为 6 度,已被许多国家作为地形图的 数学基础 UTM 投影与高斯投影的主要区别在南北格网线的比例系数上,高斯克吕格投影的中央经线投影后保持长度不变,即比例系数为 1,而 UTM 投影的

比例系数为 0.9996UTM 投影沿每一条南北格网线比例系数为常数,在东西方向 则为变数,中心格网线的比例系数为 0.9996,在南北纵行最宽部分的边缘上距 离中心点大约 363 公里,比例系数为 1.00158 高斯-克吕格投影与 UTM 投影可*似采用 Xutm=0.9996 * X 高斯, Yutm=0.9996 * Y 高斯进行坐标转换以下举例说明(基准面为 WGS84): 输入坐标(度) 高斯投影(米) UTM 投影(米) Xutm=0.9996 * X 高斯, Yutm=0.9996 * Y 高斯 纬度值(X) 32 3543600.9 3542183.5 3543600.9*0.9996 3542183.5 经度值 (Y) 121 21310996.8 311072.4 (310996.8-500000)*0.9996+500000 311072.4 注:坐标点(32,121)位于高斯投影的 21 带,高斯投影 Y 值 21310996.8 中前 两位 21 为带号;坐标点(32,121)位于 UTM 投影的 51 带,上表中 UTM 投影 的 Y 值没加带号因坐标纵轴西移了 500000 米,转换时必须将 Y 值减去 500000 乘上比例因子后再加 500000 理解:高斯投影的方法就是保持赤道和中央经线不变形,把球面摊*方法:用一 个椭圆柱套住椭球,把它投影到椭圆柱上,然后打开椭圆柱即可 注:坐标点(32,121)位于高斯投影的 21 带,高斯投影 Y 值 21310996.8 中前 两位 21 为带号;坐标点(32,121)位于 UTM 投影的 51 带,上表中 UTM 投影 的 Y 值没加带号因坐标纵轴西移了 500000 米,转换时必须将 Y 值减去 500000 乘上比例因子后再加 500000 理解:高斯投影的方法就是保持赤道和中央经线不变形,把球面摊*方法:用一 个椭圆柱套住椭球,把它投影到椭圆柱上,然后打开椭圆柱即可

这篇文章写得不错,我留下了,对我有帮助,也是尊重原作者.


相关推荐

最新更新

猜你喜欢