@math.gl/geospatial 5.0.0-alpha.9 → 5.0.0-beta.2

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

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

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

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

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

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

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

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

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

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

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

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

package/dist/projection/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/projection/experimental/grids/velocity-geotiff.js
+var velocity_geotiff_exports = {};
+__export(velocity_geotiff_exports, {
+  loadVelocityGeoTIFFGrid: () => loadVelocityGeoTIFFGrid
+});
+module.exports = __toCommonJS(velocity_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/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/projection/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/projection/experimental/grids/velocity-geotiff.cjs.map

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

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

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

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

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