Cesium地形高度获取：核心方法与实战应用

1. Cesium地形高度获取的核心概念

在三维地理信息系统中，地形高度（Terrain Height）是最基础也是最重要的数据之一。作为一名长期从事Cesium开发的工程师，我经常需要处理各种与地形高度相关的需求。理解地形高度的本质，是掌握后续所有技术方法的前提。

1.1 地形高度的精确定义

Cesium中的地形高度指的是地球表面某一点相对于平均海平面的垂直距离，也就是我们常说的海拔高度。这里需要特别注意三个关键概念的区别：

• 椭球体高度：这是基于WGS84参考椭球体计算的高度，不考虑实际地形起伏。当地形数据未加载时，Cesium默认使用这个值，通常表现为平坦的地球表面。

• 地形高度：这是叠加在椭球体上的真实地形数据，包含了山脉、峡谷等地形特征。只有加载了地形服务（如Cesium World Terrain）后，才能获取到真实的地形高度值。

• 模型高度：这是指人工添加的三维模型（如建筑物、桥梁等）相对于其所处地形的高度。

重要提示：很多初学者容易犯的错误是直接使用Cartographic.height属性，这个值在没有加载地形服务时只会返回椭球体高度。正确的做法是必须通过sampleTerrain或globe.pick等API获取真实地形高度。

1.2 Cesium中的坐标系系统

理解Cesium的坐标系系统是掌握高度获取的关键。Cesium主要使用以下几种坐标系：

1. 经纬度高程坐标系（Cartographic）

   • 包含经度(longitude)、纬度(latitude)和高度(height)三个分量
   • 角度单位为弧度(radians)
   • 高度单位为米(meters)

2. 笛卡尔空间坐标系（Cartesian3）

   • 以地心为原点的三维直角坐标系
   • 常用于空间计算和渲染

3. 屏幕坐标系（Cartesian2）

   • 二维屏幕像素坐标
   • 用于处理鼠标交互等场景

坐标系间的转换方法：

javascript复制// 经纬度转笛卡尔坐标
const cartesian = Cesium.Cartesian3.fromDegrees(lon, lat, height);

// 笛卡尔坐标转经纬度
const cartographic = Cesium.Cartographic.fromCartesian(cartesian);

// 屏幕坐标转拾取射线
const ray = viewer.camera.getPickRay(screenPosition);

2. 地形高度获取的四大核心方法

在实际项目中，根据不同的应用场景，我们需要采用不同的高度获取方法。下面我将详细介绍四种最常用的技术方案。

2.1 基于经纬度的地形采样（sampleTerrain）

这是最基础也是最可靠的高度获取方式，适用于已知经纬度坐标的情况。

2.1.1 基础实现代码

javascript复制/**
 * 获取单个点的地形高度
 * @param {Cesium.Viewer} viewer - Cesium实例
 * @param {number} lon - 经度(度)
 * @param {number} lat - 纬度(度)
 * @returns {Promise<number>} 地形高度(米)
 */
async function getTerrainHeight(viewer, lon, lat) {
  // 转换为弧度制Cartographic
  const position = Cesium.Cartographic.fromDegrees(lon, lat);
  
  // 使用最高精度采样地形
  const [sampledPosition] = await Cesium.sampleTerrainMostDetailed(
    viewer.terrainProvider,
    [position]
  );
  
  return sampledPosition.height;
}

2.1.2 批量查询优化

当需要处理大量点时，直接循环调用上述方法会导致性能问题。我们可以使用批量查询：

javascript复制async function batchGetTerrainHeights(viewer, positions) {
  // 转换为Cartographic数组
  const cartographics = positions.map(pos => 
    Cesium.Cartographic.fromDegrees(pos.lon, pos.lat)
  );
  
  // 批量采样
  const results = await Cesium.sampleTerrainMostDetailed(
    viewer.terrainProvider,
    cartographics
  );
  
  // 返回带原始坐标的结果
  return results.map((result, index) => ({
    ...positions[index],
    height: result.height
  }));
}

性能提示：对于超过100个点的批量查询，建议采用分批次处理（如每批50个点），并在批次间添加短暂延时，避免阻塞主线程。

2.2 屏幕坐标拾取（globe.pick）

在交互式应用中，我们经常需要获取鼠标位置对应的地形高度。

2.2.1 基础实现

javascript复制function getHeightFromScreen(viewer, screenPosition) {
  // 生成拾取射线
  const ray = viewer.camera.getPickRay(screenPosition);
  if (!ray) return null;
  
  // 拾取地形
  const hitPosition = viewer.scene.globe.pick(ray, viewer.scene);
  if (!hitPosition) return null;
  
  // 转换为高度
  return Cesium.Cartographic.fromCartesian(hitPosition).height;
}

// 鼠标点击示例
viewer.screenSpaceEventHandler.setInputAction(event => {
  const height = getHeightFromScreen(viewer, event.po