Cesium 已知两点坐标计算方位角

在Cesium中，计算两点之间的方位角通常用于确定从一个点到另一个点的理想直线路径相对于地球表面的朝向。这个角度通常被称为"航向角"或"真北角"。以下是一个简单的步骤来计算方位角：

1. **准备数据**：首先获取两个点的经纬度坐标。例如，设A点经度lonA，纬度latA，B点经度lonB，纬度latB。

2. **转换坐标**：将地理坐标转换为三维笛卡尔坐标，因为Cesium基于这种坐标系统。可以使用`Cesium.Ellipsoid.WGS84.cartographicToCartesian`方法：

   const pointA = Cesium.Ellipsoid.WGS84.cartographicToCartesian(new Cesium.Cartographic(lonA, latA, 0));
   const pointB = Cesium.Ellipsoid.WGS84.cartographicToCartesian(new Cesium.Cartographic(lonB, latB, 0));

3. **计算矢量**：通过减法得到从A点到B点的向量：

   const vectorAB = Cesium.Cartesian3.subtract(pointB, pointA, new Cesium.Cartesian3());

4. **找到地磁北向量**：由于地球上存在磁偏角，Cesium提供了一个方法`Cesium.Math.geodeticSurfaceNormal`来计算地心引力作用下的正常向量，这可以近似看作是地磁北向量：

   const geodeticNormal = Cesium.Ellipsoid.WGS84.geodeticSurfaceNormal(latA, lonA);

5. **计算方位角**：最后，通过向量A-B与地磁北向量的内积除以向量A-B的长度，然后取反正切值（注意结果可能是负数，所以加上π如果小于0），就可以得到方位角：

   const angle = Math.atan2(Cesium.Cartesian3.dot(geodeticNormal, vectorAB), Cesium.Cartesian3.magnitude(vectorAB)) % (2 * Math.PI);
   const bearing = angle * (180 / Math.PI);