@math.gl/geospatial 5.0.0-alpha.9 → 5.0.0-beta.1

package/dist/projection/experimental/datum.d.ts

@@ -0,0 +1,17 @@
+import type { Datum, Ellipsoid } from "./crs/types.js";
+import type { ProjectionPoint } from "./types.js";
+export type Coordinate3D = [number, number, number];
+export declare function geodeticToGeocentricInPlace(point: ProjectionPoint, ellipsoid: Ellipsoid): void;
+/** Projection retains its center error and historical near-axis canonicalization. */
+export declare function geocentricToGeodeticInPlace(point: ProjectionPoint, ellipsoid: Ellipsoid): void;
+/** Convert rotation units and scale once for the lifetime of a compiled datum stage. */
+export declare function createHelmert(values: readonly number[], inverse: boolean): DatumOperation;
+type DatumOperation = (point: ProjectionPoint) => void;
+/** Resolve identity, grid and Helmert stages once, outside coordinate loops. */
+export declare function createDatumTransform(source: Datum, destination: Datum): DatumOperation | undefined;
+/** Original per-instance dispatch informed by proj4js grid list semantics. */
+export declare function applyDatumGrids(point: ProjectionPoint, datum: Datum, inverse: boolean): void;
+export declare function geodeticToGeocentric(coordinate: Coordinate3D, ellipsoid: Ellipsoid): Coordinate3D;
+export declare function geocentricToGeodetic(coordinate: Coordinate3D, ellipsoid: Ellipsoid): Coordinate3D;
+export {};
+//# sourceMappingURL=datum.d.ts.map

package/dist/projection/experimental/datum.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"datum.d.ts","sourceRoot":"","sources":["../../../src/projection/experimental/datum.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,KAAK,EAAE,SAAS,EAAC,uBAAoB;AAClD,OAAO,KAAK,EAAC,eAAe,EAAC,mBAAgB;AAC7C,MAAM,MAAM,YAAY,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAWpD,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,eAAe,EAAE,SAAS,EAAE,SAAS,GAAG,IAAI,CAI9F;AAED,qFAAqF;AACrF,wBAAgB,2BAA2B,CAAC,KAAK,EAAE,eAAe,EAAE,SAAS,EAAE,SAAS,GAAG,IAAI,CAuB9F;AAED,wFAAwF;AACxF,wBAAgB,aAAa,CAAC,MAAM,EAAE,SAAS,MAAM,EAAE,EAAE,OAAO,EAAE,OAAO,GAAG,cAAc,CAyBzF;AAID,KAAK,cAAc,GAAG,CAAC,KAAK,EAAE,eAAe,KAAK,IAAI,CAAC;AACvD,gFAAgF;AAChF,wBAAgB,oBAAoB,CAClC,MAAM,EAAE,KAAK,EACb,WAAW,EAAE,KAAK,GACjB,cAAc,GAAG,SAAS,CAqC5B;AAuBD,8EAA8E;AAC9E,wBAAgB,eAAe,CAAC,KAAK,EAAE,eAAe,EAAE,KAAK,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,GAAG,IAAI,CA0B5F;AAED,wBAAgB,oBAAoB,CAAC,UAAU,EAAE,YAAY,EAAE,SAAS,EAAE,SAAS,GAAG,YAAY,CAIjG;AACD,wBAAgB,oBAAoB,CAAC,UAAU,EAAE,YAAY,EAAE,SAAS,EAAE,SAAS,GAAG,YAAY,CAIjG"}

package/dist/projection/experimental/datum.js

@@ -0,0 +1,160 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Helmert equations directly adapted from proj4js 2.22.0 datumUtils.js. See ../../../projection/PROJ4-LICENSE.md. Geocentric equations now delegate to the attributed @math.gl/core leaf.
+import { spheroidToCartesian, cartesianToSpheroid } from '@math.gl/core';
+const ARC_SECOND = Math.PI / (180 * 3600);
+const WGS84 = {
+    ellipsoid: {
+        semiMajorAxis: 6378137,
+        semiMinorAxis: 6356752.314245179,
+        eccentricitySquared: 0.0066943799901413165
+    },
+    towgs84: [0, 0, 0]
+};
+export function geodeticToGeocentricInPlace(point, ellipsoid) {
+    if (Math.abs(point.y) > Math.PI / 2)
+        throw new Error('Latitude outside geocentric domain');
+    if (!spheroidToCartesian(point, ellipsoid))
+        throw new Error('Geocentric coordinate must be finite');
+}
+/** Projection retains its center error and historical near-axis canonicalization. */
+export function geocentricToGeodeticInPlace(point, ellipsoid) {
+    const x = point.x, y = point.y, z = point.z;
+    const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = ellipsoid;
+    if (x === 0 && y === 0 && z === 0)
+        throw new Error('Geodetic coordinates are undefined at the Earth center');
+    if (es !== 0 &&
+        Math.abs(x) < 1e-12 * a &&
+        Math.abs(y) < 1e-12 * a &&
+        Math.hypot(x, y) < 1e-12 * a) {
+        if (!Number.isFinite(Math.hypot(x, y, z)))
+            throw new Error('Geocentric radius must be finite');
+        point.x = 0;
+        point.y = (Math.sign(z) * Math.PI) / 2;
+        point.z = Math.abs(z) - b;
+        return;
+    }
+    if (!cartesianToSpheroid(point, ellipsoid)) {
+        if (!Number.isFinite(Math.hypot(x, y, z)))
+            throw new Error('Geocentric radius must be finite');
+        throw new Error('Geocentric inverse did not converge');
+    }
+}
+/** Convert rotation units and scale once for the lifetime of a compiled datum stage. */
+export function createHelmert(values, inverse) {
+    const dx = values[0], dy = values[1], dz = values[2];
+    const rx = (values[3] || 0) * ARC_SECOND, ry = (values[4] || 0) * ARC_SECOND, rz = (values[5] || 0) * ARC_SECOND;
+    const scale = 1 + (values[6] || 0) / 1e6;
+    return inverse
+        ? point => {
+            const x = (point.x - dx) / scale, y = (point.y - dy) / scale, z = (point.z - dz) / scale;
+            point.x = x + rz * y - ry * z;
+            point.y = -rz * x + y + rx * z;
+            point.z = ry * x - rx * y + z;
+        }
+        : point => {
+            const x = point.x, y = point.y, z = point.z;
+            point.x = scale * (x - rz * y + ry * z) + dx;
+            point.y = scale * (rz * x + y - rx * z) + dy;
+            point.z = scale * (-ry * x + rx * y + z) + dz;
+        };
+}
+function shifted(datum) {
+    return Boolean(datum.towgs84?.some(value => value !== 0));
+}
+/** Resolve identity, grid and Helmert stages once, outside coordinate loops. */
+export function createDatumTransform(source, destination) {
+    // datum=none disables the entire chain, including a shifted opposite endpoint.
+    if (!source.towgs84 || !destination.towgs84)
+        return undefined;
+    if (source.grids || destination.grids) {
+        if (source.grids &&
+            destination.grids &&
+            source.ellipsoid.semiMajorAxis === destination.ellipsoid.semiMajorAxis &&
+            source.ellipsoid.eccentricitySquared === destination.ellipsoid.eccentricitySquared &&
+            source.grids.length === destination.grids.length &&
+            source.grids.every((entry, i) => entry.name === destination.grids[i].name &&
+                entry.optional === destination.grids[i].optional &&
+                entry.grid === destination.grids[i].grid))
+            return undefined;
+        const middle = createDatumTransform(source.grids ? WGS84 : source, destination.grids ? WGS84 : destination);
+        return point => {
+            if (source.grids)
+                applyDatumGrids(point, source, false);
+            middle?.(point);
+            if (destination.grids)
+                applyDatumGrids(point, destination, true);
+        };
+    }
+    if (shifted(source) || shifted(destination)) {
+        const first = convertDatum(source, WGS84), second = convertDatum(WGS84, destination);
+        return point => {
+            first?.(point);
+            second?.(point);
+        };
+    }
+    return convertDatum(source, destination);
+}
+function convertDatum(from, to) {
+    if (!from.towgs84 || !to.towgs84)
+        return undefined;
+    const a = from.ellipsoid, b = to.ellipsoid;
+    let sameParameters = true;
+    for (let i = 0; i < 7; i++)
+        if ((from.towgs84[i] || 0) !== (to.towgs84[i] || 0))
+            sameParameters = false;
+    if (a.semiMajorAxis === b.semiMajorAxis &&
+        Math.abs(a.eccentricitySquared - b.eccentricitySquared) <= 5e-11 &&
+        sameParameters)
+        return undefined;
+    const sourceShift = shifted(from) ? createHelmert(from.towgs84, false) : undefined, targetShift = shifted(to) ? createHelmert(to.towgs84, true) : undefined;
+    return point => {
+        geodeticToGeocentricInPlace(point, a);
+        sourceShift?.(point);
+        targetShift?.(point);
+        geocentricToGeodeticInPlace(point, b);
+    };
+}
+/** Original per-instance dispatch informed by proj4js grid list semantics. */
+export function applyDatumGrids(point, datum, inverse) {
+    for (const reference of datum.grids) {
+        if (reference.name === 'null')
+            return;
+        if (!reference.grid) {
+            if (reference.optional)
+                continue;
+            throw new Error('Required datum grid is not registered: ' + reference.name);
+        }
+        let found = false;
+        if (reference.grid.shiftInPlace)
+            found = reference.grid.shiftInPlace(point, inverse);
+        else {
+            const output = reference.grid.shift(point.x, point.y, inverse);
+            if (output) {
+                point.x = output[0];
+                point.y = output[1];
+                found = true;
+            }
+        }
+        if (found) {
+            if (!Number.isFinite(point.x) || !Number.isFinite(point.y) || Math.abs(point.y) > Math.PI / 2)
+                throw new Error('Datum grid produced a coordinate outside the geographic domain');
+            return;
+        }
+    }
+    throw new Error('No datum grid covers coordinate: ' + datum.grids.map(grid => grid.name).join(','));
+}
+// Scalar adapters retain the internal tuple API for existing callers.
+export function geodeticToGeocentric(coordinate, ellipsoid) {
+    const point = { x: coordinate[0], y: coordinate[1], z: coordinate[2] };
+    geodeticToGeocentricInPlace(point, ellipsoid);
+    return [point.x, point.y, point.z];
+}
+export function geocentricToGeodetic(coordinate, ellipsoid) {
+    const point = { x: coordinate[0], y: coordinate[1], z: coordinate[2] };
+    geocentricToGeodeticInPlace(point, ellipsoid);
+    return [point.x, point.y, point.z];
+}
+//# sourceMappingURL=datum.js.map

package/dist/projection/experimental/datum.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"datum.js","sourceRoot":"","sources":["../../../src/projection/experimental/datum.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,4MAA4M;AAC5M,OAAO,EAAC,mBAAmB,EAAE,mBAAmB,EAAC,MAAM,eAAe,CAAC;AAIvE,MAAM,UAAU,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;AAC1C,MAAM,KAAK,GAAU;IACnB,SAAS,EAAE;QACT,aAAa,EAAE,OAAO;QACtB,aAAa,EAAE,iBAAiB;QAChC,mBAAmB,EAAE,qBAAqB;KAC3C;IACD,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;CACnB,CAAC;AAEF,MAAM,UAAU,2BAA2B,CAAC,KAAsB,EAAE,SAAoB;IACtF,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;IAC3F,IAAI,CAAC,mBAAmB,CAAC,KAAK,EAAE,SAAS,CAAC;QACxC,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;AAC5D,CAAC;AAED,qFAAqF;AACrF,MAAM,UAAU,2BAA2B,CAAC,KAAsB,EAAE,SAAoB;IACtF,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EACf,CAAC,GAAG,KAAK,CAAC,CAAC,EACX,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACd,MAAM,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,mBAAmB,EAAE,EAAE,EAAC,GAAG,SAAS,CAAC;IAChF,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC;QAC/B,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,IACE,EAAE,KAAK,CAAC;QACR,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC;QACvB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC;QACvB,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAC5B,CAAC;QACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;QAC/F,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QACZ,KAAK,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;QACvC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,OAAO;IACT,CAAC;IACD,IAAI,CAAC,mBAAmB,CAAC,KAAK,EAAE,SAAS,CAAC,EAAE,CAAC;QAC3C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;QAC/F,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IACzD,CAAC;AACH,CAAC;AAED,wFAAwF;AACxF,MAAM,UAAU,aAAa,CAAC,MAAyB,EAAE,OAAgB;IACvE,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,EAClB,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,EACd,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IACjB,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,UAAU,EACtC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,UAAU,EAClC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,UAAU,CAAC;IACrC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;IACzC,OAAO,OAAO;QACZ,CAAC,CAAC,KAAK,CAAC,EAAE;YACN,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,EAC9B,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,EAC1B,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;YAC7B,KAAK,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YAC9B,KAAK,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YAC/B,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QAChC,CAAC;QACH,CAAC,CAAC,KAAK,CAAC,EAAE;YACN,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,EACf,CAAC,GAAG,KAAK,CAAC,CAAC,EACX,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YACd,KAAK,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAC7C,KAAK,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAC7C,KAAK,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,CAAC;AACR,CAAC;AACD,SAAS,OAAO,CAAC,KAAY;IAC3B,OAAO,OAAO,CAAC,KAAK,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC;AAC5D,CAAC;AAED,gFAAgF;AAChF,MAAM,UAAU,oBAAoB,CAClC,MAAa,EACb,WAAkB;IAElB,+EAA+E;IAC/E,IAAI,CAAC,MAAM,CAAC,OAAO,IAAI,CAAC,WAAW,CAAC,OAAO;QAAE,OAAO,SAAS,CAAC;IAC9D,IAAI,MAAM,CAAC,KAAK,IAAI,WAAW,CAAC,KAAK,EAAE,CAAC;QACtC,IACE,MAAM,CAAC,KAAK;YACZ,WAAW,CAAC,KAAK;YACjB,MAAM,CAAC,SAAS,CAAC,aAAa,KAAK,WAAW,CAAC,SAAS,CAAC,aAAa;YACtE,MAAM,CAAC,SAAS,CAAC,mBAAmB,KAAK,WAAW,CAAC,SAAS,CAAC,mBAAmB;YAClF,MAAM,CAAC,KAAK,CAAC,MAAM,KAAK,WAAW,CAAC,KAAK,CAAC,MAAM;YAChD,MAAM,CAAC,KAAK,CAAC,KAAK,CAChB,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CACX,KAAK,CAAC,IAAI,KAAK,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI;gBACxC,KAAK,CAAC,QAAQ,KAAK,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ;gBAChD,KAAK,CAAC,IAAI,KAAK,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAC3C;YAED,OAAO,SAAS,CAAC;QACnB,MAAM,MAAM,GAAG,oBAAoB,CACjC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,MAAM,EAC7B,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,WAAW,CACxC,CAAC;QACF,OAAO,KAAK,CAAC,EAAE;YACb,IAAI,MAAM,CAAC,KAAK;gBAAE,eAAe,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,CAAC,CAAC;YACxD,MAAM,EAAE,CAAC,KAAK,CAAC,CAAC;YAChB,IAAI,WAAW,CAAC,KAAK;gBAAE,eAAe,CAAC,KAAK,EAAE,WAAW,EAAE,IAAI,CAAC,CAAC;QACnE,CAAC,CAAC;IACJ,CAAC;IACD,IAAI,OAAO,CAAC,MAAM,CAAC,IAAI,OAAO,CAAC,WAAW,CAAC,EAAE,CAAC;QAC5C,MAAM,KAAK,GAAG,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,EACvC,MAAM,GAAG,YAAY,CAAC,KAAK,EAAE,WAAW,CAAC,CAAC;QAC5C,OAAO,KAAK,CAAC,EAAE;YACb,KAAK,EAAE,CAAC,KAAK,CAAC,CAAC;YACf,MAAM,EAAE,CAAC,KAAK,CAAC,CAAC;QAClB,CAAC,CAAC;IACJ,CAAC;IACD,OAAO,YAAY,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;AAC3C,CAAC;AACD,SAAS,YAAY,CAAC,IAAW,EAAE,EAAS;IAC1C,IAAI,CAAC,IAAI,CAAC,OAAO,IAAI,CAAC,EAAE,CAAC,OAAO;QAAE,OAAO,SAAS,CAAC;IACnD,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,EACtB,CAAC,GAAG,EAAE,CAAC,SAAS,CAAC;IACnB,IAAI,cAAc,GAAG,IAAI,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QACxB,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YAAE,cAAc,GAAG,KAAK,CAAC;IAC9E,IACE,CAAC,CAAC,aAAa,KAAK,CAAC,CAAC,aAAa;QACnC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,mBAAmB,GAAG,CAAC,CAAC,mBAAmB,CAAC,IAAI,KAAK;QAChE,cAAc;QAEd,OAAO,SAAS,CAAC;IACnB,MAAM,WAAW,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,IAAI,CAAC,OAAO,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,EAChF,WAAW,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,EAAE,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;IAC1E,OAAO,KAAK,CAAC,EAAE;QACb,2BAA2B,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QACtC,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC;QACrB,WAAW,EAAE,CAAC,KAAK,CAAC,CAAC;QACrB,2BAA2B,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;IACxC,CAAC,CAAC;AACJ,CAAC;AACD,8EAA8E;AAC9E,MAAM,UAAU,eAAe,CAAC,KAAsB,EAAE,KAAY,EAAE,OAAgB;IACpF,KAAK,MAAM,SAAS,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC;QACpC,IAAI,SAAS,CAAC,IAAI,KAAK,MAAM;YAAE,OAAO;QACtC,IAAI,CAAC,SAAS,CAAC,IAAI,EAAE,CAAC;YACpB,IAAI,SAAS,CAAC,QAAQ;gBAAE,SAAS;YACjC,MAAM,IAAI,KAAK,CAAC,yCAAyC,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;QAC9E,CAAC;QACD,IAAI,KAAK,GAAG,KAAK,CAAC;QAClB,IAAI,SAAS,CAAC,IAAI,CAAC,YAAY;YAAE,KAAK,GAAG,SAAS,CAAC,IAAI,CAAC,YAAY,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC;aAChF,CAAC;YACJ,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YAC/D,IAAI,MAAM,EAAE,CAAC;gBACX,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;gBACpB,KAAK,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;gBACpB,KAAK,GAAG,IAAI,CAAC;YACf,CAAC;QACH,CAAC;QACD,IAAI,KAAK,EAAE,CAAC;YACV,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC;gBAC3F,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;YACpF,OAAO;QACT,CAAC;IACH,CAAC;IACD,MAAM,IAAI,KAAK,CACb,mCAAmC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CACnF,CAAC;AACJ,CAAC;AACD,sEAAsE;AACtE,MAAM,UAAU,oBAAoB,CAAC,UAAwB,EAAE,SAAoB;IACjF,MAAM,KAAK,GAAG,EAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAC,CAAC;IACrE,2BAA2B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC9C,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AACD,MAAM,UAAU,oBAAoB,CAAC,UAAwB,EAAE,SAAoB;IACjF,MAAM,KAAK,GAAG,EAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,EAAC,CAAC;IACrE,2BAA2B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC9C,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC"}

package/dist/projection/experimental/deformation-inverse.d.ts

@@ -0,0 +1,5 @@
+import type { ProjectionPoint } from "./types.js";
+export type FixedDisplacement = (point: ProjectionPoint, source: number, target: number) => void;
+/** Source-sampled map inverse. Successful iterations reuse the caller's working point. */
+export declare function inverseDisplacement(point: ProjectionPoint, x: number, y: number, z: number, source: number, target: number, sample: FixedDisplacement): void;
+//# sourceMappingURL=deformation-inverse.d.ts.map

package/dist/projection/experimental/deformation-inverse.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"deformation-inverse.d.ts","sourceRoot":"","sources":["../../../src/projection/experimental/deformation-inverse.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,eAAe,EAAC,mBAAgB;AAC7C,MAAM,MAAM,iBAAiB,GAAG,CAAC,KAAK,EAAE,eAAe,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,KAAK,IAAI,CAAC;AACjG,0FAA0F;AAC1F,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,eAAe,EACtB,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,MAAM,EAAE,MAAM,EACd,MAAM,EAAE,MAAM,EACd,MAAM,EAAE,iBAAiB,GACxB,IAAI,CA6GN"}

package/dist/projection/experimental/deformation-inverse.js

@@ -0,0 +1,89 @@
+/** Source-sampled map inverse. Successful iterations reuse the caller's working point. */
+export function inverseDisplacement(point, x, y, z, source, target, sample) {
+    let px = x, py = y, pz = z;
+    for (let iteration = 0; iteration < 16; iteration++) {
+        point.x = px;
+        point.y = py;
+        point.z = pz;
+        sample(point, source, target);
+        const rx = px - x + point.x, ry = py - y + point.y, rz = pz - z + point.z;
+        const norm = Math.hypot(rx, ry, rz);
+        if (!Number.isFinite(norm))
+            throw new Error('Non-finite temporal inverse residual');
+        if (norm <= 1e-8) {
+            point.x = px;
+            point.y = py;
+            point.z = pz;
+            return;
+        }
+        const epsilon = Math.sqrt(Number.EPSILON) * Math.max(1, Math.abs(px), Math.abs(py), Math.abs(pz));
+        let j00 = 0, j01 = 0, j02 = 0, j10 = 0, j11 = 0, j12 = 0, j20 = 0, j21 = 0, j22 = 0;
+        for (let axis = 0; axis < 3; axis++) {
+            point.x = px + (axis === 0 ? epsilon : 0);
+            point.y = py + (axis === 1 ? epsilon : 0);
+            point.z = pz + (axis === 2 ? epsilon : 0);
+            sample(point, source, target);
+            const vx = point.x, vy = point.y, vz = point.z;
+            point.x = px - (axis === 0 ? epsilon : 0);
+            point.y = py - (axis === 1 ? epsilon : 0);
+            point.z = pz - (axis === 2 ? epsilon : 0);
+            sample(point, source, target);
+            const dx = (vx - point.x) / (2 * epsilon), dy = (vy - point.y) / (2 * epsilon), dz = (vz - point.z) / (2 * epsilon);
+            if (axis === 0) {
+                j00 = 1 + dx;
+                j10 = dy;
+                j20 = dz;
+            }
+            else if (axis === 1) {
+                j01 = dx;
+                j11 = 1 + dy;
+                j21 = dz;
+            }
+            else {
+                j02 = dx;
+                j12 = dy;
+                j22 = 1 + dz;
+            }
+        }
+        const c00 = j11 * j22 - j12 * j21, c01 = j02 * j21 - j01 * j22, c02 = j01 * j12 - j02 * j11;
+        const c10 = j12 * j20 - j10 * j22, c11 = j00 * j22 - j02 * j20, c12 = j02 * j10 - j00 * j12;
+        const c20 = j10 * j21 - j11 * j20, c21 = j01 * j20 - j00 * j21, c22 = j00 * j11 - j01 * j10;
+        const determinant = j00 * c00 + j01 * c10 + j02 * c20;
+        if (!Number.isFinite(determinant) || determinant === 0)
+            throw new Error('Singular temporal inverse Jacobian');
+        const sx = (c00 * rx + c01 * ry + c02 * rz) / determinant, sy = (c10 * rx + c11 * ry + c12 * rz) / determinant, sz = (c20 * rx + c21 * ry + c22 * rz) / determinant;
+        if (!Number.isFinite(sx) || !Number.isFinite(sy) || !Number.isFinite(sz))
+            throw new Error('Non-finite temporal inverse update');
+        let accepted = false;
+        for (let half = 0; half < 12; half++) {
+            const scale = 2 ** -half, nx = px - scale * sx, ny = py - scale * sy, nz = pz - scale * sz;
+            point.x = nx;
+            point.y = ny;
+            point.z = nz;
+            try {
+                sample(point, source, target);
+                const error = Math.hypot(nx - x + point.x, ny - y + point.y, nz - z + point.z);
+                if (Number.isFinite(error) && error < norm) {
+                    if (error <= 1e-8) {
+                        point.x = nx;
+                        point.y = ny;
+                        point.z = nz;
+                        return;
+                    }
+                    px = nx;
+                    py = ny;
+                    pz = nz;
+                    accepted = true;
+                    break;
+                }
+            }
+            catch {
+                /* Domain/coverage failures reject this trial; no extrapolation. */
+            }
+        }
+        if (!accepted)
+            throw new Error('Temporal inverse line search failed');
+    }
+    throw new Error('Temporal inverse did not converge');
+}
+//# sourceMappingURL=deformation-inverse.js.map

package/dist/projection/experimental/deformation-inverse.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"deformation-inverse.js","sourceRoot":"","sources":["../../../src/projection/experimental/deformation-inverse.ts"],"names":[],"mappings":"AAMA,0FAA0F;AAC1F,MAAM,UAAU,mBAAmB,CACjC,KAAsB,EACtB,CAAS,EACT,CAAS,EACT,CAAS,EACT,MAAc,EACd,MAAc,EACd,MAAyB;IAEzB,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,CAAC,EACN,EAAE,GAAG,CAAC,CAAC;IACT,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,EAAE,EAAE,SAAS,EAAE,EAAE,CAAC;QACpD,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;QACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;QACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;QACb,MAAM,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;QAC9B,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,EACzB,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,EACrB,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QACxB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACpC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;QACpF,IAAI,IAAI,IAAI,IAAI,EAAE,CAAC;YACjB,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,OAAO;QACT,CAAC;QACD,MAAM,OAAO,GACX,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACpF,IAAI,GAAG,GAAG,CAAC,EACT,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,EACP,GAAG,GAAG,CAAC,CAAC;QACV,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;YACpC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,MAAM,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;YAC9B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,EAChB,EAAE,GAAG,KAAK,CAAC,CAAC,EACZ,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;YACf,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,MAAM,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;YAC9B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EACvC,EAAE,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EACnC,EAAE,GAAG,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;YACtC,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;gBACf,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;gBACb,GAAG,GAAG,EAAE,CAAC;gBACT,GAAG,GAAG,EAAE,CAAC;YACX,CAAC;iBAAM,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;gBACtB,GAAG,GAAG,EAAE,CAAC;gBACT,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;gBACb,GAAG,GAAG,EAAE,CAAC;YACX,CAAC;iBAAM,CAAC;gBACN,GAAG,GAAG,EAAE,CAAC;gBACT,GAAG,GAAG,EAAE,CAAC;gBACT,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;YACf,CAAC;QACH,CAAC;QACD,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC/B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC3B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC/B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC3B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC/B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAC3B,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,WAAW,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACtD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,WAAW,KAAK,CAAC;YACpD,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;QACxD,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,WAAW,EACvD,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,WAAW,EACnD,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,WAAW,CAAC;QACtD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;YACtE,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;QACxD,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,EAAE,EAAE,IAAI,EAAE,EAAE,CAAC;YACrC,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,EACtB,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,EACpB,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,EACpB,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;YACvB,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;YACb,IAAI,CAAC;gBACH,MAAM,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;gBAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC/E,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,IAAI,EAAE,CAAC;oBAC3C,IAAI,KAAK,IAAI,IAAI,EAAE,CAAC;wBAClB,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;wBACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;wBACb,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC;wBACb,OAAO;oBACT,CAAC;oBACD,EAAE,GAAG,EAAE,CAAC;oBACR,EAAE,GAAG,EAAE,CAAC;oBACR,EAAE,GAAG,EAAE,CAAC;oBACR,QAAQ,GAAG,IAAI,CAAC;oBAChB,MAAM;gBACR,CAAC;YACH,CAAC;YAAC,MAAM,CAAC;gBACP,mEAAmE;YACrE,CAAC;QACH,CAAC;QACD,IAAI,CAAC,QAAQ;YAAE,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IACxE,CAAC;IACD,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;AACvD,CAAC"}

package/dist/projection/experimental/deformation.d.ts

@@ -0,0 +1,23 @@
+import type { ProjectionPoint } from "./types.js";
+import type { VelocityGrid } from "./grids/velocity.js";
+/** Both methods take the ORIGINAL source/target epoch pair. inverse solves forward,
+ * rather than sampling the endpoint with a negated duration. XYZ are ECEF metres.
+ * Successful calls mutate XYZ; a failure restores the input. No per-point allocations.
+ */
+export type DeformationModel = {
+    forward(point: ProjectionPoint, sourceEpoch: number, targetEpoch: number): void;
+    inverse(point: ProjectionPoint, sourceEpoch: number, targetEpoch: number): void;
+};
+export type DeformationModelOptions = {
+    readonly grid: VelocityGrid;
+    /** Inclusive, application-reviewed validity interval in decimal years. */
+    readonly epochRange: readonly [number, number];
+    /** Must match the velocity model's geographic reference frame. Default WGS84. */
+    readonly ellipsoid?: {
+        readonly semiMajorAxis: number;
+        readonly flattening: number;
+    };
+};
+/** Time-invariant ENU velocity, source-sampled in fixed XYZ space. */
+export declare function createDeformationModel(options: DeformationModelOptions): DeformationModel;
+//# sourceMappingURL=deformation.d.ts.map

package/dist/projection/experimental/deformation.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"deformation.d.ts","sourceRoot":"","sources":["../../../src/projection/experimental/deformation.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,eAAe,EAAC,mBAAgB;AAC7C,OAAO,KAAK,EAAC,YAAY,EAAC,4BAAyB;AACnD;;;GAGG;AACH,MAAM,MAAM,gBAAgB,GAAG;IAC7B,OAAO,CAAC,KAAK,EAAE,eAAe,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,GAAG,IAAI,CAAC;IAChF,OAAO,CAAC,KAAK,EAAE,eAAe,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,MAAM,GAAG,IAAI,CAAC;CACjF,CAAC;AACF,MAAM,MAAM,uBAAuB,GAAG;IACpC,QAAQ,CAAC,IAAI,EAAE,YAAY,CAAC;IAC5B,0EAA0E;IAC1E,QAAQ,CAAC,UAAU,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,iFAAiF;IACjF,QAAQ,CAAC,SAAS,CAAC,EAAE;QAAC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;QAAC,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAA;KAAC,CAAC;CACpF,CAAC;AAEF,sEAAsE;AACtE,wBAAgB,sBAAsB,CAAC,OAAO,EAAE,uBAAuB,GAAG,gBAAgB,CAKzF"}

package/dist/projection/experimental/deformation.js

@@ -0,0 +1,13 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original ENU basis rotation and inverse, inspired by PROJ's documented deformation contract; no PROJ implementation is copied or forked. Geocentric conversion reuses the attributed proj4js adaptation in datum.ts.
+import { createSpatialDeformationModel } from "./spatial-deformation.js";
+/** Time-invariant ENU velocity, source-sampled in fixed XYZ space. */
+export function createDeformationModel(options) {
+    if (!options.grid || typeof options.grid.sample !== 'function')
+        throw new Error('Prepared velocity grid required');
+    const sample = options.grid.sample.bind(options.grid);
+    return createSpatialDeformationModel(options, sample, true);
+}
+//# sourceMappingURL=deformation.js.map

package/dist/projection/experimental/deformation.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"deformation.js","sourceRoot":"","sources":["../../../src/projection/experimental/deformation.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,yOAAyO;AACzO,OAAO,EAAC,6BAA6B,EAAC,iCAA8B;AAmBpE,sEAAsE;AACtE,MAAM,UAAU,sBAAsB,CAAC,OAAgC;IACrE,IAAI,CAAC,OAAO,CAAC,IAAI,IAAI,OAAO,OAAO,CAAC,IAAI,CAAC,MAAM,KAAK,UAAU;QAC5D,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;IACrD,MAAM,MAAM,GAAG,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IACtD,OAAO,6BAA6B,CAAC,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,CAAC;AAC9D,CAAC"}

package/dist/projection/experimental/ellipsoidal-gnomonic.d.ts

@@ -0,0 +1,10 @@
+import { createGeodesicEngine } from "../../ellipsoid-geodesic.js";
+import type { KernelParameters, Point } from "./kernel.js";
+export type EllipsoidalGnomonic = {
+    engine: ReturnType<typeof createGeodesicEngine>;
+    latitude: number;
+};
+export declare function prepareEllipsoidalGnomonic(base: KernelParameters): EllipsoidalGnomonic;
+export declare function forwardEllipsoidalGnomonic(base: KernelParameters, state: EllipsoidalGnomonic, p: Point): Point;
+export declare function inverseEllipsoidalGnomonic(base: KernelParameters, state: EllipsoidalGnomonic, p: Point): Point;
+//# sourceMappingURL=ellipsoidal-gnomonic.d.ts.map

package/dist/projection/experimental/ellipsoidal-gnomonic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoidal-gnomonic.d.ts","sourceRoot":"","sources":["../../../src/projection/experimental/ellipsoidal-gnomonic.ts"],"names":[],"mappings":"AAKA,OAAO,EAAC,oBAAoB,EAAsB,oCAAiC;AACnF,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,oBAAiB;AAItD,MAAM,MAAM,mBAAmB,GAAG;IAChC,MAAM,EAAE,UAAU,CAAC,OAAO,oBAAoB,CAAC,CAAC;IAChD,QAAQ,EAAE,MAAM,CAAC;CAClB,CAAC;AAEF,wBAAgB,0BAA0B,CAAC,IAAI,EAAE,gBAAgB,GAAG,mBAAmB,CAKtF;AAED,wBAAgB,0BAA0B,CACxC,IAAI,EAAE,gBAAgB,EACtB,KAAK,EAAE,mBAAmB,EAC1B,CAAC,EAAE,KAAK,GACP,KAAK,CAcP;AAED,wBAAgB,0BAA0B,CACxC,IAAI,EAAE,gBAAgB,EACtB,KAAK,EAAE,mBAAmB,EAC1B,CAAC,EAAE,KAAK,GACP,KAAK,CAoCP"}

package/dist/projection/experimental/ellipsoidal-gnomonic.js

@@ -0,0 +1,59 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) PROJ contributors
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Ellipsoidal gnomonic dispatch and inverse Newton equations adapted from PROJ 9.5.1 src/projections/gnom.cpp (Karney geodesic formulation). See ../../../projection/PROJ-LICENSE.txt. Geodesic evaluation uses the shared ellipsoid engine; public optional polygon/line facade code stays out of projection bundles.
+import { createGeodesicEngine, nonAreaGeodesicMask } from "../../ellipsoid-geodesic.js";
+import adjust_lon from "./common/adjust_lon.js";
+const degrees = 180 / Math.PI;
+export function prepareEllipsoidalGnomonic(base) {
+    return {
+        engine: createGeodesicEngine({ semiMajorAxis: 1, semiMinorAxis: Math.sqrt(1 - base.es) }),
+        latitude: base.lat0 * degrees
+    };
+}
+export function forwardEllipsoidalGnomonic(base, state, p) {
+    const result = state.engine.Inverse(state.latitude, 0, p.y * degrees, adjust_lon(p.x - base.long0, base.over) * degrees, nonAreaGeodesicMask);
+    if (!(result.M12 > 0))
+        throw new Error('Coordinate is outside the ellipsoidal gnomonic horizon');
+    const radius = (base.a * result.m12) / result.M12;
+    const azimuth = result.azi1 / degrees;
+    p.x = base.x0 + radius * Math.sin(azimuth);
+    p.y = base.y0 + radius * Math.cos(azimuth);
+    return p;
+}
+export function inverseEllipsoidalGnomonic(base, state, p) {
+    const x = (p.x - base.x0) / base.a;
+    const y = (p.y - base.y0) / base.a;
+    const radius = Math.hypot(x, y);
+    if (radius === 0) {
+        p.x = base.long0;
+        p.y = base.lat0;
+        return p;
+    }
+    const line = state.engine.Line(state.latitude, 0, Math.atan2(x, y) * degrees, nonAreaGeodesicMask);
+    const little = radius <= 1;
+    const target = little ? radius : 1 / radius;
+    let distance = Math.atan(radius);
+    // d(m/M)/ds = 1/M². Use its reciprocal near the horizon to avoid
+    // overflow, and evaluate once more at the converged distance.
+    for (let iteration = 0; iteration < 20; iteration++) {
+        const result = line.Position(distance, nonAreaGeodesicMask);
+        const step = little
+            ? (result.m12 - target * result.M12) * result.M12
+            : (target * result.m12 - result.M12) * result.m12;
+        if (!Number.isFinite(step))
+            break;
+        distance -= step;
+        if (Math.abs(step) < 1e-14) {
+            const final = line.Position(distance, nonAreaGeodesicMask);
+            if (!(final.M12 > 0))
+                break;
+            p.x = adjust_lon(base.long0 + final.lon2 / degrees, base.over);
+            p.y = final.lat2 / degrees;
+            return p;
+        }
+    }
+    throw new Error('Ellipsoidal gnomonic inverse did not converge');
+}
+//# sourceMappingURL=ellipsoidal-gnomonic.js.map

package/dist/projection/experimental/ellipsoidal-gnomonic.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoidal-gnomonic.js","sourceRoot":"","sources":["../../../src/projection/experimental/ellipsoidal-gnomonic.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,0DAA0D;AAC1D,4DAA4D;AAC5D,yUAAyU;AACzU,OAAO,EAAC,oBAAoB,EAAE,mBAAmB,EAAC,oCAAiC;AAEnF,OAAO,UAAU,+BAA4B;AAE7C,MAAM,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;AAM9B,MAAM,UAAU,0BAA0B,CAAC,IAAsB;IAC/D,OAAO;QACL,MAAM,EAAE,oBAAoB,CAAC,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAC,CAAC;QACvF,QAAQ,EAAE,IAAI,CAAC,IAAI,GAAG,OAAO;KAC9B,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,0BAA0B,CACxC,IAAsB,EACtB,KAA0B,EAC1B,CAAQ;IAER,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC,OAAO,CACjC,KAAK,CAAC,QAAQ,EACd,CAAC,EACD,CAAC,CAAC,CAAC,GAAG,OAAO,EACb,UAAU,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,OAAO,EACjD,mBAAmB,CACpB,CAAC;IACF,IAAI,CAAC,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IACjG,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;IAClD,MAAM,OAAO,GAAG,MAAM,CAAC,IAAI,GAAG,OAAO,CAAC;IACtC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;IAC3C,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;IAC3C,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,0BAA0B,CACxC,IAAsB,EACtB,KAA0B,EAC1B,CAAQ;IAER,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IACnC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAChC,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;QACjB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACjB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QAChB,OAAO,CAAC,CAAC;IACX,CAAC;IACD,MAAM,IAAI,GAAG,KAAK,CAAC,MAAM,CAAC,IAAI,CAC5B,KAAK,CAAC,QAAQ,EACd,CAAC,EACD,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,EAC1B,mBAAmB,CACpB,CAAC;IACF,MAAM,MAAM,GAAG,MAAM,IAAI,CAAC,CAAC;IAC3B,MAAM,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;IAC5C,IAAI,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACjC,iEAAiE;IACjE,8DAA8D;IAC9D,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,EAAE,EAAE,SAAS,EAAE,EAAE,CAAC;QACpD,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,QAAQ,EAAE,mBAAmB,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,MAAM;YACjB,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG;YACjD,CAAC,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;QACpD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YAAE,MAAM;QAClC,QAAQ,IAAI,IAAI,CAAC;QACjB,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,QAAQ,EAAE,mBAAmB,CAAC,CAAC;YAC3D,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC;gBAAE,MAAM;YAC5B,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,IAAI,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;YAC/D,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,GAAG,OAAO,CAAC;YAC3B,OAAO,CAAC,CAAC;QACX,CAAC;IACH,CAAC;IACD,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;AACnE,CAAC"}

package/dist/projection/experimental/engine-options.d.ts

@@ -0,0 +1,10 @@
+import type { ProjectionEngineOptions } from "../types.js";
+import type { ProjectionPlugin } from "./types.js";
+import type { ProjectionDescriptor } from "./projection-descriptor.js";
+type Registration = ProjectionPlugin | ProjectionDescriptor;
+/** Copy registration containers while retaining owned plugin and grid data. */
+export declare function snapshotEngineOptions<P extends Registration>(options: ProjectionEngineOptions<P>): ProjectionEngineOptions<P>;
+/** Custom execution registrations can reinterpret the same input definition. */
+export declare function hasCRSOverrides(options: ProjectionEngineOptions<Registration>): boolean;
+export {};
+//# sourceMappingURL=engine-options.d.ts.map

package/dist/projection/experimental/engine-options.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"engine-options.d.ts","sourceRoot":"","sources":["../../../src/projection/experimental/engine-options.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,uBAAuB,EAAC,oBAAiB;AACtD,OAAO,KAAK,EAAC,gBAAgB,EAAC,mBAAgB;AAC9C,OAAO,KAAK,EAAC,oBAAoB,EAAC,mCAAgC;AAClE,KAAK,YAAY,GAAG,gBAAgB,GAAG,oBAAoB,CAAC;AAE5D,+EAA+E;AAC/E,wBAAgB,qBAAqB,CAAC,CAAC,SAAS,YAAY,EAC1D,OAAO,EAAE,uBAAuB,CAAC,CAAC,CAAC,GAClC,uBAAuB,CAAC,CAAC,CAAC,CAU5B;AAED,gFAAgF;AAChF,wBAAgB,eAAe,CAAC,OAAO,EAAE,uBAAuB,CAAC,YAAY,CAAC,GAAG,OAAO,CAUvF"}

package/dist/projection/experimental/engine-options.js

@@ -0,0 +1,26 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Copy registration containers while retaining owned plugin and grid data. */
+export function snapshotEngineOptions(options) {
+    return {
+        ...options,
+        projections: options.projections && [...options.projections],
+        parsers: options.parsers && [...options.parsers],
+        datumCatalogs: options.datumCatalogs && [...options.datumCatalogs],
+        aliases: options.aliases && JSON.parse(JSON.stringify(options.aliases)),
+        datumGrids: options.datumGrids && { ...options.datumGrids },
+        verticalGrids: options.verticalGrids && { ...options.verticalGrids }
+    };
+}
+/** Custom execution registrations can reinterpret the same input definition. */
+export function hasCRSOverrides(options) {
+    return Boolean(options.projections ||
+        options.parsers ||
+        options.datumCatalogs ||
+        options.aliases ||
+        options.datumGrids ||
+        options.verticalGrids ||
+        options.resolver);
+}
+//# sourceMappingURL=engine-options.js.map

package/dist/projection/experimental/engine-options.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"engine-options.js","sourceRoot":"","sources":["../../../src/projection/experimental/engine-options.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAO5D,+EAA+E;AAC/E,MAAM,UAAU,qBAAqB,CACnC,OAAmC;IAEnC,OAAO;QACL,GAAG,OAAO;QACV,WAAW,EAAE,OAAO,CAAC,WAAW,IAAI,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC;QAC5D,OAAO,EAAE,OAAO,CAAC,OAAO,IAAI,CAAC,GAAG,OAAO,CAAC,OAAO,CAAC;QAChD,aAAa,EAAE,OAAO,CAAC,aAAa,IAAI,CAAC,GAAG,OAAO,CAAC,aAAa,CAAC;QAClE,OAAO,EAAE,OAAO,CAAC,OAAO,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC;QACvE,UAAU,EAAE,OAAO,CAAC,UAAU,IAAI,EAAC,GAAG,OAAO,CAAC,UAAU,EAAC;QACzD,aAAa,EAAE,OAAO,CAAC,aAAa,IAAI,EAAC,GAAG,OAAO,CAAC,aAAa,EAAC;KACnE,CAAC;AACJ,CAAC;AAED,gFAAgF;AAChF,MAAM,UAAU,eAAe,CAAC,OAA8C;IAC5E,OAAO,OAAO,CACZ,OAAO,CAAC,WAAW;QACjB,OAAO,CAAC,OAAO;QACf,OAAO,CAAC,aAAa;QACrB,OAAO,CAAC,OAAO;QACf,OAAO,CAAC,UAAU;QAClB,OAAO,CAAC,aAAa;QACrB,OAAO,CAAC,QAAQ,CACnB,CAAC;AACJ,CAAC"}