Cesium填挖方计算实战：如何用turf.fs实现地形与模型的通用分析（附完整代码）

Cesium填挖方计算实战：如何用turf.fs实现地形与模型的通用分析（附完整代码）

在三维地理信息系统的工程实践中，填挖方计算是一个绕不开的核心课题。无论是大型基建项目的土方平衡，还是景观设计中的地形微调，精准的体积估算直接关系到成本控制与方案可行性。过去，许多开发者面对Cesium平台时，常常将地形分析和模型分析割裂开来，使用不同的技术栈，这不仅增加了开发复杂度，也带来了数据一致性的隐患。今天，我想分享一套基于turf.fs库的通用解决方案，它能优雅地统一地形表面与三维模型表面的填挖方计算。这套方法的核心，在于将三维空间问题巧妙地转化为二维平面上的几何运算，再通过空间插值与积分思想还原体积，其思路之精妙，在实际项目中已被验证具备极高的可靠性与灵活性。

1. 理解通用填挖方计算的核心挑战与解决思路

在深入代码之前，我们有必要厘清为什么需要一个“通用”的方法。传统的Cesium填挖方计算，大多依赖于SampledPositionProperty或直接对地形瓦片数据进行采样。这种方法在地形上表现尚可，但一旦遇到倾斜摄影模型、BIM模型或人工构建的三维Mesh，就立刻失灵了。因为模型内部没有连续的高程场，无法通过简单的经纬度去采样获取高度值。

核心矛盾在于： 我们需要一个统一的方法，既能处理连续的地形高程数据，又能处理离散的、任意复杂的三维模型表面。

我的解决思路源于一个经典的几何原理：任何三维表面的体积计算，都可以通过计算该表面与一个参考基准面之间所夹空间的体积来获得。关键在于，我们如何高效、准确地对这个不规则空间进行“切片”和“求和”。这里，我引入了计算几何中两个强大的工具：

• 泰森多边形（Voronoi Diagram）：用于将我们关心的区域（一个多边形范围）自动细分为无数个小的、不重叠的单元。每个单元内任意一点到其发生点（我们随机生成或按规则布设的点）的距离都比到其他发生点近。这为我们提供了天然的空间离散化方法。
• 多边形布尔运算（Intersection）：用于处理我们的分析范围与自动生成的泰森多边形之间的空间关系，精确地裁剪出落在分析范围内的每一个不规则小单元。

通过turf.fs这个强大的地理空间分析库，我们可以在二维平面（屏幕坐标或投影坐标）上轻松实现上述操作。对于每一个裁剪后的小多边形，我们将其投影回三维空间，获取其覆盖区域在真实地形或模型表面的多个采样点的高度，进而估算该小柱体的体积。最后，对所有小柱体体积进行累加和分类（填方或挖方），即得到总工程量。

注意：此方法将三维体积计算分解为大量二维多边形操作与三维点采样，其精度取决于平面细分的粒度（泰森多边形的数量）和三维采样点的密度。这是一种在计算效率与精度之间取得平衡的实用策略。

2. 构建实战开发环境与关键工具库

工欲善其事，必先利其器。让我们先搭建好开发环境，并理解几个核心工具。

2.1 项目初始化与依赖安装

假设你已有一个基于Vite、Webpack或类似构建工具的Cesium项目。我们需要安装核心依赖：

npm install cesium turf @turf/random @turf/voronoi @turf/intersect
# 或者使用 yarn
yarn add cesium turf @turf/random @turf/voronoi @turf/intersect

这里简要说明一下这几个turf模块的作用：

• @turf/random: 用于在指定边界框内生成随机点，这些点将作为泰森多边形的发生点。
• @turf/voronoi: 基于生成的点，创建泰森多边形。
• @turf/intersect: 执行多边形之间的求交运算，这是我们裁剪分析范围的关键。

2.2 核心辅助函数剖析

在实现主算法前，需要几个辅助函数来处理Cesium坐标与Turf几何对象之间的转换。这是连接Cesium三维世界与Turf二维平面分析的桥梁。

首先，一个将Cesium屏幕坐标数组(Cartesian2[])转换为Turf多边形对象的函数至关重要。Turf使用GeoJSON格式，其坐标是[经度, 纬度]或[x, y]的数组。

/**
 * 将Cesium Cartesian2数组转换为Turf Polygon几何对象
 * @param positions Cesium屏幕坐标数组
 *