CesiumCesium实现淹没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;
     });     
}

最终效果如下:

这篇好文章是转载于：学新通技术网