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

package/dist/experimental/grids/velocity-geotiff.cjs

@@ -0,0 +1,247 @@
+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/velocity-geotiff.js
+var velocity_geotiff_exports = {};
+__export(velocity_geotiff_exports, {
+  loadVelocityGeoTIFFGrid: () => loadVelocityGeoTIFFGrid
+});
+module.exports = __toCommonJS(velocity_geotiff_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/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/experimental/grids/velocity.js
+function velocityGridFromComponents(east, north, up, scale) {
+  return Object.freeze({
+    sample(longitude, latitude, output) {
+      const e = east.getOffset(longitude, latitude), n = north.getOffset(longitude, latitude), u = up.getOffset(longitude, latitude);
+      if (e === void 0 || n === void 0 || u === void 0)
+        return false;
+      const x = e * scale, y = n * scale, z = u * scale;
+      if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))
+        throw new Error("Non-finite velocity interpolation");
+      output.x = x;
+      output.y = y;
+      output.z = z;
+      return true;
+    }
+  });
+}
+
+// dist/experimental/grids/velocity-geotiff.js
+async function loadVelocityGeoTIFFGrid(tiff) {
+  const count = "images" in tiff ? tiff.images.length : await tiff.getImageCount();
+  if (!Number.isSafeInteger(count) || count < 1)
+    throw new Error("Velocity GeoTIFF requires images");
+  const images = [[], [], []];
+  const names = ["east_velocity", "north_velocity", "up_velocity"];
+  for (let index = 0; index < count; index++) {
+    const raster = "images" in tiff ? tiff.images[index] : void 0;
+    const image = raster ? void 0 : await tiff.getImage(index);
+    const metadata = raster ? raster.metadata : await image.getGDALMetadata();
+    if ((metadata == null ? void 0 : metadata["TYPE"]) !== "VELOCITY")
+      throw new Error("Velocity GeoTIFF requires TYPE=VELOCITY");
+    for (let component = 0; component < 3; component++) {
+      const band = raster == null ? void 0 : raster.bands.find((b) => b.index === component);
+      if (raster && !band)
+        throw new Error("Velocity GeoTIFF requires original ENU bands 0, 1, 2");
+      const bandMetadata = raster ? band.metadata : await image.getGDALMetadata(component);
+      if ((bandMetadata == null ? void 0 : bandMetadata["DESCRIPTION"]) !== names[component] || !["millimetres per year", "mm/year"].includes(bandMetadata == null ? void 0 : bandMetadata["UNITTYPE"]))
+        throw new Error("Velocity GeoTIFF requires explicit ENU descriptions and millimetres per year");
+      images[component].push({
+        getWidth: () => raster ? raster.width : image.getWidth(),
+        getHeight: () => raster ? raster.height : image.getHeight(),
+        getGeoKeys: () => raster ? raster.geoKeys : image.getGeoKeys(),
+        getGDALMetadata: (sample) => sample === 0 ? { ...bandMetadata, DESCRIPTION: "geoid_undulation", UNITTYPE: "metre" } : { ...metadata, TYPE: "VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL" },
+        getGDALNoData: () => {
+          var _a;
+          return raster ? raster.noData : ((_a = image.getGDALNoData) == null ? void 0 : _a.call(image)) ?? null;
+        },
+        fileDirectory: raster ? raster.fileDirectory : image.fileDirectory,
+        readRasters: async () => raster ? [band.data] : image.readRasters({ samples: [component], interleave: false })
+      });
+    }
+  }
+  const bands = await Promise.all(images.map((list) => prepareScalarGeoTIFFGrids({
+    getImageCount: async () => list.length,
+    getImage: async (index) => list[index]
+  })));
+  const grids = bands[0].map((east, index) => velocityGridFromComponents(east, bands[1][index], bands[2][index], 1e-3));
+  return Object.freeze({
+    sample(longitude, latitude, output) {
+      for (let i = grids.length - 1; i >= 0; i--)
+        if (grids[i].sample(longitude, latitude, output))
+          return true;
+      return false;
+    }
+  });
+}
+//# sourceMappingURL=velocity-geotiff.cjs.map

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

@@ -0,0 +1,6 @@
+{
+  "version": 3,
+  "sources": ["../../../src/experimental/grids/velocity-geotiff.ts", "../../../src/experimental/grids/vertical.ts", "../../../src/experimental/grids/scalar-geotiff.ts", "../../../src/experimental/grids/velocity.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;;;AC/MM,SAAU,2BACd,MACA,OACA,IACA,OAAa;AAEb,SAAO,OAAO,OAAO;IACnB,OAAO,WAAmB,UAAkB,QAAuB;AACjE,YAAM,IAAI,KAAK,UAAU,WAAW,QAAQ,GAC1C,IAAI,MAAM,UAAU,WAAW,QAAQ,GACvC,IAAI,GAAG,UAAU,WAAW,QAAQ;AACtC,UAAI,MAAM,UAAa,MAAM,UAAa,MAAM;AAAW,eAAO;AAClE,YAAM,IAAI,IAAI,OACZ,IAAI,IAAI,OACR,IAAI,IAAI;AACV,UAAI,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC,KAAK,CAAC,OAAO,SAAS,CAAC;AAClE,cAAM,IAAI,MAAM,mCAAmC;AACrD,aAAO,IAAI;AACX,aAAO,IAAI;AACX,aAAO,IAAI;AACX,aAAO;IACT;GACD;AACH;;;AH1BA,eAAsB,wBACpB,MAAmD;AAEnD,QAAM,QAAQ,YAAY,OAAO,KAAK,OAAO,SAAS,MAAM,KAAK,cAAa;AAC9E,MAAI,CAAC,OAAO,cAAc,KAAK,KAAK,QAAQ;AAC1C,UAAM,IAAI,MAAM,kCAAkC;AACpD,QAAM,SAAuC,CAAC,CAAA,GAAI,CAAA,GAAI,CAAA,CAAE;AACxD,QAAM,QAAQ,CAAC,iBAAiB,kBAAkB,aAAa;AAC/D,WAAS,QAAQ,GAAG,QAAQ,OAAO,SAAS;AAC1C,UAAM,SAAS,YAAY,OAAO,KAAK,OAAO,KAAK,IAAI;AACvD,UAAM,QAAQ,SAAS,SAAY,MAAO,KAA6B,SAAS,KAAK;AACrF,UAAM,WAAW,SAAS,OAAO,WAAW,MAAM,MAAM,gBAAe;AACvE,SAAI,qCAAW,aAAY;AACzB,YAAM,IAAI,MAAM,yCAAyC;AAC3D,aAAS,YAAY,GAAG,YAAY,GAAG,aAAa;AAClD,YAAM,OAAO,iCAAQ,MAAM,KAAK,OAAK,EAAE,UAAU;AACjD,UAAI,UAAU,CAAC;AAAM,cAAM,IAAI,MAAM,sDAAsD;AAC3F,YAAM,eAAe,SAAS,KAAK,WAAW,MAAM,MAAM,gBAAgB,SAAS;AACnF,WACE,6CAAe,oBAAmB,MAAM,SAAS,KACjD,CAAC,CAAC,wBAAwB,SAAS,EAAE,SAAS,6CAAe,WAAqB;AAElF,cAAM,IAAI,MACR,8EAA8E;AAKlF,aAAO,SAAS,EAAE,KAAK;QACrB,UAAU,MAAO,SAAS,OAAO,QAAQ,MAAM,SAAQ;QACvD,WAAW,MAAO,SAAS,OAAO,SAAS,MAAM,UAAS;QAC1D,YAAY,MAAO,SAAS,OAAO,UAAU,MAAM,WAAU;QAC7D,iBAAiB,YACf,WAAW,IACP,EAAC,GAAG,cAAc,aAAa,oBAAoB,UAAU,QAAO,IACpE,EAAC,GAAG,UAAU,MAAM,yCAAwC;QAClE,eAAe,MAAG;AAxD1B;AAwD8B,0BAAS,OAAO,WAAU,WAAM,kBAAN,mCAA2B;;QAC3E,eAAe,SAAS,OAAO,gBAAgB,MAAM;QACrD,aAAa,YACX,SAAS,CAAC,KAAK,IAAI,IAAI,MAAM,YAAY,EAAC,SAAS,CAAC,SAAS,GAAG,YAAY,MAAK,CAAC;OACrF;IACH;EACF;AACA,QAAM,QAAQ,MAAM,QAAQ,IAC1B,OAAO,IAAI,UACT,0BAA0B;IACxB,eAAe,YAAY,KAAK;IAChC,UAAU,OAAM,UAAS,KAAK,KAAK;GACpC,CAAC,CACH;AAEH,QAAM,QAAQ,MAAM,CAAC,EAAE,IAAI,CAAC,MAAM,UAChC,2BAA2B,MAAM,MAAM,CAAC,EAAE,KAAK,GAAG,MAAM,CAAC,EAAE,KAAK,GAAG,IAAK,CAAC;AAE3E,SAAO,OAAO,OAAO;IACnB,OAAO,WAAW,UAAU,QAAM;AAChC,eAAS,IAAI,MAAM,SAAS,GAAG,KAAK,GAAG;AACrC,YAAI,MAAM,CAAC,EAAE,OAAO,WAAW,UAAU,MAAM;AAAG,iBAAO;AAC3D,aAAO;IACT;GACD;AACH;",
+  "names": []
+}

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

@@ -0,0 +1,11 @@
+import type { VelocityGrid } from "./velocity.js";
+import type { VerticalGridGeoTIFF, VerticalGridGeoTIFFData, VerticalGridGeoTIFFImage } from "./vertical-geotiff.js";
+export type VelocityGridGeoTIFF = VerticalGridGeoTIFF;
+export type VelocityGridGeoTIFFImage = VerticalGridGeoTIFFImage;
+export type VelocityGridGeoTIFFData = VerticalGridGeoTIFFData;
+/** Prepare geographic, north-up GTG VELOCITY images with explicit first-three ENU
+ * bands in millimetres per year. Accuracy bands are ignored. No model/frame/epoch
+ * inference. Decoded loaders.gl rasters and structural geotiff.js objects are accepted.
+ */
+export declare function loadVelocityGeoTIFFGrid(tiff: VelocityGridGeoTIFF | VelocityGridGeoTIFFData): Promise<VelocityGrid>;
+//# sourceMappingURL=velocity-geotiff.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"velocity-geotiff.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/velocity-geotiff.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,YAAY,EAAC,sBAAmB;AAC7C,OAAO,KAAK,EACV,mBAAmB,EACnB,uBAAuB,EACvB,wBAAwB,EACzB,8BAA2B;AAC5B,MAAM,MAAM,mBAAmB,GAAG,mBAAmB,CAAC;AACtD,MAAM,MAAM,wBAAwB,GAAG,wBAAwB,CAAC;AAChE,MAAM,MAAM,uBAAuB,GAAG,uBAAuB,CAAC;AAE9D;;;GAGG;AACH,wBAAsB,uBAAuB,CAC3C,IAAI,EAAE,mBAAmB,GAAG,uBAAuB,GAClD,OAAO,CAAC,YAAY,CAAC,CA2DvB"}

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

@@ -0,0 +1,61 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adapter inspired by PROJ's Geodetic TIFF Grid specification, not a fork. Caller-owned TIFF decoding; no proj4js, PROJ or TIFF decoder code is copied here.
+import { prepareScalarGeoTIFFGrids } from "./scalar-geotiff.js";
+import { velocityGridFromComponents } from "./velocity.js";
+/** Prepare geographic, north-up GTG VELOCITY images with explicit first-three ENU
+ * bands in millimetres per year. Accuracy bands are ignored. No model/frame/epoch
+ * inference. Decoded loaders.gl rasters and structural geotiff.js objects are accepted.
+ */
+export async function loadVelocityGeoTIFFGrid(tiff) {
+    const count = 'images' in tiff ? tiff.images.length : await tiff.getImageCount();
+    if (!Number.isSafeInteger(count) || count < 1)
+        throw new Error('Velocity GeoTIFF requires images');
+    const images = [[], [], []];
+    const names = ['east_velocity', 'north_velocity', 'up_velocity'];
+    for (let index = 0; index < count; index++) {
+        const raster = 'images' in tiff ? tiff.images[index] : undefined;
+        const image = raster ? undefined : await tiff.getImage(index);
+        const metadata = raster ? raster.metadata : await image.getGDALMetadata();
+        if (metadata?.['TYPE'] !== 'VELOCITY')
+            throw new Error('Velocity GeoTIFF requires TYPE=VELOCITY');
+        for (let component = 0; component < 3; component++) {
+            const band = raster?.bands.find(b => b.index === component);
+            if (raster && !band)
+                throw new Error('Velocity GeoTIFF requires original ENU bands 0, 1, 2');
+            const bandMetadata = raster ? band.metadata : await image.getGDALMetadata(component);
+            if (bandMetadata?.['DESCRIPTION'] !== names[component] ||
+                !['millimetres per year', 'mm/year'].includes(bandMetadata?.['UNITTYPE']))
+                throw new Error('Velocity GeoTIFF requires explicit ENU descriptions and millimetres per year');
+            // Reuse scalar geometry, scale/offset, nodata, ownership and hierarchy checks.
+            // This internal view presents native millimetres as generic scalar values;
+            // conversion to metres/year occurs only after all three components are sampled.
+            images[component].push({
+                getWidth: () => (raster ? raster.width : image.getWidth()),
+                getHeight: () => (raster ? raster.height : image.getHeight()),
+                getGeoKeys: () => (raster ? raster.geoKeys : image.getGeoKeys()),
+                getGDALMetadata: sample => sample === 0
+                    ? { ...bandMetadata, DESCRIPTION: 'geoid_undulation', UNITTYPE: 'metre' }
+                    : { ...metadata, TYPE: 'VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL' },
+                getGDALNoData: () => (raster ? raster.noData : (image.getGDALNoData?.() ?? null)),
+                fileDirectory: raster ? raster.fileDirectory : image.fileDirectory,
+                readRasters: async () => raster ? [band.data] : image.readRasters({ samples: [component], interleave: false })
+            });
+        }
+    }
+    const bands = await Promise.all(images.map(list => prepareScalarGeoTIFFGrids({
+        getImageCount: async () => list.length,
+        getImage: async (index) => list[index]
+    })));
+    const grids = bands[0].map((east, index) => velocityGridFromComponents(east, bands[1][index], bands[2][index], 0.001));
+    return Object.freeze({
+        sample(longitude, latitude, output) {
+            for (let i = grids.length - 1; i >= 0; i--)
+                if (grids[i].sample(longitude, latitude, output))
+                    return true;
+            return false;
+        }
+    });
+}
+//# sourceMappingURL=velocity-geotiff.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"velocity-geotiff.js","sourceRoot":"","sources":["../../../src/experimental/grids/velocity-geotiff.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,wLAAwL;AACxL,OAAO,EAAC,yBAAyB,EAAC,4BAAyB;AAC3D,OAAO,EAAC,0BAA0B,EAAC,sBAAmB;AAWtD;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,uBAAuB,CAC3C,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,kCAAkC,CAAC,CAAC;IACtD,MAAM,MAAM,GAAiC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAC1D,MAAM,KAAK,GAAG,CAAC,eAAe,EAAE,gBAAgB,EAAE,aAAa,CAAC,CAAC;IACjE,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC3C,MAAM,MAAM,GAAG,QAAQ,IAAI,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;QACjE,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,MAAO,IAA4B,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QACvF,MAAM,QAAQ,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,eAAe,EAAE,CAAC;QAC1E,IAAI,QAAQ,EAAE,CAAC,MAAM,CAAC,KAAK,UAAU;YACnC,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;QAC7D,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE,CAAC;YACnD,MAAM,IAAI,GAAG,MAAM,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC;YAC5D,IAAI,MAAM,IAAI,CAAC,IAAI;gBAAE,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;YAC7F,MAAM,YAAY,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,KAAK,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACrF,IACE,YAAY,EAAE,CAAC,aAAa,CAAC,KAAK,KAAK,CAAC,SAAS,CAAC;gBAClD,CAAC,CAAC,sBAAsB,EAAE,SAAS,CAAC,CAAC,QAAQ,CAAC,YAAY,EAAE,CAAC,UAAU,CAAW,CAAC;gBAEnF,MAAM,IAAI,KAAK,CACb,8EAA8E,CAC/E,CAAC;YACJ,+EAA+E;YAC/E,2EAA2E;YAC3E,gFAAgF;YAChF,MAAM,CAAC,SAAS,CAAC,CAAC,IAAI,CAAC;gBACrB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,QAAQ,EAAE,CAAC;gBAC1D,SAAS,EAAE,GAAG,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,SAAS,EAAE,CAAC;gBAC7D,UAAU,EAAE,GAAG,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC;gBAChE,eAAe,EAAE,MAAM,CAAC,EAAE,CACxB,MAAM,KAAK,CAAC;oBACV,CAAC,CAAC,EAAC,GAAG,YAAY,EAAE,WAAW,EAAE,kBAAkB,EAAE,QAAQ,EAAE,OAAO,EAAC;oBACvE,CAAC,CAAC,EAAC,GAAG,QAAQ,EAAE,IAAI,EAAE,wCAAwC,EAAC;gBACnE,aAAa,EAAE,GAAG,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa,EAAE,EAAE,IAAI,IAAI,CAAC,CAAC;gBACjF,aAAa,EAAE,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,CAAC,aAAa;gBAClE,WAAW,EAAE,KAAK,IAAI,EAAE,CACtB,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,WAAW,CAAC,EAAC,OAAO,EAAE,CAAC,SAAS,CAAC,EAAE,UAAU,EAAE,KAAK,EAAC,CAAC;aACtF,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,OAAO,CAAC,GAAG,CAC7B,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAChB,yBAAyB,CAAC;QACxB,aAAa,EAAE,KAAK,IAAI,EAAE,CAAC,IAAI,CAAC,MAAM;QACtC,QAAQ,EAAE,KAAK,EAAC,KAAK,EAAC,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC;KACrC,CAAC,CACH,CACF,CAAC;IACF,MAAM,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CACzC,0BAA0B,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,KAAK,CAAC,CAC1E,CAAC;IACF,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM;YAChC,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBACxC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC;oBAAE,OAAO,IAAI,CAAC;YAChE,OAAO,KAAK,CAAC;QACf,CAAC;KACF,CAAC,CAAC;AACL,CAAC"}

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

@@ -0,0 +1,21 @@
+import type { VerticalGridOptions } from "./vertical.js";
+import type { ProjectionPoint } from "../types.js";
+import type { VerticalGrid } from "./types.js";
+/** Prepared velocities in metres per decimal year at Greenwich lon/lat radians.
+ * sample writes ENU into output.x/y/z only on success. False means no coverage.
+ */
+export type VelocityGrid = {
+    sample(longitude: number, latitude: number, output: ProjectionPoint): boolean;
+};
+export type VelocityGridOptions = Omit<VerticalGridOptions, 'offsets'> & {
+    readonly units: 'm/year' | 'mm/year';
+    /** Owned snapshots, row-major south-to-north and west-to-east. */
+    readonly east: ArrayLike<number>;
+    readonly north: ArrayLike<number>;
+    readonly up: ArrayLike<number>;
+};
+/** Internal atomic three-band sampler: never mix components from different images. */
+export declare function velocityGridFromComponents(east: VerticalGrid, north: VerticalGrid, up: VerticalGrid, scale: number): VelocityGrid;
+/** No extrapolation; nodata in any contributing component makes the sample unavailable. */
+export declare function createVelocityGrid(options: VelocityGridOptions): VelocityGrid;
+//# sourceMappingURL=velocity.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"velocity.d.ts","sourceRoot":"","sources":["../../../src/experimental/grids/velocity.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,mBAAmB,EAAC,sBAAmB;AACpD,OAAO,KAAK,EAAC,eAAe,EAAC,oBAAiB;AAC9C,OAAO,KAAK,EAAC,YAAY,EAAC,mBAAgB;AAE1C;;GAEG;AACH,MAAM,MAAM,YAAY,GAAG;IACzB,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,eAAe,GAAG,OAAO,CAAC;CAC/E,CAAC;AACF,MAAM,MAAM,mBAAmB,GAAG,IAAI,CAAC,mBAAmB,EAAE,SAAS,CAAC,GAAG;IACvE,QAAQ,CAAC,KAAK,EAAE,QAAQ,GAAG,SAAS,CAAC;IACrC,kEAAkE;IAClE,QAAQ,CAAC,IAAI,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;IACjC,QAAQ,CAAC,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;IAClC,QAAQ,CAAC,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,CAAC;CAChC,CAAC;AACF,sFAAsF;AACtF,wBAAgB,0BAA0B,CACxC,IAAI,EAAE,YAAY,EAClB,KAAK,EAAE,YAAY,EACnB,EAAE,EAAE,YAAY,EAChB,KAAK,EAAE,MAAM,GACZ,YAAY,CAkBd;AACD,2FAA2F;AAC3F,wBAAgB,kBAAkB,CAAC,OAAO,EAAE,mBAAmB,GAAG,YAAY,CAS7E"}

package/dist/experimental/grids/velocity.js

@@ -0,0 +1,29 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original ENU grid preparation inspired by PROJ's velocity grid conventions. Bilinear interpolation reuses original math.gl vertical grids; no upstream fork.
+import { createVerticalGrid } from "./vertical.js";
+/** Internal atomic three-band sampler: never mix components from different images. */
+export function velocityGridFromComponents(east, north, up, scale) {
+    return Object.freeze({
+        sample(longitude, latitude, output) {
+            const e = east.getOffset(longitude, latitude), n = north.getOffset(longitude, latitude), u = up.getOffset(longitude, latitude);
+            if (e === undefined || n === undefined || u === undefined)
+                return false;
+            const x = e * scale, y = n * scale, z = u * scale;
+            if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))
+                throw new Error('Non-finite velocity interpolation');
+            output.x = x;
+            output.y = y;
+            output.z = z;
+            return true;
+        }
+    });
+}
+/** No extrapolation; nodata in any contributing component makes the sample unavailable. */
+export function createVelocityGrid(options) {
+    if (!['m/year', 'mm/year'].includes(options.units))
+        throw new Error('Explicit velocity units required');
+    return velocityGridFromComponents(createVerticalGrid({ ...options, offsets: options.east }), createVerticalGrid({ ...options, offsets: options.north }), createVerticalGrid({ ...options, offsets: options.up }), options.units === 'mm/year' ? 0.001 : 1);
+}
+//# sourceMappingURL=velocity.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"velocity.js","sourceRoot":"","sources":["../../../src/experimental/grids/velocity.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,iLAAiL;AACjL,OAAO,EAAC,kBAAkB,EAAC,sBAAmB;AAkB9C,sFAAsF;AACtF,MAAM,UAAU,0BAA0B,CACxC,IAAkB,EAClB,KAAmB,EACnB,EAAgB,EAChB,KAAa;IAEb,OAAO,MAAM,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,SAAiB,EAAE,QAAgB,EAAE,MAAuB;YACjE,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC,EAC3C,CAAC,GAAG,KAAK,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC,EACxC,CAAC,GAAG,EAAE,CAAC,SAAS,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;YACxC,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;gBAAE,OAAO,KAAK,CAAC;YACxE,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,EACjB,CAAC,GAAG,CAAC,GAAG,KAAK,EACb,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;YAChB,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,mCAAmC,CAAC,CAAC;YACvD,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,OAAO,IAAI,CAAC;QACd,CAAC;KACF,CAAC,CAAC;AACL,CAAC;AACD,2FAA2F;AAC3F,MAAM,UAAU,kBAAkB,CAAC,OAA4B;IAC7D,IAAI,CAAC,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;IACtD,OAAO,0BAA0B,CAC/B,kBAAkB,CAAC,EAAC,GAAG,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,IAAI,EAAC,CAAC,EACvD,kBAAkB,CAAC,EAAC,GAAG,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,KAAK,EAAC,CAAC,EACxD,kBAAkB,CAAC,EAAC,GAAG,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,EAAE,EAAC,CAAC,EACrD,OAAO,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CACxC,CAAC;AACJ,CAAC"}

package/dist/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/experimental/grids/vertical-geotiff.js
+var vertical_geotiff_exports = {};
+__export(vertical_geotiff_exports, {
+  loadVerticalGeoTIFFGrid: () => loadVerticalGeoTIFFGrid
+});
+module.exports = __toCommonJS(vertical_geotiff_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/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/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/experimental/grids/vertical-geotiff.cjs.map

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