3. 矿区地形数据生成与质量检查

在矿区水系与构造解译中，DEM（数字高程模型）是所有地形因子的“母数据”。坡度、坡向、曲率、山体阴影并不是彼此独立的四张图，而是从同一份高程栅格派生出的四组地形表达：坡度回答“有多陡”，坡向回答“朝哪边”，曲率回答“坡面是凸是凹”，山体阴影回答“在给定光照下地形呈现什么明暗纹理”。它们共同服务于两个判断——水系沿何处汇聚与流动，构造（断层、线性体、岩层界线）在何处留下地形痕迹。因此本节的重点不是把四个工具各点一遍，而是理解每个因子的物理含义、合理参数，以及如何用它们反过来检验DEM本身的质量。

一、四个地形因子的含义与生成要点

坡度通常以度或百分比表示，QGIS 的“坡度”工具在输出单位中选择度即可。矿区若以缓坡丘陵为主，坡度多集中在 0°～15°；若为露天采场或山区，30° 以上坡面会明显增多。坡向以正北为 0°、顺时针增加，输出为 0°～360°，其中平坦区常被赋为 -1 或 0，需要在符号化时单独处理，否则会出现“平地朝北”的假象。曲率分为总曲率、剖面曲率和平面曲率：剖面曲率影响水流沿坡面的加速与减速，正值（凸坡）利于发散、负值（凹坡）利于汇聚；平面曲率影响水流在水平方向的辐散与收敛，对判断沟谷走向更有用。山体阴影本质是模拟光照，关键参数是方位角与高度角，常用方位角 315°（西北光）、高度角 45°，但矿区若想突出近南北向线性构造，可把方位角调到与构造走向近似垂直的方向，让断层崖在阴影中“跳出来”。

需要强调的是，这四个因子对DEM误差的敏感度不同：坡度对高程噪声最敏感，曲率次之，山体阴影最不敏感但最容易被“条带”或“台阶”暴露。因此质量检查时，应优先看坡度与曲率是否出现异常纹理。

二、DEM空洞与异常值的检查方法

DEM空洞多来自原始测绘数据的无值区、云影遮挡或拼接缝隙，表现为栅格中的 NoData 斑块。在 QGIS 中可用“栅格计算器”或“重分类”把 NoData 单独提取出来，也可用“栅格图层统计”查看有效像元比例。若空洞面积小且位于平缓区，可用“填补 NoData”工具以邻域插值修补；若空洞跨越沟谷或断层带，插值会抹平关键地形，应回到原始数据补测或换用更高精度源，而不是强行填平。

异常值包括“尖刺”（孤立极高/极低像元）和“台阶”（相邻像元高程突变形成规则条带）。尖刺在坡度图上会表现为孤立的高值亮点，在曲率图上更明显；台阶则会在山体阴影中形成平行条纹。检查时可先对DEM做一次“焦点统计”均值滤波，再用原始DEM减去滤波结果得到残差栅格，残差绝对值超过设定阈值（示例可取 3 倍标准差）的像元即为可疑点。这一步只是定位，是否修改要结合影像与实地判断：采坑边缘、陡崖、堆场边坡本身就会产生大残差，不能一律当作错误。

三、把地形因子用于水系与构造解译

水系解译时，坡度与曲率是“找汇水”的组合：在曲率图上，负值（凹形）像元连成的条带往往就是沟谷中心线；再叠加坡度图，沟谷两侧坡度较大、谷底坡度较小，边界更清楚。若把山体阴影调成半透明叠在曲率之上，沟谷的立体感会明显增强。构造解译则更依赖阴影与坡向：断层崖常表现为阴影图中一条明暗突变、走向稳定的线性带，坡向图在该带两侧出现方向反转，说明地形面被错断。此时应回到坡度图确认该带是否对应陡坡，三者一致才可初步判定为构造线性体，避免把道路切坡或采场台阶误判为断层。

一个实用的检查习惯是“因子互检”：如果阴影图上有一条醒目的线性体，但坡度图和曲率图在该位置毫无响应，那它更可能是光照方向造成的视觉假象或数据条带，而不是真实地形。反过来，若坡度图有异常高值带而阴影图平淡，则要怀疑该处存在尖刺噪声。

四、示例参数与结果判读

以某矿区 5 m 分辨率 DEM 为例（示例参数，非实测）：坡度输出单位为度，坡向 0°～360°，曲率采用剖面曲率，山体阴影方位角 315°、高度角 45°、Z 因子 1。生成后先看统计：若坡度最大值接近 90° 且分布零散，多为尖刺；若坡向在平缓区大量堆积于某一角度，说明平地处理不当。再用残差法定位可疑像元，结合正射影像逐点判断。最终交付的应是“经过检查、纹理合理、因子相互印证”的地形因子组，而不是四张各自漂亮却互相矛盾的图。