@math.gl/geospatial 5.0.0-alpha.19 → 5.0.0-alpha.20

package/dist/ellipsoid-tangent-plane.js

@@ -2,66 +2,214 @@
 // SPDX-License-Identifier: MIT AND Apache-2.0
 // SPDX-FileCopyrightText: Copyright 2011-2018 CesiumJS Contributors
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Derived from Cesium. See the repository LICENSE for upstream attribution and Apache-2.0 terms.
+// SPDX-FileComment: Derived from Cesium. Original math.gl follow-up adds explicit radial/normal inverses, packed projection and an orthonormal three-radius tangent basis. See the repository LICENSE for upstream attribution and Apache-2.0 terms.
 // This file is derived from the Cesium library under Apache 2 license
-// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
-import { Vector2, Vector3, Matrix4 } from '@math.gl/core';
+// See ../LICENSE for the retained Cesium notice and Apache-2.0 terms.
+import { Vector2, Vector3 } from '@math.gl/core';
 import { Plane } from '@math.gl/geometry/culling';
 import { Ellipsoid } from "./ellipsoid.js";
-const scratchOrigin = new Vector3();
-const scratchCart3 = new Vector3();
-const scratchEastNorthUp = new Matrix4();
-const scratchProjectPointOntoPlaneCartesian3 = new Vector3();
+import { compensatedDot3, exactNormalIntersection } from "./ellipsoid-helpers/tangent-plane-numerics.js";
 /** A two-dimensional east-north plane tangent to an ellipsoid. */
 export class EllipsoidTangentPlane {
     /**
-     * Creates a new plane tangent to the WGS84 ellipsoid at the provided origin.
+     * Creates a new plane tangent to an ellipsoid at the provided origin.
      * If origin is not on the surface of the ellipsoid, its surface projection will be used.
      *
      * @param origin A nonzero Cartesian point on or near the ellipsoid.
      * @param ellipsoid The ellipsoid to which the origin belongs. Defaults to WGS84.
      */
     constructor(origin, ellipsoid = Ellipsoid.WGS84) {
-        const originOnSurface = ellipsoid.scaleToGeodeticSurface(origin, scratchOrigin);
-        const eastNorthUp = ellipsoid.eastNorthUpToFixedFrame(originOnSurface, scratchEastNorthUp);
-        this._origin = new Vector3(originOnSurface);
-        this._xAxis = new Vector3(scratchCart3.from(eastNorthUp.getColumn(0)));
-        this._yAxis = new Vector3(scratchCart3.from(eastNorthUp.getColumn(1)));
-        const normal = new Vector3(scratchCart3.from(eastNorthUp.getColumn(2)));
-        this._plane = new Plane().fromPointNormal(this._origin, normal);
+        finiteInput(origin, 3);
+        this.radii = Object.freeze(Array.from(ellipsoid.radii));
+        if (this.radii[0] === this.radii[1] && this.radii[1] === this.radii[2])
+            this.radius = this.radii[0];
+        if (!this.radii.every(radius => radius > 0 && Number.isFinite(radius))) {
+            throw new RangeError('Tangent plane requires positive finite ellipsoid radii');
+        }
+        const originOnSurface = ellipsoid.scaleToGeodeticSurface(Array.from(origin), new Vector3());
+        if (!originOnSurface)
+            throw new RangeError('Origin has no unique ellipsoid surface point');
+        this.anchor = new Vector3(originOnSurface);
+        const normal = new Vector3(ellipsoid.geodeticSurfaceNormal(originOnSurface));
+        if (!normal.every(Number.isFinite) || !(normal.len() > 0)) {
+            throw new RangeError('Cannot resolve tangent-plane normal');
+        }
+        normal.normalize();
+        // Use the surface normal, rather than the radial position, so unequal X/Y
+        // radii still produce an orthonormal tangent basis.
+        const horizontal = Math.hypot(normal[0], normal[1]);
+        this.east =
+            horizontal > 0
+                ? new Vector3(-normal[1] / horizontal, normal[0] / horizontal, 0)
+                : new Vector3(0, 1, 0);
+        this.north = new Vector3(normal).cross(this.east).normalize();
+        this._plane = new Plane().fromPointNormal(this.anchor, normal);
     }
     /**
      * Computes the projection of the provided 3D position onto the 2D plane, along the plane normal.
      *
-     * @param cartesian The WGS84 Cartesian point to project.
+     * @param cartesian The fixed Cartesian point to project.
      * @param result The optional object in which to store the local east-north coordinates.
      * @returns The supplied result or a new Vector2.
      */
-    projectPointToNearestOnPlane(cartesian, result) {
-        if (!result)
-            result = new Vector2();
-        const intersectionPoint = this._plane.projectPointOntoPlane(cartesian, scratchProjectPointOntoPlaneCartesian3);
-        const v = intersectionPoint.subtract(this._origin);
-        const x = this._xAxis.dot(v);
-        const y = this._yAxis.dot(v);
-        result.x = x;
-        result.y = y;
-        return result;
+    projectPointToNearestOnPlane(cartesian, result = new Vector2()) {
+        finiteInput(cartesian, 3);
+        // Tangential dot products equal orthogonal projection without constructing the 3D footpoint.
+        const x = cartesian[0] - this.anchor[0];
+        const y = cartesian[1] - this.anchor[1];
+        const z = cartesian[2] - this.anchor[2];
+        const projected = writePlanePoint(this, x, y, z, result, this.radius);
+        if (!projected)
+            throw new RangeError('Tangent-plane projection overflow');
+        return projected;
+    }
+    /** Project along the outward ray from the ellipsoid center.
+     * Back-facing, horizon and nonfinite intersections return undefined without writing result.
+     */
+    projectPointRadiallyOntoPlane(cartesian, result = new Vector2()) {
+        finiteInput(cartesian, 3);
+        const scale = Math.max(Math.abs(cartesian[0]), Math.abs(cartesian[1]), Math.abs(cartesian[2]));
+        if (scale === 0)
+            return undefined;
+        const sx = cartesian[0] / scale, sy = cartesian[1] / scale, sz = cartesian[2] / scale;
+        const length = Math.hypot(sx, sy, sz);
+        const x = sx / length;
+        const y = sy / length;
+        const z = sz / length;
+        const n = this.zAxis;
+        const denominator = x * n[0] + y * n[1] + z * n[2];
+        if (denominator <= 1e-14)
+            return undefined;
+        const distance = this.anchor.dot(n) / denominator;
+        return writePlanePoint(this, x * distance - this.anchor[0], y * distance - this.anchor[1], z * distance - this.anchor[2], result, this.radius);
+    }
+    /** Lift a local east/north point to the ellipsoid. Radial is the inverse of radial
+     * projection. Normal selects the nearest-facing line intersection and can miss.
+     * Neither operation restores a source point's height. Failure leaves result untouched.
+     */
+    projectPointOntoEllipsoid(point, mode = 'radial', result = new Vector3()) {
+        finiteInput(point, 2);
+        checkMode(mode);
+        const [rx, ry, rz] = this.radii;
+        if (mode === 'radial') {
+            const scale = Math.max(Math.abs(this.anchor[0]), Math.abs(this.anchor[1]), Math.abs(this.anchor[2]), Math.abs(point[0]), Math.abs(point[1]));
+            const u = point[0] / scale, v = point[1] / scale;
+            const x = this.anchor[0] / scale + this.east[0] * u + this.north[0] * v;
+            const y = this.anchor[1] / scale + this.east[1] * u + this.north[1] * v;
+            const z = this.anchor[2] / scale + this.east[2] * u + this.north[2] * v;
+            const minimum = Math.min(rx, ry, rz);
+            const ax = x * (minimum / rx), ay = y * (minimum / ry), az = z * (minimum / rz);
+            const length = Math.hypot(ax, ay, az);
+            if (!(length > 0 && Number.isFinite(length)))
+                return undefined;
+            return writeSurfacePoint(rx * (ax / length), ry * (ay / length), rz * (az / length), result);
+        }
+        const n = this.zAxis;
+        const p = [
+            this.anchor[0] + this.east[0] * point[0] + this.north[0] * point[1],
+            this.anchor[1] + this.east[1] * point[0] + this.north[1] * point[1],
+            this.anchor[2] + this.east[2] * point[0] + this.north[2] * point[1]
+        ];
+        if (!p.every(Number.isFinite))
+            return undefined;
+        const ax = n[0] / rx, ay = n[1] / ry, az = n[2] / rz;
+        const px = p[0] / rx, py = p[1] / ry, pz = p[2] / rz;
+        const a = ax * ax + ay * ay + az * az;
+        const b = ax * px + ay * py + az * pz;
+        const c = px * px + py * py + pz * pz - 1;
+        const squaredB = b * b, ac = a * c;
+        const error = 1e-6 * (Math.abs(squaredB) + Math.abs(ac));
+        const discriminant = squaredB - ac;
+        if (!(a > 0) || !Number.isFinite(discriminant) || Math.abs(discriminant) <= error) {
+            const exact = exactNormalIntersection(this.radii, this.anchor, this.east, this.north, n, point, this.radius);
+            return exact && writeSurfacePoint(exact[0], exact[1], exact[2], result);
+        }
+        if (discriminant < 0)
+            return undefined;
+        // Stable form of the nearer root (-b + sqrt(discriminant)) / a.
+        const root = Math.sqrt(Math.max(0, discriminant));
+        const t = b > 0 ? -c / (b + root) : (-b + root) / a;
+        return writeSurfacePoint(p[0] + t * n[0], p[1] + t * n[1], p[2] + t * n[2], result);
+    }
+    /** Packed XYZ to packed east/north. Invalid radial points produce NaN pairs and a zero mask.
+     * Optional output must have exactly two components per point and must not overlap input.
+     */
+    projectPointsFlat(input, mode = 'normal', output) {
+        checkMode(mode);
+        if (input.length % 3)
+            throw new RangeError('Expected packed XYZ coordinates');
+        const count = input.length / 3;
+        const positions = output || new Float64Array(count * 2);
+        if (positions.length !== count * 2)
+            throw new RangeError('Invalid tangent-plane output length');
+        if (ArrayBuffer.isView(input) &&
+            input.buffer === positions.buffer &&
+            input.byteOffset < positions.byteOffset + positions.byteLength &&
+            positions.byteOffset < input.byteOffset + input.byteLength) {
+            throw new RangeError('Input and output buffers must not overlap');
+        }
+        // Validate before writing any output, including malformed later records.
+        for (let i = 0; i < input.length; i++) {
+            if (!Number.isFinite(input[i]))
+                throw new RangeError('Expected finite coordinates');
+        }
+        const valid = new Uint8Array(count);
+        const point = new Vector3();
+        const projected = new Vector2();
+        for (let i = 0; i < count; i++) {
+            point.set(input[i * 3], input[i * 3 + 1], input[i * 3 + 2]);
+            const value = mode === 'normal'
+                ? this.projectPointToNearestOnPlane(point, projected)
+                : this.projectPointRadiallyOntoPlane(point, projected);
+            valid[i] = value ? 1 : 0;
+            positions[i * 2] = value ? value[0] : NaN;
+            positions[i * 2 + 1] = value ? value[1] : NaN;
+        }
+        return { positions, valid };
     }
     get plane() {
         return this._plane;
     }
     get origin() {
-        return this._origin;
+        return this.anchor;
     }
     get xAxis() {
-        return this._xAxis;
+        return this.east;
     }
     get yAxis() {
-        return this._yAxis;
+        return this.north;
     }
     get zAxis() {
         return this._plane.normal;
     }
 }
+function finiteInput(point, length) {
+    if (point.length !== length)
+        throw new RangeError(`Expected ${length} finite coordinates`);
+    for (let i = 0; i < length; i++) {
+        if (!Number.isFinite(point[i]))
+            throw new RangeError(`Expected ${length} finite coordinates`);
+    }
+}
+function checkMode(mode) {
+    if (mode !== 'normal' && mode !== 'radial')
+        throw new RangeError('Invalid tangent-plane mode');
+}
+function writePlanePoint(plane, x, y, z, result, radius) {
+    let east = x * plane.xAxis[0] + y * plane.xAxis[1] + z * plane.xAxis[2];
+    let north = x * plane.yAxis[0] + y * plane.yAxis[1] + z * plane.yAxis[2];
+    if (radius !== undefined &&
+        Math.abs((east / radius) ** 2 + (north / radius) ** 2 - 1) <= 32 * Number.EPSILON) {
+        east = compensatedDot3(x, y, z, plane.xAxis[0], plane.xAxis[1], plane.xAxis[2]);
+        north = compensatedDot3(x, y, z, plane.yAxis[0], plane.yAxis[1], plane.yAxis[2]);
+    }
+    if (!Number.isFinite(east) || !Number.isFinite(north))
+        return undefined;
+    return result.set(east, north);
+}
+function writeSurfacePoint(x, y, z, result) {
+    if (![x, y, z].every(Number.isFinite))
+        return undefined;
+    return result.set(x, y, z);
+}
 //# sourceMappingURL=ellipsoid-tangent-plane.js.map

package/dist/ellipsoid-tangent-plane.js.map

@@ -1 +1 @@
-{"version":3,"file":"ellipsoid-tangent-plane.js","sourceRoot":"","sources":["../src/ellipsoid-tangent-plane.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oEAAoE;AACpE,4DAA4D;AAC5D,mHAAmH;AAEnH,sEAAsE;AACtE,4FAA4F;AAE5F,OAAO,EAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AACxD,OAAO,EAAC,KAAK,EAAC,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAC,SAAS,EAAC,uBAAoB;AAEtC,MAAM,aAAa,GAAG,IAAI,OAAO,EAAE,CAAC;AACpC,MAAM,YAAY,GAAG,IAAI,OAAO,EAAE,CAAC;AACnC,MAAM,kBAAkB,GAAG,IAAI,OAAO,EAAE,CAAC;AACzC,MAAM,sCAAsC,GAAG,IAAI,OAAO,EAAE,CAAC;AAE7D,kEAAkE;AAClE,MAAM,OAAO,qBAAqB;IAMhC;;;;;;OAMG;IACH,YAAY,MAAgB,EAAE,YAAuB,SAAS,CAAC,KAAK;QAClE,MAAM,eAAe,GAAG,SAAS,CAAC,sBAAsB,CAAC,MAAM,EAAE,aAAa,CAAC,CAAC;QAEhF,MAAM,WAAW,GAAG,SAAS,CAAC,uBAAuB,CAAC,eAAe,EAAE,kBAAkB,CAAC,CAAC;QAE3F,IAAI,CAAC,OAAO,GAAG,IAAI,OAAO,CAAC,eAAe,CAAC,CAAC;QAC5C,IAAI,CAAC,MAAM,GAAG,IAAI,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvE,IAAI,CAAC,MAAM,GAAG,IAAI,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvE,MAAM,MAAM,GAAG,IAAI,OAAO,CAAC,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAExE,IAAI,CAAC,MAAM,GAAG,IAAI,KAAK,EAAE,CAAC,eAAe,CAAC,IAAI,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;IAClE,CAAC;IAED;;;;;;OAMG;IACH,4BAA4B,CAAC,SAAkB,EAAE,MAAgB;QAC/D,IAAI,CAAC,MAAM;YAAE,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QAEpC,MAAM,iBAAiB,GAAG,IAAI,CAAC,MAAM,CAAC,qBAAqB,CACzD,SAAS,EACT,sCAAsC,CACvC,CAAC;QAEF,MAAM,CAAC,GAAG,iBAAiB,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAE7B,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC;IAC5B,CAAC;CACF"}
+{"version":3,"file":"ellipsoid-tangent-plane.js","sourceRoot":"","sources":["../src/ellipsoid-tangent-plane.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oEAAoE;AACpE,4DAA4D;AAC5D,qPAAqP;AAErP,sEAAsE;AACtE,sEAAsE;AAEtE,OAAO,EAAC,OAAO,EAAE,OAAO,EAAC,MAAM,eAAe,CAAC;AAE/C,OAAO,EAAC,KAAK,EAAC,MAAM,2BAA2B,CAAC;AAChD,OAAO,EAAC,SAAS,EAAC,uBAAoB;AACtC,OAAO,EAAC,eAAe,EAAE,uBAAuB,EAAC,sDAAmD;AASpG,kEAAkE;AAClE,MAAM,OAAO,qBAAqB;IAQhC;;;;;;OAMG;IACH,YAAY,MAA8B,EAAE,YAAuB,SAAS,CAAC,KAAK;QAChF,WAAW,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;QACxD,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;YACpE,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC9B,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC;YACvE,MAAM,IAAI,UAAU,CAAC,wDAAwD,CAAC,CAAC;QACjF,CAAC;QACD,MAAM,eAAe,GAAG,SAAS,CAAC,sBAAsB,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,IAAI,OAAO,EAAE,CAAC,CAAC;QAC5F,IAAI,CAAC,eAAe;YAAE,MAAM,IAAI,UAAU,CAAC,8CAA8C,CAAC,CAAC;QAE3F,IAAI,CAAC,MAAM,GAAG,IAAI,OAAO,CAAC,eAAe,CAAC,CAAC;QAC3C,MAAM,MAAM,GAAG,IAAI,OAAO,CAAC,SAAS,CAAC,qBAAqB,CAAC,eAAe,CAAC,CAAC,CAAC;QAC7E,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;YAC1D,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;QAC9D,CAAC;QACD,MAAM,CAAC,SAAS,EAAE,CAAC;QACnB,0EAA0E;QAC1E,oDAAoD;QACpD,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACpD,IAAI,CAAC,IAAI;YACP,UAAU,GAAG,CAAC;gBACZ,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC,CAAC,GAAG,UAAU,EAAE,CAAC,CAAC;gBACjE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3B,IAAI,CAAC,KAAK,GAAG,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,SAAS,EAAE,CAAC;QAC9D,IAAI,CAAC,MAAM,GAAG,IAAI,KAAK,EAAE,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACjE,CAAC;IAED;;;;;;OAMG;IACH,4BAA4B,CAAC,SAAiC,EAAE,MAAM,GAAG,IAAI,OAAO,EAAE;QACpF,WAAW,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QAC1B,6FAA6F;QAC7F,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,SAAS,GAAG,eAAe,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACtE,IAAI,CAAC,SAAS;YAAE,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;QAC1E,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;OAEG;IACH,6BAA6B,CAC3B,SAAiC,EACjC,MAAM,GAAG,IAAI,OAAO,EAAE;QAEtB,WAAW,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QAC1B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/F,IAAI,KAAK,KAAK,CAAC;YAAE,OAAO,SAAS,CAAC;QAClC,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,KAAK,EAC7B,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,KAAK,EACzB,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAC5B,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,MAAM,CAAC,GAAG,EAAE,GAAG,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,EAAE,GAAG,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,EAAE,GAAG,MAAM,CAAC;QACtB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,MAAM,WAAW,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnD,IAAI,WAAW,IAAI,KAAK;YAAE,OAAO,SAAS,CAAC;QAC3C,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC;QAClD,OAAO,eAAe,CACpB,IAAI,EACJ,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAC7B,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAC7B,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAC7B,MAAM,EACN,IAAI,CAAC,MAAM,CACZ,CAAC;IACJ,CAAC;IAED;;;OAGG;IACH,yBAAyB,CACvB,KAA6B,EAC7B,OAAmC,QAAQ,EAC3C,MAAM,GAAG,IAAI,OAAO,EAAE;QAEtB,WAAW,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QACtB,SAAS,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QAChC,IAAI,IAAI,KAAK,QAAQ,EAAE,CAAC;YACtB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CACpB,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,EACxB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAClB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CACnB,CAAC;YACF,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,EACxB,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACvB,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACxE,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACxE,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACxE,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,EAAE,CAAC,EAC3B,EAAE,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,EAAE,CAAC,EACvB,EAAE,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC;YAC1B,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtC,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;gBAAE,OAAO,SAAS,CAAC;YAC/D,OAAO,iBAAiB,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,MAAM,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,MAAM,CAAC,EAAE,MAAM,CAAC,CAAC;QAC/F,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrB,MAAM,CAAC,GAAG;YACR,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YACnE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YACnE,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;SACpE,CAAC;QACF,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAAE,OAAO,SAAS,CAAC;QAChD,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAClB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EACd,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACjB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAClB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EACd,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACjB,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACtC,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC1C,MAAM,QAAQ,GAAG,CAAC,GAAG,CAAC,EACpB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACb,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACzD,MAAM,YAAY,GAAG,QAAQ,GAAG,EAAE,CAAC;QACnC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,IAAI,KAAK,EAAE,CAAC;YAClF,MAAM,KAAK,GAAG,uBAAuB,CACnC,IAAI,CAAC,KAAK,EACV,IAAI,CAAC,MAAM,EACX,IAAI,CAAC,IAAI,EACT,IAAI,CAAC,KAAK,EACV,CAAC,EACD,KAAK,EACL,IAAI,CAAC,MAAM,CACZ,CAAC;YACF,OAAO,KAAK,IAAI,iBAAiB,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;QAC1E,CAAC;QACD,IAAI,YAAY,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QACvC,gEAAgE;QAChE,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC;QAClD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpD,OAAO,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACtF,CAAC;IAED;;OAEG;IACH,iBAAiB,CACf,KAA6B,EAC7B,OAAmC,QAAQ,EAC3C,MAAqB;QAErB,SAAS,CAAC,IAAI,CAAC,CAAC;QAChB,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;QAC9E,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/B,MAAM,SAAS,GAAG,MAAM,IAAI,IAAI,YAAY,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QACxD,IAAI,SAAS,CAAC,MAAM,KAAK,KAAK,GAAG,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;QAChG,IACE,WAAW,CAAC,MAAM,CAAC,KAAK,CAAC;YACzB,KAAK,CAAC,MAAM,KAAK,SAAS,CAAC,MAAM;YACjC,KAAK,CAAC,UAAU,GAAG,SAAS,CAAC,UAAU,GAAG,SAAS,CAAC,UAAU;YAC9D,SAAS,CAAC,UAAU,GAAG,KAAK,CAAC,UAAU,GAAG,KAAK,CAAC,UAAU,EAC1D,CAAC;YACD,MAAM,IAAI,UAAU,CAAC,2CAA2C,CAAC,CAAC;QACpE,CAAC;QACD,yEAAyE;QACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACtC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAAE,MAAM,IAAI,UAAU,CAAC,6BAA6B,CAAC,CAAC;QACtF,CAAC;QACD,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC;QACpC,MAAM,KAAK,GAAG,IAAI,OAAO,EAAE,CAAC;QAC5B,MAAM,SAAS,GAAG,IAAI,OAAO,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC5D,MAAM,KAAK,GACT,IAAI,KAAK,QAAQ;gBACf,CAAC,CAAC,IAAI,CAAC,4BAA4B,CAAC,KAAK,EAAE,SAAS,CAAC;gBACrD,CAAC,CAAC,IAAI,CAAC,6BAA6B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC3D,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzB,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;YAC1C,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAChD,CAAC;QACD,OAAO,EAAC,SAAS,EAAE,KAAK,EAAC,CAAC;IAC5B,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,IAAI,CAAC;IACnB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC;IAC5B,CAAC;CACF;AAED,SAAS,WAAW,CAAC,KAA6B,EAAE,MAAc;IAChE,IAAI,KAAK,CAAC,MAAM,KAAK,MAAM;QAAE,MAAM,IAAI,UAAU,CAAC,YAAY,MAAM,qBAAqB,CAAC,CAAC;IAC3F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,YAAY,MAAM,qBAAqB,CAAC,CAAC;IAChG,CAAC;AACH,CAAC;AAED,SAAS,SAAS,CAAC,IAAgC;IACjD,IAAI,IAAI,KAAK,QAAQ,IAAI,IAAI,KAAK,QAAQ;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;AACjG,CAAC;AAED,SAAS,eAAe,CACtB,KAA4B,EAC5B,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAe,EACf,MAAe;IAEf,IAAI,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,IAAI,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACzE,IACE,MAAM,KAAK,SAAS;QACpB,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,MAAM,CAAC,OAAO,EACjF,CAAC;QACD,IAAI,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAChF,KAAK,GAAG,eAAe,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;QAAE,OAAO,SAAS,CAAC;IACxE,OAAO,MAAM,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;AACjC,CAAC;AAED,SAAS,iBAAiB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,MAAe;IACzE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,SAAS,CAAC;IACxD,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAC7B,CAAC"}

package/dist/ellipsoid-transforms.d.ts

@@ -0,0 +1,12 @@
+import { Matrix4 } from '@math.gl/core';
+import type { NumericArray } from '@math.gl/core/types';
+import type { Ellipsoid } from "./ellipsoid.js";
+import type { AxisDirection } from "./ellipsoid-helpers/ellipsoid-transform.js";
+export type { AxisDirection } from "./ellipsoid-helpers/ellipsoid-transform.js";
+/** Compute an east-north-up frame at a Cartesian origin on the given ellipsoid. */
+export declare function eastNorthUpToFixedFrame(ellipsoid: Ellipsoid, origin: Readonly<NumericArray>, result?: Matrix4): Matrix4;
+export declare function eastNorthUpToFixedFrame<T extends number[]>(ellipsoid: Ellipsoid, origin: Readonly<NumericArray>, result: T): T;
+/** Compute a frame with the specified axes at a Cartesian origin. */
+export declare function localFrameToFixedFrame(ellipsoid: Ellipsoid, firstAxis: AxisDirection, secondAxis: AxisDirection, thirdAxis: AxisDirection, origin: Readonly<NumericArray>, result?: Matrix4): Matrix4;
+export declare function localFrameToFixedFrame<T extends number[]>(ellipsoid: Ellipsoid, firstAxis: AxisDirection, secondAxis: AxisDirection, thirdAxis: AxisDirection, origin: Readonly<NumericArray>, result: T): T;
+//# sourceMappingURL=ellipsoid-transforms.d.ts.map

package/dist/ellipsoid-transforms.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-transforms.d.ts","sourceRoot":"","sources":["../src/ellipsoid-transforms.ts"],"names":[],"mappings":"AAMA,OAAO,EAAC,OAAO,EAAC,MAAM,eAAe,CAAC;AACtC,OAAO,KAAK,EAAC,YAAY,EAAC,MAAM,qBAAqB,CAAC;AACtD,OAAO,KAAK,EAAC,SAAS,EAAC,uBAAoB;AAE3C,OAAO,KAAK,EAAC,aAAa,EAAC,mDAAgD;AAE3E,YAAY,EAAC,aAAa,EAAC,mDAAgD;AAE3E,mFAAmF;AACnF,wBAAgB,uBAAuB,CACrC,SAAS,EAAE,SAAS,EACpB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,CAAC,EAAE,OAAO,GACf,OAAO,CAAC;AACX,wBAAgB,uBAAuB,CAAC,CAAC,SAAS,MAAM,EAAE,EACxD,SAAS,EAAE,SAAS,EACpB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,EAAE,CAAC,GACR,CAAC,CAAC;AASL,qEAAqE;AACrE,wBAAgB,sBAAsB,CACpC,SAAS,EAAE,SAAS,EACpB,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,CAAC,EAAE,OAAO,GACf,OAAO,CAAC;AACX,wBAAgB,sBAAsB,CAAC,CAAC,SAAS,MAAM,EAAE,EACvD,SAAS,EAAE,SAAS,EACpB,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,EAAE,CAAC,GACR,CAAC,CAAC"}

package/dist/ellipsoid-transforms.js

@@ -0,0 +1,15 @@
+// math.gl
+// SPDX-License-Identifier: MIT AND Apache-2.0
+// SPDX-FileCopyrightText: Copyright 2011-2018 CesiumJS Contributors
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Extracted Cesium-derived frame adapters; see ../LICENSE for attribution.
+import { Matrix4 } from '@math.gl/core';
+import { localFrameToFixedFrame as computeLocalFrame } from "./ellipsoid-helpers/ellipsoid-transform.js";
+export function eastNorthUpToFixedFrame(ellipsoid, origin, result = new Matrix4()) {
+    return computeLocalFrame(ellipsoid, 'east', 'north', 'up', origin, result);
+}
+// eslint-disable-next-line max-params
+export function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {
+    return computeLocalFrame(ellipsoid, firstAxis, secondAxis, thirdAxis, origin, result);
+}
+//# sourceMappingURL=ellipsoid-transforms.js.map

package/dist/ellipsoid-transforms.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-transforms.js","sourceRoot":"","sources":["../src/ellipsoid-transforms.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oEAAoE;AACpE,4DAA4D;AAC5D,6FAA6F;AAE7F,OAAO,EAAC,OAAO,EAAC,MAAM,eAAe,CAAC;AAGtC,OAAO,EAAC,sBAAsB,IAAI,iBAAiB,EAAC,mDAAgD;AAgBpG,MAAM,UAAU,uBAAuB,CACrC,SAAoB,EACpB,MAA8B,EAC9B,MAAM,GAAG,IAAI,OAAO,EAAE;IAEtB,OAAO,iBAAiB,CAAC,SAAS,EAAE,MAAM,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAC7E,CAAC;AAmBD,sCAAsC;AACtC,MAAM,UAAU,sBAAsB,CACpC,SAAoB,EACpB,SAAwB,EACxB,UAAyB,EACzB,SAAwB,EACxB,MAA8B,EAC9B,MAAM,GAAG,IAAI,OAAO,EAAE;IAEtB,OAAO,iBAAiB,CAAC,SAAS,EAAE,SAAS,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AACxF,CAAC"}

package/dist/ellipsoid.d.ts

@@ -1,6 +1,5 @@
-import { Vector3, Matrix4, NumericArray } from '@math.gl/core';
+import { Vector3, NumericArray } from '@math.gl/core';
 import type { SpheroidParameters } from '@math.gl/core/types';
-import type { AxisDirection } from "./ellipsoid-helpers/ellipsoid-transform.js";
 /**
  * A quadratic surface defined in Cartesian coordinates by the equation
  * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`.  Primarily used
@@ -33,7 +32,7 @@ export declare class Ellipsoid {
     toSpheroid(): SpheroidParameters;
     /** Compares this Ellipsoid against the provided Ellipsoid componentwise */
     equals(right: Ellipsoid): boolean;
-    /** Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'. */
+    /** Creates a string representing this Ellipsoid in the format '[radii.x, radii.y, radii.z]'. */
     toString(): string;
     /** Converts the provided cartographic position to Cartesian representation.
      * The position is supplied as [longitude, latitude, height], where longitude and latitude
@@ -47,15 +46,6 @@ export declare class Ellipsoid {
      */
     cartesianToCartographic(cartesian: Readonly<NumericArray>, result: Vector3): Vector3;
     cartesianToCartographic(cartesian: Readonly<NumericArray>, result?: number[]): number[];
-    /** Computes a 4x4 transformation matrix from a reference frame with an east-north-up axes
-     * centered at the provided origin to the provided ellipsoid's fixed reference frame. */
-    eastNorthUpToFixedFrame(origin: Readonly<NumericArray>, result?: Matrix4): Matrix4;
-    eastNorthUpToFixedFrame(origin: Readonly<NumericArray>, result: number[]): number[];
-    /** Computes a 4x4 transformation matrix from a reference frame centered at
-     * the provided origin to the ellipsoid's fixed reference frame.
-     */
-    localFrameToFixedFrame(firstAxis: AxisDirection, secondAxis: AxisDirection, thirdAxis: AxisDirection, origin: Readonly<NumericArray>, result?: Matrix4): Matrix4;
-    localFrameToFixedFrame<_Matrix4T>(firstAxis: AxisDirection, secondAxis: AxisDirection, thirdAxis: AxisDirection, origin: Readonly<NumericArray>, result: number[]): number[];
     /** Computes the unit vector directed from the center of this ellipsoid toward
      * the provided Cartesian position. */
     geocentricSurfaceNormal(cartesian: number[], result?: number[]): number[];

package/dist/ellipsoid.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ellipsoid.d.ts","sourceRoot":"","sources":["../src/ellipsoid.ts"],"names":[],"mappings":"AAUA,OAAO,EAAC,OAAO,EAAE,OAAO,EAA8B,YAAY,EAAC,MAAM,eAAe,CAAC;AACzF,OAAO,KAAK,EAAC,kBAAkB,EAAC,MAAM,qBAAqB,CAAC;AAW5D,OAAO,KAAK,EAAC,aAAa,EAAC,mDAAgD;AAY3E;;;;GAIG;AACH,qBAAa,SAAS;IACpB,+DAA+D;IAC/D,MAAM,CAAC,QAAQ,CAAC,KAAK,EAAE,SAAS,CAAiE;IAEjG;;OAEG;IACH,MAAM,CAAC,YAAY,CAAC,UAAU,EAAE,kBAAkB,GAAG,SAAS;IAM9D,QAAQ,CAAC,KAAK,EAAE,OAAO,CAAC;IACxB,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC;IAC/B,QAAQ,CAAC,gBAAgB,EAAE,OAAO,CAAC;IACnC,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC;IAC/B,QAAQ,CAAC,mBAAmB,EAAE,OAAO,CAAC;IACtC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,sBAAsB,EAAE,MAAM,CAAuB;IAC9D,QAAQ,CAAC,oBAAoB,EAAE,MAAM,CAAC;IACtC,OAAO,CAAC,QAAQ,CAAC,gBAAgB,CAAC,CAAmB;IAErD,4FAA4F;gBAChF,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;;IAkD3C;;;OAGG;IACH,UAAU,IAAI,kBAAkB;IAOhC,2EAA2E;IAC3E,MAAM,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAIjC,gGAAgG;IAChG,QAAQ,IAAI,MAAM;IAIlB;;;OAGG;IACH,uBAAuB,CAAC,YAAY,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IACzE,uBAAuB,CAAC,YAAY,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IA4B5E;;;OAGG;IACH,uBAAuB,CAAC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IACpF,uBAAuB,CAAC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAiDvF;4FACwF;IACxF,uBAAuB,CAAC,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAClF,uBAAuB,CAAC,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAMnF;;OAEG;IACH,sBAAsB,CACpB,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,CAAC,EAAE,OAAO,GACf,OAAO;IACV,sBAAsB,CAAC,SAAS,EAC9B,SAAS,EAAE,aAAa,EACxB,UAAU,EAAE,aAAa,EACzB,SAAS,EAAE,aAAa,EACxB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,EAAE,MAAM,EAAE,GACf,MAAM,EAAE;IAcX;0CACsC;IACtC,uBAAuB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IACzE,uBAAuB,CAAC,QAAQ,EAAE,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,QAAQ,GAAG,QAAQ;IAKlF;;OAEG;IACH,iCAAiC,CAAC,QAAQ,EACxC,YAAY,EAAE,QAAQ,CAAC,YAAY,CAAC,EACpC,MAAM,EAAE,QAAQ,GACf,QAAQ;IACX,iCAAiC,CAAC,YAAY,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAgBnE,yGAAyG;IACzG,qBAAqB,CAAC,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,SAAS,GAAG,SAAS;IACnF,qBAAqB,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAKpD;;0EAEsE;IACtE,sBAAsB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAIxE;yDACqD;IACrD,wBAAwB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IAmBrF;kEAC8D;IAC9D,8BAA8B,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IAI1F;4DACwD;IACxD,gCAAgC,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IAI5F,wFAAwF;IACxF,qCAAqC,CACnC,QAAQ,EAAE,MAAM,EAAE,EAClB,MAAM,GAAE,MAAU,EAClB,MAAM,GAAE,MAAM,EAAc,GAC3B,MAAM,EAAE;CAcZ"}
+{"version":3,"file":"ellipsoid.d.ts","sourceRoot":"","sources":["../src/ellipsoid.ts"],"names":[],"mappings":"AAUA,OAAO,EAAC,OAAO,EAAsB,YAAY,EAAS,MAAM,eAAe,CAAC;AAChF,OAAO,KAAK,EAAC,kBAAkB,EAAC,MAAM,qBAAqB,CAAC;AAsB5D;;;;GAIG;AACH,qBAAa,SAAS;IACpB,+DAA+D;IAC/D,MAAM,CAAC,QAAQ,CAAC,KAAK,EAAE,SAAS,CAAiE;IAEjG;;OAEG;IACH,MAAM,CAAC,YAAY,CAAC,UAAU,EAAE,kBAAkB,GAAG,SAAS;IAM9D,QAAQ,CAAC,KAAK,EAAE,OAAO,CAAC;IACxB,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC;IAC/B,QAAQ,CAAC,gBAAgB,EAAE,OAAO,CAAC;IACnC,QAAQ,CAAC,YAAY,EAAE,OAAO,CAAC;IAC/B,QAAQ,CAAC,mBAAmB,EAAE,OAAO,CAAC;IACtC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,sBAAsB,EAAE,MAAM,CAAuB;IAC9D,QAAQ,CAAC,oBAAoB,EAAE,MAAM,CAAC;IACtC,OAAO,CAAC,QAAQ,CAAC,gBAAgB,CAAC,CAAmB;IAErD,4FAA4F;gBAChF,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;;IAyD3C;;;OAGG;IACH,UAAU,IAAI,kBAAkB;IAOhC,2EAA2E;IAC3E,MAAM,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAIjC,gGAAgG;IAChG,QAAQ,IAAI,MAAM;IAIlB;;;OAGG;IACH,uBAAuB,CAAC,YAAY,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IACzE,uBAAuB,CAAC,YAAY,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IA4C5E;;;OAGG;IACH,uBAAuB,CAAC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IACpF,uBAAuB,CAAC,SAAS,EAAE,QAAQ,CAAC,YAAY,CAAC,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAsEvF;0CACsC;IACtC,uBAAuB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IACzE,uBAAuB,CAAC,QAAQ,EAAE,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,QAAQ,GAAG,QAAQ;IAOlF;;OAEG;IACH,iCAAiC,CAAC,QAAQ,EACxC,YAAY,EAAE,QAAQ,CAAC,YAAY,CAAC,EACpC,MAAM,EAAE,QAAQ,GACf,QAAQ;IACX,iCAAiC,CAAC,YAAY,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAuBnE,yGAAyG;IACzG,qBAAqB,CAAC,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,SAAS,GAAG,SAAS;IACnF,qBAAqB,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IA6CpD;;0EAEsE;IACtE,sBAAsB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAIxE;yDACqD;IACrD,wBAAwB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IA4BrF;kEAC8D;IAC9D,8BAA8B,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IAM1F;4DACwD;IACxD,gCAAgC,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,GAAE,MAAM,EAAc,GAAG,MAAM,EAAE;IAM5F,wFAAwF;IACxF,qCAAqC,CACnC,QAAQ,EAAE,MAAM,EAAE,EAClB,MAAM,GAAE,MAAU,EAClB,MAAM,GAAE,MAAM,EAAc,GAC3B,MAAM,EAAE;CAqBZ"}