尧图网站建设 
简介这是一份围绕 ArcGIS 平台进行土方量计算实用性论证的技术文档面向测绘、工程算量及地理信息从业者旨在解决传统 CASS 方法在三角网调整中费时费力、可视化不足的痛点。文档从计算原理入手说明土方精度依赖数字高程模型DEM的构建质量介绍了 DEM 的常见数据结构并演示了通过两期 DEM 叠加获取填挖分界线、统计开挖与回填体积的计算流程同时讨论了构建 TIN 的两种常见建模方式强调结合实际地形选择合适的插值模式以提升地面仿真度。作者通过多个工程将 ArcGIS 与 CASS 的计算结果进行对比验证了 ArcGIS 在保证准确性的同时具备三维直观、自动化程度高、计算速度快、边界可灵活拆分等实用优势。资源为 1 个 PDF 文件大小约 183KB已有 86 人学习适合需要了解 ArcGIS 土方计算可行性及操作思路的工程技术人员下载参考。1. 土方量计算这件事ArcGIS 的实用性该重新估值《ArcGIS 计算土方的实用性论证》这几年频繁出现在投标方案、技术交底和竣工结算里背景很现实无人机航测、GPS 碎部点、多期地形数据越来越常见传统 CAD 平台算土方开始吃力。ArcGIS 的看家本领是表面建模与空间统计这与土方量的数学定义严丝合缝所以它不是 CASS 的替代品而是一条能互相验证的路径。真正的实用性不在“能不能算”而在“算完能不能交代清楚”三角网模型、栅格高差、体积统计的每一步都能存档回查。两期土方对比、大范围场平、多地块同时核量正好是这套能力的舒适区。测量员、施工与造价技术人员都能照下面的流程走通拿到与 CASS 互相校验的体积数也能看清哪些参数会把结果带偏。2. 从高程点到填挖方ArcGIS 土方计算的原理与数据准备2.1 土方量的本质是两张曲面之间的体积差先立住数学底子。不管断面法、方格网法还是三角网法土方量本质上都是对高差场做积分把设计面或完工面高程减去现状面高程得到每个位置的填挖深度再乘上对应面积并累加。离散到栅格里就是 V Σ(ΔZ × S)S 是单个像元面积。ArcGIS 的 TIN 和栅格 DEM只是同一个积分在两种数据结构上的实现不存在“谁算得准、谁算不准”的玄学只存在离散化方式不同带来的精度差。业务上要区分两种口径。一期土方是“现状面 vs 设计高程面”设计面往往是恒定标高或带坡度的平面两期土方是“施工前实测面 vs 施工后实测面”两期均为外业实测数据正好对应 CASS 两期土方计算里最常出现的场景开工前收原地面完工后收竣工面中间差值就是结算依据。ArcGIS 对两种口径的处理手段一致只是数据来源不同这也是它能被独立论证实用性的前提。2.2 离散点、等高线与 CAD三种数据进 ArcGIS 的路线2.2.1 离散点从 Excel 进入要素类外业碎部点通常是 CSV 或 XLSXArcGIS 里最稳的入口是 XY 表转点。界面操作为 文件 → 添加数据 → XY 点数据但批量处理时不如一句工具命令可控import arcpy arcpy.env.workspace rD:\earthwork\data.gdb arcpy.management.XYTableToPoint( rD:\earthwork\survey_before.csv, survey_before, X, Y, Z, coordinate_systemarcpy.SpatialReference( CGCS2000 / 3-degree Gauss-Kruger CM 108E))这段代码把 CSV 的 X、Y、Z 三列转成点要素类关键参数有三个坐标字段要显式指定不要把默认经纬度当成平面坐标坐标系必须和外业投影一致否则后面所有面积与体积都错Z 字段在这里作为属性写入供建 TIN 时调用。转完先打开属性表扫一遍 Z 值范围确认高程单位是米还是厘米这一步能省掉一整轮排错。2.2.2 等高线与 CAD 转 SHP 的三个动作设计院给的 DWG 常有现成等高线用 ArcGIS 将 CAD 转化为 SHP 文件是常规操作但直接全层导入会带进大量无关实体。我一般只导出 Polyline 图层的等高线并只保留高程字段非空的要素。三个动作缺一不可先确认 CAD 里等高线高程是存在几何 Z 值还是扩展属性里前者可直接用后者需用要素转点再挂接属性再用投影工具把要素统一到土方计算的坐标系再做 TIN避免单位混乱最后明确等高线内插出的表面属于平滑后的近似精度受等高距限制只能用于方案阶段初算不能作结算依据。2.3 TIN 还是栅格两种表面模型的选型边界选型决定后续工具走向。TIN 由实测点直接构成不规则三角网节点高程等于原始观测值栅格 DEM 是规则矩阵每个像元值来自对 TIN 或插值面的重采样。二者各有一个明显的长板直接对比对比项TIN 三角网栅格 DEM数据结构不规则三角面节点密度自适应规则矩阵像元大小固定高程保真度保留全部实测点受像元尺寸与内插方法影响大范围运算面片多速度慢矩阵运算快适合批量体积计算工具Surface Volume、Polygon Volume栅格计算器、Cut/Fill典型土方场景小范围改造、有陡坎断裂线多期对比、场平、跨地块核算工程里做两期土方对比我几乎总是转成栅格计算因为栅格相减和分区统计都是一条命令的事但测区有陡坎、河堤这类线性突变时直接用 TIN 更稳断裂线不会被栅格化过程抹平。2.4 坐标系统与单位Z 因子错一位方量差一个量级土方计算必须在投影坐标系下进行平面单位是米。地理坐标系的度做不了面积积分算出的体积没有任何意义。很多项目拿到影像或 CAD 后不确认坐标系直接开算我的习惯是先检查属性里的线性单位再做定义投影或投影栅格。如果涉及在 ArcGIS 里给影像改坐标系统用 Project Raster 重采样到目标投影对 DEM 重采样方法优先选双线性或三次卷积最邻近会留下块状台阶。Z 值单位是另一个高发坑。CAD 高程常按毫米或厘米绘制点表的 Z 又是米混在一起体积差 10 倍甚至 100 倍。Surface Volume、Cut/Fill、TinRaster 都带 Z factor 参数X/Y 单位米、Z 单位米时取 1.0Z 是厘米时取 0.01Z 是毫米时取 0.001。把 Z factor 和坐标系一起写进脚本注释或日志是这个场景最便宜的保险。3. 用 ArcGIS 跑通两期土方计算的最小流程3.1 第一步从两期点集生成 TIN 与 DEM 栅格流程起点是两次外业测量得到的点要素类 survey_before 与 survey_after。先用 Create TIN 建表面ArcMap 10.2 到 10.8 的界面路径是 ArcToolbox → 3D Analyst Tools → Data Management → TIN → Create TINArcGIS Pro 3.x 的同名工具也在 3D Analyst 工具箱里。脚本方式更可控import arcpy arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.overwriteOutput True sr arcpy.SpatialReference(CGCS2000 / 3-degree Gauss-Kruger CM 108E) arcpy.ddd.CreateTIN( tin_before, sr, [[survey_before, Shape.Z, masspoints], [site_boundary, Shape, hardclip]]) arcpy.ddd.TinRaster( tin_before, dem_before, FLOAT, LINEAR, CELLSIZE 2, 1)CreateTIN 的输入列表每项是 [要素类, 高程字段, 表面要素类型]。高程字段用 Shape.Z 表示直接读三维几何的 Z 值masspoints 把点作为质量点参与构网hardclip 把边界作为硬裁剪边防止三角网外扩。TinRaster 的关键参数里LINEAR 是线性内插CELLSIZE 2 把像元设为 2 米末尾的 1 是 Z factor。同样代码对 survey_after 重复一遍得到 dem_after。像元尺寸不是随手填的它要跟测点密度匹配。粗略经验平均点距 3~5 米时取 1~2 米点距 10 米或更疏时取 2~5 米。像元过小不会带来更多精度只放大点间插值噪声像元过大则平滑真实地形。第 4 章会给定量结论。提示两期 DEM 必须使用同一坐标系、同一 CELLSIZE、同一范围栅格计算器相减时才不会出现条带错位。3.2 第二步栅格计算器做填挖分离并累计体积拿到 dem_before 和 dem_after核心计算只有一行完工面减现状面。界面操作是 Spatial Analyst → Map Algebra → Raster Calculator表达式为dem_after - dem_before。正值表示完工面更高需要填方负值表示完工面更低需要挖方。屏幕上红色区是要挖的、蓝色区是要填的但视觉图不能直接结算体积需要用像元面积乘累计高差import arcpy from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.cellSize dem_before arcpy.env.snapRaster dem_before arcpy.env.mask site_boundary arcpy.env.overwriteOutput True before Raster(dem_before) after Raster(dem_after) diff after - before fill_ras Con(diff 0, diff, 0) # 填方完工面高于现状面 cut_ras Con(diff 0, -diff, 0) # 挖方取正便于累计 zone Con(IsNull(diff), 0, 1) # 全测区单值分区 ZonalStatisticsAsTable(zone, Value, fill_ras, stat_fill, SUM) ZonalStatisticsAsTable(zone, Value, cut_ras, stat_cut, SUM) def read_sum(tbl): with arcpy.da.SearchCursor(tbl, [SUM]) as cur: return [row[0] for row in cur][0] or 0 cell_area 2.0 ** 2 # 像元边长平方单位 m² fill_vol read_sum(stat_fill) * cell_area cut_vol read_sum(stat_cut) * cell_area print(f填方 {fill_vol:.2f} m3挖方 {cut_vol:.2f} m3 f填挖相抵 {fill_vol - cut_vol:.2f} m3)这段脚本有四个环境变量要说明。cellSize 和 snapRaster 都指向基准栅格保证两期 DEM 像元严格对齐否则相减出现错位条纹mask 指向计算边界杜绝矩形外推ZonalStatisticsAsTable 的 SUM 统计的是高差之和乘上像元面积才是体积。填挖相抵为正值表示整体需要借方负值表示整体需要外运这个净方量在土方平衡里常被单独要求。3.3 第三步Cut/Fill 与 Surface Volume 交叉验证栅格计算器是基准结果但不能只靠一条公式对业主交代。工程上习惯再用 3D Analyst 的 Cut/Fill 工具验一遍位置在 ArcToolbox → 3D Analyst Tools → Raster Surface → Cut/Fill输入 before 面与 after 面输出栅格属性表自带 Volume 和 Area 字段。要注意 Cut/Fill 输出值的正负约定在不同版本里并不完全一致多数版本正值对应挖方、负值对应填方稳妥做法是拿它跟 3.2 节的结果做绝对值比对不依赖符号。一期土方现状面 vs 设计高程面用 Surface Volume 更直接工具在 3D Analyst Tools → Functional Surface → Surface Volume输入表面选现状 DEM基准面高程填设计标高。方向选 ABOVE 得到设计标高以上的实体体积也就是需要挖走的方量选 BELOW 得到设计标高以下的空腔体积也就是需要填进来的方量。这个工具的 Z factor 参数和 2.4 节的问题直接挂钩。三套数若相差超过 2%优先查像元尺寸、范围和 NoData 区域而不是怀疑算法。3.4 与 CASS 方格网结果核对的口径CASS 是测绘结算的行业惯例ArcGIS 结果最终要拿到会上和它对账。两者不是两种互相矛盾的方法CASS 方格网法和 ArcGIS 栅格法在数学上同源都是在规则网格交点上取高差、乘网格面积累加差别在网格怎么生成、交点高程怎么内插。核对维度CASS 方格网ArcGIS 栅格网格生成手输网格尺寸 1m/5m/10mCELLSIZE 参数可任意高程获取对三角网内插交点对 TIN 线性内插到像元边界处理方格压界处按面积折减mask 掩膜或按掩膜提取体积累计逐格统计分区统计 SUM 乘像元面积可回溯性图纸与表格栅格与脚本全链路存档对账有三条铁律计算范围与 CASS 收方范围完全一致网格尺寸对齐CASS 用 5 米格网时 CELLSIZE 也取 5 米允许 1%~3% 的偏差超过 5% 先查边界和坐标系而不是改算法参数。把栅格相减的高差图和 CASS 方格填挖值叠在一起偏差集中在哪个区块一目了然这是 ArcGIS 在结算争议里最有说服力的输出。4. 精度与数据污染ArcGIS 土方计算实用化绕不开的三个问题4.1 像元尺寸的敏感性测试体积跟着分辨率漂移栅格法最容易被忽略的参数是 CELLSIZE。同一个 TIN用 0.5、1、2、5、10 米分别转栅格再算体积结果不会相等而且差值有规律像元越小越逼近 TIN 真值但计算量与文件尺寸非线性增长像元太大时山脊沟谷被抹平高差被平均掉填挖体积整体向小偏移陡坡地段更明显。平缓地形从 1 米换到 5 米可能只差 0.5%高差大的场平项目能差到 5% 以上。实操里我习惯做一组敏感性测试循环改 CELLSIZE 并打印体积import arcpy from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.snapRaster dem_before arcpy.env.mask site_boundary for cs in [0.5, 1, 2, 5, 10]: arcpy.env.cellSize cs diff Raster(dem_after) - Raster(dem_before) zone Con(IsNull(diff), 0, 1) ZonalStatisticsAsTable(zone, Value, diff, fin_memory/s_{cs}, SUM) with arcpy.da.SearchCursor(fin_memory/s_{cs}, [SUM]) as cur: s [r[0] for r in cur][0] or 0 vol s * (cs ** 2) print(fCELLSIZE{cs:4} 体积{vol:15.2f} m3)判断标准是收敛性相邻两档体积变化小于 1%就认为该分辨率可接受如果从 1 米到 2 米还在跳 3%说明当前点密度撑不起 1 米网格再加密像元只是渲染出的假精度。记录这组数据的另一个好处是结算报告里可以写“网格尺寸经敏感性测试确定”这句话在审计时很有分量。4.2 边界不闭合范围线没画好方量一定偏大栅格天生是矩形不做掩膜时计算范围会把边界外区域全包进来。边界外像元高程来自 TIN 外推往往是坡脚或低洼处的插值值面积乘高差就是白送出去的土方。最常见的现场是 ArcGIS 算出来比 CASS 多 10%查到最后是边界线没闭合或没参与掩膜。标准做法分三步启动编辑器目标图层选边界面要素沿红线或权属线绘制闭合多边形需要时用追踪工具沿已有线划。保存编辑后用面转线再线转面回检闭合性或者用拓扑规则检查“不能有悬挂点”。进入体积计算前设置 arcpy.env.mask或用 Spatial Analyst → Extraction → Extract by Mask按掩膜提取把两期 DEM 先裁一遍。边界只解决栅格参与范围问题。如果边界取自测绘范围而非设计范围ArcGIS 和 CASS 永远对不上账这一点要在数据准备阶段和甲方书面确认不是技术问题是口径问题。注意内插方法与参数一旦选定必须随成果一起留档否则复核方无法复现你的体积数实用性论证就失去了根基。4.3 内插方法差异TIN、IDW、克里金差多少从离散点生成连续表面不止一种方法TIN 的线性内插、IDW 的反距离加权、克里金的统计最优估计会得到不同 DEM进而产生不同土方量。差异程度取决于地形起伏与点密度平缓地形通常小于 1%山脊沟谷密集区域可达 2%~5%。内插方法适用场景常用参数体积差异经验值TIN 线性内插实测点密、有线性地物断裂线参与构网基准值IDW点分布均匀、局部平稳power2搜索半径可变与 TIN 差 0~2%克里金点稀疏、有明显空间趋势半变异函数模型与 TIN 差 1~5%土方结算需要可复核我推荐以 TIN 为基准因为它的节点就是实测点本身不经过统计假设IDW 和克里金适合方案阶段的快速估算用了就必须把功率、搜索半径、半变异模型写进原始记录。更严格的项目可以同时用 TIN 和克里金各算一套差异落在 2% 内再定案这也是“实用性论证”里最站得住的论据。4.4 坏点、重复点与 Z 值突变外业数据里总有粗差点树顶、车顶、仪器误操作。这些点不剔除TIN 会在局部拉出“火山口”填挖分布图上一片异常红蓝体积失真。第一步做两件规定动作对点要素按 Z 值排序检查最高最低值是否落在测区的合理地形范围内再用删除相同项处理完全重复的记录arcpy.management.DeleteIdentical(survey_before, [Shape, Z])DeleteIdentical 的字段列表只放 Shape 和 Z意思是位置与高程完全相同的记录只保留一条避免同一点被重复加权。第二步更严格用克里金插值做交叉验证把每个点的实测 Z 与预测 Z 相减得到残差残差绝对值超过 3 倍标准差的点单独导出成图层回现场复核或直接剔除。这两道工序做完数据污染问题基本能压到结算争议线以下。5. 从土方量算出来到交出去ArcGIS 批量出表、出图与数据交换的落地技巧5.1 多地块批量土方统计的 Python 模板一个场平项目常有十几个地块逐个开界面跑栅格计算器效率太低。用循环把边界要素名作为参数一次输出全部分区填挖方对照表import arcpy, csv from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.cellSize 2 arcpy.env.snapRaster dem_before arcpy.env.overwriteOutput True def calc(poly_name): arcpy.env.mask poly_name _bd # 各地块边界面 diff Raster(dem_after) - Raster(dem_before) fill_ras Con(diff 0, diff, 0) cut_ras Con(diff 0, -diff, 0) zone Con(IsNull(diff), 0, 1) ZonalStatisticsAsTable(zone, Value, fill_ras, in_memory/f, SUM) ZonalStatisticsAsTable(zone, Value, cut_ras, in_memory/c, SUM) def s(t): with arcpy.da.SearchCursor(t, [SUM]) as cur: return [r[0] for r in cur][0] or 0 return s(in_memory/f) * 4.0, s(in_memory/c) * 4.0 rows [] for dk in [DK01, DK02, DK03]: vf, vc calc(dk) rows.append([dk, round(vf, 2), round(vc, 2)]) with open(rD:\earthwork\volumes.csv, w, newline) as f: w csv.writer(f) w.writerow([地块, 填方_m3, 挖方_m3]) w.writerows(rows)每次循环切换 mask等于用不同边界裁剪两期 DEM像元面积是 2 × 2 4SUM 乘 4.0 即体积。这个模板在 ArcGIS 10.2 到 10.8 及 ArcGIS Pro 3.x 都能跑差别只在 SpatialReference 字符串写法。5.2 填挖方分布图与断面复核体积数字之外结算会议更认图。把 diff 栅格用分类色带渲染负值暖色表示挖方、正值冷色表示填方叠加边界、坐标网格和图例就是一张能直接进报告附件的填挖方分布图。断面复核推荐 ArcGIS Pro 的 Stack Profile 工具在中线上生成剖面把现状线与设计线的高差画出来与 CASS 断面法成果做逐交点对比ArcMap 老用户可用 3D Analyst 的剖面功能。剖面图能暴露“总方量对但局部高差错了”的隐蔽问题这类问题在纯数字汇总里根本看不出来。5.3 成果转出导出 Excel 与转 CAD 时保住精度土方成果最终要进 Excel 表格和 CAD 图纸转出环节最容易掉精度。ArcGIS 导出 Excel 用 Conversion Tools → Table to Excel一步到位导出前先检查小数位数ArcGIS 里高程字段常见的小数点前不显示 0 只是显示格式问题右键字段 → 数字格式设为 0.000 即可正常显示和导出数值本身不受影响。CAD 方向有两个高频问题一是转 CAD 后地类注记和符号全部丢失因为 Export to CAD 只带几何和属性不带符号渲染符号要先转注记或靠种子图模板实现二是三维点转 CAD 多段线时要确认 Z 值写进了实体几何否则到 CASS 里又变回二维点前功尽弃。转出后立刻用 CAD 打开抽查三个标高是最快的验收动作。至于把边界和高程点转成 DWG 再给 CASS 补算一次通常不用改任何参数就能对齐这一步相当于在交付前加了最后一道双保险。本文还有配套的精品资源点击获取 ArcGIS土方量计算实用指南：原理、流程与精度控制
简介这是一份围绕 ArcGIS 平台进行土方量计算实用性论证的技术文档面向测绘、工程算量及地理信息从业者旨在解决传统 CASS 方法在三角网调整中费时费力、可视化不足的痛点。文档从计算原理入手说明土方精度依赖数字高程模型DEM的构建质量介绍了 DEM 的常见数据结构并演示了通过两期 DEM 叠加获取填挖分界线、统计开挖与回填体积的计算流程同时讨论了构建 TIN 的两种常见建模方式强调结合实际地形选择合适的插值模式以提升地面仿真度。作者通过多个工程将 ArcGIS 与 CASS 的计算结果进行对比验证了 ArcGIS 在保证准确性的同时具备三维直观、自动化程度高、计算速度快、边界可灵活拆分等实用优势。资源为 1 个 PDF 文件大小约 183KB已有 86 人学习适合需要了解 ArcGIS 土方计算可行性及操作思路的工程技术人员下载参考。1. 土方量计算这件事ArcGIS 的实用性该重新估值《ArcGIS 计算土方的实用性论证》这几年频繁出现在投标方案、技术交底和竣工结算里背景很现实无人机航测、GPS 碎部点、多期地形数据越来越常见传统 CAD 平台算土方开始吃力。ArcGIS 的看家本领是表面建模与空间统计这与土方量的数学定义严丝合缝所以它不是 CASS 的替代品而是一条能互相验证的路径。真正的实用性不在“能不能算”而在“算完能不能交代清楚”三角网模型、栅格高差、体积统计的每一步都能存档回查。两期土方对比、大范围场平、多地块同时核量正好是这套能力的舒适区。测量员、施工与造价技术人员都能照下面的流程走通拿到与 CASS 互相校验的体积数也能看清哪些参数会把结果带偏。2. 从高程点到填挖方ArcGIS 土方计算的原理与数据准备2.1 土方量的本质是两张曲面之间的体积差先立住数学底子。不管断面法、方格网法还是三角网法土方量本质上都是对高差场做积分把设计面或完工面高程减去现状面高程得到每个位置的填挖深度再乘上对应面积并累加。离散到栅格里就是 V Σ(ΔZ × S)S 是单个像元面积。ArcGIS 的 TIN 和栅格 DEM只是同一个积分在两种数据结构上的实现不存在“谁算得准、谁算不准”的玄学只存在离散化方式不同带来的精度差。业务上要区分两种口径。一期土方是“现状面 vs 设计高程面”设计面往往是恒定标高或带坡度的平面两期土方是“施工前实测面 vs 施工后实测面”两期均为外业实测数据正好对应 CASS 两期土方计算里最常出现的场景开工前收原地面完工后收竣工面中间差值就是结算依据。ArcGIS 对两种口径的处理手段一致只是数据来源不同这也是它能被独立论证实用性的前提。2.2 离散点、等高线与 CAD三种数据进 ArcGIS 的路线2.2.1 离散点从 Excel 进入要素类外业碎部点通常是 CSV 或 XLSXArcGIS 里最稳的入口是 XY 表转点。界面操作为 文件 → 添加数据 → XY 点数据但批量处理时不如一句工具命令可控import arcpy arcpy.env.workspace rD:\earthwork\data.gdb arcpy.management.XYTableToPoint( rD:\earthwork\survey_before.csv, survey_before, X, Y, Z, coordinate_systemarcpy.SpatialReference( CGCS2000 / 3-degree Gauss-Kruger CM 108E))这段代码把 CSV 的 X、Y、Z 三列转成点要素类关键参数有三个坐标字段要显式指定不要把默认经纬度当成平面坐标坐标系必须和外业投影一致否则后面所有面积与体积都错Z 字段在这里作为属性写入供建 TIN 时调用。转完先打开属性表扫一遍 Z 值范围确认高程单位是米还是厘米这一步能省掉一整轮排错。2.2.2 等高线与 CAD 转 SHP 的三个动作设计院给的 DWG 常有现成等高线用 ArcGIS 将 CAD 转化为 SHP 文件是常规操作但直接全层导入会带进大量无关实体。我一般只导出 Polyline 图层的等高线并只保留高程字段非空的要素。三个动作缺一不可先确认 CAD 里等高线高程是存在几何 Z 值还是扩展属性里前者可直接用后者需用要素转点再挂接属性再用投影工具把要素统一到土方计算的坐标系再做 TIN避免单位混乱最后明确等高线内插出的表面属于平滑后的近似精度受等高距限制只能用于方案阶段初算不能作结算依据。2.3 TIN 还是栅格两种表面模型的选型边界选型决定后续工具走向。TIN 由实测点直接构成不规则三角网节点高程等于原始观测值栅格 DEM 是规则矩阵每个像元值来自对 TIN 或插值面的重采样。二者各有一个明显的长板直接对比对比项TIN 三角网栅格 DEM数据结构不规则三角面节点密度自适应规则矩阵像元大小固定高程保真度保留全部实测点受像元尺寸与内插方法影响大范围运算面片多速度慢矩阵运算快适合批量体积计算工具Surface Volume、Polygon Volume栅格计算器、Cut/Fill典型土方场景小范围改造、有陡坎断裂线多期对比、场平、跨地块核算工程里做两期土方对比我几乎总是转成栅格计算因为栅格相减和分区统计都是一条命令的事但测区有陡坎、河堤这类线性突变时直接用 TIN 更稳断裂线不会被栅格化过程抹平。2.4 坐标系统与单位Z 因子错一位方量差一个量级土方计算必须在投影坐标系下进行平面单位是米。地理坐标系的度做不了面积积分算出的体积没有任何意义。很多项目拿到影像或 CAD 后不确认坐标系直接开算我的习惯是先检查属性里的线性单位再做定义投影或投影栅格。如果涉及在 ArcGIS 里给影像改坐标系统用 Project Raster 重采样到目标投影对 DEM 重采样方法优先选双线性或三次卷积最邻近会留下块状台阶。Z 值单位是另一个高发坑。CAD 高程常按毫米或厘米绘制点表的 Z 又是米混在一起体积差 10 倍甚至 100 倍。Surface Volume、Cut/Fill、TinRaster 都带 Z factor 参数X/Y 单位米、Z 单位米时取 1.0Z 是厘米时取 0.01Z 是毫米时取 0.001。把 Z factor 和坐标系一起写进脚本注释或日志是这个场景最便宜的保险。3. 用 ArcGIS 跑通两期土方计算的最小流程3.1 第一步从两期点集生成 TIN 与 DEM 栅格流程起点是两次外业测量得到的点要素类 survey_before 与 survey_after。先用 Create TIN 建表面ArcMap 10.2 到 10.8 的界面路径是 ArcToolbox → 3D Analyst Tools → Data Management → TIN → Create TINArcGIS Pro 3.x 的同名工具也在 3D Analyst 工具箱里。脚本方式更可控import arcpy arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.overwriteOutput True sr arcpy.SpatialReference(CGCS2000 / 3-degree Gauss-Kruger CM 108E) arcpy.ddd.CreateTIN( tin_before, sr, [[survey_before, Shape.Z, masspoints], [site_boundary, Shape, hardclip]]) arcpy.ddd.TinRaster( tin_before, dem_before, FLOAT, LINEAR, CELLSIZE 2, 1)CreateTIN 的输入列表每项是 [要素类, 高程字段, 表面要素类型]。高程字段用 Shape.Z 表示直接读三维几何的 Z 值masspoints 把点作为质量点参与构网hardclip 把边界作为硬裁剪边防止三角网外扩。TinRaster 的关键参数里LINEAR 是线性内插CELLSIZE 2 把像元设为 2 米末尾的 1 是 Z factor。同样代码对 survey_after 重复一遍得到 dem_after。像元尺寸不是随手填的它要跟测点密度匹配。粗略经验平均点距 3~5 米时取 1~2 米点距 10 米或更疏时取 2~5 米。像元过小不会带来更多精度只放大点间插值噪声像元过大则平滑真实地形。第 4 章会给定量结论。提示两期 DEM 必须使用同一坐标系、同一 CELLSIZE、同一范围栅格计算器相减时才不会出现条带错位。3.2 第二步栅格计算器做填挖分离并累计体积拿到 dem_before 和 dem_after核心计算只有一行完工面减现状面。界面操作是 Spatial Analyst → Map Algebra → Raster Calculator表达式为dem_after - dem_before。正值表示完工面更高需要填方负值表示完工面更低需要挖方。屏幕上红色区是要挖的、蓝色区是要填的但视觉图不能直接结算体积需要用像元面积乘累计高差import arcpy from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.cellSize dem_before arcpy.env.snapRaster dem_before arcpy.env.mask site_boundary arcpy.env.overwriteOutput True before Raster(dem_before) after Raster(dem_after) diff after - before fill_ras Con(diff 0, diff, 0) # 填方完工面高于现状面 cut_ras Con(diff 0, -diff, 0) # 挖方取正便于累计 zone Con(IsNull(diff), 0, 1) # 全测区单值分区 ZonalStatisticsAsTable(zone, Value, fill_ras, stat_fill, SUM) ZonalStatisticsAsTable(zone, Value, cut_ras, stat_cut, SUM) def read_sum(tbl): with arcpy.da.SearchCursor(tbl, [SUM]) as cur: return [row[0] for row in cur][0] or 0 cell_area 2.0 ** 2 # 像元边长平方单位 m² fill_vol read_sum(stat_fill) * cell_area cut_vol read_sum(stat_cut) * cell_area print(f填方 {fill_vol:.2f} m3挖方 {cut_vol:.2f} m3 f填挖相抵 {fill_vol - cut_vol:.2f} m3)这段脚本有四个环境变量要说明。cellSize 和 snapRaster 都指向基准栅格保证两期 DEM 像元严格对齐否则相减出现错位条纹mask 指向计算边界杜绝矩形外推ZonalStatisticsAsTable 的 SUM 统计的是高差之和乘上像元面积才是体积。填挖相抵为正值表示整体需要借方负值表示整体需要外运这个净方量在土方平衡里常被单独要求。3.3 第三步Cut/Fill 与 Surface Volume 交叉验证栅格计算器是基准结果但不能只靠一条公式对业主交代。工程上习惯再用 3D Analyst 的 Cut/Fill 工具验一遍位置在 ArcToolbox → 3D Analyst Tools → Raster Surface → Cut/Fill输入 before 面与 after 面输出栅格属性表自带 Volume 和 Area 字段。要注意 Cut/Fill 输出值的正负约定在不同版本里并不完全一致多数版本正值对应挖方、负值对应填方稳妥做法是拿它跟 3.2 节的结果做绝对值比对不依赖符号。一期土方现状面 vs 设计高程面用 Surface Volume 更直接工具在 3D Analyst Tools → Functional Surface → Surface Volume输入表面选现状 DEM基准面高程填设计标高。方向选 ABOVE 得到设计标高以上的实体体积也就是需要挖走的方量选 BELOW 得到设计标高以下的空腔体积也就是需要填进来的方量。这个工具的 Z factor 参数和 2.4 节的问题直接挂钩。三套数若相差超过 2%优先查像元尺寸、范围和 NoData 区域而不是怀疑算法。3.4 与 CASS 方格网结果核对的口径CASS 是测绘结算的行业惯例ArcGIS 结果最终要拿到会上和它对账。两者不是两种互相矛盾的方法CASS 方格网法和 ArcGIS 栅格法在数学上同源都是在规则网格交点上取高差、乘网格面积累加差别在网格怎么生成、交点高程怎么内插。核对维度CASS 方格网ArcGIS 栅格网格生成手输网格尺寸 1m/5m/10mCELLSIZE 参数可任意高程获取对三角网内插交点对 TIN 线性内插到像元边界处理方格压界处按面积折减mask 掩膜或按掩膜提取体积累计逐格统计分区统计 SUM 乘像元面积可回溯性图纸与表格栅格与脚本全链路存档对账有三条铁律计算范围与 CASS 收方范围完全一致网格尺寸对齐CASS 用 5 米格网时 CELLSIZE 也取 5 米允许 1%~3% 的偏差超过 5% 先查边界和坐标系而不是改算法参数。把栅格相减的高差图和 CASS 方格填挖值叠在一起偏差集中在哪个区块一目了然这是 ArcGIS 在结算争议里最有说服力的输出。4. 精度与数据污染ArcGIS 土方计算实用化绕不开的三个问题4.1 像元尺寸的敏感性测试体积跟着分辨率漂移栅格法最容易被忽略的参数是 CELLSIZE。同一个 TIN用 0.5、1、2、5、10 米分别转栅格再算体积结果不会相等而且差值有规律像元越小越逼近 TIN 真值但计算量与文件尺寸非线性增长像元太大时山脊沟谷被抹平高差被平均掉填挖体积整体向小偏移陡坡地段更明显。平缓地形从 1 米换到 5 米可能只差 0.5%高差大的场平项目能差到 5% 以上。实操里我习惯做一组敏感性测试循环改 CELLSIZE 并打印体积import arcpy from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.snapRaster dem_before arcpy.env.mask site_boundary for cs in [0.5, 1, 2, 5, 10]: arcpy.env.cellSize cs diff Raster(dem_after) - Raster(dem_before) zone Con(IsNull(diff), 0, 1) ZonalStatisticsAsTable(zone, Value, diff, fin_memory/s_{cs}, SUM) with arcpy.da.SearchCursor(fin_memory/s_{cs}, [SUM]) as cur: s [r[0] for r in cur][0] or 0 vol s * (cs ** 2) print(fCELLSIZE{cs:4} 体积{vol:15.2f} m3)判断标准是收敛性相邻两档体积变化小于 1%就认为该分辨率可接受如果从 1 米到 2 米还在跳 3%说明当前点密度撑不起 1 米网格再加密像元只是渲染出的假精度。记录这组数据的另一个好处是结算报告里可以写“网格尺寸经敏感性测试确定”这句话在审计时很有分量。4.2 边界不闭合范围线没画好方量一定偏大栅格天生是矩形不做掩膜时计算范围会把边界外区域全包进来。边界外像元高程来自 TIN 外推往往是坡脚或低洼处的插值值面积乘高差就是白送出去的土方。最常见的现场是 ArcGIS 算出来比 CASS 多 10%查到最后是边界线没闭合或没参与掩膜。标准做法分三步启动编辑器目标图层选边界面要素沿红线或权属线绘制闭合多边形需要时用追踪工具沿已有线划。保存编辑后用面转线再线转面回检闭合性或者用拓扑规则检查“不能有悬挂点”。进入体积计算前设置 arcpy.env.mask或用 Spatial Analyst → Extraction → Extract by Mask按掩膜提取把两期 DEM 先裁一遍。边界只解决栅格参与范围问题。如果边界取自测绘范围而非设计范围ArcGIS 和 CASS 永远对不上账这一点要在数据准备阶段和甲方书面确认不是技术问题是口径问题。注意内插方法与参数一旦选定必须随成果一起留档否则复核方无法复现你的体积数实用性论证就失去了根基。4.3 内插方法差异TIN、IDW、克里金差多少从离散点生成连续表面不止一种方法TIN 的线性内插、IDW 的反距离加权、克里金的统计最优估计会得到不同 DEM进而产生不同土方量。差异程度取决于地形起伏与点密度平缓地形通常小于 1%山脊沟谷密集区域可达 2%~5%。内插方法适用场景常用参数体积差异经验值TIN 线性内插实测点密、有线性地物断裂线参与构网基准值IDW点分布均匀、局部平稳power2搜索半径可变与 TIN 差 0~2%克里金点稀疏、有明显空间趋势半变异函数模型与 TIN 差 1~5%土方结算需要可复核我推荐以 TIN 为基准因为它的节点就是实测点本身不经过统计假设IDW 和克里金适合方案阶段的快速估算用了就必须把功率、搜索半径、半变异模型写进原始记录。更严格的项目可以同时用 TIN 和克里金各算一套差异落在 2% 内再定案这也是“实用性论证”里最站得住的论据。4.4 坏点、重复点与 Z 值突变外业数据里总有粗差点树顶、车顶、仪器误操作。这些点不剔除TIN 会在局部拉出“火山口”填挖分布图上一片异常红蓝体积失真。第一步做两件规定动作对点要素按 Z 值排序检查最高最低值是否落在测区的合理地形范围内再用删除相同项处理完全重复的记录arcpy.management.DeleteIdentical(survey_before, [Shape, Z])DeleteIdentical 的字段列表只放 Shape 和 Z意思是位置与高程完全相同的记录只保留一条避免同一点被重复加权。第二步更严格用克里金插值做交叉验证把每个点的实测 Z 与预测 Z 相减得到残差残差绝对值超过 3 倍标准差的点单独导出成图层回现场复核或直接剔除。这两道工序做完数据污染问题基本能压到结算争议线以下。5. 从土方量算出来到交出去ArcGIS 批量出表、出图与数据交换的落地技巧5.1 多地块批量土方统计的 Python 模板一个场平项目常有十几个地块逐个开界面跑栅格计算器效率太低。用循环把边界要素名作为参数一次输出全部分区填挖方对照表import arcpy, csv from arcpy.sa import * arcpy.env.workspace rD:\earthwork\data.gdb arcpy.env.cellSize 2 arcpy.env.snapRaster dem_before arcpy.env.overwriteOutput True def calc(poly_name): arcpy.env.mask poly_name _bd # 各地块边界面 diff Raster(dem_after) - Raster(dem_before) fill_ras Con(diff 0, diff, 0) cut_ras Con(diff 0, -diff, 0) zone Con(IsNull(diff), 0, 1) ZonalStatisticsAsTable(zone, Value, fill_ras, in_memory/f, SUM) ZonalStatisticsAsTable(zone, Value, cut_ras, in_memory/c, SUM) def s(t): with arcpy.da.SearchCursor(t, [SUM]) as cur: return [r[0] for r in cur][0] or 0 return s(in_memory/f) * 4.0, s(in_memory/c) * 4.0 rows [] for dk in [DK01, DK02, DK03]: vf, vc calc(dk) rows.append([dk, round(vf, 2), round(vc, 2)]) with open(rD:\earthwork\volumes.csv, w, newline) as f: w csv.writer(f) w.writerow([地块, 填方_m3, 挖方_m3]) w.writerows(rows)每次循环切换 mask等于用不同边界裁剪两期 DEM像元面积是 2 × 2 4SUM 乘 4.0 即体积。这个模板在 ArcGIS 10.2 到 10.8 及 ArcGIS Pro 3.x 都能跑差别只在 SpatialReference 字符串写法。5.2 填挖方分布图与断面复核体积数字之外结算会议更认图。把 diff 栅格用分类色带渲染负值暖色表示挖方、正值冷色表示填方叠加边界、坐标网格和图例就是一张能直接进报告附件的填挖方分布图。断面复核推荐 ArcGIS Pro 的 Stack Profile 工具在中线上生成剖面把现状线与设计线的高差画出来与 CASS 断面法成果做逐交点对比ArcMap 老用户可用 3D Analyst 的剖面功能。剖面图能暴露“总方量对但局部高差错了”的隐蔽问题这类问题在纯数字汇总里根本看不出来。5.3 成果转出导出 Excel 与转 CAD 时保住精度土方成果最终要进 Excel 表格和 CAD 图纸转出环节最容易掉精度。ArcGIS 导出 Excel 用 Conversion Tools → Table to Excel一步到位导出前先检查小数位数ArcGIS 里高程字段常见的小数点前不显示 0 只是显示格式问题右键字段 → 数字格式设为 0.000 即可正常显示和导出数值本身不受影响。CAD 方向有两个高频问题一是转 CAD 后地类注记和符号全部丢失因为 Export to CAD 只带几何和属性不带符号渲染符号要先转注记或靠种子图模板实现二是三维点转 CAD 多段线时要确认 Z 值写进了实体几何否则到 CASS 里又变回二维点前功尽弃。转出后立刻用 CAD 打开抽查三个标高是最快的验收动作。至于把边界和高程点转成 DWG 再给 CASS 补算一次通常不用改任何参数就能对齐这一步相当于在交付前加了最后一道双保险。本文还有配套的精品资源点击获取