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

package/dist/lng-lat-rectangle.js

@@ -0,0 +1,54 @@
+// 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 {
+    /**
+     * 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, south, east, north) {
+        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, result) {
+        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() {
+        let east = this.east;
+        const west = this.west;
+        if (east < west) {
+            east += _MathUtils.TWO_PI;
+        }
+        return east - west;
+    }
+}
+//# sourceMappingURL=lng-lat-rectangle.js.map

package/dist/lng-lat-rectangle.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lng-lat-rectangle.js","sourceRoot":"","sources":["../src/lng-lat-rectangle.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,8CAA8C;AAC9C,oCAAoC;AAEpC,sEAAsE;AACtE,4FAA4F;AAE5F,OAAO,EAAC,OAAO,EAAE,cAAc,EAAE,UAAU,EAAC,MAAM,eAAe,CAAC;AAElE,gFAAgF;AAChF,MAAM,OAAO,eAAe;IAM1B;;;;;;;OAOG;IACH,YAAY,IAAY,EAAE,KAAa,EAAE,IAAY,EAAE,KAAa;QAClE,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACrB,CAAC;IAED;;;;;;OAMG;IACH,MAAM,CAAC,MAAM,CAAC,SAA0B,EAAE,MAAgB;QACxD,IAAI,CAAC,MAAM;YAAE,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QAEpC,IAAI,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC;QAC1B,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC;QAE5B,IAAI,IAAI,GAAG,IAAI,EAAE,CAAC;YAChB,IAAI,IAAI,UAAU,CAAC,MAAM,CAAC;QAC5B,CAAC;QAED,MAAM,SAAS,GAAG,cAAc,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,mBAAmB,CAAC,CAAC;QAC3E,MAAM,QAAQ,GAAG,CAAC,SAAS,CAAC,KAAK,GAAG,SAAS,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QAE3D,MAAM,CAAC,CAAC,GAAG,SAAS,CAAC;QACrB,MAAM,CAAC,CAAC,GAAG,QAAQ,CAAC;QACpB,MAAM,CAAC,CAAC,GAAG,GAAG,CAAC;QACf,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,IAAI,KAAK;QACP,IAAI,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;QACrB,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;QAEvB,IAAI,IAAI,GAAG,IAAI,EAAE,CAAC;YAChB,IAAI,IAAI,UAAU,CAAC,MAAM,CAAC;QAC5B,CAAC;QACD,OAAO,IAAI,GAAG,IAAI,CAAC;IACrB,CAAC;CACF"}

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

@@ -0,0 +1,11 @@
+import { OrientedBoundingBox } from '@math.gl/culling';
+/**
+ * Computes an OrientedBoundingBox that bounds a region on the surface of the WGS84 ellipsoid.
+ * There are no guarantees about the orientation of the bounding box.
+ *
+ * @param {number[]} region The cartographic region ([west, south, east, north, minimum height, maximum height])
+ *                          on the surface of the ellipsoid.
+ * @returns {OrientedBoundingBox} The modified result parameter or a new OrientedBoundingBox instance if none was provided.
+ */
+export declare function makeOBBfromRegion(region: number[]): 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":"AAQA,OAAO,EAAC,mBAAmB,EAAQ,MAAM,kBAAkB,CAAC;AAuC5D;;;;;;;GAOG;AAEH,wBAAgB,iBAAiB,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,mBAAmB,CA6MvE"}

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

@@ -0,0 +1,173 @@
+// 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 { Vector2, Vector3, degrees, _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 scratchMinY = new Vector3();
+const scratchZ = new Vector3();
+const VECTOR3_UNIT_X = new Vector3(1, 0, 0);
+const VECTOR3_UNIT_Z = new Vector3(0, 0, 1);
+/**
+ * Computes an OrientedBoundingBox that bounds a region on the surface of the WGS84 ellipsoid.
+ * There are no guarantees about the orientation of the bounding box.
+ *
+ * @param {number[]} region The cartographic region ([west, south, east, north, minimum height, maximum height])
+ *                          on the surface of the ellipsoid.
+ * @returns {OrientedBoundingBox} The modified result parameter or a new OrientedBoundingBox instance if none was provided.
+ */
+// eslint-disable-next-line max-statements
+export function makeOBBfromRegion(region) {
+    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.WGS84.cartographicToCartesian(tangentPointCartographic);
+        const ellipsoidTangentPlane = new EllipsoidTangentPlane(tangentPoint);
+        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.WGS84.cartographicToCartesian(perimeterCartographicNC, scratchPerimeterCartesianNC);
+        let perimeterCartesianNW = Ellipsoid.WGS84.cartographicToCartesian(perimeterCartographicNW, scratchPerimeterCartesianNW);
+        const perimeterCartesianCW = Ellipsoid.WGS84.cartographicToCartesian(perimeterCartographicCW, scratchPerimeterCartesianCW);
+        let perimeterCartesianSW = Ellipsoid.WGS84.cartographicToCartesian(perimeterCartographicSW, scratchPerimeterCartesianSW);
+        const perimeterCartesianSC = Ellipsoid.WGS84.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.WGS84.cartographicToCartesian(perimeterCartographicNW, scratchPerimeterCartesianNW);
+        perimeterCartesianSW = Ellipsoid.WGS84.cartographicToCartesian(perimeterCartographicSW, scratchPerimeterCartesianSW);
+        plane = ellipsoidTangentPlane.plane;
+        minZ = Math.min(plane.getPointDistance(perimeterCartesianNW), plane.getPointDistance(perimeterCartesianSW));
+        maxZ = maximumHeight;
+        return fromPlaneExtents(ellipsoidTangentPlane.origin, ellipsoidTangentPlane.xAxis, ellipsoidTangentPlane.yAxis, ellipsoidTangentPlane.zAxis, minX, maxX, minY, maxY, minZ, maxZ, obb);
+    }
+    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.WGS84.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.WGS84.cartographicToCartesian([degrees(lonCenter + _MathUtils.PI_OVER_TWO), latitudeNearestToEquator, maximumHeight], scratchHorizonCartesian);
+    const projectedPoint = plane.projectPointOntoPlane(horizonCartesian, scratchHorizonProjected);
+    maxX = projectedPoint.dot(planeXAxis);
+    minX = -maxX;
+    maxY = Ellipsoid.WGS84.cartographicToCartesian([0.0, northDeg, fullyBelowEquator ? minimumHeight : maximumHeight], scratchMaxY).z;
+    minY = Ellipsoid.WGS84.cartographicToCartesian([0.0, southDeg, fullyAboveEquator ? minimumHeight : maximumHeight], scratchMinY).z;
+    const farZ = Ellipsoid.WGS84.cartographicToCartesian([eastDeg, latitudeNearestToEquator, maximumHeight], scratchZ);
+    minZ = plane.getPointDistance(farZ);
+    maxZ = 0.0;
+    return fromPlaneExtents(planeOrigin, planeXAxis, planeYAxis, planeNormal, minX, maxX, minY, maxY, minZ, maxZ, obb);
+}
+const scratchScale = new Vector3();
+/** helper function for OrientedBoundingBox.fromRegion() */
+// 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);
+    let centerOffset = scratchOffset;
+    centerOffset.x = (minimumX + maximumX) / 2.0;
+    centerOffset.y = (minimumY + maximumY) / 2.0;
+    centerOffset.z = (minimumZ + maximumZ) / 2.0;
+    centerOffset = halfAxes.multiplyByVector(centerOffset, scratchOffset);
+    const scale = scratchScale;
+    scale.x = (maximumX - minimumX) / 2.0;
+    scale.y = (maximumY - minimumY) / 2.0;
+    scale.z = (maximumZ - minimumZ) / 2.0;
+    planeOrigin = center.copy(planeOrigin);
+    result.center = planeOrigin.add(centerOffset);
+    result.halfAxes = halfAxes.multiplyByScale(scale, halfAxes);
+    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,EAAE,OAAO,EAAE,OAAO,EAAE,UAAU,EAAC,MAAM,eAAe,CAAC;AACpE,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,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;AAE5C;;;;;;;GAOG;AACH,0CAA0C;AAC1C,MAAM,UAAU,iBAAiB,CAAC,MAAgB;IAChD,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,KAAK,CAAC,uBAAuB,CAAC,wBAAwB,CAAC,CAAC;QACvF,MAAM,qBAAqB,GAAG,IAAI,qBAAqB,CAAC,YAAY,CAAC,CAAC;QAEtE,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,KAAK,CAAC,uBAAuB,CAClE,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,IAAI,oBAAoB,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAChE,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAClE,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,IAAI,oBAAoB,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAChE,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,MAAM,oBAAoB,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAClE,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,KAAK,CAAC,uBAAuB,CAC5D,uBAAuB,EACvB,2BAA2B,CAC5B,CAAC;QACF,oBAAoB,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAC5D,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,gBAAgB,CACrB,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,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,KAAK,CAAC,uBAAuB,CACzD,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,KAAK,CAAC,uBAAuB,CAC9D,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,IAAI,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAC5C,CAAC,GAAG,EAAE,QAAQ,EAAE,iBAAiB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC,EAClE,WAAW,CACZ,CAAC,CAAC,CAAC;IACJ,IAAI,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAC5C,CAAC,GAAG,EAAE,QAAQ,EAAE,iBAAiB,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,aAAa,CAAC,EAClE,WAAW,CACZ,CAAC,CAAC,CAAC;IAEJ,MAAM,IAAI,GAAG,SAAS,CAAC,KAAK,CAAC,uBAAuB,CAClD,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,gBAAgB,CACrB,WAAW,EACX,UAAU,EACV,UAAU,EACV,WAAW,EACX,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,GAAG,CACJ,CAAC;AACJ,CAAC;AAED,MAAM,YAAY,GAAG,IAAI,OAAO,EAAE,CAAC;AAEnC,2DAA2D;AAC3D,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,IAAI,YAAY,GAAG,aAAa,CAAC;IACjC,YAAY,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC7C,YAAY,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC7C,YAAY,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAC7C,YAAY,GAAG,QAAQ,CAAC,gBAAgB,CAAC,YAAY,EAAE,aAAa,CAAC,CAAC;IAEtE,MAAM,KAAK,GAAG,YAAY,CAAC;IAC3B,KAAK,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IACtC,KAAK,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IACtC,KAAK,CAAC,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC;IAEtC,WAAW,GAAG,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IACvC,MAAM,CAAC,MAAM,GAAG,WAAW,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IAC9C,MAAM,CAAC,QAAQ,GAAG,QAAQ,CAAC,eAAe,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;IAE5D,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/package.json

@@ -6,7 +6,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "4.1.0",
+  "version": "4.2.0-alpha.3",
   "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.3",
+    "@math.gl/culling": "4.2.0-alpha.3",
+    "@math.gl/types": "4.2.0-alpha.3"
   },
-  "gitHead": "76820f54e2a9034b42bd24ffeb381cc4e01ad46e"
+  "gitHead": "4238aa95fb38da43a65c6e5f2b8343a875666d92"
 }

package/src/ellipsoid-tangent-plane.ts

@@ -0,0 +1,98 @@
+// 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, Ray} from '@math.gl/culling';
+import {Ellipsoid} from './ellipsoid';
+
+const scratchOrigin = new Vector3();
+const scratchCart3 = new Vector3();
+const scratchEastNorthUp = new Matrix4();
+const scratchPlane = new Plane();
+
+const scratchProjectPointOntoPlaneRay = new Ray();
+const scratchProjectPointOntoPlaneCartesian3 = new Vector3();
+const scratchDirection = new Vector3();
+
+/** A plane tangent to the WGS84 ellipsoid at the provided origin */
+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, it's surface projection will be used.
+   *
+   * @param {Cartesian3} origin The point on the surface of the ellipsoid where the tangent plane touches.
+   */
+  constructor(origin: number[]) {
+    origin = Ellipsoid.WGS84.scaleToGeodeticSurface(origin, scratchOrigin);
+
+    const eastNorthUp = Ellipsoid.WGS84.eastNorthUpToFixedFrame(origin, scratchEastNorthUp);
+
+    this._origin = origin as Vector3;
+    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 = scratchPlane.fromPointNormal(origin, normal);
+  }
+
+  /**
+   * Computes the projection of the provided 3D position onto the 2D plane, along the plane normal.
+   *
+   * @param {Vector3} cartesian The point to project.
+   * @param {Vector2} [result] The object onto which to store the result.
+   * @returns {Vector2} The modified result parameter or a new Cartesian2 instance if none was provided.
+   */
+  projectPointToNearestOnPlane(cartesian: Vector3, result?: Vector2): Vector2 {
+    if (!result) result = new Vector2();
+
+    const plane = this._plane;
+
+    const ray = scratchProjectPointOntoPlaneRay;
+    scratchProjectPointOntoPlaneRay.origin = cartesian;
+    scratchProjectPointOntoPlaneRay.direction = scratchDirection.copy(plane.normal);
+
+    let intersectionPoint = plane.intersectWithRay(ray, scratchProjectPointOntoPlaneCartesian3);
+
+    if (!intersectionPoint) {
+      ray.direction = ray.direction.negate();
+      intersectionPoint = plane.intersectWithRay(ray, 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() {
+    return this._plane;
+  }
+
+  get origin() {
+    return this._origin;
+  }
+
+  get xAxis() {
+    return this._xAxis;
+  }
+
+  get yAxis() {
+    return this._yAxis;
+  }
+
+  get zAxis() {
+    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();

package/src/index.ts

@@ -2,5 +2,7 @@
 // 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 {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() {
+    let east = this.east;
+    const west = this.west;
+
+    if (east < west) {
+      east += _MathUtils.TWO_PI;
+    }
+    return east - west;
+  }
+}