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

package/dist/tile-matrix.d.ts

@@ -0,0 +1,107 @@
+/** Bounds in the matrix coordinate system: minimum X/Y followed by maximum X/Y. */
+export type TileMatrixBounds = readonly [number, number, number, number];
+/** Coalescence of contiguous base-grid columns over an inclusive band of rows. */
+export type VariableMatrixWidth = {
+    /** Number of base columns per tile; an integer greater than one dividing matrix width. */
+    readonly coalesce: number;
+    /** First row in this band, in the matrix's native row direction. */
+    readonly minTileRow: number;
+    /** Last row in this band, inclusive. */
+    readonly maxTileRow: number;
+};
+/** A rectangular tile grid at one resolution, independent of service protocols. */
+export type TileMatrix = {
+    /** Opaque identifier; it is not a numeric zoom or array index. */
+    readonly id: string;
+    /** X coordinate units per pixel, positive and finite; also used for Y unless overridden. */
+    readonly resolution: number;
+    /** Y coordinate units per pixel, positive and finite; defaults to resolution. */
+    readonly resolutionY?: number;
+    /** Origin in X/Y order in the matrix coordinate system. */
+    readonly origin: readonly [number, number];
+    /** Origin corner; columns increase in X, rows decrease/increase in Y for top/bottom left. */
+    readonly cornerOfOrigin: 'topLeft' | 'bottomLeft';
+    /** Positive integer tile width and height in pixels. */
+    readonly tileSize: readonly [number, number];
+    /** Positive integer base-grid width and height in tiles, before coalescence. */
+    readonly matrixSize: readonly [number, number];
+    /** Optional disjoint row bands; omitted rows have coalescence one. Order is immaterial. */
+    readonly variableMatrixWidths?: readonly VariableMatrixWidth[];
+};
+/** Matrices sharing a coordinate system; no ordering or quadtree relationship is assumed. */
+export type TileMatrixSet = {
+    /** Opaque coordinate-system identifier. Utilities do not transform coordinates. */
+    readonly crs: string;
+    /** Matrices with unique identifiers. */
+    readonly matrices: readonly TileMatrix[];
+};
+/** Inclusive column and row limits for a rectangular subset of a matrix. */
+export type TileMatrixLimits = {
+    /** First available column. */
+    readonly minTileColumn: number;
+    /** Last available column. */
+    readonly maxTileColumn: number;
+    /** First available row. */
+    readonly minTileRow: number;
+    /** Last available row. */
+    readonly maxTileRow: number;
+};
+/** Inclusive base-grid coverage identified by native matrix ID; omitted levels have no coverage. */
+export type TileMatrixSetLimits = readonly (TileMatrixLimits & {
+    /** Native identifier of the matrix addressed by this coverage rectangle. */
+    readonly tileMatrix: string;
+})[];
+/** An exact row-band range of canonical columns, enumerated using columnStep. */
+export type TileMatrixRowRange = TileMatrixLimits & {
+    /** Increment between canonical columns; rows increment by one. */
+    readonly columnStep: number;
+};
+/** Rejects malformed or numerically unrepresentable matrix geometry. */
+export declare function validateTileMatrix(matrix: TileMatrix): void;
+/**
+ * Validates complete geometry, a nonempty opaque CRS identifier, unique matrix IDs and optional
+ * identified coverage. Empty sets and empty coverage lists are valid. Does not mutate inputs,
+ * interpret the CRS, or require ordered levels or quadtree relationships.
+ */
+export declare function validateTileMatrixSet(matrixSet: TileMatrixSet, limits?: TileMatrixSetLimits): void;
+/** Validates complete matrix geometry and inclusive base-grid coverage without mutating either. */
+export declare function validateTileMatrixLimits(matrix: TileMatrix, limits: TileMatrixLimits): void;
+/** Tests integer indices against the matrix and optional inclusive coverage limits. Does not wrap. */
+export declare function isTileIndexInRange(matrix: TileMatrix, column: number, row: number, limits?: TileMatrixLimits): boolean;
+/**
+ * Resolves integer base-grid indices to a canonical column/row, or null for noninteger,
+ * nonfinite, outside-matrix or unavailable indices. A coalesced tile is available when any of
+ * its aliases intersects the inclusive base-grid limits; the canonical column may precede
+ * minTileColumn. Validates geometry and limits even for an invalid query. Does not wrap or mutate.
+ */
+export declare function resolveTileIndex(matrix: TileMatrix, column: number, row: number, limits?: TileMatrixLimits): [number, number] | null;
+/** Returns the entire tile footprint, including column aliases. Throws for out-of-range indices. */
+export declare function getTileBounds(matrix: TileMatrix, column: number, row: number): TileMatrixBounds;
+/**
+ * Finds the containing tile without clamping or wrapping. Tile intervals are half-open in column
+ * and row direction: the origin edge belongs to the matrix, the opposite outer edge does not.
+ * Returns a canonical column for coalesced rows, or null for nonfinite/outside coordinates.
+ */
+export declare function getTileIndex(matrix: TileMatrix, coordinate: readonly [number, number]): [number, number] | null;
+/**
+ * Returns a rectangular range for a regular matrix. Variable-width matrices require getTileRanges
+ * so callers cannot accidentally enumerate aliases as separate tiles.
+ */
+export declare function getTileRange(matrix: TileMatrix, bounds: TileMatrixBounds, limits?: TileMatrixLimits): TileMatrixLimits | null;
+/**
+ * Returns exact canonical-column ranges, ordered by native row, split at coalescence changes.
+ * Columns advance by columnStep, rows by one. Limits use base-grid indices: a tile is included
+ * when any of its aliases falls within those limits. Canonical columns may precede minTileColumn.
+ * Work and output scale with row bands rather than the number of tiles or rows.
+ */
+export declare function getTileRanges(matrix: TileMatrix, bounds: TileMatrixBounds, limits?: TileMatrixLimits): TileMatrixRowRange[];
+/**
+ * Selects the coarsest matrix with resolution at or below the positive finite target, or the
+ * finest available matrix when all levels are coarser. Target units must match the matrix set's
+ * coordinate units per pixel. Returns null for an empty set; ties retain the first matrix.
+ * Uses nominal uncoalesced resolution; selection is independent of row.
+ * Uses the larger axis resolution for anisotropic pixels, so both axes meet the target.
+ * Validates the complete set, including its CRS identifier. Does not sort or mutate the set.
+ */
+export declare function selectTileMatrix(matrixSet: TileMatrixSet, targetResolution: number): TileMatrix | null;
+//# sourceMappingURL=tile-matrix.d.ts.map

package/dist/tile-matrix.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"tile-matrix.d.ts","sourceRoot":"","sources":["../src/tile-matrix.ts"],"names":[],"mappings":"AAIA,mFAAmF;AACnF,MAAM,MAAM,gBAAgB,GAAG,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AAEzE,kFAAkF;AAClF,MAAM,MAAM,mBAAmB,GAAG;IAChC,0FAA0F;IAC1F,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,oEAAoE;IACpE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;IAC5B,wCAAwC;IACxC,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B,CAAC;AAEF,mFAAmF;AACnF,MAAM,MAAM,UAAU,GAAG;IACvB,kEAAkE;IAClE,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,4FAA4F;IAC5F,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;IAC5B,iFAAiF;IACjF,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;IAC9B,2DAA2D;IAC3D,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC3C,6FAA6F;IAC7F,QAAQ,CAAC,cAAc,EAAE,SAAS,GAAG,YAAY,CAAC;IAClD,wDAAwD;IACxD,QAAQ,CAAC,QAAQ,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7C,gFAAgF;IAChF,QAAQ,CAAC,UAAU,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,2FAA2F;IAC3F,QAAQ,CAAC,oBAAoB,CAAC,EAAE,SAAS,mBAAmB,EAAE,CAAC;CAChE,CAAC;AAEF,6FAA6F;AAC7F,MAAM,MAAM,aAAa,GAAG;IAC1B,mFAAmF;IACnF,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,wCAAwC;IACxC,QAAQ,CAAC,QAAQ,EAAE,SAAS,UAAU,EAAE,CAAC;CAC1C,CAAC;AAEF,4EAA4E;AAC5E,MAAM,MAAM,gBAAgB,GAAG;IAC7B,8BAA8B;IAC9B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,6BAA6B;IAC7B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,2BAA2B;IAC3B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;IAC5B,0BAA0B;IAC1B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B,CAAC;AAEF,oGAAoG;AACpG,MAAM,MAAM,mBAAmB,GAAG,SAAS,CAAC,gBAAgB,GAAG;IAC7D,4EAA4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B,CAAC,EAAE,CAAC;AAEL,iFAAiF;AACjF,MAAM,MAAM,kBAAkB,GAAG,gBAAgB,GAAG;IAClD,kEAAkE;IAClE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B,CAAC;AAEF,wEAAwE;AACxE,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,UAAU,GAAG,IAAI,CAE3D;AAED;;;;GAIG;AACH,wBAAgB,qBAAqB,CACnC,SAAS,EAAE,aAAa,EACxB,MAAM,CAAC,EAAE,mBAAmB,GAC3B,IAAI,CAgCN;AAED,mGAAmG;AACnG,wBAAgB,wBAAwB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,gBAAgB,GAAG,IAAI,CAG3F;AAgED,sGAAsG;AACtG,wBAAgB,kBAAkB,CAChC,MAAM,EAAE,UAAU,EAClB,MAAM,EAAE,MAAM,EACd,GAAG,EAAE,MAAM,EACX,MAAM,CAAC,EAAE,gBAAgB,GACxB,OAAO,CAWT;AAcD;;;;;GAKG;AACH,wBAAgB,gBAAgB,CAC9B,MAAM,EAAE,UAAU,EAClB,MAAM,EAAE,MAAM,EACd,GAAG,EAAE,MAAM,EACX,MAAM,CAAC,EAAE,gBAAgB,GACxB,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,IAAI,CAgBzB;AAED,oGAAoG;AACpG,wBAAgB,aAAa,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,gBAAgB,CAmB/F;AAED;;;;GAIG;AACH,wBAAgB,YAAY,CAC1B,MAAM,EAAE,UAAU,EAClB,UAAU,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,GACpC,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,IAAI,CAqBzB;AA4HD;;;GAGG;AACH,wBAAgB,YAAY,CAC1B,MAAM,EAAE,UAAU,EAClB,MAAM,EAAE,gBAAgB,EACxB,MAAM,CAAC,EAAE,gBAAgB,GACxB,gBAAgB,GAAG,IAAI,CAOzB;AAUD;;;;;GAKG;AACH,wBAAgB,aAAa,CAC3B,MAAM,EAAE,UAAU,EAClB,MAAM,EAAE,gBAAgB,EACxB,MAAM,CAAC,EAAE,gBAAgB,GACxB,kBAAkB,EAAE,CAgEtB;AAsCD;;;;;;;GAOG;AACH,wBAAgB,gBAAgB,CAC9B,SAAS,EAAE,aAAa,EACxB,gBAAgB,EAAE,MAAM,GACvB,UAAU,GAAG,IAAI,CAqBnB"}

package/dist/tile-matrix.js

@@ -0,0 +1,383 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+/** Rejects malformed or numerically unrepresentable matrix geometry. */
+export function validateTileMatrix(matrix) {
+    validateTileMatrixGeometry(matrix);
+}
+/**
+ * Validates complete geometry, a nonempty opaque CRS identifier, unique matrix IDs and optional
+ * identified coverage. Empty sets and empty coverage lists are valid. Does not mutate inputs,
+ * interpret the CRS, or require ordered levels or quadtree relationships.
+ */
+export function validateTileMatrixSet(matrixSet, limits) {
+    if (!matrixSet ||
+        typeof matrixSet.crs !== 'string' ||
+        !matrixSet.crs.trim() ||
+        !Array.isArray(matrixSet.matrices)) {
+        throw new RangeError('Invalid tile matrix set');
+    }
+    const matrices = new Map();
+    for (const matrix of matrixSet.matrices) {
+        validateTileMatrix(matrix);
+        if (matrices.has(matrix.id))
+            throw new RangeError('Tile matrix identifiers must be unique');
+        matrices.set(matrix.id, matrix);
+    }
+    if (limits === undefined)
+        return;
+    if (!Array.isArray(limits))
+        throw new RangeError('Tile matrix set limits require an array');
+    const identifiers = new Set();
+    for (const limit of limits) {
+        if (!limit ||
+            typeof limit.tileMatrix !== 'string' ||
+            !limit.tileMatrix ||
+            identifiers.has(limit.tileMatrix)) {
+            throw new RangeError('Tile matrix set limits require unique matrix identifiers');
+        }
+        const matrix = matrices.get(limit.tileMatrix);
+        if (!matrix)
+            throw new RangeError('Tile matrix limits reference an unknown matrix');
+        validateTileMatrixCoverage(matrix, limit);
+        identifiers.add(limit.tileMatrix);
+    }
+}
+/** Validates complete matrix geometry and inclusive base-grid coverage without mutating either. */
+export function validateTileMatrixLimits(matrix, limits) {
+    validateTileMatrix(matrix);
+    validateTileMatrixCoverage(matrix, limits);
+}
+/** Validates geometry and returns a sorted copy of the declared coalescence bands. */
+function validateTileMatrixGeometry(matrix) {
+    if (!matrix ||
+        typeof matrix.id !== 'string' ||
+        !Array.isArray(matrix.origin) ||
+        !Array.isArray(matrix.tileSize) ||
+        !Array.isArray(matrix.matrixSize) ||
+        (matrix.variableMatrixWidths !== undefined && !Array.isArray(matrix.variableMatrixWidths)) ||
+        !matrix.id ||
+        matrix.origin.length !== 2 ||
+        matrix.tileSize.length !== 2 ||
+        matrix.matrixSize.length !== 2 ||
+        !Number.isFinite(matrix.resolution) ||
+        matrix.resolution <= 0 ||
+        (matrix.resolutionY !== undefined &&
+            (!Number.isFinite(matrix.resolutionY) || matrix.resolutionY <= 0)) ||
+        !matrix.origin.every(Number.isFinite) ||
+        ![...matrix.tileSize, ...matrix.matrixSize].every(value => Number.isSafeInteger(value) && value > 0) ||
+        !['topLeft', 'bottomLeft'].includes(matrix.cornerOfOrigin)) {
+        throw new RangeError('Invalid tile matrix geometry');
+    }
+    const bands = [...(matrix.variableMatrixWidths ?? [])];
+    if (bands.some(band => !band))
+        throw new RangeError('Invalid variable matrix width');
+    bands.sort((first, second) => first.minTileRow - second.minTileRow);
+    let previousRow = -1;
+    for (const band of bands) {
+        if (!Number.isSafeInteger(band.coalesce) ||
+            band.coalesce <= 1 ||
+            matrix.matrixSize[0] % band.coalesce !== 0 ||
+            !Number.isSafeInteger(band.minTileRow) ||
+            !Number.isSafeInteger(band.maxTileRow) ||
+            band.minTileRow <= previousRow ||
+            band.maxTileRow < band.minTileRow ||
+            band.maxTileRow >= matrix.matrixSize[1])
+            throw new RangeError('Invalid variable matrix width');
+        previousRow = band.maxTileRow;
+    }
+    for (let axis = 0; axis < 2; axis++) {
+        const resolution = axis === 1 ? (matrix.resolutionY ?? matrix.resolution) : matrix.resolution;
+        const tileSpan = resolution * matrix.tileSize[axis];
+        const span = tileSpan * matrix.matrixSize[axis];
+        const direction = axis === 1 && matrix.cornerOfOrigin === 'topLeft' ? -1 : 1;
+        const end = matrix.origin[axis] + direction * span;
+        if (!Number.isFinite(span) ||
+            span <= 0 ||
+            !Number.isFinite(end) ||
+            matrix.origin[axis] + direction * tileSpan === matrix.origin[axis] ||
+            end - direction * tileSpan === end) {
+            throw new RangeError('Tile matrix extent must be finite');
+        }
+    }
+    return bands;
+}
+/** Tests integer indices against the matrix and optional inclusive coverage limits. Does not wrap. */
+export function isTileIndexInRange(matrix, column, row, limits) {
+    validateTileMatrix(matrix);
+    if (limits)
+        validateTileMatrixCoverage(matrix, limits);
+    return (isTileIndexWithinMatrix(matrix, column, row) &&
+        (!limits ||
+            (column >= limits.minTileColumn &&
+                column <= limits.maxTileColumn &&
+                row >= limits.minTileRow &&
+                row <= limits.maxTileRow)));
+}
+/** Tests base-grid indices against already validated dimensions. */
+function isTileIndexWithinMatrix(matrix, column, row) {
+    return (Number.isSafeInteger(column) &&
+        Number.isSafeInteger(row) &&
+        column >= 0 &&
+        row >= 0 &&
+        column < matrix.matrixSize[0] &&
+        row < matrix.matrixSize[1]);
+}
+/**
+ * Resolves integer base-grid indices to a canonical column/row, or null for noninteger,
+ * nonfinite, outside-matrix or unavailable indices. A coalesced tile is available when any of
+ * its aliases intersects the inclusive base-grid limits; the canonical column may precede
+ * minTileColumn. Validates geometry and limits even for an invalid query. Does not wrap or mutate.
+ */
+export function resolveTileIndex(matrix, column, row, limits) {
+    validateTileMatrix(matrix);
+    if (limits !== undefined)
+        validateTileMatrixCoverage(matrix, limits);
+    if (!isTileIndexWithinMatrix(matrix, column, row))
+        return null;
+    const coalesce = getRowCoalesce(matrix, row);
+    const canonicalColumn = Math.floor(column / coalesce) * coalesce;
+    if (limits &&
+        (row < limits.minTileRow ||
+            row > limits.maxTileRow ||
+            canonicalColumn > limits.maxTileColumn ||
+            canonicalColumn + coalesce - 1 < limits.minTileColumn)) {
+        return null;
+    }
+    return [canonicalColumn, row];
+}
+/** Returns the entire tile footprint, including column aliases. Throws for out-of-range indices. */
+export function getTileBounds(matrix, column, row) {
+    if (!isTileIndexInRange(matrix, column, row))
+        throw new RangeError('Tile index outside matrix');
+    const width = matrix.resolution * matrix.tileSize[0];
+    const height = (matrix.resolutionY ?? matrix.resolution) * matrix.tileSize[1];
+    const coalesce = getRowCoalesce(matrix, row);
+    const canonicalColumn = Math.floor(column / coalesce) * coalesce;
+    const minimumX = matrix.origin[0] + canonicalColumn * width;
+    const minimumY = matrix.cornerOfOrigin === 'topLeft'
+        ? matrix.origin[1] - (row + 1) * height
+        : matrix.origin[1] + row * height;
+    return [
+        minimumX,
+        minimumY,
+        matrix.origin[0] + (canonicalColumn + coalesce) * width,
+        matrix.cornerOfOrigin === 'topLeft'
+            ? matrix.origin[1] - row * height
+            : matrix.origin[1] + (row + 1) * height
+    ];
+}
+/**
+ * Finds the containing tile without clamping or wrapping. Tile intervals are half-open in column
+ * and row direction: the origin edge belongs to the matrix, the opposite outer edge does not.
+ * Returns a canonical column for coalesced rows, or null for nonfinite/outside coordinates.
+ */
+export function getTileIndex(matrix, coordinate) {
+    validateTileMatrix(matrix);
+    if (!coordinate.every(Number.isFinite))
+        return null;
+    const row = getTileAxisIndex(matrix.origin[1], matrix.cornerOfOrigin === 'topLeft' ? -1 : 1, (matrix.resolutionY ?? matrix.resolution) * matrix.tileSize[1], matrix.matrixSize[1], coordinate[1]);
+    if (row === null)
+        return null;
+    const coalesce = getRowCoalesce(matrix, row);
+    const column = getTileAxisIndex(matrix.origin[0], 1, matrix.resolution * matrix.tileSize[0], matrix.matrixSize[0], coordinate[0], coalesce);
+    return column === null ? null : [column, row];
+}
+/** Finds half-open edge ownership using the same arithmetic as bounds, without division drift. */
+function getTileAxisIndex(origin, direction, tileSpan, matrixSize, coordinate, columnStep = 1) {
+    const end = origin + direction * (matrixSize * tileSpan);
+    if (direction === 1
+        ? coordinate < origin || coordinate >= end
+        : coordinate > origin || coordinate <= end) {
+        return null;
+    }
+    return ((findTileAxisBoundary(origin, direction, tileSpan, matrixSize, coordinate, false, columnStep) -
+        1) *
+        columnStep);
+}
+/** Finds the first edge beyond a coordinate, optionally including an equal edge. */
+function findTileAxisBoundary(origin, direction, tileSpan, matrixSize, coordinate, includeEqualEdge, columnStep = 1) {
+    // Binary search avoids enumerating edges or assuming division is accurate near an integer.
+    let lower = 0;
+    let upper = matrixSize / columnStep;
+    while (lower < upper) {
+        const middle = lower + Math.floor((upper - lower) / 2);
+        const edge = origin + direction * (middle * columnStep * tileSpan);
+        const precedesCoordinate = direction === 1
+            ? includeEqualEdge
+                ? edge < coordinate
+                : edge <= coordinate
+            : includeEqualEdge
+                ? edge > coordinate
+                : edge >= coordinate;
+        if (precedesCoordinate)
+            lower = middle + 1;
+        else
+            upper = middle;
+    }
+    return lower;
+}
+/** Validates inclusive coverage limits against the dimensions of an already validated matrix. */
+function validateTileMatrixCoverage(matrix, limits) {
+    if (!limits)
+        throw new RangeError('Invalid tile matrix limits');
+    const values = [limits.minTileColumn, limits.maxTileColumn, limits.minTileRow, limits.maxTileRow];
+    if (!values.every(value => Number.isSafeInteger(value) && value >= 0) ||
+        limits.minTileColumn > limits.maxTileColumn ||
+        limits.minTileRow > limits.maxTileRow ||
+        limits.maxTileColumn >= matrix.matrixSize[0] ||
+        limits.maxTileRow >= matrix.matrixSize[1]) {
+        throw new RangeError('Invalid tile matrix limits');
+    }
+}
+/**
+ * Returns a compact inclusive tile range intersecting an X/Y extent, clipped to the matrix and
+ * optional coverage limits, or null for no intersection. Positive-area queries exclude tiles
+ * that only touch an edge. Zero-width/height queries use point edge ownership on that axis.
+ * Bounds must be finite and ordered; no wrapping or coordinate transformation is performed.
+ */
+function getRegularTileRange(matrix, bounds, limits) {
+    validateTileQueryBounds(bounds);
+    const columns = getTileAxisRange(matrix.origin[0], 1, matrix.resolution * matrix.tileSize[0], matrix.matrixSize[0], bounds[0], bounds[2]);
+    const topLeft = matrix.cornerOfOrigin === 'topLeft';
+    const rows = getTileAxisRange(matrix.origin[1], topLeft ? -1 : 1, (matrix.resolutionY ?? matrix.resolution) * matrix.tileSize[1], matrix.matrixSize[1], topLeft ? bounds[3] : bounds[1], topLeft ? bounds[1] : bounds[3]);
+    if (!columns || !rows)
+        return null;
+    const range = {
+        minTileColumn: Math.max(columns[0], limits?.minTileColumn ?? 0),
+        maxTileColumn: Math.min(columns[1], limits?.maxTileColumn ?? matrix.matrixSize[0] - 1),
+        minTileRow: Math.max(rows[0], limits?.minTileRow ?? 0),
+        maxTileRow: Math.min(rows[1], limits?.maxTileRow ?? matrix.matrixSize[1] - 1)
+    };
+    return range.minTileColumn > range.maxTileColumn || range.minTileRow > range.maxTileRow
+        ? null
+        : range;
+}
+/** Rejects nonfinite or unordered query bounds independently of matrix validation. */
+function validateTileQueryBounds(bounds) {
+    if (bounds.length !== 4 ||
+        !bounds.every(Number.isFinite) ||
+        bounds[0] > bounds[2] ||
+        bounds[1] > bounds[3]) {
+        throw new RangeError('Tile query bounds must be finite and ordered');
+    }
+}
+/**
+ * Returns a rectangular range for a regular matrix. Variable-width matrices require getTileRanges
+ * so callers cannot accidentally enumerate aliases as separate tiles.
+ */
+export function getTileRange(matrix, bounds, limits) {
+    validateTileMatrix(matrix);
+    if (matrix.variableMatrixWidths?.length) {
+        throw new RangeError('Use getTileRanges for variable-width matrices');
+    }
+    if (limits)
+        validateTileMatrixCoverage(matrix, limits);
+    return getRegularTileRange(matrix, bounds, limits);
+}
+/** Returns the coalescence factor for an already validated row. */
+function getRowCoalesce(matrix, row) {
+    return (matrix.variableMatrixWidths?.find(band => row >= band.minTileRow && row <= band.maxTileRow)
+        ?.coalesce ?? 1);
+}
+/**
+ * Returns exact canonical-column ranges, ordered by native row, split at coalescence changes.
+ * Columns advance by columnStep, rows by one. Limits use base-grid indices: a tile is included
+ * when any of its aliases falls within those limits. Canonical columns may precede minTileColumn.
+ * Work and output scale with row bands rather than the number of tiles or rows.
+ */
+export function getTileRanges(matrix, bounds, limits) {
+    const declaredBands = validateTileMatrixGeometry(matrix);
+    if (limits)
+        validateTileMatrixCoverage(matrix, limits);
+    validateTileQueryBounds(bounds);
+    const bands = [];
+    let nextRow = 0;
+    for (const band of declaredBands) {
+        if (nextRow < band.minTileRow) {
+            bands.push({ coalesce: 1, minTileRow: nextRow, maxTileRow: band.minTileRow - 1 });
+        }
+        bands.push(band);
+        nextRow = band.maxTileRow + 1;
+    }
+    if (nextRow < matrix.matrixSize[1]) {
+        bands.push({ coalesce: 1, minTileRow: nextRow, maxTileRow: matrix.matrixSize[1] - 1 });
+    }
+    const topLeft = matrix.cornerOfOrigin === 'topLeft';
+    const rows = getTileAxisRange(matrix.origin[1], topLeft ? -1 : 1, (matrix.resolutionY ?? matrix.resolution) * matrix.tileSize[1], matrix.matrixSize[1], topLeft ? bounds[3] : bounds[1], topLeft ? bounds[1] : bounds[3]);
+    if (!rows)
+        return [];
+    const ranges = [];
+    for (const band of bands) {
+        const minimumRow = Math.max(band.minTileRow, rows[0], limits?.minTileRow ?? 0);
+        const maximumRow = Math.min(band.maxTileRow, rows[1], limits?.maxTileRow ?? matrix.matrixSize[1] - 1);
+        if (minimumRow > maximumRow)
+            continue;
+        const columns = getTileAxisRange(matrix.origin[0], 1, matrix.resolution * matrix.tileSize[0], matrix.matrixSize[0], bounds[0], bounds[2], band.coalesce);
+        if (!columns)
+            continue;
+        const minimumColumn = Math.max(columns[0], Math.floor((limits?.minTileColumn ?? 0) / band.coalesce) * band.coalesce);
+        const maximumColumn = Math.min(columns[1], Math.floor((limits?.maxTileColumn ?? matrix.matrixSize[0] - 1) / band.coalesce) *
+            band.coalesce);
+        if (minimumColumn <= maximumColumn)
+            ranges.push({
+                minTileColumn: minimumColumn,
+                maxTileColumn: maximumColumn,
+                minTileRow: minimumRow,
+                maxTileRow: maximumRow,
+                columnStep: band.coalesce
+            });
+    }
+    return ranges;
+}
+/** Finds an inclusive range along one directed axis, clipping before comparing computed edges. */
+function getTileAxisRange(origin, direction, tileSpan, matrixSize, start, end, columnStep = 1) {
+    if (start === end) {
+        const index = getTileAxisIndex(origin, direction, tileSpan, matrixSize, start, columnStep);
+        return index === null ? null : [index, index];
+    }
+    const matrixEnd = origin + direction * (matrixSize * tileSpan);
+    const clippedStart = direction === 1 ? Math.max(start, origin) : Math.min(start, origin);
+    const clippedEnd = direction === 1 ? Math.min(end, matrixEnd) : Math.max(end, matrixEnd);
+    if (direction === 1 ? clippedStart >= clippedEnd : clippedStart <= clippedEnd)
+        return null;
+    return [
+        (findTileAxisBoundary(origin, direction, tileSpan, matrixSize, clippedStart, false, columnStep) -
+            1) *
+            columnStep,
+        (findTileAxisBoundary(origin, direction, tileSpan, matrixSize, clippedEnd, true, columnStep) -
+            1) *
+            columnStep
+    ];
+}
+/**
+ * Selects the coarsest matrix with resolution at or below the positive finite target, or the
+ * finest available matrix when all levels are coarser. Target units must match the matrix set's
+ * coordinate units per pixel. Returns null for an empty set; ties retain the first matrix.
+ * Uses nominal uncoalesced resolution; selection is independent of row.
+ * Uses the larger axis resolution for anisotropic pixels, so both axes meet the target.
+ * Validates the complete set, including its CRS identifier. Does not sort or mutate the set.
+ */
+export function selectTileMatrix(matrixSet, targetResolution) {
+    if (!Number.isFinite(targetResolution) || targetResolution <= 0) {
+        throw new RangeError('Target resolution must be positive and finite');
+    }
+    let selected = null;
+    let finest = null;
+    let selectedResolution = -Infinity;
+    let finestResolution = Infinity;
+    validateTileMatrixSet(matrixSet);
+    for (const matrix of matrixSet.matrices) {
+        const resolution = Math.max(matrix.resolution, matrix.resolutionY ?? matrix.resolution);
+        if (!finest || resolution < finestResolution) {
+            finest = matrix;
+            finestResolution = resolution;
+        }
+        if (resolution <= targetResolution && (!selected || resolution > selectedResolution)) {
+            selected = matrix;
+            selectedResolution = resolution;
+        }
+    }
+    return selected ?? finest;
+}
+//# sourceMappingURL=tile-matrix.js.map

package/dist/tile-matrix.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"tile-matrix.js","sourceRoot":"","sources":["../src/tile-matrix.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAmE5D,wEAAwE;AACxE,MAAM,UAAU,kBAAkB,CAAC,MAAkB;IACnD,0BAA0B,CAAC,MAAM,CAAC,CAAC;AACrC,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,qBAAqB,CACnC,SAAwB,EACxB,MAA4B;IAE5B,IACE,CAAC,SAAS;QACV,OAAO,SAAS,CAAC,GAAG,KAAK,QAAQ;QACjC,CAAC,SAAS,CAAC,GAAG,CAAC,IAAI,EAAE;QACrB,CAAC,KAAK,CAAC,OAAO,CAAC,SAAS,CAAC,QAAQ,CAAC,EAClC,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,yBAAyB,CAAC,CAAC;IAClD,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAsB,CAAC;IAC/C,KAAK,MAAM,MAAM,IAAI,SAAS,CAAC,QAAQ,EAAE,CAAC;QACxC,kBAAkB,CAAC,MAAM,CAAC,CAAC;QAC3B,IAAI,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,wCAAwC,CAAC,CAAC;QAC5F,QAAQ,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,EAAE,MAAM,CAAC,CAAC;IAClC,CAAC;IACD,IAAI,MAAM,KAAK,SAAS;QAAE,OAAO;IACjC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,yCAAyC,CAAC,CAAC;IAC5F,MAAM,WAAW,GAAG,IAAI,GAAG,EAAU,CAAC;IACtC,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,IACE,CAAC,KAAK;YACN,OAAO,KAAK,CAAC,UAAU,KAAK,QAAQ;YACpC,CAAC,KAAK,CAAC,UAAU;YACjB,WAAW,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,EACjC,CAAC;YACD,MAAM,IAAI,UAAU,CAAC,0DAA0D,CAAC,CAAC;QACnF,CAAC;QACD,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC;QAC9C,IAAI,CAAC,MAAM;YAAE,MAAM,IAAI,UAAU,CAAC,gDAAgD,CAAC,CAAC;QACpF,0BAA0B,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAC1C,WAAW,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC;IACpC,CAAC;AACH,CAAC;AAED,mGAAmG;AACnG,MAAM,UAAU,wBAAwB,CAAC,MAAkB,EAAE,MAAwB;IACnF,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC3B,0BAA0B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;AAC7C,CAAC;AAED,sFAAsF;AACtF,SAAS,0BAA0B,CAAC,MAAkB;IACpD,IACE,CAAC,MAAM;QACP,OAAO,MAAM,CAAC,EAAE,KAAK,QAAQ;QAC7B,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC;QAC7B,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC;QAC/B,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC;QACjC,CAAC,MAAM,CAAC,oBAAoB,KAAK,SAAS,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,oBAAoB,CAAC,CAAC;QAC1F,CAAC,MAAM,CAAC,EAAE;QACV,MAAM,CAAC,MAAM,CAAC,MAAM,KAAK,CAAC;QAC1B,MAAM,CAAC,QAAQ,CAAC,MAAM,KAAK,CAAC;QAC5B,MAAM,CAAC,UAAU,CAAC,MAAM,KAAK,CAAC;QAC9B,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,UAAU,CAAC;QACnC,MAAM,CAAC,UAAU,IAAI,CAAC;QACtB,CAAC,MAAM,CAAC,WAAW,KAAK,SAAS;YAC/B,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,MAAM,CAAC,WAAW,IAAI,CAAC,CAAC,CAAC;QACpE,CAAC,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QACrC,CAAC,CAAC,GAAG,MAAM,CAAC,QAAQ,EAAE,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,KAAK,CAC/C,KAAK,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC,CAClD;QACD,CAAC,CAAC,SAAS,EAAE,YAAY,CAAC,CAAC,QAAQ,CAAC,MAAM,CAAC,cAAc,CAAC,EAC1D,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,8BAA8B,CAAC,CAAC;IACvD,CAAC;IACD,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,oBAAoB,IAAI,EAAE,CAAC,CAAC,CAAC;IACvD,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,+BAA+B,CAAC,CAAC;IACrF,KAAK,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,MAAM,EAAE,EAAE,CAAC,KAAK,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC;IACpE,IAAI,WAAW,GAAG,CAAC,CAAC,CAAC;IACrB,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IACE,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,QAAQ,CAAC;YACpC,IAAI,CAAC,QAAQ,IAAI,CAAC;YAClB,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,KAAK,CAAC;YAC1C,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,UAAU,CAAC;YACtC,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,UAAU,CAAC;YACtC,IAAI,CAAC,UAAU,IAAI,WAAW;YAC9B,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,UAAU;YACjC,IAAI,CAAC,UAAU,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;YAEvC,MAAM,IAAI,UAAU,CAAC,+BAA+B,CAAC,CAAC;QACxD,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC;IAChC,CAAC;IACD,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QACpC,MAAM,UAAU,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,UAAU,CAAC;QAC9F,MAAM,QAAQ,GAAG,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACpD,MAAM,IAAI,GAAG,QAAQ,GAAG,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;QAChD,MAAM,SAAS,GAAG,IAAI,KAAK,CAAC,IAAI,MAAM,CAAC,cAAc,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7E,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,SAAS,GAAG,IAAI,CAAC;QACnD,IACE,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;YACtB,IAAI,IAAI,CAAC;YACT,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;YACrB,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,SAAS,GAAG,QAAQ,KAAK,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC;YAClE,GAAG,GAAG,SAAS,GAAG,QAAQ,KAAK,GAAG,EAClC,CAAC;YACD,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,sGAAsG;AACtG,MAAM,UAAU,kBAAkB,CAChC,MAAkB,EAClB,MAAc,EACd,GAAW,EACX,MAAyB;IAEzB,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC3B,IAAI,MAAM;QAAE,0BAA0B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACvD,OAAO,CACL,uBAAuB,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,CAAC;QAC5C,CAAC,CAAC,MAAM;YACN,CAAC,MAAM,IAAI,MAAM,CAAC,aAAa;gBAC7B,MAAM,IAAI,MAAM,CAAC,aAAa;gBAC9B,GAAG,IAAI,MAAM,CAAC,UAAU;gBACxB,GAAG,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAC/B,CAAC;AACJ,CAAC;AAED,oEAAoE;AACpE,SAAS,uBAAuB,CAAC,MAAkB,EAAE,MAAc,EAAE,GAAW;IAC9E,OAAO,CACL,MAAM,CAAC,aAAa,CAAC,MAAM,CAAC;QAC5B,MAAM,CAAC,aAAa,CAAC,GAAG,CAAC;QACzB,MAAM,IAAI,CAAC;QACX,GAAG,IAAI,CAAC;QACR,MAAM,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;QAC7B,GAAG,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAC3B,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,gBAAgB,CAC9B,MAAkB,EAClB,MAAc,EACd,GAAW,EACX,MAAyB;IAEzB,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC3B,IAAI,MAAM,KAAK,SAAS;QAAE,0BAA0B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACrE,IAAI,CAAC,uBAAuB,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,CAAC;QAAE,OAAO,IAAI,CAAC;IAC/D,MAAM,QAAQ,GAAG,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAC7C,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,QAAQ,CAAC;IACjE,IACE,MAAM;QACN,CAAC,GAAG,GAAG,MAAM,CAAC,UAAU;YACtB,GAAG,GAAG,MAAM,CAAC,UAAU;YACvB,eAAe,GAAG,MAAM,CAAC,aAAa;YACtC,eAAe,GAAG,QAAQ,GAAG,CAAC,GAAG,MAAM,CAAC,aAAa,CAAC,EACxD,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,CAAC,eAAe,EAAE,GAAG,CAAC,CAAC;AAChC,CAAC;AAED,oGAAoG;AACpG,MAAM,UAAU,aAAa,CAAC,MAAkB,EAAE,MAAc,EAAE,GAAW;IAC3E,IAAI,CAAC,kBAAkB,CAAC,MAAM,EAAE,MAAM,EAAE,GAAG,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,2BAA2B,CAAC,CAAC;IAChG,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACrD,MAAM,MAAM,GAAG,CAAC,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IAC9E,MAAM,QAAQ,GAAG,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAC7C,MAAM,eAAe,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,QAAQ,CAAC;IACjE,MAAM,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,eAAe,GAAG,KAAK,CAAC;IAC5D,MAAM,QAAQ,GACZ,MAAM,CAAC,cAAc,KAAK,SAAS;QACjC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM;QACvC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,MAAM,CAAC;IACtC,OAAO;QACL,QAAQ;QACR,QAAQ;QACR,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,eAAe,GAAG,QAAQ,CAAC,GAAG,KAAK;QACvD,MAAM,CAAC,cAAc,KAAK,SAAS;YACjC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,MAAM;YACjC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM;KAC1C,CAAC;AACJ,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,YAAY,CAC1B,MAAkB,EAClB,UAAqC;IAErC,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC3B,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,IAAI,CAAC;IACpD,MAAM,GAAG,GAAG,gBAAgB,CAC1B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,MAAM,CAAC,cAAc,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAC5C,CAAC,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAC9D,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,UAAU,CAAC,CAAC,CAAC,CACd,CAAC;IACF,IAAI,GAAG,KAAK,IAAI;QAAE,OAAO,IAAI,CAAC;IAC9B,MAAM,QAAQ,GAAG,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAC7C,MAAM,MAAM,GAAG,gBAAgB,CAC7B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,CAAC,EACD,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EACtC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,UAAU,CAAC,CAAC,CAAC,EACb,QAAQ,CACT,CAAC;IACF,OAAO,MAAM,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;AAChD,CAAC;AAED,kGAAkG;AAClG,SAAS,gBAAgB,CACvB,MAAc,EACd,SAAiB,EACjB,QAAgB,EAChB,UAAkB,EAClB,UAAkB,EAClB,UAAU,GAAG,CAAC;IAEd,MAAM,GAAG,GAAG,MAAM,GAAG,SAAS,GAAG,CAAC,UAAU,GAAG,QAAQ,CAAC,CAAC;IACzD,IACE,SAAS,KAAK,CAAC;QACb,CAAC,CAAC,UAAU,GAAG,MAAM,IAAI,UAAU,IAAI,GAAG;QAC1C,CAAC,CAAC,UAAU,GAAG,MAAM,IAAI,UAAU,IAAI,GAAG,EAC5C,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,OAAO,CACL,CAAC,oBAAoB,CAAC,MAAM,EAAE,SAAS,EAAE,QAAQ,EAAE,UAAU,EAAE,UAAU,EAAE,KAAK,EAAE,UAAU,CAAC;QAC3F,CAAC,CAAC;QACJ,UAAU,CACX,CAAC;AACJ,CAAC;AAED,oFAAoF;AACpF,SAAS,oBAAoB,CAC3B,MAAc,EACd,SAAiB,EACjB,QAAgB,EAChB,UAAkB,EAClB,UAAkB,EAClB,gBAAyB,EACzB,UAAU,GAAG,CAAC;IAEd,2FAA2F;IAC3F,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,KAAK,GAAG,UAAU,GAAG,UAAU,CAAC;IACpC,OAAO,KAAK,GAAG,KAAK,EAAE,CAAC;QACrB,MAAM,MAAM,GAAG,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACvD,MAAM,IAAI,GAAG,MAAM,GAAG,SAAS,GAAG,CAAC,MAAM,GAAG,UAAU,GAAG,QAAQ,CAAC,CAAC;QACnE,MAAM,kBAAkB,GACtB,SAAS,KAAK,CAAC;YACb,CAAC,CAAC,gBAAgB;gBAChB,CAAC,CAAC,IAAI,GAAG,UAAU;gBACnB,CAAC,CAAC,IAAI,IAAI,UAAU;YACtB,CAAC,CAAC,gBAAgB;gBAChB,CAAC,CAAC,IAAI,GAAG,UAAU;gBACnB,CAAC,CAAC,IAAI,IAAI,UAAU,CAAC;QAC3B,IAAI,kBAAkB;YAAE,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;;YACtC,KAAK,GAAG,MAAM,CAAC;IACtB,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,iGAAiG;AACjG,SAAS,0BAA0B,CAAC,MAAkB,EAAE,MAAwB;IAC9E,IAAI,CAAC,MAAM;QAAE,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;IAChE,MAAM,MAAM,GAAG,CAAC,MAAM,CAAC,aAAa,EAAE,MAAM,CAAC,aAAa,EAAE,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC;IAClG,IACE,CAAC,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC;QACjE,MAAM,CAAC,aAAa,GAAG,MAAM,CAAC,aAAa;QAC3C,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,UAAU;QACrC,MAAM,CAAC,aAAa,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;QAC5C,MAAM,CAAC,UAAU,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACzC,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,4BAA4B,CAAC,CAAC;IACrD,CAAC;AACH,CAAC;AAED;;;;;GAKG;AACH,SAAS,mBAAmB,CAC1B,MAAkB,EAClB,MAAwB,EACxB,MAAyB;IAEzB,uBAAuB,CAAC,MAAM,CAAC,CAAC;IAChC,MAAM,OAAO,GAAG,gBAAgB,CAC9B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,CAAC,EACD,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EACtC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,MAAM,CAAC,CAAC,CAAC,EACT,MAAM,CAAC,CAAC,CAAC,CACV,CAAC;IACF,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,KAAK,SAAS,CAAC;IACpD,MAAM,IAAI,GAAG,gBAAgB,CAC3B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAChB,CAAC,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAC9D,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAC/B,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAChC,CAAC;IACF,IAAI,CAAC,OAAO,IAAI,CAAC,IAAI;QAAE,OAAO,IAAI,CAAC;IACnC,MAAM,KAAK,GAAG;QACZ,aAAa,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,aAAa,IAAI,CAAC,CAAC;QAC/D,aAAa,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,aAAa,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACtF,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,IAAI,CAAC,CAAC;QACtD,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;KAC9E,CAAC;IACF,OAAO,KAAK,CAAC,aAAa,GAAG,KAAK,CAAC,aAAa,IAAI,KAAK,CAAC,UAAU,GAAG,KAAK,CAAC,UAAU;QACrF,CAAC,CAAC,IAAI;QACN,CAAC,CAAC,KAAK,CAAC;AACZ,CAAC;AAED,sFAAsF;AACtF,SAAS,uBAAuB,CAAC,MAAwB;IACvD,IACE,MAAM,CAAC,MAAM,KAAK,CAAC;QACnB,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QACrB,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EACrB,CAAC;QACD,MAAM,IAAI,UAAU,CAAC,8CAA8C,CAAC,CAAC;IACvE,CAAC;AACH,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,YAAY,CAC1B,MAAkB,EAClB,MAAwB,EACxB,MAAyB;IAEzB,kBAAkB,CAAC,MAAM,CAAC,CAAC;IAC3B,IAAI,MAAM,CAAC,oBAAoB,EAAE,MAAM,EAAE,CAAC;QACxC,MAAM,IAAI,UAAU,CAAC,+CAA+C,CAAC,CAAC;IACxE,CAAC;IACD,IAAI,MAAM;QAAE,0BAA0B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACvD,OAAO,mBAAmB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;AACrD,CAAC;AAED,mEAAmE;AACnE,SAAS,cAAc,CAAC,MAAkB,EAAE,GAAW;IACrD,OAAO,CACL,MAAM,CAAC,oBAAoB,EAAE,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,IAAI,CAAC,UAAU,IAAI,GAAG,IAAI,IAAI,CAAC,UAAU,CAAC;QACzF,EAAE,QAAQ,IAAI,CAAC,CAClB,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,aAAa,CAC3B,MAAkB,EAClB,MAAwB,EACxB,MAAyB;IAEzB,MAAM,aAAa,GAAG,0BAA0B,CAAC,MAAM,CAAC,CAAC;IACzD,IAAI,MAAM;QAAE,0BAA0B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACvD,uBAAuB,CAAC,MAAM,CAAC,CAAC;IAChC,MAAM,KAAK,GAA0B,EAAE,CAAC;IACxC,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,KAAK,MAAM,IAAI,IAAI,aAAa,EAAE,CAAC;QACjC,IAAI,OAAO,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC;YAC9B,KAAK,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,CAAC,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,IAAI,CAAC,UAAU,GAAG,CAAC,EAAC,CAAC,CAAC;QAClF,CAAC;QACD,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACjB,OAAO,GAAG,IAAI,CAAC,UAAU,GAAG,CAAC,CAAC;IAChC,CAAC;IACD,IAAI,OAAO,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC;QACnC,KAAK,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,CAAC,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,EAAC,CAAC,CAAC;IACvF,CAAC;IACD,MAAM,OAAO,GAAG,MAAM,CAAC,cAAc,KAAK,SAAS,CAAC;IACpD,MAAM,IAAI,GAAG,gBAAgB,CAC3B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAChB,CAAC,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAC9D,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAC/B,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAChC,CAAC;IACF,IAAI,CAAC,IAAI;QAAE,OAAO,EAAE,CAAC;IACrB,MAAM,MAAM,GAAyB,EAAE,CAAC;IACxC,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,UAAU,IAAI,CAAC,CAAC,CAAC;QAC/E,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CACzB,IAAI,CAAC,UAAU,EACf,IAAI,CAAC,CAAC,CAAC,EACP,MAAM,EAAE,UAAU,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAC/C,CAAC;QACF,IAAI,UAAU,GAAG,UAAU;YAAE,SAAS;QACtC,MAAM,OAAO,GAAG,gBAAgB,CAC9B,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAChB,CAAC,EACD,MAAM,CAAC,UAAU,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EACtC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EACpB,MAAM,CAAC,CAAC,CAAC,EACT,MAAM,CAAC,CAAC,CAAC,EACT,IAAI,CAAC,QAAQ,CACd,CAAC;QACF,IAAI,CAAC,OAAO;YAAE,SAAS;QACvB,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,OAAO,CAAC,CAAC,CAAC,EACV,IAAI,CAAC,KAAK,CAAC,CAAC,MAAM,EAAE,aAAa,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,QAAQ,CACzE,CAAC;QACF,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAC5B,OAAO,CAAC,CAAC,CAAC,EACV,IAAI,CAAC,KAAK,CAAC,CAAC,MAAM,EAAE,aAAa,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;YAC7E,IAAI,CAAC,QAAQ,CAChB,CAAC;QACF,IAAI,aAAa,IAAI,aAAa;YAChC,MAAM,CAAC,IAAI,CAAC;gBACV,aAAa,EAAE,aAAa;gBAC5B,aAAa,EAAE,aAAa;gBAC5B,UAAU,EAAE,UAAU;gBACtB,UAAU,EAAE,UAAU;gBACtB,UAAU,EAAE,IAAI,CAAC,QAAQ;aAC1B,CAAC,CAAC;IACP,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,kGAAkG;AAClG,SAAS,gBAAgB,CACvB,MAAc,EACd,SAAiB,EACjB,QAAgB,EAChB,UAAkB,EAClB,KAAa,EACb,GAAW,EACX,UAAU,GAAG,CAAC;IAEd,IAAI,KAAK,KAAK,GAAG,EAAE,CAAC;QAClB,MAAM,KAAK,GAAG,gBAAgB,CAAC,MAAM,EAAE,SAAS,EAAE,QAAQ,EAAE,UAAU,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC;QAC3F,OAAO,KAAK,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IAChD,CAAC;IACD,MAAM,SAAS,GAAG,MAAM,GAAG,SAAS,GAAG,CAAC,UAAU,GAAG,QAAQ,CAAC,CAAC;IAC/D,MAAM,YAAY,GAAG,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;IACzF,MAAM,UAAU,GAAG,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;IACzF,IAAI,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,YAAY,IAAI,UAAU,CAAC,CAAC,CAAC,YAAY,IAAI,UAAU;QAAE,OAAO,IAAI,CAAC;IAC3F,OAAO;QACL,CAAC,oBAAoB,CACnB,MAAM,EACN,SAAS,EACT,QAAQ,EACR,UAAU,EACV,YAAY,EACZ,KAAK,EACL,UAAU,CACX;YACC,CAAC,CAAC;YACF,UAAU;QACZ,CAAC,oBAAoB,CAAC,MAAM,EAAE,SAAS,EAAE,QAAQ,EAAE,UAAU,EAAE,UAAU,EAAE,IAAI,EAAE,UAAU,CAAC;YAC1F,CAAC,CAAC;YACF,UAAU;KACb,CAAC;AACJ,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,gBAAgB,CAC9B,SAAwB,EACxB,gBAAwB;IAExB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,gBAAgB,IAAI,CAAC,EAAE,CAAC;QAChE,MAAM,IAAI,UAAU,CAAC,+CAA+C,CAAC,CAAC;IACxE,CAAC;IACD,IAAI,QAAQ,GAAsB,IAAI,CAAC;IACvC,IAAI,MAAM,GAAsB,IAAI,CAAC;IACrC,IAAI,kBAAkB,GAAG,CAAC,QAAQ,CAAC;IACnC,IAAI,gBAAgB,GAAG,QAAQ,CAAC;IAChC,qBAAqB,CAAC,SAAS,CAAC,CAAC;IACjC,KAAK,MAAM,MAAM,IAAI,SAAS,CAAC,QAAQ,EAAE,CAAC;QACxC,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,EAAE,MAAM,CAAC,WAAW,IAAI,MAAM,CAAC,UAAU,CAAC,CAAC;QACxF,IAAI,CAAC,MAAM,IAAI,UAAU,GAAG,gBAAgB,EAAE,CAAC;YAC7C,MAAM,GAAG,MAAM,CAAC;YAChB,gBAAgB,GAAG,UAAU,CAAC;QAChC,CAAC;QACD,IAAI,UAAU,IAAI,gBAAgB,IAAI,CAAC,CAAC,QAAQ,IAAI,UAAU,GAAG,kBAAkB,CAAC,EAAE,CAAC;YACrF,QAAQ,GAAG,MAAM,CAAC;YAClB,kBAAkB,GAAG,UAAU,CAAC;QAClC,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,IAAI,MAAM,CAAC;AAC5B,CAAC"}

package/dist/tiles.d.ts

@@ -0,0 +1,5 @@
+export { getGeographicTile, getGeographicTileBounds, getGeographicTileRanges } from "./geographic-tiles.js";
+export type { GeographicTile, GeographicTileRange } from "./geographic-tiles.js";
+export { validateTileMatrix, validateTileMatrixSet, validateTileMatrixLimits, resolveTileIndex, isTileIndexInRange, getTileBounds, getTileIndex, getTileRange, getTileRanges, selectTileMatrix } from "./tile-matrix.js";
+export type { TileMatrix, TileMatrixSet, TileMatrixLimits, TileMatrixSetLimits, TileMatrixBounds, VariableMatrixWidth, TileMatrixRowRange } from "./tile-matrix.js";
+//# sourceMappingURL=tiles.d.ts.map

package/dist/tiles.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"tiles.d.ts","sourceRoot":"","sources":["../src/tiles.ts"],"names":[],"mappings":"AAIA,OAAO,EACL,iBAAiB,EACjB,uBAAuB,EACvB,uBAAuB,EACxB,8BAA2B;AAC5B,YAAY,EAAC,cAAc,EAAE,mBAAmB,EAAC,8BAA2B;AAE5E,OAAO,EACL,kBAAkB,EAClB,qBAAqB,EACrB,wBAAwB,EACxB,gBAAgB,EAChB,kBAAkB,EAClB,aAAa,EACb,YAAY,EACZ,YAAY,EACZ,aAAa,EACb,gBAAgB,EACjB,yBAAsB;AACvB,YAAY,EACV,UAAU,EACV,aAAa,EACb,gBAAgB,EAChB,mBAAmB,EACnB,gBAAgB,EAChB,mBAAmB,EACnB,kBAAkB,EACnB,yBAAsB"}

package/dist/tiles.js

@@ -0,0 +1,6 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+export { getGeographicTile, getGeographicTileBounds, getGeographicTileRanges } from "./geographic-tiles.js";
+export { validateTileMatrix, validateTileMatrixSet, validateTileMatrixLimits, resolveTileIndex, isTileIndexInRange, getTileBounds, getTileIndex, getTileRange, getTileRanges, selectTileMatrix } from "./tile-matrix.js";
+//# sourceMappingURL=tiles.js.map

package/dist/tiles.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"tiles.js","sourceRoot":"","sources":["../src/tiles.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAE5D,OAAO,EACL,iBAAiB,EACjB,uBAAuB,EACvB,uBAAuB,EACxB,8BAA2B;AAG5B,OAAO,EACL,kBAAkB,EAClB,qBAAqB,EACrB,wBAAwB,EACxB,gBAAgB,EAChB,kBAAkB,EAClB,aAAa,EACb,YAAY,EACZ,YAAY,EACZ,aAAa,EACb,gBAAgB,EACjB,yBAAsB"}

package/dist/triaxial/integrate.d.ts

@@ -0,0 +1,2 @@
+export declare function integrateGeodesic(initial: readonly number[], q: readonly number[], distance: number): number[];
+//# sourceMappingURL=integrate.d.ts.map

package/dist/triaxial/integrate.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"integrate.d.ts","sourceRoot":"","sources":["../../src/triaxial/integrate.ts"],"names":[],"mappings":"AAwBA,wBAAgB,iBAAiB,CAC/B,OAAO,EAAE,SAAS,MAAM,EAAE,EAC1B,CAAC,EAAE,SAAS,MAAM,EAAE,EACpB,QAAQ,EAAE,MAAM,GACf,MAAM,EAAE,CAqDV"}

package/dist/triaxial/integrate.js

@@ -0,0 +1,76 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adaptive Dormand-Prince integration of the constrained Cartesian geodesic equation; mathematical coefficients, no third-party implementation copied.
+const C = [0, 1 / 5, 3 / 10, 4 / 5, 8 / 9, 1, 1];
+const A = [
+    [],
+    [1 / 5],
+    [3 / 40, 9 / 40],
+    [44 / 45, -56 / 15, 32 / 9],
+    [19372 / 6561, -25360 / 2187, 64448 / 6561, -212 / 729],
+    [9017 / 3168, -355 / 33, 46732 / 5247, 49 / 176, -5103 / 18656],
+    [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84]
+];
+const B = [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84, 0];
+const ERROR = [
+    35 / 384 - 5179 / 57600,
+    0,
+    500 / 1113 - 7571 / 16695,
+    125 / 192 - 393 / 640,
+    -2187 / 6784 + 92097 / 339200,
+    11 / 84 - 187 / 2100,
+    -1 / 40
+];
+export function integrateGeodesic(initial, q, distance) {
+    let state = Array.from(initial), done = 0, step = Math.sign(distance) * Math.min(0.01, Math.abs(distance));
+    const derivative = (y) => {
+        const normal = [q[0] * y[0], q[1] * y[1], q[2] * y[2]];
+        const curvature = (q[0] * y[3] ** 2 + q[1] * y[4] ** 2 + q[2] * y[5] ** 2) /
+            (normal[0] ** 2 + normal[1] ** 2 + normal[2] ** 2);
+        return [
+            y[3],
+            y[4],
+            y[5],
+            -curvature * normal[0],
+            -curvature * normal[1],
+            -curvature * normal[2]
+        ];
+    };
+    for (let attempt = 0; attempt < 50000; attempt++) {
+        if (done === distance)
+            return state;
+        if (Math.abs(step) > Math.abs(distance - done))
+            step = distance - done;
+        const stages = [];
+        for (let j = 0; j < C.length; j++) {
+            const y = state.map((v, i) => v + step * A[j].reduce((sum, a, k) => sum + a * stages[k][i], 0));
+            stages.push(derivative(y));
+        }
+        const next = state.map((v, i) => v + step * B.reduce((sum, b, j) => sum + b * stages[j][i], 0));
+        let error = 0;
+        for (let i = 0; i < 6; i++)
+            error = Math.max(error, Math.abs(step * ERROR.reduce((sum, e, j) => sum + e * stages[j][i], 0)));
+        const tolerance = 2e-15;
+        if (error <= tolerance) {
+            state = next;
+            done += step;
+            // Restore the two constraints and unit speed only at accepted steps.
+            const surface = Math.sqrt(q[0] * state[0] ** 2 + q[1] * state[1] ** 2 + q[2] * state[2] ** 2);
+            for (let i = 0; i < 3; i++)
+                state[i] /= surface;
+            const normal = state.slice(0, 3).map((x, i) => q[i] * x), nn = normal.reduce((s, x) => s + x * x, 0);
+            const dot = normal.reduce((s, x, i) => s + x * state[i + 3], 0) / nn;
+            for (let i = 0; i < 3; i++)
+                state[i + 3] -= normal[i] * dot;
+            const speed = Math.hypot(state[3], state[4], state[5]);
+            for (let i = 3; i < 6; i++)
+                state[i] /= speed;
+        }
+        step *= Math.max(0.2, Math.min(4, 0.9 * (tolerance / Math.max(error, 1e-300)) ** 0.2));
+        if (done + step === done)
+            throw new RangeError('Triaxial integration step underflow');
+    }
+    throw new RangeError('Triaxial integration work limit exceeded');
+}
+//# sourceMappingURL=integrate.js.map

package/dist/triaxial/integrate.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"integrate.js","sourceRoot":"","sources":["../../src/triaxial/integrate.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,kLAAkL;AAClL,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACjD,MAAM,CAAC,GAAG;IACR,EAAE;IACF,CAAC,CAAC,GAAG,CAAC,CAAC;IACP,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC;IAChB,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;IAC3B,CAAC,KAAK,GAAG,IAAI,EAAE,CAAC,KAAK,GAAG,IAAI,EAAE,KAAK,GAAG,IAAI,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC;IACvD,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,KAAK,GAAG,IAAI,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,KAAK,CAAC;IAC/D,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,CAAC;CAC5D,CAAC;AACF,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,IAAI,GAAG,IAAI,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;AACzE,MAAM,KAAK,GAAG;IACZ,EAAE,GAAG,GAAG,GAAG,IAAI,GAAG,KAAK;IACvB,CAAC;IACD,GAAG,GAAG,IAAI,GAAG,IAAI,GAAG,KAAK;IACzB,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG;IACrB,CAAC,IAAI,GAAG,IAAI,GAAG,KAAK,GAAG,MAAM;IAC7B,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,IAAI;IACpB,CAAC,CAAC,GAAG,EAAE;CACR,CAAC;AACF,MAAM,UAAU,iBAAiB,CAC/B,OAA0B,EAC1B,CAAoB,EACpB,QAAgB;IAEhB,IAAI,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,EAC7B,IAAI,GAAG,CAAC,EACR,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;IAClE,MAAM,UAAU,GAAG,CAAC,CAAoB,EAAE,EAAE;QAC1C,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvD,MAAM,SAAS,GACb,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YACxD,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;QACrD,OAAO;YACL,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,CAAC,CAAC,CAAC;YACJ,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;YACtB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;YACtB,CAAC,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC;SACvB,CAAC;IACJ,CAAC,CAAC;IACF,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,KAAK,EAAE,OAAO,EAAE,EAAE,CAAC;QACjD,IAAI,IAAI,KAAK,QAAQ;YAAE,OAAO,KAAK,CAAC;QACpC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC;YAAE,IAAI,GAAG,QAAQ,GAAG,IAAI,CAAC;QACvE,MAAM,MAAM,GAAe,EAAE,CAAC;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CACjB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAC3E,CAAC;YACF,MAAM,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAChG,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,KAAK,GAAG,IAAI,CAAC,GAAG,CACd,KAAK,EACL,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CACxE,CAAC;QACJ,MAAM,SAAS,GAAG,KAAK,CAAC;QACxB,IAAI,KAAK,IAAI,SAAS,EAAE,CAAC;YACvB,KAAK,GAAG,IAAI,CAAC;YACb,IAAI,IAAI,IAAI,CAAC;YACb,qEAAqE;YACrE,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;YAC9F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC;YAChD,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EACtD,EAAE,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7C,MAAM,GAAG,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;YACrE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAC5D,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;QAChD,CAAC;QACD,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QACvF,IAAI,IAAI,GAAG,IAAI,KAAK,IAAI;YAAE,MAAM,IAAI,UAAU,CAAC,qCAAqC,CAAC,CAAC;IACxF,CAAC;IACD,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,CAAC;AACnE,CAAC"}