@developmentseed/deck.gl-raster 0.6.0-alpha.1 → 0.6.0-beta.1

package/dist/raster-tile-layer/raster-tile-layer.js

@@ -0,0 +1,175 @@
+import { CompositeLayer } from "@deck.gl/core";
+import { TileLayer } from "@deck.gl/geo-layers";
+import { renderDebugTileOutline } from "../layer-utils.js";
+import { RasterLayer } from "../raster-layer.js";
+import { RasterTileset2D } from "../raster-tileset/index.js";
+import { TILE_SIZE, WEB_MERCATOR_TO_WORLD_SCALE } from "./constants.js";
+const defaultProps = {
+    ...TileLayer.defaultProps,
+    maxError: 0.125,
+    debug: false,
+    debugOpacity: 0.5,
+};
+/**
+ * Base layer that renders a tiled raster source driven by a generic
+ * {@link TilesetDescriptor}.
+ *
+ * Usable directly (provide `tilesetDescriptor`, `getTileData`, and `renderTile`
+ * as props) or as a base class (override the protected `_tilesetDescriptor`,
+ * `_getTileDataCallback`, `_renderTileCallback` accessors to source them from
+ * state).
+ *
+ * The generic `ExtraProps` parameter lets a subclass redeclare any of the
+ * overridable fields with a domain-specific signature (e.g. `COGLayer`'s
+ * `getTileData(image, options)`).
+ */
+export class RasterTileLayer extends CompositeLayer {
+    static layerName = "RasterTileLayer";
+    static defaultProps = defaultProps;
+    /**
+     * The currently effective {@link TilesetDescriptor}.
+     *
+     * Subclasses override this to return a descriptor built from their own
+     * async-parsed state. Returns `undefined` while the source is still
+     * loading; `renderLayers()` returns `null` in that case.
+     *
+     * The inline cast to `RasterTileLayerProps<DataT>` is required because
+     * `tilesetDescriptor` is declared on `ExtraProps`, not on the base's
+     * `RasterTileLayerBaseProps`. For direct use the default `ExtraProps`
+     * brings it in; for subclass use this method is overridden and the cast
+     * is never reached.
+     */
+    _tilesetDescriptor() {
+        return this.props
+            .tilesetDescriptor;
+    }
+    /**
+     * The currently effective tile-fetch callback.
+     *
+     * Subclasses override this to adapt their user-facing `getTileData`
+     * signature into the base's `(tile, options) => Promise<DataT>` shape.
+     * Returns `undefined` when the callback is not yet available.
+     */
+    _getTileDataCallback() {
+        return this.props.getTileData;
+    }
+    /**
+     * The currently effective per-tile render callback.
+     *
+     * Subclasses override this to thread their user-facing `renderTile` and
+     * any inferred default. Returns `undefined` when no callback is available.
+     */
+    _renderTileCallback() {
+        return this.props.renderTile;
+    }
+    renderLayers() {
+        const descriptor = this._tilesetDescriptor();
+        const getTileData = this._getTileDataCallback();
+        const renderTile = this._renderTileCallback();
+        if (!descriptor || !getTileData || !renderTile) {
+            return null;
+        }
+        return this._renderTileLayer(descriptor, getTileData, renderTile);
+    }
+    _renderTileLayer(descriptor, getTileData, renderTile) {
+        class TilesetFactory extends RasterTileset2D {
+            constructor(opts) {
+                super(opts, descriptor);
+            }
+        }
+        const { tileSize, zoomOffset, maxZoom, minZoom, extent, debounceTime, maxCacheSize, maxCacheByteSize, maxRequests, refinementStrategy, updateTriggers, } = this.props;
+        return new TileLayer({
+            id: `raster-tile-layer-${this.id}`,
+            TilesetClass: TilesetFactory,
+            getTileData: (tile) => this._wrapGetTileData(tile, getTileData),
+            renderSubLayers: (props) => this._renderSubLayers(props, descriptor, renderTile),
+            updateTriggers: {
+                renderSubLayers: updateTriggers?.renderTile,
+            },
+            tileSize,
+            zoomOffset,
+            maxZoom,
+            minZoom,
+            extent,
+            debounceTime,
+            maxCacheSize,
+            maxCacheByteSize,
+            maxRequests,
+            refinementStrategy,
+        });
+    }
+    async _wrapGetTileData(tile, getTileData) {
+        const { signal: tileSignal } = tile;
+        const userSignal = this.props.signal;
+        const signal = userSignal && tileSignal
+            ? AbortSignal.any([userSignal, tileSignal])
+            : (userSignal ?? tileSignal);
+        const options = {
+            device: this.context.device,
+            signal,
+        };
+        return getTileData(tile, options);
+    }
+    _renderSubLayers(props, descriptor, renderTile) {
+        const { maxError, debug, debugOpacity } = this.props;
+        const tile = props.tile;
+        const layers = [];
+        if (debug) {
+            layers.push(...renderDebugTileOutline(`${this.id}-${tile.id}-bounds`, tile, descriptor.projectTo4326));
+        }
+        if (!props.data) {
+            return layers;
+        }
+        const { x, y, z } = tile.index;
+        const level = descriptor.levels[z];
+        if (!level) {
+            return layers;
+        }
+        const { forwardTransform, inverseTransform } = level.tileTransform(x, y);
+        const { image, renderPipeline } = renderTile(props.data);
+        const { width, height } = props.data;
+        const isGlobe = this.context.viewport.resolution !== undefined;
+        const reprojectionFns = isGlobe
+            ? {
+                forwardTransform,
+                inverseTransform,
+                forwardReproject: descriptor.projectTo4326,
+                inverseReproject: descriptor.projectFrom4326,
+            }
+            : {
+                forwardTransform,
+                inverseTransform,
+                forwardReproject: descriptor.projectTo3857,
+                inverseReproject: descriptor.projectFrom3857,
+            };
+        const deckProjectionProps = isGlobe
+            ? {}
+            : {
+                coordinateSystem: "cartesian",
+                coordinateOrigin: [TILE_SIZE / 2, TILE_SIZE / 2, 0],
+                // biome-ignore format: array
+                modelMatrix: [
+                    WEB_MERCATOR_TO_WORLD_SCALE, 0, 0, 0,
+                    0, WEB_MERCATOR_TO_WORLD_SCALE, 0, 0,
+                    0, 0, 1, 0,
+                    0, 0, 0, 1,
+                ],
+            };
+        const rasterLayer = new RasterLayer(this.getSubLayerProps({
+            id: `${props.id}-raster`,
+            width,
+            height,
+            // Passing `image: undefined` explicitly would trip isAsyncPropLoading
+            // and cause a transient black flash (see issue #376).
+            ...(image !== undefined && { image }),
+            renderPipeline,
+            maxError,
+            reprojectionFns,
+            debug,
+            debugOpacity,
+            ...deckProjectionProps,
+        }));
+        return [rasterLayer, ...layers];
+    }
+}
+//# sourceMappingURL=raster-tile-layer.js.map

package/dist/raster-tile-layer/raster-tile-layer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"raster-tile-layer.js","sourceRoot":"","sources":["../../src/raster-tile-layer/raster-tile-layer.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAO/C,OAAO,EAAE,SAAS,EAAE,MAAM,qBAAqB,CAAC;AAGhD,OAAO,EAAE,sBAAsB,EAAE,MAAM,mBAAmB,CAAC;AAE3D,OAAO,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC;AAEjD,OAAO,EAAE,eAAe,EAAE,MAAM,4BAA4B,CAAC;AAE7D,OAAO,EAAE,SAAS,EAAE,2BAA2B,EAAE,MAAM,gBAAgB,CAAC;AA8GxE,MAAM,YAAY,GAAuC;IACvD,GAAG,SAAS,CAAC,YAAY;IACzB,QAAQ,EAAE,KAAK;IACf,KAAK,EAAE,KAAK;IACZ,YAAY,EAAE,GAAG;CAClB,CAAC;AAyBF;;;;;;;;;;;;GAYG;AACH,MAAM,OAAO,eAGX,SAAQ,cAA4D;IACpE,MAAM,CAAU,SAAS,GAAG,iBAAiB,CAAC;IAC9C,MAAM,CAAU,YAAY,GAAG,YAAY,CAAC;IAE5C;;;;;;;;;;;;OAYG;IACO,kBAAkB;QAC1B,OAAQ,IAAI,CAAC,KAAgD;aAC1D,iBAAiB,CAAC;IACvB,CAAC;IAED;;;;;;OAMG;IACO,oBAAoB;QAC5B,OAAQ,IAAI,CAAC,KAAgD,CAAC,WAAW,CAAC;IAC5E,CAAC;IAED;;;;;OAKG;IACO,mBAAmB;QAC3B,OAAQ,IAAI,CAAC,KAAgD,CAAC,UAAU,CAAC;IAC3E,CAAC;IAEQ,YAAY;QACnB,MAAM,UAAU,GAAG,IAAI,CAAC,kBAAkB,EAAE,CAAC;QAC7C,MAAM,WAAW,GAAG,IAAI,CAAC,oBAAoB,EAAE,CAAC;QAChD,MAAM,UAAU,GAAG,IAAI,CAAC,mBAAmB,EAAE,CAAC;QAE9C,IAAI,CAAC,UAAU,IAAI,CAAC,WAAW,IAAI,CAAC,UAAU,EAAE,CAAC;YAC/C,OAAO,IAAI,CAAC;QACd,CAAC;QAED,OAAO,IAAI,CAAC,gBAAgB,CAAC,UAAU,EAAE,WAAW,EAAE,UAAU,CAAC,CAAC;IACpE,CAAC;IAEO,gBAAgB,CACtB,UAA6B,EAC7B,WAAoE,EACpE,UAAkE;QAElE,MAAM,cAAe,SAAQ,eAAe;YAC1C,YAAY,IAAoB;gBAC9B,KAAK,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC;YAC1B,CAAC;SACF;QAED,MAAM,EACJ,QAAQ,EACR,UAAU,EACV,OAAO,EACP,OAAO,EACP,MAAM,EACN,YAAY,EACZ,YAAY,EACZ,gBAAgB,EAChB,WAAW,EACX,kBAAkB,EAClB,cAAc,GACf,GAAG,IAAI,CAAC,KAAK,CAAC;QAEf,OAAO,IAAI,SAAS,CAAQ;YAC1B,EAAE,EAAE,qBAAqB,IAAI,CAAC,EAAE,EAAE;YAClC,YAAY,EAAE,cAAc;YAC5B,WAAW,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,WAAW,CAAC;YAC/D,eAAe,EAAE,CAAC,KAAK,EAAE,EAAE,CACzB,IAAI,CAAC,gBAAgB,CACnB,KAKC,EACD,UAAU,EACV,UAAU,CACX;YACH,cAAc,EAAE;gBACd,eAAe,EAAE,cAAc,EAAE,UAAU;aAC5C;YACD,QAAQ;YACR,UAAU;YACV,OAAO;YACP,OAAO;YACP,MAAM;YACN,YAAY;YACZ,YAAY;YACZ,gBAAgB;YAChB,WAAW;YACX,kBAAkB;SACnB,CAAC,CAAC;IACL,CAAC;IAEO,KAAK,CAAC,gBAAgB,CAC5B,IAAmB,EACnB,WAAoE;QAEpE,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,GAAG,IAAI,CAAC;QACpC,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC;QACrC,MAAM,MAAM,GACV,UAAU,IAAI,UAAU;YACtB,CAAC,CAAC,WAAW,CAAC,GAAG,CAAC,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;YAC3C,CAAC,CAAC,CAAC,UAAU,IAAI,UAAU,CAAC,CAAC;QACjC,MAAM,OAAO,GAAuB;YAClC,MAAM,EAAE,IAAI,CAAC,OAAO,CAAC,MAAM;YAC3B,MAAM;SACP,CAAC;QACF,OAAO,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IACpC,CAAC;IAEO,gBAAgB,CACtB,KAKC,EACD,UAA6B,EAC7B,UAAkE;QAElE,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,YAAY,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QACrD,MAAM,IAAI,GAAG,KAAK,CAAC,IAA0C,CAAC;QAE9D,MAAM,MAAM,GAAY,EAAE,CAAC;QAC3B,IAAI,KAAK,EAAE,CAAC;YACV,MAAM,CAAC,IAAI,CACT,GAAG,sBAAsB,CACvB,GAAG,IAAI,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,SAAS,EAC9B,IAAI,EACJ,UAAU,CAAC,aAAa,CACzB,CACF,CAAC;QACJ,CAAC;QACD,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;YAChB,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QAC/B,MAAM,KAAK,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,KAAK,EAAE,CAAC;YACX,OAAO,MAAM,CAAC;QAChB,CAAC;QACD,MAAM,EAAE,gBAAgB,EAAE,gBAAgB,EAAE,GAAG,KAAK,CAAC,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACzE,MAAM,EAAE,KAAK,EAAE,cAAc,EAAE,GAAG,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACzD,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC;QAErC,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,UAAU,KAAK,SAAS,CAAC;QAC/D,MAAM,eAAe,GAAoB,OAAO;YAC9C,CAAC,CAAC;gBACE,gBAAgB;gBAChB,gBAAgB;gBAChB,gBAAgB,EAAE,UAAU,CAAC,aAAa;gBAC1C,gBAAgB,EAAE,UAAU,CAAC,eAAe;aAC7C;YACH,CAAC,CAAC;gBACE,gBAAgB;gBAChB,gBAAgB;gBAChB,gBAAgB,EAAE,UAAU,CAAC,aAAa;gBAC1C,gBAAgB,EAAE,UAAU,CAAC,eAAe;aAC7C,CAAC;QACN,MAAM,mBAAmB,GAAwB,OAAO;YACtD,CAAC,CAAC,EAAE;YACJ,CAAC,CAAC;gBACE,gBAAgB,EAAE,WAAW;gBAC7B,gBAAgB,EAAE,CAAC,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,CAAC,CAAC;gBACnD,6BAA6B;gBAC7B,WAAW,EAAE;oBACX,2BAA2B,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;oBACpC,CAAC,EAAE,2BAA2B,EAAE,CAAC,EAAE,CAAC;oBACpC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;oBACV,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;iBACX;aACF,CAAC;QAEN,MAAM,WAAW,GAAG,IAAI,WAAW,CACjC,IAAI,CAAC,gBAAgB,CAAC;YACpB,EAAE,EAAE,GAAG,KAAK,CAAC,EAAE,SAAS;YACxB,KAAK;YACL,MAAM;YACN,sEAAsE;YACtE,sDAAsD;YACtD,GAAG,CAAC,KAAK,KAAK,SAAS,IAAI,EAAE,KAAK,EAAE,CAAC;YACrC,cAAc;YACd,QAAQ;YACR,eAAe;YACf,KAAK;YACL,YAAY;YACZ,GAAG,mBAAmB;SACvB,CAAC,CACH,CAAC;QACF,OAAO,CAAC,WAAW,EAAE,GAAG,MAAM,CAAC,CAAC;IAClC,CAAC"}

package/dist/raster-tileset/affine-tileset-level.d.ts

@@ -0,0 +1,56 @@
+import type { Affine } from "@developmentseed/affine";
+import type { TilesetLevel } from "./tileset-interface.js";
+import type { Bounds, Corners } from "./types.js";
+/**
+ * Constructor options for {@link AffineTilesetLevel}.
+ */
+export interface AffineTilesetLevelOptions {
+    /** Pixel → CRS affine transform for this resolution level. */
+    affine: Affine;
+    /** Full level width, in pixels. */
+    arrayWidth: number;
+    /** Full level height, in pixels. */
+    arrayHeight: number;
+    /** Tile width, in pixels. */
+    tileWidth: number;
+    /** Tile height, in pixels. */
+    tileHeight: number;
+    /** Meters per CRS unit (1 for metric CRSes, ≈111000 for WGS84). */
+    mpu: number;
+}
+/**
+ * A {@link TilesetLevel} described by a single affine transform plus tile and
+ * array sizes.
+ *
+ * This handles axis-aligned, rotated, skewed, and non-square-pixel grids
+ * uniformly. Sources that fit this shape (tiled GeoTIFF overviews, GeoZarr
+ * multiscales) can construct one of these per resolution level instead of
+ * implementing {@link TilesetLevel} manually.
+ */
+export declare class AffineTilesetLevel implements TilesetLevel {
+    readonly tileWidth: number;
+    readonly tileHeight: number;
+    readonly matrixWidth: number;
+    readonly matrixHeight: number;
+    readonly metersPerPixel: number;
+    /**
+     * Source-CRS bounding box of the level's array `[minX, minY, maxX, maxY]`.
+     * Computed from the affine applied to the four array corners.
+     */
+    readonly projectedBounds: Bounds;
+    private readonly _affine;
+    private readonly _invAffine;
+    constructor(options: AffineTilesetLevelOptions);
+    projectedTileCorners(col: number, row: number): Corners;
+    tileTransform(col: number, row: number): {
+        forwardTransform: (x: number, y: number) => [number, number];
+        inverseTransform: (x: number, y: number) => [number, number];
+    };
+    crsBoundsToTileRange(projectedMinX: number, projectedMinY: number, projectedMaxX: number, projectedMaxY: number): {
+        minCol: number;
+        maxCol: number;
+        minRow: number;
+        maxRow: number;
+    };
+}
+//# sourceMappingURL=affine-tileset-level.d.ts.map

package/dist/raster-tileset/affine-tileset-level.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine-tileset-level.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/affine-tileset-level.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,yBAAyB,CAAC;AAEtD,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,YAAY,CAAC;AAElD;;GAEG;AACH,MAAM,WAAW,yBAAyB;IACxC,8DAA8D;IAC9D,MAAM,EAAE,MAAM,CAAC;IACf,mCAAmC;IACnC,UAAU,EAAE,MAAM,CAAC;IACnB,oCAAoC;IACpC,WAAW,EAAE,MAAM,CAAC;IACpB,6BAA6B;IAC7B,SAAS,EAAE,MAAM,CAAC;IAClB,8BAA8B;IAC9B,UAAU,EAAE,MAAM,CAAC;IACnB,mEAAmE;IACnE,GAAG,EAAE,MAAM,CAAC;CACb;AAED;;;;;;;;GAQG;AACH,qBAAa,kBAAmB,YAAW,YAAY;IACrD,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;IAC5B,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,QAAQ,CAAC,cAAc,EAAE,MAAM,CAAC;IAChC;;;OAGG;IACH,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IAEjC,OAAO,CAAC,QAAQ,CAAC,OAAO,CAAS;IACjC,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAS;gBAExB,OAAO,EAAE,yBAAyB;IA6B9C,oBAAoB,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,OAAO;IAavD,aAAa,CACX,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM,GACV;QACD,gBAAgB,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAC7D,gBAAgB,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;KAC9D;IAaD,oBAAoB,CAClB,aAAa,EAAE,MAAM,EACrB,aAAa,EAAE,MAAM,EACrB,aAAa,EAAE,MAAM,EACrB,aAAa,EAAE,MAAM,GACpB;QAAE,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAA;KAAE;CAmCtE"}

package/dist/raster-tileset/affine-tileset-level.js

@@ -0,0 +1,101 @@
+import * as affine from "@developmentseed/affine";
+/**
+ * A {@link TilesetLevel} described by a single affine transform plus tile and
+ * array sizes.
+ *
+ * This handles axis-aligned, rotated, skewed, and non-square-pixel grids
+ * uniformly. Sources that fit this shape (tiled GeoTIFF overviews, GeoZarr
+ * multiscales) can construct one of these per resolution level instead of
+ * implementing {@link TilesetLevel} manually.
+ */
+export class AffineTilesetLevel {
+    tileWidth;
+    tileHeight;
+    matrixWidth;
+    matrixHeight;
+    metersPerPixel;
+    /**
+     * Source-CRS bounding box of the level's array `[minX, minY, maxX, maxY]`.
+     * Computed from the affine applied to the four array corners.
+     */
+    projectedBounds;
+    _affine;
+    _invAffine;
+    constructor(options) {
+        this._affine = options.affine;
+        this._invAffine = affine.invert(options.affine);
+        this.tileWidth = options.tileWidth;
+        this.tileHeight = options.tileHeight;
+        this.matrixWidth = Math.ceil(options.arrayWidth / options.tileWidth);
+        this.matrixHeight = Math.ceil(options.arrayHeight / options.tileHeight);
+        // Geometric mean of the two pixel-edge scales handles non-square pixels.
+        const a = affine.a(options.affine);
+        const e = affine.e(options.affine);
+        this.metersPerPixel = Math.sqrt(Math.abs(a * e)) * options.mpu;
+        const corners = [
+            affine.apply(options.affine, 0, 0),
+            affine.apply(options.affine, options.arrayWidth, 0),
+            affine.apply(options.affine, 0, options.arrayHeight),
+            affine.apply(options.affine, options.arrayWidth, options.arrayHeight),
+        ];
+        const xs = corners.map(([x]) => x);
+        const ys = corners.map(([, y]) => y);
+        this.projectedBounds = [
+            Math.min(...xs),
+            Math.min(...ys),
+            Math.max(...xs),
+            Math.max(...ys),
+        ];
+    }
+    projectedTileCorners(col, row) {
+        const tw = this.tileWidth;
+        const th = this.tileHeight;
+        const af = this._affine;
+        return {
+            topLeft: affine.apply(af, col * tw, row * th),
+            topRight: affine.apply(af, (col + 1) * tw, row * th),
+            bottomLeft: affine.apply(af, col * tw, (row + 1) * th),
+            bottomRight: affine.apply(af, (col + 1) * tw, (row + 1) * th),
+        };
+    }
+    tileTransform(col, row) {
+        const tileOffset = affine.translation(col * this.tileWidth, row * this.tileHeight);
+        const tileAffine = affine.compose(this._affine, tileOffset);
+        const invTileAffine = affine.invert(tileAffine);
+        return {
+            forwardTransform: (x, y) => affine.apply(tileAffine, x, y),
+            inverseTransform: (x, y) => affine.apply(invTileAffine, x, y),
+        };
+    }
+    crsBoundsToTileRange(projectedMinX, projectedMinY, projectedMaxX, projectedMaxY) {
+        // Map all four CRS corners through the inverse affine to get pixel coords,
+        // then take the bbox in pixel space (handles rotation correctly).
+        const inv = this._invAffine;
+        const pixelCorners = [
+            affine.apply(inv, projectedMinX, projectedMinY),
+            affine.apply(inv, projectedMaxX, projectedMinY),
+            affine.apply(inv, projectedMinX, projectedMaxY),
+            affine.apply(inv, projectedMaxX, projectedMaxY),
+        ];
+        const xs = pixelCorners.map(([px]) => px);
+        const ys = pixelCorners.map(([, py]) => py);
+        const pixMinX = Math.min(...xs);
+        const pixMaxX = Math.max(...xs);
+        const pixMinY = Math.min(...ys);
+        const pixMaxY = Math.max(...ys);
+        const tw = this.tileWidth;
+        const th = this.tileHeight;
+        const maxColIdx = this.matrixWidth - 1;
+        const maxRowIdx = this.matrixHeight - 1;
+        // Asymmetric clamping: only clamp minCol/minRow up from below and
+        // maxCol/maxRow down from above. If the bbox lies entirely outside the
+        // array, this produces an empty range (min > max) so the consumer's
+        // `for (i = min; i <= max; i++)` loop does nothing.
+        const minCol = Math.max(0, Math.floor(pixMinX / tw));
+        const maxCol = Math.min(maxColIdx, Math.floor(pixMaxX / tw));
+        const minRow = Math.max(0, Math.floor(pixMinY / th));
+        const maxRow = Math.min(maxRowIdx, Math.floor(pixMaxY / th));
+        return { minCol, maxCol, minRow, maxRow };
+    }
+}
+//# sourceMappingURL=affine-tileset-level.js.map

package/dist/raster-tileset/affine-tileset-level.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine-tileset-level.js","sourceRoot":"","sources":["../../src/raster-tileset/affine-tileset-level.ts"],"names":[],"mappings":"AACA,OAAO,KAAK,MAAM,MAAM,yBAAyB,CAAC;AAsBlD;;;;;;;;GAQG;AACH,MAAM,OAAO,kBAAkB;IACpB,SAAS,CAAS;IAClB,UAAU,CAAS;IACnB,WAAW,CAAS;IACpB,YAAY,CAAS;IACrB,cAAc,CAAS;IAChC;;;OAGG;IACM,eAAe,CAAS;IAEhB,OAAO,CAAS;IAChB,UAAU,CAAS;IAEpC,YAAY,OAAkC;QAC5C,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC,MAAM,CAAC;QAC9B,IAAI,CAAC,UAAU,GAAG,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAChD,IAAI,CAAC,SAAS,GAAG,OAAO,CAAC,SAAS,CAAC;QACnC,IAAI,CAAC,UAAU,GAAG,OAAO,CAAC,UAAU,CAAC;QACrC,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,UAAU,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC;QACrE,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,WAAW,GAAG,OAAO,CAAC,UAAU,CAAC,CAAC;QAExE,yEAAyE;QACzE,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QACnC,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC;QAE/D,MAAM,OAAO,GAAG;YACd,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,CAAC;YAClC,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,EAAE,OAAO,CAAC,UAAU,EAAE,CAAC,CAAC;YACnD,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,OAAO,CAAC,WAAW,CAAC;YACpD,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,EAAE,OAAO,CAAC,UAAU,EAAE,OAAO,CAAC,WAAW,CAAC;SACtE,CAAC;QACF,MAAM,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACrC,IAAI,CAAC,eAAe,GAAG;YACrB,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YACf,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YACf,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;YACf,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;SAChB,CAAC;IACJ,CAAC;IAED,oBAAoB,CAAC,GAAW,EAAE,GAAW;QAC3C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC;QAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC;QAExB,OAAO;YACL,OAAO,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC;YAC7C,QAAQ,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,CAAC;YACpD,UAAU,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YACtD,WAAW,EAAE,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;SAC9D,CAAC;IACJ,CAAC;IAED,aAAa,CACX,GAAW,EACX,GAAW;QAKX,MAAM,UAAU,GAAG,MAAM,CAAC,WAAW,CACnC,GAAG,GAAG,IAAI,CAAC,SAAS,EACpB,GAAG,GAAG,IAAI,CAAC,UAAU,CACtB,CAAC;QACF,MAAM,UAAU,GAAG,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,OAAO,EAAE,UAAU,CAAC,CAAC;QAC5D,MAAM,aAAa,GAAG,MAAM,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC;QAChD,OAAO;YACL,gBAAgB,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,EAAE,CAAC,CAAC;YAC1D,gBAAgB,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,CAAC;SAC9D,CAAC;IACJ,CAAC;IAED,oBAAoB,CAClB,aAAqB,EACrB,aAAqB,EACrB,aAAqB,EACrB,aAAqB;QAErB,2EAA2E;QAC3E,kEAAkE;QAClE,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC;QAC5B,MAAM,YAAY,GAAG;YACnB,MAAM,CAAC,KAAK,CAAC,GAAG,EAAE,aAAa,EAAE,aAAa,CAAC;YAC/C,MAAM,CAAC,KAAK,CAAC,GAAG,EAAE,aAAa,EAAE,aAAa,CAAC;YAC/C,MAAM,CAAC,KAAK,CAAC,GAAG,EAAE,aAAa,EAAE,aAAa,CAAC;YAC/C,MAAM,CAAC,KAAK,CAAC,GAAG,EAAE,aAAa,EAAE,aAAa,CAAC;SAChD,CAAC;QAEF,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;QAE5C,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAChC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAEhC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC;QAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC;QAC3B,MAAM,SAAS,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,CAAC;QACvC,MAAM,SAAS,GAAG,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC;QAExC,kEAAkE;QAClE,uEAAuE;QACvE,oEAAoE;QACpE,oDAAoD;QACpD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC;QACrD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC;QAC7D,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC;QACrD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC;QAE7D,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC;IAC5C,CAAC;CACF"}

package/dist/raster-tileset/affine-tileset.d.ts

@@ -0,0 +1,34 @@
+import type { AffineTilesetLevel } from "./affine-tileset-level.js";
+import type { TilesetDescriptor } from "./tileset-interface.js";
+import type { Bounds, ProjectionFunction } from "./types.js";
+/**
+ * Constructor options for {@link AffineTileset}.
+ */
+export interface AffineTilesetOptions {
+    /** Resolution levels, ordered coarsest first. */
+    levels: AffineTilesetLevel[];
+    /** Forward projection function from source CRS to EPSG:3857. */
+    projectTo3857: ProjectionFunction;
+    /** Inverse projection function from EPSG:3857 to source CRS. */
+    projectFrom3857: ProjectionFunction;
+    /** Forward projection function from source CRS to EPSG:4326. */
+    projectTo4326: ProjectionFunction;
+    /** Inverse projection function from EPSG:4326 to source CRS. */
+    projectFrom4326: ProjectionFunction;
+}
+/**
+ * A {@link TilesetDescriptor} backed by per-level affine transforms.
+ *
+ * Derives `projectedBounds` from the coarsest level's array. Everything else
+ * is passed through from the constructor options.
+ */
+export declare class AffineTileset implements TilesetDescriptor {
+    readonly levels: AffineTilesetLevel[];
+    readonly projectTo3857: ProjectionFunction;
+    readonly projectFrom3857: ProjectionFunction;
+    readonly projectTo4326: ProjectionFunction;
+    readonly projectFrom4326: ProjectionFunction;
+    readonly projectedBounds: Bounds;
+    constructor(options: AffineTilesetOptions);
+}
+//# sourceMappingURL=affine-tileset.d.ts.map

package/dist/raster-tileset/affine-tileset.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine-tileset.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/affine-tileset.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AACpE,OAAO,KAAK,EAAE,iBAAiB,EAAE,MAAM,wBAAwB,CAAC;AAChE,OAAO,KAAK,EAAE,MAAM,EAAE,kBAAkB,EAAE,MAAM,YAAY,CAAC;AAE7D;;GAEG;AACH,MAAM,WAAW,oBAAoB;IACnC,iDAAiD;IACjD,MAAM,EAAE,kBAAkB,EAAE,CAAC;IAC7B,gEAAgE;IAChE,aAAa,EAAE,kBAAkB,CAAC;IAClC,gEAAgE;IAChE,eAAe,EAAE,kBAAkB,CAAC;IACpC,gEAAgE;IAChE,aAAa,EAAE,kBAAkB,CAAC;IAClC,gEAAgE;IAChE,eAAe,EAAE,kBAAkB,CAAC;CACrC;AAED;;;;;GAKG;AACH,qBAAa,aAAc,YAAW,iBAAiB;IACrD,QAAQ,CAAC,MAAM,EAAE,kBAAkB,EAAE,CAAC;IACtC,QAAQ,CAAC,aAAa,EAAE,kBAAkB,CAAC;IAC3C,QAAQ,CAAC,eAAe,EAAE,kBAAkB,CAAC;IAC7C,QAAQ,CAAC,aAAa,EAAE,kBAAkB,CAAC;IAC3C,QAAQ,CAAC,eAAe,EAAE,kBAAkB,CAAC;IAC7C,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;gBAErB,OAAO,EAAE,oBAAoB;CAW1C"}

package/dist/raster-tileset/affine-tileset.js

@@ -0,0 +1,26 @@
+/**
+ * A {@link TilesetDescriptor} backed by per-level affine transforms.
+ *
+ * Derives `projectedBounds` from the coarsest level's array. Everything else
+ * is passed through from the constructor options.
+ */
+export class AffineTileset {
+    levels;
+    projectTo3857;
+    projectFrom3857;
+    projectTo4326;
+    projectFrom4326;
+    projectedBounds;
+    constructor(options) {
+        if (options.levels.length === 0) {
+            throw new Error("AffineTileset requires at least one level");
+        }
+        this.levels = options.levels;
+        this.projectTo3857 = options.projectTo3857;
+        this.projectFrom3857 = options.projectFrom3857;
+        this.projectTo4326 = options.projectTo4326;
+        this.projectFrom4326 = options.projectFrom4326;
+        this.projectedBounds = options.levels[0].projectedBounds;
+    }
+}
+//# sourceMappingURL=affine-tileset.js.map

package/dist/raster-tileset/affine-tileset.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"affine-tileset.js","sourceRoot":"","sources":["../../src/raster-tileset/affine-tileset.ts"],"names":[],"mappings":"AAoBA;;;;;GAKG;AACH,MAAM,OAAO,aAAa;IACf,MAAM,CAAuB;IAC7B,aAAa,CAAqB;IAClC,eAAe,CAAqB;IACpC,aAAa,CAAqB;IAClC,eAAe,CAAqB;IACpC,eAAe,CAAS;IAEjC,YAAY,OAA6B;QACvC,IAAI,OAAO,CAAC,MAAM,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAChC,MAAM,IAAI,KAAK,CAAC,2CAA2C,CAAC,CAAC;QAC/D,CAAC;QACD,IAAI,CAAC,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;QAC7B,IAAI,CAAC,aAAa,GAAG,OAAO,CAAC,aAAa,CAAC;QAC3C,IAAI,CAAC,eAAe,GAAG,OAAO,CAAC,eAAe,CAAC;QAC/C,IAAI,CAAC,aAAa,GAAG,OAAO,CAAC,aAAa,CAAC;QAC3C,IAAI,CAAC,eAAe,GAAG,OAAO,CAAC,eAAe,CAAC;QAC/C,IAAI,CAAC,eAAe,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAE,CAAC,eAAe,CAAC;IAC5D,CAAC;CACF"}

package/dist/raster-tileset/index.d.ts

@@ -1,3 +1,7 @@
+export type { AffineTilesetOptions } from "./affine-tileset.js";
+export { AffineTileset } from "./affine-tileset.js";
+export type { AffineTilesetLevelOptions } from "./affine-tileset-level.js";
+export { AffineTilesetLevel } from "./affine-tileset-level.js";
 export type { TileMetadata } from "./raster-tileset-2d.js";
 export { RasterTileset2D } from "./raster-tileset-2d.js";
 export { TileMatrixSetAdaptor } from "./tile-matrix-set.js";

package/dist/raster-tileset/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/index.ts"],"names":[],"mappings":"AAAA,YAAY,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AACzD,OAAO,EAAE,oBAAoB,EAAE,MAAM,sBAAsB,CAAC;AAC5D,YAAY,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC9E,YAAY,EACV,MAAM,EACN,YAAY,EACZ,OAAO,EACP,kBAAkB,GACnB,MAAM,YAAY,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/index.ts"],"names":[],"mappings":"AAAA,YAAY,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAChE,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AACpD,YAAY,EAAE,yBAAyB,EAAE,MAAM,2BAA2B,CAAC;AAC3E,OAAO,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AAC/D,YAAY,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AACzD,OAAO,EAAE,oBAAoB,EAAE,MAAM,sBAAsB,CAAC;AAC5D,YAAY,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC9E,YAAY,EACV,MAAM,EACN,YAAY,EACZ,OAAO,EACP,kBAAkB,GACnB,MAAM,YAAY,CAAC"}

package/dist/raster-tileset/index.js

@@ -1,3 +1,5 @@
+export { AffineTileset } from "./affine-tileset.js";
+export { AffineTilesetLevel } from "./affine-tileset-level.js";
 export { RasterTileset2D } from "./raster-tileset-2d.js";
 export { TileMatrixSetAdaptor } from "./tile-matrix-set.js";
 //# sourceMappingURL=index.js.map

package/dist/raster-tileset/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/raster-tileset/index.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AACzD,OAAO,EAAE,oBAAoB,EAAE,MAAM,sBAAsB,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/raster-tileset/index.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAEpD,OAAO,EAAE,kBAAkB,EAAE,MAAM,2BAA2B,CAAC;AAE/D,OAAO,EAAE,eAAe,EAAE,MAAM,wBAAwB,CAAC;AACzD,OAAO,EAAE,oBAAoB,EAAE,MAAM,sBAAsB,CAAC"}

package/dist/raster-tileset/raster-tile-traversal.d.ts

@@ -139,6 +139,28 @@ export declare class RasterTileNode {
      */
     private _getGenericBoundingVolume;
 }
+/**
+ * Build the list of root (z=0) `RasterTileNode`s for the traversal.
+ *
+ * Small root matrices (≤ {@link MAX_ROOT_TILES_NO_CULL}) are enumerated
+ * directly — traditional pyramids with a 1×1 or 4×5 root grid skip any
+ * projection work and keep bit-identical behavior to the pre-optimization
+ * traversal.
+ *
+ * Large root matrices are culled to the intersection of the dataset extent
+ * (`datasetWgs84Bounds`) and the viewport's WGS84 bounds, projected into
+ * the source CRS via `transformBounds` (which densifies the edges so a
+ * curving projection doesn't escape the 4-corner hull). If the viewport
+ * and dataset don't overlap, an empty array is returned and the rest of
+ * the traversal short-circuits.
+ *
+ * Exported for unit testing.
+ */
+export declare function createRootTiles(opts: {
+    descriptor: TilesetDescriptor;
+    viewport: Pick<Viewport, "getBounds">;
+    datasetWgs84Bounds: Bounds;
+}): RasterTileNode[];
 /**
  * Get tile indices visible in viewport.
  *

package/dist/raster-tileset/raster-tile-traversal.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"raster-tile-traversal.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/raster-tile-traversal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;GAiBG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAE9C,OAAO,KAAK,EAAE,mBAAmB,EAAE,MAAM,kBAAkB,CAAC;AAC5D,OAAO,EACL,aAAa,EAGd,MAAM,kBAAkB,CAAC;AAG1B,OAAO,KAAK,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC9E,OAAO,KAAK,EACV,MAAM,EAIN,SAAS,EACT,MAAM,EACP,MAAM,YAAY,CAAC;AAiEpB;;;;;;;;;;GAUG;AACH,qBAAa,cAAc;IACzB,yBAAyB;IACzB,CAAC,EAAE,MAAM,CAAC;IAEV,0BAA0B;IAC1B,CAAC,EAAE,MAAM,CAAC;IAEV,uDAAuD;IACvD,CAAC,EAAE,MAAM,CAAC;IAEV,OAAO,CAAC,UAAU,CAAoB;IAEtC;;;;;OAKG;IACH,OAAO,CAAC,YAAY,CAAC,CAAU;IAE/B;;;;OAIG;IACH,OAAO,CAAC,QAAQ,CAAC,CAAU;IAE3B,4CAA4C;IAC5C,OAAO,CAAC,SAAS,CAAC,CAA0B;IAE5C;;;;OAIG;IACH,OAAO,CAAC,eAAe,CAAC,CAItB;gBAGA,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,UAAU,EAAE,EAAE;QAAE,UAAU,EAAE,iBAAiB,CAAA;KAAE;IAQnD,kDAAkD;IAClD,IAAI,KAAK,IAAI,YAAY,CAExB;IAED;;;;;;;OAOG;IACH,IAAI,QAAQ,IAAI,cAAc,EAAE,GAAG,IAAI,CA+BtC;IAED;;;;;;;;;;;;;;;;;OAiBG;IACH,MAAM,CAAC,MAAM,EAAE;QACb,QAAQ,EAAE,QAAQ,CAAC;QAEnB,OAAO,EAAE,CAAC,CAAC,GAAG,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,CAAC,GAAG,IAAI,CAAC;QAE9C,aAAa,EAAE,aAAa,CAAC;QAE7B,eAAe,EAAE,MAAM,CAAC;QACxB,wCAAwC;QACxC,IAAI,EAAE,MAAM,CAAC;QACb,sCAAsC;QACtC,IAAI,CAAC,EAAE,MAAM,CAAC;QACd,wCAAwC;QACxC,MAAM,CAAC,EAAE,MAAM,CAAC;KACjB,GAAG,OAAO;IAiFX;;;;;;OAMG;IACH,WAAW,CAAC,MAAM,GAAE,cAAc,EAAO,GAAG,cAAc,EAAE;IAY5D;;;;;;OAMG;IACH,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,GAAG,OAAO;IAUhE;;;;;;OAMG;IACH,iBAAiB,CACf,MAAM,EAAE,MAAM,EACd,OAAO,EAAE,CAAC,CAAC,GAAG,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,CAAC,GAAG,IAAI,GAC5C;QAAE,cAAc,EAAE,mBAAmB,CAAC;QAAC,iBAAiB,EAAE,MAAM,CAAA;KAAE;IAgCrE;;;;OAIG;IACH,OAAO,CAAC,yBAAyB;CAiDlC;AA4GD;;;;;;;GAOG;AACH,wBAAgB,cAAc,CAC5B,UAAU,EAAE,iBAAiB,EAC7B,IAAI,EAAE;IACJ,QAAQ,EAAE,QAAQ,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,WAAW,EAAE,MAAM,CAAC;CACrB,GACA,SAAS,EAAE,CA4Fb"}
+{"version":3,"file":"raster-tile-traversal.d.ts","sourceRoot":"","sources":["../../src/raster-tileset/raster-tile-traversal.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;GAiBG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,eAAe,CAAC;AAG9C,OAAO,KAAK,EAAE,mBAAmB,EAAE,MAAM,kBAAkB,CAAC;AAC5D,OAAO,EACL,aAAa,EAGd,MAAM,kBAAkB,CAAC;AAG1B,OAAO,KAAK,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,wBAAwB,CAAC;AAC9E,OAAO,KAAK,EACV,MAAM,EAIN,SAAS,EACT,MAAM,EACP,MAAM,YAAY,CAAC;AAiEpB;;;;;;;;;;GAUG;AACH,qBAAa,cAAc;IACzB,yBAAyB;IACzB,CAAC,EAAE,MAAM,CAAC;IAEV,0BAA0B;IAC1B,CAAC,EAAE,MAAM,CAAC;IAEV,uDAAuD;IACvD,CAAC,EAAE,MAAM,CAAC;IAEV,OAAO,CAAC,UAAU,CAAoB;IAEtC;;;;;OAKG;IACH,OAAO,CAAC,YAAY,CAAC,CAAU;IAE/B;;;;OAIG;IACH,OAAO,CAAC,QAAQ,CAAC,CAAU;IAE3B,4CAA4C;IAC5C,OAAO,CAAC,SAAS,CAAC,CAA0B;IAE5C;;;;OAIG;IACH,OAAO,CAAC,eAAe,CAAC,CAItB;gBAGA,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,CAAC,EAAE,MAAM,EACT,EAAE,UAAU,EAAE,EAAE;QAAE,UAAU,EAAE,iBAAiB,CAAA;KAAE;IAQnD,kDAAkD;IAClD,IAAI,KAAK,IAAI,YAAY,CAExB;IAED;;;;;;;OAOG;IACH,IAAI,QAAQ,IAAI,cAAc,EAAE,GAAG,IAAI,CA+BtC;IAED;;;;;;;;;;;;;;;;;OAiBG;IACH,MAAM,CAAC,MAAM,EAAE;QACb,QAAQ,EAAE,QAAQ,CAAC;QAEnB,OAAO,EAAE,CAAC,CAAC,GAAG,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,CAAC,GAAG,IAAI,CAAC;QAE9C,aAAa,EAAE,aAAa,CAAC;QAE7B,eAAe,EAAE,MAAM,CAAC;QACxB,wCAAwC;QACxC,IAAI,EAAE,MAAM,CAAC;QACb,sCAAsC;QACtC,IAAI,CAAC,EAAE,MAAM,CAAC;QACd,wCAAwC;QACxC,MAAM,CAAC,EAAE,MAAM,CAAC;KACjB,GAAG,OAAO;IAiFX;;;;;;OAMG;IACH,WAAW,CAAC,MAAM,GAAE,cAAc,EAAO,GAAG,cAAc,EAAE;IAY5D;;;;;;OAMG;IACH,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,iBAAiB,EAAE,MAAM,GAAG,OAAO;IAUhE;;;;;;OAMG;IACH,iBAAiB,CACf,MAAM,EAAE,MAAM,EACd,OAAO,EAAE,CAAC,CAAC,GAAG,EAAE,MAAM,EAAE,KAAK,MAAM,EAAE,CAAC,GAAG,IAAI,GAC5C;QAAE,cAAc,EAAE,mBAAmB,CAAC;QAAC,iBAAiB,EAAE,MAAM,CAAA;KAAE;IAgCrE;;;;OAIG;IACH,OAAO,CAAC,yBAAyB;CAyDlC;AAsHD;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,eAAe,CAAC,IAAI,EAAE;IACpC,UAAU,EAAE,iBAAiB,CAAC;IAC9B,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,WAAW,CAAC,CAAC;IACtC,kBAAkB,EAAE,MAAM,CAAC;CAC5B,GAAG,cAAc,EAAE,CA4CnB;AAED;;;;;;;GAOG;AACH,wBAAgB,cAAc,CAC5B,UAAU,EAAE,iBAAiB,EAC7B,IAAI,EAAE;IACJ,QAAQ,EAAE,QAAQ,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;IACb,MAAM,EAAE,MAAM,GAAG,IAAI,CAAC;IACtB,WAAW,EAAE,MAAM,CAAC;CACrB,GACA,SAAS,EAAE,CAqFb"}

package/dist/raster-tileset/raster-tile-traversal.js

@@ -17,6 +17,7 @@
  * both OGC TileMatrixSet grids and Zarr multiscale pyramids.
  */
 import { _GlobeViewport, assert } from "@deck.gl/core";
+import { transformBounds } from "@developmentseed/proj";
 import { CullingVolume, makeOrientedBoundingBoxFromPoints, Plane, } from "@math.gl/culling";
 import { lngLatToWorld, worldToLngLat } from "@math.gl/web-mercator";
 /**
@@ -316,14 +317,18 @@ export class RasterTileNode {
         let maxX = Number.NEGATIVE_INFINITY;
         let maxY = Number.NEGATIVE_INFINITY;
         for (const [x, y] of commonSpacePositions) {
-            if (x < minX)
+            if (x < minX) {
                 minX = x;
-            if (y < minY)
+            }
+            if (y < minY) {
                 minY = y;
-            if (x > maxX)
+            }
+            if (x > maxX) {
                 maxX = x;
-            if (y > maxY)
+            }
+            if (y > maxY) {
                 maxY = y;
+            }
         }
         const commonSpaceBounds = [minX, minY, maxX, maxY];
         return {
@@ -404,6 +409,68 @@ function rescaleEPSG3857ToCommonSpace([x, y]) {
         (clampedY / EPSG_3857_CIRCUMFERENCE + 0.5) * TILE_SIZE,
     ];
 }
+/**
+ * Above this root-tile count, `createRootTiles` culls to the viewport
+ * before instantiation. Below it, every root tile is created and downstream
+ * frustum culling filters the unused ones. Typical OGC pyramids have 1–a
+ * few dozen tiles at z=0, so they stay on the unchanged path. Large
+ * single-level zarr descriptors (e.g. AEF mosaic: ~15000 × 7000 ≈ 100M root
+ * tiles) must take the culled path or instantiation hangs the page.
+ */
+const MAX_ROOT_TILES_NO_CULL = 100;
+/**
+ * Build the list of root (z=0) `RasterTileNode`s for the traversal.
+ *
+ * Small root matrices (≤ {@link MAX_ROOT_TILES_NO_CULL}) are enumerated
+ * directly — traditional pyramids with a 1×1 or 4×5 root grid skip any
+ * projection work and keep bit-identical behavior to the pre-optimization
+ * traversal.
+ *
+ * Large root matrices are culled to the intersection of the dataset extent
+ * (`datasetWgs84Bounds`) and the viewport's WGS84 bounds, projected into
+ * the source CRS via `transformBounds` (which densifies the edges so a
+ * curving projection doesn't escape the 4-corner hull). If the viewport
+ * and dataset don't overlap, an empty array is returned and the rest of
+ * the traversal short-circuits.
+ *
+ * Exported for unit testing.
+ */
+export function createRootTiles(opts) {
+    const { descriptor, viewport, datasetWgs84Bounds } = opts;
+    const rootLevel = descriptor.levels[0];
+    const roots = [];
+    const rootTileCount = rootLevel.matrixWidth * rootLevel.matrixHeight;
+    if (rootTileCount <= MAX_ROOT_TILES_NO_CULL) {
+        // Small root matrix → enumerate every tile; downstream frustum culling
+        // handles the small amount of waste.
+        for (let y = 0; y < rootLevel.matrixHeight; y++) {
+            for (let x = 0; x < rootLevel.matrixWidth; x++) {
+                roots.push(new RasterTileNode(x, y, 0, { descriptor }));
+            }
+        }
+        return roots;
+    }
+    // Large root matrix → intersect dataset extent with viewport, project
+    // to source CRS, use the root level's tile-range helper.
+    const vpBounds = viewport.getBounds();
+    const cullBounds = [
+        Math.max(datasetWgs84Bounds[0], vpBounds[0]),
+        Math.max(datasetWgs84Bounds[1], vpBounds[1]),
+        Math.min(datasetWgs84Bounds[2], vpBounds[2]),
+        Math.min(datasetWgs84Bounds[3], vpBounds[3]),
+    ];
+    if (cullBounds[0] > cullBounds[2] || cullBounds[1] > cullBounds[3]) {
+        return roots;
+    }
+    const [minX, minY, maxX, maxY] = transformBounds(descriptor.projectFrom4326, cullBounds[0], cullBounds[1], cullBounds[2], cullBounds[3]);
+    const rootRange = rootLevel.crsBoundsToTileRange(minX, minY, maxX, maxY);
+    for (let y = rootRange.minRow; y <= rootRange.maxRow; y++) {
+        for (let x = rootRange.minCol; x <= rootRange.maxCol; x++) {
+            roots.push(new RasterTileNode(x, y, 0, { descriptor }));
+        }
+    }
+    return roots;
+}
 /**
  * Get tile indices visible in viewport.
  *
@@ -460,17 +527,11 @@ export function getTileIndices(descriptor, opts) {
         topRight[0],
         topRight[1],
     ];
-    // Start from coarsest level
-    const rootLevel = descriptor.levels[0];
-    // Create root tiles at coarsest level.
-    // In contrary to OSM tiling, we might have more than one tile at the
-    // coarsest level (z=0).
-    const roots = [];
-    for (let y = 0; y < rootLevel.matrixHeight; y++) {
-        for (let x = 0; x < rootLevel.matrixWidth; x++) {
-            roots.push(new RasterTileNode(x, y, 0, { descriptor }));
-        }
-    }
+    const roots = createRootTiles({
+        descriptor,
+        viewport,
+        datasetWgs84Bounds: wgs84Bounds,
+    });
     // Traverse and update visibility
     const traversalParams = {
         viewport,