同一套数据CASS10.1和7.1算出来土方差了4000方？

这几天，遇到生态圈的小伙伴拿着一套很有代表性实例数据来咨询：同一套数据，计算范围和坐标数据、设计参数都相同，用CASS方格网法来算土方。CASS10.1和7.1算出的结果相差将近4000方。

是软件bug，还是操作错误？

下面是计算结果的截图，数值上看，差别是挺大的。

CASS7.1计算结果

CASS10.1计算结果

填方差=1972.1-2751.9=-779.8（方）

挖方差=24390.7-20477.9=3192.8（方）

❔

发现问题

拿到实例数据之后，按照小伙伴描述的操作步骤，分别用CASS7.1和10.1计算了一遍。得到的结果和上面的截图接近，但也发现了几个比较大的问题：

1

计算范围比原始地形数据范围，相差过大

下图所示，计算范围线，远远小于参与计算的原始地形数据范围。这种情况，会影响方格网的高程点内插值。

2

计算时没有指定方格网起点

可以指定范围线节点（本次测试选择左下角点）。

3

计算范围线附近的空白区域，无实测高程点

下图所示，计算范围线附近，有空白区域中，并没有实测的高程点。这片区域的计算结果，要仔细核对。

❗

解决问题

1

裁剪数据

⭕依据计算范围，向外偏移一定距离（命令: _offset）

⭕裁剪图形

执行命令: PLSJ，将多余的高程点区域删除。

⭕提取高程点数据

执行命令: LINKSJX1，将上述高程点提取生成坐标文件（*.dat）

2

指定方格网起点

在CASS7.1输入命令fgwtf_getpt，依据命令行提示，指定计算范围线左下角点作为起点，输入设计参数进行方格网土方计算。

在CASS10.1中，在命令行输入fgwtf，按下图所示指定起点。

3

检查空白区域的内插高程点

因为范围线附近的空白区域，无实测高程点，在生成的方格网上检查下图的内插高程点（黄色）和附近的高程点，是否匹配。

同样的操作，在7.1和10.1中各操作一次，得到的计算结果相差个位数。

💭

一句话

总结

这个实例很典型，一次踩了三个坑：数据范围、方格网起点，空白区域内插点。最大的坑就是，由于选了过大的计算数据，导致空白区域高程点内插值偏差，最后的计算结果误差，就没法估计了。还是那句话，不管用什么方法算土方，自己要检查关键点的设计高程、内插高程是否正确哈！

- END -