@math.gl/geospatial 4.1.0 → 4.2.0-alpha.5

package/dist/make-obb-from-region.d.ts

@@ -0,0 +1,28 @@
+import { Matrix4 } from '@math.gl/core';
+import { OrientedBoundingBox } from '@math.gl/culling';
+import { Ellipsoid } from "./ellipsoid.js";
+export type GeodeticRegion = readonly [
+    west: number,
+    south: number,
+    east: number,
+    north: number,
+    minimumHeight: number,
+    maximumHeight: number
+];
+export type MakeOBBFromRegionOptions = {
+    /** Longitude and latitude units. Defaults to radians, as in 3D Tiles regions. */
+    units?: 'radians' | 'degrees';
+    /** Affine transform from ellipsoid-fixed coordinates into the caller's world space. */
+    transform?: Matrix4;
+};
+/**
+ * Creates a conservative oriented bounding box for a geodetic region.
+ *
+ * Longitudes and latitudes are radians by default. An east longitude smaller than west
+ * denotes the directed eastward interval crossing the antimeridian; equal longitudes denote
+ * zero width unless the original endpoints differ by a full turn.
+ * Heights use the ellipsoid's linear unit. Inputs and the supplied ellipsoid and transform are
+ * never mutated.
+ */
+export declare function makeOBBFromRegion(region: readonly number[], ellipsoid?: Ellipsoid, options?: MakeOBBFromRegionOptions): OrientedBoundingBox;
+//# sourceMappingURL=make-obb-from-region.d.ts.map

package/dist/make-obb-from-region.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"make-obb-from-region.d.ts","sourceRoot":"","sources":["../src/make-obb-from-region.ts"],"names":[],"mappings":"AAOA,OAAO,EAAU,OAAO,EAAiD,MAAM,eAAe,CAAC;AAC/F,OAAO,EAAC,mBAAmB,EAAQ,MAAM,kBAAkB,CAAC;AAI5D,OAAO,EAAC,SAAS,EAAC,uBAAoB;AAkCtC,MAAM,MAAM,cAAc,GAAG,SAAS;IACpC,IAAI,EAAE,MAAM;IACZ,KAAK,EAAE,MAAM;IACb,IAAI,EAAE,MAAM;IACZ,KAAK,EAAE,MAAM;IACb,aAAa,EAAE,MAAM;IACrB,aAAa,EAAE,MAAM;CACtB,CAAC;AAEF,MAAM,MAAM,wBAAwB,GAAG;IACrC,iFAAiF;IACjF,KAAK,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC;IAC9B,uFAAuF;IACvF,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB,CAAC;AAEF;;;;;;;;GAQG;AAEH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,SAAS,MAAM,EAAE,EACzB,SAAS,GAAE,SAA2B,EACtC,OAAO,GAAE,wBAA6B,GACrC,mBAAmB,CAmCrB"}

package/dist/make-obb-from-region.js

@@ -0,0 +1,246 @@
+// math.gl
+// SPDX-License-Identifier: MIT and Apache-2.0
+// Copyright (c) vis.gl contributors
+// 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
+import { Matrix3, Vector2, Vector3, degrees, radians, _MathUtils } from '@math.gl/core';
+import { OrientedBoundingBox, Plane } from '@math.gl/culling';
+// @ts-ignore
+// eslint-disable-next-line import/no-unresolved
+import { LngLatRectangle } from "./lng-lat-rectangle.js";
+import { Ellipsoid } from "./ellipsoid.js";
+import { EllipsoidTangentPlane } from "./ellipsoid-tangent-plane.js";
+const scratchOffset = new Vector3();
+const scratchTangentPoint = new Vector3();
+const scratchPerimeterCartographicNC = new Vector3();
+const scratchPerimeterCartographicNW = new Vector3();
+const scratchPerimeterCartographicCW = new Vector3();
+const scratchPerimeterCartographicSW = new Vector3();
+const scratchPerimeterCartographicSC = new Vector3();
+const scratchPerimeterCartesianNC = new Vector3();
+const scratchPerimeterCartesianNW = new Vector3();
+const scratchPerimeterCartesianCW = new Vector3();
+const scratchPerimeterCartesianSW = new Vector3();
+const scratchPerimeterCartesianSC = new Vector3();
+const scratchPerimeterProjectedNC = new Vector2();
+const scratchPerimeterProjectedNW = new Vector2();
+const scratchPerimeterProjectedCW = new Vector2();
+const scratchPerimeterProjectedSW = new Vector2();
+const scratchPerimeterProjectedSC = new Vector2();
+const scratchPlane = new Plane();
+const scratchPlaneOrigin = new Vector3();
+const scratchPlaneNormal = new Vector3();
+const scratchPlaneXAxis = new Vector3();
+const scratchHorizonCartesian = new Vector3();
+const scratchHorizonProjected = new Vector2();
+const scratchMaxY = new Vector3();
+const scratchZ = new Vector3();
+const VECTOR3_UNIT_X = new Vector3(1, 0, 0);
+const VECTOR3_UNIT_Z = new Vector3(0, 0, 1);
+/**
+ * Creates a conservative oriented bounding box for a geodetic region.
+ *
+ * Longitudes and latitudes are radians by default. An east longitude smaller than west
+ * denotes the directed eastward interval crossing the antimeridian; equal longitudes denote
+ * zero width unless the original endpoints differ by a full turn.
+ * Heights use the ellipsoid's linear unit. Inputs and the supplied ellipsoid and transform are
+ * never mutated.
+ */
+// eslint-disable-next-line max-statements
+export function makeOBBFromRegion(region, ellipsoid = Ellipsoid.WGS84, options = {}) {
+    if (!region || region.length !== 6 || region.some(value => !Number.isFinite(value))) {
+        throw new Error('makeOBBFromRegion: region must contain exactly six finite numbers');
+    }
+    const unitScale = options.units === 'degrees' ? radians(1) : 1;
+    const [rawWest, rawSouth, rawEast, rawNorth, minimumHeight, maximumHeight] = region;
+    if (minimumHeight > maximumHeight) {
+        throw new Error('makeOBBFromRegion: minimumHeight must not exceed maximumHeight');
+    }
+    const south = rawSouth * unitScale;
+    const north = rawNorth * unitScale;
+    const west = normalizeLongitude(rawWest * unitScale);
+    let east = normalizeLongitude(rawEast * unitScale);
+    // Preserve a complete directed turn, which would otherwise collapse when both endpoints
+    // are normalized to the same longitude. Reversed boundaries retain Cesium/3D Tiles semantics.
+    const directedSpan = rawEast * unitScale - rawWest * unitScale;
+    if (directedSpan >= _MathUtils.TWO_PI - 1e-12) {
+        east = west + _MathUtils.TWO_PI;
+    }
+    const latitudeTolerance = 1e-12;
+    if (south > north ||
+        south < -_MathUtils.PI_OVER_TWO - latitudeTolerance ||
+        north > _MathUtils.PI_OVER_TWO + latitudeTolerance) {
+        throw new Error('makeOBBFromRegion: latitude must be within [-π/2, π/2] and south ≤ north');
+    }
+    // Avoid passing a tiny out-of-range value to the ellipsoid conversion at a boundary.
+    const clampedSouth = Math.max(-_MathUtils.PI_OVER_TWO, Math.min(_MathUtils.PI_OVER_TWO, south));
+    const clampedNorth = Math.max(-_MathUtils.PI_OVER_TWO, Math.min(_MathUtils.PI_OVER_TWO, north));
+    return makeOBBFromNormalizedRegion([west, clampedSouth, east, clampedNorth, minimumHeight, maximumHeight], ellipsoid, options.transform);
+}
+function normalizeLongitude(value) {
+    const normalized = ((((value + Math.PI) % _MathUtils.TWO_PI) + _MathUtils.TWO_PI) % _MathUtils.TWO_PI) - Math.PI;
+    return normalized === -Math.PI && value > 0 ? Math.PI : normalized;
+}
+function makeOBBFromNormalizedRegion(region, ellipsoid, transform) {
+    const obb = new OrientedBoundingBox();
+    const [west, south, east, north, minimumHeight, maximumHeight] = region;
+    const northDeg = degrees(north);
+    const southDeg = degrees(south);
+    let maxX;
+    let maxY;
+    let maxZ;
+    let minX;
+    let minY;
+    let minZ;
+    let plane;
+    const rectangle = new LngLatRectangle(west, south, east, north);
+    const tangentPoint = LngLatRectangle.center(rectangle, scratchTangentPoint);
+    const tangentPointCartographic = new Vector3([
+        degrees(tangentPoint.x),
+        degrees(tangentPoint.y),
+        0.0
+    ]);
+    const lonCenter = tangentPoint.x;
+    const lonCenterDeg = tangentPointCartographic.x;
+    if (rectangle.width <= _MathUtils.PI) {
+        const westDeg = degrees(west);
+        const tangentPoint = ellipsoid.cartographicToCartesian(tangentPointCartographic);
+        const ellipsoidTangentPlane = new EllipsoidTangentPlane(tangentPoint, ellipsoid);
+        const latCenter = southDeg < 0.0 && northDeg > 0.0 ? 0.0 : tangentPointCartographic.y;
+        const perimeterCartographicNC = scratchPerimeterCartographicNC.copy([
+            lonCenterDeg,
+            northDeg,
+            maximumHeight
+        ]);
+        const perimeterCartographicNW = scratchPerimeterCartographicNW.copy([
+            westDeg,
+            northDeg,
+            maximumHeight
+        ]);
+        const perimeterCartographicCW = scratchPerimeterCartographicCW.copy([
+            westDeg,
+            latCenter,
+            maximumHeight
+        ]);
+        const perimeterCartographicSW = scratchPerimeterCartographicSW.copy([
+            westDeg,
+            southDeg,
+            maximumHeight
+        ]);
+        const perimeterCartographicSC = scratchPerimeterCartographicSC.copy([
+            lonCenterDeg,
+            southDeg,
+            maximumHeight
+        ]);
+        const perimeterCartesianNC = ellipsoid.cartographicToCartesian(perimeterCartographicNC, scratchPerimeterCartesianNC);
+        let perimeterCartesianNW = ellipsoid.cartographicToCartesian(perimeterCartographicNW, scratchPerimeterCartesianNW);
+        const perimeterCartesianCW = ellipsoid.cartographicToCartesian(perimeterCartographicCW, scratchPerimeterCartesianCW);
+        let perimeterCartesianSW = ellipsoid.cartographicToCartesian(perimeterCartographicSW, scratchPerimeterCartesianSW);
+        const perimeterCartesianSC = ellipsoid.cartographicToCartesian(perimeterCartographicSC, scratchPerimeterCartesianSC);
+        const perimeterProjectedNC = ellipsoidTangentPlane.projectPointToNearestOnPlane(perimeterCartesianNC, scratchPerimeterProjectedNC);
+        const perimeterProjectedNW = ellipsoidTangentPlane.projectPointToNearestOnPlane(perimeterCartesianNW, scratchPerimeterProjectedNW);
+        const perimeterProjectedCW = ellipsoidTangentPlane.projectPointToNearestOnPlane(perimeterCartesianCW, scratchPerimeterProjectedCW);
+        const perimeterProjectedSW = ellipsoidTangentPlane.projectPointToNearestOnPlane(perimeterCartesianSW, scratchPerimeterProjectedSW);
+        const perimeterProjectedSC = ellipsoidTangentPlane.projectPointToNearestOnPlane(perimeterCartesianSC, scratchPerimeterProjectedSC);
+        minX = Math.min(perimeterProjectedNW.x, perimeterProjectedCW.x, perimeterProjectedSW.x);
+        maxX = -minX;
+        maxY = Math.max(perimeterProjectedNW.y, perimeterProjectedNC.y);
+        minY = Math.min(perimeterProjectedSW.y, perimeterProjectedSC.y);
+        perimeterCartographicNW.z = perimeterCartographicSW.z = minimumHeight;
+        perimeterCartesianNW = ellipsoid.cartographicToCartesian(perimeterCartographicNW, scratchPerimeterCartesianNW);
+        perimeterCartesianSW = ellipsoid.cartographicToCartesian(perimeterCartographicSW, scratchPerimeterCartesianSW);
+        plane = ellipsoidTangentPlane.plane;
+        minZ = Math.min(plane.getPointDistance(perimeterCartesianNW), plane.getPointDistance(perimeterCartesianSW));
+        maxZ = maximumHeight;
+        return applyTransform(fromPlaneExtents(ellipsoidTangentPlane.origin, ellipsoidTangentPlane.xAxis, ellipsoidTangentPlane.yAxis, ellipsoidTangentPlane.zAxis, minX, maxX, minY, maxY, minZ, maxZ, obb), transform);
+    }
+    const eastDeg = degrees(east);
+    const fullyAboveEquator = south > 0.0;
+    const fullyBelowEquator = north < 0.0;
+    const latitudeNearestToEquator = fullyAboveEquator
+        ? southDeg
+        : fullyBelowEquator
+            ? northDeg
+            : 0.0;
+    const planeOrigin = ellipsoid.cartographicToCartesian([lonCenterDeg, latitudeNearestToEquator, maximumHeight], scratchPlaneOrigin);
+    planeOrigin.z = 0.0;
+    const isPole = Math.abs(planeOrigin.x) < _MathUtils.EPSILON10 &&
+        Math.abs(planeOrigin.y) < _MathUtils.EPSILON10;
+    const planeNormal = !isPole ? scratchPlaneNormal.copy(planeOrigin).normalize() : VECTOR3_UNIT_X;
+    const planeYAxis = VECTOR3_UNIT_Z;
+    const planeXAxis = scratchPlaneXAxis.copy(planeNormal).cross(planeYAxis);
+    plane = scratchPlane.fromPointNormal(planeOrigin, planeNormal);
+    const horizonCartesian = ellipsoid.cartographicToCartesian([degrees(lonCenter + _MathUtils.PI_OVER_TWO), latitudeNearestToEquator, maximumHeight], scratchHorizonCartesian);
+    const projectedPoint = plane.projectPointOntoPlane(horizonCartesian, scratchHorizonProjected);
+    maxX = projectedPoint.dot(planeXAxis);
+    minX = -maxX;
+    const northHeight = fullyBelowEquator ? minimumHeight : maximumHeight;
+    const southHeight = fullyAboveEquator ? minimumHeight : maximumHeight;
+    maxY = latitudeZExtent(ellipsoid, northDeg, northHeight, true);
+    minY = latitudeZExtent(ellipsoid, southDeg, southHeight, false);
+    const farZ = ellipsoid.cartographicToCartesian([eastDeg, latitudeNearestToEquator, maximumHeight], scratchZ);
+    minZ = plane.getPointDistance(farZ);
+    maxZ = 0.0;
+    return applyTransform(fromPlaneExtents(planeOrigin, planeXAxis, planeYAxis, planeNormal, minX, maxX, minY, maxY, minZ, maxZ, obb), transform);
+}
+function applyTransform(obb, transform) {
+    if (!transform)
+        return obb;
+    // A general affine transform turns an OBB into a parallelepiped. Refit its eight
+    // corners to an axis-aligned OBB so culling methods that assume orthogonal axes remain safe.
+    const axes = [0, 1, 2].map(column => obb.halfAxes.getColumn(column, new Vector3()));
+    const center = transform.transformAsPoint(obb.center);
+    const minimum = [Infinity, Infinity, Infinity];
+    const maximum = [-Infinity, -Infinity, -Infinity];
+    for (let mask = 0; mask < 8; mask++) {
+        const corner = new Vector3(center);
+        for (let axis = 0; axis < 3; axis++) {
+            const transformedAxis = transform.transformAsVector(axes[axis]);
+            corner.add(new Vector3(transformedAxis).scale((mask & (1 << axis)) === 0 ? -1 : 1));
+        }
+        for (let component = 0; component < 3; component++) {
+            minimum[component] = Math.min(minimum[component], corner[component]);
+            maximum[component] = Math.max(maximum[component], corner[component]);
+        }
+    }
+    obb.center.copy(minimum).add(maximum).scale(0.5);
+    const halfSize = new Vector3(maximum).subtract(minimum).scale(0.5);
+    obb.halfAxes = new Matrix3().set(halfSize.x, 0, 0, 0, halfSize.y, 0, 0, 0, halfSize.z);
+    return obb;
+}
+function latitudeZExtent(ellipsoid, latitude, height, maximum) {
+    let extent = maximum ? -Infinity : Infinity;
+    // A triaxial ellipsoid's fixed-latitude extrema occur at a principal longitude.
+    for (const longitude of [0, 90, 180, 270]) {
+        const z = ellipsoid.cartographicToCartesian([longitude, latitude, height], scratchMaxY).z;
+        extent = maximum ? Math.max(extent, z) : Math.min(extent, z);
+    }
+    return extent;
+}
+/** Helper function for makeOBBFromRegion(). */
+// eslint-disable-next-line max-params
+function fromPlaneExtents(planeOrigin, planeXAxis, planeYAxis, planeZAxis, minimumX, maximumX, minimumY, maximumY, minimumZ, maximumZ, result) {
+    const center = result.center;
+    const halfAxes = result.halfAxes;
+    halfAxes.setColumn(0, planeXAxis);
+    halfAxes.setColumn(1, planeYAxis);
+    halfAxes.setColumn(2, planeZAxis);
+    const centerOffset = scratchOffset
+        .set((minimumX + maximumX) / 2.0, (minimumY + maximumY) / 2.0, (minimumZ + maximumZ) / 2.0)
+        .transformByMatrix3(halfAxes);
+    const scaleX = (maximumX - minimumX) / 2.0;
+    const scaleY = (maximumY - minimumY) / 2.0;
+    const scaleZ = (maximumZ - minimumZ) / 2.0;
+    halfAxes[0] *= scaleX;
+    halfAxes[1] *= scaleX;
+    halfAxes[2] *= scaleX;
+    halfAxes[3] *= scaleY;
+    halfAxes[4] *= scaleY;
+    halfAxes[5] *= scaleY;
+    halfAxes[6] *= scaleZ;
+    halfAxes[7] *= scaleZ;
+    halfAxes[8] *= scaleZ;
+    center.copy(planeOrigin).add(centerOffset);
+    return result;
+}
+//# sourceMappingURL=make-obb-from-region.js.map

package/dist/make-obb-from-region.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"make-obb-from-region.js","sourceRoot":"","sources":["../src/make-obb-from-region.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oCAAoC;AAEpC,2EAA2E;AAC3E,4FAA4F;AAE5F,OAAO,EAAC,OAAO,EAAW,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,UAAU,EAAC,MAAM,eAAe,CAAC;AAC/F,OAAO,EAAC,mBAAmB,EAAE,KAAK,EAAC,MAAM,kBAAkB,CAAC;AAC5D,aAAa;AACb,gDAAgD;AAChD,OAAO,EAAC,eAAe,EAAC,+BAA4B;AACpD,OAAO,EAAC,SAAS,EAAC,uBAAoB;AACtC,OAAO,EAAC,qBAAqB,EAAC,qCAAkC;AAEhE,MAAM,aAAa,GAAG,IAAI,OAAO,EAAE,CAAC;AAEpC,MAAM,mBAAmB,GAAG,IAAI,OAAO,EAAE,CAAC;AAC1C,MAAM,8BAA8B,GAAG,IAAI,OAAO,EAAE,CAAC;AACrD,MAAM,8BAA8B,GAAG,IAAI,OAAO,EAAE,CAAC;AACrD,MAAM,8BAA8B,GAAG,IAAI,OAAO,EAAE,CAAC;AACrD,MAAM,8BAA8B,GAAG,IAAI,OAAO,EAAE,CAAC;AACrD,MAAM,8BAA8B,GAAG,IAAI,OAAO,EAAE,CAAC;AACrD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAClD,MAAM,2BAA2B,GAAG,IAAI,OAAO,EAAE,CAAC;AAElD,MAAM,YAAY,GAAG,IAAI,KAAK,EAAE,CAAC;AACjC,MAAM,kBAAkB,GAAG,IAAI,OAAO,EAAE,CAAC;AACzC,MAAM,kBAAkB,GAAG,IAAI,OAAO,EAAE,CAAC;AACzC,MAAM,iBAAiB,GAAG,IAAI,OAAO,EAAE,CAAC;AACxC,MAAM,uBAAuB,GAAG,IAAI,OAAO,EAAE,CAAC;AAC9C,MAAM,uBAAuB,GAAG,IAAI,OAAO,EAAE,CAAC;AAC9C,MAAM,WAAW,GAAG,IAAI,OAAO,EAAE,CAAC;AAClC,MAAM,QAAQ,GAAG,IAAI,OAAO,EAAE,CAAC;AAE/B,MAAM,cAAc,GAAG,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAC5C,MAAM,cAAc,GAAG,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAkB5C;;;;;;;;GAQG;AACH,0CAA0C;AAC1C,MAAM,UAAU,iBAAiB,CAC/B,MAAyB,EACzB,YAAuB,SAAS,CAAC,KAAK,EACtC,UAAoC,EAAE;IAEtC,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;QACpF,MAAM,IAAI,KAAK,CAAC,mEAAmE,CAAC,CAAC;IACvF,CAAC;IACD,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/D,MAAM,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,aAAa,EAAE,aAAa,CAAC,GAAG,MAAM,CAAC;IACpF,IAAI,aAAa,GAAG,aAAa,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;IACpF,CAAC;IACD,MAAM,KAAK,GAAG,QAAQ,GAAG,SAAS,CAAC;IACnC,MAAM,KAAK,GAAG,QAAQ,GAAG,SAAS,CAAC;IACnC,MAAM,IAAI,GAAG,kBAAkB,CAAC,OAAO,GAAG,SAAS,CAAC,CAAC;IACrD,IAAI,IAAI,GAAG,kBAAkB,CAAC,OAAO,GAAG,SAAS,CAAC,CAAC;IACnD,wFAAwF;IACxF,8FAA8F;IAC9F,MAAM,YAAY,GAAG,OAAO,GAAG,SAAS,GAAG,OAAO,GAAG,SAAS,CAAC;IAC/D,IAAI,YAAY,IAAI,UAAU,CAAC,MAAM,GAAG,KAAK,EAAE,CAAC;QAC9C,IAAI,GAAG,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC;IAClC,CAAC;IACD,MAAM,iBAAiB,GAAG,KAAK,CAAC;IAChC,IACE,KAAK,GAAG,KAAK;QACb,KAAK,GAAG,CAAC,UAAU,CAAC,WAAW,GAAG,iBAAiB;QACnD,KAAK,GAAG,UAAU,CAAC,WAAW,GAAG,iBAAiB,EAClD,CAAC;QACD,MAAM,IAAI,KAAK,CAAC,0EAA0E,CAAC,CAAC;IAC9F,CAAC;IACD,qFAAqF;IACrF,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,WAAW,EAAE,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,CAAC;IAChG,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,WAAW,EAAE,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,WAAW,EAAE,KAAK,CAAC,CAAC,CAAC;IAChG,OAAO,2BAA2B,CAChC,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,YAAY,EAAE,aAAa,EAAE,aAAa,CAAC,EACtE,SAAS,EACT,OAAO,CAAC,SAAS,CAClB,CAAC;AACJ,CAAC;AAED,SAAS,kBAAkB,CAAC,KAAa;IACvC,MAAM,UAAU,GACd,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IAChG,OAAO,UAAU,KAAK,CAAC,IAAI,CAAC,EAAE,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC;AACrE,CAAC;AAED,SAAS,2BAA2B,CAClC,MAAsB,EACtB,SAAoB,EACpB,SAAmB;IAEnB,MAAM,GAAG,GAAG,IAAI,mBAAmB,EAAE,CAAC;IAEtC,MAAM,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,aAAa,EAAE,aAAa,CAAC,GAAG,MAAM,CAAC;IAExE,MAAM,QAAQ,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC;IAChC,MAAM,QAAQ,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC;IAEhC,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,KAAY,CAAC;IAEjB,MAAM,SAAS,GAAG,IAAI,eAAe,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;IAChE,MAAM,YAAY,GAAG,eAAe,CAAC,MAAM,CAAC,SAAS,EAAE,mBAAmB,CAAC,CAAC;IAC5E,MAAM,wBAAwB,GAAG,IAAI,OAAO,CAAC;QAC3C,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC;QACvB,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC;QACvB,GAAG;KACJ,CAAC,CAAC;IAEH,MAAM,SAAS,GAAG,YAAY,CAAC,CAAC,CAAC;IACjC,MAAM,YAAY,GAAG,wBAAwB,CAAC,CAAC,CAAC;IAEhD,IAAI,SAAS,CAAC,KAAK,IAAI,UAAU,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;QAE9B,MAAM,YAAY,GAAG,SAAS,CAAC,uBAAuB,CAAC,wBAAwB,CAAC,CAAC;QACjF,MAAM,qBAAqB,GAAG,IAAI,qBAAqB,CAAC,YAAY,EAAE,SAAS,CAAC,CAAC;QAEjF,MAAM,SAAS,GAAG,QAAQ,GAAG,GAAG,IAAI,QAAQ,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,wBAAwB,CAAC,CAAC,CAAC;QAEtF,MAAM,uBAAuB,GAAG,8BAA8B,CAAC,IAAI,CAAC;YAClE,YAAY;YACZ,QAAQ;YACR,aAAa;SACd,CAAC,CAAC;QACH,MAAM,uBAAuB,GAAG,8BAA8B,CAAC,IAAI,CAAC;YAClE,OAAO;YACP,QAAQ;YACR,aAAa;SACd,CAAC,CAAC;QACH,MAAM,uBAAuB,GAAG,8BAA8B,CAAC,IAAI,CAAC;YAClE,OAAO;YACP,SAAS;YACT,aAAa;SACd,CAAC,CAAC;QACH,MAAM,uBAAuB,GAAG,8BAA8B,CAAC,IAAI,CAAC;YAClE,OAAO;YACP,QAAQ;YACR,aAAa;SACd,CAAC,CAAC;QACH,MAAM,uBAAuB,GAAG,8BAA8B,CAAC,IAAI,CAAC;YAClE,YAAY;YACZ,QAAQ;YACR,aAAa;SACd,CAAC,CAAC;QAEH,MAAM,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CAC5D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,IAAI,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CAC1D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CAC5D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,IAAI,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CAC1D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CAC5D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QAEF,MAAM,oBAAoB,GAAG,qBAAqB,CAAC,4BAA4B,CAC7E,oBAAoB,EACpB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,qBAAqB,CAAC,4BAA4B,CAC7E,oBAAoB,EACpB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,qBAAqB,CAAC,4BAA4B,CAC7E,oBAAoB,EACpB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,qBAAqB,CAAC,4BAA4B,CAC7E,oBAAoB,EACpB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,qBAAqB,CAAC,4BAA4B,CAC7E,oBAAoB,EACpB,2BAA2B,CAC5B,CAAC;QAEF,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,oBAAoB,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,CAAC,CAAC;QACxF,IAAI,GAAG,CAAC,IAAI,CAAC;QAEb,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,oBAAoB,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,CAAC,CAAC;QAChE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,oBAAoB,CAAC,CAAC,EAAE,oBAAoB,CAAC,CAAC,CAAC,CAAC;QAEhE,uBAAuB,CAAC,CAAC,GAAG,uBAAuB,CAAC,CAAC,GAAG,aAAa,CAAC;QACtE,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CACtD,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,oBAAoB,GAAG,SAAS,CAAC,uBAAuB,CACtD,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QAEF,KAAK,GAAG,qBAAqB,CAAC,KAAK,CAAC;QACpC,IAAI,GAAG,IAAI,CAAC,GAAG,CACb,KAAK,CAAC,gBAAgB,CAAC,oBAAoB,CAAC,EAC5C,KAAK,CAAC,gBAAgB,CAAC,oBAAoB,CAAC,CAC7C,CAAC;QACF,IAAI,GAAG,aAAa,CAAC;QAErB,OAAO,cAAc,CACnB,gBAAgB,CACd,qBAAqB,CAAC,MAAM,EAC5B,qBAAqB,CAAC,KAAK,EAC3B,qBAAqB,CAAC,KAAK,EAC3B,qBAAqB,CAAC,KAAK,EAC3B,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,GAAG,CACJ,EACD,SAAS,CACV,CAAC;IACJ,CAAC;IAED,MAAM,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE9B,MAAM,iBAAiB,GAAG,KAAK,GAAG,GAAG,CAAC;IACtC,MAAM,iBAAiB,GAAG,KAAK,GAAG,GAAG,CAAC;IACtC,MAAM,wBAAwB,GAAG,iBAAiB;QAChD,CAAC,CAAC,QAAQ;QACV,CAAC,CAAC,iBAAiB;YACjB,CAAC,CAAC,QAAQ;YACV,CAAC,CAAC,GAAG,CAAC;IAEV,MAAM,WAAW,GAAG,SAAS,CAAC,uBAAuB,CACnD,CAAC,YAAY,EAAE,wBAAwB,EAAE,aAAa,CAAC,EACvD,kBAAkB,CACnB,CAAC;IACF,WAAW,CAAC,CAAC,GAAG,GAAG,CAAC;IACpB,MAAM,MAAM,GACV,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,SAAS;QAC9C,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,SAAS,CAAC;IACjD,MAAM,WAAW,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,kBAAkB,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,cAAc,CAAC;IAChG,MAAM,UAAU,GAAG,cAAc,CAAC;IAClC,MAAM,UAAU,GAAG,iBAAiB,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC;IACzE,KAAK,GAAG,YAAY,CAAC,eAAe,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;IAE/D,MAAM,gBAAgB,GAAG,SAAS,CAAC,uBAAuB,CACxD,CAAC,OAAO,CAAC,SAAS,GAAG,UAAU,CAAC,WAAW,CAAC,EAAE,wBAAwB,EAAE,aAAa,CAAC,EACtF,uBAAuB,CACxB,CAAC;IACF,MAAM,cAAc,GAAG,KAAK,CAAC,qBAAqB,CAChD,gBAAgB,EAChB,uBAAuB,CACb,CAAC;IACb,IAAI,GAAG,cAAc,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACtC,IAAI,GAAG,CAAC,IAAI,CAAC;IAEb,MAAM,WAAW,GAAG,iBAAiB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC;IACtE,MAAM,WAAW,GAAG,iBAAiB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC;IACtE,IAAI,GAAG,eAAe,CAAC,SAAS,EAAE,QAAQ,EAAE,WAAW,EAAE,IAAI,CAAC,CAAC;IAC/D,IAAI,GAAG,eAAe,CAAC,SAAS,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,CAAC,CAAC;IAEhE,MAAM,IAAI,GAAG,SAAS,CAAC,uBAAuB,CAC5C,CAAC,OAAO,EAAE,wBAAwB,EAAE,aAAa,CAAC,EAClD,QAAQ,CACT,CAAC;IAEF,IAAI,GAAG,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC;IACpC,IAAI,GAAG,GAAG,CAAC;IAEX,OAAO,cAAc,CACnB,gBAAgB,CACd,WAAW,EACX,UAAU,EACV,UAAU,EACV,WAAW,EACX,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,GAAG,CACJ,EACD,SAAS,CACV,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CAAC,GAAwB,EAAE,SAAmB;IACnE,IAAI,CAAC,SAAS;QAAE,OAAO,GAAG,CAAC;IAC3B,iFAAiF;IACjF,6FAA6F;IAC7F,MAAM,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,QAAQ,CAAC,SAAS,CAAC,MAAM,EAAE,IAAI,OAAO,EAAE,CAAC,CAAC,CAAC;IACpF,MAAM,MAAM,GAAG,SAAS,CAAC,gBAAgB,CAAC,GAAG,CAAC,MAAM,CAAa,CAAC;IAClE,MAAM,OAAO,GAAG,CAAC,QAAQ,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAC/C,MAAM,OAAO,GAAG,CAAC,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,CAAC;IAClD,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC;QACnC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;YACpC,MAAM,eAAe,GAAG,SAAS,CAAC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAa,CAAC;YAC5E,MAAM,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,eAAe,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACtF,CAAC;QACD,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE,CAAC;YACnD,OAAO,CAAC,SAAS,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;YACrE,OAAO,CAAC,SAAS,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;QACvE,CAAC;IACH,CAAC;IACD,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACjD,MAAM,QAAQ,GAAG,IAAI,OAAO,CAAC,OAAO,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACnE,GAAG,CAAC,QAAQ,GAAG,IAAI,OAAO,EAAE,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC;IACvF,OAAO,GAAG,CAAC;AACb,CAAC;AAED,SAAS,eAAe,CACtB,SAAoB,EACpB,QAAgB,EAChB,MAAc,EACd,OAAgB;IAEhB,IAAI,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC;IAC5C,gFAAgF;IAChF,KAAK,MAAM,SAAS,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;QAC1C,MAAM,CAAC,GAAG,SAAS,CAAC,uBAAuB,CAAC,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC;QAC1F,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,+CAA+C;AAC/C,sCAAsC;AACtC,SAAS,gBAAgB,CACvB,WAAoB,EACpB,UAAmB,EACnB,UAAmB,EACnB,UAAmB,EACnB,QAAgB,EAChB,QAAgB,EAChB,QAAgB,EAChB,QAAgB,EAChB,QAAgB,EAChB,QAAgB,EAChB,MAA2B;IAE3B,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;IAC7B,MAAM,QAAQ,GAAG,MAAM,CAAC,QAAQ,CAAC;IACjC,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IAClC,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IAClC,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IAElC,MAAM,YAAY,GAAG,aAAa;SAC/B,GAAG,CAAC,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,EAAE,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,EAAE,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;SAC1F,kBAAkB,CAAC,QAAQ,CAAC,CAAC;IAEhC,MAAM,MAAM,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC3C,MAAM,MAAM,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC3C,MAAM,MAAM,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC3C,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACtB,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IAEtB,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IAE3C,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/package.json

@@ -6,7 +6,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "4.1.0",
+  "version": "4.2.0-alpha.5",
   "keywords": [
     "webgl",
     "javascript",
@@ -39,8 +39,9 @@
     "src"
   ],
   "dependencies": {
-    "@math.gl/core": "4.1.0",
-    "@math.gl/types": "4.1.0"
+    "@math.gl/core": "4.2.0-alpha.5",
+    "@math.gl/culling": "4.2.0-alpha.5",
+    "@math.gl/types": "4.2.0-alpha.5"
   },
-  "gitHead": "76820f54e2a9034b42bd24ffeb381cc4e01ad46e"
+  "gitHead": "9c13785761867748e2edb167b5f2749dce803039"
 }

package/src/ellipsoid-tangent-plane.ts

@@ -0,0 +1,87 @@
+// math.gl
+// SPDX-License-Identifier: MIT and Apache-2.0
+// Copyright (c) vis.gl contributors
+
+// 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';
+import {Plane} from '@math.gl/culling';
+import {Ellipsoid} from './ellipsoid';
+
+const scratchOrigin = new Vector3();
+const scratchCart3 = new Vector3();
+const scratchEastNorthUp = new Matrix4();
+const scratchProjectPointOntoPlaneCartesian3 = new Vector3();
+
+/** A two-dimensional east-north plane tangent to an ellipsoid. */
+export class EllipsoidTangentPlane {
+  private _origin: Vector3;
+  private _xAxis: Vector3;
+  private _yAxis: Vector3;
+  private _plane: Plane;
+
+  /**
+   * Creates a new plane tangent to the WGS84 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: number[], ellipsoid: 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);
+  }
+
+  /**
+   * 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 result The optional object in which to store the local east-north coordinates.
+   * @returns The supplied result or a new Vector2.
+   */
+  projectPointToNearestOnPlane(cartesian: Vector3, result?: Vector2): Vector2 {
+    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;
+  }
+
+  get plane(): Plane {
+    return this._plane;
+  }
+
+  get origin(): Vector3 {
+    return this._origin;
+  }
+
+  get xAxis(): Vector3 {
+    return this._xAxis;
+  }
+
+  get yAxis(): Vector3 {
+    return this._yAxis;
+  }
+
+  get zAxis(): Vector3 {
+    return this._plane.normal;
+  }
+}

package/src/{ellipsoid/ellipsoid.ts → ellipsoid.ts}

@@ -8,12 +8,12 @@
 /* eslint-disable */
 import {Vector3, Matrix4, assert, equals, _MathUtils, NumericArray, vec3} from '@math.gl/core';
 
-import {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from '../constants';
-import {fromCartographicToRadians, toCartographicFromRadians} from '../type-utils';
+import {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from './constants';
+import {fromCartographicToRadians, toCartographicFromRadians} from './type-utils';
 
-import type {AxisDirection} from './helpers/ellipsoid-transform';
-import {localFrameToFixedFrame} from './helpers/ellipsoid-transform';
-import {scaleToGeodeticSurface} from './helpers/scale-to-geodetic-surface';
+import type {AxisDirection} from './ellipsoid-helpers/ellipsoid-transform';
+import {localFrameToFixedFrame} from './ellipsoid-helpers/ellipsoid-transform';
+import {scaleToGeodeticSurface} from './ellipsoid-helpers/scale-to-geodetic-surface';
 
 const scratchVector = new Vector3();
 const scratchNormal = new Vector3();
@@ -89,7 +89,10 @@ export class Ellipsoid {
     return this.radii.toString();
   }
 
-  /** Converts the provided cartographic to Cartesian representation. */
+  /** Converts the provided cartographic position to Cartesian representation.
+   * The position is supplied as [longitude, latitude, height], where longitude and latitude
+   * are in degrees and height is in meters above the ellipsoid.
+   */
   cartographicToCartesian(cartographic: number[], result: Vector3): Vector3;
   cartographicToCartesian(cartographic: number[], result?: number[]): number[];
 
@@ -111,8 +114,10 @@ export class Ellipsoid {
     return k.to(result);
   }
 
-  /** Converts the provided cartesian to cartographic (lng/lat/z) representation.
-   * The cartesian is undefined at the center of the ellipsoid. */
+  /** Converts the provided Cartesian position to [longitude, latitude, height],
+   * where longitude and latitude are in degrees and height is in meters above the ellipsoid.
+   * The cartographic position is undefined at the center of the ellipsoid.
+   */
   cartesianToCartographic(cartesian: Readonly<NumericArray>, result: Vector3): Vector3;
   cartesianToCartographic(cartesian: Readonly<NumericArray>, result?: number[]): number[];
 
@@ -155,7 +160,7 @@ export class Ellipsoid {
     origin: Readonly<NumericArray>,
     result?: Matrix4
   ): Matrix4;
-  localFrameToFixedFrame<Matrix4T>(
+  localFrameToFixedFrame<_Matrix4T>(
     firstAxis: AxisDirection,
     secondAxis: AxisDirection,
     thirdAxis: AxisDirection,
@@ -183,7 +188,9 @@ export class Ellipsoid {
     return scratchVector.from(cartesian).normalize().to(result);
   }
 
-  /** Computes the normal of the plane tangent to the surface of the ellipsoid at provided position. */
+  /** Computes the normal of the plane tangent to the surface of the ellipsoid at the provided
+   * [longitude, latitude, height], where longitude and latitude are in degrees.
+   */
   geodeticSurfaceNormalCartographic<NumArray>(
     cartographic: Readonly<NumericArray>,
     result: NumArray

package/src/index.ts

@@ -2,5 +2,9 @@
 // SPDX-License-Identifier: MIT
 // Copyright (c) vis.gl contributors
 
-export {Ellipsoid} from './ellipsoid/ellipsoid';
+export {Ellipsoid} from './ellipsoid';
+export {EllipsoidTangentPlane} from './ellipsoid-tangent-plane';
+export {LngLatRectangle} from './lng-lat-rectangle';
+export {makeOBBFromRegion} from './make-obb-from-region';
+export type {GeodeticRegion, MakeOBBFromRegionOptions} from './make-obb-from-region';
 export {isWGS84} from './type-utils';

package/src/lng-lat-rectangle.ts

@@ -0,0 +1,67 @@
+// math.gl
+// SPDX-License-Identifier: MIT and Apache-2.0
+// Copyright (c) vis.gl contributors
+
+// 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 {Vector3, normalizeAngle, _MathUtils} from '@math.gl/core';
+
+/** A two dimensional region specified as longitude and latitude coordinates. */
+export class LngLatRectangle {
+  west: number;
+  south: number;
+  east: number;
+  north: number;
+
+  /**
+   * Creates a new two dimensional region specified as longitude and latitude coordinates.
+   *
+   * @param west The westernmost longitude, in radians, in the range [-Pi, Pi].
+   * @param south The southernmost latitude, in radians, in the range [-Pi/2, Pi/2].
+   * @param east The easternmost longitude, in radians, in the range [-Pi, Pi].
+   * @param north The northernmost latitude, in radians, in the range [-Pi/2, Pi/2].
+   */
+  constructor(west: number, south: number, east: number, north: number) {
+    this.west = west;
+    this.south = south;
+    this.east = east;
+    this.north = north;
+  }
+
+  /**
+   * Computes the center of a rectangle.
+   *
+   * @param rectangle The rectangle for which to find the center
+   * @param  [result] The object onto which to store the result.
+   * @returns  The modified result parameter or a new Cartographic instance if none was provided.
+   */
+  static center(rectangle: LngLatRectangle, result?: Vector3): Vector3 {
+    if (!result) result = new Vector3();
+
+    let east = rectangle.east;
+    const west = rectangle.west;
+
+    if (east < west) {
+      east += _MathUtils.TWO_PI;
+    }
+
+    const longitude = normalizeAngle((west + east) * 0.5, 'negative-pi-to-pi');
+    const latitude = (rectangle.south + rectangle.north) * 0.5;
+
+    result.x = longitude;
+    result.y = latitude;
+    result.z = 0.0;
+    return result;
+  }
+
+  get width(): number {
+    let east = this.east;
+    const west = this.west;
+
+    if (east < west) {
+      east += _MathUtils.TWO_PI;
+    }
+    return east - west;
+  }
+}