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

package/dist/experimental/grids/geotiff.cjs

@@ -0,0 +1,180 @@
+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/geotiff.js
+var geotiff_exports = {};
+__export(geotiff_exports, {
+  loadGeoTIFFGrid: () => loadGeoTIFFGrid
+});
+module.exports = __toCommonJS(geotiff_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/geotiff.js
+async function loadGeoTIFFGrid(tiff) {
+  var _a;
+  const count = await tiff.getImageCount();
+  if (!Number.isSafeInteger(count) || count < 1)
+    throw new Error("GeoTIFF requires at least one image");
+  const subgrids = [];
+  for (let n = count - 1; n >= 0; n--) {
+    const image = await tiff.getImage(n), width = image.getWidth(), height = image.getHeight();
+    if (![width, height].every((v) => Number.isSafeInteger(v) && v >= 2))
+      throw new Error("Invalid GeoTIFF grid dimensions");
+    const bbox = image.getBoundingBox();
+    const directory = image.fileDirectory;
+    const scale = "getValue" in directory ? directory.getValue("ModelPixelScale") : directory.ModelPixelScale;
+    if (!scale || typeof scale !== "object" || !("0" in scale) || !("1" in scale))
+      throw new Error("GeoTIFF requires ModelPixelScale");
+    const sx = Number(scale[0]), sy = Number(scale[1]);
+    if (bbox.length !== 4 || !bbox.every(Number.isFinite) || !(sx > 0 && sy > 0 && Number.isFinite(sx) && Number.isFinite(sy)))
+      throw new Error("Invalid GeoTIFF geographic grid geometry");
+    const rasters = await image.readRasters(), latitude = rasters[0], longitude = rasters[1];
+    if (!latitude || !longitude || typeof latitude === "number" || typeof longitude === "number" || latitude.length !== width * height || longitude.length !== width * height)
+      throw new Error("GeoTIFF requires latitude and longitude offset bands");
+    const nodata = (_a = image.getGDALNoData) == null ? void 0 : _a.call(image);
+    const shifts = new Float64Array(width * height * 2);
+    let node = 0;
+    for (let row = height - 1; row >= 0; row--)
+      for (let col = width - 1; col >= 0; col--) {
+        const index = row * width + col, lat = latitude[index], lon = longitude[index];
+        const missing = lat === nodata || lon === nodata || !Number.isFinite(lat) || !Number.isFinite(lon);
+        shifts[node++] = missing ? NaN : -lon * Math.PI / 648e3;
+        shifts[node++] = missing ? NaN : lat * Math.PI / 648e3;
+      }
+    subgrids.push({
+      origin: [
+        -(bbox[0] + (width - 1) * sx) * Math.PI / 180,
+        (bbox[3] - (height - 1) * sy) * Math.PI / 180
+      ],
+      step: [sx * Math.PI / 180, sy * Math.PI / 180],
+      size: [width, height],
+      shifts
+    });
+  }
+  return createDatumGrid(subgrids);
+}
+//# sourceMappingURL=geotiff.cjs.map

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

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../src/experimental/grids/geotiff.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;;;AFlHA,eAAsB,gBAAgB,MAAsB;AA1B5D;AA2BE,QAAM,QAAQ,MAAM,KAAK,cAAa;AACtC,MAAI,CAAC,OAAO,cAAc,KAAK,KAAK,QAAQ;AAC1C,UAAM,IAAI,MAAM,qCAAqC;AACvD,QAAM,WAAsB,CAAA;AAC5B,WAAS,IAAI,QAAQ,GAAG,KAAK,GAAG,KAAK;AACnC,UAAM,QAAQ,MAAM,KAAK,SAAS,CAAC,GACjC,QAAQ,MAAM,SAAQ,GACtB,SAAS,MAAM,UAAS;AAC1B,QAAI,CAAC,CAAC,OAAO,MAAM,EAAE,MAAM,OAAK,OAAO,cAAc,CAAC,KAAK,KAAK,CAAC;AAC/D,YAAM,IAAI,MAAM,iCAAiC;AACnD,UAAM,OAAO,MAAM,eAAc;AACjC,UAAM,YAAY,MAAM;AACxB,UAAM,QACJ,cAAc,YAAY,UAAU,SAAS,iBAAiB,IAAI,UAAU;AAC9E,QAAI,CAAC,SAAS,OAAO,UAAU,YAAY,EAAE,OAAO,UAAU,EAAE,OAAO;AACrE,YAAM,IAAI,MAAM,kCAAkC;AACpD,UAAM,KAAK,OAAO,MAAM,CAAC,CAAC,GACxB,KAAK,OAAO,MAAM,CAAC,CAAC;AACtB,QACE,KAAK,WAAW,KAChB,CAAC,KAAK,MAAM,OAAO,QAAQ,KAC3B,EAAE,KAAK,KAAK,KAAK,KAAK,OAAO,SAAS,EAAE,KAAK,OAAO,SAAS,EAAE;AAE/D,YAAM,IAAI,MAAM,0CAA0C;AAC5D,UAAM,UAAU,MAAM,MAAM,YAAW,GACrC,WAAW,QAAQ,CAAC,GACpB,YAAY,QAAQ,CAAC;AACvB,QACE,CAAC,YACD,CAAC,aACD,OAAO,aAAa,YACpB,OAAO,cAAc,YACrB,SAAS,WAAW,QAAQ,UAC5B,UAAU,WAAW,QAAQ;AAE7B,YAAM,IAAI,MAAM,sDAAsD;AACxE,UAAM,UAAS,WAAM,kBAAN;AACf,UAAM,SAAS,IAAI,aAAa,QAAQ,SAAS,CAAC;AAClD,QAAI,OAAO;AACX,aAAS,MAAM,SAAS,GAAG,OAAO,GAAG;AACnC,eAAS,MAAM,QAAQ,GAAG,OAAO,GAAG,OAAO;AACzC,cAAM,QAAQ,MAAM,QAAQ,KAC1B,MAAM,SAAS,KAAK,GACpB,MAAM,UAAU,KAAK;AACvB,cAAM,UACJ,QAAQ,UAAU,QAAQ,UAAU,CAAC,OAAO,SAAS,GAAG,KAAK,CAAC,OAAO,SAAS,GAAG;AACnF,eAAO,MAAM,IAAI,UAAU,MAAO,CAAC,MAAM,KAAK,KAAM;AACpD,eAAO,MAAM,IAAI,UAAU,MAAO,MAAM,KAAK,KAAM;MACrD;AACF,aAAS,KAAK;MACZ,QAAQ;QACL,EAAE,KAAK,CAAC,KAAK,QAAQ,KAAK,MAAM,KAAK,KAAM;SAC1C,KAAK,CAAC,KAAK,SAAS,KAAK,MAAM,KAAK,KAAM;;MAE9C,MAAM,CAAE,KAAK,KAAK,KAAM,KAAM,KAAK,KAAK,KAAM,GAAG;MACjD,MAAM,CAAC,OAAO,MAAM;MACpB;KACD;EACH;AACA,SAAO,gBAAgB,QAAQ;AACjC;",
+  "names": []
+}

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

@@ -0,0 +1,23 @@
+import type { DatumGrid } from "./types.js";
+/** Structural subset supported by geotiff.js v2/v3; callers own file/network decoding. */
+export type DatumGridGeoTIFFImage = {
+    getWidth(): number;
+    getHeight(): number;
+    getBoundingBox(): number[];
+    fileDirectory: {
+        ModelPixelScale?: ArrayLike<number>;
+    } | {
+        getValue(name: 'ModelPixelScale'): unknown;
+    };
+    readRasters(): Promise<ArrayLike<number | ArrayLike<number>>>;
+    getGDALNoData?(): number | null;
+};
+export type DatumGridGeoTIFF = {
+    getImageCount(): Promise<number>;
+    getImage(index: number): Promise<DatumGridGeoTIFFImage>;
+};
+/** Load horizontal latitude/longitude offset bands in arcseconds on a geographic degree grid.
+ * Images are tried last-to-first, matching proj4js's child-before-parent GeoTIFF ordering.
+ */
+export declare function loadGeoTIFFGrid(tiff: DatumGridGeoTIFF): Promise<DatumGrid>;
+//# sourceMappingURL=geotiff.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/geotiff.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,SAAS,EAAU,mBAAgB;AAEhD,0FAA0F;AAC1F,MAAM,MAAM,qBAAqB,GAAG;IAClC,QAAQ,IAAI,MAAM,CAAC;IACnB,SAAS,IAAI,MAAM,CAAC;IACpB,cAAc,IAAI,MAAM,EAAE,CAAC;IAC3B,aAAa,EACT;QAAC,eAAe,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,CAAA;KAAC,GACrC;QAAC,QAAQ,CAAC,IAAI,EAAE,iBAAiB,GAAG,OAAO,CAAA;KAAC,CAAC;IACjD,WAAW,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC9D,aAAa,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC;CACjC,CAAC;AACF,MAAM,MAAM,gBAAgB,GAAG;IAC7B,aAAa,IAAI,OAAO,CAAC,MAAM,CAAC,CAAC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,GAAG,OAAO,CAAC,qBAAqB,CAAC,CAAC;CACzD,CAAC;AACF;;GAEG;AACH,wBAAsB,eAAe,CAAC,IAAI,EAAE,gBAAgB,GAAG,OAAO,CAAC,SAAS,CAAC,CA6DhF"}

package/dist/experimental/grids/geotiff.js

@@ -0,0 +1,59 @@
+// 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: GeoTIFF node orientation directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../PROJ4-LICENSE.md. Modified for explicit async loading, validation and nodata handling; no geotiff runtime import.
+import { createDatumGrid } from "./grid.js";
+/** Load horizontal latitude/longitude offset bands in arcseconds on a geographic degree grid.
+ * Images are tried last-to-first, matching proj4js's child-before-parent GeoTIFF ordering.
+ */
+export async function loadGeoTIFFGrid(tiff) {
+    const count = await tiff.getImageCount();
+    if (!Number.isSafeInteger(count) || count < 1)
+        throw new Error('GeoTIFF requires at least one image');
+    const subgrids = [];
+    for (let n = count - 1; n >= 0; n--) {
+        const image = await tiff.getImage(n), width = image.getWidth(), height = image.getHeight();
+        if (![width, height].every(v => Number.isSafeInteger(v) && v >= 2))
+            throw new Error('Invalid GeoTIFF grid dimensions');
+        const bbox = image.getBoundingBox();
+        const directory = image.fileDirectory;
+        const scale = 'getValue' in directory ? directory.getValue('ModelPixelScale') : directory.ModelPixelScale;
+        if (!scale || typeof scale !== 'object' || !('0' in scale) || !('1' in scale))
+            throw new Error('GeoTIFF requires ModelPixelScale');
+        const sx = Number(scale[0]), sy = Number(scale[1]);
+        if (bbox.length !== 4 ||
+            !bbox.every(Number.isFinite) ||
+            !(sx > 0 && sy > 0 && Number.isFinite(sx) && Number.isFinite(sy)))
+            throw new Error('Invalid GeoTIFF geographic grid geometry');
+        const rasters = await image.readRasters(), latitude = rasters[0], longitude = rasters[1];
+        if (!latitude ||
+            !longitude ||
+            typeof latitude === 'number' ||
+            typeof longitude === 'number' ||
+            latitude.length !== width * height ||
+            longitude.length !== width * height)
+            throw new Error('GeoTIFF requires latitude and longitude offset bands');
+        const nodata = image.getGDALNoData?.();
+        const shifts = new Float64Array(width * height * 2);
+        let node = 0;
+        for (let row = height - 1; row >= 0; row--)
+            for (let col = width - 1; col >= 0; col--) {
+                const index = row * width + col, lat = latitude[index], lon = longitude[index];
+                const missing = lat === nodata || lon === nodata || !Number.isFinite(lat) || !Number.isFinite(lon);
+                shifts[node++] = missing ? NaN : (-lon * Math.PI) / 648000;
+                shifts[node++] = missing ? NaN : (lat * Math.PI) / 648000;
+            }
+        subgrids.push({
+            origin: [
+                (-(bbox[0] + (width - 1) * sx) * Math.PI) / 180,
+                ((bbox[3] - (height - 1) * sy) * Math.PI) / 180
+            ],
+            step: [(sx * Math.PI) / 180, (sy * Math.PI) / 180],
+            size: [width, height],
+            shifts
+        });
+    }
+    return createDatumGrid(subgrids);
+}
+//# sourceMappingURL=geotiff.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff.js","sourceRoot":"","sources":["../../../src/experimental/grids/geotiff.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,4DAA4D;AAC5D,iOAAiO;AACjO,OAAO,EAAC,eAAe,EAAC,kBAAe;AAkBvC;;GAEG;AACH,MAAM,CAAC,KAAK,UAAU,eAAe,CAAC,IAAsB;IAC1D,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,aAAa,EAAE,CAAC;IACzC,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC;QAC3C,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IACzD,MAAM,QAAQ,GAAc,EAAE,CAAC;IAC/B,KAAK,IAAI,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAClC,KAAK,GAAG,KAAK,CAAC,QAAQ,EAAE,EACxB,MAAM,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAChE,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;QACrD,MAAM,IAAI,GAAG,KAAK,CAAC,cAAc,EAAE,CAAC;QACpC,MAAM,SAAS,GAAG,KAAK,CAAC,aAAa,CAAC;QACtC,MAAM,KAAK,GACT,UAAU,IAAI,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,QAAQ,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,eAAe,CAAC;QAC9F,IAAI,CAAC,KAAK,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,CAAC,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,KAAK,CAAC;YAC3E,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;QACtD,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EACzB,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACxB,IACE,IAAI,CAAC,MAAM,KAAK,CAAC;YACjB,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAC5B,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;YAEjE,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;QAC9D,MAAM,OAAO,GAAG,MAAM,KAAK,CAAC,WAAW,EAAE,EACvC,QAAQ,GAAG,OAAO,CAAC,CAAC,CAAC,EACrB,SAAS,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACzB,IACE,CAAC,QAAQ;YACT,CAAC,SAAS;YACV,OAAO,QAAQ,KAAK,QAAQ;YAC5B,OAAO,SAAS,KAAK,QAAQ;YAC7B,QAAQ,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM;YAClC,SAAS,CAAC,MAAM,KAAK,KAAK,GAAG,MAAM;YAEnC,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,MAAM,MAAM,GAAG,KAAK,CAAC,aAAa,EAAE,EAAE,CAAC;QACvC,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC;QACpD,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,GAAG,GAAG,MAAM,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE;YACxC,KAAK,IAAI,GAAG,GAAG,KAAK,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;gBAC1C,MAAM,KAAK,GAAG,GAAG,GAAG,KAAK,GAAG,GAAG,EAC7B,GAAG,GAAG,QAAQ,CAAC,KAAK,CAAC,EACrB,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;gBACzB,MAAM,OAAO,GACX,GAAG,KAAK,MAAM,IAAI,GAAG,KAAK,MAAM,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;gBACrF,MAAM,CAAC,IAAI,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC;gBAC3D,MAAM,CAAC,IAAI,EAAE,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC;YAC5D,CAAC;QACH,QAAQ,CAAC,IAAI,CAAC;YACZ,MAAM,EAAE;gBACN,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG;gBAC/C,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG;aAChD;YACD,IAAI,EAAE,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YAClD,IAAI,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC;YACrB,MAAM;SACP,CAAC,CAAC;IACL,CAAC;IACD,OAAO,eAAe,CAAC,QAAQ,CAAC,CAAC;AACnC,CAAC"}

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

@@ -0,0 +1,4 @@
+import type { DatumGrid, Subgrid } from "./types.js";
+/** Readers own their decoded arrays; keep the snapshot private behind a frozen executable grid. */
+export declare function createDatumGrid(subgrids: Subgrid[]): DatumGrid;
+//# sourceMappingURL=grid.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"grid.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/grid.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,SAAS,EAAE,OAAO,EAAC,mBAAgB;AAwFhD,mGAAmG;AACnG,wBAAgB,eAAe,CAAC,QAAQ,EAAE,OAAO,EAAE,GAAG,SAAS,CA8C9D"}

package/dist/experimental/grids/grid.js

@@ -0,0 +1,111 @@
+import { wrapLongitude } from "../parameters.js";
+const 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;
+    }
+    // Seed shifted boundary nodes within the measured displacement envelope.
+    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) {
+            // Never accept the clamped extrapolation outside actual source coverage.
+            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');
+}
+/** Readers own their decoded arrays; keep the snapshot private behind a frozen executable grid. */
+export 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;
+    };
+    // This private point never escapes or invokes application hooks. A returned
+    // tuple owns its values; the next call resets both coordinates even after failure.
+    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] : undefined;
+        }
+    });
+}
+//# sourceMappingURL=grid.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"grid.js","sourceRoot":"","sources":["../../../src/experimental/grids/grid.ts"],"names":[],"mappings":"AAMA,OAAO,EAAC,aAAa,EAAC,yBAAsB;AAE5C,MAAM,SAAS,GAAG,KAAK,CAAC;AAIxB,SAAS,WAAW,CAClB,IAAqB,EACrB,CAAS,EACT,CAAS,EACT,MAAa,EACb,KAAK,GAAG,KAAK;IAEb,IAAI,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAC3C,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EACxB,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,IACE,CAAC,KAAK;QACN,CAAC,GAAG,GAAG,CAAC,KAAK,IAAI,GAAG,GAAG,CAAC,KAAK,IAAI,GAAG,GAAG,KAAK,GAAG,CAAC,GAAG,KAAK,IAAI,GAAG,GAAG,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;QAErF,OAAO,KAAK,CAAC;IACf,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAC5C,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAC7C,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,EAC5C,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5C,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,EACf,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IACd,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,CAAC,CAAC;IACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,EACjB,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,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;QACzD,IAAI,MAAM,KAAK,CAAC;YAAE,SAAS;QAC3B,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QAC9D,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,EAC/B,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAClC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YAAE,OAAO,KAAK,CAAC;QACvE,EAAE,IAAI,MAAM,GAAG,MAAM,CAAC;QACtB,EAAE,IAAI,MAAM,GAAG,MAAM,CAAC;IACxB,CAAC;IACD,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;IACd,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC;IACd,OAAO,IAAI,CAAC;AACd,CAAC;AACD,SAAS,YAAY,CACnB,IAAqB,EACrB,CAAS,EACT,CAAS,EACT,OAAgB,EAChB,MAAa;IAEb,IAAI,CAAC,OAAO,EAAE,CAAC;QACb,IAAI,CAAC,WAAW,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC;YAAE,OAAO,KAAK,CAAC;QACnD,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC;QACd,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC;QACd,OAAO,IAAI,CAAC;IACd,CAAC;IACD,yEAAyE;IACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,IACE,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,SAAS;YACrD,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,SAAS;YAEzF,OAAO,KAAK,CAAC;IACjB,CAAC;IACD,IAAI,CAAC,WAAW,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC;QAAE,OAAO,KAAK,CAAC;IACzD,IAAI,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,EACvB,MAAM,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACxB,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,EAAE,EAAE,SAAS,EAAE,EAAE,CAAC;QACpD,IAAI,CAAC,WAAW,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC;YAAE,OAAO,KAAK,CAAC;QACnE,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,EAC9B,EAAE,GAAG,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;QAC7B,MAAM,IAAI,EAAE,CAAC;QACb,MAAM,IAAI,EAAE,CAAC;QACb,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,IAAI,SAAS,EAAE,CAAC;YACtD,yEAAyE;YACzE,IAAI,CAAC,WAAW,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;gBAAE,OAAO,KAAK,CAAC;YAC7D,MAAM,CAAC,CAAC,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;YACjC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,CAAC;QACd,CAAC;IACH,CAAC;IACD,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;AAC/D,CAAC;AACD,mGAAmG;AACnG,MAAM,UAAU,eAAe,CAAC,QAAmB;IACjD,IAAI,CAAC,QAAQ,CAAC,MAAM;QAAE,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAClF,MAAM,QAAQ,GAAsB,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;QACtD,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;QAClC,IACE,CAAC,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC9D,KAAK,GAAG,MAAM,GAAG,CAAC,KAAK,IAAI,CAAC,MAAM,CAAC,MAAM;YACzC,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YACnC,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YAElD,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;QAClE,MAAM,QAAQ,GAAqB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC1C,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC;YACxD,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;YACjC,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;gBACxB,QAAQ,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;QACzE,CAAC;QACD,OAAO,EAAC,GAAG,IAAI,EAAE,QAAQ,EAAC,CAAC;IAC7B,CAAC,CAAC,CAAC;IACH,MAAM,YAAY,GAAG,CAAC,KAAY,EAAE,OAAgB,EAAW,EAAE;QAC/D,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,EACvB,QAAQ,GAAG,KAAK,CAAC,CAAC,CAAC;QACrB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YAC3D,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;QACrD,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE,CAAC;YAC5B,IAAI,YAAY,CAAC,IAAI,EAAE,CAAC,SAAS,EAAE,QAAQ,EAAE,OAAO,EAAE,KAAK,CAAC,EAAE,CAAC;gBAC7D,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;gBACnB,OAAO,IAAI,CAAC;YACd,CAAC;QACH,CAAC;QACD,KAAK,CAAC,CAAC,GAAG,SAAS,CAAC;QACpB,KAAK,CAAC,CAAC,GAAG,QAAQ,CAAC;QACnB,OAAO,KAAK,CAAC;IACf,CAAC,CAAC;IACF,4EAA4E;IAC5E,mFAAmF;IACnF,MAAM,YAAY,GAAU,EAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAC,CAAC;IAC7C,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,YAAY,EAAE,QAAQ,CAAC,MAAM;QAC7B,YAAY;QACZ,KAAK,CAAC,SAAiB,EAAE,QAAgB,EAAE,OAAgB;YACzD,YAAY,CAAC,CAAC,GAAG,SAAS,CAAC;YAC3B,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC;YAC1B,OAAO,YAAY,CAAC,YAAY,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;QAC5F,CAAC;KACF,CAAC,CAAC;AACL,CAAC"}

package/dist/experimental/grids/gtx.cjs

@@ -0,0 +1,91 @@
+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/gtx.js
+var gtx_exports = {};
+__export(gtx_exports, {
+  parseGTXGrid: () => parseGTXGrid
+});
+module.exports = __toCommonJS(gtx_exports);
+
+// dist/experimental/grids/vertical.js
+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;
+    }
+  });
+}
+
+// dist/experimental/grids/gtx.js
+function parseGTXGrid(buffer) {
+  const data = new DataView(buffer);
+  if (data.byteLength < 40)
+    throw new Error("Truncated GTX header");
+  const height = data.getInt32(32), width = data.getInt32(36);
+  if (width < 2 || height < 2 || data.byteLength !== 40 + width * height * 4)
+    throw new Error("Invalid GTX dimensions or data length");
+  const offsets = new Float32Array(width * height);
+  for (let i = 0; i < offsets.length; i++) {
+    const value = data.getFloat32(40 + i * 4);
+    offsets[i] = Math.abs(value) > 1e3 ? NaN : value;
+  }
+  return createVerticalGrid({
+    origin: [data.getFloat64(8), data.getFloat64(0)],
+    step: [data.getFloat64(24), data.getFloat64(16)],
+    size: [width, height],
+    offsets,
+    noData: Math.fround(-88.8888)
+  });
+}
+//# sourceMappingURL=gtx.cjs.map

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

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../src/experimental/grids/gtx.ts", "../../../src/experimental/grids/vertical.ts"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;;;;;;;ACoBM,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;;;ADtEM,SAAU,aAAa,QAAmB;AAC9C,QAAM,OAAO,IAAI,SAAS,MAAM;AAChC,MAAI,KAAK,aAAa;AAAI,UAAM,IAAI,MAAM,sBAAsB;AAChE,QAAM,SAAS,KAAK,SAAS,EAAE,GAC7B,QAAQ,KAAK,SAAS,EAAE;AAC1B,MAAI,QAAQ,KAAK,SAAS,KAAK,KAAK,eAAe,KAAK,QAAQ,SAAS;AACvE,UAAM,IAAI,MAAM,uCAAuC;AACzD,QAAM,UAAU,IAAI,aAAa,QAAQ,MAAM;AAC/C,WAAS,IAAI,GAAG,IAAI,QAAQ,QAAQ,KAAK;AACvC,UAAM,QAAQ,KAAK,WAAW,KAAK,IAAI,CAAC;AACxC,YAAQ,CAAC,IAAI,KAAK,IAAI,KAAK,IAAI,MAAO,MAAM;EAC9C;AACA,SAAO,mBAAmB;IACxB,QAAQ,CAAC,KAAK,WAAW,CAAC,GAAG,KAAK,WAAW,CAAC,CAAC;IAC/C,MAAM,CAAC,KAAK,WAAW,EAAE,GAAG,KAAK,WAAW,EAAE,CAAC;IAC/C,MAAM,CAAC,OAAO,MAAM;IACpB;IACA,QAAQ,KAAK,OAAO,QAAQ;GAC7B;AACH;",
+  "names": []
+}

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

@@ -0,0 +1,7 @@
+import type { VerticalGrid } from "./types.js";
+/** Read a big-endian GTX grid (float32 metre offsets) into an owned snapshot.
+ * The conventional -88.8888 float32 sentinel and offsets outside ±1000 metres are nodata, matching PROJ GTX semantics.
+ * Fetching/decompression remain caller-owned; the prepared grid is synchronous.
+ */
+export declare function parseGTXGrid(buffer: ArrayBuffer): VerticalGrid;
+//# sourceMappingURL=gtx.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"gtx.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/gtx.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,YAAY,EAAC,mBAAgB;AAE1C;;;GAGG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,WAAW,GAAG,YAAY,CAmB9D"}

package/dist/experimental/grids/gtx.js

@@ -0,0 +1,30 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original reader for the NOAA GTX format documented by GDAL and PROJ. No upstream implementation is copied. Independent tests use PROJ 9.5.1.
+import { createVerticalGrid } from "./vertical.js";
+/** Read a big-endian GTX grid (float32 metre offsets) into an owned snapshot.
+ * The conventional -88.8888 float32 sentinel and offsets outside ±1000 metres are nodata, matching PROJ GTX semantics.
+ * Fetching/decompression remain caller-owned; the prepared grid is synchronous.
+ */
+export function parseGTXGrid(buffer) {
+    const data = new DataView(buffer);
+    if (data.byteLength < 40)
+        throw new Error('Truncated GTX header');
+    const height = data.getInt32(32), width = data.getInt32(36);
+    if (width < 2 || height < 2 || data.byteLength !== 40 + width * height * 4)
+        throw new Error('Invalid GTX dimensions or data length');
+    const offsets = new Float32Array(width * height);
+    for (let i = 0; i < offsets.length; i++) {
+        const value = data.getFloat32(40 + i * 4);
+        offsets[i] = Math.abs(value) > 1000 ? NaN : value;
+    }
+    return createVerticalGrid({
+        origin: [data.getFloat64(8), data.getFloat64(0)],
+        step: [data.getFloat64(24), data.getFloat64(16)],
+        size: [width, height],
+        offsets,
+        noData: Math.fround(-88.8888)
+    });
+}
+//# sourceMappingURL=gtx.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"gtx.js","sourceRoot":"","sources":["../../../src/experimental/grids/gtx.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,iKAAiK;AACjK,OAAO,EAAC,kBAAkB,EAAC,sBAAmB;AAG9C;;;GAGG;AACH,MAAM,UAAU,YAAY,CAAC,MAAmB;IAC9C,MAAM,IAAI,GAAG,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC;IAClC,IAAI,IAAI,CAAC,UAAU,GAAG,EAAE;QAAE,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;IAClE,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,EAC9B,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IAC5B,IAAI,KAAK,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,IAAI,IAAI,CAAC,UAAU,KAAK,EAAE,GAAG,KAAK,GAAG,MAAM,GAAG,CAAC;QACxE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;IAC3D,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;IACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxC,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1C,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC;IACpD,CAAC;IACD,OAAO,kBAAkB,CAAC;QACxB,MAAM,EAAE,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QAChD,IAAI,EAAE,CAAC,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAChD,IAAI,EAAE,CAAC,KAAK,EAAE,MAAM,CAAC;QACrB,OAAO;QACP,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,OAAO,CAAC;KAC9B,CAAC,CAAC;AACL,CAAC"}