Cesium与气体扩散模型实战：从原理到可视化实现

1次阅读

没有评论

共计 2204 个字符，预计需要花费 6 分钟才能阅读完成。

在环境科学领域，气体扩散模型（Gas Dispersion Model）被广泛用于模拟有毒或危险气体在空气中的传播过程。典型场景包括：

化工企业泄漏应急响应：快速预测氯气、氨气等危险气体的扩散范围，辅助制定疏散方案
核事故后果评估：模拟放射性物质的扩散路径和浓度分布
城市空气质量监测：分析 PM2.5 等污染物的传输规律

这些应用对可视化提出了特殊要求：需要在高精度地理底图上叠加动态扩散效果，并保持实时交互能力。

原生支持 WGS84 坐标系，完美匹配地理数据
内置地形、影像服务，简化地理场景构建
时间动态系统（TimeDynamicPointCloud）适合时序数据

需要自定义着色器实现特殊效果时更灵活
非地理参考的纯三维场景性能更好
更轻量级的 WebGL 封装

对于气体扩散可视化这种强地理属性的应用，Cesium 通常是更优选择。

扩散浓度计算公式：

C(x,y,z) = \frac{Q}{2\pi u \sigma_y \sigma_z} \exp\left(-\frac{y^2}{2\sigma_y^2}\right) \left[\exp\left(-\frac{(z-H)^2}{2\sigma_z^2}\right) + \exp\left(-\frac{(z+H)^2}{2\sigma_z^2}\right)\right]

实际实现时可简化为：
1. 忽略垂直扩散项（z 轴）
2. 使用查表法预计算 σ 值
3. 将连续模型离散为粒子集合

const particleSystem = viewer.scene.primitives.add(
  new Cesium.ParticleSystem({
    image: 'path/to/particle.png',  // 使用圆形透明渐变贴图
    imageSize: new Cesium.Cartesian2(15, 15), // 像素尺寸
    startScale: 1.0,
    endScale: 0.1,  // 粒子消散时缩小
    lifetime: 10.0, // 粒子存活时间(秒)
    emitter: new Cesium.CircleEmitter(100.0), // 圆形发射区域
    emissionRate: 50.0, // 每秒发射粒子数
    modelMatrix: Cesium.Matrix4.IDENTITY,
    emitterModelMatrix: computeWindMatrix(windDirection) // 风速影响矩阵
  })
);

// 地理坐标转场景坐标
function geoToScene(position: Cartographic) {
  const cartesian = Cesium.Cartesian3.fromRadians(
    position.longitude,
    position.latitude,
    position.height
  );
  return Cesium.Matrix4.fromTranslation(cartesian);
}

// 风速影响矩阵
function computeWindMatrix(windVec: Cartesian3) {
  const rotation = Cesium.Quaternion.fromHeadingPitchRoll(
    new Cesium.HeadingPitchRoll(Math.atan2(windVec.y, windVec.x),
      0,
      0
    )
  );
  return Cesium.Matrix4.fromTranslationRotationScale(
    Cesium.Cartesian3.ZERO,
    rotation,
    Cesium.Cartesian3.ONE
  );
}

根据视距动态调整粒子密度（LOD 技术）
使用视锥体剔除（Frustum Culling）不可见区域的粒子
设置最大粒子数上限（通常 5000-10000）

// 主线程
const worker = new Worker('diffusion-worker.js');
worker.postMessage({originPos: [longitude, latitude],
  windData: currentWindField 
});

// Worker 线程
onmessage = (e) => {const positions = simulateDiffusion(e.data);
  postMessage(positions);
};