@math.gl/geospatial 3.5.3 → 3.6.0-alpha.1

package/dist/constants.d.ts

@@ -0,0 +1,12 @@
+export declare const WGS84_RADIUS_X = 6378137;
+export declare const WGS84_RADIUS_Y = 6378137;
+export declare const WGS84_RADIUS_Z = 6356752.314245179;
+export declare const WGS84_CONSTANTS: {
+    radii: number[];
+    radiiSquared: number[];
+    oneOverRadii: number[];
+    oneOverRadiiSquared: number[];
+    maximumRadius: number;
+    centerToleranceSquared: number;
+};
+//# sourceMappingURL=constants.d.ts.map

package/dist/constants.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"constants.d.ts","sourceRoot":"","sources":["../src/constants.ts"],"names":[],"mappings":"AAGA,eAAO,MAAM,cAAc,UAAY,CAAC;AACxC,eAAO,MAAM,cAAc,UAAY,CAAC;AACxC,eAAO,MAAM,cAAc,oBAAqB,CAAC;AAIjD,eAAO,MAAM,eAAe;;;;;;;CAe3B,CAAC"}

package/dist/constants.js

@@ -0,0 +1,22 @@
+// This file is derived from the Cesium math library under Apache 2 license
+// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
+export const WGS84_RADIUS_X = 6378137.0;
+export const WGS84_RADIUS_Y = 6378137.0;
+export const WGS84_RADIUS_Z = 6356752.3142451793;
+// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`
+export const WGS84_CONSTANTS = {
+    radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],
+    radiiSquared: [
+        WGS84_RADIUS_X * WGS84_RADIUS_X,
+        WGS84_RADIUS_Y * WGS84_RADIUS_Y,
+        WGS84_RADIUS_Z * WGS84_RADIUS_Z
+    ],
+    oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],
+    oneOverRadiiSquared: [
+        1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),
+        1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),
+        1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)
+    ],
+    maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),
+    centerToleranceSquared: 1e-1 // EPSILON1;
+};

package/dist/ellipsoid/ellipsoid.d.ts

@@ -0,0 +1,68 @@
+import { Vector3, Matrix4, NumericArray } from '@math.gl/core';
+/**
+ * A quadratic surface defined in Cartesian coordinates by the equation
+ * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`.  Primarily used
+ * to represent the shape of planetary bodies.
+ */
+export default class Ellipsoid {
+    /** An Ellipsoid instance initialized to the WGS84 standard. */
+    static readonly WGS84: Ellipsoid;
+    readonly radii: Vector3;
+    readonly radiiSquared: Vector3;
+    readonly radiiToTheFourth: Vector3;
+    readonly oneOverRadii: Vector3;
+    readonly oneOverRadiiSquared: Vector3;
+    readonly minimumRadius: number;
+    readonly maximumRadius: number;
+    readonly centerToleranceSquared: number;
+    readonly squaredXOverSquaredZ: number;
+    /** Creates an Ellipsoid from a Cartesian specifying the radii in x, y, and z directions. */
+    constructor(x: number, y: number, z: number);
+    constructor();
+    /** 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)'. */
+    toString(): string;
+    /** Converts the provided cartographic to Cartesian representation. */
+    cartographicToCartesian(cartographic: number[], result: Vector3): Vector3;
+    cartographicToCartesian(cartographic: number[], result?: number[]): number[];
+    /** Converts the provided cartesian to cartographic (lng/lat/z) representation.
+     * The cartesian is undefined at the center of the ellipsoid. */
+    cartesianToCartographic(cartesian: number[], result: Vector3): Vector3;
+    cartesianToCartographic(cartesian: number[], 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: number[], result?: Matrix4): Matrix4;
+    eastNorthUpToFixedFrame(origin: number[], 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: string, secondAxis: string, thirdAxis: string, origin: Readonly<NumericArray>, result?: Matrix4): Matrix4;
+    localFrameToFixedFrame<Matrix4T>(firstAxis: string, secondAxis: string, thirdAxis: string, 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[];
+    geocentricSurfaceNormal<NumArray>(cartesian: number[], result: NumArray): NumArray;
+    /** Computes the normal of the plane tangent to the surface of the ellipsoid at provided position. */
+    geodeticSurfaceNormalCartographic<NumArray>(cartographic: number[], result: NumArray): NumArray;
+    geodeticSurfaceNormalCartographic(cartographic: number[]): number[];
+    /** Computes the normal of the plane tangent to the surface of the ellipsoid at the provided position. */
+    geodeticSurfaceNormal<NumArrayT>(cartesian: number[], result: NumArrayT): NumArrayT;
+    geodeticSurfaceNormal(cartesian: number[]): number[];
+    /** Scales the provided Cartesian position along the geodetic surface normal
+     * so that it is on the surface of this ellipsoid.  If the position is
+     * at the center of the ellipsoid, this function returns undefined. */
+    scaleToGeodeticSurface(cartesian: number[], result?: number[]): number[];
+    /** Scales the provided Cartesian position along the geocentric surface normal
+     * so that it is on the surface of this ellipsoid. */
+    scaleToGeocentricSurface(cartesian: number[], result?: number[]): number[];
+    /** Transforms a Cartesian X, Y, Z position to the ellipsoid-scaled space by multiplying
+     * its components by the result of `Ellipsoid#oneOverRadii` */
+    transformPositionToScaledSpace(position: number[], result?: number[]): number[];
+    /** Transforms a Cartesian X, Y, Z position from the ellipsoid-scaled space by multiplying
+     * its components by the result of `Ellipsoid#radii`. */
+    transformPositionFromScaledSpace(position: number[], result?: number[]): number[];
+    /** Computes a point which is the intersection of the surface normal with the z-axis. */
+    getSurfaceNormalIntersectionWithZAxis(position: number[], buffer?: number, result?: number[]): number[];
+}
+//# sourceMappingURL=ellipsoid.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid.d.ts","sourceRoot":"","sources":["../../src/ellipsoid/ellipsoid.ts"],"names":[],"mappings":"AAIA,OAAO,EACL,OAAO,EACP,OAAO,EAMP,YAAY,EACb,MAAM,eAAe,CAAC;AAkBvB;;;;GAIG;AACH,MAAM,CAAC,OAAO,OAAO,SAAS;IAC5B,+DAA+D;IAC/D,MAAM,CAAC,QAAQ,CAAC,KAAK,YAAiE;IAEtF,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;IAEtC,4FAA4F;gBAChF,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;;IAqC3C,2EAA2E;IAC3E,MAAM,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAIjC,gGAAgG;IAChG,QAAQ,IAAI,MAAM;IAIlB,sEAAsE;IACtE,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;IAoB5E;oEACgE;IAChE,uBAAuB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,OAAO,GAAG,OAAO;IACtE,uBAAuB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAsBzE;4FACwF;IACxF,uBAAuB,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IACpE,uBAAuB,CAAC,MAAM,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAMrE;;OAEG;IACH,sBAAsB,CACpB,SAAS,EAAE,MAAM,EACjB,UAAU,EAAE,MAAM,EAClB,SAAS,EAAE,MAAM,EACjB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,CAAC,EAAE,OAAO,GACf,OAAO;IACV,sBAAsB,CAAC,QAAQ,EAC7B,SAAS,EAAE,MAAM,EACjB,UAAU,EAAE,MAAM,EAClB,SAAS,EAAE,MAAM,EACjB,MAAM,EAAE,QAAQ,CAAC,YAAY,CAAC,EAC9B,MAAM,EAAE,MAAM,EAAE,GACf,MAAM,EAAE;IAQX;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,qGAAqG;IACrG,iCAAiC,CAAC,QAAQ,EAAE,YAAY,EAAE,MAAM,EAAE,EAAE,MAAM,EAAE,QAAQ,GAAG,QAAQ;IAC/F,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;IAKxE;yDACqD;IACrD,wBAAwB,CAAC,SAAS,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAoB1E;kEAC8D;IAC9D,8BAA8B,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAK/E;4DACwD;IACxD,gCAAgC,CAAC,QAAQ,EAAE,MAAM,EAAE,EAAE,MAAM,CAAC,EAAE,MAAM,EAAE,GAAG,MAAM,EAAE;IAKjF,wFAAwF;IACxF,qCAAqC,CACnC,QAAQ,EAAE,MAAM,EAAE,EAClB,MAAM,CAAC,EAAE,MAAM,EACf,MAAM,CAAC,EAAE,MAAM,EAAE,GAChB,MAAM,EAAE;CAeZ"}

package/dist/ellipsoid/ellipsoid.js

@@ -0,0 +1,132 @@
+// This file is derived from the Cesium math library under Apache 2 license
+// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
+/* eslint-disable */
+import { Vector3, Matrix4, assert, equals, _MathUtils } from '@math.gl/core';
+import * as vec3 from 'gl-matrix/vec3';
+import { WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z } from '../constants';
+import { fromCartographicToRadians, toCartographicFromRadians } from '../type-utils';
+import scaleToGeodeticSurface from './helpers/scale-to-geodetic-surface';
+import localFrameToFixedFrame from './helpers/ellipsoid-transform';
+const scratchVector = new Vector3();
+const scratchNormal = new Vector3();
+const scratchK = new Vector3();
+const scratchPosition = new Vector3();
+const scratchHeight = new Vector3();
+const scratchCartesian = new Vector3();
+let wgs84;
+/**
+ * A quadratic surface defined in Cartesian coordinates by the equation
+ * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`.  Primarily used
+ * to represent the shape of planetary bodies.
+ */
+export default class Ellipsoid {
+    constructor(x = 0.0, y = 0.0, z = 0.0) {
+        this.centerToleranceSquared = _MathUtils.EPSILON1;
+        assert(x >= 0.0);
+        assert(y >= 0.0);
+        assert(z >= 0.0);
+        this.radii = new Vector3(x, y, z);
+        this.radiiSquared = new Vector3(x * x, y * y, z * z);
+        this.radiiToTheFourth = new Vector3(x * x * x * x, y * y * y * y, z * z * z * z);
+        this.oneOverRadii = new Vector3(x === 0.0 ? 0.0 : 1.0 / x, y === 0.0 ? 0.0 : 1.0 / y, z === 0.0 ? 0.0 : 1.0 / z);
+        this.oneOverRadiiSquared = new Vector3(x === 0.0 ? 0.0 : 1.0 / (x * x), y === 0.0 ? 0.0 : 1.0 / (y * y), z === 0.0 ? 0.0 : 1.0 / (z * z));
+        this.minimumRadius = Math.min(x, y, z);
+        this.maximumRadius = Math.max(x, y, z);
+        if (this.radiiSquared.z !== 0) {
+            this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;
+        }
+        Object.freeze(this);
+    }
+    /** Compares this Ellipsoid against the provided Ellipsoid componentwise */
+    equals(right) {
+        return this === right || Boolean(right && this.radii.equals(right.radii));
+    }
+    /** Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'. */
+    toString() {
+        return this.radii.toString();
+    }
+    cartographicToCartesian(cartographic, result = [0, 0, 0]) {
+        const normal = scratchNormal;
+        const k = scratchK;
+        const [, , height] = cartographic;
+        this.geodeticSurfaceNormalCartographic(cartographic, normal);
+        k.copy(this.radiiSquared).scale(normal);
+        const gamma = Math.sqrt(normal.dot(k));
+        k.scale(1 / gamma);
+        normal.scale(height);
+        k.add(normal);
+        return k.to(result);
+    }
+    cartesianToCartographic(cartesian, result = [0, 0, 0]) {
+        scratchCartesian.from(cartesian);
+        const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);
+        if (!point) {
+            return undefined;
+        }
+        const normal = this.geodeticSurfaceNormal(point, scratchNormal);
+        const h = scratchHeight;
+        h.copy(scratchCartesian).subtract(point);
+        const longitude = Math.atan2(normal.y, normal.x);
+        const latitude = Math.asin(normal.z);
+        const height = Math.sign(vec3.dot(h, scratchCartesian)) * vec3.length(h);
+        return toCartographicFromRadians([longitude, latitude, height], result);
+    }
+    eastNorthUpToFixedFrame(origin, result = new Matrix4()) {
+        return localFrameToFixedFrame(this, 'east', 'north', 'up', origin, result);
+    }
+    // Computes a 4x4 transformation matrix from a reference frame centered at
+    // the provided origin to the ellipsoid's fixed reference frame.
+    localFrameToFixedFrame(firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {
+        return localFrameToFixedFrame(this, firstAxis, secondAxis, thirdAxis, origin, result);
+    }
+    geocentricSurfaceNormal(cartesian, result = [0, 0, 0]) {
+        return scratchVector.from(cartesian).normalize().to(result);
+    }
+    geodeticSurfaceNormalCartographic(cartographic, result = [0, 0, 0]) {
+        const cartographicVectorRadians = fromCartographicToRadians(cartographic);
+        const longitude = cartographicVectorRadians[0];
+        const latitude = cartographicVectorRadians[1];
+        const cosLatitude = Math.cos(latitude);
+        scratchVector
+            .set(cosLatitude * Math.cos(longitude), cosLatitude * Math.sin(longitude), Math.sin(latitude))
+            .normalize();
+        return scratchVector.to(result);
+    }
+    geodeticSurfaceNormal(cartesian, result = [0, 0, 0]) {
+        return scratchVector.from(cartesian).scale(this.oneOverRadiiSquared).normalize().to(result);
+    }
+    scaleToGeodeticSurface(cartesian, result) {
+        return scaleToGeodeticSurface(cartesian, this, result);
+    }
+    scaleToGeocentricSurface(cartesian, result = [0, 0, 0]) {
+        scratchPosition.from(cartesian);
+        const positionX = scratchPosition.x;
+        const positionY = scratchPosition.y;
+        const positionZ = scratchPosition.z;
+        const oneOverRadiiSquared = this.oneOverRadiiSquared;
+        const beta = 1.0 /
+            Math.sqrt(positionX * positionX * oneOverRadiiSquared.x +
+                positionY * positionY * oneOverRadiiSquared.y +
+                positionZ * positionZ * oneOverRadiiSquared.z);
+        return scratchPosition.multiplyScalar(beta).to(result);
+    }
+    transformPositionToScaledSpace(position, result = [0, 0, 0]) {
+        return scratchPosition.from(position).scale(this.oneOverRadii).to(result);
+    }
+    transformPositionFromScaledSpace(position, result = [0, 0, 0]) {
+        return scratchPosition.from(position).scale(this.radii).to(result);
+    }
+    getSurfaceNormalIntersectionWithZAxis(position, buffer = 0.0, result = [0, 0, 0]) {
+        // Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)
+        assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));
+        assert(this.radii.z > 0);
+        scratchPosition.from(position);
+        const z = scratchPosition.z * (1 - this.squaredXOverSquaredZ);
+        if (Math.abs(z) >= this.radii.z - buffer) {
+            return undefined;
+        }
+        return scratchPosition.set(0.0, 0.0, z).to(result);
+    }
+}
+/** An Ellipsoid instance initialized to the WGS84 standard. */
+Ellipsoid.WGS84 = new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);

package/dist/ellipsoid/helpers/ellipsoid-transform.d.ts

@@ -0,0 +1,2 @@
+export default function localFrameToFixedFrame(ellipsoid: any, firstAxis: any, secondAxis: any, thirdAxis: any, cartesianOrigin: any, result: any): any;
+//# sourceMappingURL=ellipsoid-transform.d.ts.map

package/dist/ellipsoid/helpers/ellipsoid-transform.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-transform.d.ts","sourceRoot":"","sources":["../../../src/ellipsoid/helpers/ellipsoid-transform.ts"],"names":[],"mappings":"AAuEA,MAAM,CAAC,OAAO,UAAU,sBAAsB,CAC5C,SAAS,KAAA,EACT,SAAS,KAAA,EACT,UAAU,KAAA,EACV,SAAS,KAAA,EACT,eAAe,KAAA,EACf,MAAM,KAAA,OAwEP"}

package/dist/ellipsoid/helpers/ellipsoid-transform.js

@@ -0,0 +1,124 @@
+import { Vector3, assert, equals as equalsEpsilon } from '@math.gl/core';
+const EPSILON14 = 1e-14;
+const scratchOrigin = new Vector3();
+// Caclulate third axis from given two axii
+const VECTOR_PRODUCT_LOCAL_FRAME = {
+    up: {
+        south: 'east',
+        north: 'west',
+        west: 'south',
+        east: 'north'
+    },
+    down: {
+        south: 'west',
+        north: 'east',
+        west: 'north',
+        east: 'south'
+    },
+    south: {
+        up: 'west',
+        down: 'east',
+        west: 'down',
+        east: 'up'
+    },
+    north: {
+        up: 'east',
+        down: 'west',
+        west: 'up',
+        east: 'down'
+    },
+    west: {
+        up: 'north',
+        down: 'south',
+        north: 'down',
+        south: 'up'
+    },
+    east: {
+        up: 'south',
+        down: 'north',
+        north: 'up',
+        south: 'down'
+    }
+};
+const degeneratePositionLocalFrame = {
+    north: [-1, 0, 0],
+    east: [0, 1, 0],
+    up: [0, 0, 1],
+    south: [1, 0, 0],
+    west: [0, -1, 0],
+    down: [0, 0, -1]
+};
+const scratchAxisVectors = {
+    east: new Vector3(),
+    north: new Vector3(),
+    up: new Vector3(),
+    west: new Vector3(),
+    south: new Vector3(),
+    down: new Vector3()
+};
+const scratchVector1 = new Vector3();
+const scratchVector2 = new Vector3();
+const scratchVector3 = new Vector3();
+// Computes a 4x4 transformation matrix from a reference frame
+// centered at the provided origin to the provided ellipsoid's fixed reference frame.
+// eslint-disable-next-line max-statements, max-params, complexity
+export default function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAxis, cartesianOrigin, result) {
+    const thirdAxisInferred = VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] && VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis];
+    // firstAxis and secondAxis must be east, north, up, west, south or down.');
+    assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));
+    let firstAxisVector;
+    let secondAxisVector;
+    let thirdAxisVector;
+    const origin = scratchOrigin.copy(cartesianOrigin);
+    // If x and y are zero, assume origin is at a pole, which is a special case.
+    const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);
+    if (atPole) {
+        // Look up axis value and adjust
+        const sign = Math.sign(origin.z);
+        firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);
+        if (firstAxis !== 'east' && firstAxis !== 'west') {
+            firstAxisVector.scale(sign);
+        }
+        secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);
+        if (secondAxis !== 'east' && secondAxis !== 'west') {
+            secondAxisVector.scale(sign);
+        }
+        thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);
+        if (thirdAxis !== 'east' && thirdAxis !== 'west') {
+            thirdAxisVector.scale(sign);
+        }
+    }
+    else {
+        // Calculate all axis
+        const { up, east, north } = scratchAxisVectors;
+        east.set(-origin.y, origin.x, 0.0).normalize();
+        ellipsoid.geodeticSurfaceNormal(origin, up);
+        north.copy(up).cross(east);
+        const { down, west, south } = scratchAxisVectors;
+        down.copy(up).scale(-1);
+        west.copy(east).scale(-1);
+        south.copy(north).scale(-1);
+        // Pick three axis based on desired orientation
+        firstAxisVector = scratchAxisVectors[firstAxis];
+        secondAxisVector = scratchAxisVectors[secondAxis];
+        thirdAxisVector = scratchAxisVectors[thirdAxis];
+    }
+    // TODO - assuming the result is column-major
+    result[0] = firstAxisVector.x;
+    result[1] = firstAxisVector.y;
+    result[2] = firstAxisVector.z;
+    result[3] = 0.0;
+    result[4] = secondAxisVector.x;
+    result[5] = secondAxisVector.y;
+    result[6] = secondAxisVector.z;
+    result[7] = 0.0;
+    result[8] = thirdAxisVector.x;
+    result[9] = thirdAxisVector.y;
+    result[10] = thirdAxisVector.z;
+    result[11] = 0.0;
+    result[12] = origin.x;
+    result[13] = origin.y;
+    result[14] = origin.z;
+    result[15] = 1.0;
+    return result;
+}

package/dist/ellipsoid/helpers/scale-to-geodetic-surface.d.ts

@@ -0,0 +1,3 @@
+import { Vector3 } from '@math.gl/core';
+export default function scaleToGeodeticSurface(cartesian: any, ellipsoid: any, result?: Vector3): any;
+//# sourceMappingURL=scale-to-geodetic-surface.d.ts.map

package/dist/ellipsoid/helpers/scale-to-geodetic-surface.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"scale-to-geodetic-surface.d.ts","sourceRoot":"","sources":["../../../src/ellipsoid/helpers/scale-to-geodetic-surface.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,OAAO,EAAqB,MAAM,eAAe,CAAC;AAU1D,MAAM,CAAC,OAAO,UAAU,sBAAsB,CAAC,SAAS,KAAA,EAAE,SAAS,KAAA,EAAE,MAAM,UAAgB,OAuF1F"}

package/dist/ellipsoid/helpers/scale-to-geodetic-surface.js

@@ -0,0 +1,70 @@
+/* eslint-disable */
+import { Vector3, _MathUtils } from '@math.gl/core';
+const scratchVector = new Vector3();
+const scaleToGeodeticSurfaceIntersection = new Vector3();
+const scaleToGeodeticSurfaceGradient = new Vector3();
+// Scales the provided Cartesian position along the geodetic surface normal
+// so that it is on the surface of this ellipsoid.  If the position is
+// at the center of the ellipsoid, this function returns undefined.
+export default function scaleToGeodeticSurface(cartesian, ellipsoid, result = new Vector3()) {
+    const { oneOverRadii, oneOverRadiiSquared, centerToleranceSquared } = ellipsoid;
+    scratchVector.from(cartesian);
+    const positionX = cartesian.x;
+    const positionY = cartesian.y;
+    const positionZ = cartesian.z;
+    const oneOverRadiiX = oneOverRadii.x;
+    const oneOverRadiiY = oneOverRadii.y;
+    const oneOverRadiiZ = oneOverRadii.z;
+    const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;
+    const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;
+    const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;
+    // Compute the squared ellipsoid norm.
+    const squaredNorm = x2 + y2 + z2;
+    const ratio = Math.sqrt(1.0 / squaredNorm);
+    // When very close to center or at center
+    if (!Number.isFinite(ratio)) {
+        return undefined;
+    }
+    // As an initial approximation, assume that the radial intersection is the projection point.
+    const intersection = scaleToGeodeticSurfaceIntersection;
+    intersection.copy(cartesian).scale(ratio);
+    // If the position is near the center, the iteration will not converge.
+    if (squaredNorm < centerToleranceSquared) {
+        return intersection.to(result);
+    }
+    const oneOverRadiiSquaredX = oneOverRadiiSquared.x;
+    const oneOverRadiiSquaredY = oneOverRadiiSquared.y;
+    const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;
+    // Use the gradient at the intersection point in place of the true unit normal.
+    // The difference in magnitude will be absorbed in the multiplier.
+    const gradient = scaleToGeodeticSurfaceGradient;
+    gradient.set(intersection.x * oneOverRadiiSquaredX * 2.0, intersection.y * oneOverRadiiSquaredY * 2.0, intersection.z * oneOverRadiiSquaredZ * 2.0);
+    // Compute the initial guess at the normal vector multiplier, lambda.
+    let lambda = ((1.0 - ratio) * cartesian.len()) / (0.5 * gradient.len());
+    let correction = 0.0;
+    let xMultiplier;
+    let yMultiplier;
+    let zMultiplier;
+    let func;
+    do {
+        lambda -= correction;
+        xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);
+        yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);
+        zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);
+        const xMultiplier2 = xMultiplier * xMultiplier;
+        const yMultiplier2 = yMultiplier * yMultiplier;
+        const zMultiplier2 = zMultiplier * zMultiplier;
+        const xMultiplier3 = xMultiplier2 * xMultiplier;
+        const yMultiplier3 = yMultiplier2 * yMultiplier;
+        const zMultiplier3 = zMultiplier2 * zMultiplier;
+        func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;
+        // "denominator" here refers to the use of this expression in the velocity and acceleration
+        // computations in the sections to follow.
+        const denominator = x2 * xMultiplier3 * oneOverRadiiSquaredX +
+            y2 * yMultiplier3 * oneOverRadiiSquaredY +
+            z2 * zMultiplier3 * oneOverRadiiSquaredZ;
+        const derivative = -2.0 * denominator;
+        correction = func / derivative;
+    } while (Math.abs(func) > _MathUtils.EPSILON12);
+    return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);
+}

package/dist/es5/constants.js

@@ -4,13 +4,13 @@ Object.defineProperty(exports, "__esModule", {
   value: true
 });
 exports.WGS84_CONSTANTS = exports.WGS84_RADIUS_Z = exports.WGS84_RADIUS_Y = exports.WGS84_RADIUS_X = void 0;
-var WGS84_RADIUS_X = 6378137.0;
+const WGS84_RADIUS_X = 6378137.0;
 exports.WGS84_RADIUS_X = WGS84_RADIUS_X;
-var WGS84_RADIUS_Y = 6378137.0;
+const WGS84_RADIUS_Y = 6378137.0;
 exports.WGS84_RADIUS_Y = WGS84_RADIUS_Y;
-var WGS84_RADIUS_Z = 6356752.3142451793;
+const WGS84_RADIUS_Z = 6356752.3142451793;
 exports.WGS84_RADIUS_Z = WGS84_RADIUS_Z;
-var WGS84_CONSTANTS = {
+const WGS84_CONSTANTS = {
   radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],
   radiiSquared: [WGS84_RADIUS_X * WGS84_RADIUS_X, WGS84_RADIUS_Y * WGS84_RADIUS_Y, WGS84_RADIUS_Z * WGS84_RADIUS_Z],
   oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],

package/dist/es5/constants.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../src/constants.js"],"names":["WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","WGS84_CONSTANTS","radii","radiiSquared","oneOverRadii","oneOverRadiiSquared","maximumRadius","Math","max","centerToleranceSquared"],"mappings":";;;;;;AAGO,IAAMA,cAAc,GAAG,SAAvB;;AACA,IAAMC,cAAc,GAAG,SAAvB;;AACA,IAAMC,cAAc,GAAG,kBAAvB;;AAIA,IAAMC,eAAe,GAAG;AAC7BC,EAAAA,KAAK,EAAE,CAACJ,cAAD,EAAiBC,cAAjB,EAAiCC,cAAjC,CADsB;AAE7BG,EAAAA,YAAY,EAAE,CACZL,cAAc,GAAGA,cADL,EAEZC,cAAc,GAAGA,cAFL,EAGZC,cAAc,GAAGA,cAHL,CAFe;AAO7BI,EAAAA,YAAY,EAAE,CAAC,MAAMN,cAAP,EAAuB,MAAMC,cAA7B,EAA6C,MAAMC,cAAnD,CAPe;AAQ7BK,EAAAA,mBAAmB,EAAE,CACnB,OAAOP,cAAc,GAAGA,cAAxB,CADmB,EAEnB,OAAOC,cAAc,GAAGA,cAAxB,CAFmB,EAGnB,OAAOC,cAAc,GAAGA,cAAxB,CAHmB,CARQ;AAa7BM,EAAAA,aAAa,EAAEC,IAAI,CAACC,GAAL,CAASV,cAAT,EAAyBC,cAAzB,EAAyCC,cAAzC,CAbc;AAc7BS,EAAAA,sBAAsB,EAAE;AAdK,CAAxB","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nexport const WGS84_RADIUS_X = 6378137.0;\nexport const WGS84_RADIUS_Y = 6378137.0;\nexport const WGS84_RADIUS_Z = 6356752.3142451793;\n\n// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`\n\nexport const WGS84_CONSTANTS = {\n  radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],\n  radiiSquared: [\n    WGS84_RADIUS_X * WGS84_RADIUS_X,\n    WGS84_RADIUS_Y * WGS84_RADIUS_Y,\n    WGS84_RADIUS_Z * WGS84_RADIUS_Z\n  ],\n  oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],\n  oneOverRadiiSquared: [\n    1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),\n    1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),\n    1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)\n  ],\n  maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),\n  centerToleranceSquared: 1e-1 // EPSILON1;\n};\n"],"file":"constants.js"}
+{"version":3,"sources":["../../src/constants.ts"],"names":["WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","WGS84_CONSTANTS","radii","radiiSquared","oneOverRadii","oneOverRadiiSquared","maximumRadius","Math","max","centerToleranceSquared"],"mappings":";;;;;;AAGO,MAAMA,cAAc,GAAG,SAAvB;;AACA,MAAMC,cAAc,GAAG,SAAvB;;AACA,MAAMC,cAAc,GAAG,kBAAvB;;AAIA,MAAMC,eAAe,GAAG;AAC7BC,EAAAA,KAAK,EAAE,CAACJ,cAAD,EAAiBC,cAAjB,EAAiCC,cAAjC,CADsB;AAE7BG,EAAAA,YAAY,EAAE,CACZL,cAAc,GAAGA,cADL,EAEZC,cAAc,GAAGA,cAFL,EAGZC,cAAc,GAAGA,cAHL,CAFe;AAO7BI,EAAAA,YAAY,EAAE,CAAC,MAAMN,cAAP,EAAuB,MAAMC,cAA7B,EAA6C,MAAMC,cAAnD,CAPe;AAQ7BK,EAAAA,mBAAmB,EAAE,CACnB,OAAOP,cAAc,GAAGA,cAAxB,CADmB,EAEnB,OAAOC,cAAc,GAAGA,cAAxB,CAFmB,EAGnB,OAAOC,cAAc,GAAGA,cAAxB,CAHmB,CARQ;AAa7BM,EAAAA,aAAa,EAAEC,IAAI,CAACC,GAAL,CAASV,cAAT,EAAyBC,cAAzB,EAAyCC,cAAzC,CAbc;AAc7BS,EAAAA,sBAAsB,EAAE;AAdK,CAAxB","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nexport const WGS84_RADIUS_X = 6378137.0;\nexport const WGS84_RADIUS_Y = 6378137.0;\nexport const WGS84_RADIUS_Z = 6356752.3142451793;\n\n// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`\n\nexport const WGS84_CONSTANTS = {\n  radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],\n  radiiSquared: [\n    WGS84_RADIUS_X * WGS84_RADIUS_X,\n    WGS84_RADIUS_Y * WGS84_RADIUS_Y,\n    WGS84_RADIUS_Z * WGS84_RADIUS_Z\n  ],\n  oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],\n  oneOverRadiiSquared: [\n    1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),\n    1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),\n    1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)\n  ],\n  maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),\n  centerToleranceSquared: 1e-1 // EPSILON1;\n};\n"],"file":"constants.js"}