cesium-xs-sdk 1.0.0

package/src/utils/GeoCalcUtil.js

@@ -0,0 +1,95 @@
+/**
+ * 地理计算工具（距离、面积、高程）
+ */
+export class GeoCalcUtil {
+    /**
+     * 计算两点地理坐标距离（米）
+     * @param {Array<number>} point1 [lng, lat, height]
+     * @param {Array<number>} point2 [lng, lat, height]
+     * @returns {number} 距离
+     */
+    static calcDistance(point1, point2) {
+        const coord1 = Cesium.Cartesian3.fromDegrees(...point1);
+        const coord2 = Cesium.Cartesian3.fromDegrees(...point2);
+        return Cesium.Cartesian3.distance(coord1, coord2);
+    }
+
+    /**
+     * 计算多边形面积（平方米）
+     * 使用球面多边形面积公式（基于 Cesium 几何顶点）
+     * @param {Array<Array<number>>} points 地理坐标点集 [lon, lat, height?]
+     * @returns {number} 面积
+     */
+    static calcArea(points) {
+        if (!points || points.length < 3) return 0;
+
+        const positions = points.map(p => Cesium.Cartesian3.fromDegrees(p[0], p[1], p[2] || 0));
+        const polygon = new Cesium.PolygonGeometry({
+            polygonHierarchy: new Cesium.PolygonHierarchy(positions)
+        });
+        const geometry = Cesium.PolygonGeometry.createGeometry(polygon);
+
+        if (!geometry || !geometry.attributes || !geometry.attributes.position) {
+            return 0;
+        }
+
+        const values = geometry.attributes.position.values;
+        const vertexPositions = [];
+        for (let i = 0; i < values.length; i += 3) {
+            vertexPositions.push(new Cesium.Cartesian3(values[i], values[i + 1], values[i + 2]));
+        }
+
+        // 球面多边形面积（基于单位球上的多边形面积公式）
+        const area = GeoCalcUtil.#sphericalPolygonArea(vertexPositions);
+        return Math.abs(area);
+    }
+
+    /**
+     * 单位球面上多边形的有向面积（平方米）
+     * 采用局部切平面投影法：以多边形重心为原点建立切平面，投影后计算二维多边形面积。
+     * 对于常规 GIS 多边形精度足够；极大规模多边形可后续替换为更严格的球面 excess 算法。
+     * @param {Cesium.Cartesian3[]} positions 世界坐标顶点
+     * @returns {number} 面积（平方米）
+     */
+    static #sphericalPolygonArea(positions) {
+        const ellipsoid = Cesium.Ellipsoid.WGS84;
+        const n = positions.length;
+        if (n < 3) return 0;
+
+        // 计算近似重心（归一化到椭球面）
+        let center = new Cesium.Cartesian3(0, 0, 0);
+        for (const p of positions) {
+            center = Cesium.Cartesian3.add(center, p, center);
+        }
+        center = Cesium.Cartesian3.normalize(center, new Cesium.Cartesian3());
+        center = ellipsoid.scaleToGeodeticSurface(center, new Cesium.Cartesian3()) || center;
+
+        // 在重心处建立东北天（ENU）局部坐标系
+        const tangentPlane = new Cesium.EllipsoidTangentPlane(center, ellipsoid);
+
+        // 投影到切平面
+        const projected = positions.map(p => tangentPlane.projectPointOntoPlane(p, new Cesium.Cartesian2()));
+
+        // 计算二维多边形面积（Shoelace 公式）
+        let area2 = 0;
+        for (let i = 0; i < n; i++) {
+            const p1 = projected[i];
+            const p2 = projected[(i + 1) % n];
+            area2 += p1.x * p2.y - p2.x * p1.y;
+        }
+
+        return Math.abs(area2) * 0.5;
+    }
+
+    /**
+     * 获取指定地理坐标的高程
+     * @param {Cesium.Viewer} viewer Cesium 实例
+     * @param {Array<number>} coord [lng, lat]
+     * @returns {Promise<number>} 高程值
+     */
+    static async getHeight(viewer, coord) {
+        const cartographic = Cesium.Cartographic.fromDegrees(coord[0], coord[1]);
+        const results = await Cesium.sampleTerrainMostDetailed(viewer.terrainProvider, [cartographic]);
+        return results[0]?.height ?? 0;
+    }
+}