在地形分析与工程地质调查中，计算数字地形模型（DTM）的体积（即填挖方量）是一项核心任务。本文将详细介绍如何在 QGIS 环境下利用 SAGA Next Gen 插件中的 Raster Volume 工具快速、准确地完成这一计算。
1. 核心工具介绍
为了计算 DTM 图层的体积，推荐使用 SAGA Next Gen 提供的 Raster Volume 算法。
- 算法名称：
Raster Volume
- 算法 ID：
sagang:rastervolume
- 功能定义：该工具对应 SAGA 软件中的
Grid | Calculus → Grid Volume 模块，专门用于计算数字高程模型（DEM/DTM）在指定基准面（Base Level）上方或下方的累积体积。
2. 操作实操步骤
2.1 定位工具
启动 QGIS，打开右侧的 工具箱（Processing Toolbox），在搜索栏输入 Raster Volume 或 sagang:rastervolume 即可找到该工具。
2.2 参数配置
双击打开工具，重点配置以下参数：
- Grid (栅格层)：选择需要计算的 DTM 图层。
- Base Level (基准面高度)：
- 输入作为参照的高程值。
- 场景举例：若项目仅需关注高程 920 米以上的土方量，则将此处设为 920.0。
- Method (计算方法)：根据工程需求选择不同的体积统计模式：
- [0] Count Only Above Base Level： 只统计基准面以上的体积（填方）。
- [1] Count Only Below Base Level： 只统计基准面以下的体积（挖方，但不取绝对值，可能以负数形式输出）。
- [2] Subtract Volumes Below Base Level： 用基准面以上的体积减去基准面以下的体积（上方‑下方），适合计算净体积。
- [3] Add Volumes Below Base Level： 将上方体积与下方体积的绝对值相加，得到总移动土方量（挖方+填方）。
2.3 执行与结果获取
点击 Run（运行）。计算完成后，结果通常不会直接生成新图层，而是显示在 日志 窗口中。
例如：Grid Volume: Volume: 2136510.319504
提示：除 SAGA 工具外，QGIS 原生工具 Raster surface volume (native:rastersurfacevolume) 同样具备类似功能，可作为备选方案。
3. 关键注意事项与精度分析
在实际操作中，为确保体积计算结果的准确性，需关注以下三点：
3.1 坐标系
必须使用正确的投影坐标系。
投影坐标系是以 “米” 作为单位的。这意味着在 QGIS 软件内部，一个像素的宽度或长度是以米来计算的。
如果使用的是 EPSG:4326 (WGS 84)，单位是 “度”。在度量坐标系下计算体积，结果将是“平方度 * 米”，这在物理上没有任何意义。
3.2 单位
像素尺寸：若像素尺寸为 0.5 * 0.5 m，则单个像元的面积为 0.25 m²，体积计算结果自动以 立方米（m³） 呈现。
3.3 区域限制
Raster Volume 工具默认对整个图层进行计算。如果只需要计算特定区域（如某个施工区）的土方量，建议在运行工具前：
- 使用 裁剪栅格 (Clip Raster by Mask Layer) 工具，利用多边形将 DTM 裁剪为目标区域。
- 将裁剪后的结果作为输入图层。
3.4 结果解读
正值：该区域在基准面以上的体积，即 填方量。
负值：该区域在基准面以下的体积，即 挖方量。
代数和：正负相抵后的 净土方量，正值表示总体需要填方，负值表示总体需要挖方。
更多问题，欢迎留言。转载须注明出处。