3.3 栅格分析与数字地形建模
3.3 栅格分析与数字地形建模
本节摘要：栅格分析是加工链的"煨汤"工序——慢工出深味。本节从 DEM 出发，讲坡度、坡向、汇流等地形因子的计算逻辑，再讲重分类与叠加打分，最终组装成一张场地适宜性评价图。这是国土、规划、新能源项目里复用率最高的分析链。
从 DEM 到适宜性一张图
某新能源公司要在山区县选光伏电站场址，条件有四条：坡度小于 15 度、朝向偏南、远离林区、距道路 500 米以内。这四条没有一条能靠矢量叠加直接解决——坡度与坡向是连续表面，必须走栅格链：DEM 算地形因子，逐条重分类打分，最后加权合成一张适宜性图。
图：栅格地形分析的加工链
一、地形因子：邻域运算的三个代表作
坡度、坡向、汇流累积量是地形分析的三件套，它们共同的计算逻辑是"逐格看邻域"：
# 地形因子计算操作与输出 操作：地形分析菜单 坡度分析，输入 DEM，输出单位选度 输出示例： 坡度范围 0 至 63.4 度 直方图特征：0至8度占全县面积 46%，8至15度占 27%，15度以上占 27% 坡向计算： 输出 8 方位分类（北 东北 东 南东 南 南西 西 北西） 平地格网单独标记（坡度近零处坡向无意义） 汇流累积量： 每格的值 = 汇入该格的上游格数，河道即高汇流带 阈值 5000 格以上者标记为河网，与实测水系比对吻合度 # 检查动作：坡度直方图与县域地貌对照，山区县若 8 度以下占八成， # 先怀疑 DEM 高程基准或格距配置有误
汇流累积量值得多说一句：它是"假想下雨"模拟出来的——每个格子降一格水，水往低处流，下游格子累计上游来水。阈值一调，河网密度就变，所以河网提取的阈值必须用实测水系校准，不是拍脑袋定的。
二、重分类与打分：把连续值翻译成业务语言
业务听不懂"坡度 12.7 度"，只听得懂"满足条件"或"得几分"。重分类就是这个翻译：
# 光伏选址四因子打分表（1 为适宜，0 为排除，中间态给分） 因子 分级规则 得分 权重 坡度 小于8度 3 0.3 8至15度 2 大于15度 0（否决） 坡向 南、南西、南东 3 0.3 东、西 1 北向 0（否决） 林地 非林地 1 0.2 林地 0（否决） 道路可达 距道路小于500米 1 0.2 500至1000米 0.5 大于1000米 0 # 设计要点：坡度坡向是连续改善因子，用等级分； # 林地是红线因子，一票否决，不给中间分
三、地图代数：逐格合成最终答案
打分表定好后，加权合成靠地图代数完成——本质是"多张栅格当多列数据，写一条列间公式"：
# 概念演示：适宜性合成的地图代数表达（与具体平台语法略有出入，逻辑一致） # 各图层为同格网同范围的重分类结果 suitability = ( 0.3 * slope_score # 坡度得分栅格 + 0.3 * aspect_score # 坡向得分栅格（权重同为0.3，见打分表） + 0.2 * forest_score # 林地因子栅格 + 0.2 * road_score # 道路可达栅格 ) # 约束一票否决：任一否决因子为 0 的格子直接清零 suitability = suitability * (slope_score > 0) * (aspect_score > 0) * (forest_score > 0) print("合成完成，理论值域 0 至 2.4") # 输出示例：高适宜区（得分大于2.0）面积 6.8 平方公里，共 17 个斑块 # 其中最大的连续斑块 2.3 平方公里，位于县境东南
两个工程细节决定成败。格网对齐：所有参与合成的栅格必须同格距、同范围、同对齐方式，差一个像素都会导致错位运算，合成前用重采样统一到最粗的那张。权重敏感性：权重是主观判断，改一次权重结论可能翻盘，规范做法是做敏感性检验——权重上下浮动两成，看高适宜区面积变化是否平稳，不平稳就说明结论对权重过敏，需要更多依据。
结果落地：从得分栅格到决策图斑
合成出的得分栅格还不能直接交给决策者——"每格 2.13 分"没人看得懂。落地三步：分级切片把连续得分切成三五档，高适宜、适宜、一般、限制、禁建；矢量化把高适宜档转成面图斑，附上面积与位置属性；成果核算统计各档面积与占比，标出前三大斑块作为备选地块。切片的断点不是等间距乱切，要贴合得分直方图的天然断点，避免把一个连续的山坡切成两半。
# 适宜性成果落地报告骨架（可复用） 一 分级结果 五档面积表 高适宜占比 直方图断点说明 二 备选地块 前三大斑块 位置 面积 主导因子 得分 三 否决清单 被红线因子清零的区域分布 说明否决依据 四 敏感性 权重浮动两成 高适宜面积波动区间 五 数据说明 DEM来源与格距 各因子图层时点 已知局限 # 心法：报告的第五部分最能建立信任 敢写局限的分析才敢被使用
常见问题两则
坡度分析结果与实地感受对不上？ 先查 DEM 的来源与格距——用光滑处理过的地形产品算坡度会系统性偏缓，山地项目尤其明显；再查投影，用未投影的经纬度 DEM 直接算坡度，单位错位会把坡度放大或缩小一个量级。汇流量算出来河网断断续续？ 先查洼地填充是否执行——DEM 里的假洼地会把水流"扣住"，不填洼的汇流结果必然断流；填洼之后如果仍断，查无值区（湖泊边界处的空洞）是否被当成了自由出流口。
# 地形链排障对照（三分钟自查） 症状 先查 次查 坡度整体偏缓 DEM光滑产品 格距过粗 坡度出现异常极值 经纬度未投影 高程单位错（英尺混入） 河网断续 洼地未填充 无值区处理 汇流量量级离谱 出流口设置 流向算法参数 # 四条症状覆盖地形链九成求助工单
⚠️ 常见坑：拿 30 米格距的免费 DEM 做 1 比 1 万尺度的坡度评价。格距粗于所需精度的两倍时，坡度被系统性平滑，陡坡会被算缓——选址结论"安全"得虚假。精度不够时，要么换高精度 DEM，要么明确降级使用并写入报告。
💡 关键直觉：栅格链是"先分后合"——每个因子独立成层、独立打分，最后才合成。任何一步想偷懒把两个因子混在一起算，错误就再也无法定位到具体因子。
本节要点回顾
- 三件套：坡度坡向汇流共享"逐格看邻域"的运算逻辑，直方图是质检窗口
- 阈值校准：河网提取阈值必须用实测水系标定，不能默认
- 打分表：连续因子给等级分，红线因子一票否决，两类规则不能混
- 地图代数：同格网对齐是前提，权重敏感性检验是结论可信度的体检
- 精度上限：格距决定能分辨的地物尺度，先问数据再问算法
连续表面的世界走完，下一节给加工链加上第三根轴与时间轴：三维与时空分析。