Cesium专栏-填挖方分析（附源码下载）

Cesium 是一款面向三维地球和地图的，世界级的JavaScript开源产品。它提供了基于JavaScript语言的开发包，方便用户快速搭建一款零插件的虚拟地球Web应用，并在性能，精度，渲染质量以及多平台，易用性上都有高质量的保证。

关于填挖方分析

---

从百度百科查找到的关于填方和挖方的定义如下：

填方：填方指的是路基表面高于原地面时，从原地面填筑至路基表面部分的土石体积。
挖方：指的是路基表面低于原地面时，从原地面至路基表面挖去部分的土石体积。

在简书里，也找到一篇 基于GIS的填挖方平衡分析 的文章，并给出了在ArcGIS软件中的实操。

本文着重要介绍的是在Cesium中如何计算填方与挖方，其中会借鉴小专栏里的另外一篇文章 Cesium专栏-剖面分析 中关于地形高度采样的内容。

需要注意的是：本文使用的计算方式只是一种研究行方案，用于拓展思路，并不是最准确的计算方式，请评估后慎用！

具体做法

---

1.定义地形服务，绘制施工区域

• // 在使用地形服务前，请设置好 Cesium.Ion.defaultAccessToken
• var terrainProvider = Cesium.createWorldTerrain({
• requestWaterMask: true,
• requestVertexNormals: true
• });
•  
• // 定义填挖方的基准面高度
• var excavateHeight = 2300; // 开挖高度
• var buryHeight = 6000; // 填埋高度
•  
• // 定义施工区域
• var scope = [
• Cesium.Cartesian3.fromDegrees( 99, 29),
• Cesium.Cartesian3.fromDegrees( 100, 29),
• Cesium.Cartesian3.fromDegrees( 100, 30),
• Cesium.Cartesian3.fromDegrees( 99, 30)
• ];

2.剖分施工区域，计算每个区域的面积与填挖高度

核心思想：计算 开挖/填埋 的 开挖量/填方量 的 核心思想就是 剖分 微积分

• // 设置剖分最小单元 0.01°
• var subdivisionCell = 0.01; // 剖分精度自定义
• // 存储所有的剖分矩形
• var subRectangles = [];
• for (var i = 99; i <= 100; i = i + subdivisionCell) {
• for (var j = 29; j <= 30; j = j + subdivisionCell) {
• var subRectangle = new Cesium.Rectangle(
• Cesium.Math.toRadians(i),
• Cesium.Math.toRadians(j),
• Cesium.Math.toRadians(i + subdivisionCell),
• Cesium.Math.toRadians(j + subdivisionCell)
• );
• subRectangles.push(subRectangle);
• }
• }

扫描二维码关注公众号，回复： 7059436 查看本文章

• // 计算每个矩形的中心点作为这个矩形的代表
• var subRectanglesCenterPoints = [];
• subRectangles.forEach( subRectangle => {
• var centerPoint = Cesium.Cartographic.fromRadians((subRectangle.west + subRectangle.east) / 2, (subRectangle
• .north +
• subRectangle.south) / 2);
• subRectanglesCenterPoints.push(centerPoint);
• });

• // 采样每个中心点到达地表的高度
• var promise = Cesium.sampleTerrainMostDetailed(terrainProvider, subRectanglesCenterPoints);
• Cesium.when(promise, function (updatedPositions) {
• // 所有高度
• var heights = [];
• updatedPositions.forEach( point => {
• heights.push(point.height);
• });
• });

3.计算填挖方

更多详情见下面链接文章：

GIS之家小专栏此文章：Cesium专栏-填挖方分析（附源码下载）

文章提供源码，对本专栏感兴趣的话，可以关注一波