@math.gl/projection 0.0.0-stage → 5.0.0-alpha.12

package/dist/experimental/grids/vertical-geotiff.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical-geotiff.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/vertical-geotiff.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,YAAY,EAAC,mBAAgB;AAC1C,OAAO,KAAK,EAAC,mBAAmB,EAAE,uBAAuB,EAAC,4BAAyB;AACnF,YAAY,EACV,mBAAmB,EACnB,wBAAwB,EACxB,uBAAuB,EACxB,4BAAyB;AAE1B,2FAA2F;AAC3F,wBAAsB,uBAAuB,CAC3C,IAAI,EAAE,mBAAmB,GAAG,uBAAuB,GAClD,OAAO,CAAC,YAAY,CAAC,CAWvB"}

package/dist/experimental/grids/vertical-geotiff.js

@@ -0,0 +1,20 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adapter for PROJ's GTG conventions; no upstream decoder code is copied.
+import { prepareScalarGeoTIFFGrids } from "./scalar-geotiff.js";
+/** Prepare owned geoid grids; later nested images take precedence with nodata fallback. */
+export async function loadVerticalGeoTIFFGrid(tiff) {
+    const grids = await prepareScalarGeoTIFFGrids(tiff);
+    return Object.freeze({
+        getOffset(longitude, latitude) {
+            for (let i = grids.length - 1; i >= 0; i--) {
+                const value = grids[i].getOffset(longitude, latitude);
+                if (value !== undefined)
+                    return value;
+            }
+            return undefined;
+        }
+    });
+}
+//# sourceMappingURL=vertical-geotiff.js.map

package/dist/experimental/grids/vertical-geotiff.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical-geotiff.js","sourceRoot":"","sources":["../../../src/experimental/grids/vertical-geotiff.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,qGAAqG;AACrG,OAAO,EAAC,yBAAyB,EAAC,4BAAyB;AAS3D,2FAA2F;AAC3F,MAAM,CAAC,KAAK,UAAU,uBAAuB,CAC3C,IAAmD;IAEnD,MAAM,KAAK,GAAG,MAAM,yBAAyB,CAAC,IAAI,CAAC,CAAC;IACpD,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,SAAS,CAAC,SAAiB,EAAE,QAAgB;YAC3C,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3C,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;gBACtD,IAAI,KAAK,KAAK,SAAS;oBAAE,OAAO,KAAK,CAAC;YACxC,CAAC;YACD,OAAO,SAAS,CAAC;QACnB,CAAC;KACF,CAAC,CAAC;AACL,CAAC"}

package/dist/experimental/grids/vertical.cjs

@@ -0,0 +1,80 @@
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// dist/experimental/grids/vertical.js
+var vertical_exports = {};
+__export(vertical_exports, {
+  createGeoidGrid: () => createGeoidGrid,
+  createVerticalGrid: () => createVerticalGrid
+});
+module.exports = __toCommonJS(vertical_exports);
+function createVerticalGrid(options) {
+  const [west, south] = options.origin;
+  const [dx, dy] = options.step;
+  const [width, height] = options.size;
+  const east = west + (width - 1) * dx;
+  const north = south + (height - 1) * dy;
+  if (![west, south, dx, dy, east, north].every(Number.isFinite) || ![width, height].every((n) => Number.isSafeInteger(n) && n >= 2) || !(dx > 0 && dy > 0) || east - west > 360 || south < -90 || north > 90 || width * height !== options.offsets.length)
+    throw new Error("Invalid vertical grid geometry or node count");
+  const noData = options.noData;
+  const values = Float64Array.from(options.offsets, (value) => value === noData || !Number.isFinite(value) ? NaN : value);
+  return Object.freeze({
+    getOffset(longitude, latitude) {
+      if (!Number.isFinite(longitude) || !Number.isFinite(latitude))
+        throw new Error("Vertical grid coordinates must be finite");
+      let x = longitude * (180 / Math.PI);
+      const epsilon = 1e-10;
+      if (x < west - epsilon * dx || x > east + epsilon * dx)
+        x += 360 * Math.round(((west + east) / 2 - x) / 360);
+      let col = (x - west) / dx;
+      let row = (latitude * (180 / Math.PI) - south) / dy;
+      if (!Number.isFinite(col) || !Number.isFinite(row) || col < -epsilon || row < -epsilon || col > width - 1 + epsilon || row > height - 1 + epsilon)
+        return void 0;
+      col = Math.max(0, Math.min(width - 1, col));
+      row = Math.max(0, Math.min(height - 1, row));
+      const i = Math.min(width - 2, Math.floor(col));
+      const j = Math.min(height - 2, Math.floor(row));
+      const u = col - i, v = row - j;
+      let offset = 0;
+      for (let n = 0; n < 4; n++) {
+        const right = n % 2, upper = n > 1;
+        const weight = (right ? u : 1 - u) * (upper ? v : 1 - v);
+        if (!weight)
+          continue;
+        const value = values[(j + (upper ? 1 : 0)) * width + i + right];
+        if (!Number.isFinite(value))
+          return void 0;
+        offset += weight * value;
+      }
+      return offset;
+    }
+  });
+}
+function createGeoidGrid(geoid) {
+  if (!geoid || typeof geoid.getHeight !== "function")
+    throw new Error("Invalid geoid model");
+  return Object.freeze({
+    getOffset(longitude, latitude) {
+      if (!Number.isFinite(longitude) || !Number.isFinite(latitude))
+        throw new Error("Vertical grid coordinates must be finite");
+      const offset = geoid.getHeight(latitude * (180 / Math.PI), longitude * (180 / Math.PI));
+      return Number.isFinite(offset) ? offset : void 0;
+    }
+  });
+}
+//# sourceMappingURL=vertical.cjs.map

package/dist/experimental/grids/vertical.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../src/experimental/grids/vertical.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAmBA;;;;;;AACM,SAAU,mBAAmB,SAA4B;AAC7D,QAAM,CAAC,MAAM,KAAK,IAAI,QAAQ;AAC9B,QAAM,CAAC,IAAI,EAAE,IAAI,QAAQ;AACzB,QAAM,CAAC,OAAO,MAAM,IAAI,QAAQ;AAChC,QAAM,OAAO,QAAQ,QAAQ,KAAK;AAClC,QAAM,QAAQ,SAAS,SAAS,KAAK;AACrC,MACE,CAAC,CAAC,MAAM,OAAO,IAAI,IAAI,MAAM,KAAK,EAAE,MAAM,OAAO,QAAQ,KACzD,CAAC,CAAC,OAAO,MAAM,EAAE,MAAM,OAAK,OAAO,cAAc,CAAC,KAAK,KAAK,CAAC,KAC7D,EAAE,KAAK,KAAK,KAAK,MACjB,OAAO,OAAO,OACd,QAAQ,OACR,QAAQ,MACR,QAAQ,WAAW,QAAQ,QAAQ;AAEnC,UAAM,IAAI,MAAM,8CAA8C;AAChE,QAAM,SAAS,QAAQ;AACvB,QAAM,SAAS,aAAa,KAAK,QAAQ,SAAS,WAChD,UAAU,UAAU,CAAC,OAAO,SAAS,KAAK,IAAI,MAAM,KAAK;AAE3D,SAAO,OAAO,OAAO;IACnB,UAAU,WAAmB,UAAgB;AAC3C,UAAI,CAAC,OAAO,SAAS,SAAS,KAAK,CAAC,OAAO,SAAS,QAAQ;AAC1D,cAAM,IAAI,MAAM,0CAA0C;AAG5D,UAAI,IAAI,aAAa,MAAM,KAAK;AAChC,YAAM,UAAU;AAEhB,UAAI,IAAI,OAAO,UAAU,MAAM,IAAI,OAAO,UAAU;AAClD,aAAK,MAAM,KAAK,QAAQ,OAAO,QAAQ,IAAI,KAAK,GAAG;AACrD,UAAI,OAAO,IAAI,QAAQ;AACvB,UAAI,OAAO,YAAY,MAAM,KAAK,MAAM,SAAS;AACjD,UACE,CAAC,OAAO,SAAS,GAAG,KACpB,CAAC,OAAO,SAAS,GAAG,KACpB,MAAM,CAAC,WACP,MAAM,CAAC,WACP,MAAM,QAAQ,IAAI,WAClB,MAAM,SAAS,IAAI;AAEnB,eAAO;AACT,YAAM,KAAK,IAAI,GAAG,KAAK,IAAI,QAAQ,GAAG,GAAG,CAAC;AAC1C,YAAM,KAAK,IAAI,GAAG,KAAK,IAAI,SAAS,GAAG,GAAG,CAAC;AAC3C,YAAM,IAAI,KAAK,IAAI,QAAQ,GAAG,KAAK,MAAM,GAAG,CAAC;AAC7C,YAAM,IAAI,KAAK,IAAI,SAAS,GAAG,KAAK,MAAM,GAAG,CAAC;AAC9C,YAAM,IAAI,MAAM,GACd,IAAI,MAAM;AACZ,UAAI,SAAS;AACb,eAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,cAAM,QAAQ,IAAI,GAChB,QAAQ,IAAI;AACd,cAAM,UAAU,QAAQ,IAAI,IAAI,MAAM,QAAQ,IAAI,IAAI;AACtD,YAAI,CAAC;AAAQ;AACb,cAAM,QAAQ,QAAQ,KAAK,QAAQ,IAAI,MAAM,QAAQ,IAAI,KAAK;AAC9D,YAAI,CAAC,OAAO,SAAS,KAAK;AAAG,iBAAO;AACpC,kBAAU,SAAS;MACrB;AACA,aAAO;IACT;GACD;AACH;AAKM,SAAU,gBAAgB,OAE/B;AACC,MAAI,CAAC,SAAS,OAAO,MAAM,cAAc;AAAY,UAAM,IAAI,MAAM,qBAAqB;AAC1F,SAAO,OAAO,OAAO;IACnB,UAAU,WAAmB,UAAgB;AAC3C,UAAI,CAAC,OAAO,SAAS,SAAS,KAAK,CAAC,OAAO,SAAS,QAAQ;AAC1D,cAAM,IAAI,MAAM,0CAA0C;AAC5D,YAAM,SAAS,MAAM,UAAU,YAAY,MAAM,KAAK,KAAK,aAAa,MAAM,KAAK,GAAG;AACtF,aAAO,OAAO,SAAS,MAAM,IAAI,SAAS;IAC5C;GACD;AACH;",
+  "names": []
+}

package/dist/experimental/grids/vertical.d.ts

@@ -0,0 +1,22 @@
+import type { VerticalGrid } from "./types.js";
+export type { VerticalGrid, VerticalGridCollection } from "./types.js";
+export type VerticalGridOptions = {
+    /** Southwest node, [longitude, latitude] in degrees. */
+    origin: readonly [number, number];
+    /** Positive [longitude, latitude] node spacing in degrees. */
+    step: readonly [number, number];
+    size: readonly [number, number];
+    /** Row-major south-to-north, west-to-east geoid undulations in metres. Copied. */
+    offsets: ArrayLike<number>;
+    /** Non-finite values are also treated as nodata. */
+    noData?: number;
+};
+/** Prepare an owned, synchronous grid with inclusive edges and no extrapolation. */
+export declare function createVerticalGrid(options: VerticalGridOptions): VerticalGrid;
+/** Adapt an already loaded @math.gl/geoid Geoid without importing its implementation.
+ * The supplied model owns its data and interpolation choice; do not mutate it while in use.
+ */
+export declare function createGeoidGrid(geoid: {
+    getHeight(latitude: number, longitude: number): number;
+}): VerticalGrid;
+//# sourceMappingURL=vertical.d.ts.map

package/dist/experimental/grids/vertical.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/vertical.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,YAAY,EAAC,mBAAgB;AAC1C,YAAY,EAAC,YAAY,EAAE,sBAAsB,EAAC,mBAAgB;AAElE,MAAM,MAAM,mBAAmB,GAAG;IAChC,wDAAwD;IACxD,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAClC,8DAA8D;IAC9D,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAChC,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAChC,kFAAkF;IAClF,OAAO,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;IAC3B,oDAAoD;IACpD,MAAM,CAAC,EAAE,MAAM,CAAC;CACjB,CAAC;AAEF,oFAAoF;AACpF,wBAAgB,kBAAkB,CAAC,OAAO,EAAE,mBAAmB,GAAG,YAAY,CA6D7E;AAED;;GAEG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE;IACrC,SAAS,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM,CAAC;CACxD,GAAG,YAAY,CAUf"}

package/dist/experimental/grids/vertical.js

@@ -0,0 +1,73 @@
+/** Prepare an owned, synchronous grid with inclusive edges and no extrapolation. */
+export function createVerticalGrid(options) {
+    const [west, south] = options.origin;
+    const [dx, dy] = options.step;
+    const [width, height] = options.size;
+    const east = west + (width - 1) * dx;
+    const north = south + (height - 1) * dy;
+    if (![west, south, dx, dy, east, north].every(Number.isFinite) ||
+        ![width, height].every(n => Number.isSafeInteger(n) && n >= 2) ||
+        !(dx > 0 && dy > 0) ||
+        east - west > 360 ||
+        south < -90 ||
+        north > 90 ||
+        width * height !== options.offsets.length)
+        throw new Error('Invalid vertical grid geometry or node count');
+    const noData = options.noData;
+    const values = Float64Array.from(options.offsets, value => value === noData || !Number.isFinite(value) ? NaN : value);
+    return Object.freeze({
+        getOffset(longitude, latitude) {
+            if (!Number.isFinite(longitude) || !Number.isFinite(latitude))
+                throw new Error('Vertical grid coordinates must be finite');
+            // Choose the equivalent longitude nearest the grid's centre, including
+            // grids expressed in 0..360 degrees or crossing the antimeridian.
+            let x = longitude * (180 / Math.PI);
+            const epsilon = 1e-10;
+            // Preserve either explicit edge of a grid with duplicate seam nodes.
+            if (x < west - epsilon * dx || x > east + epsilon * dx)
+                x += 360 * Math.round(((west + east) / 2 - x) / 360);
+            let col = (x - west) / dx;
+            let row = (latitude * (180 / Math.PI) - south) / dy;
+            if (!Number.isFinite(col) ||
+                !Number.isFinite(row) ||
+                col < -epsilon ||
+                row < -epsilon ||
+                col > width - 1 + epsilon ||
+                row > height - 1 + epsilon)
+                return undefined;
+            col = Math.max(0, Math.min(width - 1, col));
+            row = Math.max(0, Math.min(height - 1, row));
+            const i = Math.min(width - 2, Math.floor(col));
+            const j = Math.min(height - 2, Math.floor(row));
+            const u = col - i, v = row - j;
+            let offset = 0;
+            for (let n = 0; n < 4; n++) {
+                const right = n % 2, upper = n > 1;
+                const weight = (right ? u : 1 - u) * (upper ? v : 1 - v);
+                if (!weight)
+                    continue;
+                const value = values[(j + (upper ? 1 : 0)) * width + i + right];
+                if (!Number.isFinite(value))
+                    return undefined;
+                offset += weight * value;
+            }
+            return offset;
+        }
+    });
+}
+/** Adapt an already loaded @math.gl/geoid Geoid without importing its implementation.
+ * The supplied model owns its data and interpolation choice; do not mutate it while in use.
+ */
+export function createGeoidGrid(geoid) {
+    if (!geoid || typeof geoid.getHeight !== 'function')
+        throw new Error('Invalid geoid model');
+    return Object.freeze({
+        getOffset(longitude, latitude) {
+            if (!Number.isFinite(longitude) || !Number.isFinite(latitude))
+                throw new Error('Vertical grid coordinates must be finite');
+            const offset = geoid.getHeight(latitude * (180 / Math.PI), longitude * (180 / Math.PI));
+            return Number.isFinite(offset) ? offset : undefined;
+        }
+    });
+}
+//# sourceMappingURL=vertical.js.map

package/dist/experimental/grids/vertical.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical.js","sourceRoot":"","sources":["../../../src/experimental/grids/vertical.ts"],"names":[],"mappings":"AAmBA,oFAAoF;AACpF,MAAM,UAAU,kBAAkB,CAAC,OAA4B;IAC7D,MAAM,CAAC,IAAI,EAAE,KAAK,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;IACrC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC;IAC9B,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC;IACrC,MAAM,IAAI,GAAG,IAAI,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IACrC,MAAM,KAAK,GAAG,KAAK,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IACxC,IACE,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAC1D,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9D,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;QACnB,IAAI,GAAG,IAAI,GAAG,GAAG;QACjB,KAAK,GAAG,CAAC,EAAE;QACX,KAAK,GAAG,EAAE;QACV,KAAK,GAAG,MAAM,KAAK,OAAO,CAAC,OAAO,CAAC,MAAM;QAEzC,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAClE,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAC9B,MAAM,MAAM,GAAG,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,EAAE,KAAK,CAAC,EAAE,CACxD,KAAK,KAAK,MAAM,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAC1D,CAAC;IACF,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,SAAS,CAAC,SAAiB,EAAE,QAAgB;YAC3C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;gBAC3D,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,uEAAuE;YACvE,kEAAkE;YAClE,IAAI,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;YACpC,MAAM,OAAO,GAAG,KAAK,CAAC;YACtB,qEAAqE;YACrE,IAAI,CAAC,GAAG,IAAI,GAAG,OAAO,GAAG,EAAE,IAAI,CAAC,GAAG,IAAI,GAAG,OAAO,GAAG,EAAE;gBACpD,CAAC,IAAI,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;YACvD,IAAI,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;YAC1B,IAAI,GAAG,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC;YACpD,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;gBACrB,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;gBACrB,GAAG,GAAG,CAAC,OAAO;gBACd,GAAG,GAAG,CAAC,OAAO;gBACd,GAAG,GAAG,KAAK,GAAG,CAAC,GAAG,OAAO;gBACzB,GAAG,GAAG,MAAM,GAAG,CAAC,GAAG,OAAO;gBAE1B,OAAO,SAAS,CAAC;YACnB,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAC5C,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAC7C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;YAChD,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,EACf,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YACd,IAAI,MAAM,GAAG,CAAC,CAAC;YACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,EACjB,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;gBAChB,MAAM,MAAM,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACzD,IAAI,CAAC,MAAM;oBAAE,SAAS;gBACtB,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;gBAChE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;oBAAE,OAAO,SAAS,CAAC;gBAC9C,MAAM,IAAI,MAAM,GAAG,KAAK,CAAC;YAC3B,CAAC;YACD,OAAO,MAAM,CAAC;QAChB,CAAC;KACF,CAAC,CAAC;AACL,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,eAAe,CAAC,KAE/B;IACC,IAAI,CAAC,KAAK,IAAI,OAAO,KAAK,CAAC,SAAS,KAAK,UAAU;QAAE,MAAM,IAAI,KAAK,CAAC,qBAAqB,CAAC,CAAC;IAC5F,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,SAAS,CAAC,SAAiB,EAAE,QAAgB;YAC3C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;gBAC3D,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,MAAM,MAAM,GAAG,KAAK,CAAC,SAAS,CAAC,QAAQ,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,SAAS,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;YACxF,OAAO,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC;QACtD,CAAC;KACF,CAAC,CAAC;AACL,CAAC"}

package/dist/experimental/helmert-flat.d.ts

@@ -0,0 +1,4 @@
+import type { NumericFlatOperation } from "./numeric-flat.js";
+/** Setup-owned coefficients: small-angle rx/ry/rz or row-major exact matrix, then dx/dy/dz/scale. */
+export declare function createHelmertFlat(coefficients: readonly number[], inverse: boolean, exact: boolean): NumericFlatOperation;
+//# sourceMappingURL=helmert-flat.d.ts.map

package/dist/experimental/helmert-flat.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"helmert-flat.d.ts","sourceRoot":"","sources":["../../src/experimental/helmert-flat.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,oBAAoB,EAAC,0BAAuB;AAEzD,qGAAqG;AACrG,wBAAgB,iBAAiB,CAC/B,YAAY,EAAE,SAAS,MAAM,EAAE,EAC/B,OAAO,EAAE,OAAO,EAChB,KAAK,EAAE,OAAO,GACb,oBAAoB,CA8DtB"}

package/dist/experimental/helmert-flat.js

@@ -0,0 +1,53 @@
+/** Setup-owned coefficients: small-angle rx/ry/rz or row-major exact matrix, then dx/dy/dz/scale. */
+export function createHelmertFlat(coefficients, inverse, exact) {
+    const dx = coefficients[9], dy = coefficients[10], dz = coefficients[11], scale = coefficients[12];
+    const a = coefficients[0], b = coefficients[1], c = coefficients[2], d = coefficients[3], e = coefficients[4], f = coefficients[5], g = coefficients[6], h = coefficients[7], j = coefficients[8];
+    return (coordinates, dimension, epochs) => {
+        const float32 = coordinates instanceof Float32Array;
+        const epochBuffer = typeof epochs === 'number' ? undefined : epochs;
+        for (let offset = 0, record = 0; offset < coordinates.length; offset += dimension, record++) {
+            let x = coordinates[offset], y = coordinates[offset + 1], z = coordinates[offset + 2];
+            if (epochBuffer && !Number.isFinite(epochBuffer[record]))
+                throw new Error('Coordinate epoch must be a finite decimal year');
+            if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))
+                throw new Error('Pipeline coordinate must be finite');
+            if (inverse) {
+                x = (x - dx) / scale;
+                y = (y - dy) / scale;
+                z = (z - dz) / scale;
+            }
+            let outX, outY, outZ;
+            if (exact) {
+                if (inverse) {
+                    outX = a * x + d * y + g * z;
+                    outY = b * x + e * y + h * z;
+                    outZ = c * x + f * y + j * z;
+                }
+                else {
+                    outX = scale * (a * x + b * y + c * z) + dx;
+                    outY = scale * (d * x + e * y + f * z) + dy;
+                    outZ = scale * (g * x + h * y + j * z) + dz;
+                }
+            }
+            else if (inverse) {
+                outX = x + c * y - b * z;
+                outY = -c * x + y + a * z;
+                outZ = b * x - a * y + z;
+            }
+            else {
+                outX = scale * (x - c * y + b * z) + dx;
+                outY = scale * (c * x + y - a * z) + dy;
+                outZ = scale * (-b * x + a * y + z) + dz;
+            }
+            if (!Number.isFinite(outX) || !Number.isFinite(outY) || !Number.isFinite(outZ))
+                throw new Error('Pipeline coordinate must be finite');
+            if (float32 &&
+                Math.max(Math.abs(outX), Math.abs(outY), Math.abs(outZ)) > 3.4028234663852886e38)
+                throw new Error('Pipeline output exceeds Float32 range');
+            coordinates[offset] = outX;
+            coordinates[offset + 1] = outY;
+            coordinates[offset + 2] = outZ;
+        }
+    };
+}
+//# sourceMappingURL=helmert-flat.js.map

package/dist/experimental/helmert-flat.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"helmert-flat.js","sourceRoot":"","sources":["../../src/experimental/helmert-flat.ts"],"names":[],"mappings":"AAQA,qGAAqG;AACrG,MAAM,UAAU,iBAAiB,CAC/B,YAA+B,EAC/B,OAAgB,EAChB,KAAc;IAEd,MAAM,EAAE,GAAG,YAAY,CAAC,CAAC,CAAC,EACxB,EAAE,GAAG,YAAY,CAAC,EAAE,CAAC,EACrB,EAAE,GAAG,YAAY,CAAC,EAAE,CAAC,EACrB,KAAK,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACvB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EACnB,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;IACtB,OAAO,CAAC,WAAW,EAAE,SAAS,EAAE,MAAM,EAAE,EAAE;QACxC,MAAM,OAAO,GAAG,WAAW,YAAY,YAAY,CAAC;QACpD,MAAM,WAAW,GAAG,OAAO,MAAM,KAAK,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAM,CAAC;QACpE,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,WAAW,CAAC,MAAM,EAAE,MAAM,IAAI,SAAS,EAAE,MAAM,EAAE,EAAE,CAAC;YAC5F,IAAI,CAAC,GAAG,WAAW,CAAC,MAAM,CAAC,EACzB,CAAC,GAAG,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,EAC3B,CAAC,GAAG,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;YAC9B,IAAI,WAAW,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC;gBACtD,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;YACpE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;gBACnE,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;YACxD,IAAI,OAAO,EAAE,CAAC;gBACZ,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;gBACrB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;gBACrB,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC;YACvB,CAAC;YACD,IAAI,IAAY,EAAE,IAAY,EAAE,IAAY,CAAC;YAC7C,IAAI,KAAK,EAAE,CAAC;gBACV,IAAI,OAAO,EAAE,CAAC;oBACZ,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAC7B,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAC7B,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC/B,CAAC;qBAAM,CAAC;oBACN,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;oBAC5C,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;oBAC5C,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;gBAC9C,CAAC;YACH,CAAC;iBAAM,IAAI,OAAO,EAAE,CAAC;gBACnB,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACzB,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC1B,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC3B,CAAC;iBAAM,CAAC;gBACN,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;gBACxC,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;gBACxC,IAAI,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAC3C,CAAC;YACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;gBAC5E,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;YACxD,IACE,OAAO;gBACP,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,qBAAqB;gBAEhF,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;YAC3D,WAAW,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC;YAC3B,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;YAC/B,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;QACjC,CAAC;IACH,CAAC,CAAC;AACJ,CAAC"}