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

package/dist/experimental/grids/ntv2.cjs

@@ -0,0 +1,181 @@
+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/ntv2.js
+var ntv2_exports = {};
+__export(ntv2_exports, {
+  parseNTv2Grid: () => parseNTv2Grid
+});
+module.exports = __toCommonJS(ntv2_exports);
+
+// dist/experimental/parameters.js
+var DEGREES_TO_RADIANS = Math.PI / 180;
+function wrapLongitude(longitude) {
+  if (Math.abs(longitude) <= Math.PI)
+    return longitude;
+  const circle = 2 * Math.PI;
+  return ((longitude + Math.PI) % circle + circle) % circle - Math.PI;
+}
+
+// dist/experimental/grids/grid.js
+var TOLERANCE = 1e-12;
+function interpolate(grid, x, y, output, clamp = false) {
+  let col = (x - grid.origin[0]) / grid.step[0], row = (y - grid.origin[1]) / grid.step[1];
+  const width = grid.size[0], height = grid.size[1];
+  if (!clamp && (col < -1e-10 || row < -1e-10 || col > width - 1 + 1e-10 || row > height - 1 + 1e-10))
+    return false;
+  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)), j = Math.min(height - 2, Math.floor(row));
+  const u = col - i, v = row - j;
+  let dx = 0, dy = 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 === 0)
+      continue;
+    const index = ((j + (upper ? 1 : 0)) * width + i + right) * 2;
+    const xShift = grid.shifts[index], yShift = grid.shifts[index + 1];
+    if (!Number.isFinite(xShift) || !Number.isFinite(yShift))
+      return false;
+    dx += weight * xShift;
+    dy += weight * yShift;
+  }
+  output.x = dx;
+  output.y = dy;
+  return true;
+}
+function shiftSubgrid(grid, x, y, inverse, output) {
+  if (!inverse) {
+    if (!interpolate(grid, x, y, output))
+      return false;
+    output.x += x;
+    output.y += y;
+    return true;
+  }
+  for (let i = 0; i < 2; i++) {
+    const value = i === 0 ? x : y;
+    if (value < grid.origin[i] - grid.maxShift[i] - TOLERANCE || value > grid.origin[i] + (grid.size[i] - 1) * grid.step[i] + grid.maxShift[i] + TOLERANCE)
+      return false;
+  }
+  if (!interpolate(grid, x, y, output, true))
+    return false;
+  let guessX = x - output.x, guessY = y - output.y;
+  for (let iteration = 0; iteration < 20; iteration++) {
+    if (!interpolate(grid, guessX, guessY, output, true))
+      return false;
+    const dx = x - guessX - output.x, dy = y - guessY - output.y;
+    guessX += dx;
+    guessY += dy;
+    if (Math.max(Math.abs(dx), Math.abs(dy)) <= TOLERANCE) {
+      if (!interpolate(grid, guessX, guessY, output))
+        return false;
+      output.x = wrapLongitude(guessX);
+      output.y = guessY;
+      return true;
+    }
+  }
+  throw new Error("Inverse datum grid shift did not converge");
+}
+function createDatumGrid(subgrids) {
+  if (!subgrids.length)
+    throw new Error("Datum grid requires at least one subgrid");
+  const prepared = subgrids.map((grid) => {
+    const [width, height] = grid.size;
+    if (![width, height].every((n) => Number.isSafeInteger(n) && n >= 2) || width * height * 2 !== grid.shifts.length || !grid.origin.every(Number.isFinite) || !grid.step.every((n) => Number.isFinite(n) && n > 0))
+      throw new Error("Invalid datum subgrid geometry or node count");
+    const maxShift = [0, 0];
+    for (let index = 0; index < grid.shifts.length; index++) {
+      const value = grid.shifts[index];
+      if (Number.isFinite(value))
+        maxShift[index % 2] = Math.max(maxShift[index % 2], Math.abs(value));
+    }
+    return { ...grid, maxShift };
+  });
+  const shiftInPlace = (point, inverse) => {
+    const longitude = point.x, latitude = point.y;
+    if (!Number.isFinite(longitude) || !Number.isFinite(latitude))
+      throw new Error("Grid coordinates must be finite");
+    for (const grid of prepared) {
+      if (shiftSubgrid(grid, -longitude, latitude, inverse, point)) {
+        point.x = -point.x;
+        return true;
+      }
+    }
+    point.x = longitude;
+    point.y = latitude;
+    return false;
+  };
+  const tupleScratch = { x: NaN, y: NaN };
+  return Object.freeze({
+    subgridCount: prepared.length,
+    shiftInPlace,
+    shift(longitude, latitude, inverse) {
+      tupleScratch.x = longitude;
+      tupleScratch.y = latitude;
+      return shiftInPlace(tupleScratch, inverse) ? [tupleScratch.x, tupleScratch.y] : void 0;
+    }
+  });
+}
+
+// dist/experimental/grids/ntv2.js
+var ARC_SECOND = Math.PI / 648e3;
+function parseNTv2Grid(data, options = {}) {
+  const view = new DataView(data);
+  if (view.byteLength < 176)
+    throw new Error("Truncated NTv2 overview header");
+  const little = view.getInt32(8, true) === 11;
+  if (!little && view.getInt32(8, false) !== 11)
+    throw new Error("Invalid NTv2 header or byte order");
+  const integer = (offset2) => view.getInt32(offset2, little);
+  const real = (offset2) => view.getFloat64(offset2, little);
+  const label = (offset2) => String.fromCharCode(...new Uint8Array(data, offset2, 8)).replace(/\0/g, "").trim();
+  if (label(0) !== "NUM_OREC" || integer(24) !== 11 || label(56) !== "SECONDS")
+    throw new Error("Unsupported NTv2 header or shift units");
+  const count = integer(40), stride = options.includeErrorFields === false ? 8 : 16;
+  if (count < 1 || count > Math.floor((view.byteLength - 176) / (176 + 4 * stride)))
+    throw new Error("Invalid NTv2 subgrid count");
+  let offset = 176;
+  const subgrids = [];
+  for (let n = 0; n < count; n++) {
+    if (offset + 176 > view.byteLength)
+      throw new Error("Truncated NTv2 subgrid header");
+    const lowerLat = real(offset + 72), upperLat = real(offset + 88), lowerLon = real(offset + 104), upperLon = real(offset + 120);
+    const latStep = real(offset + 136), lonStep = real(offset + 152), nodes = integer(offset + 168);
+    const width = 1 + (upperLon - lowerLon) / lonStep, height = 1 + (upperLat - lowerLat) / latStep;
+    if (![width, height].every((v) => Number.isFinite(v) && v >= 2 && Math.abs(v - Math.round(v)) < 1e-7) || !(lonStep > 0 && latStep > 0) || nodes !== Math.round(width) * Math.round(height))
+      throw new Error("Invalid NTv2 subgrid dimensions");
+    const start = offset + 176, end = start + nodes * stride;
+    if (end > view.byteLength)
+      throw new Error("Truncated NTv2 node records");
+    const shifts = new Float64Array(nodes * 2);
+    for (let i = 0; i < nodes; i++) {
+      shifts[i * 2] = view.getFloat32(start + i * stride + 4, little) * ARC_SECOND;
+      shifts[i * 2 + 1] = view.getFloat32(start + i * stride, little) * ARC_SECOND;
+    }
+    subgrids.push({
+      origin: [lowerLon * ARC_SECOND, lowerLat * ARC_SECOND],
+      step: [lonStep * ARC_SECOND, latStep * ARC_SECOND],
+      size: [Math.round(width), Math.round(height)],
+      shifts
+    });
+    offset = end;
+  }
+  return createDatumGrid(subgrids);
+}
+//# sourceMappingURL=ntv2.cjs.map

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

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../src/experimental/grids/ntv2.ts", "../../../src/experimental/parameters.ts", "../../../src/experimental/grids/grid.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;;;;;;;ACMO,IAAM,qBAAqB,KAAK,KAAK;AAwBtC,SAAU,cAAc,WAAiB;AAC7C,MAAI,KAAK,IAAI,SAAS,KAAK,KAAK;AAAI,WAAO;AAC3C,QAAM,SAAS,IAAI,KAAK;AACxB,WAAW,YAAY,KAAK,MAAM,SAAU,UAAU,SAAU,KAAK;AACvE;;;AC1BA,IAAM,YAAY;AAIlB,SAAS,YACP,MACA,GACA,GACA,QACA,QAAQ,OAAK;AAEb,MAAI,OAAO,IAAI,KAAK,OAAO,CAAC,KAAK,KAAK,KAAK,CAAC,GAC1C,OAAO,IAAI,KAAK,OAAO,CAAC,KAAK,KAAK,KAAK,CAAC;AAC1C,QAAM,QAAQ,KAAK,KAAK,CAAC,GACvB,SAAS,KAAK,KAAK,CAAC;AACtB,MACE,CAAC,UACA,MAAM,UAAU,MAAM,UAAU,MAAM,QAAQ,IAAI,SAAS,MAAM,SAAS,IAAI;AAE/E,WAAO;AACT,QAAM,KAAK,IAAI,GAAG,KAAK,IAAI,QAAQ,GAAG,GAAG,CAAC;AAC1C,QAAM,KAAK,IAAI,GAAG,KAAK,IAAI,SAAS,GAAG,GAAG,CAAC;AAC3C,QAAM,IAAI,KAAK,IAAI,QAAQ,GAAG,KAAK,MAAM,GAAG,CAAC,GAC3C,IAAI,KAAK,IAAI,SAAS,GAAG,KAAK,MAAM,GAAG,CAAC;AAC1C,QAAM,IAAI,MAAM,GACd,IAAI,MAAM;AACZ,MAAI,KAAK,GACP,KAAK;AACP,WAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,UAAM,QAAQ,IAAI,GAChB,QAAQ,IAAI;AACd,UAAM,UAAU,QAAQ,IAAI,IAAI,MAAM,QAAQ,IAAI,IAAI;AACtD,QAAI,WAAW;AAAG;AAClB,UAAM,UAAU,KAAK,QAAQ,IAAI,MAAM,QAAQ,IAAI,SAAS;AAC5D,UAAM,SAAS,KAAK,OAAO,KAAK,GAC9B,SAAS,KAAK,OAAO,QAAQ,CAAC;AAChC,QAAI,CAAC,OAAO,SAAS,MAAM,KAAK,CAAC,OAAO,SAAS,MAAM;AAAG,aAAO;AACjE,UAAM,SAAS;AACf,UAAM,SAAS;EACjB;AACA,SAAO,IAAI;AACX,SAAO,IAAI;AACX,SAAO;AACT;AACA,SAAS,aACP,MACA,GACA,GACA,SACA,QAAa;AAEb,MAAI,CAAC,SAAS;AACZ,QAAI,CAAC,YAAY,MAAM,GAAG,GAAG,MAAM;AAAG,aAAO;AAC7C,WAAO,KAAK;AACZ,WAAO,KAAK;AACZ,WAAO;EACT;AAEA,WAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,UAAM,QAAQ,MAAM,IAAI,IAAI;AAC5B,QACE,QAAQ,KAAK,OAAO,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI,aAC5C,QAAQ,KAAK,OAAO,CAAC,KAAK,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,KAAK,CAAC,IAAI,KAAK,SAAS,CAAC,IAAI;AAEhF,aAAO;EACX;AACA,MAAI,CAAC,YAAY,MAAM,GAAG,GAAG,QAAQ,IAAI;AAAG,WAAO;AACnD,MAAI,SAAS,IAAI,OAAO,GACtB,SAAS,IAAI,OAAO;AACtB,WAAS,YAAY,GAAG,YAAY,IAAI,aAAa;AACnD,QAAI,CAAC,YAAY,MAAM,QAAQ,QAAQ,QAAQ,IAAI;AAAG,aAAO;AAC7D,UAAM,KAAK,IAAI,SAAS,OAAO,GAC7B,KAAK,IAAI,SAAS,OAAO;AAC3B,cAAU;AACV,cAAU;AACV,QAAI,KAAK,IAAI,KAAK,IAAI,EAAE,GAAG,KAAK,IAAI,EAAE,CAAC,KAAK,WAAW;AAErD,UAAI,CAAC,YAAY,MAAM,QAAQ,QAAQ,MAAM;AAAG,eAAO;AACvD,aAAO,IAAI,cAAc,MAAM;AAC/B,aAAO,IAAI;AACX,aAAO;IACT;EACF;AACA,QAAM,IAAI,MAAM,2CAA2C;AAC7D;AAEM,SAAU,gBAAgB,UAAmB;AACjD,MAAI,CAAC,SAAS;AAAQ,UAAM,IAAI,MAAM,0CAA0C;AAChF,QAAM,WAA8B,SAAS,IAAI,UAAO;AACtD,UAAM,CAAC,OAAO,MAAM,IAAI,KAAK;AAC7B,QACE,CAAC,CAAC,OAAO,MAAM,EAAE,MAAM,OAAK,OAAO,cAAc,CAAC,KAAK,KAAK,CAAC,KAC7D,QAAQ,SAAS,MAAM,KAAK,OAAO,UACnC,CAAC,KAAK,OAAO,MAAM,OAAO,QAAQ,KAClC,CAAC,KAAK,KAAK,MAAM,OAAK,OAAO,SAAS,CAAC,KAAK,IAAI,CAAC;AAEjD,YAAM,IAAI,MAAM,8CAA8C;AAChE,UAAM,WAA6B,CAAC,GAAG,CAAC;AACxC,aAAS,QAAQ,GAAG,QAAQ,KAAK,OAAO,QAAQ,SAAS;AACvD,YAAM,QAAQ,KAAK,OAAO,KAAK;AAC/B,UAAI,OAAO,SAAS,KAAK;AACvB,iBAAS,QAAQ,CAAC,IAAI,KAAK,IAAI,SAAS,QAAQ,CAAC,GAAG,KAAK,IAAI,KAAK,CAAC;IACvE;AACA,WAAO,EAAC,GAAG,MAAM,SAAQ;EAC3B,CAAC;AACD,QAAM,eAAe,CAAC,OAAc,YAA6B;AAC/D,UAAM,YAAY,MAAM,GACtB,WAAW,MAAM;AACnB,QAAI,CAAC,OAAO,SAAS,SAAS,KAAK,CAAC,OAAO,SAAS,QAAQ;AAC1D,YAAM,IAAI,MAAM,iCAAiC;AACnD,eAAW,QAAQ,UAAU;AAC3B,UAAI,aAAa,MAAM,CAAC,WAAW,UAAU,SAAS,KAAK,GAAG;AAC5D,cAAM,IAAI,CAAC,MAAM;AACjB,eAAO;MACT;IACF;AACA,UAAM,IAAI;AACV,UAAM,IAAI;AACV,WAAO;EACT;AAGA,QAAM,eAAsB,EAAC,GAAG,KAAK,GAAG,IAAG;AAC3C,SAAO,OAAO,OAAO;IACnB,cAAc,SAAS;IACvB;IACA,MAAM,WAAmB,UAAkB,SAAgB;AACzD,mBAAa,IAAI;AACjB,mBAAa,IAAI;AACjB,aAAO,aAAa,cAAc,OAAO,IAAI,CAAC,aAAa,GAAG,aAAa,CAAC,IAAI;IAClF;GACD;AACH;;;AFhIA,IAAM,aAAa,KAAK,KAAK;AAEvB,SAAU,cAAc,MAAmB,UAA2B,CAAA,GAAE;AAC5E,QAAM,OAAO,IAAI,SAAS,IAAI;AAC9B,MAAI,KAAK,aAAa;AAAK,UAAM,IAAI,MAAM,gCAAgC;AAC3E,QAAM,SAAS,KAAK,SAAS,GAAG,IAAI,MAAM;AAC1C,MAAI,CAAC,UAAU,KAAK,SAAS,GAAG,KAAK,MAAM;AACzC,UAAM,IAAI,MAAM,mCAAmC;AACrD,QAAM,UAAU,CAACA,YAAmB,KAAK,SAASA,SAAQ,MAAM;AAChE,QAAM,OAAO,CAACA,YAAmB,KAAK,WAAWA,SAAQ,MAAM;AAC/D,QAAM,QAAQ,CAACA,YACb,OAAO,aAAa,GAAG,IAAI,WAAW,MAAMA,SAAQ,CAAC,CAAC,EACnD,QAAQ,OAAO,EAAE,EACjB,KAAI;AACT,MAAI,MAAM,CAAC,MAAM,cAAc,QAAQ,EAAE,MAAM,MAAM,MAAM,EAAE,MAAM;AACjE,UAAM,IAAI,MAAM,wCAAwC;AAC1D,QAAM,QAAQ,QAAQ,EAAE,GACtB,SAAS,QAAQ,uBAAuB,QAAQ,IAAI;AACtD,MAAI,QAAQ,KAAK,QAAQ,KAAK,OAAO,KAAK,aAAa,QAAQ,MAAM,IAAI,OAAO;AAC9E,UAAM,IAAI,MAAM,4BAA4B;AAC9C,MAAI,SAAS;AACb,QAAM,WAAsB,CAAA;AAC5B,WAAS,IAAI,GAAG,IAAI,OAAO,KAAK;AAC9B,QAAI,SAAS,MAAM,KAAK;AAAY,YAAM,IAAI,MAAM,+BAA+B;AACnF,UAAM,WAAW,KAAK,SAAS,EAAE,GAC/B,WAAW,KAAK,SAAS,EAAE,GAC3B,WAAW,KAAK,SAAS,GAAG,GAC5B,WAAW,KAAK,SAAS,GAAG;AAC9B,UAAM,UAAU,KAAK,SAAS,GAAG,GAC/B,UAAU,KAAK,SAAS,GAAG,GAC3B,QAAQ,QAAQ,SAAS,GAAG;AAC9B,UAAM,QAAQ,KAAK,WAAW,YAAY,SACxC,SAAS,KAAK,WAAW,YAAY;AACvC,QACE,CAAC,CAAC,OAAO,MAAM,EAAE,MACf,OAAK,OAAO,SAAS,CAAC,KAAK,KAAK,KAAK,KAAK,IAAI,IAAI,KAAK,MAAM,CAAC,CAAC,IAAI,IAAI,KAEzE,EAAE,UAAU,KAAK,UAAU,MAC3B,UAAU,KAAK,MAAM,KAAK,IAAI,KAAK,MAAM,MAAM;AAE/C,YAAM,IAAI,MAAM,iCAAiC;AACnD,UAAM,QAAQ,SAAS,KACrB,MAAM,QAAQ,QAAQ;AACxB,QAAI,MAAM,KAAK;AAAY,YAAM,IAAI,MAAM,6BAA6B;AACxE,UAAM,SAAS,IAAI,aAAa,QAAQ,CAAC;AACzC,aAAS,IAAI,GAAG,IAAI,OAAO,KAAK;AAC9B,aAAO,IAAI,CAAC,IAAI,KAAK,WAAW,QAAQ,IAAI,SAAS,GAAG,MAAM,IAAI;AAClE,aAAO,IAAI,IAAI,CAAC,IAAI,KAAK,WAAW,QAAQ,IAAI,QAAQ,MAAM,IAAI;IACpE;AACA,aAAS,KAAK;MACZ,QAAQ,CAAC,WAAW,YAAY,WAAW,UAAU;MACrD,MAAM,CAAC,UAAU,YAAY,UAAU,UAAU;MACjD,MAAM,CAAC,KAAK,MAAM,KAAK,GAAG,KAAK,MAAM,MAAM,CAAC;MAC5C;KACD;AACD,aAAS;EACX;AACA,SAAO,gBAAgB,QAAQ;AACjC;",
+  "names": ["offset"]
+}

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

@@ -0,0 +1,8 @@
+import type { DatumGrid } from "./types.js";
+export type NTv2GridOptions = {
+    /** False selects compact 8-byte nodes; true (default) selects 16-byte nodes with accuracy fields. */
+    includeErrorFields?: boolean;
+};
+/** Decode an owned grid snapshot. The input buffer may be changed or released after this returns. */
+export declare function parseNTv2Grid(data: ArrayBuffer, options?: NTv2GridOptions): DatumGrid;
+//# sourceMappingURL=ntv2.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"ntv2.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/ntv2.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,SAAS,EAAU,mBAAgB;AAEhD,MAAM,MAAM,eAAe,GAAG;IAC5B,qGAAqG;IACrG,kBAAkB,CAAC,EAAE,OAAO,CAAC;CAC9B,CAAC;AAEF,qGAAqG;AACrG,wBAAgB,aAAa,CAAC,IAAI,EAAE,WAAW,EAAE,OAAO,GAAE,eAAoB,GAAG,SAAS,CAwDzF"}

package/dist/experimental/grids/ntv2.js

@@ -0,0 +1,56 @@
+// 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: NTv2 decoding directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../PROJ4-LICENSE.md. Modified for explicit ownership, bounds/geometry validation and no global registration.
+import { createDatumGrid } from "./grid.js";
+const ARC_SECOND = Math.PI / 648000;
+/** Decode an owned grid snapshot. The input buffer may be changed or released after this returns. */
+export function parseNTv2Grid(data, options = {}) {
+    const view = new DataView(data);
+    if (view.byteLength < 176)
+        throw new Error('Truncated NTv2 overview header');
+    const little = view.getInt32(8, true) === 11;
+    if (!little && view.getInt32(8, false) !== 11)
+        throw new Error('Invalid NTv2 header or byte order');
+    const integer = (offset) => view.getInt32(offset, little);
+    const real = (offset) => view.getFloat64(offset, little);
+    const label = (offset) => String.fromCharCode(...new Uint8Array(data, offset, 8))
+        .replace(/\0/g, '')
+        .trim();
+    if (label(0) !== 'NUM_OREC' || integer(24) !== 11 || label(56) !== 'SECONDS')
+        throw new Error('Unsupported NTv2 header or shift units');
+    const count = integer(40), stride = options.includeErrorFields === false ? 8 : 16;
+    if (count < 1 || count > Math.floor((view.byteLength - 176) / (176 + 4 * stride)))
+        throw new Error('Invalid NTv2 subgrid count');
+    let offset = 176;
+    const subgrids = [];
+    for (let n = 0; n < count; n++) {
+        if (offset + 176 > view.byteLength)
+            throw new Error('Truncated NTv2 subgrid header');
+        const lowerLat = real(offset + 72), upperLat = real(offset + 88), lowerLon = real(offset + 104), upperLon = real(offset + 120);
+        const latStep = real(offset + 136), lonStep = real(offset + 152), nodes = integer(offset + 168);
+        const width = 1 + (upperLon - lowerLon) / lonStep, height = 1 + (upperLat - lowerLat) / latStep;
+        if (![width, height].every(v => Number.isFinite(v) && v >= 2 && Math.abs(v - Math.round(v)) < 1e-7) ||
+            !(lonStep > 0 && latStep > 0) ||
+            nodes !== Math.round(width) * Math.round(height))
+            throw new Error('Invalid NTv2 subgrid dimensions');
+        const start = offset + 176, end = start + nodes * stride;
+        if (end > view.byteLength)
+            throw new Error('Truncated NTv2 node records');
+        const shifts = new Float64Array(nodes * 2);
+        for (let i = 0; i < nodes; i++) {
+            shifts[i * 2] = view.getFloat32(start + i * stride + 4, little) * ARC_SECOND;
+            shifts[i * 2 + 1] = view.getFloat32(start + i * stride, little) * ARC_SECOND;
+        }
+        subgrids.push({
+            origin: [lowerLon * ARC_SECOND, lowerLat * ARC_SECOND],
+            step: [lonStep * ARC_SECOND, latStep * ARC_SECOND],
+            size: [Math.round(width), Math.round(height)],
+            shifts
+        });
+        offset = end;
+    }
+    return createDatumGrid(subgrids);
+}
+//# sourceMappingURL=ntv2.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"ntv2.js","sourceRoot":"","sources":["../../../src/experimental/grids/ntv2.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,8MAA8M;AAC9M,OAAO,EAAC,eAAe,EAAC,kBAAe;AAOvC,MAAM,UAAU,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC;AACpC,qGAAqG;AACrG,MAAM,UAAU,aAAa,CAAC,IAAiB,EAAE,UAA2B,EAAE;IAC5E,MAAM,IAAI,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,CAAC;IAChC,IAAI,IAAI,CAAC,UAAU,GAAG,GAAG;QAAE,MAAM,IAAI,KAAK,CAAC,gCAAgC,CAAC,CAAC;IAC7E,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC;IAC7C,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,KAAK,EAAE;QAC3C,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;IACvD,MAAM,OAAO,GAAG,CAAC,MAAc,EAAE,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAClE,MAAM,IAAI,GAAG,CAAC,MAAc,EAAE,EAAE,CAAC,IAAI,CAAC,UAAU,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACjE,MAAM,KAAK,GAAG,CAAC,MAAc,EAAE,EAAE,CAC/B,MAAM,CAAC,YAAY,CAAC,GAAG,IAAI,UAAU,CAAC,IAAI,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;SACpD,OAAO,CAAC,KAAK,EAAE,EAAE,CAAC;SAClB,IAAI,EAAE,CAAC;IACZ,IAAI,KAAK,CAAC,CAAC,CAAC,KAAK,UAAU,IAAI,OAAO,CAAC,EAAE,CAAC,KAAK,EAAE,IAAI,KAAK,CAAC,EAAE,CAAC,KAAK,SAAS;QAC1E,MAAM,IAAI,KAAK,CAAC,wCAAwC,CAAC,CAAC;IAC5D,MAAM,KAAK,GAAG,OAAO,CAAC,EAAE,CAAC,EACvB,MAAM,GAAG,OAAO,CAAC,kBAAkB,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IACzD,IAAI,KAAK,GAAG,CAAC,IAAI,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;QAC/E,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;IAChD,IAAI,MAAM,GAAG,GAAG,CAAC;IACjB,MAAM,QAAQ,GAAc,EAAE,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,UAAU;YAAE,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;QACrF,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,GAAG,EAAE,CAAC,EAChC,QAAQ,GAAG,IAAI,CAAC,MAAM,GAAG,EAAE,CAAC,EAC5B,QAAQ,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC,EAC7B,QAAQ,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC,EAChC,OAAO,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC,EAC5B,KAAK,GAAG,OAAO,CAAC,MAAM,GAAG,GAAG,CAAC,CAAC;QAChC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,OAAO,EAC/C,MAAM,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,OAAO,CAAC;QAC/C,IACE,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,KAAK,CACpB,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CACxE;YACD,CAAC,CAAC,OAAO,GAAG,CAAC,IAAI,OAAO,GAAG,CAAC,CAAC;YAC7B,KAAK,KAAK,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;YAEhD,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;QACrD,MAAM,KAAK,GAAG,MAAM,GAAG,GAAG,EACxB,GAAG,GAAG,KAAK,GAAG,KAAK,GAAG,MAAM,CAAC;QAC/B,IAAI,GAAG,GAAG,IAAI,CAAC,UAAU;YAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,CAAC,CAAC;QAC1E,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,EAAE,MAAM,CAAC,GAAG,UAAU,CAAC;YAC7E,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,GAAG,CAAC,GAAG,MAAM,EAAE,MAAM,CAAC,GAAG,UAAU,CAAC;QAC/E,CAAC;QACD,QAAQ,CAAC,IAAI,CAAC;YACZ,MAAM,EAAE,CAAC,QAAQ,GAAG,UAAU,EAAE,QAAQ,GAAG,UAAU,CAAC;YACtD,IAAI,EAAE,CAAC,OAAO,GAAG,UAAU,EAAE,OAAO,GAAG,UAAU,CAAC;YAClD,IAAI,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;YAC7C,MAAM;SACP,CAAC,CAAC;QACH,MAAM,GAAG,GAAG,CAAC;IACf,CAAC;IACD,OAAO,eAAe,CAAC,QAAQ,CAAC,CAAC;AACnC,CAAC"}

package/dist/experimental/grids/scalar-geotiff.d.ts

@@ -0,0 +1,58 @@
+import type { VerticalGrid } from "./types.js";
+type Metadata = Record<string, unknown> | null;
+type DirectoryTag = 'ModelPixelScale' | 'ModelTiepoint' | 'ModelTransformation' | 'NewSubfileType';
+/** Structural subset of geotiff.js v2/v3. Metadata may be synchronous or asynchronous. */
+export type VerticalGridGeoTIFFImage = {
+    getWidth(): number;
+    getHeight(): number;
+    getGeoKeys(): Record<string, unknown> | null;
+    getGDALMetadata(sample?: number | null): Metadata | Promise<Metadata>;
+    getGDALNoData?(): number | null;
+    fileDirectory: Partial<Record<DirectoryTag, unknown>> | {
+        getValue(name: DirectoryTag): unknown;
+    };
+    readRasters(options: {
+        samples: number[];
+        interleave: false;
+    }): Promise<ArrayLike<number | ArrayLike<number>>>;
+};
+export type VerticalGridGeoTIFF = {
+    getImageCount(): Promise<number>;
+    getImage(index: number): Promise<VerticalGridGeoTIFFImage>;
+};
+/** Plain decoded raster input, structurally compatible with loaders.gl GeoTIFFRasterData. */
+export type VerticalGridGeoTIFFData = {
+    /** Images in original file order, with parents before nested children. */
+    images: readonly {
+        /** Original dimensions in pixels. */
+        width: number;
+        /** Original dimensions in pixels. */
+        height: number;
+        /** Raw bands; original band zero must be present. */
+        bands: readonly {
+            /** Original zero-based sample index, retained when selecting bands. */
+            index: number;
+            /** Native, unscaled samples in TIFF row order. */
+            data: ArrayLike<number>;
+            /** Per-band GDAL metadata. */
+            metadata: Metadata;
+        }[];
+        /** Decoded CRS and pixel-registration GeoKeys. */
+        geoKeys: Metadata;
+        /** Image-level GDAL metadata. */
+        metadata: Metadata;
+        /** Declared nodata value, before scale/offset. */
+        noData: number | null;
+        /** Original geometry and image-kind tags. */
+        fileDirectory: Partial<Record<DirectoryTag, unknown>>;
+    }[];
+};
+/** Prepare metre geoid undulations from geographic, north-up PROJ Geodetic TIFF grids.
+ * The first band must declare geoid_undulation; dataset TYPE must be
+ * VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL. PixelIsPoint and PixelIsArea are supported.
+ * Later nested images take precedence. Overviews, ambiguous overlaps, rotations,
+ * non-metre bands and other operation types are rejected rather than inferred.
+ */
+export declare function prepareScalarGeoTIFFGrids(tiff: VerticalGridGeoTIFF | VerticalGridGeoTIFFData): Promise<VerticalGrid[]>;
+export {};
+//# sourceMappingURL=scalar-geotiff.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"scalar-geotiff.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/scalar-geotiff.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,YAAY,EAAC,mBAAgB;AAE1C,KAAK,QAAQ,GAAG,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC;AAC/C,KAAK,YAAY,GAAG,iBAAiB,GAAG,eAAe,GAAG,qBAAqB,GAAG,gBAAgB,CAAC;AACnG,0FAA0F;AAC1F,MAAM,MAAM,wBAAwB,GAAG;IACrC,QAAQ,IAAI,MAAM,CAAC;IACnB,SAAS,IAAI,MAAM,CAAC;IACpB,UAAU,IAAI,MAAM,CAAC,MAAM,EAAE,OAAO,CAAC,GAAG,IAAI,CAAC;IAC7C,eAAe,CAAC,MAAM,CAAC,EAAE,MAAM,GAAG,IAAI,GAAG,QAAQ,GAAG,OAAO,CAAC,QAAQ,CAAC,CAAC;IACtE,aAAa,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC;IAChC,aAAa,EAAE,OAAO,CAAC,MAAM,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC,GAAG;QAAC,QAAQ,CAAC,IAAI,EAAE,YAAY,GAAG,OAAO,CAAA;KAAC,CAAC;IAChG,WAAW,CAAC,OAAO,EAAE;QACnB,OAAO,EAAE,MAAM,EAAE,CAAC;QAClB,UAAU,EAAE,KAAK,CAAC;KACnB,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;CACpD,CAAC;AACF,MAAM,MAAM,mBAAmB,GAAG;IAChC,aAAa,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,GAAG,OAAO,CAAC,wBAAwB,CAAC,CAAC;CAC5D,CAAC;AAEF,6FAA6F;AAC7F,MAAM,MAAM,uBAAuB,GAAG;IACpC,0EAA0E;IAC1E,MAAM,EAAE,SAAS;QACf,qCAAqC;QACrC,KAAK,EAAE,MAAM,CAAC;QACd,qCAAqC;QACrC,MAAM,EAAE,MAAM,CAAC;QACf,qDAAqD;QACrD,KAAK,EAAE,SAAS;YACd,uEAAuE;YACvE,KAAK,EAAE,MAAM,CAAC;YACd,kDAAkD;YAClD,IAAI,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;YACxB,8BAA8B;YAC9B,QAAQ,EAAE,QAAQ,CAAC;SACpB,EAAE,CAAC;QACJ,kDAAkD;QAClD,OAAO,EAAE,QAAQ,CAAC;QAClB,iCAAiC;QACjC,QAAQ,EAAE,QAAQ,CAAC;QACnB,kDAAkD;QAClD,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;QACtB,6CAA6C;QAC7C,aAAa,EAAE,OAAO,CAAC,MAAM,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC,CAAC;KACvD,EAAE,CAAC;CACL,CAAC;AAoCF;;;;;GAKG;AACH,wBAAsB,yBAAyB,CAC7C,IAAI,EAAE,mBAAmB,GAAG,uBAAuB,GAClD,OAAO,CAAC,YAAY,EAAE,CAAC,CAqIzB"}

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

@@ -0,0 +1,147 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adapter for PROJ's Geodetic TIFF Grid specification, not a proj4js/PROJ fork. Decoding is caller-owned; no TIFF library is imported at runtime.
+import { createVerticalGrid } from "./vertical.js";
+/** Adapts plain numeric rasters without importing or retaining a TIFF decoder. */
+function getRasterImage(data, index) {
+    const image = data.images[index];
+    const band = image.bands.find(candidate => candidate.index === 0);
+    if (!band)
+        throw new Error('Vertical GeoTIFF requires original band zero');
+    return {
+        getWidth: () => image.width,
+        getHeight: () => image.height,
+        getGeoKeys: () => image.geoKeys,
+        getGDALMetadata: sample => (sample === 0 ? band.metadata : image.metadata),
+        getGDALNoData: () => image.noData,
+        fileDirectory: image.fileDirectory,
+        readRasters: async () => [band.data]
+    };
+}
+function number(value, fallback) {
+    if (value === undefined)
+        return fallback;
+    if ((typeof value !== 'number' && typeof value !== 'string') ||
+        (typeof value === 'string' && !value.trim()) ||
+        !Number.isFinite(Number(value)))
+        throw new Error('Invalid vertical GeoTIFF scale or offset');
+    return Number(value);
+}
+function vector(value, length) {
+    if (!value || typeof value !== 'object' || !('length' in value) || value.length !== length)
+        throw new Error('Vertical GeoTIFF requires one tiepoint and pixel scale');
+    const result = Array.from(value);
+    if (!result.every(Number.isFinite))
+        throw new Error('Invalid vertical GeoTIFF geometry');
+    return result;
+}
+/** Prepare metre geoid undulations from geographic, north-up PROJ Geodetic TIFF grids.
+ * The first band must declare geoid_undulation; dataset TYPE must be
+ * VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL. PixelIsPoint and PixelIsArea are supported.
+ * Later nested images take precedence. Overviews, ambiguous overlaps, rotations,
+ * non-metre bands and other operation types are rejected rather than inferred.
+ */
+export async function prepareScalarGeoTIFFGrids(tiff) {
+    const count = 'images' in tiff ? tiff.images.length : await tiff.getImageCount();
+    if (!Number.isSafeInteger(count) || count < 1)
+        throw new Error('Vertical GeoTIFF requires at least one image');
+    const prepared = [];
+    for (let n = 0; n < count; n++) {
+        const image = 'images' in tiff ? getRasterImage(tiff, n) : await tiff.getImage(n);
+        const directory = image.fileDirectory;
+        const tag = (name) => 'getValue' in directory ? directory.getValue(name) : directory[name];
+        if (tag('NewSubfileType') !== undefined && tag('NewSubfileType') !== 0)
+            throw new Error('Vertical GeoTIFF overviews and masks are unsupported');
+        if (tag('ModelTransformation') !== undefined)
+            throw new Error('Vertical GeoTIFF requires axis-aligned tiepoint/scale geometry');
+        const keys = image.getGeoKeys();
+        if (keys?.['GTModelTypeGeoKey'] !== 2 ||
+            (keys['GeogAngularUnitsGeoKey'] !== undefined && keys['GeogAngularUnitsGeoKey'] !== 9102))
+            throw new Error('Vertical GeoTIFF requires a geographic degree grid');
+        if ((keys['GeogPrimeMeridianGeoKey'] !== undefined && keys['GeogPrimeMeridianGeoKey'] !== 8901) ||
+            (keys['GeogPrimeMeridianLongGeoKey'] !== undefined &&
+                keys['GeogPrimeMeridianLongGeoKey'] !== 0))
+            throw new Error('Vertical GeoTIFF requires Greenwich longitudes');
+        const rasterType = keys['GTRasterTypeGeoKey'];
+        if (rasterType !== 1 && rasterType !== 2)
+            throw new Error('Vertical GeoTIFF requires explicit PixelIsPoint or PixelIsArea');
+        const metadata = await image.getGDALMetadata();
+        const band = await image.getGDALMetadata(0);
+        if (metadata?.['TYPE'] !== 'VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL' ||
+            band?.['DESCRIPTION'] !== 'geoid_undulation')
+            throw new Error('Vertical GeoTIFF requires a geographic-to-vertical geoid_undulation band');
+        if (band['UNITTYPE'] !== undefined && band['UNITTYPE'] !== 'metre')
+            throw new Error('Vertical GeoTIFF geoid undulations must be in metres');
+        if (metadata['interpolation_method'] !== undefined &&
+            metadata['interpolation_method'] !== 'bilinear')
+            throw new Error('Vertical GeoTIFF supports only bilinear interpolation');
+        const scale = number(band['SCALE'], 1), offset = number(band['OFFSET'], 0);
+        const [dx, dy] = vector(tag('ModelPixelScale'), 3);
+        const [i, j, k, x, y, z] = vector(tag('ModelTiepoint'), 6);
+        const width = image.getWidth(), height = image.getHeight();
+        if (!(dx > 0 && dy > 0) ||
+            k !== 0 ||
+            z !== 0 ||
+            ![width, height].every(v => Number.isSafeInteger(v) && v >= 2))
+            throw new Error('Invalid vertical GeoTIFF geometry or dimensions');
+        const half = rasterType === 1 ? 0.5 : 0;
+        const west = x + (half - i) * dx;
+        const north = y - (half - j) * dy;
+        const south = north - (height - 1) * dy, east = west + (width - 1) * dx;
+        if (![west, south, east, north].every(Number.isFinite) ||
+            !(east > west && north > south) ||
+            east - west > 360 ||
+            south < -90 ||
+            north > 90 ||
+            !Number.isSafeInteger(width * height))
+            throw new Error('Invalid vertical GeoTIFF geographic extent');
+        // Require ordered nested or disjoint images, so TIFF ordering cannot silently
+        // replace a fine child with a coarse parent or an ambiguous overlapping grid.
+        for (const other of prepared) {
+            const overlaps = west < other.east && east > other.west && south < other.north && north > other.south;
+            if (overlaps &&
+                !(west >= other.west &&
+                    east <= other.east &&
+                    south >= other.south &&
+                    north <= other.north &&
+                    dx <= other.dx &&
+                    dy <= other.dy))
+                throw new Error('Vertical GeoTIFF images must be ordered parent before nested child or disjoint');
+        }
+        const rasters = await image.readRasters({ samples: [0], interleave: false });
+        const values = rasters[0];
+        if (!values || typeof values === 'number' || values.length !== width * height)
+            throw new Error('Vertical GeoTIFF requires a complete non-interleaved height band');
+        const declaredNoData = image.getGDALNoData?.();
+        // GDAL_NODATA is text; compare it in the decoded band's numeric precision.
+        const nodata = values instanceof Float32Array && typeof declaredNoData === 'number'
+            ? Math.fround(declaredNoData)
+            : declaredNoData;
+        const offsets = new Float64Array(width * height);
+        for (let row = 0; row < height; row++)
+            for (let col = 0; col < width; col++) {
+                const raw = values[row * width + col];
+                const value = raw === nodata || !Number.isFinite(raw) ? NaN : raw * scale + offset;
+                if (Number.isFinite(raw) && raw !== nodata && !Number.isFinite(value))
+                    throw new Error('Vertical GeoTIFF scale/offset produced a non-finite height');
+                offsets[(height - 1 - row) * width + col] = value;
+            }
+        prepared.push({
+            grid: createVerticalGrid({
+                origin: [west, south],
+                step: [dx, dy],
+                size: [width, height],
+                offsets
+            }),
+            west,
+            south,
+            east,
+            north,
+            dx,
+            dy
+        });
+    }
+    return prepared.map(image => image.grid);
+}
+//# sourceMappingURL=scalar-geotiff.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"scalar-geotiff.js","sourceRoot":"","sources":["../../../src/experimental/grids/scalar-geotiff.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,6KAA6K;AAC7K,OAAO,EAAC,kBAAkB,EAAC,sBAAmB;AAmD9C,kFAAkF;AAClF,SAAS,cAAc,CAAC,IAA6B,EAAE,KAAa;IAClE,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;IACjC,MAAM,IAAI,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,EAAE,CAAC,SAAS,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC;IAClE,IAAI,CAAC,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAC3E,OAAO;QACL,QAAQ,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,KAAK;QAC3B,SAAS,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,MAAM;QAC7B,UAAU,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,OAAO;QAC/B,eAAe,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,QAAQ,CAAC;QAC1E,aAAa,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,MAAM;QACjC,aAAa,EAAE,KAAK,CAAC,aAAa;QAClC,WAAW,EAAE,KAAK,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;KACrC,CAAC;AACJ,CAAC;AAED,SAAS,MAAM,CAAC,KAAc,EAAE,QAAgB;IAC9C,IAAI,KAAK,KAAK,SAAS;QAAE,OAAO,QAAQ,CAAC;IACzC,IACE,CAAC,OAAO,KAAK,KAAK,QAAQ,IAAI,OAAO,KAAK,KAAK,QAAQ,CAAC;QACxD,CAAC,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;QAC5C,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;QAE/B,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,OAAO,MAAM,CAAC,KAAK,CAAC,CAAC;AACvB,CAAC;AACD,SAAS,MAAM,CAAC,KAAc,EAAE,MAAc;IAC5C,IAAI,CAAC,KAAK,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,CAAC,QAAQ,IAAI,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,MAAM;QACxF,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,KAA0B,CAAC,CAAC;IACtD,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;IACzF,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;;GAKG;AACH,MAAM,CAAC,KAAK,UAAU,yBAAyB,CAC7C,IAAmD;IAEnD,MAAM,KAAK,GAAG,QAAQ,IAAI,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,IAAI,CAAC,aAAa,EAAE,CAAC;IACjF,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC;QAC3C,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;IAClE,MAAM,QAAQ,GAQR,EAAE,CAAC;IACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,KAAK,GAAG,QAAQ,IAAI,IAAI,CAAC,CAAC,CAAC,cAAc,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAClF,MAAM,SAAS,GAAG,KAAK,CAAC,aAAa,CAAC;QACtC,MAAM,GAAG,GAAG,CAAC,IAAkB,EAAE,EAAE,CACjC,UAAU,IAAI,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC;QACvE,IAAI,GAAG,CAAC,gBAAgB,CAAC,KAAK,SAAS,IAAI,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACpE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,IAAI,GAAG,CAAC,qBAAqB,CAAC,KAAK,SAAS;YAC1C,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;QACpF,MAAM,IAAI,GAAG,KAAK,CAAC,UAAU,EAAE,CAAC;QAChC,IACE,IAAI,EAAE,CAAC,mBAAmB,CAAC,KAAK,CAAC;YACjC,CAAC,IAAI,CAAC,wBAAwB,CAAC,KAAK,SAAS,IAAI,IAAI,CAAC,wBAAwB,CAAC,KAAK,IAAI,CAAC;YAEzF,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;QACxE,IACE,CAAC,IAAI,CAAC,yBAAyB,CAAC,KAAK,SAAS,IAAI,IAAI,CAAC,yBAAyB,CAAC,KAAK,IAAI,CAAC;YAC3F,CAAC,IAAI,CAAC,6BAA6B,CAAC,KAAK,SAAS;gBAChD,IAAI,CAAC,6BAA6B,CAAC,KAAK,CAAC,CAAC;YAE5C,MAAM,IAAI,KAAK,CAAC,gDAAgD,CAAC,CAAC;QACpE,MAAM,UAAU,GAAG,IAAI,CAAC,oBAAoB,CAAC,CAAC;QAC9C,IAAI,UAAU,KAAK,CAAC,IAAI,UAAU,KAAK,CAAC;YACtC,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;QACpF,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,eAAe,EAAE,CAAC;QAC/C,MAAM,IAAI,GAAG,MAAM,KAAK,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;QAC5C,IACE,QAAQ,EAAE,CAAC,MAAM,CAAC,KAAK,wCAAwC;YAC/D,IAAI,EAAE,CAAC,aAAa,CAAC,KAAK,kBAAkB;YAE5C,MAAM,IAAI,KAAK,CAAC,0EAA0E,CAAC,CAAC;QAC9F,IAAI,IAAI,CAAC,UAAU,CAAC,KAAK,SAAS,IAAI,IAAI,CAAC,UAAU,CAAC,KAAK,OAAO;YAChE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,IACE,QAAQ,CAAC,sBAAsB,CAAC,KAAK,SAAS;YAC9C,QAAQ,CAAC,sBAAsB,CAAC,KAAK,UAAU;YAE/C,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,EACpC,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;QACrC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,iBAAiB,CAAC,EAAE,CAAC,CAAC,CAAC;QACnD,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,eAAe,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3D,MAAM,KAAK,GAAG,KAAK,CAAC,QAAQ,EAAE,EAC5B,MAAM,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;QAC7B,IACE,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;YACnB,CAAC,KAAK,CAAC;YACP,CAAC,KAAK,CAAC;YACP,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAE9D,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;QACrE,MAAM,IAAI,GAAG,UAAU,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACjC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QAClC,MAAM,KAAK,GAAG,KAAK,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,EACrC,IAAI,GAAG,IAAI,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACjC,IACE,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAClD,CAAC,CAAC,IAAI,GAAG,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;YAC/B,IAAI,GAAG,IAAI,GAAG,GAAG;YACjB,KAAK,GAAG,CAAC,EAAE;YACX,KAAK,GAAG,EAAE;YACV,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,GAAG,MAAM,CAAC;YAErC,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;QAChE,8EAA8E;QAC9E,8EAA8E;QAC9E,KAAK,MAAM,KAAK,IAAI,QAAQ,EAAE,CAAC;YAC7B,MAAM,QAAQ,GACZ,IAAI,GAAG,KAAK,CAAC,IAAI,IAAI,IAAI,GAAG,KAAK,CAAC,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC,KAAK,IAAI,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;YACvF,IACE,QAAQ;gBACR,CAAC,CACC,IAAI,IAAI,KAAK,CAAC,IAAI;oBAClB,IAAI,IAAI,KAAK,CAAC,IAAI;oBAClB,KAAK,IAAI,KAAK,CAAC,KAAK;oBACpB,KAAK,IAAI,KAAK,CAAC,KAAK;oBACpB,EAAE,IAAI,KAAK,CAAC,EAAE;oBACd,EAAE,IAAI,KAAK,CAAC,EAAE,CACf;gBAED,MAAM,IAAI,KAAK,CACb,gFAAgF,CACjF,CAAC;QACN,CAAC;QACD,MAAM,OAAO,GAAG,MAAM,KAAK,CAAC,WAAW,CAAC,EAAC,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,UAAU,EAAE,KAAK,EAAC,CAAC,CAAC;QAC3E,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QAC1B,IAAI,CAAC,MAAM,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,MAAM,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM;YAC3E,MAAM,IAAI,KAAK,CAAC,kEAAkE,CAAC,CAAC;QACtF,MAAM,cAAc,GAAG,KAAK,CAAC,aAAa,EAAE,EAAE,CAAC;QAC/C,2EAA2E;QAC3E,MAAM,MAAM,GACV,MAAM,YAAY,YAAY,IAAI,OAAO,cAAc,KAAK,QAAQ;YAClE,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC;YAC7B,CAAC,CAAC,cAAc,CAAC;QACrB,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;QACjD,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,GAAG,EAAE;YACnC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC;gBACrC,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,GAAG,KAAK,GAAG,GAAG,CAAC,CAAC;gBACtC,MAAM,KAAK,GAAG,GAAG,KAAK,MAAM,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,KAAK,GAAG,MAAM,CAAC;gBACnF,IAAI,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,MAAM,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;oBACnE,MAAM,IAAI,KAAK,CAAC,4DAA4D,CAAC,CAAC;gBAChF,OAAO,CAAC,CAAC,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC;YACpD,CAAC;QACH,QAAQ,CAAC,IAAI,CAAC;YACZ,IAAI,EAAE,kBAAkB,CAAC;gBACvB,MAAM,EAAE,CAAC,IAAI,EAAE,KAAK,CAAC;gBACrB,IAAI,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;gBACd,IAAI,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC;gBACrB,OAAO;aACR,CAAC;YACF,IAAI;YACJ,KAAK;YACL,IAAI;YACJ,KAAK;YACL,EAAE;YACF,EAAE;SACH,CAAC,CAAC;IACL,CAAC;IACD,OAAO,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;AAC3C,CAAC"}

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

@@ -0,0 +1,35 @@
+/** Prepared horizontal datum grid. Inputs/outputs are east-positive longitude/latitude radians. */
+export type DatumGrid = {
+    readonly subgridCount: number;
+    /** Optional mutable hook. On false, x/y must remain unchanged; always preserve z. */
+    shiftInPlace?(point: {
+        x: number;
+        y: number;
+    }, inverse: boolean): boolean;
+    /** undefined means no usable subgrid covers this coordinate; other failures throw. */
+    shift(longitude: number, latitude: number, inverse: boolean): [number, number] | undefined;
+};
+/** Per-instance registration. Readers produce grids before synchronous construction. */
+export type DatumGridCollection = Readonly<Record<string, DatumGrid>>;
+export type DatumGridReference = {
+    readonly name: string;
+    readonly optional: boolean;
+    readonly grid?: DatumGrid;
+};
+/** Internal west-positive node layout, matching the pinned proj4js grid convention. */
+export type Subgrid = {
+    readonly origin: readonly [number, number];
+    readonly step: readonly [number, number];
+    readonly size: readonly [number, number];
+    /** Reader-owned west-positive longitude/latitude shift pairs in radians. */
+    readonly shifts: Float64Array;
+};
+/** Geoid undulation N in metres: ellipsoidal height h = gravity-related height H + N.
+ * Inputs are east-positive Greenwich longitude and latitude in radians, in the
+ * horizontal datum of the CRS carrying +geoidgrids. Undefined means no coverage.
+ */
+export type VerticalGrid = {
+    getOffset(longitude: number, latitude: number): number | undefined;
+};
+export type VerticalGridCollection = Readonly<Record<string, VerticalGrid>>;
+//# sourceMappingURL=types.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"types.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/types.ts"],"names":[],"mappings":"AAGA,mGAAmG;AACnG,MAAM,MAAM,SAAS,GAAG;IACtB,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,qFAAqF;IACrF,YAAY,CAAC,CAAC,KAAK,EAAE;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAA;KAAC,EAAE,OAAO,EAAE,OAAO,GAAG,OAAO,CAAC;IACxE,sFAAsF;IACtF,KAAK,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,SAAS,CAAC;CAC5F,CAAC;AACF,wFAAwF;AACxF,MAAM,MAAM,mBAAmB,GAAG,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,CAAC;AACtE,MAAM,MAAM,kBAAkB,GAAG;IAC/B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,QAAQ,EAAE,OAAO,CAAC;IAC3B,QAAQ,CAAC,IAAI,CAAC,EAAE,SAAS,CAAC;CAC3B,CAAC;AACF,uFAAuF;AACvF,MAAM,MAAM,OAAO,GAAG;IACpB,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3C,QAAQ,CAAC,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACzC,QAAQ,CAAC,IAAI,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACzC,4EAA4E;IAC5E,QAAQ,CAAC,MAAM,EAAE,YAAY,CAAC;CAC/B,CAAC;AAEF;;;GAGG;AACH,MAAM,MAAM,YAAY,GAAG;IACzB,SAAS,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,MAAM,GAAG,SAAS,CAAC;CACpE,CAAC;AACF,MAAM,MAAM,sBAAsB,GAAG,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,YAAY,CAAC,CAAC,CAAC"}

package/dist/experimental/grids/types.js

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=types.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"types.js","sourceRoot":"","sources":["../../../src/experimental/grids/types.ts"],"names":[],"mappings":""}