ArcGIS软件中的坐标系及数据类型的概念

1.坐标系基础知识

想简明的介绍坐标系的概念和它在ArcGIS软件中的应用，基础知识部分参考了《地理信息系统教程》（汤国安等著）。

1.1地球空间参考

对于规划类的同鞋们来说，学习使用GIS软件的过程中经常会碰到各种坐标系的问题，对于各种坐标系（北京54、西安80、国家2000等等）傻傻分不清楚。但是进行空间分析或则其他操作又要把空间数据纳入到相同的空间参考（坐标系）才行。今天我们就来探讨一下GIS和坐标系的关系还有它的应用。
首先提一个问题，我们用ArcMap软件大多时候是画地图，地图是和现实世界是一一对应的，但是我们身处的世界是立体的、空间的、三维的；但我们画的图确是扁平的、平面的、二维的。从三维世界到二维地图的这个过程是怎么实现的呢？这里就需要了解一些测绘的知识，有了这些基础呢方便理解后面应用时的操作。当然你如果对这些基础知识不感兴趣的话，也可以跳过，直接看操作部分就行了，但是还是建议看看基础知识。

1.1.1地球形状与地球椭球体

既然我们想研究地球做地图，那我们嘚让它方便计算，方便我们进行各种测量。
根据大地测量学很古很古的前辈大佬们的研究成果来看：因为地球表面是非常凹凸不平的（不规律），为了方便研究他们将地球模拟抽象成了4种球体。
第一个球：地球的自然表面
因为地球表面是非常不圆润、非常不规则（支楞八叉的）的表面，非常复杂，很难用来数学建模，在很多方面用地球的自然表面来考虑问题呢，很困难，不好使。但是不好使怎么办呢，那就让它规则点，圆溜点呗
                                图片来自bing搜索

第二个球：大地水准面
前面说了用地球自然表面很难搞，但是前人提出了一个相对抽象的面----大地水准面。这个面是一个理想状态下的假想的面，因为地球大部分被海水覆盖，前人假设在只考虑重力的情况下当海水不再流动，当海水不在流动，即海水的表面与重力方向处处正交，大家学过物理知道当只考虑重力，海水处处就只有向下的运动趋势，而不会流动。这样海面是圆润了，陆地呢，不处理了么（抛弃了么），显然不会，沿着刚才设想的静止的海水面延伸呗，这样这个连续、闭合的面就是大地水准面。注意地球上有很多重力等位面（水准面），只有与静止海平面重合的面才叫大地水准面。然后前人又引入了一个概念：把大地水准面包围的形体叫水准椭球也称大地体，可以理解把一个梨分为梨皮和梨剩下的部分。好了，现在第二个球介绍差不多了，我们把地球的自然表面提出了一点想象，让它圆润了一点，但是呢从结果看又没完全圆润，这是因为由于地球的内部质量分布不均匀，使得地球表面上各点的铅垂线方向产生不规则变化，因而大地水准面实际是一个复杂不规则的曲面，并不是严格的数学曲面，直接拿它来应用解决实际问题还是挺难搞的。但是大地水准面有明确的物理意义,在地球上实际存在,很大程度上能反应地球的真实情况,从而被沿用至今。

                                图片来自bing搜索
第三个球：地球椭球面
现在来看第三个球（终于敲到这个球了），前面说了我们对地球表面进行了一次“圆润”的过程，从地球表面到了大地水准面。但是用大地水准面来进行实际测量之类的工作还是很困难的，前辈们又提出了一个与大地体非常接近的旋转椭球，提出这个概念后前辈们觉得，哎这是个好用的地球模型，是个规则的数学曲面，方便计算和应用。所以在大地测量和GIS应用中一般都选择一个旋转椭球作为地球的模型，称为地球椭球公式就不复制上来了，有兴趣自己去查查哈。
不知道你注意到上面“一般都选择一个旋转椭球作为地球的模型”这句话没有，是不是说旋转椭球有很多个呢？地球椭球并不是一个任意的旋转椭球体，只有与水准椭球一致起来的旋转椭球才能用作地球椭球。但是这个“一致起来”怎么理解？按照我想的是二者比较吻合，即旋转椭球比较贴合大地水准面，这样用这个旋转椭球测量的数据才能贴近现实中的世界，精度也更高，当然这是我个人的理解不知道对不对（其实前面好多都是个人见解，不知道对不对，希望看到错误的小伙伴留言指出，以便及时指正。doge)。当然这个地球椭球也不是随随便便就确的这涉及非常复杂的大地测量学内容，大多数国家都有专门确定这个的机构，我就不展开了（主要是不会，狗头保命）。但(跳过的从这里开始)我们要知道的是不同的限制条件和不同的研究方法确定的地球椭球多种多样（这也就决定了为什么会有那么多不同的地理空间坐标系统），下面列出目前世界地图以及我国不同时期所采用的的地球椭球。

椭球名称	创立年代
WGS-84	1984
1975年国际椭球（西安80坐标系采用）	1975
海福特（hayford））（中国1953年以前采用）	1910
克拉索夫斯基（ красовский ）（北京54坐标系采用）	1940

有了地球椭球，在实际建立地理空间坐标系时，还需要指定一个大地基准面将这个椭球体与大地体（忘了啥叫大地水准面的你点一下这三个字）联系起来，在大地测量学中称之为椭球定位。按照我的理解这个过程的作用就是将这个椭球和你要研究的区域进行“配准”让坐标系更加精准。

1.1.2坐标系统

地理空间坐标系统通常分为球面坐标系统（GIS中的地理坐标系）和平面坐标系统（GIS中的投影坐标系），从二者的名称来看它两个的关系就是一个是三维的一个是二维的，平面（投影）坐标系是按照球面（地理）坐标与平面（投影）坐标之间一一对应的关系把三维坐标转绘到平面。

1.2空间数据投影

地球椭球体在数学上讲也是一个不能展开成平面的曲面，要把这样一个曲面表现成平面，就会发生变形。这种变形是不可避免的，变形又细分为长度变形、角度变形、和面积变形。对应的有等角投影（角度不会变形），等积投影（面积不变形可以用来计算面积）和等距投影（距离不变形）。

1.2.1常用的地图投影

高斯-克吕格投影（gauss kruger）
投影类型：横轴切圆柱等角投影
变形特点：角度不变形长度会变形
分带说明：分为6°带和3°带，6°带是从0°子午线（某条经线）起西向东每隔经差6°为一投影带，全球共分为60个分带。三°带与6°带类似只是每隔3°为一带。
横轴墨卡托投影（UTM）
投影类型：横轴割圆柱等角投影
变形特点：与高斯-克吕格相同
其他的投影因为书上没讲到我也不想查（大写的懒），而且这些基础内容感觉好像挺多的了(虽然还是省略了不少内容)，读起来比较累，影响阅读体验就敲到这打住。

2.ArcMap中的坐标系的应用

接下来就是今天讨论的重点了坐标系在ArcMap（后面简称GIS）中的应用，以上的基础内容皆是为此做准备。在GIS中我们用到坐标系的地方就是坐标系的投影与转换（注意：这是两个不同的概念）。投影 是针对地理坐标系将三维的球转绘成二维的平面图，投影坐标系是地理坐标系在投影这种“加工方法”下产生的产品类似父与子的关系；转换是将两个不同的地理坐标系之间进行转换。结合我在规划工作中的经历来介绍这两个操作的不同。
情景一：将有坐标的CAD文件导入到GIS中并叠加到影像上

GIS可以直接打开CAD的dwg文件，为了方便处理我们首先把cad数据转换成shapefile格式的数据，选择合适的位置保存但是因为dwg文件没有坐标系文件所以gis会蹦出这个提示框，提示你数据缺少坐标系信息，没有正确的坐标系GIS就不能把cad数据对应到它的坐标所在的位置，尽管cad数据是有坐标的，那我们应该怎么处理呢。最快捷的办法是问给你数据的那个人或则单位，问他知不知道坐标系是什么。如果能问到那是最好，直接用工具箱里的定义投影（定义投影只用来将未知或不正确的投影定义为正确的投影）工具就行了
但是我们可能会遇到给数据的那个人也不清楚坐标系的情况，这是后怎么办呢。不要忘了我们的CAD数据是有坐标的虽然不是有坐标系，我们可以从坐标里获取一定的信息。
可以看到右下角红色框框里的坐标整数部分是：6位，7位这种形式。这种形式代表了什么呢？代表了它是投影坐标系下的坐标，因为地理坐标系是由经纬线构成的，单位是度，投影坐标系是由地理坐标系按照某种方法投影出来的，投影坐标系有坐标原点，单位是米所以它的坐标的整数部分会很长，可以是6位，7位或则8位，7位这种两种组合。6位，7位是我们现在遇到的情况，那8位，7位是什么情况呢？其实是在前一种情况的X（东西向）坐标前加上带号像（36123456.001,1234567.001）这样。而且一般来说我遇到这种情况可以判断为西安80坐标系或则北京54坐标系或则国家2000坐标系（cgcs2000），首先判断西安80坐标系，因为这个坐标系在cgcs2000没普遍使用前用的最多的就是西安80坐标系，如果西安80不对的话试试国家2000对不对，最后在试北京54坐标系因为这个坐标系比较有年代感现在不是很常用。理论有了我们现在开始操作
在用定义投影这个工具时需要一个指定一个坐标系参数，我们这里先用西安80坐标系来尝试一下（我这个CAD数据在导入前已经问到了它是西安80坐标系但具体是哪个分带不清楚）。

我们点开投影坐标系的文件夹会发现有这么多先坐标系，选哪个呢？这时候我们需要借助场外工具的帮助了，对于规划行业来说这个cad数据你一定是知道它的实际位置的（我们现实生活中它所对应的位置），比如我这里的数据是贵州省黄平县谷陇镇岩门司村的cad数据，现在我们在浏览器里搜索天地图贵州进入它的在线地图服务

然后搜索框里搜索黄平县谷陇镇我们可以看到它的经纬度坐标是：经度:108.05 纬度:26.83

因为中华人民共和国成立后，我国大中比例尺地形图一律规定采用以克拉索夫斯基椭球体元素计算的高斯-克吕格投影。而上图是高斯投影的分带计算公式。根据获取的坐标我们来计算应该选择的分带带号，直接计算3度分带就行了，因为1:10000以下的小比例尺地图采用6度分带，像我们一个村子大小的地方应该属于大比例尺用3度分带来算（不知道这样想有没有问题，但是后面用这样的思路来搞还是正确的）。n=108/3=36，所以选择36分带。计算了分带之后再怎么处理呢？如果我们的cad数据坐标是带分带号的那种（8位，7位的形式）我们直接选哪个有代号数字的投影名称就行了，比如下面我们选中的这个

在投影的名字后面的数字就是带号36，但是我们的数据貌似是不带带号的那种这怎么办呢？再介绍一下3度分带投影坐标系的命名方式（以下面两个投影坐标系为例子讲解：）你就知道怎么办了。

第一种Xian_1980_3_Degree_GK_Zone_36：它代表的含义是西安80地理坐标系进行3度分带的高斯投影下的第36度带横坐标前加带号。Xian_1980代表西安80坐标系（简称西安80）；3_Degree代表3度分带；
GK是Gauss_Kruger（高斯克吕格）的缩写；Zone_36代表地36度带区域；
第二种Xian_1980_3_Degree_GK_CM_108E 解释一下CM_108E的意思其他的和第一种没区别CM_108E就是Central_Meridian（中央子午线）是东经108度横坐标前不加带号；
细心的同学可能已经发现其实下面的框框里是有具体的说明的，而且这两种投影其实是一回事只是第一种投影会在X坐标前加上带号，第二种不加带号，就是这个区别，它俩的中央子午线都是108。然后我用啦这么多篇幅去解释，不知道是不是有点多余。
西安80投影文件夹里还有这么一种比较短的投影它代表6度分带

现在在了解这个文件夹里的各种投影后我们就知道选什么了，36度带坐标不加带号那就选Xian_1980_3_Degree_GK_CM_108E这个，因为我们知道要选第36度带而且是坐标不带带号的那种，前面那两个一个是带带号的，一个不带那就选不带的呗。
然后把我们定义投影过的数据和遥感影像同时添加进来看看效果
可以看到是叠起来的，达到了预期想要达到的效果。
虽然影像和刚转换的数据坐标系不一样但还叠起来了，这是因为ArcGIS本身有一个动态投影机制，就是只要我们是同一个地方的地理要素，只要每个数据都定义了正确的坐标系，那么无论各要素是采用了什么坐标系，都是可以直接叠加对齐的，这也就是上文讲到坐标变换并不是必须的原因，就是说我们在做分析或者可视化的时候，并不一定都要变换为同一个坐标系，只要定义了正确的坐标系注意此处我多次强调“正确”，这个很重要。），一般都是能自动对齐的。

先敲到这剩下的后面再敲。所以，席，你看懂了么？ 不懂我们再讨论哈

=============================================================================

3.ArcGIS中的数据介绍

敲不动了，后面有空在写了