这几天，遇到生态圈的小伙伴拿着一套很有代表性实例数据来咨询：同一套数据，计算范围和坐标数据、设计参数都相同，用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 -
