Cesium JS椭球切平面计算与法向应用问题咨询

核心问题原因 #

1. 法向计算逻辑错误：大地法向（geodetic normal）才是椭球表面切平面的正确法向，地心法向为地心到点的连线方向，仅在赤道、极点位置与大地法向重合，你得到的两个法向完全一致说明你误将地心法向当成了大地法向使用。
2. Cesium.EllipsoidTangentPlane参数错误：该构造函数要求传入的原点为无高度的椭球表面笛卡尔坐标，如果传入带高程的点击坐标、或坐标转换逻辑错误，就会出现平面悬浮、不相切的问题。

正确实现步骤 #

1. 预处理点击坐标，得到椭球表面点 #

// clickCartesian为点击事件获取的三维笛卡尔坐标
const cartographic = Cesium.Cartographic.fromCartesian(clickCartesian);
// 重置高度为0，确保坐标落在WGS84椭球表面
cartographic.height = 0;
const surfaceCartesian = Cesium.Ellipsoid.WGS84.cartographicToCartesian(cartographic);

2. 获取正确的大地法向 #

// 调用Cesium内置接口获取椭球表面大地法向，为切平面的法向
const geodeticNormal = Cesium.Ellipsoid.WGS84.geodeticSurfaceNormal(surfaceCartesian);

3. 正确使用EllipsoidTangentPlane实现点投影 #

// 传入椭球表面坐标，初始化切平面
const tangentPlane = new Cesium.EllipsoidTangentPlane(surfaceCartesian, Cesium.Ellipsoid.WGS84);
// 将其他三维世界坐标点投影到切平面，返回平面二维坐标
const projectedPoints = tangentPlane.projectPoints(yourPointList);

4. 切平面可视化实现（如果需要） #

viewer.entities.add({
  // 平面锚点必须为椭球表面坐标
  position: surfaceCartesian,
  plane: {
    // 法向为大地法向，平面到锚点距离为0
    plane: new Cesium.Plane(geodeticNormal, 0),
    // 平面尺寸，单位为米
    dimensions: new Cesium.Cartesian2(200000, 200000),
    material: Cesium.Color.BLUE.withAlpha(0.4)
  }
});

可选：获取东北天局部坐标系的切平面XY轴 #

如果你需要切平面的X轴为东向、Y轴为北向，直接调用坐标转换接口生成局部变换矩阵即可：

// 生成东北天坐标系到世界坐标系的变换矩阵
const enuMatrix = Cesium.Transforms.eastNorthUpToFixedFrame(surfaceCartesian);
// 提取东向、北向、法向向量
const eastAxis = Cesium.Matrix4.getColumn(enuMatrix, 0, new Cesium.Cartesian3());
const northAxis = Cesium.Matrix4.getColumn(enuMatrix, 1, new Cesium.Cartesian3());
const upAxis = Cesium.Matrix4.getColumn(enuMatrix, 2, new Cesium.Cartesian3());

内容的提问来源于stack exchange，提问作者Ilia Shevelev