4. 矿区坡度分级与重分类

坡度栅格生成之后，每个像元都带有一个连续的坡度值（单位通常为度）。这个连续值本身信息量很大，但对矿区工程判断来说并不直观：0.3° 和 1.2° 在数值上不同，在工程意义上却都属于"几乎平坦"；而 24° 和 26° 可能分别落在"可通行"与"需防护"的两侧。重分类（Reclassify）的任务，就是把连续坡度值按矿区地形特征与工程需求切分为若干有明确工程含义的等级，并输出一张分级栅格，再据此制作坡度分级图。

本节的核心不是"点一下重分类工具"，而是三件事：分级阈值从哪里来、分级结果如何验证、分级图如何表达才不误导读者。

一、分级阈值必须由工程用途反推

矿区坡度分级没有全国统一的唯一标准，常见做法是结合《矿山开采设计规范》类文件、道路与排土场设计坡度要求、以及本矿区实测地形统计来确定。一个可参考的示例分级方案（示例参数，需按实际矿区调整）如下：

阈值不是拍脑袋定的。实操中应先对坡度栅格做一次统计：在 QGIS 中打开栅格图层的"属性 → 信息"，查看最小值、最大值与均值；再用"栅格 → 杂项 → 栅格信息"或"处理工具箱 → 栅格分析 → 栅格图层统计"获取分位数。若矿区 80% 的面积集中在 3°–12°，那么把 5° 和 15° 作为切分点就比均匀等分更合理，因为它把主要地形差异区分开了。反之，若矿区包含大面积陡崖，35° 以上的分级就应再细分，否则急坡区会被压缩成一个无法指导施工的大色块。

二、在 QGIS 中执行重分类

推荐使用处理工具箱中的"栅格分析 → 重分类 → 重分类栅格"（Reclassify by Table），它比栅格计算器更可控，且能直接生成整型分级栅格。关键设置如下：

• 输入栅格：上一步生成的坡度栅格（如 slope_deg.tif），确认其单位为度而非百分比。
• 重分类表：逐行填写"最小值、最大值、新值"。注意区间边界约定——QGIS 默认按左闭右开处理，即 [min, max)。因此 5–15 应写为 5、15、2，而 0–5 写为 0、5、1，避免 5° 被重复计入两级。
• 输出数据类型：选择 Int16 或 Byte，分级值是整数编码，不需要浮点。
• 使用 NoData：勾选，使矿区边界外的空白区不参与分级，避免出现"0 级"伪色块。

若使用栅格计算器，等价表达式为嵌套条件，例如：

("slope_deg@1" >= 0  AND "slope_deg@1" < 5)  * 1 +
("slope_deg@1" >= 5  AND "slope_deg@1" < 15) * 2 +
("slope_deg@1" >= 15 AND "slope_deg@1" < 25) * 3 +
("slope_deg@1" >= 25 AND "slope_deg@1" < 35) * 4 +
("slope_deg@1" >= 35) * 5

该写法依赖各条件互斥，若边界处理不当会出现叠加。因此更推荐重分类表方式，边界规则由工具统一保证。

三、分级结果的验证与制图表达

重分类完成后不能直接出图，必须先验证。第一步做面积统计：用"处理工具箱 → 栅格分析 → 栅格图层唯一值报告"统计各等级像元数，换算为面积（像元数 × 像元面积），检查是否与坡度统计的分位数预期一致。若"急坡"面积异常偏大，多半是坡度单位误用了百分比，或 NoData 未生效。第二步做空间抽查：在陡坡与平坡交界处叠加等高线，确认分级边界与地形转折吻合，而不是沿栅格行列出现锯齿状伪边界。

制图表达上，分级栅格应使用"调色板/唯一值"渲染，按等级赋色，并在图例中同时标注等级名称与坡度区间，例如"陡坡 25°–35°"。仅写"4 级"对读图者没有意义。底图可叠加矿界、道路、排土场边界与等高线，使分级色块与工程对象产生对应关系。下图展示了从连续坡度到分级栅格、再到工程判读的完整关系。

最后需要强调一个容易被忽略的细节：坡度分级图是"工程判读图"，不是"地形美观图"。配色应从低到高形成明确的冷暖递进（如绿→黄→橙→红），并在图例中写清区间与工程含义。若矿区存在采空塌陷区、排土场台阶等特殊对象，可在分级图之上叠加其边界，让读者一眼看出哪些工程对象落在陡坡或急坡区。分级完成后，建议将重分类表、像元面积统计与最终图件一并归档，作为后续山体阴影制图与边坡分析的基础数据。