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

package/dist/terrain/contours.js

@@ -0,0 +1,133 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Contour the field's fixed-diagonal piecewise-linear surface, not its bilinear interpolant.
+ * Below vs at-or-above is the tie rule. Flat triangles and point-only touches emit no segment.
+ * Output is an indexed segment graph; no path stitching, simplification or polygon filling.
+ */
+export function createTerrainContours(field, levels, options = {}) {
+    const maxSamples = options.maxSamples ?? 1100000, maxLevels = options.maxLevels ?? 256;
+    const maxVertices = options.maxVertices ?? 2000000, maxSegments = options.maxSegments ?? 1000000;
+    const maxTriangleTests = options.maxTriangleTests ?? 10000000;
+    for (const limit of [maxSamples, maxLevels, maxVertices, maxSegments, maxTriangleTests])
+        if (!Number.isSafeInteger(limit) || limit < 0)
+            throw new Error('Invalid terrain contour capacity');
+    const count = field.width * field.height;
+    if (count > maxSamples ||
+        count > 0xffffffff ||
+        maxVertices > 0xffffffff ||
+        maxSegments > 0xffffffff)
+        throw new Error('Terrain contour sample or index capacity exceeded');
+    if (!Number.isSafeInteger(levels.length) || levels.length < 0 || levels.length > maxLevels)
+        throw new Error('Terrain contour level capacity exceeded');
+    for (let i = 0; i < levels.length; i++)
+        if (!Number.isFinite(levels[i]) || (i > 0 && levels[i] <= levels[i - 1]))
+            throw new Error('Terrain contour levels must be finite and strictly increasing');
+    if (levels.length * 2 * (field.width - 1) * (field.height - 1) > maxTriangleTests)
+        throw new Error('Terrain contour triangle-test budget exceeded');
+    const origin = options.origin ?? [field.origin[0], field.origin[1], 0];
+    if (origin.length !== 3 ||
+        !Number.isFinite(origin[0]) ||
+        !Number.isFinite(origin[1]) ||
+        !Number.isFinite(origin[2]))
+        throw new Error('Terrain contour origin must be finite XYZ');
+    const outputLevels = Float64Array.from(levels), positions = [], coordinates = [], segments = [], levelIndices = [];
+    const shift = field.registration === 'cell' ? 0.5 : 0;
+    for (let levelIndex = 0; levelIndex < outputLevels.length; levelIndex++) {
+        const level = outputLevels[levelIndex], vertices = new Map(), edges = new Set();
+        const height = (i) => field.getHeight(i % field.width, Math.floor(i / field.width));
+        const intersection = (a, b) => {
+            if (a > b)
+                [a, b] = [b, a];
+            const exact = height(a) === level ? a : height(b) === level ? b : undefined;
+            return exact === undefined
+                ? { a, b, key: `e${a}/${b}` }
+                : { a: exact, b: exact, key: `v${exact}` };
+        };
+        const vertex = ({ a, b, key }) => {
+            const existing = vertices.get(key);
+            if (existing !== undefined)
+                return existing;
+            if (positions.length / 3 >= maxVertices)
+                throw new Error('Terrain contour vertex capacity exceeded');
+            let fraction = 0;
+            if (a !== b) {
+                const ha = height(a), hb = height(b), delta = hb - ha, offset = level - ha;
+                if (Number.isFinite(delta) && Number.isFinite(offset))
+                    fraction = offset / delta;
+                else {
+                    const scale = Math.max(Math.abs(ha), Math.abs(hb), Math.abs(level));
+                    fraction = (level / scale - ha / scale) / (hb / scale - ha / scale);
+                }
+                if (!Number.isFinite(fraction) || fraction <= 0 || fraction >= 1)
+                    throw new Error('Terrain contour intersection loses precision');
+            }
+            const ax = a % field.width, ay = Math.floor(a / field.width);
+            const x = ax + ((b % field.width) - ax) * fraction, y = ay + (Math.floor(b / field.width) - ay) * fraction;
+            const px = field.origin[0] -
+                origin[0] +
+                (x + shift) * field.columnStep[0] +
+                (y + shift) * field.rowStep[0];
+            const py = field.origin[1] -
+                origin[1] +
+                (x + shift) * field.columnStep[1] +
+                (y + shift) * field.rowStep[1];
+            const pz = level - origin[2];
+            if (!Number.isFinite(px) || !Number.isFinite(py) || !Number.isFinite(pz))
+                throw new Error('Terrain contour coordinate overflow');
+            const index = positions.length / 3;
+            positions.push(px, py, pz);
+            coordinates.push(x, y);
+            vertices.set(key, index);
+            return index;
+        };
+        const triangle = (a, b, c) => {
+            const ha = height(a), hb = height(b), hc = height(c);
+            if (ha === undefined || hb === undefined || hc === undefined)
+                return;
+            const hits = [];
+            if (ha < level !== hb < level)
+                hits.push(intersection(a, b));
+            if (hb < level !== hc < level)
+                hits.push(intersection(b, c));
+            if (hc < level !== ha < level)
+                hits.push(intersection(c, a));
+            if (hits.length !== 2 || hits[0].key === hits[1].key)
+                return;
+            const va = vertex(hits[0]), vb = vertex(hits[1]);
+            const key = va < vb ? `${va}/${vb}` : `${vb}/${va}`;
+            if (edges.has(key))
+                return;
+            if (segments.length / 2 >= maxSegments)
+                throw new Error('Terrain contour segment capacity exceeded');
+            if (positions[3 * va] === positions[3 * vb] &&
+                positions[3 * va + 1] === positions[3 * vb + 1])
+                throw new Error('Terrain contour segment collapses; use a local origin');
+            edges.add(key);
+            segments.push(va, vb);
+            levelIndices.push(levelIndex);
+        };
+        for (let row = 0; row < field.height - 1; row++)
+            for (let column = 0; column < field.width - 1; column++) {
+                const a = row * field.width + column, b = a + 1, c = a + field.width, d = c + 1;
+                if (field.diagonal === '00-11') {
+                    triangle(a, b, d);
+                    triangle(a, d, c);
+                }
+                else {
+                    triangle(a, b, c);
+                    triangle(b, d, c);
+                }
+            }
+    }
+    return {
+        positions: Float64Array.from(positions),
+        gridCoordinates: Float64Array.from(coordinates),
+        segments: Uint32Array.from(segments),
+        levelIndices: Uint32Array.from(levelIndices),
+        levels: outputLevels,
+        origin: Object.freeze([...origin]),
+        heightReference: field.heightReference
+    };
+}
+//# sourceMappingURL=contours.js.map

package/dist/terrain/contours.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"contours.js","sourceRoot":"","sources":["../../src/terrain/contours.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AA2B5D;;;GAGG;AACH,MAAM,UAAU,qBAAqB,CACnC,KAAkB,EAClB,MAAyB,EACzB,UAAiC,EAAE;IAEnC,MAAM,UAAU,GAAG,OAAO,CAAC,UAAU,IAAI,OAAS,EAChD,SAAS,GAAG,OAAO,CAAC,SAAS,IAAI,GAAG,CAAC;IACvC,MAAM,WAAW,GAAG,OAAO,CAAC,WAAW,IAAI,OAAS,EAClD,WAAW,GAAG,OAAO,CAAC,WAAW,IAAI,OAAS,CAAC;IACjD,MAAM,gBAAgB,GAAG,OAAO,CAAC,gBAAgB,IAAI,QAAU,CAAC;IAChE,KAAK,MAAM,KAAK,IAAI,CAAC,UAAU,EAAE,SAAS,EAAE,WAAW,EAAE,WAAW,EAAE,gBAAgB,CAAC;QACrF,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC;YAC3C,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;IACxD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC;IACzC,IACE,KAAK,GAAG,UAAU;QAClB,KAAK,GAAG,UAAU;QAClB,WAAW,GAAG,UAAU;QACxB,WAAW,GAAG,UAAU;QAExB,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,CAAC,MAAM,GAAG,SAAS;QACxF,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,EAAE;QACpC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACtE,MAAM,IAAI,KAAK,CAAC,+DAA+D,CAAC,CAAC;IACrF,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,gBAAgB;QAC/E,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACvE,IACE,MAAM,CAAC,MAAM,KAAK,CAAC;QACnB,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAE3B,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,MAAM,YAAY,GAAG,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,EAC5C,SAAS,GAAa,EAAE,EACxB,WAAW,GAAa,EAAE,EAC1B,QAAQ,GAAa,EAAE,EACvB,YAAY,GAAa,EAAE,CAAC;IAC9B,MAAM,KAAK,GAAG,KAAK,CAAC,YAAY,KAAK,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACtD,KAAK,IAAI,UAAU,GAAG,CAAC,EAAE,UAAU,GAAG,YAAY,CAAC,MAAM,EAAE,UAAU,EAAE,EAAE,CAAC;QACxE,MAAM,KAAK,GAAG,YAAY,CAAC,UAAU,CAAC,EACpC,QAAQ,GAAG,IAAI,GAAG,EAAkB,EACpC,KAAK,GAAG,IAAI,GAAG,EAAU,CAAC;QAC5B,MAAM,MAAM,GAAG,CAAC,CAAS,EAAsB,EAAE,CAC/C,KAAK,CAAC,SAAS,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;QAChE,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,CAAS,EAAuC,EAAE;YACjF,IAAI,CAAC,GAAG,CAAC;gBAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;YAC5E,OAAO,KAAK,KAAK,SAAS;gBACxB,CAAC,CAAC,EAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,EAAE,EAAC;gBAC3B,CAAC,CAAC,EAAC,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,KAAK,EAAE,EAAC,CAAC;QAC7C,CAAC,CAAC;QACF,MAAM,MAAM,GAAG,CAAC,EAAC,CAAC,EAAE,CAAC,EAAE,GAAG,EAAsC,EAAU,EAAE;YAC1E,MAAM,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACnC,IAAI,QAAQ,KAAK,SAAS;gBAAE,OAAO,QAAQ,CAAC;YAC5C,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC,IAAI,WAAW;gBACrC,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;YAC9D,IAAI,QAAQ,GAAG,CAAC,CAAC;YACjB,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;gBACZ,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAE,EACnB,EAAE,GAAG,MAAM,CAAC,CAAC,CAAE,EACf,KAAK,GAAG,EAAE,GAAG,EAAE,EACf,MAAM,GAAG,KAAK,GAAG,EAAE,CAAC;gBACtB,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;oBAAE,QAAQ,GAAG,MAAM,GAAG,KAAK,CAAC;qBAC5E,CAAC;oBACJ,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;oBACpE,QAAQ,GAAG,CAAC,KAAK,GAAG,KAAK,GAAG,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;gBACtE,CAAC;gBACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,IAAI,QAAQ,IAAI,CAAC;oBAC9D,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;YACpE,CAAC;YACD,MAAM,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,EACxB,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC;YACnC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,QAAQ,EAChD,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC,GAAG,QAAQ,CAAC;YACzD,MAAM,EAAE,GACN,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;gBACf,MAAM,CAAC,CAAC,CAAC;gBACT,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC;gBACjC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YACjC,MAAM,EAAE,GACN,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC;gBACf,MAAM,CAAC,CAAC,CAAC;gBACT,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC;gBACjC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YACjC,MAAM,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YAC7B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;gBACtE,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;YACzD,MAAM,KAAK,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YACnC,SAAS,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YAC3B,WAAW,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACvB,QAAQ,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;YACzB,OAAO,KAAK,CAAC;QACf,CAAC,CAAC;QACF,MAAM,QAAQ,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAQ,EAAE;YACzD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,EAClB,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,EACd,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;YACjB,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS;gBAAE,OAAO;YACrE,MAAM,IAAI,GAAsC,EAAE,CAAC;YACnD,IAAI,EAAE,GAAG,KAAK,KAAK,EAAE,GAAG,KAAK;gBAAE,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7D,IAAI,EAAE,GAAG,KAAK,KAAK,EAAE,GAAG,KAAK;gBAAE,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7D,IAAI,EAAE,GAAG,KAAK,KAAK,EAAE,GAAG,KAAK;gBAAE,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7D,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG;gBAAE,OAAO;YAC7D,MAAM,EAAE,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EACxB,EAAE,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;YACpD,IAAI,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC;gBAAE,OAAO;YAC3B,IAAI,QAAQ,CAAC,MAAM,GAAG,CAAC,IAAI,WAAW;gBACpC,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;YAC/D,IACE,SAAS,CAAC,CAAC,GAAG,EAAE,CAAC,KAAK,SAAS,CAAC,CAAC,GAAG,EAAE,CAAC;gBACvC,SAAS,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;gBAE/C,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;YAC3E,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACf,QAAQ,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACtB,YAAY,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;QAChC,CAAC,CAAC;QACF,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE;YAC7C,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;gBACxD,MAAM,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,KAAK,GAAG,MAAM,EAClC,CAAC,GAAG,CAAC,GAAG,CAAC,EACT,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,KAAK,EACnB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;gBACZ,IAAI,KAAK,CAAC,QAAQ,KAAK,OAAO,EAAE,CAAC;oBAC/B,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;oBAClB,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACpB,CAAC;qBAAM,CAAC;oBACN,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;oBAClB,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACpB,CAAC;YACH,CAAC;IACL,CAAC;IACD,OAAO;QACL,SAAS,EAAE,YAAY,CAAC,IAAI,CAAC,SAAS,CAAC;QACvC,eAAe,EAAE,YAAY,CAAC,IAAI,CAAC,WAAW,CAAC;QAC/C,QAAQ,EAAE,WAAW,CAAC,IAAI,CAAC,QAAQ,CAAC;QACpC,YAAY,EAAE,WAAW,CAAC,IAAI,CAAC,YAAY,CAAC;QAC5C,MAAM,EAAE,YAAY;QACpB,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAA6B,CAAC;QAC9D,eAAe,EAAE,KAAK,CAAC,eAAe;KACvC,CAAC;AACJ,CAAC"}

package/dist/terrain/delatin.d.ts

@@ -0,0 +1,12 @@
+import type { HeightField } from "./height-field.js";
+import type { AdaptiveTerrainMesh, AdaptiveTerrainMeshOptions } from "./adaptive-mesh-types.js";
+export type { AdaptiveTerrainMesh } from "./adaptive-mesh-types.js";
+export type DelatinTerrainMeshOptions = AdaptiveTerrainMeshOptions & {
+    /** Maximum point insertions; defaults to source sample count minus the four corners. */
+    maxRefinements?: number;
+};
+/** Optional greedy Delaunay terrain meshing for complete rectangular grids.
+ * Backend ownership lives here so loaders can retain decoding, padding, UVs, skirts and schema work.
+ */
+export declare function createDelatinTerrainMesh(field: HeightField, options: DelatinTerrainMeshOptions): AdaptiveTerrainMesh;
+//# sourceMappingURL=delatin.d.ts.map

package/dist/terrain/delatin.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"delatin.d.ts","sourceRoot":"","sources":["../../src/terrain/delatin.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,WAAW,EAAC,0BAAuB;AAIhD,OAAO,KAAK,EAAC,mBAAmB,EAAE,0BAA0B,EAAC,iCAA8B;AAC3F,YAAY,EAAC,mBAAmB,EAAC,iCAA8B;AAE/D,MAAM,MAAM,yBAAyB,GAAG,0BAA0B,GAAG;IACnE,wFAAwF;IACxF,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB,CAAC;AAEF;;GAEG;AACH,wBAAgB,wBAAwB,CACtC,KAAK,EAAE,WAAW,EAClB,OAAO,EAAE,yBAAyB,GACjC,mBAAmB,CA0DrB"}

package/dist/terrain/delatin.js

@@ -0,0 +1,69 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import Delatin from "./vendor/delatin.js";
+import { assembleTerrainMesh } from "./mesh-data.js";
+import { createTerrainMesh } from "./mesh.js";
+import { measureTerrainMeshError } from "./mesh-error.js";
+/** Optional greedy Delaunay terrain meshing for complete rectangular grids.
+ * Backend ownership lives here so loaders can retain decoding, padding, UVs, skirts and schema work.
+ */
+export function createDelatinTerrainMesh(field, options) {
+    const count = field.width * field.height;
+    const maxSamples = options.maxSamples ?? 1100000;
+    const maxTriangles = options.maxTriangles ?? 2 * (field.width - 1) * (field.height - 1);
+    const maxRefinements = options.maxRefinements ?? count - 4;
+    if (!Number.isFinite(options.maxError) || options.maxError < 0)
+        throw new Error('Delatin terrain maxError must be finite and nonnegative');
+    for (const limit of [maxSamples, maxTriangles, maxRefinements])
+        if (!Number.isSafeInteger(limit) || limit < 0)
+            throw new Error('Invalid Delatin terrain capacity');
+    if (count > maxSamples || count > 0xffffffff)
+        throw new Error('Delatin terrain sample capacity exceeded');
+    if (field.validCount !== count)
+        throw new Error('Delatin terrain requires complete source data; use a regular mesh for holes');
+    if (maxTriangles < 2)
+        throw new Error('Delatin terrain triangle capacity exceeded');
+    if (options.fullEdges !== undefined)
+        throw new Error('Delatin does not support fullEdges; use MARTINI for forced boundary refinement');
+    const { heights } = field.copyData();
+    const tin = new Delatin(heights, field.width, field.height);
+    let threshold = options.maxError, refinements = 0;
+    for (let pass = 0; pass <= 8; pass++) {
+        if (pass < 8) {
+            while (true) {
+                const error = tin.getMaxError();
+                if (!Number.isFinite(error))
+                    throw new Error('Delatin error arithmetic overflow');
+                if (error <= threshold)
+                    break;
+                if (refinements >= maxRefinements)
+                    throw new Error('Delatin terrain refinement budget exceeded');
+                tin.refine();
+                refinements++;
+                if (tin.triangles.length / 3 > maxTriangles)
+                    throw new Error('Delatin terrain triangle capacity exceeded');
+            }
+        }
+        if (pass === 8 && 2 * (field.width - 1) * (field.height - 1) > maxTriangles)
+            throw new Error('Delatin terrain fallback triangle capacity exceeded');
+        const indices = Uint32Array.from(tin.triangles);
+        for (let i = 0; i < indices.length; i += 3)
+            [indices[i + 1], indices[i + 2]] = [indices[i + 2], indices[i + 1]];
+        const mesh = pass === 8
+            ? createTerrainMesh(field, options)
+            : assembleTerrainMesh(field, Uint32Array.from(tin.coords), indices, options);
+        const error = measureTerrainMeshError(field, mesh, { maxSampleTests: options.maxSampleTests });
+        if (!error.uncoveredSampleCount && error.maxSampleError <= options.maxError)
+            return {
+                ...mesh,
+                error,
+                requestedMaxError: options.maxError,
+                refinementPasses: pass + 1,
+                tessellator: pass === 8 ? 'regular' : 'delatin'
+            };
+        threshold /= 2;
+    }
+    throw new Error('Delatin terrain could not meet the requested sample error');
+}
+//# sourceMappingURL=delatin.js.map

package/dist/terrain/delatin.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"delatin.js","sourceRoot":"","sources":["../../src/terrain/delatin.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,OAAO,4BAAyB;AAEvC,OAAO,EAAC,mBAAmB,EAAC,uBAAoB;AAChD,OAAO,EAAC,iBAAiB,EAAC,kBAAe;AACzC,OAAO,EAAC,uBAAuB,EAAC,wBAAqB;AASrD;;GAEG;AACH,MAAM,UAAU,wBAAwB,CACtC,KAAkB,EAClB,OAAkC;IAElC,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC;IACzC,MAAM,UAAU,GAAG,OAAO,CAAC,UAAU,IAAI,OAAS,CAAC;IACnD,MAAM,YAAY,GAAG,OAAO,CAAC,YAAY,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACxF,MAAM,cAAc,GAAG,OAAO,CAAC,cAAc,IAAI,KAAK,GAAG,CAAC,CAAC;IAC3D,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,QAAQ,CAAC,IAAI,OAAO,CAAC,QAAQ,GAAG,CAAC;QAC5D,MAAM,IAAI,KAAK,CAAC,yDAAyD,CAAC,CAAC;IAC7E,KAAK,MAAM,KAAK,IAAI,CAAC,UAAU,EAAE,YAAY,EAAE,cAAc,CAAC;QAC5D,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC;YAC3C,MAAM,IAAI,KAAK,CAAC,kCAAkC,CAAC,CAAC;IACxD,IAAI,KAAK,GAAG,UAAU,IAAI,KAAK,GAAG,UAAU;QAC1C,MAAM,IAAI,KAAK,CAAC,0CAA0C,CAAC,CAAC;IAC9D,IAAI,KAAK,CAAC,UAAU,KAAK,KAAK;QAC5B,MAAM,IAAI,KAAK,CAAC,6EAA6E,CAAC,CAAC;IACjG,IAAI,YAAY,GAAG,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;IACpF,IAAK,OAAiC,CAAC,SAAS,KAAK,SAAS;QAC5D,MAAM,IAAI,KAAK,CACb,gFAAgF,CACjF,CAAC;IACJ,MAAM,EAAC,OAAO,EAAC,GAAG,KAAK,CAAC,QAAQ,EAAE,CAAC;IACnC,MAAM,GAAG,GAAG,IAAI,OAAO,CAAC,OAAO,EAAE,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;IAC5D,IAAI,SAAS,GAAG,OAAO,CAAC,QAAQ,EAC9B,WAAW,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QACrC,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;YACb,OAAO,IAAI,EAAE,CAAC;gBACZ,MAAM,KAAK,GAAG,GAAG,CAAC,WAAW,EAAE,CAAC;gBAChC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;oBAAE,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;gBAClF,IAAI,KAAK,IAAI,SAAS;oBAAE,MAAM;gBAC9B,IAAI,WAAW,IAAI,cAAc;oBAC/B,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;gBAChE,GAAG,CAAC,MAAM,EAAE,CAAC;gBACb,WAAW,EAAE,CAAC;gBACd,IAAI,GAAG,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,GAAG,YAAY;oBACzC,MAAM,IAAI,KAAK,CAAC,4CAA4C,CAAC,CAAC;YAClE,CAAC;QACH,CAAC;QACD,IAAI,IAAI,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,YAAY;YACzE,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,CAAC;QACzE,MAAM,OAAO,GAAG,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QAChD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC;YACxC,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,IAAI,GACR,IAAI,KAAK,CAAC;YACR,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,OAAO,CAAC;YACnC,CAAC,CAAC,mBAAmB,CAAC,KAAK,EAAE,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QACjF,MAAM,KAAK,GAAG,uBAAuB,CAAC,KAAK,EAAE,IAAI,EAAE,EAAC,cAAc,EAAE,OAAO,CAAC,cAAc,EAAC,CAAC,CAAC;QAC7F,IAAI,CAAC,KAAK,CAAC,oBAAoB,IAAI,KAAK,CAAC,cAAc,IAAI,OAAO,CAAC,QAAQ;YACzE,OAAO;gBACL,GAAG,IAAI;gBACP,KAAK;gBACL,iBAAiB,EAAE,OAAO,CAAC,QAAQ;gBACnC,gBAAgB,EAAE,IAAI,GAAG,CAAC;gBAC1B,WAAW,EAAE,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS;aAChD,CAAC;QACJ,SAAS,IAAI,CAAC,CAAC;IACjB,CAAC;IACD,MAAM,IAAI,KAAK,CAAC,2DAA2D,CAAC,CAAC;AAC/E,CAAC"}

package/dist/terrain/derivatives.d.ts

@@ -0,0 +1,23 @@
+import type { HeightField } from "./height-field.js";
+export type TerrainGradientMethod = 'horn' | 'zevenbergen-thorne';
+export type TerrainGradientGrid = {
+    width: number;
+    height: number;
+    /** Packed world-frame dz/dx, dz/dy; zero where status is not VALID. */
+    gradients: Float64Array;
+    statuses: Uint8Array;
+    method: TerrainGradientMethod;
+    origin: readonly [number, number];
+    columnStep: readonly [number, number];
+    rowStep: readonly [number, number];
+    registration: HeightField['registration'];
+    heightReference: string;
+};
+/** Raster stencil derivatives, distinct from derivatives of an interpolation patch.
+ * Both methods require all nine samples; no extrapolation or nodata filling occurs.
+ */
+export declare function createTerrainGradientGrid(field: HeightField, options?: {
+    method?: TerrainGradientMethod;
+    maxSamples?: number;
+}): TerrainGradientGrid;
+//# sourceMappingURL=derivatives.d.ts.map

package/dist/terrain/derivatives.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"derivatives.d.ts","sourceRoot":"","sources":["../../src/terrain/derivatives.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,WAAW,EAAC,0BAAuB;AAGhD,MAAM,MAAM,qBAAqB,GAAG,MAAM,GAAG,oBAAoB,CAAC;AAClE,MAAM,MAAM,mBAAmB,GAAG;IAChC,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IACf,uEAAuE;IACvE,SAAS,EAAE,YAAY,CAAC;IACxB,QAAQ,EAAE,UAAU,CAAC;IACrB,MAAM,EAAE,qBAAqB,CAAC;IAC9B,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAClC,UAAU,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACtC,OAAO,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACnC,YAAY,EAAE,WAAW,CAAC,cAAc,CAAC,CAAC;IAC1C,eAAe,EAAE,MAAM,CAAC;CACzB,CAAC;AAEF;;GAEG;AACH,wBAAgB,yBAAyB,CACvC,KAAK,EAAE,WAAW,EAClB,OAAO,GAAE;IAAC,MAAM,CAAC,EAAE,qBAAqB,CAAC;IAAC,UAAU,CAAC,EAAE,MAAM,CAAA;CAAM,GAClE,mBAAmB,CAkErB"}

package/dist/terrain/derivatives.js

@@ -0,0 +1,71 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { TerrainSampleStatus as Status } from "./types.js";
+/** Raster stencil derivatives, distinct from derivatives of an interpolation patch.
+ * Both methods require all nine samples; no extrapolation or nodata filling occurs.
+ */
+export function createTerrainGradientGrid(field, options = {}) {
+    const method = options.method ?? 'horn';
+    if (method !== 'horn' && method !== 'zevenbergen-thorne')
+        throw new Error('Unknown terrain gradient method');
+    const count = field.width * field.height, limit = options.maxSamples ?? 1100000;
+    if (!Number.isSafeInteger(limit) || limit < 0 || count > limit)
+        throw new Error('Terrain gradient sample capacity exceeded');
+    const gradients = new Float64Array(count * 2), statuses = new Uint8Array(count).fill(Status.OUTSIDE_DOMAIN);
+    const [a, c] = field.columnStep, [b, d] = field.rowStep;
+    const determinant = a * d - b * c;
+    const ix = d / determinant, iy = -c / determinant, jx = -b / determinant, jy = a / determinant;
+    const window = new Float64Array(9);
+    for (let row = 1; row < field.height - 1; row++)
+        for (let column = 1; column < field.width - 1; column++) {
+            const index = row * field.width + column;
+            let valid = true;
+            for (let y = -1; y <= 1; y++)
+                for (let x = -1; x <= 1; x++) {
+                    const value = field.getHeight(column + x, row + y);
+                    if (value === undefined)
+                        valid = false;
+                    window[(y + 1) * 3 + x + 1] = value ?? 0;
+                }
+            if (!valid) {
+                statuses[index] = Status.NO_DATA;
+                continue;
+            }
+            const dc = method === 'horn'
+                ? difference(window[2], window[0], 8) +
+                    difference(window[5], window[3], 4) +
+                    difference(window[8], window[6], 8)
+                : difference(window[5], window[3], 2);
+            const dr = method === 'horn'
+                ? difference(window[6], window[0], 8) +
+                    difference(window[7], window[1], 4) +
+                    difference(window[8], window[2], 8)
+                : difference(window[7], window[1], 2);
+            const gx = dc * ix + dr * iy, gy = dc * jx + dr * jy;
+            if (!Number.isFinite(gx) || !Number.isFinite(gy)) {
+                statuses[index] = Status.NUMERIC_ERROR;
+                continue;
+            }
+            gradients[2 * index] = gx;
+            gradients[2 * index + 1] = gy;
+            statuses[index] = Status.VALID;
+        }
+    return {
+        width: field.width,
+        height: field.height,
+        gradients,
+        statuses,
+        method,
+        origin: field.origin,
+        columnStep: field.columnStep,
+        rowStep: field.rowStep,
+        registration: field.registration,
+        heightReference: field.heightReference
+    };
+}
+function difference(a, b, divisor) {
+    const delta = a - b;
+    return Number.isFinite(delta) ? delta / divisor : a / divisor - b / divisor;
+}
+//# sourceMappingURL=derivatives.js.map

package/dist/terrain/derivatives.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"derivatives.js","sourceRoot":"","sources":["../../src/terrain/derivatives.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAG5D,OAAO,EAAC,mBAAmB,IAAI,MAAM,EAAC,mBAAgB;AAiBtD;;GAEG;AACH,MAAM,UAAU,yBAAyB,CACvC,KAAkB,EAClB,UAAiE,EAAE;IAEnE,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,MAAM,CAAC;IACxC,IAAI,MAAM,KAAK,MAAM,IAAI,MAAM,KAAK,oBAAoB;QACtD,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;IACrD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,MAAM,EACtC,KAAK,GAAG,OAAO,CAAC,UAAU,IAAI,OAAS,CAAC;IAC1C,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC,IAAI,KAAK,GAAG,KAAK;QAC5D,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;IAC/D,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,CAAC,CAAC,EAC3C,QAAQ,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC;IAC/D,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,UAAU,EAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC;IACzB,MAAM,WAAW,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClC,MAAM,EAAE,GAAG,CAAC,GAAG,WAAW,EACxB,EAAE,GAAG,CAAC,CAAC,GAAG,WAAW,EACrB,EAAE,GAAG,CAAC,CAAC,GAAG,WAAW,EACrB,EAAE,GAAG,CAAC,GAAG,WAAW,CAAC;IACvB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IACnC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE;QAC7C,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC;YACxD,MAAM,KAAK,GAAG,GAAG,GAAG,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC;YACzC,IAAI,KAAK,GAAG,IAAI,CAAC;YACjB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC7B,MAAM,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;oBACnD,IAAI,KAAK,KAAK,SAAS;wBAAE,KAAK,GAAG,KAAK,CAAC;oBACvC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,IAAI,CAAC,CAAC;gBAC3C,CAAC;YACH,IAAI,CAAC,KAAK,EAAE,CAAC;gBACX,QAAQ,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC;gBACjC,SAAS;YACX,CAAC;YACD,MAAM,EAAE,GACN,MAAM,KAAK,MAAM;gBACf,CAAC,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;oBACnC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;oBACnC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACrC,CAAC,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC1C,MAAM,EAAE,GACN,MAAM,KAAK,MAAM;gBACf,CAAC,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;oBACnC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;oBACnC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACrC,CAAC,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC1C,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAC1B,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YACzB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;gBACjD,QAAQ,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,aAAa,CAAC;gBACvC,SAAS;YACX,CAAC;YACD,SAAS,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC;YAC1B,SAAS,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YAC9B,QAAQ,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC;QACjC,CAAC;IACH,OAAO;QACL,KAAK,EAAE,KAAK,CAAC,KAAK;QAClB,MAAM,EAAE,KAAK,CAAC,MAAM;QACpB,SAAS;QACT,QAAQ;QACR,MAAM;QACN,MAAM,EAAE,KAAK,CAAC,MAAM;QACpB,UAAU,EAAE,KAAK,CAAC,UAAU;QAC5B,OAAO,EAAE,KAAK,CAAC,OAAO;QACtB,YAAY,EAAE,KAAK,CAAC,YAAY;QAChC,eAAe,EAAE,KAAK,CAAC,eAAe;KACvC,CAAC;AACJ,CAAC;AAED,SAAS,UAAU,CAAC,CAAS,EAAE,CAAS,EAAE,OAAe;IACvD,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;IACpB,OAAO,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO,CAAC;AAC9E,CAAC"}

package/dist/terrain/geographic.d.ts

@@ -0,0 +1,21 @@
+import { type HeightFieldProps, type TerrainSample, type TerrainSampleOptions } from "./types.js";
+export type GeographicHeightFieldProps = Omit<HeightFieldProps, 'origin' | 'columnStep' | 'rowStep'> & {
+    /** Unwrapped west, south, east, north in degrees. East must exceed west. */
+    bounds: readonly [number, number, number, number];
+    projection?: 'geographic' | 'web-mercator';
+    /** Raster row direction. Defaults to north to south. */
+    rowDirection?: 'north-to-south' | 'south-to-north';
+};
+/** Longitude/latitude queries with WGS84 surface-metric east/north derivatives. */
+export declare class GeographicHeightField {
+    readonly bounds: readonly [number, number, number, number];
+    readonly projection: 'geographic' | 'web-mercator';
+    readonly heightReference: string;
+    private readonly field;
+    constructor(props: GeographicHeightFieldProps);
+    /** Longitude must use the bounds' unwrapped interval; no implicit wrapping or clamping. */
+    sample(longitude: number, latitude: number, options?: TerrainSampleOptions, result?: TerrainSample): TerrainSample;
+    private projectX;
+    private projectY;
+}
+//# sourceMappingURL=geographic.d.ts.map

package/dist/terrain/geographic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geographic.d.ts","sourceRoot":"","sources":["../../src/terrain/geographic.ts"],"names":[],"mappings":"AAKA,OAAO,EAEL,KAAK,gBAAgB,EACrB,KAAK,aAAa,EAClB,KAAK,oBAAoB,EAC1B,mBAAgB;AAOjB,MAAM,MAAM,0BAA0B,GAAG,IAAI,CAC3C,gBAAgB,EAChB,QAAQ,GAAG,YAAY,GAAG,SAAS,CACpC,GAAG;IACF,4EAA4E;IAC5E,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAClD,UAAU,CAAC,EAAE,YAAY,GAAG,cAAc,CAAC;IAC3C,wDAAwD;IACxD,YAAY,CAAC,EAAE,gBAAgB,GAAG,gBAAgB,CAAC;CACpD,CAAC;AAEF,mFAAmF;AACnF,qBAAa,qBAAqB;IAChC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3D,QAAQ,CAAC,UAAU,EAAE,YAAY,GAAG,cAAc,CAAC;IACnD,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IACjC,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAc;gBAExB,KAAK,EAAE,0BAA0B;IAuC7C,2FAA2F;IAC3F,MAAM,CACJ,SAAS,EAAE,MAAM,EACjB,QAAQ,EAAE,MAAM,EAChB,OAAO,CAAC,EAAE,oBAAoB,EAC9B,MAAM,CAAC,EAAE,aAAa,GACrB,aAAa;IAyChB,OAAO,CAAC,QAAQ;IAIhB,OAAO,CAAC,QAAQ;CAKjB"}

package/dist/terrain/geographic.js

@@ -0,0 +1,91 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { HeightField } from "./height-field.js";
+import { TerrainSampleStatus as Status } from "./types.js";
+const RADIANS = Math.PI / 180;
+const RADIUS = 6378137;
+const ECCENTRICITY_SQUARED = 6.6943799901413165e-3;
+const MAX_MERCATOR_LATITUDE = Math.atan(Math.sinh(Math.PI)) / RADIANS;
+/** Longitude/latitude queries with WGS84 surface-metric east/north derivatives. */
+export class GeographicHeightField {
+    constructor(props) {
+        const { bounds } = props;
+        this.projection = props.projection ?? 'geographic';
+        this.heightReference = props.heightReference ?? 'unknown';
+        const limit = this.projection === 'web-mercator' ? MAX_MERCATOR_LATITUDE : 90;
+        if (bounds.length !== 4 ||
+            !bounds.every(Number.isFinite) ||
+            bounds[2] <= bounds[0] ||
+            bounds[2] - bounds[0] > 360 ||
+            bounds[3] <= bounds[1] ||
+            bounds[1] < -limit ||
+            bounds[3] > limit ||
+            (this.projection === 'geographic' && (bounds[1] === -90 || bounds[3] === 90))) {
+            throw new Error('Terrain bounds must be increasing, span at most 360 degrees and exclude the projection limits');
+        }
+        if (this.projection !== 'geographic' && this.projection !== 'web-mercator')
+            throw new Error('Unknown terrain projection');
+        const direction = props.rowDirection ?? 'north-to-south';
+        if (direction !== 'north-to-south' && direction !== 'south-to-north')
+            throw new Error('Unknown terrain row direction');
+        this.bounds = Object.freeze([...bounds]);
+        const west = this.projectX(bounds[0]);
+        const east = this.projectX(bounds[2]);
+        const south = this.projectY(bounds[1]);
+        const north = this.projectY(bounds[3]);
+        const cell = props.registration === 'cell';
+        const descending = direction === 'north-to-south';
+        this.field = new HeightField({
+            ...props,
+            origin: [west, descending ? north : south],
+            columnStep: [(east - west) / (props.width - (cell ? 0 : 1)), 0],
+            rowStep: [0, (descending ? south - north : north - south) / (props.height - (cell ? 0 : 1))]
+        });
+    }
+    /** Longitude must use the bounds' unwrapped interval; no implicit wrapping or clamping. */
+    sample(longitude, latitude, options, result) {
+        const valid = Number.isFinite(longitude) && Number.isFinite(latitude);
+        const inside = valid &&
+            longitude >= this.bounds[0] &&
+            longitude <= this.bounds[2] &&
+            latitude >= this.bounds[1] &&
+            latitude <= this.bounds[3];
+        const sample = this.field.sample(inside ? this.projectX(longitude) : NaN, inside ? this.projectY(latitude) : NaN, options, result);
+        if (!inside) {
+            sample.status = sample.gradientStatus = valid ? Status.OUTSIDE_DOMAIN : Status.INVALID_INPUT;
+            return sample;
+        }
+        if (sample.gradientStatus !== Status.VALID)
+            return sample;
+        const phi = latitude * RADIANS;
+        const denominator = 1 - ECCENTRICITY_SQUARED * Math.sin(phi) ** 2;
+        const n = RADIUS / Math.sqrt(denominator);
+        const m = (RADIUS * (1 - ECCENTRICITY_SQUARED)) / denominator ** 1.5;
+        const mercator = this.projection === 'web-mercator';
+        const gx = (sample.gradientX * (mercator ? RADIUS : 1 / RADIANS)) / (n * Math.cos(phi));
+        const gy = (sample.gradientY * (mercator ? RADIUS / Math.cos(phi) : 1 / RADIANS)) / m;
+        if (!Number.isFinite(gx) || !Number.isFinite(gy)) {
+            sample.gradientStatus = Status.NUMERIC_ERROR;
+            sample.gradientX = sample.gradientY = sample.normalX = sample.normalY = sample.normalZ = 0;
+            return sample;
+        }
+        sample.gradientX = gx;
+        sample.gradientY = gy;
+        const scale = Math.max(1, Math.abs(gx), Math.abs(gy));
+        const length = Math.hypot(gx / scale, gy / scale, 1 / scale);
+        sample.normalX = -gx / scale / length;
+        sample.normalY = -gy / scale / length;
+        sample.normalZ = 1 / scale / length;
+        return sample;
+    }
+    projectX(longitude) {
+        return this.projection === 'web-mercator' ? RADIUS * longitude * RADIANS : longitude;
+    }
+    projectY(latitude) {
+        return this.projection === 'web-mercator'
+            ? RADIUS * Math.asinh(Math.tan(latitude * RADIANS))
+            : latitude;
+    }
+}
+//# sourceMappingURL=geographic.js.map

package/dist/terrain/geographic.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geographic.js","sourceRoot":"","sources":["../../src/terrain/geographic.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EAAC,WAAW,EAAC,0BAAuB;AAC3C,OAAO,EACL,mBAAmB,IAAI,MAAM,EAI9B,mBAAgB;AAEjB,MAAM,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;AAC9B,MAAM,MAAM,GAAG,OAAO,CAAC;AACvB,MAAM,oBAAoB,GAAG,qBAAqB,CAAC;AACnD,MAAM,qBAAqB,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC;AAatE,mFAAmF;AACnF,MAAM,OAAO,qBAAqB;IAMhC,YAAY,KAAiC;QAC3C,MAAM,EAAC,MAAM,EAAC,GAAG,KAAK,CAAC;QACvB,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC,UAAU,IAAI,YAAY,CAAC;QACnD,IAAI,CAAC,eAAe,GAAG,KAAK,CAAC,eAAe,IAAI,SAAS,CAAC;QAC1D,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,KAAK,cAAc,CAAC,CAAC,CAAC,qBAAqB,CAAC,CAAC,CAAC,EAAE,CAAC;QAC9E,IACE,MAAM,CAAC,MAAM,KAAK,CAAC;YACnB,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAC9B,MAAM,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC;YACtB,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG;YAC3B,MAAM,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC;YACtB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK;YAClB,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK;YACjB,CAAC,IAAI,CAAC,UAAU,KAAK,YAAY,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,IAAI,MAAM,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,EAC7E,CAAC;YACD,MAAM,IAAI,KAAK,CACb,+FAA+F,CAChG,CAAC;QACJ,CAAC;QACD,IAAI,IAAI,CAAC,UAAU,KAAK,YAAY,IAAI,IAAI,CAAC,UAAU,KAAK,cAAc;YACxE,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;QAChD,MAAM,SAAS,GAAG,KAAK,CAAC,YAAY,IAAI,gBAAgB,CAAC;QACzD,IAAI,SAAS,KAAK,gBAAgB,IAAI,SAAS,KAAK,gBAAgB;YAClE,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;QACnD,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAkB,CAAC;QAC1D,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,KAAK,CAAC,YAAY,KAAK,MAAM,CAAC;QAC3C,MAAM,UAAU,GAAG,SAAS,KAAK,gBAAgB,CAAC;QAClD,IAAI,CAAC,KAAK,GAAG,IAAI,WAAW,CAAC;YAC3B,GAAG,KAAK;YACR,MAAM,EAAE,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC;YAC1C,UAAU,EAAE,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC/D,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC7F,CAAC,CAAC;IACL,CAAC;IAED,2FAA2F;IAC3F,MAAM,CACJ,SAAiB,EACjB,QAAgB,EAChB,OAA8B,EAC9B,MAAsB;QAEtB,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;QACtE,MAAM,MAAM,GACV,KAAK;YACL,SAAS,IAAI,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;YAC3B,SAAS,IAAI,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;YAC3B,QAAQ,IAAI,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;YAC1B,QAAQ,IAAI,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAC9B,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,EACvC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,EACtC,OAAO,EACP,MAAM,CACP,CAAC;QACF,IAAI,CAAC,MAAM,EAAE,CAAC;YACZ,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,cAAc,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,CAAC,MAAM,CAAC,aAAa,CAAC;YAC7F,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,IAAI,MAAM,CAAC,cAAc,KAAK,MAAM,CAAC,KAAK;YAAE,OAAO,MAAM,CAAC;QAC1D,MAAM,GAAG,GAAG,QAAQ,GAAG,OAAO,CAAC;QAC/B,MAAM,WAAW,GAAG,CAAC,GAAG,oBAAoB,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAClE,MAAM,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;QAC1C,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,oBAAoB,CAAC,CAAC,GAAG,WAAW,IAAI,GAAG,CAAC;QACrE,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,KAAK,cAAc,CAAC;QACpD,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,SAAS,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACxF,MAAM,EAAE,GAAG,CAAC,MAAM,CAAC,SAAS,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC;QACtF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;YACjD,MAAM,CAAC,cAAc,GAAG,MAAM,CAAC,aAAa,CAAC;YAC7C,MAAM,CAAC,SAAS,GAAG,MAAM,CAAC,SAAS,GAAG,MAAM,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,GAAG,MAAM,CAAC,OAAO,GAAG,CAAC,CAAC;YAC3F,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;QACtB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;QACtB,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;QAC7D,MAAM,CAAC,OAAO,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;QACtC,MAAM,CAAC,OAAO,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;QACtC,MAAM,CAAC,OAAO,GAAG,CAAC,GAAG,KAAK,GAAG,MAAM,CAAC;QACpC,OAAO,MAAM,CAAC;IAChB,CAAC;IAEO,QAAQ,CAAC,SAAiB;QAChC,OAAO,IAAI,CAAC,UAAU,KAAK,cAAc,CAAC,CAAC,CAAC,MAAM,GAAG,SAAS,GAAG,OAAO,CAAC,CAAC,CAAC,SAAS,CAAC;IACvF,CAAC;IAEO,QAAQ,CAAC,QAAgB;QAC/B,OAAO,IAAI,CAAC,UAAU,KAAK,cAAc;YACvC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,OAAO,CAAC,CAAC;YACnD,CAAC,CAAC,QAAQ,CAAC;IACf,CAAC;CACF"}

package/dist/terrain/geometry.d.ts

@@ -0,0 +1,14 @@
+import { Geometry } from '@math.gl/geometry';
+import type { HeightField } from "./height-field.js";
+export type HeightFieldGeometryOptions = {
+    id?: string;
+    /** World-space XYZ origin subtracted before storing positions. Defaults to [grid origin X, Y, 0]. */
+    origin?: readonly [number, number, number];
+    precision?: 'float32' | 'float64';
+};
+/** Z-up regular terrain mesh. Triangles and triangle sampling use the same cell diagonal. */
+export declare class HeightFieldGeometry extends Geometry {
+    readonly origin: readonly [number, number, number];
+    constructor(field: HeightField, options?: HeightFieldGeometryOptions);
+}
+//# sourceMappingURL=geometry.d.ts.map

package/dist/terrain/geometry.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geometry.d.ts","sourceRoot":"","sources":["../../src/terrain/geometry.ts"],"names":[],"mappings":"AAIA,OAAO,EAAC,QAAQ,EAAC,MAAM,mBAAmB,CAAC;AAC3C,OAAO,KAAK,EAAC,WAAW,EAAC,0BAAuB;AAEhD,MAAM,MAAM,0BAA0B,GAAG;IACvC,EAAE,CAAC,EAAE,MAAM,CAAC;IACZ,qGAAqG;IACrG,MAAM,CAAC,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3C,SAAS,CAAC,EAAE,SAAS,GAAG,SAAS,CAAC;CACnC,CAAC;AAEF,6FAA6F;AAC7F,qBAAa,mBAAoB,SAAQ,QAAQ;IAC/C,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;gBAEvC,KAAK,EAAE,WAAW,EAAE,OAAO,GAAE,0BAA+B;CA6GzE"}

package/dist/terrain/geometry.js

@@ -0,0 +1,116 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+import { Geometry } from '@math.gl/geometry';
+/** Z-up regular terrain mesh. Triangles and triangle sampling use the same cell diagonal. */
+export class HeightFieldGeometry extends Geometry {
+    constructor(field, options = {}) {
+        const origin = options.origin ?? [field.origin[0], field.origin[1], 0];
+        if (origin.length !== 3 || !origin.every(Number.isFinite))
+            throw new Error('Mesh origin must be a finite XYZ triple');
+        const precision = options.precision ?? 'float32';
+        if (precision !== 'float32' && precision !== 'float64')
+            throw new Error('Unknown terrain mesh precision');
+        const count = field.width * field.height;
+        const positions = precision === 'float64' ? new Float64Array(count * 3) : new Float32Array(count * 3);
+        const normals = new Float32Array(count * 3);
+        const accumulatedNormals = new Float64Array(count * 3);
+        const uv = new Float32Array(count * 2);
+        const valid = new Uint8Array(count);
+        const indices = count > 0xffff
+            ? new Uint32Array((field.width - 1) * (field.height - 1) * 6)
+            : new Uint16Array((field.width - 1) * (field.height - 1) * 6);
+        const shift = field.registration === 'cell' ? 0.5 : 0;
+        for (let row = 0; row < field.height; row++) {
+            for (let column = 0; column < field.width; column++) {
+                const i = row * field.width + column;
+                const height = field.getHeight(column, row);
+                valid[i] = height === undefined ? 0 : 1;
+                positions[3 * i] =
+                    field.origin[0] -
+                        origin[0] +
+                        (column + shift) * field.columnStep[0] +
+                        (row + shift) * field.rowStep[0];
+                positions[3 * i + 1] =
+                    field.origin[1] -
+                        origin[1] +
+                        (column + shift) * field.columnStep[1] +
+                        (row + shift) * field.rowStep[1];
+                positions[3 * i + 2] = height === undefined ? 0 : height - origin[2];
+                uv[2 * i] =
+                    (column + shift) / (field.registration === 'cell' ? field.width : field.width - 1);
+                uv[2 * i + 1] =
+                    (row + shift) / (field.registration === 'cell' ? field.height : field.height - 1);
+            }
+        }
+        if (!positions.every(Number.isFinite))
+            throw new Error('Terrain mesh positions overflow the selected precision');
+        const reverse = field.columnStep[0] * field.rowStep[1] - field.columnStep[1] * field.rowStep[0] < 0;
+        let indexCount = 0;
+        const addTriangle = (a, b, c) => {
+            if (!valid[a] || !valid[b] || !valid[c])
+                return;
+            if (reverse)
+                [b, c] = [c, b];
+            indices[indexCount++] = a;
+            indices[indexCount++] = b;
+            indices[indexCount++] = c;
+            const abx = positions[3 * b] - positions[3 * a];
+            const aby = positions[3 * b + 1] - positions[3 * a + 1];
+            const abz = positions[3 * b + 2] - positions[3 * a + 2];
+            const acx = positions[3 * c] - positions[3 * a];
+            const acy = positions[3 * c + 1] - positions[3 * a + 1];
+            const acz = positions[3 * c + 2] - positions[3 * a + 2];
+            const nx = aby * acz - abz * acy;
+            const ny = abz * acx - abx * acz;
+            const nz = abx * acy - aby * acx;
+            if (!Number.isFinite(nx) || !Number.isFinite(ny) || !Number.isFinite(nz) || nz <= 0) {
+                throw new Error('Terrain triangle degenerates or overflows at the selected precision; use a local origin');
+            }
+            for (const i of [a, b, c]) {
+                accumulatedNormals[3 * i] += nx;
+                accumulatedNormals[3 * i + 1] += ny;
+                accumulatedNormals[3 * i + 2] += nz;
+            }
+        };
+        for (let row = 0; row < field.height - 1; row++) {
+            for (let column = 0; column < field.width - 1; column++) {
+                const a = row * field.width + column;
+                const b = a + 1;
+                const c = a + field.width;
+                const d = c + 1;
+                if (field.diagonal === '00-11') {
+                    addTriangle(a, b, d);
+                    addTriangle(a, d, c);
+                }
+                else {
+                    addTriangle(a, b, c);
+                    addTriangle(b, d, c);
+                }
+            }
+        }
+        for (let i = 0; i < count; i++) {
+            const nx = accumulatedNormals[3 * i];
+            const ny = accumulatedNormals[3 * i + 1];
+            const nz = accumulatedNormals[3 * i + 2];
+            const length = Math.hypot(nx, ny, nz);
+            if (!Number.isFinite(length))
+                throw new Error('Terrain normal accumulation overflow');
+            normals[3 * i] = length ? nx / length : 0;
+            normals[3 * i + 1] = length ? ny / length : 0;
+            normals[3 * i + 2] = length ? nz / length : 1;
+        }
+        super({
+            id: options.id,
+            topology: 'triangle-list',
+            indices: indexCount === indices.length ? indices : indices.slice(0, indexCount),
+            attributes: {
+                POSITION: { size: 3, value: positions },
+                NORMAL: { size: 3, value: normals },
+                TEXCOORD_0: { size: 2, value: uv }
+            }
+        });
+        this.origin = Object.freeze([...origin]);
+    }
+}
+//# sourceMappingURL=geometry.js.map