GIS中关于投影的几点心得

写在开始

以下的内容或是看资料、视频，或是借鉴前人自学而来，仅作参考，如有错误，敬请谅解，欢迎批评指正！

投影变换的几个概念

之前一直搞不懂正解变换和反解变换，恰好看见了相关的介绍，就简单记录一下好了。

正解变换：通过建立一种投影变换为另一种投影的严密或近似的解析关系式，直接由一种投影的数字化坐标x、y变换到另一种投影的直角坐标X、Y。

反解变换：由一种投影的坐标反解出地理坐标（x、y→B、L），然后再将地理坐标代入另一种投影的坐标公式（B、L→X、Y），从而实现由一种投影的坐标到另一种投影坐标的变换（x、y→X、Y）。

数值变换：根据两种投影在变换区内的若干同名数字化点，采用插值法、有限差分法、最小二乘法、有限元法、待定系数法等，从而实现由一种投影的坐标到另一种投影坐标的变换。

GIS中的几个投影工具介绍

GIS中ArcToolBox->Data Management Tools ->Projections and Transformations中提供了几个工具：

1. Define Projection
2. Feature -> Project
3. Raster -> Project Raster
4. Create Custom Geographic Transformation

若数据没有任何空间参考信息，那么需要先用工具1给数据定义一个投影之后，再用工具2和工具3进行投影变化，其中工具2对矢量数据进行操作，工具3对栅格数据进行操作。

Create Custom Geographic Transformation
创建自定义地理坐标变换
可创建一种变换方法，用于在两个地理坐标系或基准面直接对数据进行转换。对于任何参数要求进行地理变换的工具，都可使用此工具的输出作为变换方法。

几个投影变换情景

在进行坐标转换时，工具上有一个选项，Geographic Transformation
(可选)，有时候发现这个参数必须填，有时候又可以省路，这是为什么呢？

1.不涉及地理坐标系变换的坐标变换，不需要参数。

如：从GCS_Xian_1980变换到Xian_1980_3_Degree_GK_CM_120E，这个过程不涉及地理坐标变换，因此无需参数。

2.涉及到地理坐标系的变化，并且GIS中有转换方法，这个参数是必须的，可以直接从列表中选择。

如：从GCS_Beijing_1954转换为GCS_WGS_1984坐标系。

3.涉及到地理坐标系的变化，这个参数是必须的，但是列表中没有，

需要自己定义一个。 如：从GCS_Beijing_1954转换为GCS_Xian_1980坐标系。

自定义采用前面的Create Custom Geographic Transformation
创建自定义地理坐标变换工具进行。

例子

如何将两个平面投影互相转换？

以下载的图层为例，将墨卡托平面投影坐标系（WGS 1984 UTM ZONE 50N）转换为北京54坐标系

首先定义一个投影，如图选择工具箱，Data Management Tools/Projections and Transformations/Create Custom Geografic Transformation。按照下图设置之后，点击ok即可。

定义完成之后，选择工具箱中投影工具条，Data Management Tools/Projections and Transformations/Raster/Project Raster，打开工具条。

根据图示方法进行转换，最后点击ok完成。

加载投影后的图层，查看属性，可以看到坐标系已经换了。

如何将WGS84坐标系投影到北京坐标系？

同理先定义一个转换，然后再进行投影。
1.定义转换

2.进行投影

将北京坐标系投影到WGS1984
此时不用定义转换，可以直接进行投影，利用project raster，坐标系投影方法的选择如下：