gdal坐标转换总结

首先，在进行坐标转换之前，有必要先了解一下有关坐标系的几个基本概念。

地理坐标系（Geographic Coordinate Systems）

地理坐标系是一个球面的坐标系统，以经纬度为单位，它由椭球体和大地基准面两个部分组成。

椭球体（spheroid）

我们要将地理信息以球面坐标系的方式表达，首先需要找到一个可以量化计算的椭球体。一个椭球体的确定需要以下参数：长半轴、短半轴、偏心率，其中偏心率可根据长短半轴计算得到。

例如，WGS84椭球的参数如下：

Spheroid(椭球名):"WGS_84";
Semimajor Axis(长半轴):6378137
Semimajor Axis(长半轴):6356752.3142
Inverse Flattening(扁率):1/298.2572236

大地基准面(datum)

有了椭球体以后，还需要一个大地基准面将这个椭球定位。

大地基准面（Geodetic datum），设计为最密合部份或全部大地水准面的数学模式。它由椭球体本身及椭球体和地表上一点(原点)之间的关系来定义。此关系能以 6个量来定义，通常是大地纬度、大地经度、原点高度、原点垂线偏差之两分量及原点至某点的大地方位角。

同一个椭球面，不同的地区由于关心的位置不同，需要最大限度的贴合自己的那一部分，因而大地基准面就会不同。

有了Spheroid和Datum两个基本条件，便可以确定一个地理坐标系统。

投影坐标系

将球面坐标转化为平面坐标的过程称为投影。因此，投影坐标系实质上是在地理坐标系的基础上通过投影得到的。投影坐标系其单位通常为m。

例如我国常用的高斯-克吕格投影,其通常是按6度和3度分带投影，1:2.5万－1:50万比例尺地形图采用经差6度分带，1:1万比例尺的地形图采用经差3度分带。具体分带法是：6度分带从本初子午线开始，按经差6度为一个投影带自西向东划分，全球共分60个投影带，带号分别为1－60；3度投影带是从东经1度30秒经线开始，按经差3度为一个投影带自西向东划分，全球共分120个投影带。为了便于地形图的测量作业，在高斯-克吕格投影带内布置了平面直角坐标系统，具体方法是，规定中央经线为X轴，赤道为Y轴，中央经线与赤道交点为坐标原点，x值在北半球为正，南半球为负，y值在中央经线以东为正，中央经线以西为负。由于我国疆域均在北半球，x值均为正值，为了避免y值出现负值，规定各投影带的坐标纵轴均西移500km，中央经线上原横坐标值由0变为500km。为了方便带间点位的区分，可以在每个点位横坐标y值的百千米位数前加上所在带号，如20带内A点的坐标可以表示为YA=20 745 921.8m。

ArcGIS中常用的投影坐标系命名方式

• Beijing 1954 3 Degree GK CM 117E

  北京54坐标系 3度分带投影 高斯克吕格 中央经线东经117度 横坐标前加带号

• Beijing 1954 3 Degree GK Zone 25

  北京54坐标系 3度分带投影 高斯克吕格 带号25 横坐标前加带号

• Beijing 1954 GK Zone 13

  北京54坐标系 6度分带投影 高斯克吕格 带号13 横坐标前加带号

• Beijing 1954 GK Zone 13N

  北京54坐标系 6度分带投影 高斯克吕格 带号13 横坐标前不加带号

ArcGIS动态投影

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

坐标系的通用描述方法

WKT

WKT，全程 well know text，是一种文本标记语言，用于表示矢量几何对象、空间参照系统及空间参照系统之间的转换，该格式由开放地理空间联盟(OGC)制定。

例如，WGS84地理坐标系的OGC WTK定义如下：

GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4326"]]

ESRI使用的WTK描述与OGC标准的WTK存在区别，因此，有时候会需要用到ESRI WTK，例如，WGS84地理坐标系的ESRI WTK如下：

GEOGCS["GCS_WGS_1984",
    DATUM["D_WGS_1984",
        SPHEROID["WGS_1984",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["Degree",0.017453292519943295]]

在arcgis中，.prj文件中存储的即是要素的坐标系的WKT文本表示，可用记事本打开查阅。

EPSG

EPSG的英文全称是European Petroleum Survey Group，中文名称为欧洲石油调查组