0.序

在工程设计中，经常需要根据设计方案和现状地形之间进行方量计算，尤其关注方量变化的区域，哪些区域需要填方，哪些区域需要挖方，并依据此进行方量的平衡。

在流域管理中，尤其是湿地、三角洲等容易淤积或者河岸容易被冲刷的河床需要经常关注水下地形的变化，根据变化制定合理的管理策略。到底是需要进行清淤处理还是河床加固。

本文的重点内容是基于GlobalMapper进行2期地形的变化区域提取，根据区域计算每个区域的方量变化（填方、挖方）。

1.随机地形的生成（可跳过）

示例数据只有一期河底地形，为了模拟地形的变化，需要制作一个变化后的现状地形。那就需要在现状地形的基础上做一些地形的变化

当然地形的生成也经常用于解决倾斜模型、地形的融合。可以自己绘制一个指定范围的地形，也可用过一些离散高程点来生产一个地形数据。

本过程是通过对离散高程点进行随机数的一个变化。

新建一个变化字段：

简单粗暴，就按照随机数生成一个指定范围的值，先正的，再负的。这样填方、挖方都有了。

通过rand函数，生成0-0.1之间的值。

也可以在原来值的基础上加。

生成新的地形如下：

2.地形求差

用新地形和老地形的高程值进行相减，求差。

选择减法

ps：也可用于计算两份地形的最大值、最小值、平均值

结果如下：

3.填方、挖方范围提取

这里采用的方法是根据地形差，大于0的就是需要挖方的区域，小于0的就是需要填方的区域

设定地形的有效值：

提取边界：

挖方范围提取结果如下：

方量计算：

选中所有的面对象

可以把计算结果直接导出为一个csv文件。或者把用于计算方量的范围文件导出为shp。

结果如下图：