做ArcGIS实验做到坐标系这一章很多人会突然觉得脑子不够用。数据加载进去显示不出来两个图层明明看着是同一片区域却隔着十万八千里算面积算出来跟实际差了几个数量级这些问题几乎都绕不开坐标系。这次实验六的名字就叫ArcGIS坐标系总结听起来像是个收尾总结实际上它是把前面所有实验里隐含的那条线单独拎出来摊开讲。我的感受是坐标系这关过不了后面的缓冲区分析、叠加分析、栅格计算全都站不住脚因为你所有的操作对象如果不在同一个空间参照下算出来的东西就是无源之水。这篇内容我会按自己在实验里实际操作的顺序来讲先解释坐标系到底在解决什么问题再把常见坐标系摆在一起比较接着把ArcGIS里查看、定义、转换的核心操作一步步走通最后集中说几个最容易翻车的坑和排查思路。不管你是第一次做这个实验的同学还是工作中被坐标系搞到怀疑人生的从业者跟着这篇把逻辑理一遍很多问题当场就能看明白。1. 坐标系到底解决什么问题先想明白再动手1.1 地理坐标系用经纬度描述位置的规则地理坐标系Geographic Coordinate SystemGCS用经纬度来表达一个点在地球表面的位置。经纬度本身是角度单位经度从本初子午线往东往西各0到180度纬度从赤道往南北各0到90度。这套系统看起来简单但它有一个隐含问题地球并不是一个规则的球体而是一个赤道鼓起、两极稍扁的椭球体你要描述某个点在哪里必须先规定你把这颗星球当成什么形状。这就引出了椭球体和基准面的概念。椭球体Ellipsoid/Spheroid是一个用长半轴和扁率定义的数学模型用来近似地球形状基准面Datum则在椭球体的基础上规定了椭球中心、旋转轴方向以及和大地水准面之间的拟合关系。也就是说同样的经纬度坐标放在不同的基准面下对应的物理位置是不同的。举个日常能感受的例子同一个经纬度值在WGS84基准面下和在西安80基准面下实际位置可能相差几十米甚至上百米。这个差异在电子地图上看不出来但一到了需要精确边界、地籍测量、工程放线的场景就是天壤之别。做实验的时候很多同学看到图层属性里写着一串类似GCS_WGS_1984的字符就以为坐标系都一样这恰恰是后面所有混乱的起点。地理坐标系只是第一层它告诉你用什么经纬度规则但还没有解决怎么把地球表面画到平面上的问题。1.2 投影坐标系为了让平面计算有意义投影坐标系Projected Coordinate SystemPCS是在地理坐标系的基础上通过某种数学变换把椭球面上的经纬度换算成平面上的X、Y坐标。为什么要多此一举因为计算机绘图、测量面积、计算距离、做拓扑分析都必须在平面直角坐标系下才能快速又合理地进行。你在球面上量两点之间的距离要用大圆公式但把这一片区域投影到平面之后直接用勾股定理就能算了效率完全不同。常见投影方式包括高斯-克吕格投影Gauss-Kruger、横轴墨卡托投影UTM、墨卡托投影、兰伯特等角圆锥投影、阿尔伯斯等积圆锥投影等等。每种投影都有自己的取舍有的保证角度不变等角有的保证面积不变等积有的在某个方向上距离变形最小。没有任何一种投影能同时保持角度、面积、距离全部不变这是数学上已经被证明的结论。所以选择哪种投影坐标系本质上取决于你要拿数据做什么。投影坐标系在ArcGIS里通常写作某某投影坐标系_基于某地理坐标系的形式比如WGS_1984_UTM_Zone_50N含义就是底下的地理基准是WGS84投影方式是UTM第50带这个坐标系属于北半球。看到这种名字你能立刻读出三件事基准面是什么、投影方式是什么、适用于哪个区域。1.3 椭球体、基准面和坐标系三层概念别混为一谈我在实验里给学生讲的时候喜欢把这三层比作做衣服的过程。椭球体是先画一个人的标准体型底模基准面是把底模往真人身上套、调整到某个地区最合身的那一次裁剪坐标系则是最终用尺子量出来的一套标注规则。GCS里面包含了基准面和角度单位PCS里面包含了GCS、投影方法、线性单位和一堆投影参数。数据属性里那一长串名称其实已经把这三层信息都压缩进去了。这三个概念一旦混淆就会出现下面这些常见错误有人把北京54和西安80当成两种投影方式其实它们是两套不同基准面下的地理坐标系统有人觉得WGS84只是经纬度坐标系忽略了它也有对应的UTM投影坐标系还有人一看到投影两个字就以为必须把数据转成平面结果把本来已经是投影坐标系的图层又强行投影了一遍数值全乱。搞清楚这三层的关系之后ArcGIS里的很多操作逻辑就变得非常顺理成章定义投影是在告诉软件这套数据用的是哪套底模和尺子投影转换是把一套标注规则换算成另一套标注规则的过程。接下来我们看看具体有哪些坐标系可选。2. 常用坐标系选型不同场景下的正确选择2.1 全球层面的 WGS84 与 UTMWGS84全称World Geodetic System 1984是全球定位系统GPS直接采用的地理坐标系也是目前国际上应用最广泛、被互联网地图和导航设备默认采用的标准。它的椭球参数是长半轴6378137米扁率倒数298.257223563。大多数手持GPS设备导出的坐标默认都是WGS84经纬度。UTMUniversal Transverse Mercator投影则是建立在WGS84等基准之上的常用投影体系它把全球从西经180度开始按每6度一带来划分共60带每带中央经线采用比例因子0.9996来减小边缘变形北半球还有500000米的假东偏移。在ArcGIS里新建要素类的时候只要选择WGS84地理坐标系或对应的UTM投影带全球通用的数据基本都能对接上。这里有一个实操细节值得记下来如果数据范围是跨带的大尺度区域比如跨越好几个UTM带就不建议硬套某个单带投影而要考虑使用适合大范围的投影方式或者干脆保留经纬度进行分析。UTM每一带在中央经线附近变形很小但到了带的边缘距离变形可能达到千分之一左右对精度有要求的分析就得小心了。2.2 国内数据常见的北京54、西安80和CGCS2000国内数据接触最多的坐标系有三个关键词北京54、西安80、CGCS2000。北京54坐标系采用克拉索夫斯基椭球长半轴6378245米扁率倒数298.3始建于上世纪五十年代是一套参心坐标系椭球中心与地心不重合。西安80坐标系采用IAG75椭球长半轴6378140米扁率倒数298.257同样属于参心坐标系但是椭球定位更符合我国大陆的实际情况。CGCS20002000国家大地坐标系则是地心坐标系椭球长半轴6378137米扁率倒数298.257222101从全国范围看它和WGS84的差异通常很小大部分情况下可以忽略但在高精度应用里仍要按实际需求处理。从实验和实际工作的角度我最常给出的建议是手头的数据若涉及历史图件先确认它是哪一套基准面再谈后续分析新采集的数据尽可能统一到CGCS2000因为它和GPS/WGS84的兼容性最好如果在混用北京54和西安80的数据必须明确做基准面转换不能直接把坐标数值当成同一套体系来用。北京54和WGS84之间的平面位置差异在中低纬度地区普遍有几十米到上百米不做转换直接叠加边界错位一眼就能看出来。2.3 互联网底图数据里的坐标系偏移现象第三个需要了解的现实问题是互联网地图和卫星影像的坐标系并不总是直接对标WGS84。很多商用在线底图由于自身服务策略的原因会采用经过偏移处理的标准这导致你在ArcGIS里看到的影像内容和它标注的坐标系名称之间并不严格对应。具体表现为把野外采集的GPS点加载到ArcGIS里叠加在线影像发现点全部偏到马路对面甚至几百米外。很多同学第一反应是投影搞错了实际上可能是底图数据的坐标偏移问题。做实验的时候我建议尽量使用不依赖在线图源的数据或者使用带有明确坐标系元数据、并且和原始数据坐标系一致的本地数据这样能少踩很多莫名其妙的坑。如果实验内容确实需要影像底图可以把在线底图只当作参考背景而不作为精确分析的对象一切定量计算都用本地矢量数据完成。3. 实验六的核心操作查看坐标系、定义投影、坐标转换3.1 如何快速判断一个图层当前处于什么坐标系在ArcGIS里查看一个数据的坐标系方法很简单右键图层打开属性Properties切到源Source选项卡就能看到完整的坐标系描述。如果是矢量要素类坐标系信息还会写明线性单位是什么、投影参数有哪些如果是栅格数据除了坐标系还会显示像元大小和范围。我比较常用的是把要素类和栅格都拖进地图然后到图层属性→源里核对三样东西第一底下的地理坐标系名称第二如果是投影坐标系投影方法和中央经线、比例因子等参数第三单位是不是米。光是单位一条就能筛掉很多问题。比如数据源的线性单位如果显示为Degree而你在算面积时却期望得到平方米结果肯定会差很远。ArcCatalog里同样可以看到这些信息而且ArcCatalog更适合批量查看多个数据集的坐标系状况一边是列表模式一边在预览里看坐标范围。实验里如果给你一批数据让你判断坐标系最快的办法就是用ArcCatalog或属性面板挨个看把坐标系名称、基准面、投影类型记录下来整理成一个小表格。这个表格会成为你后面所有操作的出发点。3.2 定义投影和投影的区别一个改标签一个做变换这一步是整个实验里最容易混淆的地方。ArcToolbox里有两个名字极像的工具一个叫定义投影Define Projection在数据管理工具→投影和变换下另一个叫投影Project也在同一个工具集下。它们的区别是本质性的定义投影只是给数据贴一个正确的坐标系标签并不改变数据里任何坐标数值投影则是把坐标值从一套坐标系真正换算到另一套坐标系坐标数值会发生变化。什么时候用定义投影最常见的情形是你拿到一份数据形状位置都对但它缺少坐标系信息属性里显示未知或No Spatial Reference同时你确定它实际对应哪个坐标系这时用定义投影把它标回去。另一种情形是数据本身有坐标系信息但写错了需要在确认无误的前提下把它改回来。什么时候用投影Project当你有两套不同坐标系的数据要叠加分析并且它们是同一物理区域但基准面不同或者基准面相同但投影方式不同时就需要用投影工具把其中一个转换到另一个的坐标系下。记住一个判断口诀坐标值没变只是标签错了→定义投影坐标值需要换算→投影。把这两个工具搞反是所有坐标系实验事故中最常见的一种没有之一。3.3 矢量数据投影转换的完整操作流程矢量数据的投影转换标准操作是用工具箱里的Project工具。具体流程是打开ArcToolbox→数据管理工具→投影和变换→要素→投影然后在对话框里选择输入要素再点输出坐标系旁边的图标选择一个目标坐标系或者浏览现有数据来直接导入它的坐标系。如果两个坐标系属于不同基准面对话框里会要求选择地理变换方法和参数如果是相同基准面下的不同投影则不需要地理变换直接可以完成换算。这里要说一下地理变换参数。从WGS84转到CGCS2000很多情况下可以不选变换或者使用简化参数但从北京54转到WGS84或者西安80转到CGCS2000就必须提供合适的变换方法比如三参数、七参数、区域拟合等多种方法。参数从哪里来对于实验教学通常由老师直接给出一套参数对于实际生产项目一般由测绘部门提供覆盖该区域的转换参数。拿不到准确参数却强行转换结果比不转换还糟这一点务必记清楚。做完转换之后数据会生成一个新的要素类。我建议用源选项卡核对一遍新数据的坐标系名称再用测量工具在图上量一段已知距离来验证。如果转换前和转换后测出来的距离差不多说明整个流程是对的如果量出来差了一倍或者单位不对多半是目标坐标系选成了经纬度或者线性单位没设置对。3.4 栅格数据的投影转换与重采样栅格数据和矢量不太一样因为投影转换会改变像元的位置和形状必须配合重采样方法才能生成新栅格。ArcGIS里的工具叫做栅格投影Project Raster在数据管理工具→投影和变换→栅格下。界面里除了输入栅格和目标坐标系之外还会要求选择重采样技术常见的有最邻近Nearest Neighbor、双线性Bilinear、三次卷积Cubic Convolution三种。我的经验是分类数据如土地利用类型、植被分类用最邻近因为它不会在像元之间插值出新的类别值能保持类别完整连续型数据如高程、气温、降雨量用双线性或三次卷积得到的表面更平滑。还有一个细节栅格重投影之后输出像元大小可能发生变化需要手动设定目标像元大小否则ArcGIS会按默认方式算出一个新值不同软件之间的默认逻辑还不太一样。投影转换后的栅格还会遇到一个现象数据范围会比原来大一圈因为栅格旋转后边缘会出现空白像元这些像元会被赋予NoData值。在后面的分析中如果这些空白区域影响到了统计结果可以考虑先裁剪再分析或者把范围设置得更紧。这些细节在实验报告里写清楚是加分项。3.5 动态投影的坑看起来对了不等于真的对了ArcGIS默认开启动态投影On-the-fly Projection也就是在地图上显示数据时软件会临时把不同坐标系的数据统一重投影到当前数据框的坐标系下让你能够叠在一起看。这个功能非常方便但也带来一个很大的认知陷阱动态投影只是显示层面的转换它并不改变数据本身的坐标系和坐标值。做过实验的同学应该见过这种场景加载不同来源的两份数据数据框坐标系是WGS84软件自动把两份数据都动态投影到WGS84经纬度下显示看起来位置也大致对。于是很多同学就以为数据已经转换好了实际上原始数据还是它们各自原来的坐标系如果直接对它们做空间分析、算面积、量距离软件会在后台临时重投影但结果未必符合你的预期。我一直强调的做法是在正式分析之前先把所有参与分析的图层统一到同一个投影坐标系下并另存为新的数据确保底层数据真的是同一套坐标系。怎么判断底层数据是否已经统一右键每个图层看源选项卡里写的是什么坐标系只要两个图层的坐标系描述不完全一致就别急着开分析工具。动态投影适合快速看图不适合作为分析依据这条要刻在脑子里。4. 实验中最常翻车的场景与完整排查链路4.1 两个图层明明都是WGS84却对不上位置这是我在实验里遇到频率最高的问题。一份要素类显示坐标系是WGS84另一份影像坐标系也是WGS84但叠加之后点跑到影像边缘完全对不上。这时候很多人第一反应是软件坏了其实大多数情况跟软件无关。你仔细看两个数据的坐标系发现一个是GCS_WGS_1984另一个是WGS_1984_Web_Mercator虽然都带WGS84字样但前者是地理坐标系经纬度后者是投影坐标系以米为单位的Web墨卡托。由于Web墨卡托把全球纬度限制在约±85.06度且经过特殊的球面墨卡托模型处理坐标数值和真实经纬度有系统性的差异所以直接叠加必然错位。只要把其中一个数据转换到和另一个一致的坐标系或者把数据框坐标系设置成一个通用的投影坐标系并让两者都动态投影问题立刻就能解决。另一个隐藏原因是基准面虽然同名但描述文本不同。比如有的数据写GCS_WGS_1984有的写成D_WGS_1984这通常是同一回事但如果一个是GCS_Beijing_1954却被人为改名为GCS_WGS_1984用定义投影随便改的那么就算名字一样实地的位置也不会对齐。这种情况下你得回到数据源确认真实基准而不是相信标签。4.2 面积、距离计算结果离谱实验里经常要做面积统计。某个区域用某种坐标系算出来面积是几千平方千米换个坐标系算出来变成几十万而且单位变成平方度之类的怪东西这就是目标坐标系选成了经纬度或者干脆在数据框是地理坐标系的情况下直接用几何计算面积。ArcGIS里如果数据框是GCS经纬度几何面积工具返回的单位就不是平方米而是平方度数值自然千奇百怪。解决思路很清晰分析之前把数据转换到合适的投影坐标系。投影选择上如果分析区域是全国范围的用兰伯特等积圆锥或阿尔伯斯等积投影比较合适因为它们保证面积不变如果分析区域是某个省或地市用高斯-克吕格3度带或6度带都可以但要注意选择离区域最近的中央经线。高斯-克吕格投影不是等积投影它的面积变形在离中央经线远的地方会变大跨带分析时误差可能比较明显。我还见过一个有趣的案例A同学用高斯-克吕格投影算面积结果发现比真实值大了将近1%查了很久才发现是选错了投影带——数据本身落在21带他却选了23带的坐标系定义中央经线差了6度变形值完全不同。所以做面积之前先查一下数据落在哪个带、对应哪条中央经线非常有必要。4.3 数据完全没有坐标系信息第三种高频翻车场景拿到的矢量数据在ArcGIS里显示未知坐标系Spatial Reference: Unknown加载时软件还会给出一个警告提示缺少空间参考信息。很多同学顺手点了确定然后开始在未知坐标系的地图上做分析最后得到的结果完全没法解释。正确做法是先停下来搞清楚这份数据本来的坐标系是什么。判断依据一般有三个来源第一数据来源方提供的元数据说明第二同类数据的默认坐标系比如某单位以前的图纸都是西安80第三数据自身的坐标范围和单位。比如一个数据的范围X大约是500000到540000Y大约是3300000到3400000而且看起来单位是米大致能推断它可能是某投影带下的高斯-克吕格坐标如果X范围是100到120、Y范围是20到40那八成是经纬度。推断出坐标系之后用定义投影工具把它标回去再加载到地图里和已知地物核对如果吻合就说明判断正确。要注意的是定义投影一旦执行数据的XYZ坐标值本身不会变它只是给这份数据补上了坐标系标签。如果标签补错了后面所有的转换都会基于一个错误前提所以定义投影前最好备份原始数据并且把推断依据写进实验报告里。5. 一次完整的实验六实操记录5.1 实验数据准备为了把前面的理论串起来我在这里模拟一份实验六的典型数据组合一份某市中海拔区域的边界矢量数据用北京54/高斯-克吕格6度带表示单位米一份同区域的遥感影像用WGS84/UTM投影表示一批野外采集的GPS点WGS84经纬度。实验目标是把三份数据统一到同一个坐标系下完成面积统计和距离测量。我把数据命名为模拟项目X的boundary、image、gps_points三个图层。注意如果你在真实环境里遇到类似数据先不要急着动手第一步永远是把三份数据的原始坐标系信息截图存下来这是实验报告最重要的依据之一。5.2 分步实验操作记录第一步用ArcCatalog或属性面板查看三份数据的坐标系描述分别记录如下boundary是Beijing_1954_GK_Zone_21N北京54基准面、高斯-克吕格21带、北半球image是WGS_1984_UTM_Zone_50NWGS84基准面、UTM50带gps_points是GCS_WGS_1984WGS84地理坐标系、经纬度。第二步因为gps_points是经纬度其余两份数据是投影坐标如果直接把三份数据加进来叠加靠动态投影虽然能显示但统一性不足。所以我把boundary作为基准数据先把它设置成数据框的坐标系然后把image和gps_points分别转换到boundary的坐标系上。这一步我用Project工具完成对gps_points定义目标坐标系Beijing_1954_GK_Zone_21N由于基准面从WGS84变成北京54需要指定地理变换方法这里按教学数据给出的区域三参数进行变换。对image我在输入栅格和目标坐标系之间选择相同的三参数变换同时把重采样设为双线性像元大小保持与原始影像一致。第三步投影转换完成后把三份转换后的数据统一加载进地图。在图层属性→源里确认三者坐标系完全一致然后叠加检查GPS点应当落在边界内影像底图应该和boundary边界重合边缘没有大的错位。如果出现超过一两个像元的系统性偏移就要回头检查地理变换参数是否选择正确。第四步在统一坐标系下做两个验证实验一是用几何计算工具统计boundary的面积二是用测量工具量取两个GPS点间的距离。将结果与按原始数据分别计算的结果做对比验证投影转换是否影响了数值。如果面积结果出现异常检查数据框坐标系是否是投影坐标系而不是经纬度。这四步做完实验六的核心内容就基本完整了。整套操作的核心其实只有一句话先确认数据各自是什么坐标系再统一目标坐标系最后在统一坐标系下做分析和验证。5.3 结果对比与常见偏差速查表为了写报告方便我把实验里常见的对比结果整理成一个表格方便你对照自己的输出结果判断问题出在哪一步。下表列出的是在相同区域、不同坐标系/不同操作下面积和距离计算结果的典型差异方向场景计算工具条件结果表现原因数据框为WGS84经纬度直接算面积得到平方度数量级异常单位不是米不能在球面经纬度下直接算平面面积数据框为北京54高斯投影投影前计算面积正确单位米且投影带匹配区域数据框投影带选错选成相邻带面积偏大或偏小约0.1%-1%高斯-克吕格离中央经线越远变形越大北京54数据直接命名WGS84不转换硬叠加位置偏移几十米以上基准面不同需地理变换参数栅格重投影用最邻近对分类数据类别保持边缘有锯齿重采样不插值栅格重投影用三次卷积对高程数据表面平滑可能改变极值插值产生中间值这张表我建议原样贴进实验报告再配合每一行的截图说明基本能一眼看出你确实理解了坐标系的核心逻辑。如果你在实际操作中得到了表中没列出的异常结果最值得优先检查的两件事是数据框的坐标系是不是投影坐标系以及参与分析的数据是否真的已经被转换到了同一坐标系。6. 做完实验想补充的三点经验第一个经验是坐标系实验不像其他实验那样有固定唯一的标准答案它更考验你能不能判断自己当前看到的结果是否合理。每次操作完我都习惯给自己提三个问题单位对不对范围是不是变了和已知的参照物对不对得上三个问题都通过才敢说转换成功。这三个问题看似简单但绝大多数坐标系事故最后都能回溯到其中某一个没答上来。第二个经验是养成随手记录坐标系信息的习惯。拿到任何一份数据先把坐标系描述复制到实验文档里最好连投影参数也存下来。这个习惯在实验阶段看起来有点多余等到了做多源数据融合或者处理测绘成果时你会发现它帮你省下大把排查时间。很多时候问题不是转换算法不懂而是从一开始就没人知道某份数据的坐标系到底是谁定义的元数据缺失是比任何工具错误都要棘手的难题。第三个经验是多利用ArcGIS自带的工具文档来学习而不是死记硬背工具名。把鼠标悬停在一个工具上看它对参数的解释尤其是输入和输出的单位说明很多困惑其实文档里早就写明白了。坐标系这个主题初学觉得绕主要是因为概念有嵌套关系但只要把椭球体、基准面、地理坐标系、投影坐标系这四层在脑子里理顺配合一次完整的实验操作之后再遇到更复杂的投影问题都不过是换个区域、换套参数重复同样流程而已。