地形分析实战指南：从DEM选型到坡度水文分析

做地形分析这些年我最大的感受是很多人卡住的不是软件操作而是不知道手里的结果到底在算什么、算得对不对。上周还有个做规划的朋友拿着一个坡度图问我为什么同一片地方在ArcGIS里算出来是15度换到QGIS里就变成了17度。这个问题听着简单但真要讲清楚得把DEM、算法、投影、分辨率这些底层的东西全翻出来。所以这篇我干脆把地形分析这条线上的核心问题从头捋一遍从DEM的选型到坡度坡向、山体阴影、视域分析、水文分析再到等高线和剖面最后说说实战里经常踩的坑。内容偏实操原理也会讲适合刚入门GIS的学生、做资源环境和规划项目的同行以及打算做GIS二次开发的工程师参考。1. 一口吃不下地形分析先从DEM数据这件事说起地形分析看起来是一堆工具的组合但真正决定分析质量上限的其实是你手里那份底图数据。很多人一上来就打开Spatial Analyst工具开始算坡度、算坡向结果算出来的图自己都不敢信问题往往出在DEM本身就不可靠。1.1 DEM的三种形态选错了一步步都白搭DEM数字高程模型在GIS里有三种常见表达方式栅格GRID、不规则三角网TIN和等高线。做地形分析绝大多数工具都要求输入栅格DEM因为坡度、坡向、水文分析这些算法本质都是在规则格网上做邻域运算格网结构最简单、计算效率最高。TIN的优势在于它能在地形变化剧烈的地方加密三角面比如陡崖、山谷用更少的存储表达更精细的地形细节但它的数据结构复杂直接拿来做邻域分析很费劲通常是把TIN转成栅格后再用。等高线是大家最熟悉的但等高线本身不是连续表面直接拿等高线做坡度分析需要先插值成栅格插值方法选不好锯齿、平顶、伪洼地全跑出来了。我平时接项目如果甲方只给等高线我会先仔细看等高距和地形特征点山顶、鞍部、洼底有没有标全。等高距大、地形又破碎的话插值出来的DEM基本只能做宏观展示不能做精细的水文分析。1.2 分辨率、坐标系和精度这三个参数决定了你的分析可信度先说分辨率。DEM栅格的分辨率通俗讲就是每个格子代表地面上多大的范围。30米分辨率的DEM一个格子对应30乘30米的面积你要是拿它去做一个几百米宽的光伏场址的坡度分级那勉强够用但如果要用它去规划一段5米宽的排水沟结果就没有实际意义了。分辨率和地形特征尺度不匹配是地形分析结果失真的第一大原因。再说坐标系。很多人在做地形分析时用的是WGS84经纬度坐标单位是度而坡度、坡向、面积、长度这些计算依赖平面单位米。如果在经纬度坐标下直接算坡度算法会把经度和纬度方向的格网距离都当成“度”来算算出来的数值完全没有物理意义。所以做地形分析之前第一步一定是把DEM投影到合适的平面坐标系比如高斯-克里格投影或者UTM山区项目我更倾向于用当地的中央经线投影变形小。精度这块容易被忽略。DEM的垂直精度哪怕只有几米的误差对坡度的影响在平坦地区会被急剧放大。平地上差2米可能就把2度的坡算成8度。所以做项目前最好收集一些实测高程控制点来验证别迷信所谓的公开数据。1.3 公开DEM数据源怎么选我的建议目前国内能用的公开DEM主要有SRTM、ASTER GDEM、ALOS AW3D30以及国内的天地图有时也会提供地形服务注用于配合在线底图展示用数据质检和投影转换还是建议拿到本地做。SRTM覆盖全球分辨率约30米1弧秒垂直精度相对均衡适合大区域水文分析和地形宏观特征提取。ASTER GDEM分辨率也是30米但局部区域存在一些异常坑洞和条纹做精细的坡度分析前要先做平滑和填洼处理。ALOS AW3D30同样30米但在山区表现更细腻很多做山区项目的同行喜欢用它。如果条件允许用无人机航测或LiDAR生成的高分辨率DEM做小范围精细地形分析效果是公开数据完全比不了的。只是数据处理量大流程长需要权衡项目成本。2. 坡度与坡向人人都会做但你真的知道它算的是什么吗坡度坡向是地形分析里出场率最高的两个工具操作起来就是点一下的事但很多人并不清楚结果背后的算法逻辑更不知道不同的算法之间竟然能差出好几度。2.1 坡度算法之间的差别ArcGIS和QGIS到底谁准ArcGIS的坡度工具有两种输出单位度和百分比。度的取值范围是0到90百分比是0到无穷大垂直陡崖接近无限。在做土地利用评价、农业适宜性分析时我一般用度在做地表径流和土壤侵蚀相关的研究时百分比更常见因为百分比直接反映“水平走100米上升了多少米”和径流速度、侵蚀能量的物理含义对应得上。算法层面ArcGIS默认用的是一种基于3x3窗口的三阶反距离平方权法也叫Horn算法。它会给中心格网周围8个格网分配权重距离越近权重越大再通过南北和东西方向的高程差分算出坡度。QGIS里如果选“默认”或“Horn”结果基本一致但如果你在QGIS里选了其他算法比如简单的一阶差分结果就会有偏差。我前面提到的同一片地ArcGIS算出15度、QGIS算出17度的例子大概率就是两个软件用的默认算法不一样或者坐标系没对齐。解决的方法很简单做之前固定好投影坐标系然后在两边选择相同的坡度算法再对比结果。2.2 坡向的环形数据陷阱坡向输出的是0到360度的方向角0度或360度表示正北90度表示正东以此类推。坑在于0度和360度本质上是同一个方向但在栅格数据里它们是两个不同的数值。如果你对坡向做重分类比如想把“北坡”划为337.5度到22.5度区间直方图上北坡的像元会被劈成两半一个靠近0度一个靠近359度统计结果就会莫名其妙少了一大截。处理办法有两个一是用三角函数对坡向做变换比如计算cos(坡向)和sin(坡向)把环形数据转成线性。二是重分类时先把坡向减1让北坡落在整数区间中间再按新范围分类。做山地生态研究时北坡南坡的日照、水分差异巨大这个细节会直接影响后续的植被适宜性结论。2.3 一个看实际情况的案例坡度坡向的组合用法光是坡度或者光是坡向信息量都有限组合起来才是地形分析的高频场景。举个例子某个区域想做茶园选址茶树怕涝喜荫但又不能完全无光照。这时候我们可以做三件事生成坡度图筛选出坡度小于25度太陡不方便采摘和保水生成坡向图优先选北坡和东坡日照温和、蒸发少再叠加海拔带避开高海拔冻害区。三个条件叠起来做栅格计算几分钟就能出一张候选区分布图。这类分析在农业区划、林地方向规划、山区道路选线里都非常常用。关键是阈值不能拍脑袋定要结合当地作物的实际习性、地方规程和现场踏勘结果来调整。3. 山体阴影与视域分析让地形“开口说话”的关键技术如果说坡度坡向是用数字描述地形那山体阴影和视域分析就是用“光和视线”去理解地形。这两项技术在制图和工程选址里都很实用但参数设置不当结果会完全不一样。3.1 山体阴影的参数与制图搭配山体阴影Hillshade模拟给定太阳方位角和高度角时地表的光照与阴影分布。它的输出灰度值范围是0到2550代表全黑阴影255代表全亮。ArcGIS和QGIS里都要设置两个参数太阳方位角默认315度即西北方向和太阳高度角默认45度。制图时一个常见的做法是把山体阴影放到DEM或坡度图下面图层混合模式调成“正片叠底”或“叠加”透明度设到50%左右这样既能让地形有立体感又不会遮盖底下的专题颜色。ArcGIS里这个技巧非常实用出图效果好学过之后基本就离不开。但要注意山体阴影的视觉感受高度依赖光照方向。如果整片区域的山脊走向恰好和光线方向平行山脊可能看起来不立体或者阴影被过度拉长。这时候可以多试几个方位角找出最有立体感的组合。做峡谷或山地制图的时候我一般会生成4个不同角度的山体阴影再在权重合成后使用能有效避免单一方向光照造成的地形细节丢失。3.2 视域分析不是所有看到的地方都能看到视域分析Viewshed解决的问题是站在某个观察点上周围哪些地方能被看到。听起来简单实际上参数很多。基础参数包括观察点高度OFFSETA、目标点高度OFFSETB、观察半径内半径和外半径以及地形曲率是否参与计算。默认状态下OFFSETA和OFFSETB都是0含义是站在地表裸地上看地面上的一个点这在现实里几乎不存在。做基站规划时OFFSETA要填入天线挂高比如35米做景观视廊评估时你关心的是建筑上部和山脊线的通视OFFSETA可能就要填楼层高度。这里有几个实用经验。第一观察半径不要设太大复杂山地地形下超出一定距离的“可见性”计算误差急剧上升因为大气折射、植被遮挡这些因素都没法被简单的通视算法考虑进去。第二栅格分辨率对视域结果影响很大分辨率越粗“能看见”的像元越容易被高估尤其是细长形的地物。第三如果你有建筑物或植被高度数据应该预先加到DEM上形成数字表面模型DSM再做视域分析否则算出来的是“地理上可見”而不是“实际上可見”。3.3 通视分析的工程价值通信铁塔选址、森林防火瞭望塔选址、景观视廊控制本质都是视域分析问题。做瞭望塔选址时我通常的做法是先生成多种候选点位的视域范围图再统计每个塔能覆盖多少重点防护面积最后用栅格叠加的方式选出覆盖率最高的一组点位组合着手优化点位布局。另外一个容易忽略的点是视域分析的结果受观察点高度影响极大。同样一个山头塔高30米和60米可视范围可能相差近一倍。所以做选址评估时一定要把塔高当作变量去跑多个场景再对比覆盖差异而不是拿着一个固定高度就出结论。4. 水文分析从DEM自动提取流域的完整链路地形分析里最复杂的、也最容易翻车的就是水文分析。很多人在平缓地区提取河道提取出来的河网七扭八歪或者干脆断成一截一截的。要搞明白为什么得从填洼、流向、流量累积这几个基础步骤说起。4.1 填洼Fill看似简单却是成败关键天然DEM里有很多洼地有的是真实的比如喀斯特漏斗、冰川湖盆有的则是数据误差造成的伪洼地。水文分析想得到连续的水流路径就必须把它认为的“洼地”填平。ArcGIS里的填洼工具就是把低于周围格网的洼地高程抬升到水流能流出去的水平。坑就在这填洼是无差别处理的它分不清真洼地还是伪洼地。你要是跑到四川、贵州那种喀斯特地貌区用默认参数填洼真实的天坑、溶蚀漏斗会被一个不剩地填平提取出来的河网就和实际水系完全对不上。处理办法是控制填洼阈值或者使用“有条件填洼”的工具比如只填深度小于某个值的洼地保留真实的深洼地。ArcGIS的Fill工具可以设置z_limit参数QGIS里也有类似功能。我在喀斯特地区做水文分析一般会先把填洼深度栅格算出来人工看一眼哪些区域是高填洼聚集区再决定要不要手动处理。4.2 流向与流量累积两条核心算法搞清楚就赢了一半流向Flow Direction目前最常用的是D8算法原理是假设每个格网的水只能流向周围8个格网中坡度下降最陡的那个方向用2的幂次编码记录1、2、4、8、16、32、64、128。D8算法简单高效适合坡面流和河道提取但在平缓地区和地形平坦的冲积平原上它会产生大量平行流动河网提取效果差。更精细的还有D∞算法它允许水流在相邻格网之间按比例分流适合模拟扩散性质的地表漫流但计算量大日常项目用得少。流量累积Flow Accumulation是把流向结果累加得到每个格网汇入的水量准确说是“积累的上游格网数量”。流量累积值高的格网就是潜在的河道位置。很多新手直接用默认阈值提取河网成果稀疏或密集得离谱原因就在没去仔细看流量累积的数值分布。4.3 河网阈值怎么定、分水岭怎么出提取河网时一个关键参数是“累积流量阈值”——当像元的流量累积超过这个值才把它归为河道。阈值低河网密集、支流多阈值高河网稀疏、只剩主水道。确定阈值没有万能公式我常用的方法是做一个流量累积值的双对数分布图找曲线拐点。拐点附近的流量值就是比较合理的阈值起始参考。如果时间紧直接试几个阈值比如总像元数的0.5%、1%、2%然后叠加遥感影像和高精度地形对比提取出的河道走向和实际河流的吻合度。分水岭Watershed则是在流向和已提取河网的基础上把每条河道的汇水区边界划分出来。操作时关键一步是把河网栅格按连接性等级矢量化栅格河网链接再以这些链接的每个交汇点或出水口为“倾泻点”生成对应的子流域多边形。4.4 提取完河网还要补哪些功课不少人在提完河网后直接就拿去用等到叠加实地照片才发现河网位置偏差很大。原因很简单30米分辨率的DEM提取出的河网本质是一条条“格网路径”不是真实河流的精确中心线。尤其是在沟谷狭窄、地形陡峭的山区提取的河道可能偏到山腰上。所以拿到河网后我通常会做两件事一是对照高分辨率影像人工核查关键河段必要时手动修改二是做平滑和简化用消除伪节点、简化错误的算法让河网更干净。做流域面积统计时还要特别注意边界是否闭合河网提取不完整常常导致流域面积明显偏小。5. 等高线与地形剖面传统制图手段在现代项目中的新玩法说到地形分析很多人脑子里全是栅格计算却忽略了等高线和剖面图这些“老工具”。实际上等高线生成和剖面分析在今天依然是地形表达和工程设计中不可替代的环节。5.1 等高线生成参数与平滑技巧从DEM生成等高线在ArcGIS里非常快但直接生成的等高线往往呈锯齿状如果不做处理打印出来后很难看。锯齿的来源还是栅格结构——等高线在像元边缘转弯时只能走直角或对角线。处理建议生成等高线后先做一次简化Simplify再选用合适的平滑算法比如PAEK算法做适度平滑。平滑容差根据地形的破碎程度调整山地一般用10到30米平原地区甚至可以更小。平滑太狠会把微地形特征抹掉太平缓的地方等高线过于密集又会让图面乱成一团。实际工作中我还会把等高线按高程分段设置不同的颜色和线宽比如低海拔区域用浅棕色细线高海拔区域用深棕色粗线这样出图层次分明比单纯一条黑线好看很多。5.2 剖面线从A点到B点的地形起伏密码剖面分析是沿着你画的一条线生成地形起伏的二维断面图。它的应用场景太多了进行输水管道和道路选线时用来判断挖填方量进行景观规划时用来分析视线遮挡进行山地运动路径设计时用来评估坡度变化。操作非常简单ArcGIS里用3D Analyst或Spatial Analyst的剖面工具先画一条线然后基于DEM生成剖面图。但这个简单的过程里有两个容易忽略的点。第一剖面线的位置和长度直接影响分析结论线放在不同位置可能的挖填方量差异巨大第二DEM分辨率越粗剖面细节越少山脊和山谷容易在剖面上被“磨平”。所以做精细剖面分析时尽量使用高分辨率DEM。5.3 剖面分析在工程规划里的实际应用拿道路选线举例。两个山头之间修一条连接路最省成本的路线未必是直线。借助剖面分析可以沿着几条候选路线分别生成剖面图比较最大纵坡、累计爬升高度和挖填方量。哪条线路上山段短、坡度均匀哪条线路能避开陡崖在剖面上看得一清二楚。我做过一个景区游步道规划项目当时提了三条候选线路凭感觉根本没法判断哪条走起来舒服。生成剖面图后一条线在400米距离内有连续3段超过25度的陡坡立刻就被淘汰了。剩下两条再结合景观视域分析选出视廊最好的一条。这就是“老工具”和新分析结合的价值。6. 那些年我们一起踩过的地形分析坑地形分析工具用多了总会遇到一些让人抓狂的问题。这一节把我踩过、也看别人踩过的常见坑集中说一下大概率能帮你省下几个晚上的排查时间。6.1 Web墨卡托下做地形分析结果基本没法看很多从在线地图上抓下来的DEM或高程服务默认都是Web墨卡托投影EPSG:3857或经纬度坐标EPSG:4326。在这个坐标下直接做坡度、坡向结果是失真的因为坡度计算依赖水平距离而经纬度网在不同纬度上的实际距离是不一样的。我见过有人拿EPSG:4326的DEM直接提取河网提取出来的水系像被揉碎的纸条几乎无法使用。正确的做法是先把DEM投影到当地合适的投影坐标系比如高斯-克吕格CGCS2000 3度带或6度带或UTM对应分带再做后续分析。做完分析如果要叠加在线地图显示结果转回地理坐标或用自带的动态投影机制显示分析本身必须在平面坐标下完成。6.2 填洼过度、河网断裂这些典型问题怎么排查填洼过度通常表现为提取的河网走向和自然地形明显矛盾或者原来有湖泊、坑塘的区域被完全填平河网变成一条光秃秃的直线。遇到这种情况先回看填洼深度分布如果某个区域填洼深度特别大就要考虑是不是把真实洼地也填了。河网断裂的常见原因有两个一是流向分析前没填洼水流被伪洼地截断二是河网提取阈值太高细支流被过滤掉导致主河道在某些位置“断流”。排查思路其实很简单在ArcGIS里把填洼前后的DEM做一下对比再逐级降低阈值试提取看河网是在哪个阈值下开始变连贯的。6.3 软件之间算法差异导致结论对不上怎么办ArcGIS、QGIS、SAGA GIS等软件在坡度、流向、汇流累积等算法实现上存在差异。比如SAGA GIS在提取水文分析时有多种流向算法可选结果自然不同QGIS某些版本在处理大范围DEM时默认会用到不同的重采样方式也会让高程统计出现微差别。这种差异不一定是谁错了而是算法实现细节不同。处理思路项目开始前先统一软件、统一算法、统一坐标系并在报告里写清楚数据处理流程。不同软件的结果如果需要对比也必须保证使用了相同的表面建模方法、相同的填洼阈值、相同的流向算法。否则互相之间的差异根本说明不了任何问题。6.4 出图和文件管理的几个建议分析做完不算完成果能不能交出去很重要。做地形分析出图我一般会保存两类工程文件一类是ArcGIS里的地图文档/工程文件另一类是导出的高分辨率图片和图层包。图层包的好处是可以把DEM、坡度图、河网矢量等一起打包发送给别人后对方不需要有原始数据也能打开查看避免“拷贝的GIS地图包跑到C盘”这种默认存储路径问题——在软件的选项设置里把默认保存路径改成项目目录能省不少事。如果要把地形分析结果转给CAD用户记得先投影到目标坐标系再用“导出为CAD”功能并在转换前检查单位和坐标偏移否则到CAD里容易出现坐标对不上的尴尬局面。面积平差、图例换行、新建shp文件这些基础操作平时没事多过几遍真正做项目的时候效率会高很多。我自己做地形分析项目时最后一步永远是回到现场或影像上做验证。再精细的分析没有实地佐证都只能算预判。地形分析的价值恰恰在于它能帮我们在大范围区域内快速圈出重点、排除风险而不是替代我们做最终决策。把数据准备、算法选择、阈值调参和成果验证这几个环节都做扎实了地形分析才能真正成为项目里值得信赖的工具。