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

package/dist/projection/experimental/grids/vertical-geotiff.cjs

@@ -0,0 +1,195 @@
+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/projection/experimental/grids/vertical-geotiff.js
+var vertical_geotiff_exports = {};
+__export(vertical_geotiff_exports, {
+  loadVerticalGeoTIFFGrid: () => loadVerticalGeoTIFFGrid
+});
+module.exports = __toCommonJS(vertical_geotiff_exports);
+
+// dist/projection/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/projection/experimental/grids/scalar-geotiff.js
+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 === void 0)
+    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;
+}
+async function prepareScalarGeoTIFFGrids(tiff) {
+  var _a;
+  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") !== void 0 && tag("NewSubfileType") !== 0)
+      throw new Error("Vertical GeoTIFF overviews and masks are unsupported");
+    if (tag("ModelTransformation") !== void 0)
+      throw new Error("Vertical GeoTIFF requires axis-aligned tiepoint/scale geometry");
+    const keys = image.getGeoKeys();
+    if ((keys == null ? void 0 : keys["GTModelTypeGeoKey"]) !== 2 || keys["GeogAngularUnitsGeoKey"] !== void 0 && keys["GeogAngularUnitsGeoKey"] !== 9102)
+      throw new Error("Vertical GeoTIFF requires a geographic degree grid");
+    if (keys["GeogPrimeMeridianGeoKey"] !== void 0 && keys["GeogPrimeMeridianGeoKey"] !== 8901 || keys["GeogPrimeMeridianLongGeoKey"] !== void 0 && 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 == null ? void 0 : metadata["TYPE"]) !== "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL" || (band == null ? void 0 : band["DESCRIPTION"]) !== "geoid_undulation")
+      throw new Error("Vertical GeoTIFF requires a geographic-to-vertical geoid_undulation band");
+    if (band["UNITTYPE"] !== void 0 && band["UNITTYPE"] !== "metre")
+      throw new Error("Vertical GeoTIFF geoid undulations must be in metres");
+    if (metadata["interpolation_method"] !== void 0 && 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");
+    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 = (_a = image.getGDALNoData) == null ? void 0 : _a.call(image);
+    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);
+}
+
+// dist/projection/experimental/grids/vertical-geotiff.js
+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 !== void 0)
+          return value;
+      }
+      return void 0;
+    }
+  });
+}
+//# sourceMappingURL=vertical-geotiff.cjs.map

package/dist/projection/experimental/grids/vertical-geotiff.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/projection/experimental/grids/vertical-geotiff.ts", "../../../../src/projection/experimental/grids/vertical.ts", "../../../../src/projection/experimental/grids/scalar-geotiff.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;;;ACzBA,SAAS,eAAe,MAA+B,OAAa;AAClE,QAAM,QAAQ,KAAK,OAAO,KAAK;AAC/B,QAAM,OAAO,MAAM,MAAM,KAAK,eAAa,UAAU,UAAU,CAAC;AAChE,MAAI,CAAC;AAAM,UAAM,IAAI,MAAM,8CAA8C;AACzE,SAAO;IACL,UAAU,MAAM,MAAM;IACtB,WAAW,MAAM,MAAM;IACvB,YAAY,MAAM,MAAM;IACxB,iBAAiB,YAAW,WAAW,IAAI,KAAK,WAAW,MAAM;IACjE,eAAe,MAAM,MAAM;IAC3B,eAAe,MAAM;IACrB,aAAa,YAAY,CAAC,KAAK,IAAI;;AAEvC;AAEA,SAAS,OAAO,OAAgB,UAAgB;AAC9C,MAAI,UAAU;AAAW,WAAO;AAChC,MACG,OAAO,UAAU,YAAY,OAAO,UAAU,YAC9C,OAAO,UAAU,YAAY,CAAC,MAAM,KAAI,KACzC,CAAC,OAAO,SAAS,OAAO,KAAK,CAAC;AAE9B,UAAM,IAAI,MAAM,0CAA0C;AAC5D,SAAO,OAAO,KAAK;AACrB;AACA,SAAS,OAAO,OAAgB,QAAc;AAC5C,MAAI,CAAC,SAAS,OAAO,UAAU,YAAY,EAAE,YAAY,UAAU,MAAM,WAAW;AAClF,UAAM,IAAI,MAAM,wDAAwD;AAC1E,QAAM,SAAS,MAAM,KAAK,KAA0B;AACpD,MAAI,CAAC,OAAO,MAAM,OAAO,QAAQ;AAAG,UAAM,IAAI,MAAM,mCAAmC;AACvF,SAAO;AACT;AAQA,eAAsB,0BACpB,MAAmD;AAhGrD;AAkGE,QAAM,QAAQ,YAAY,OAAO,KAAK,OAAO,SAAS,MAAM,KAAK,cAAa;AAC9E,MAAI,CAAC,OAAO,cAAc,KAAK,KAAK,QAAQ;AAC1C,UAAM,IAAI,MAAM,8CAA8C;AAChE,QAAM,WAQA,CAAA;AACN,WAAS,IAAI,GAAG,IAAI,OAAO,KAAK;AAC9B,UAAM,QAAQ,YAAY,OAAO,eAAe,MAAM,CAAC,IAAI,MAAM,KAAK,SAAS,CAAC;AAChF,UAAM,YAAY,MAAM;AACxB,UAAM,MAAM,CAAC,SACX,cAAc,YAAY,UAAU,SAAS,IAAI,IAAI,UAAU,IAAI;AACrE,QAAI,IAAI,gBAAgB,MAAM,UAAa,IAAI,gBAAgB,MAAM;AACnE,YAAM,IAAI,MAAM,sDAAsD;AACxE,QAAI,IAAI,qBAAqB,MAAM;AACjC,YAAM,IAAI,MAAM,gEAAgE;AAClF,UAAM,OAAO,MAAM,WAAU;AAC7B,SACE,6BAAO,0BAAyB,KAC/B,KAAK,wBAAwB,MAAM,UAAa,KAAK,wBAAwB,MAAM;AAEpF,YAAM,IAAI,MAAM,oDAAoD;AACtE,QACG,KAAK,yBAAyB,MAAM,UAAa,KAAK,yBAAyB,MAAM,QACrF,KAAK,6BAA6B,MAAM,UACvC,KAAK,6BAA6B,MAAM;AAE1C,YAAM,IAAI,MAAM,gDAAgD;AAClE,UAAM,aAAa,KAAK,oBAAoB;AAC5C,QAAI,eAAe,KAAK,eAAe;AACrC,YAAM,IAAI,MAAM,gEAAgE;AAClF,UAAM,WAAW,MAAM,MAAM,gBAAe;AAC5C,UAAM,OAAO,MAAM,MAAM,gBAAgB,CAAC;AAC1C,SACE,qCAAW,aAAY,6CACvB,6BAAO,oBAAmB;AAE1B,YAAM,IAAI,MAAM,0EAA0E;AAC5F,QAAI,KAAK,UAAU,MAAM,UAAa,KAAK,UAAU,MAAM;AACzD,YAAM,IAAI,MAAM,sDAAsD;AACxE,QACE,SAAS,sBAAsB,MAAM,UACrC,SAAS,sBAAsB,MAAM;AAErC,YAAM,IAAI,MAAM,uDAAuD;AACzE,UAAM,QAAQ,OAAO,KAAK,OAAO,GAAG,CAAC,GACnC,SAAS,OAAO,KAAK,QAAQ,GAAG,CAAC;AACnC,UAAM,CAAC,IAAI,EAAE,IAAI,OAAO,IAAI,iBAAiB,GAAG,CAAC;AACjD,UAAM,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,OAAO,IAAI,eAAe,GAAG,CAAC;AACzD,UAAM,QAAQ,MAAM,SAAQ,GAC1B,SAAS,MAAM,UAAS;AAC1B,QACE,EAAE,KAAK,KAAK,KAAK,MACjB,MAAM,KACN,MAAM,KACN,CAAC,CAAC,OAAO,MAAM,EAAE,MAAM,OAAK,OAAO,cAAc,CAAC,KAAK,KAAK,CAAC;AAE7D,YAAM,IAAI,MAAM,iDAAiD;AACnE,UAAM,OAAO,eAAe,IAAI,MAAM;AACtC,UAAM,OAAO,KAAK,OAAO,KAAK;AAC9B,UAAM,QAAQ,KAAK,OAAO,KAAK;AAC/B,UAAM,QAAQ,SAAS,SAAS,KAAK,IACnC,OAAO,QAAQ,QAAQ,KAAK;AAC9B,QACE,CAAC,CAAC,MAAM,OAAO,MAAM,KAAK,EAAE,MAAM,OAAO,QAAQ,KACjD,EAAE,OAAO,QAAQ,QAAQ,UACzB,OAAO,OAAO,OACd,QAAQ,OACR,QAAQ,MACR,CAAC,OAAO,cAAc,QAAQ,MAAM;AAEpC,YAAM,IAAI,MAAM,4CAA4C;AAG9D,eAAW,SAAS,UAAU;AAC5B,YAAM,WACJ,OAAO,MAAM,QAAQ,OAAO,MAAM,QAAQ,QAAQ,MAAM,SAAS,QAAQ,MAAM;AACjF,UACE,YACA,EACE,QAAQ,MAAM,QACd,QAAQ,MAAM,QACd,SAAS,MAAM,SACf,SAAS,MAAM,SACf,MAAM,MAAM,MACZ,MAAM,MAAM;AAGd,cAAM,IAAI,MACR,gFAAgF;IAEtF;AACA,UAAM,UAAU,MAAM,MAAM,YAAY,EAAC,SAAS,CAAC,CAAC,GAAG,YAAY,MAAK,CAAC;AACzE,UAAM,SAAS,QAAQ,CAAC;AACxB,QAAI,CAAC,UAAU,OAAO,WAAW,YAAY,OAAO,WAAW,QAAQ;AACrE,YAAM,IAAI,MAAM,kEAAkE;AACpF,UAAM,kBAAiB,WAAM,kBAAN;AAEvB,UAAM,SACJ,kBAAkB,gBAAgB,OAAO,mBAAmB,WACxD,KAAK,OAAO,cAAc,IAC1B;AACN,UAAM,UAAU,IAAI,aAAa,QAAQ,MAAM;AAC/C,aAAS,MAAM,GAAG,MAAM,QAAQ;AAC9B,eAAS,MAAM,GAAG,MAAM,OAAO,OAAO;AACpC,cAAM,MAAM,OAAO,MAAM,QAAQ,GAAG;AACpC,cAAM,QAAQ,QAAQ,UAAU,CAAC,OAAO,SAAS,GAAG,IAAI,MAAM,MAAM,QAAQ;AAC5E,YAAI,OAAO,SAAS,GAAG,KAAK,QAAQ,UAAU,CAAC,OAAO,SAAS,KAAK;AAClE,gBAAM,IAAI,MAAM,4DAA4D;AAC9E,iBAAS,SAAS,IAAI,OAAO,QAAQ,GAAG,IAAI;MAC9C;AACF,aAAS,KAAK;MACZ,MAAM,mBAAmB;QACvB,QAAQ,CAAC,MAAM,KAAK;QACpB,MAAM,CAAC,IAAI,EAAE;QACb,MAAM,CAAC,OAAO,MAAM;QACpB;OACD;MACD;MACA;MACA;MACA;MACA;MACA;KACD;EACH;AACA,SAAO,SAAS,IAAI,WAAS,MAAM,IAAI;AACzC;;;AFxNA,eAAsB,wBACpB,MAAmD;AAEnD,QAAM,QAAQ,MAAM,0BAA0B,IAAI;AAClD,SAAO,OAAO,OAAO;IACnB,UAAU,WAAmB,UAAgB;AAC3C,eAAS,IAAI,MAAM,SAAS,GAAG,KAAK,GAAG,KAAK;AAC1C,cAAM,QAAQ,MAAM,CAAC,EAAE,UAAU,WAAW,QAAQ;AACpD,YAAI,UAAU;AAAW,iBAAO;MAClC;AACA,aAAO;IACT;GACD;AACH;",
+  "names": []
+}

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

@@ -0,0 +1,6 @@
+import type { VerticalGrid } from "./types.js";
+import type { VerticalGridGeoTIFF, VerticalGridGeoTIFFData } from "./scalar-geotiff.js";
+export type { VerticalGridGeoTIFF, VerticalGridGeoTIFFImage, VerticalGridGeoTIFFData } from "./scalar-geotiff.js";
+/** Prepare owned geoid grids; later nested images take precedence with nodata fallback. */
+export declare function loadVerticalGeoTIFFGrid(tiff: VerticalGridGeoTIFF | VerticalGridGeoTIFFData): Promise<VerticalGrid>;
+//# sourceMappingURL=vertical-geotiff.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"vertical-geotiff.d.ts","sourceRoot":"","sources":["../../../../src/projection/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/projection/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/projection/experimental/grids/vertical-geotiff.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical-geotiff.js","sourceRoot":"","sources":["../../../../src/projection/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/projection/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/projection/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/projection/experimental/grids/vertical.cjs.map

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/projection/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/projection/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/earth 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/projection/experimental/grids/vertical.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical.d.ts","sourceRoot":"","sources":["../../../../src/projection/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/projection/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/earth 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/projection/experimental/grids/vertical.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"vertical.js","sourceRoot":"","sources":["../../../../src/projection/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/projection/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/projection/experimental/helmert-flat.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"helmert-flat.d.ts","sourceRoot":"","sources":["../../../src/projection/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/projection/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/projection/experimental/helmert-flat.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"helmert-flat.js","sourceRoot":"","sources":["../../../src/projection/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"}