【Cesium】Cesium实现淹没分析（CallbackProperty方式）

一、淹没分析简介

淹没分析是指根据某区域的地形以及洪水流量速度，动态模拟该地形区域水位上涨过程。**该功能适用于地形起伏较大的区域。**目前Cesium实现淹没分析功能主要利用polygon的extruedHeight属性，在地形区域拾取几个点划定洪水淹没区，并规定洪水上涨限定高度，洪水从划定区域的最低处以polygon的方式逐渐上涨，达到限定高度后停止上涨，淹没过程结束。

二、淹没分析实现

2.1 淹没核心代码

//淹没分析
function floodAnalysis(positions,targertWaterHeight,waterHeight){console.log(targertWaterHeight);console.log(waterHeight);viewer.entities.add({polygon: {hierarchy: new Cesium.PolygonHierarchy(positions),extrudedHeight: new Cesium.CallbackProperty(function(){waterHeight += 0.5;if (waterHeight > targertWaterHeight) {waterHeight = targertWaterHeight;}return waterHeight;}, false),material: new Cesium.Color.fromBytes(64, 157, 253, 150),perPositionHeight: true,}})
}

其中positions为限定区域端点坐标数组，targetWaterHeight为洪水上涨限定高度，waterHeight为洪水起始高度。这里通过CallbackProperty来动态拉伸polygon。CallbackProperty是一个值，其值通过回调函数来决定，从它被创建开始，它就会不断地调用回调函数。

2.2 获取范围内地形高程极值

        由于不确定地形高程的最低点，若extruedHeight从0开始增加，则很可能从地底处开始上涨（因为地形高度一般大于0），所以要确定范围地形的最低点，才能很好的实现淹没分析效果
        目前Cesium没有提供相关方法来获取范围地形的高度极值，参考
Cesium获取绘制范围内地形高程进行淹没分析
这篇文章提供的思路，先构建范围地形的外接矩形，再构建矩形范围的格网，获取每个格网点处的高程，即可粗略计算洪水范围内地形的极值(采用了第三方库turf.js来构建网格)：

//创建polygon外接矩形，并生成点格网，返回所有格网点坐标
function buildPolygonGrid(positions){const tempPoints = []for (let i = 0; i < positions.length; i++) {cid;var cartographic = ellipsoid.cartesianToCartographic(positions[i]);var lat = Cesium.Math.toDegrees(cartographic.latitude);var lng = Cesium.Math.toDegrees(cartographic.longitude);tempPoints.push([lng, lat]);}//生成外接矩形var line = turf.lineString(tempPoints);var bbox = turf.bbox(line);var bboxPolygon = turf.bboxPolygon(bbox);var area = turf.area(bboxPolygon);//生成格网//计算网格点之间的距离，尽量保证范围内有1万个左右格网点。var cellSide = Math.sqrt(area/1000000)/100;var options = { units: 'kilometers'};var grid = turf.pointGrid(bbox, cellSide, options);const gridPositions = [];grid.features.forEach(f => {gridPositions .push(Cesium.Cartographic.fromDegrees(f.geometry.coordinates[0], f.geometry.coordinates[1]));})const promise = Cesium.sampleTerrainMostDetailed(terrain, gridPositions);let maxHeight = 0;let minHeight = 10000.0;Promise.resolve(promise).then(function(updatedPositions) {for(let i=0;i<updatedPositions.length;i++){let height = updatedPositions[i].height;//获取格网点处地形高度minHeight = height < minHeight ? height:minHeight;maxHeight = height > maxHeight ? height:maxHeight;}document.getElementById("maxHeight").value = maxHeight;document.getElementById("minHeight").value = minHeight;});     
}

最终效果如下: