@math.gl/geospatial 5.0.0-alpha.19 → 5.0.0-alpha.20

package/dist/discontinuous/elastic/geometry.d.ts

@@ -0,0 +1,22 @@
+import type { Position } from "../math.js";
+import type { MapSamplingOptions } from "../map-layout.js";
+export type Polygons = Position[][][];
+export declare const rectangle: (w: number, s: number, e: number, n: number) => Position[][];
+export declare function onSegment(p: readonly number[], a: readonly number[], b: readonly number[]): boolean;
+export declare function inRing(p: readonly number[], ring: readonly (readonly number[])[], includeBoundary?: boolean): boolean;
+export declare const contains: (polygons: Polygons, p: readonly number[]) => boolean;
+export declare const lerp: (a: readonly number[], b: readonly number[], t: number) => Position;
+/** Split a sampled source segment at every domain edge, including narrow lobes. */
+export declare function clipSegment(a: Position, b: Position, polygons: Polygons): Position[][];
+export declare function sampling(options: MapSamplingOptions): {
+    step: number;
+    tolerance: number;
+};
+export declare function sample(line: Position[], project: (p: Position) => Position | null, options: MapSamplingOptions, cellStep: number): {
+    source: Position[];
+    xy: Position[];
+};
+/** Latitude buckets accelerate repeated coverage queries without approximating
+ * the published boundary. Narrow cuts retain every original segment. */
+export declare function indexDomain(polygons: Polygons): (p: readonly number[]) => boolean;
+//# sourceMappingURL=geometry.d.ts.map

package/dist/discontinuous/elastic/geometry.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geometry.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/elastic/geometry.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAC,QAAQ,EAAC,mBAAgB;AACtC,OAAO,KAAK,EAAC,kBAAkB,EAAC,yBAAsB;AACtD,MAAM,MAAM,QAAQ,GAAG,QAAQ,EAAE,EAAE,EAAE,CAAC;AACtC,eAAO,MAAM,SAAS,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,QAAQ,EAAE,EAQhF,CAAC;AACF,wBAAgB,SAAS,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,WASzF;AACD,wBAAgB,MAAM,CACpB,CAAC,EAAE,SAAS,MAAM,EAAE,EACpB,IAAI,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,EACpC,eAAe,UAAO,WAcvB;AACD,eAAO,MAAM,QAAQ,GAAI,UAAU,QAAQ,EAAE,GAAG,SAAS,MAAM,EAAE,YAC8B,CAAC;AAChG,eAAO,MAAM,IAAI,GAAI,GAAG,SAAS,MAAM,EAAE,EAAE,GAAG,SAAS,MAAM,EAAE,EAAE,GAAG,MAAM,KAAG,QAG5E,CAAC;AACF,mFAAmF;AACnF,wBAAgB,WAAW,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,GAAG,QAAQ,EAAE,EAAE,CAuBtF;AACD,wBAAgB,QAAQ,CAAC,OAAO,EAAE,kBAAkB;;;EAYnD;AACD,wBAAgB,MAAM,CACpB,IAAI,EAAE,QAAQ,EAAE,EAChB,OAAO,EAAE,CAAC,CAAC,EAAE,QAAQ,KAAK,QAAQ,GAAG,IAAI,EACzC,OAAO,EAAE,kBAAkB,EAC3B,QAAQ,EAAE,MAAM;;;EAkCjB;AAED;wEACwE;AACxE,wBAAgB,WAAW,CAAC,QAAQ,EAAE,QAAQ,GAAG,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,OAAO,CA0CjF"}

package/dist/discontinuous/elastic/geometry.js

@@ -0,0 +1,131 @@
+export const rectangle = (w, s, e, n) => [
+    [
+        [w, s],
+        [e, s],
+        [e, n],
+        [w, n],
+        [w, s]
+    ]
+];
+export function onSegment(p, a, b) {
+    const dx = b[0] - a[0], dy = b[1] - a[1], l = Math.hypot(dx, dy);
+    return (l > 0 &&
+        Math.abs((p[0] - a[0]) * dy - (p[1] - a[1]) * dx) <= 1e-8 * l &&
+        (p[0] - a[0]) * (p[0] - b[0]) + (p[1] - a[1]) * (p[1] - b[1]) <= 1e-8);
+}
+export function inRing(p, ring, includeBoundary = true) {
+    let inside = false;
+    for (let i = 0, j = ring.length - 1; i < ring.length; j = i++) {
+        const a = ring[j], b = ring[i];
+        if (onSegment(p, a, b))
+            return includeBoundary;
+        if (a[1] > p[1] !== b[1] > p[1] &&
+            p[0] < ((b[0] - a[0]) * (p[1] - a[1])) / (b[1] - a[1]) + a[0])
+            inside = !inside;
+    }
+    return inside;
+}
+export const contains = (polygons, p) => polygons.some(rings => inRing(p, rings[0]) && !rings.slice(1).some(r => inRing(p, r, false)));
+export const lerp = (a, b, t) => [
+    a[0] + (b[0] - a[0]) * t,
+    a[1] + (b[1] - a[1]) * t
+];
+/** Split a sampled source segment at every domain edge, including narrow lobes. */
+export function clipSegment(a, b, polygons) {
+    const ts = [0, 1], dx = b[0] - a[0], dy = b[1] - a[1];
+    for (const rings of polygons)
+        for (const ring of rings)
+            for (let i = 1; i < ring.length; i++) {
+                const c = ring[i - 1], d = ring[i], ex = d[0] - c[0], ey = d[1] - c[1], det = dx * ey - dy * ex;
+                if (Math.abs(det) < 1e-14)
+                    continue;
+                const t = ((c[0] - a[0]) * ey - (c[1] - a[1]) * ex) / det;
+                const u = ((c[0] - a[0]) * dy - (c[1] - a[1]) * dx) / det;
+                if (t > 0 && t < 1 && u >= -1e-10 && u <= 1 + 1e-10)
+                    ts.push(t);
+            }
+    ts.sort((a, b) => a - b);
+    const result = [];
+    for (let i = 1; i < ts.length; i++)
+        if (ts[i] - ts[i - 1] > 1e-12 && contains(polygons, lerp(a, b, (ts[i] + ts[i - 1]) / 2)))
+            result.push([lerp(a, b, ts[i - 1]), lerp(a, b, ts[i])]);
+    return result;
+}
+export function sampling(options) {
+    const step = options.maxAngularStep ?? 2, tolerance = options.tolerance ?? 1e-4;
+    if (!Number.isFinite(step) ||
+        step <= 0 ||
+        step > 90 ||
+        !Number.isFinite(tolerance) ||
+        tolerance <= 0)
+        throw new Error('Invalid map sampling options');
+    return { step, tolerance };
+}
+export function sample(line, project, options, cellStep) {
+    const { step, tolerance } = sampling(options), source = [], xy = [];
+    const emit = (p, q) => {
+        if (source.length >= 100000)
+            throw new Error('Elastic sampling limit exceeded');
+        source.push(p);
+        xy.push(q);
+    };
+    const at = (p) => {
+        const q = project(p);
+        if (!q)
+            throw new Error('Elastic grid is undefined inside its coverage');
+        return q;
+    };
+    const visit = (a, b, qa, qb, depth) => {
+        const mid = lerp(a, b, 0.5), qm = at(mid), chord = lerp(qa, qb, 0.5);
+        if (Math.hypot(b[0] - a[0], b[1] - a[1]) > Math.min(step, cellStep / 2) ||
+            Math.hypot(qm[0] - chord[0], qm[1] - chord[1]) > tolerance) {
+            if (depth >= 24)
+                throw new Error('Elastic sampling limit exceeded');
+            visit(a, mid, qa, qm, depth + 1);
+            visit(mid, b, qm, qb, depth + 1);
+        }
+        else
+            emit(b, qb);
+    };
+    if (line.length) {
+        emit(line[0], at(line[0]));
+        for (let i = 1; i < line.length; i++)
+            visit(line[i - 1], line[i], at(line[i - 1]), at(line[i]), 0);
+    }
+    return { source, xy };
+}
+/** Latitude buckets accelerate repeated coverage queries without approximating
+ * the published boundary. Narrow cuts retain every original segment. */
+export function indexDomain(polygons) {
+    const compiled = polygons.map(rings => rings.map(ring => {
+        let minX = Infinity, maxX = -Infinity, minY = Infinity, maxY = -Infinity;
+        for (const [x, y] of ring) {
+            minX = Math.min(minX, x);
+            maxX = Math.max(maxX, x);
+            minY = Math.min(minY, y);
+            maxY = Math.max(maxY, y);
+        }
+        const bins = Array.from({ length: 64 }, () => []);
+        const index = (y) => Math.max(0, Math.min(63, Math.floor(((y + 90) / 180) * 64)));
+        for (let i = 1; i < ring.length; i++) {
+            const a = ring[i - 1], b = ring[i];
+            for (let j = index(Math.min(a[1], b[1]) - 1e-8); j <= index(Math.max(a[1], b[1]) + 1e-8); j++)
+                bins[j].push([a, b]);
+        }
+        return (p, inclusive) => {
+            if (p[0] < minX - 1e-8 || p[0] > maxX + 1e-8 || p[1] < minY - 1e-8 || p[1] > maxY + 1e-8)
+                return false;
+            let inside = false;
+            for (const [a, b] of bins[index(p[1])]) {
+                if (onSegment(p, a, b))
+                    return inclusive;
+                if (a[1] > p[1] !== b[1] > p[1] &&
+                    p[0] < ((b[0] - a[0]) * (p[1] - a[1])) / (b[1] - a[1]) + a[0])
+                    inside = !inside;
+            }
+            return inside;
+        };
+    }));
+    return p => compiled.some(rings => rings[0](p, true) && !rings.slice(1).some(r => r(p, false)));
+}
+//# sourceMappingURL=geometry.js.map

package/dist/discontinuous/elastic/geometry.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geometry.js","sourceRoot":"","sources":["../../../src/discontinuous/elastic/geometry.ts"],"names":[],"mappings":"AAMA,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAgB,EAAE,CAAC;IACrF;QACE,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;KACP;CACF,CAAC;AACF,MAAM,UAAU,SAAS,CAAC,CAAoB,EAAE,CAAoB,EAAE,CAAoB;IACxF,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EACpB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IACzB,OAAO,CACL,CAAC,GAAG,CAAC;QACL,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,IAAI,GAAG,CAAC;QAC7D,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CACtE,CAAC;AACJ,CAAC;AACD,MAAM,UAAU,MAAM,CACpB,CAAoB,EACpB,IAAoC,EACpC,eAAe,GAAG,IAAI;IAEtB,IAAI,MAAM,GAAG,KAAK,CAAC;IACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC;QAC9D,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EACf,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACd,IAAI,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YAAE,OAAO,eAAe,CAAC;QAC/C,IACE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAE7D,MAAM,GAAG,CAAC,MAAM,CAAC;IACrB,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AACD,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,QAAkB,EAAE,CAAoB,EAAE,EAAE,CACnE,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;AAChG,MAAM,CAAC,MAAM,IAAI,GAAG,CAAC,CAAoB,EAAE,CAAoB,EAAE,CAAS,EAAY,EAAE,CAAC;IACvF,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACxB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;CACzB,CAAC;AACF,mFAAmF;AACnF,MAAM,UAAU,WAAW,CAAC,CAAW,EAAE,CAAW,EAAE,QAAkB;IACtE,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EACf,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACnB,KAAK,MAAM,KAAK,IAAI,QAAQ;QAC1B,KAAK,MAAM,IAAI,IAAI,KAAK;YACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,EACnB,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EACX,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;gBAC1B,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK;oBAAE,SAAS;gBACpC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;gBAC1D,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC;gBAC1D,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,CAAC,GAAG,KAAK;oBAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAClE,CAAC;IACL,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACzB,MAAM,MAAM,GAAiB,EAAE,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;QAChC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,IAAI,QAAQ,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACtF,MAAM,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,OAAO,MAAM,CAAC;AAChB,CAAC;AACD,MAAM,UAAU,QAAQ,CAAC,OAA2B;IAClD,MAAM,IAAI,GAAG,OAAO,CAAC,cAAc,IAAI,CAAC,EACtC,SAAS,GAAG,OAAO,CAAC,SAAS,IAAI,IAAI,CAAC;IACxC,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QACtB,IAAI,IAAI,CAAC;QACT,IAAI,GAAG,EAAE;QACT,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC;QAC3B,SAAS,IAAI,CAAC;QAEd,MAAM,IAAI,KAAK,CAAC,8BAA8B,CAAC,CAAC;IAClD,OAAO,EAAC,IAAI,EAAE,SAAS,EAAC,CAAC;AAC3B,CAAC;AACD,MAAM,UAAU,MAAM,CACpB,IAAgB,EAChB,OAAyC,EACzC,OAA2B,EAC3B,QAAgB;IAEhB,MAAM,EAAC,IAAI,EAAE,SAAS,EAAC,GAAG,QAAQ,CAAC,OAAO,CAAC,EACzC,MAAM,GAAe,EAAE,EACvB,EAAE,GAAe,EAAE,CAAC;IACtB,MAAM,IAAI,GAAG,CAAC,CAAW,EAAE,CAAW,EAAE,EAAE;QACxC,IAAI,MAAM,CAAC,MAAM,IAAI,MAAM;YAAE,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;QAChF,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACb,CAAC,CAAC;IACF,MAAM,EAAE,GAAG,CAAC,CAAW,EAAE,EAAE;QACzB,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACrB,IAAI,CAAC,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;QACzE,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;IACF,MAAM,KAAK,GAAG,CAAC,CAAW,EAAE,CAAW,EAAE,EAAY,EAAE,EAAY,EAAE,KAAa,EAAE,EAAE;QACpF,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,EACzB,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EACZ,KAAK,GAAG,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAC5B,IACE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,QAAQ,GAAG,CAAC,CAAC;YACnE,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,SAAS,EAC1D,CAAC;YACD,IAAI,KAAK,IAAI,EAAE;gBAAE,MAAM,IAAI,KAAK,CAAC,iCAAiC,CAAC,CAAC;YACpE,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;YACjC,KAAK,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;QACnC,CAAC;;YAAM,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IACrB,CAAC,CAAC;IACF,IAAI,IAAI,CAAC,MAAM,EAAE,CAAC;QAChB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE;YAClC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACjE,CAAC;IACD,OAAO,EAAC,MAAM,EAAE,EAAE,EAAC,CAAC;AACtB,CAAC;AAED;wEACwE;AACxE,MAAM,UAAU,WAAW,CAAC,QAAkB;IAC5C,MAAM,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CACpC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;QACf,IAAI,IAAI,GAAG,QAAQ,EACjB,IAAI,GAAG,CAAC,QAAQ,EAChB,IAAI,GAAG,QAAQ,EACf,IAAI,GAAG,CAAC,QAAQ,CAAC;QACnB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,IAAI,EAAE,CAAC;YAC1B,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAC3B,CAAC;QACD,MAAM,IAAI,GAAmB,KAAK,CAAC,IAAI,CAAC,EAAC,MAAM,EAAE,EAAE,EAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QAChE,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,EACnB,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACd,KACE,IAAI,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAC1C,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EACvC,CAAC,EAAE;gBAEH,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACzB,CAAC;QACD,OAAO,CAAC,CAAoB,EAAE,SAAkB,EAAE,EAAE;YAClD,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI;gBACtF,OAAO,KAAK,CAAC;YACf,IAAI,MAAM,GAAG,KAAK,CAAC;YACnB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACvC,IAAI,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBAAE,OAAO,SAAS,CAAC;gBACzC,IACE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBAE7D,MAAM,GAAG,CAAC,MAAM,CAAC;YACrB,CAAC;YACD,OAAO,MAAM,CAAC;QAChB,CAAC,CAAC;IACJ,CAAC,CAAC,CACH,CAAC;IACF,OAAO,CAAC,CAAC,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;AAClG,CAAC"}

package/dist/discontinuous/elastic/grid.d.ts

@@ -0,0 +1,15 @@
+import type { ElasticSection } from "./data.js";
+import type { Position } from "../math.js";
+/** Bilinear forward interpolation and cell-local inversion. No interpolation
+ * across missing nodes or section cuts; inverse candidates are residual checked. */
+export declare class ElasticGrid {
+    readonly section: ElasticSection;
+    private scale;
+    private cells;
+    constructor(section: ElasticSection, scale: number);
+    project([lon, lat]: readonly number[]): Position | null;
+    /** Every interpolation node needed for this geographic rectangle is defined. */
+    covers([w, s, e, n]: readonly number[]): boolean;
+    inverse(target: readonly number[]): Position[];
+}
+//# sourceMappingURL=grid.d.ts.map

package/dist/discontinuous/elastic/grid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"grid.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/elastic/grid.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAC,cAAc,EAAC,kBAAe;AAC3C,OAAO,KAAK,EAAC,QAAQ,EAAC,mBAAgB;AAGtC;oFACoF;AACpF,qBAAa,WAAW;IAGpB,QAAQ,CAAC,OAAO,EAAE,cAAc;IAChC,OAAO,CAAC,KAAK;IAHf,OAAO,CAAC,KAAK,CAAc;gBAEhB,OAAO,EAAE,cAAc,EACxB,KAAK,EAAE,MAAM;IA8CvB,OAAO,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,SAAS,MAAM,EAAE,GAAG,QAAQ,GAAG,IAAI;IAyBvD,gFAAgF;IAChF,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,GAAG,OAAO;IAUhD,OAAO,CAAC,MAAM,EAAE,SAAS,MAAM,EAAE,GAAG,QAAQ,EAAE;CAsB/C"}

package/dist/discontinuous/elastic/grid.js

@@ -0,0 +1,177 @@
+/** Bilinear forward interpolation and cell-local inversion. No interpolation
+ * across missing nodes or section cuts; inverse candidates are residual checked. */
+export class ElasticGrid {
+    constructor(section, scale) {
+        this.section = section;
+        this.scale = scale;
+        this.cells = [];
+        if (!Number.isFinite(scale) || scale <= 0)
+            throw new Error('Elastic output scale is not representable');
+        const { rows, columns, values } = section;
+        if (!Number.isInteger(rows) ||
+            !Number.isInteger(columns) ||
+            rows < 2 ||
+            columns < 2 ||
+            rows * columns > 1000000 ||
+            values.length !== rows * columns * 2 ||
+            values.some(v => v !== null && !Number.isFinite(v)))
+            throw new Error('Invalid Elastic grid');
+        for (let i = 0; i < values.length; i += 2)
+            if ((values[i] === null) !== (values[i + 1] === null))
+                throw new Error('Incomplete Elastic node');
+        for (let y = 0; y < rows - 1; y++)
+            for (let x = 0; x < columns - 1; x++) {
+                const indices = [
+                    y * columns + x,
+                    y * columns + x + 1,
+                    (y + 1) * columns + x,
+                    (y + 1) * columns + x + 1
+                ];
+                if (indices.some(i => values[2 * i] === null))
+                    continue;
+                const p = indices.map(i => [values[2 * i] * scale, values[2 * i + 1] * scale]);
+                if (!p.flat().every(Number.isFinite))
+                    throw new Error('Elastic output coordinates are not representable');
+                this.cells.push({
+                    x,
+                    y,
+                    p,
+                    bounds: [
+                        Math.min(...p.map(v => v[0])),
+                        Math.min(...p.map(v => v[1])),
+                        Math.max(...p.map(v => v[0])),
+                        Math.max(...p.map(v => v[1]))
+                    ]
+                });
+            }
+        if (!this.cells.length)
+            throw new Error('Elastic section has no defined cells');
+    }
+    project([lon, lat]) {
+        const { rows, columns, values } = this.section;
+        const x = ((lon + 180) / 360) * (columns - 1), y = ((lat + 90) / 180) * (rows - 1);
+        if (x < -1e-9 || x > columns - 1 + 1e-9 || y < -1e-9 || y > rows - 1 + 1e-9)
+            return null;
+        const i = Math.max(0, Math.min(columns - 2, Math.floor(x))), j = Math.max(0, Math.min(rows - 2, Math.floor(y)));
+        const u = x - i, v = y - j;
+        const indices = [
+            j * columns + i,
+            j * columns + i + 1,
+            (j + 1) * columns + i,
+            (j + 1) * columns + i + 1
+        ];
+        const weights = [(1 - u) * (1 - v), u * (1 - v), (1 - u) * v, u * v];
+        const result = [0, 0];
+        for (let k = 0; k < 4; k++) {
+            if (Math.abs(weights[k]) < 1e-14)
+                continue;
+            if (values[indices[k] * 2] === null)
+                return null;
+            for (let axis = 0; axis < 2; axis++)
+                result[axis] += weights[k] * values[indices[k] * 2 + axis] * this.scale;
+        }
+        return result;
+    }
+    /** Every interpolation node needed for this geographic rectangle is defined. */
+    covers([w, s, e, n]) {
+        const { rows, columns, values } = this.section;
+        const left = Math.max(0, Math.floor(((w + 180) / 360) * (columns - 1))), right = Math.min(columns - 1, Math.ceil(((e + 180) / 360) * (columns - 1)));
+        const bottom = Math.max(0, Math.floor(((s + 90) / 180) * (rows - 1))), top = Math.min(rows - 1, Math.ceil(((n + 90) / 180) * (rows - 1)));
+        for (let y = bottom; y <= top; y++)
+            for (let x = left; x <= right; x++)
+                if (values[2 * (y * columns + x)] === null)
+                    return false;
+        return true;
+    }
+    inverse(target) {
+        const results = [];
+        const epsilon = Math.max(this.scale * 1e-5, Math.hypot(...target) * 1e-10);
+        for (const { x, y, p, bounds: b } of this.cells) {
+            if (target[0] < b[0] - epsilon ||
+                target[0] > b[2] + epsilon ||
+                target[1] < b[1] - epsilon ||
+                target[1] > b[3] + epsilon)
+                continue;
+            for (const [u, v] of invertCell(p, target)) {
+                const point = [
+                    -180 + (360 * (x + u)) / (this.section.columns - 1),
+                    -90 + (180 * (y + v)) / (this.section.rows - 1)
+                ];
+                if (!results.some(r => Math.hypot(r[0] - point[0], r[1] - point[1]) < 1e-7))
+                    results.push(point);
+            }
+        }
+        return results;
+    }
+}
+const cross2 = (a, b) => a[0] * b[1] - a[1] * b[0];
+/** Eliminate either parameter to solve the bilinear patch's quadratic. Solving
+ * both directions also handles affine cells and a vanishing elimination axis.
+ * Normalize locally so output radius does not change inverse conditioning. */
+function invertCell(p, target) {
+    const scale = Math.max(...p.slice(1).flatMap(q => q.map((v, i) => Math.abs(v - p[0][i]))));
+    if (!(scale > 0))
+        return target.every((v, i) => v === p[0][i])
+            ? [
+                [0, 0],
+                [1, 0],
+                [0, 1],
+                [1, 1],
+                [0.5, 0.5]
+            ]
+            : [];
+    const b = p[1].map((v, i) => (v - p[0][i]) / scale), c = p[2].map((v, i) => (v - p[0][i]) / scale);
+    const d = p[3].map((v, i) => (v - p[0][i]) / scale - b[i] - c[i]);
+    const h = target.map((v, i) => (v - p[0][i]) / scale);
+    const roots = (a, b, c) => {
+        if (Math.abs(a) < 1e-14)
+            return Math.abs(b) < 1e-14 ? [] : [-c / b];
+        const discriminant = b * b - 4 * a * c;
+        if (discriminant < -1e-14)
+            return [];
+        const q = -0.5 * (b + (b < 0 ? -1 : 1) * Math.sqrt(Math.max(0, discriminant)));
+        return q === 0 ? [-b / (2 * a)] : [q / a, c / q];
+    };
+    const output = [];
+    const accept = (xy) => {
+        if (!xy.every(v => Number.isFinite(v) && v >= -1e-9 && v <= 1 + 1e-9))
+            return;
+        const [x, y] = xy.map(v => Math.max(0, Math.min(1, v)));
+        if (Math.hypot(...h.map((v, i) => b[i] * x + c[i] * y + d[i] * x * y - v)) > 1e-9)
+            return;
+        if (!output.some(p => Math.hypot(p[0] - x, p[1] - y) < 1e-8))
+            output.push([x, y]);
+    };
+    // Elimination loses solutions when a whole edge collapses (notably poles).
+    // Test edge restrictions directly, retaining residual-checked representatives
+    // of a collapsed edge rather than treating its zero denominator as no solution.
+    const corners = [
+        [0, 0],
+        [1, 0],
+        [1, 1],
+        [0, 1]
+    ];
+    const value = ([x, y]) => b.map((v, i) => v * x + c[i] * y + d[i] * x * y);
+    for (let i = 0; i < 4; i++) {
+        const start = corners[i], end = corners[(i + 1) % 4], a = value(start), z = value(end);
+        const delta = z.map((v, j) => v - a[j]);
+        const axis = Math.abs(delta[0]) > Math.abs(delta[1]) ? 0 : 1;
+        const parameters = Math.abs(delta[axis]) < 1e-14 ? [0, 0.5, 1] : [(h[axis] - a[axis]) / delta[axis]];
+        for (const t of parameters)
+            accept([start[0] + t * (end[0] - start[0]), start[1] + t * (end[1] - start[1])]);
+    }
+    for (const swap of [false, true]) {
+        const first = swap ? c : b, second = swap ? b : c;
+        for (const u of roots(-cross2(first, d), cross2(h, d) - cross2(first, second), cross2(h, second))) {
+            const denominator = second.map((v, i) => v + d[i] * u), numerator = h.map((v, i) => v - first[i] * u);
+            const axis = Math.abs(denominator[0]) > Math.abs(denominator[1]) ? 0 : 1;
+            if (Math.abs(denominator[axis]) < 1e-14)
+                continue;
+            const v = numerator[axis] / denominator[axis];
+            const xy = swap ? [v, u] : [u, v];
+            accept(xy);
+        }
+    }
+    return output;
+}
+//# sourceMappingURL=grid.js.map

package/dist/discontinuous/elastic/grid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"grid.js","sourceRoot":"","sources":["../../../src/discontinuous/elastic/grid.ts"],"names":[],"mappings":"AAOA;oFACoF;AACpF,MAAM,OAAO,WAAW;IAEtB,YACW,OAAuB,EACxB,KAAa;QADZ,YAAO,GAAP,OAAO,CAAgB;QACxB,UAAK,GAAL,KAAK,CAAQ;QAHf,UAAK,GAAW,EAAE,CAAC;QAKzB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;YACvC,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;QAC/D,MAAM,EAAC,IAAI,EAAE,OAAO,EAAE,MAAM,EAAC,GAAG,OAAO,CAAC;QACxC,IACE,CAAC,MAAM,CAAC,SAAS,CAAC,IAAI,CAAC;YACvB,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC;YAC1B,IAAI,GAAG,CAAC;YACR,OAAO,GAAG,CAAC;YACX,IAAI,GAAG,OAAO,GAAG,OAAS;YAC1B,MAAM,CAAC,MAAM,KAAK,IAAI,GAAG,OAAO,GAAG,CAAC;YACpC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAEnD,MAAM,IAAI,KAAK,CAAC,sBAAsB,CAAC,CAAC;QAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC;YACvC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,IAAI,CAAC;gBACnD,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;QAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE;YAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACrC,MAAM,OAAO,GAAG;oBACd,CAAC,GAAG,OAAO,GAAG,CAAC;oBACf,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;oBACnB,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC;oBACrB,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;iBAC1B,CAAC;gBACF,IAAI,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,IAAI,CAAC;oBAAE,SAAS;gBACxD,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CACnB,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,KAAK,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,KAAK,CAAa,CACtE,CAAC;gBACF,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;oBAClC,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;gBACtE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC;oBACd,CAAC;oBACD,CAAC;oBACD,CAAC;oBACD,MAAM,EAAE;wBACN,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBAC7B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBAC7B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBAC7B,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;qBAC9B;iBACF,CAAC,CAAC;YACL,CAAC;QACH,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM;YAAE,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;IAClF,CAAC;IACD,OAAO,CAAC,CAAC,GAAG,EAAE,GAAG,CAAoB;QACnC,MAAM,EAAC,IAAI,EAAE,OAAO,EAAE,MAAM,EAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,EAC3C,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QACtC,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI;YAAE,OAAO,IAAI,CAAC;QACzF,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EACzD,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACrD,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EACb,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,OAAO,GAAG;YACd,CAAC,GAAG,OAAO,GAAG,CAAC;YACf,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;YACnB,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC;YACrB,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;SAC1B,CAAC;QACF,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QACrE,MAAM,MAAM,GAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK;gBAAE,SAAS;YAC3C,IAAI,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,IAAI;gBAAE,OAAO,IAAI,CAAC;YACjD,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE;gBACjC,MAAM,CAAC,IAAI,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QAC7E,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,gFAAgF;IAChF,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAoB;QACpC,MAAM,EAAC,IAAI,EAAE,OAAO,EAAE,MAAM,EAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QAC7C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,EACrE,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9E,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,EACnE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,KAAK,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAChC,KAAK,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE;gBAAE,IAAI,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI;oBAAE,OAAO,KAAK,CAAC;QAC/F,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,CAAC,MAAyB;QAC/B,MAAM,OAAO,GAAe,EAAE,CAAC;QAC/B,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC;QAC3E,KAAK,MAAM,EAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAC,IAAI,IAAI,CAAC,KAAK,EAAE,CAAC;YAC9C,IACE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;gBAC1B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;gBAC1B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;gBAC1B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO;gBAE1B,SAAS;YACX,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,MAAM,CAAC,EAAE,CAAC;gBAC3C,MAAM,KAAK,GAAa;oBACtB,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,GAAG,CAAC,CAAC;oBACnD,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,GAAG,CAAC,CAAC;iBAChD,CAAC;gBACF,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;oBACzE,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACxB,CAAC;QACH,CAAC;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;CACF;AAED,MAAM,MAAM,GAAG,CAAC,CAAoB,EAAE,CAAoB,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AACzF;;8EAE8E;AAC9E,SAAS,UAAU,CAAC,CAAa,EAAE,MAAyB;IAC1D,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3F,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QACd,OAAO,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,CAAC,CAAC;gBACE,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,GAAG,EAAE,GAAG,CAAC;aACX;YACH,CAAC,CAAC,EAAE,CAAC;IACT,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EACjD,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IAChD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACtD,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE;QAChD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK;YAAE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACpE,MAAM,YAAY,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvC,IAAI,YAAY,GAAG,CAAC,KAAK;YAAE,OAAO,EAAE,CAAC;QACrC,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/E,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IACnD,CAAC,CAAC;IACF,MAAM,MAAM,GAAe,EAAE,CAAC;IAC9B,MAAM,MAAM,GAAG,CAAC,EAAY,EAAE,EAAE;QAC9B,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;YAAE,OAAO;QAC9E,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACxD,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI;YAAE,OAAO;QAC1F,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;YAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACpF,CAAC,CAAC;IACF,2EAA2E;IAC3E,8EAA8E;IAC9E,gFAAgF;IAChF,MAAM,OAAO,GAAe;QAC1B,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;KACP,CAAC;IACF,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAW,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACrF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,EACtB,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAC1B,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,EAChB,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7D,MAAM,UAAU,GACd,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;QACpF,KAAK,MAAM,CAAC,IAAI,UAAU;YACxB,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrF,CAAC;IACD,KAAK,MAAM,IAAI,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,EAAE,CAAC;QACjC,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EACxB,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxB,KAAK,MAAM,CAAC,IAAI,KAAK,CACnB,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,CAAC,EACjB,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,EACpC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAClB,EAAE,CAAC;YACF,MAAM,WAAW,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EACpD,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAChD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzE,IAAI,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK;gBAAE,SAAS;YAClD,MAAM,CAAC,GAAG,SAAS,CAAC,IAAI,CAAC,GAAG,WAAW,CAAC,IAAI,CAAC,CAAC;YAC9C,MAAM,EAAE,GAAa,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5C,MAAM,CAAC,EAAE,CAAC,CAAC;QACb,CAAC;IACH,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/dist/discontinuous/elastic/map.d.ts

@@ -0,0 +1,8 @@
+import type { MapSurface } from "../map-layout.js";
+import type { ElasticMapData } from "./data.js";
+/** Compile a sampled Elastic dataset. Output units are radius-scaled, matching
+ * other MapEngine maps. Dataset tables remain explicit, optional imports. */
+export declare function createElasticMap(data: ElasticMapData, options?: {
+    radius?: number;
+}): MapSurface;
+//# sourceMappingURL=map.d.ts.map

package/dist/discontinuous/elastic/map.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"map.d.ts","sourceRoot":"","sources":["../../../src/discontinuous/elastic/map.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAQV,UAAU,EACX,yBAAsB;AACvB,OAAO,KAAK,EAAC,cAAc,EAAC,kBAAe;AAkB3C;6EAC6E;AAC7E,wBAAgB,gBAAgB,CAC9B,IAAI,EAAE,cAAc,EACpB,OAAO,GAAE;IAAC,MAAM,CAAC,EAAE,MAAM,CAAA;CAAM,GAC9B,UAAU,CAEZ"}