地图投影变换详解

《地图投影变换详解》由会员分享，可在线阅读，更多相关《地图投影变换详解（112页珍藏版）》请在金锄头文库上搜索。

1、QQ:83079879 E-mail：地图投影投影变换地图配准预备知识（复习）既然GIS是解决地理位置（未知）怎样从GIS的角度去认识地球？GIS应该是从坐标系统开始的，然后带 领我们一步步向前GIS中的坐标系定义是GIS系统的基础 ，正确定义GIS系统的坐标系非常重要 。GIS中的坐标系定义由基准面和地图 投影两组参数确定，而基准面的定义 则由特定椭球体及其对应的转换参数 确定，因此欲正确定义GIS系统坐标系 ，首先必须弄清地球椭球体(Ellipsoid) 、大地基准面(Datum)及地图投影 (Projection)三者的基本概念及它们之 间的关系。 大地水准面是由静止海水面并向大陆延 伸

2、所形成的不规则的封闭曲面。它是 重力等位面，即物体沿该面运动时， 重力不做功（如水在这个面上是不会 流动的） 地球椭球体只不过是一个具有长半轴， 短半轴和变率的椭球体，可以任意放置 的，它没有为我们规定度量的起点，所 以就有基准面的产生，而基准面就是规 定了度量标准。基准面是利用特定椭球体对特定地区地球表 面的逼近，因此每个国家或地区均有各自的 基准面，我们通常称谓的北京54坐标系、西 安80坐标系实际上指的是我国的两个大地基 准面.椭球体与基准面之间的关系是一对多的 关系，也就是基准面是在椭球体基础上建立 的，但椭球体不能代表基准面，同样的椭球 体能定义不同的基准面，一般意义上基准面 与参考

3、椭球体是同一个概念。地面上选一点P，由P点投影到大地水准面P0点,使P0上的椭球面与大地 水准面相切， 此时过P0点的铅垂线与P0点的椭球面法线重合，切点P0称为 大地原点。同时要使旋转椭球短轴与地球短轴相平行（不要求重合），达 到本国范围内的大地水准面与椭球面十分接近，该椭球面称为参考椭球面 。我国大地原点选在我国中部陕西省泾阳县永乐镇。P P地球表面地球表面大地水准面大地水准面N NS SWWE EP P0 0垂线与法线垂线与法线 重合重合参考椭球体参考椭球体参考椭球面与我国大地原点地球椭球的参数可用a(长半径)、b（短半径）及（扁率）表示。扁率为 1979年国际大地测量与地球物理联合会推

4、荐的地球椭球参数a=6378140m， b=6356755.3m, =1:298.257。旋转椭球面是数学表面，可用如下的公式表示：按一定的规则将旋转椭球与大地体套合在一起，这项工作称椭球定位。定位时采用椭球中心与地球质心重合，椭球短轴与地球短轴重合，椭球与全球大地水准面差距的平方和最小，这样的椭球称总地球椭球。 总地球椭球与参考椭球（水准面）的区别对于地理坐标，只需要确定两 个参数，即椭球体和大地基准 面。 （为什么？）有了基准面我们就可以操作了，我们的一切 都是基于基准面的，为了说明这个，我还看 了一下相关资料。ENVI中的坐标定义文件存放在 HOMEITTIDL70productsenv

5、i45map _proj 文件夹下，三个文件记录了坐标信息 ：ellipse.txt 椭球体参数文件datum.txt 基准面参数文件map_proj.txt 坐标系参数文件在ENVI中自定义坐标系分三步：定义椭球 体、基准面和定义坐标参数第一步、添加椭球体语法为 ,。这里 将“Krassovsky，6378245.0，6356863.0” 因为Envi下面有这个椭球体，我没有加。 第二步、添加基准面语法为,。这里将“Beijing-54, Krassovsky, -12, -113, -41” datum.txt末端。第三步、定义坐标在ENVI任何用到投影坐标的功能模块中 都可以新建坐标系(

6、在任何地图投影选择对话 框中，点击“New”按钮。)，这里我们选择 Map-Customize Map Projection，如图所 示，将相应的参数添加，这里添加的参数如 图所示。注：投影类型选择Transverse Mercator，Scale factor填写0.9996，与 Gauss-Kruger等同。False easting中如果把 带号，即39500000，得到的坐标就带有带 号。填写相关参数后，在map_proj.txt 坐标系参数文 件中看到，完毕。ERDAS中添加椭球体和基准面F:Program FilesLeica GeosystemsGeospatial Imagin

7、g 9.2etcspheroid.tab 记住一个椭球体可以有多个基准面基本语法为： “椭球名称” “椭球序号”椭球体长半轴 椭球体短半轴 “椭球名称” 0 0 0 0 0 0 0 “基准面名称1” dx1 dy1 dz1 rx1 rz1 ds1 “基准面名称2” dx2 dy2 dz2 rx2 rz1 ds2 . 其中：“基准面名称” dx dy dz rx rz ds中，dx、dy、dz是x、y、z3个轴对 于WGS84基准点的平移参数，单位为 m。rx、ry、rz是x、y、z、3个轴对 于WGS84基准点的旋 转参数，单位为 rad。Ds是对于WGS84基准点的比例 因子。 在更多的情况

8、下椭球的基准面是基于它本 身的。这时假定椭球的中心点是与没有 经过任何平移或旋转的WGS84的基准 面相重合，即这时椭球基准面的7个参 数均 为0，即这时椭球基准面的7个参数 均为0。我国在使用克拉索夫斯基椭球 和IAG75椭球时就是用椭球体本身为基 准 在spheroid.tab文件末尾加入如下语 句即可， “IAG 75“ 75 6378140 6356755.2882 “xian 80” 0 0 0 0 0 0 0 在Viewer中打开图像数据，Utility- layer info，在projection info 栏中可 以看到目前的数据投影信息还不完整。 点击edit菜单中的cha

9、nge map model ，在弹出窗口中将unite参数设为 meters，projection，参数设为 Tansverse Mercator。接下来再点击 edit菜单中的Add/Change projection ，在弹出对话框中将原始投影参数添加 进去。 投影变换想想同一基准面的投影变换？不同基准面之间的变换？在ArcGIS Desktop中进行三参数或七参数 精确投影转换 ArcGIS中定义的投影转换方法，在对数据 的空间信息要求较高的工程中往往不能适用 ，有比较明显的偏差。在项目的前期数据准 备工作中，需要进行更加精确的三参数或七 参数投影转换。下面介绍两种办法来在ArcGIS

10、Desktop中进行这种转换。方法1：在ArcMap中进行动态转换 假设原投影坐标系统为Xian80坐标系统，本例选择为系统预设 的Projected Coordinate SystemsGauss KrugerXian 1980Xian 1980 GK Zone 20投影，中央经线为117度，要转 换成Beijing 1954Beijing 1954 GK Zone 20N。在ArcMap中加载了图层之后，打开View-Data Frame Properties对话框，显示当前的投影坐标系统为Xian 1980 GK Zone 20，在下面的选择坐标系统框中选择Beijing 1954 GK

11、 Zone 20N，在右边有一个按钮为Transformations.动态投影(ArcMap)所谓动态投影指，ArcMap中的Data 的空 间参考或是说坐标系统是默认为第一加载到当 前工作区的那个文件的坐标系统，后加入的数 据，如果和当前工作区坐标系统不相同，则 ArcMap会自动做投影变换，把后加入的数据投 影变换到当前坐标系统下显示！但此时数据文 件所存储的数据并没有改变，只是显示形态上 的变化！因此叫动态投影！表现这一点最明显 的例子就是，在Export Data时，会让你选择 是按this layers source data（数据源的坐标 系统导出），还是按照the Data （当

12、前数据框 架的坐标系统）导出数据！（注意在某些叠加 分析的时候会出错？）返回首先理解地理坐标系（Geographic coordinate system），Geographic coordinate system直译为 地理坐标系统， 是以经纬度为地图的存储单位的。很明显 ，Geographic coordinate syst em是球面 坐标系统。我们要将地球上的数字化信息 存放到球面坐标系统上，如何进行操作 呢 ？地球是一个不规则的椭球，如何将数据 信息以科学的方法存放到椭球上？这必然 要求 我们找到这样的一个椭球体。这样的 椭球体具有特点：可以量化计算的。具有 长半轴，短 半轴，偏心率。

13、以下几行便是 Krasovsky_1940椭球及其相应参数。 Spheroid: Krasovsky_1940 Semimajor Axis: 6378245.000000000000000000 Semiminor Axis: 6356863.018773047300000000 Inverse Flattening（扁率）: 298.300000000000010000 然而有了这个椭球体以后还不够，还需要一 个大地基准面将这个椭球定位。在坐标系统描 述中，可以看到有这么一行： Datum: D_Beijing_1954 表示，大地基准面是D_Beijing_1954。点击打开一个投影转换

14、对话框，可以在 对话框中看到Convert from和Into表 明了我们想从什么坐标系统转换到什 么坐标系统。在下方的using下拉框右边，点击 New.，新建一个投影转换公式，在 Method下拉框中可以选择一系列转换 方法，其中有一些是三参数的，有一 些是七参数的，然后在参数表中输入 各个转换参数。输入完毕以后，点击OK，回到之前的 投影转换对话框，再点击OK，就完成 了对当前地图的动态投影转换。这时还 没有对图层文件本身的投影进行转换， 要转换图层文件本身的投影，再使用数 据导出，导出时选择投影为当前地图的 投影即可。方法2：对于有大量图层需要进行投影转换时，这种手工操 作的办法显得比

15、较繁琐，每次都需要设置参数。可 以只定义一次投影转换公式，而在此后的转换中引用此投影转换公式即可。这种方法需要在ArcTools 中进行操作。在Data Management ToolsProjections and Tranformations下，有 Create Custom Geographic Transformation命令 。打开这个命令，选择输入和输出的投影，可 以是系统自带的也可以是自己设置的，选择 转换方法，与方法种介绍的类似，可选择 三参数或者七参数，然后输入各个参数指。 通过为这个投影转换公式指定一个名称，可 以在以后的操作中直接引用此公式而不用重 复输入各个参数了。点击

16、OK生成这个投影转 换公式。在方法一里面，我们是动态的改变了地图的投影 ，然后通过数据导出的办法将要转换投影的图层重新生成的。在这里，我们可以直接使用Data Management ToolsProjections and Tranformations下的Project命令，生成转换后 的图层文件，Project命令分别位于Feature和 Raster目录下，分别针对于矢量和栅格数据。在这 个命令中，在指定了输入的图层后，Input Coordinate System自动的识别出了输入的投影， 需要用户指定输出的投影，如果两者与之前定义投 影转换公式的输入和输入投影的话，在下面的Geographic Transformation下拉框中会出现之前 定义的公式名称，直接选择即可使用。这些参数怎么获取？ （1）直接去买 （2）如果有同一个点在不同椭球体下的 坐标可以在mapgis中实现，mapgis可 以反算出这几个参数。附加内容：高斯平面直角坐标