@basemaps/cli-raster 8.0.0

package/build/cogify/covering/cutline.js

@@ -0,0 +1,109 @@
+import { Epsg, EpsgCode, Projection } from '@basemaps/geo';
+import { fsa } from '@basemaps/shared';
+import { Area, featuresToMultiPolygon, intersection, toFeatureCollection, toFeatureMultiPolygon, } from '@linzjs/geojson';
+export async function loadCutline(path) {
+    const buf = await fsa.read(path);
+    if (path.pathname.endsWith('.geojson') || path.pathname.endsWith('.json')) {
+        const data = JSON.parse(buf.toString());
+        const projection = Epsg.parseCode(data.crs?.properties?.name ?? '') ?? EpsgCode.Wgs84;
+        const polygon = featuresToMultiPolygon(data.features, true).coordinates;
+        return { polygon, projection };
+    }
+    throw new Error('Unknown cutline type: ' + path.href);
+}
+export class CutlineOptimizer {
+    constructor(path, cutline, blend, tileMatrix) {
+        Object.defineProperty(this, "path", {
+            enumerable: true,
+            configurable: true,
+            writable: true,
+            value: void 0
+        });
+        Object.defineProperty(this, "cutline", {
+            enumerable: true,
+            configurable: true,
+            writable: true,
+            value: void 0
+        });
+        Object.defineProperty(this, "blend", {
+            enumerable: true,
+            configurable: true,
+            writable: true,
+            value: void 0
+        });
+        Object.defineProperty(this, "tileMatrix", {
+            enumerable: true,
+            configurable: true,
+            writable: true,
+            value: void 0
+        });
+        this.path = path;
+        this.cutline = cutline;
+        this.blend = blend;
+        this.tileMatrix = tileMatrix;
+    }
+    /** Cut a polygon to the cutline */
+    cut(poly) {
+        if (this.cutline == null)
+            return poly;
+        // TODO buffer the cutline by this.blend so blended tiles will be included
+        // although it is somewhat wasteful to make a cog just for a blended cut
+        return intersection(poly, this.cutline);
+    }
+    /** Optimize a cutline for a tile, if it doesn't intersect then ignore the cutline */
+    optimize(tile) {
+        if (this.cutline == null)
+            return null;
+        const tileBounds = this.tileMatrix.tileToSourceBounds(tile);
+        const tileScale = this.tileMatrix.pixelScale(tile.z);
+        // Expand the tile bounds to ensure optimized bounds are not clipping and blended.
+        const scaleAmount = (tileBounds.width + tileScale * (CutlineOptimizer.TilePixelPadding + this.blend) * 2) / tileBounds.width;
+        const scaledBounds = tileBounds.scaleFromCenter(scaleAmount);
+        const optimized = intersection(this.cutline, scaledBounds.toPolygon());
+        // Check how much of a change was the optimized cutline
+        const optimizedArea = Area.multiPolygon(optimized);
+        const scaledArea = scaledBounds.width * scaledBounds.height;
+        const areaPercent = 1 - optimizedArea / scaledArea;
+        // if area is bigger than 0.0001% then assume the cutline was applied
+        if (areaPercent < 1e-4)
+            return null;
+        const feature = toFeatureCollection([toFeatureMultiPolygon(optimized)]);
+        feature.crs = { type: 'name', properties: { name: this.tileMatrix.projection.toUrn() } };
+        return feature;
+    }
+    /**
+     * Convert cutline to geojson FeatureCollection
+     *
+     * @returns GeoJSON FeatureCollection if cutline exists, null otherwise
+     */
+    toGeoJson() {
+        if (this.cutline == null)
+            return null;
+        const feature = toFeatureCollection([toFeatureMultiPolygon(this.cutline)]);
+        feature.crs = { type: 'name', properties: { name: this.tileMatrix.projection.toUrn() } };
+        return feature;
+    }
+    static async load(path, blend, tileMatrix) {
+        if (path == null)
+            return new CutlineOptimizer(null, null, blend, tileMatrix);
+        const cut = await loadCutline(path);
+        if (cut.projection === tileMatrix.projection.code)
+            return new CutlineOptimizer(path, cut.polygon, blend, tileMatrix);
+        const projected = Projection.get(cut.projection).projectMultipolygon(cut.polygon, Projection.get(tileMatrix.projection));
+        return new CutlineOptimizer(path, projected, blend, tileMatrix);
+    }
+    /** Create a cutline optimizer from a STAC cutline link */
+    static async loadFromLink(l, tileMatrix) {
+        if (l == null)
+            return new CutlineOptimizer(null, null, 0, tileMatrix);
+        return this.load(new URL(l.href), l.blend, tileMatrix);
+    }
+}
+/** Pad tiles by 200m when optimizing cutlines so blending will not be affected */
+Object.defineProperty(CutlineOptimizer, "TilePixelPadding", {
+    enumerable: true,
+    configurable: true,
+    writable: true,
+    value: 200
+});
+//# sourceMappingURL=cutline.js.map

package/build/cogify/covering/cutline.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cutline.js","sourceRoot":"","sources":["../../../src/cogify/covering/cutline.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,UAAU,EAAuB,MAAM,eAAe,CAAC;AAChF,OAAO,EAAE,GAAG,EAAE,MAAM,kBAAkB,CAAC;AACvC,OAAO,EACL,IAAI,EAEJ,sBAAsB,EACtB,YAAY,EAEZ,mBAAmB,EACnB,qBAAqB,GACtB,MAAM,iBAAiB,CAAC;AAIzB,MAAM,CAAC,KAAK,UAAU,WAAW,CAAC,IAAS;IACzC,MAAM,GAAG,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAEjC,IAAI,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,UAAU,CAAC,IAAI,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,OAAO,CAAC,EAAE,CAAC;QAC1E,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,QAAQ,EAAE,CAA6B,CAAC;QACpE,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,UAAU,EAAE,IAAI,IAAI,EAAE,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC;QACtF,MAAM,OAAO,GAAG,sBAAsB,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC,WAA2B,CAAC;QACxF,OAAO,EAAE,OAAO,EAAE,UAAU,EAAE,CAAC;IACjC,CAAC;IAED,MAAM,IAAI,KAAK,CAAC,wBAAwB,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;AACxD,CAAC;AAED,MAAM,OAAO,gBAAgB;IAQ3B,YAAY,IAAgB,EAAE,OAA4B,EAAE,KAAa,EAAE,UAAyB;QALpG;;;;;WAAiB;QACjB;;;;;WAA6B;QAC7B;;;;;WAAc;QACd;;;;;WAA0B;QAGxB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,OAAO,GAAG,OAAO,CAAC;QACvB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;IAC/B,CAAC;IAED,mCAAmC;IACnC,GAAG,CAAC,IAAkB;QACpB,IAAI,IAAI,CAAC,OAAO,IAAI,IAAI;YAAE,OAAO,IAAI,CAAC;QACtC,0EAA0E;QAC1E,wEAAwE;QACxE,OAAO,YAAY,CAAC,IAAI,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;IAC1C,CAAC;IAED,qFAAqF;IACrF,QAAQ,CAAC,IAAU;QACjB,IAAI,IAAI,CAAC,OAAO,IAAI,IAAI;YAAE,OAAO,IAAI,CAAC;QACtC,MAAM,UAAU,GAAG,IAAI,CAAC,UAAU,CAAC,kBAAkB,CAAC,IAAI,CAAC,CAAC;QAC5D,MAAM,SAAS,GAAG,IAAI,CAAC,UAAU,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAErD,kFAAkF;QAClF,MAAM,WAAW,GACf,CAAC,UAAU,CAAC,KAAK,GAAG,SAAS,GAAG,CAAC,gBAAgB,CAAC,gBAAgB,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC;QAE3G,MAAM,YAAY,GAAG,UAAU,CAAC,eAAe,CAAC,WAAW,CAAC,CAAC;QAE7D,MAAM,SAAS,GAAG,YAAY,CAAC,IAAI,CAAC,OAAO,EAAE,YAAY,CAAC,SAAS,EAAE,CAAC,CAAC;QAEvE,uDAAuD;QACvD,MAAM,aAAa,GAAG,IAAI,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;QACnD,MAAM,UAAU,GAAG,YAAY,CAAC,KAAK,GAAG,YAAY,CAAC,MAAM,CAAC;QAC5D,MAAM,WAAW,GAAG,CAAC,GAAG,aAAa,GAAG,UAAU,CAAC;QAEnD,qEAAqE;QACrE,IAAI,WAAW,GAAG,IAAI;YAAE,OAAO,IAAI,CAAC;QAEpC,MAAM,OAAO,GAAG,mBAAmB,CAAC,CAAC,qBAAqB,CAAC,SAAS,CAAC,CAAC,CAA6B,CAAC;QACpG,OAAO,CAAC,GAAG,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,UAAU,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,EAAE,EAAE,EAAE,CAAC;QACzF,OAAO,OAAO,CAAC;IACjB,CAAC;IAED;;;;OAIG;IACH,SAAS;QACP,IAAI,IAAI,CAAC,OAAO,IAAI,IAAI;YAAE,OAAO,IAAI,CAAC;QACtC,MAAM,OAAO,GAAG,mBAAmB,CAAC,CAAC,qBAAqB,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAA6B,CAAC;QACvG,OAAO,CAAC,GAAG,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,UAAU,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,UAAU,CAAC,UAAU,CAAC,KAAK,EAAE,EAAE,EAAE,CAAC;QACzF,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAqB,EAAE,KAAa,EAAE,UAAyB;QAC/E,IAAI,IAAI,IAAI,IAAI;YAAE,OAAO,IAAI,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC;QAC7E,MAAM,GAAG,GAAG,MAAM,WAAW,CAAC,IAAI,CAAC,CAAC;QAEpC,IAAI,GAAG,CAAC,UAAU,KAAK,UAAU,CAAC,UAAU,CAAC,IAAI;YAC/C,OAAO,IAAI,gBAAgB,CAAC,IAAI,EAAE,GAAG,CAAC,OAAO,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC;QAEpE,MAAM,SAAS,GAAG,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,mBAAmB,CAClE,GAAG,CAAC,OAAO,EACX,UAAU,CAAC,GAAG,CAAC,UAAU,CAAC,UAAU,CAAC,CACtB,CAAC;QAElB,OAAO,IAAI,gBAAgB,CAAC,IAAI,EAAE,SAAS,EAAE,KAAK,EAAE,UAAU,CAAC,CAAC;IAClE,CAAC;IAED,0DAA0D;IAC1D,MAAM,CAAC,KAAK,CAAC,YAAY,CAAC,CAA2B,EAAE,UAAyB;QAC9E,IAAI,CAAC,IAAI,IAAI;YAAE,OAAO,IAAI,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QACtE,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;IACzD,CAAC;;AAhFD,kFAAkF;AAC3E;;;;WAAmB,GAAG;EAAN,CAAO"}

package/build/cogify/covering/tile.cover.d.ts

@@ -0,0 +1,38 @@
+import { Rgba } from '@basemaps/config';
+import { ConfigImageryTiff } from '@basemaps/config-loader';
+import { TileMatrixSet } from '@basemaps/geo';
+import { LogType } from '@basemaps/shared';
+import { Metrics } from '@linzjs/metrics';
+import { CogifyStacCollection, CogifyStacItem } from '../stac.js';
+import { CutlineOptimizer } from './cutline.js';
+export interface TileCoverContext {
+    /** Unique id for the covering */
+    id: string;
+    /** List of imagery to cover */
+    imagery: ConfigImageryTiff;
+    /** Cutline to apply */
+    cutline: CutlineOptimizer;
+    /** Output tile matrix */
+    tileMatrix: TileMatrixSet;
+    /** Optional metrics provider to track how long actions take */
+    metrics?: Metrics;
+    /** Optional logger to trace covering creation */
+    logger?: LogType;
+    /** GDAL configuration preset */
+    preset: string;
+    /** Optional color with which to replace all transparent COG pixels */
+    background?: Rgba;
+    /**
+     * Override the base zoom to store the output COGS as
+     */
+    targetZoomOffset?: number;
+}
+export interface TileCoverResult {
+    /** Stac collection for the imagery */
+    collection: CogifyStacCollection;
+    /** tiles to create */
+    items: CogifyStacItem[];
+    /** GeoJSON features of all the source imagery */
+    source: GeoJSON.FeatureCollection;
+}
+export declare function createTileCover(ctx: TileCoverContext): Promise<TileCoverResult>;

package/build/cogify/covering/tile.cover.js

@@ -0,0 +1,198 @@
+import { Bounds, Projection, ProjectionLoader, TileId } from '@basemaps/geo';
+import { fsa, urlToString } from '@basemaps/shared';
+import { CliDate, CliInfo } from '@basemaps/shared/build/cli/info.js';
+import { intersection, toFeatureCollection, union } from '@linzjs/geojson';
+import { Presets } from '../../preset.js';
+import { createFileStats } from '../stac.js';
+import { createCovering } from './covering.js';
+function getDateTime(ctx) {
+    const interval = ctx.imagery.collection?.extent?.temporal?.interval?.[0];
+    if (interval)
+        return { start: interval[0], end: interval[1] };
+    // TODO should we guess datetime
+    return { start: null, end: null };
+}
+/**
+ * Find a zoom level in the target tile matrix that is at least as good as resolution as `resolution`
+ *
+ * for some imagery sets this can mean a large difference in the source to target for example 1M resolution source imagery
+ * is resampled to 0.59m (z18) rather than the closer 1.19m (z17) `targetZoomOffset` can adjust the resolution
+ * to a nearer value using `-1` will convert a 1M target from 0.59m (z18) to 1.19m (z17)
+ *
+ * This can greatly reduce the size of output tiles
+ *
+ * @param tileMatrix target tile tile matrix to use
+ * @param resolution target resolution
+ * @param targetZoomOffset override the target base zoom, for instance -1 turns a z18 into z17
+ * @returns
+ */
+function getTargetBaseZoom(tileMatrix, resolution, targetZoomOffset) {
+    if (targetZoomOffset == null)
+        return Projection.getTiffResZoom(tileMatrix, resolution);
+    return Projection.getTiffResZoom(tileMatrix, resolution) + targetZoomOffset;
+}
+export async function createTileCover(ctx) {
+    // Ensure we have the projection loaded for the source imagery
+    await ProjectionLoader.load(ctx.imagery.projection);
+    // Find the zoom level that is at least as good as the source imagery
+    const targetBaseZoom = getTargetBaseZoom(ctx.tileMatrix, ctx.imagery.gsd, ctx.targetZoomOffset);
+    // The base zoom is 256x256 pixels at its resolution, we are trying to find a image that is <32k pixels wide/high
+    // zooming out 7 levels converts a 256x256 image into 32k x 32k image
+    // 256 * 2 ** 7 = 32,768 - 256x256 tile
+    // 512 * 2 ** 6 = 32,768 - 512x512 tile
+    const optimalCoveringZoom = Math.max(1, targetBaseZoom - 7); // z12 from z19
+    ctx.logger?.debug({ targetBaseZoom, cogOverZoom: optimalCoveringZoom }, 'Imagery:ZoomLevel');
+    const sourceBounds = projectPolygon(polygonFromBounds(ctx.imagery.files), ctx.imagery.projection, ctx.tileMatrix.projection.code);
+    ctx.logger?.debug('Cutline:Apply');
+    ctx.metrics?.start('cutline:apply');
+    const dateTime = getDateTime(ctx);
+    // Convert the source imagery to a geojson
+    const sourceGeoJson = ctx.imagery.files.map((file) => {
+        return Projection.get(ctx.imagery.projection).boundsToGeoJsonFeature(file, { url: file.name });
+    });
+    ctx.logger?.info({ zoom: optimalCoveringZoom }, 'Imagery:Covering:Start');
+    // Cover the source imagery in tiles
+    ctx.metrics?.start('covering:create');
+    const covering = createCovering({
+        tileMatrix: ctx.tileMatrix,
+        source: sourceBounds,
+        targetZoom: optimalCoveringZoom,
+        baseZoom: targetBaseZoom,
+        cutline: ctx.cutline,
+        metrics: ctx.metrics,
+        logger: ctx.logger,
+    });
+    ctx.metrics?.end('covering:create');
+    if (covering.length === 0)
+        throw new Error('Unable to create tile covering, no tiles created.');
+    ctx.logger?.info({ tiles: covering.length, zoom: optimalCoveringZoom }, 'Imagery:Covering:Created');
+    const imageryBounds = ctx.imagery.files.map((f) => {
+        const polygon = Bounds.fromJson(f).toPolygon();
+        return { ...f, polygon };
+    });
+    ctx.metrics?.start('covering:polygon');
+    const items = [];
+    for (const tile of covering) {
+        const bounds = ctx.tileMatrix.tileToSourceBounds(tile);
+        // Scale the tile bounds slightly to ensure we get all relevant imagery
+        const scaledBounds = bounds.scaleFromCenter(1.05);
+        const tileBounds = Projection.get(ctx.tileMatrix).projectMultipolygon([scaledBounds.toPolygon()], Projection.get(ctx.imagery.projection));
+        const source = imageryBounds.filter((f) => intersection(tileBounds, f.polygon).length > 0);
+        const feature = Projection.get(ctx.tileMatrix).boundsToGeoJsonFeature(bounds);
+        const tileId = TileId.fromTile(tile);
+        const item = {
+            id: `${ctx.id}/${tileId}`,
+            type: 'Feature',
+            collection: ctx.id,
+            stac_version: '1.0.0',
+            stac_extensions: [],
+            geometry: feature.geometry,
+            bbox: Projection.get(ctx.tileMatrix).boundsToWgs84BoundingBox(bounds),
+            links: [
+                { href: `./${tileId}.json`, rel: 'self' },
+                { href: './collection.json', rel: 'collection' },
+                { href: './collection.json', rel: 'parent' },
+            ],
+            properties: {
+                datetime: dateTime.start ? null : CliDate,
+                start_datetime: dateTime.start ?? undefined,
+                end_datetime: dateTime.end ?? undefined,
+                'proj:epsg': ctx.tileMatrix.projection.code,
+                'linz_basemaps:options': {
+                    preset: ctx.preset,
+                    ...Presets[ctx.preset].options,
+                    tile,
+                    tileMatrix: ctx.tileMatrix.identifier,
+                    sourceEpsg: ctx.imagery.projection,
+                    zoomLevel: targetBaseZoom,
+                },
+                'linz_basemaps:generated': {
+                    package: CliInfo.package,
+                    hash: CliInfo.hash,
+                    version: CliInfo.version,
+                    datetime: CliDate,
+                },
+            },
+            assets: {},
+        };
+        // Add the background color if it exists
+        if (ctx.background)
+            item.properties['linz_basemaps:options'].background = ctx.background;
+        // Add the source imagery as a STAC Link
+        for (const src of source) {
+            const srcLink = {
+                href: new URL(src.name, ctx.imagery.url).href,
+                rel: 'linz_basemaps:source',
+                type: 'image/tiff; application=geotiff;',
+                'linz_basemaps:source_height': src.height,
+                'linz_basemaps:source_width': src.width,
+            };
+            item.links.push(srcLink);
+        }
+        // Add the cutline as a STAC Link if it exists
+        if (ctx.cutline.path) {
+            const cutLink = {
+                href: urlToString(ctx.cutline.path),
+                rel: 'linz_basemaps:cutline',
+                blend: ctx.cutline.blend,
+            };
+            item.links.push(cutLink);
+        }
+        items.push(item);
+    }
+    ctx.metrics?.end('covering:polygon');
+    const collection = {
+        id: ctx.id,
+        type: 'Collection',
+        stac_version: '1.0.0',
+        stac_extensions: [],
+        license: ctx.imagery.collection?.license ?? 'CC-BY-4.0',
+        title: ctx.imagery.title,
+        description: ctx.imagery.collection?.description ?? 'Missing source STAC',
+        providers: ctx.imagery.collection?.providers,
+        extent: {
+            spatial: { bbox: [Projection.get(ctx.imagery.projection).boundsToWgs84BoundingBox(ctx.imagery.bounds)] },
+            // Default  the temporal time today if no times were found as it is required for STAC
+            temporal: { interval: dateTime.start ? [[dateTime.start, dateTime.end]] : [[CliDate, null]] },
+        },
+        links: items.map((item) => {
+            const tile = item.properties['linz_basemaps:options'].tile;
+            if (tile == null)
+                throw new Error('Tile missing from item');
+            const tileId = TileId.fromTile(tile);
+            return { href: `./${tileId}.json`, rel: 'item', type: 'application/json' };
+        }),
+    };
+    // Add a self link to the links
+    collection.links.unshift({ rel: 'self', href: './collection.json', type: 'application/json' });
+    // Include a link back to the source collection
+    if (ctx.imagery.collection) {
+        const target = new URL('collection.json', ctx.imagery.url);
+        const stac = await fsa.read(target);
+        collection.links.push({
+            rel: 'linz_basemaps:source_collection',
+            href: target.href,
+            type: 'application/json',
+            ...createFileStats(stac),
+        });
+    }
+    return { collection, items, source: toFeatureCollection(sourceGeoJson) };
+}
+/** Convert a list of bounding boxes into a multipolygon */
+function polygonFromBounds(bounds, scale = 1 + 1e-8) {
+    const srcPoly = [];
+    // merge imagery bounds
+    for (const image of bounds) {
+        srcPoly.push(Bounds.fromJson(image).scaleFromCenter(scale).toPolygon());
+    }
+    return union(srcPoly);
+}
+/** project a polygon from the source projection to a target projection  */
+function projectPolygon(p, sourceProjection, targetProjection) {
+    if (sourceProjection === targetProjection)
+        return p;
+    const sourceProj = Projection.get(sourceProjection);
+    const targetProj = Projection.get(targetProjection);
+    return sourceProj.projectMultipolygon(p, targetProj);
+}
+//# sourceMappingURL=tile.cover.js.map

package/build/cogify/covering/tile.cover.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"tile.cover.js","sourceRoot":"","sources":["../../../src/cogify/covering/tile.cover.ts"],"names":[],"mappings":"AAEA,OAAO,EAAe,MAAM,EAAY,UAAU,EAAE,gBAAgB,EAAE,MAAM,EAAiB,MAAM,eAAe,CAAC;AACnH,OAAO,EAAE,GAAG,EAAW,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAC7D,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,oCAAoC,CAAC;AACtE,OAAO,EAAE,YAAY,EAAgB,mBAAmB,EAAE,KAAK,EAAE,MAAM,iBAAiB,CAAC;AAIzF,OAAO,EAAE,OAAO,EAAE,MAAM,iBAAiB,CAAC;AAC1C,OAAO,EAA6E,eAAe,EAAE,MAAM,YAAY,CAAC;AACxH,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAkC/C,SAAS,WAAW,CAAC,GAAqB;IACxC,MAAM,QAAQ,GAAG,GAAG,CAAC,OAAO,CAAC,UAAU,EAAE,MAAM,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC;IAEzE,IAAI,QAAQ;QAAE,OAAO,EAAE,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;IAE9D,gCAAgC;IAChC,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC;AACpC,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,SAAS,iBAAiB,CAAC,UAAyB,EAAE,UAAkB,EAAE,gBAAyB;IACjG,IAAI,gBAAgB,IAAI,IAAI;QAAE,OAAO,UAAU,CAAC,cAAc,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;IACvF,OAAO,UAAU,CAAC,cAAc,CAAC,UAAU,EAAE,UAAU,CAAC,GAAG,gBAAgB,CAAC;AAC9E,CAAC;AAED,MAAM,CAAC,KAAK,UAAU,eAAe,CAAC,GAAqB;IACzD,8DAA8D;IAC9D,MAAM,gBAAgB,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC;IAEpD,qEAAqE;IACrE,MAAM,cAAc,GAAG,iBAAiB,CAAC,GAAG,CAAC,UAAU,EAAE,GAAG,CAAC,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,gBAAgB,CAAC,CAAC;IAEhG,iHAAiH;IACjH,qEAAqE;IACrE,uCAAuC;IACvC,uCAAuC;IACvC,MAAM,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,cAAc,GAAG,CAAC,CAAC,CAAC,CAAC,eAAe;IAC5E,GAAG,CAAC,MAAM,EAAE,KAAK,CAAC,EAAE,cAAc,EAAE,WAAW,EAAE,mBAAmB,EAAE,EAAE,mBAAmB,CAAC,CAAC;IAE7F,MAAM,YAAY,GAAG,cAAc,CACjC,iBAAiB,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,EACpC,GAAG,CAAC,OAAO,CAAC,UAAU,EACtB,GAAG,CAAC,UAAU,CAAC,UAAU,CAAC,IAAI,CAC/B,CAAC;IACF,GAAG,CAAC,MAAM,EAAE,KAAK,CAAC,eAAe,CAAC,CAAC;IACnC,GAAG,CAAC,OAAO,EAAE,KAAK,CAAC,eAAe,CAAC,CAAC;IAEpC,MAAM,QAAQ,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC;IAElC,0CAA0C;IAC1C,MAAM,aAAa,GAAG,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE;QACnD,OAAO,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,sBAAsB,CAAC,IAAI,EAAE,EAAE,GAAG,EAAE,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;IACjG,CAAC,CAAC,CAAC;IACH,GAAG,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,IAAI,EAAE,mBAAmB,EAAE,EAAE,wBAAwB,CAAC,CAAC;IAE1E,oCAAoC;IACpC,GAAG,CAAC,OAAO,EAAE,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,cAAc,CAAC;QAC9B,UAAU,EAAE,GAAG,CAAC,UAAU;QAC1B,MAAM,EAAE,YAAY;QACpB,UAAU,EAAE,mBAAmB;QAC/B,QAAQ,EAAE,cAAc;QACxB,OAAO,EAAE,GAAG,CAAC,OAAO;QACpB,OAAO,EAAE,GAAG,CAAC,OAAO;QACpB,MAAM,EAAE,GAAG,CAAC,MAAM;KACnB,CAAC,CAAC;IACH,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,iBAAiB,CAAC,CAAC;IAEpC,IAAI,QAAQ,CAAC,MAAM,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IAChG,GAAG,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,mBAAmB,EAAE,EAAE,0BAA0B,CAAC,CAAC;IAEpG,MAAM,aAAa,GAAG,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;QAChD,MAAM,OAAO,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QAC/C,OAAO,EAAE,GAAG,CAAC,EAAE,OAAO,EAAE,CAAC;IAC3B,CAAC,CAAC,CAAC;IAEH,GAAG,CAAC,OAAO,EAAE,KAAK,CAAC,kBAAkB,CAAC,CAAC;IACvC,MAAM,KAAK,GAAqB,EAAE,CAAC;IACnC,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,GAAG,CAAC,UAAU,CAAC,kBAAkB,CAAC,IAAI,CAAC,CAAC;QAEvD,uEAAuE;QACvE,MAAM,YAAY,GAAG,MAAM,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC;QAClD,MAAM,UAAU,GAAG,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,mBAAmB,CACnE,CAAC,YAAY,CAAC,SAAS,EAAE,CAAC,EAC1B,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,CACvB,CAAC;QAElB,MAAM,MAAM,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,YAAY,CAAC,UAAU,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAE3F,MAAM,OAAO,GAAG,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,sBAAsB,CAAC,MAAM,CAAC,CAAC;QAE9E,MAAM,MAAM,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QAErC,MAAM,IAAI,GAAmB;YAC3B,EAAE,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,MAAM,EAAE;YACzB,IAAI,EAAE,SAAS;YACf,UAAU,EAAE,GAAG,CAAC,EAAE;YAClB,YAAY,EAAE,OAAO;YACrB,eAAe,EAAE,EAAE;YACnB,QAAQ,EAAE,OAAO,CAAC,QAA0B;YAC5C,IAAI,EAAE,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,wBAAwB,CAAC,MAAM,CAAC;YACrE,KAAK,EAAE;gBACL,EAAE,IAAI,EAAE,KAAK,MAAM,OAAO,EAAE,GAAG,EAAE,MAAM,EAAE;gBACzC,EAAE,IAAI,EAAE,mBAAmB,EAAE,GAAG,EAAE,YAAY,EAAE;gBAChD,EAAE,IAAI,EAAE,mBAAmB,EAAE,GAAG,EAAE,QAAQ,EAAE;aAC7C;YACD,UAAU,EAAE;gBACV,QAAQ,EAAE,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,OAAO;gBACzC,cAAc,EAAE,QAAQ,CAAC,KAAK,IAAI,SAAS;gBAC3C,YAAY,EAAE,QAAQ,CAAC,GAAG,IAAI,SAAS;gBACvC,WAAW,EAAE,GAAG,CAAC,UAAU,CAAC,UAAU,CAAC,IAAI;gBAC3C,uBAAuB,EAAE;oBACvB,MAAM,EAAE,GAAG,CAAC,MAAM;oBAClB,GAAG,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,OAAO;oBAC9B,IAAI;oBACJ,UAAU,EAAE,GAAG,CAAC,UAAU,CAAC,UAAU;oBACrC,UAAU,EAAE,GAAG,CAAC,OAAO,CAAC,UAAU;oBAClC,SAAS,EAAE,cAAc;iBAC1B;gBACD,yBAAyB,EAAE;oBACzB,OAAO,EAAE,OAAO,CAAC,OAAO;oBACxB,IAAI,EAAE,OAAO,CAAC,IAAI;oBAClB,OAAO,EAAE,OAAO,CAAC,OAAO;oBACxB,QAAQ,EAAE,OAAO;iBAClB;aACF;YACD,MAAM,EAAE,EAAE;SACX,CAAC;QAEF,wCAAwC;QACxC,IAAI,GAAG,CAAC,UAAU;YAAE,IAAI,CAAC,UAAU,CAAC,uBAAuB,CAAC,CAAC,UAAU,GAAG,GAAG,CAAC,UAAU,CAAC;QAEzF,wCAAwC;QACxC,KAAK,MAAM,GAAG,IAAI,MAAM,EAAE,CAAC;YACzB,MAAM,OAAO,GAAqB;gBAChC,IAAI,EAAE,IAAI,GAAG,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,IAAI;gBAC7C,GAAG,EAAE,sBAAsB;gBAC3B,IAAI,EAAE,kCAAkC;gBACxC,6BAA6B,EAAE,GAAG,CAAC,MAAM;gBACzC,4BAA4B,EAAE,GAAG,CAAC,KAAK;aACxC,CAAC;YACF,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC3B,CAAC;QAED,8CAA8C;QAC9C,IAAI,GAAG,CAAC,OAAO,CAAC,IAAI,EAAE,CAAC;YACrB,MAAM,OAAO,GAAsB;gBACjC,IAAI,EAAE,WAAW,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC;gBACnC,GAAG,EAAE,uBAAuB;gBAC5B,KAAK,EAAE,GAAG,CAAC,OAAO,CAAC,KAAK;aACzB,CAAC;YACF,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAC3B,CAAC;QAED,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACnB,CAAC;IACD,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,kBAAkB,CAAC,CAAC;IAErC,MAAM,UAAU,GAAyB;QACvC,EAAE,EAAE,GAAG,CAAC,EAAE;QACV,IAAI,EAAE,YAAY;QAClB,YAAY,EAAE,OAAO;QACrB,eAAe,EAAE,EAAE;QACnB,OAAO,EAAE,GAAG,CAAC,OAAO,CAAC,UAAU,EAAE,OAAO,IAAI,WAAW;QACvD,KAAK,EAAE,GAAG,CAAC,OAAO,CAAC,KAAK;QACxB,WAAW,EAAE,GAAG,CAAC,OAAO,CAAC,UAAU,EAAE,WAAW,IAAI,qBAAqB;QACzE,SAAS,EAAE,GAAG,CAAC,OAAO,CAAC,UAAU,EAAE,SAAS;QAC5C,MAAM,EAAE;YACN,OAAO,EAAE,EAAE,IAAI,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,wBAAwB,CAAC,GAAG,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE;YACxG,qFAAqF;YACrF,QAAQ,EAAE,EAAE,QAAQ,EAAE,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC,EAAE;SAC9F;QACD,KAAK,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE;YACxB,MAAM,IAAI,GAAG,IAAI,CAAC,UAAU,CAAC,uBAAuB,CAAC,CAAC,IAAI,CAAC;YAC3D,IAAI,IAAI,IAAI,IAAI;gBAAE,MAAM,IAAI,KAAK,CAAC,wBAAwB,CAAC,CAAC;YAC5D,MAAM,MAAM,GAAG,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;YACrC,OAAO,EAAE,IAAI,EAAE,KAAK,MAAM,OAAO,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,kBAAkB,EAAE,CAAC;QAC7E,CAAC,CAAC;KACH,CAAC;IAEF,+BAA+B;IAC/B,UAAU,CAAC,KAAK,CAAC,OAAO,CAAC,EAAE,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,mBAAmB,EAAE,IAAI,EAAE,kBAAkB,EAAE,CAAC,CAAC;IAC/F,+CAA+C;IAC/C,IAAI,GAAG,CAAC,OAAO,CAAC,UAAU,EAAE,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,iBAAiB,EAAE,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;QAC3D,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACpC,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC;YACpB,GAAG,EAAE,iCAAiC;YACtC,IAAI,EAAE,MAAM,CAAC,IAAI;YACjB,IAAI,EAAE,kBAAkB;YACxB,GAAG,eAAe,CAAC,IAAI,CAAC;SACzB,CAAC,CAAC;IACL,CAAC;IAED,OAAO,EAAE,UAAU,EAAE,KAAK,EAAE,MAAM,EAAE,mBAAmB,CAAC,aAAa,CAAC,EAAE,CAAC;AAC3E,CAAC;AAED,2DAA2D;AAC3D,SAAS,iBAAiB,CAAC,MAAqB,EAAE,QAAgB,CAAC,GAAG,IAAI;IACxE,MAAM,OAAO,GAAiB,EAAE,CAAC;IAEjC,uBAAuB;IACvB,KAAK,MAAM,KAAK,IAAI,MAAM,EAAE,CAAC;QAC3B,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,eAAe,CAAC,KAAK,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC;IAC1E,CAAC;IACD,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC;AACxB,CAAC;AAED,2EAA2E;AAC3E,SAAS,cAAc,CAAC,CAAe,EAAE,gBAA0B,EAAE,gBAA0B;IAC7F,IAAI,gBAAgB,KAAK,gBAAgB;QAAE,OAAO,CAAC,CAAC;IACpD,MAAM,UAAU,GAAG,UAAU,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC;IACpD,MAAM,UAAU,GAAG,UAAU,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC;IACpD,OAAO,UAAU,CAAC,mBAAmB,CAAC,CAAC,EAAE,UAAU,CAAiB,CAAC;AACvE,CAAC"}

package/build/cogify/gdal/gdal.command.d.ts

@@ -0,0 +1,27 @@
+import { Rgba } from '@basemaps/config';
+import { EpsgCode } from '@basemaps/geo';
+import { CogifyCreationOptions } from '../stac.js';
+import { GdalCommand } from './gdal.runner.js';
+export declare function gdalBuildVrt(targetVrt: URL, source: URL[], addalpha?: boolean): GdalCommand;
+export declare function gdalBuildVrtWarp(targetVrt: URL, sourceVrt: URL, sourceProjection: EpsgCode, cutline: {
+    url: URL | null;
+    blend: number;
+}, opt: CogifyCreationOptions): GdalCommand;
+export declare function gdalBuildCog(targetTiff: URL, sourceVrt: URL, opt: CogifyCreationOptions): GdalCommand;
+/**
+ * Creates an empty tiff where all pixel values are set to the given color.
+ * Used to force a background so that there are no empty pixels in the final COG.
+ *
+ * @param targetTiff the file path and name for the created tiff
+ * @param color the color to set all pixel values
+ * @param opt a CogifyCreationOptions object
+ *
+ * @returns a 'gdal_create' GdalCommand object
+ */
+export declare function gdalCreate(targetTiff: URL, color: Rgba, opt: CogifyCreationOptions): GdalCommand;
+/**
+ * Build a topographic mapsheet cog
+ *
+ * This is specific configuration to LINZ's topo50 and 250 mapsheets
+ */
+export declare function gdalBuildTopoRasterCommands(targetTiff: URL, sourceVrt: URL, opt: CogifyCreationOptions, width: number, height: number): GdalCommand;

package/build/cogify/gdal/gdal.command.js

@@ -0,0 +1,175 @@
+import { Epsg, TileMatrixSets } from '@basemaps/geo';
+import { urlToString } from '@basemaps/shared';
+import { Presets } from '../../preset.js';
+const isPowerOfTwo = (x) => (x & (x - 1)) === 0;
+/**
+ * Topographic mapsheets are rendered generally at 1:600 dpi,
+ *
+ * A topo50 mapsheet 1:600dpi does not perfectly align to full pixels and approximatly 1.7 pixels
+ * of empty space is rendered at the edge of every mapsheet.
+ *
+ */
+const DefaultTrimPixelRight = 1.7; // 1.7 pixels to trim from the right side of the topo raster imagery
+export function gdalBuildVrt(targetVrt, source, addalpha) {
+    if (source.length === 0)
+        throw new Error('No source files given for :' + targetVrt.href);
+    return {
+        output: targetVrt,
+        command: 'gdalbuildvrt',
+        args: [addalpha ? ['-addalpha'] : undefined, urlToString(targetVrt), ...source.map(urlToString)]
+            .filter((f) => f != null)
+            .flat()
+            .map(String),
+    };
+}
+export function gdalBuildVrtWarp(targetVrt, sourceVrt, sourceProjection, cutline, opt) {
+    const tileMatrix = TileMatrixSets.find(opt.tileMatrix);
+    if (tileMatrix == null)
+        throw new Error('Unable to find tileMatrix: ' + opt.tileMatrix);
+    const targetResolution = tileMatrix.pixelScale(opt.zoomLevel);
+    return {
+        output: targetVrt,
+        command: 'gdalwarp',
+        args: [
+            ['-of', 'vrt'], // Output as a VRT
+            // ['-co', 'compress=lzw'],
+            // ['-co', 'bigtiff=yes'],
+            '-multi', // Mutithread IO
+            ['-wo', 'NUM_THREADS=ALL_CPUS'], // Multithread the warp
+            ['-s_srs', Epsg.get(sourceProjection).toEpsgString()], // Source EPSG
+            ['-t_srs', tileMatrix.projection.toEpsgString()], // Target EPSG
+            ['-tr', targetResolution, targetResolution],
+            opt.warpResampling ? ['-r', opt.warpResampling] : undefined,
+            cutline.url ? ['-cutline', urlToString(cutline.url), '-cblend', cutline.blend] : undefined,
+            urlToString(sourceVrt),
+            urlToString(targetVrt),
+        ]
+            .filter((f) => f != null)
+            .flat()
+            .map(String),
+    };
+}
+export function gdalBuildCog(targetTiff, sourceVrt, opt) {
+    const cfg = { ...Presets[opt.preset], ...opt };
+    const tileMatrix = TileMatrixSets.find(cfg.tileMatrix);
+    if (tileMatrix == null)
+        throw new Error('Unable to find tileMatrix: ' + cfg.tileMatrix);
+    const bounds = tileMatrix.tileToSourceBounds(cfg.tile);
+    const tileExtent = [
+        Math.min(bounds.x, bounds.right),
+        Math.min(bounds.y, bounds.bottom),
+        Math.max(bounds.x, bounds.right),
+        Math.max(bounds.y, bounds.bottom),
+    ];
+    const targetResolution = tileMatrix.pixelScale(cfg.zoomLevel);
+    return {
+        command: 'gdal_translate',
+        output: targetTiff,
+        args: [
+            ['-of', 'COG'],
+            ['-co', 'NUM_THREADS=ALL_CPUS'], // Use all CPUS
+            ['--config', 'GDAL_NUM_THREADS', 'all_cpus'], // Also required to NUM_THREADS till gdal 3.7.x
+            ['-co', 'BIGTIFF=IF_NEEDED'], // BigTiff is somewhat slower and most (All?) of the COGS should be well below 4GB
+            ['-co', 'ADD_ALPHA=YES'],
+            /**
+             *  GDAL will recompress existing overviews if they exist which will compound
+             *  any lossly compression on the overview, so compute new overviews instead
+             */
+            ['-co', 'OVERVIEWS=IGNORE_EXISTING'],
+            ['-co', `BLOCKSIZE=${cfg.blockSize}`],
+            // ['-co', 'RESAMPLING=cubic'],
+            ['-co', `WARP_RESAMPLING=${cfg.warpResampling}`],
+            ['-co', `OVERVIEW_RESAMPLING=${cfg.overviewResampling}`],
+            ['-co', `COMPRESS=${cfg.compression}`],
+            cfg.quality ? ['-co', `QUALITY=${cfg.quality}`] : undefined,
+            cfg.maxZError ? ['-co', `MAX_Z_ERROR=${cfg.maxZError}`] : undefined,
+            cfg.maxZErrorOverview ? ['-co', `MAX_Z_ERROR_OVERVIEW=${cfg.maxZErrorOverview}`] : undefined,
+            ['-co', 'SPARSE_OK=YES'],
+            ['-co', `TARGET_SRS=${tileMatrix.projection.toEpsgString()}`],
+            ['-co', `EXTENT=${tileExtent.join(',')},`],
+            ['-tr', targetResolution, targetResolution],
+            urlToString(sourceVrt),
+            urlToString(targetTiff),
+        ]
+            .filter((f) => f != null)
+            .flat()
+            .map(String),
+    };
+}
+/**
+ * Creates an empty tiff where all pixel values are set to the given color.
+ * Used to force a background so that there are no empty pixels in the final COG.
+ *
+ * @param targetTiff the file path and name for the created tiff
+ * @param color the color to set all pixel values
+ * @param opt a CogifyCreationOptions object
+ *
+ * @returns a 'gdal_create' GdalCommand object
+ */
+export function gdalCreate(targetTiff, color, opt) {
+    const cfg = { ...Presets[opt.preset], ...opt };
+    const tileMatrix = TileMatrixSets.find(cfg.tileMatrix);
+    if (tileMatrix == null)
+        throw new Error('Unable to find tileMatrix: ' + cfg.tileMatrix);
+    const bounds = tileMatrix.tileToSourceBounds(cfg.tile);
+    const pixelScale = tileMatrix.pixelScale(cfg.zoomLevel);
+    const size = Math.round(bounds.width / pixelScale);
+    // if the value of 'size' is not a power of 2
+    if (!isPowerOfTwo(size))
+        throw new Error('Size did not compute to a power of 2');
+    return {
+        command: 'gdal_create',
+        output: targetTiff,
+        args: [
+            ['-of', 'GTiff'],
+            ['-outsize', size, size], // set the size to match that of the final COG
+            ['-bands', '4'],
+            ['-burn', `${color.r} ${color.g} ${color.b} ${color.alpha}`], // set all pixel values to the given color
+            ['-a_srs', tileMatrix.projection.toEpsgString()],
+            ['-a_ullr', bounds.x, bounds.bottom, bounds.right, bounds.y],
+            ['-co', 'COMPRESS=LZW'],
+            urlToString(targetTiff),
+        ]
+            .filter((f) => f != null)
+            .flat()
+            .map(String),
+    };
+}
+/**
+ * Build a topographic mapsheet cog
+ *
+ * This is specific configuration to LINZ's topo50 and 250 mapsheets
+ */
+export function gdalBuildTopoRasterCommands(targetTiff, sourceVrt, opt, width, height) {
+    const command = {
+        command: 'gdal_translate',
+        output: targetTiff,
+        args: [
+            ['-q'], // Supress non-error output
+            ['-stats'], // Force stats (re)computation
+            ['-of', 'COG'], // Output format
+            ['-srcwin', '0', '0', `${width - DefaultTrimPixelRight}`, `${height}`],
+            ['-a_srs', `EPSG:${opt.sourceEpsg}`],
+            // https://gdal.org/en/latest/drivers/raster/cog.html#creation-options
+            ['-co', 'BIGTIFF=NO'],
+            ['-co', 'BLOCKSIZE=512'],
+            ['-co', 'COMPRESS=WEBP'],
+            ['-co', 'NUM_THREADS=ALL_CPUS'], // Use all CPUS
+            ['-co', 'OVERVIEW_COMPRESS=WEBP'],
+            ['-co', 'OVERVIEWS=IGNORE_EXISTING'],
+            ['-co', 'OVERVIEW_QUALITY=90'],
+            ['-co', 'OVERVIEW_RESAMPLING=LANCZOS'],
+            ['-co', 'QUALITY=100'],
+            ['-co', 'SPARSE_OK=TRUE'], // Allow for sparse writes
+            // https://gdal.org/en/latest/drivers/raster/cog.html#reprojection-related-creation-options
+            ['-co', 'ADD_ALPHA=YES'],
+            urlToString(sourceVrt),
+            urlToString(targetTiff),
+        ]
+            .filter((f) => f != null)
+            .flat()
+            .map(String),
+    };
+    return command;
+}
+//# sourceMappingURL=gdal.command.js.map

package/build/cogify/gdal/gdal.command.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"gdal.command.js","sourceRoot":"","sources":["../../../src/cogify/gdal/gdal.command.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,IAAI,EAAY,cAAc,EAAE,MAAM,eAAe,CAAC;AAC/D,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AAE/C,OAAO,EAAE,OAAO,EAAE,MAAM,iBAAiB,CAAC;AAI1C,MAAM,YAAY,GAAG,CAAC,CAAS,EAAW,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;AAEjE;;;;;;GAMG;AACH,MAAM,qBAAqB,GAAG,GAAG,CAAC,CAAC,oEAAoE;AAEvG,MAAM,UAAU,YAAY,CAAC,SAAc,EAAE,MAAa,EAAE,QAAkB;IAC5E,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;IACzF,OAAO;QACL,MAAM,EAAE,SAAS;QACjB,OAAO,EAAE,cAAc;QACvB,IAAI,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,WAAW,CAAC,SAAS,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,WAAW,CAAC,CAAC;aAC7F,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC;aACxB,IAAI,EAAE;aACN,GAAG,CAAC,MAAM,CAAC;KACf,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,gBAAgB,CAC9B,SAAc,EACd,SAAc,EACd,gBAA0B,EAC1B,OAA2C,EAC3C,GAA0B;IAE1B,MAAM,UAAU,GAAG,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACvD,IAAI,UAAU,IAAI,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC;IACxF,MAAM,gBAAgB,GAAG,UAAU,CAAC,UAAU,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IAE9D,OAAO;QACL,MAAM,EAAE,SAAS;QACjB,OAAO,EAAE,UAAU;QACnB,IAAI,EAAE;YACJ,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,kBAAkB;YAClC,2BAA2B;YAC3B,0BAA0B;YAC1B,QAAQ,EAAE,gBAAgB;YAC1B,CAAC,KAAK,EAAE,sBAAsB,CAAC,EAAE,uBAAuB;YACxD,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,gBAAgB,CAAC,CAAC,YAAY,EAAE,CAAC,EAAE,cAAc;YACrE,CAAC,QAAQ,EAAE,UAAU,CAAC,UAAU,CAAC,YAAY,EAAE,CAAC,EAAE,cAAc;YAChE,CAAC,KAAK,EAAE,gBAAgB,EAAE,gBAAgB,CAAC;YAC3C,GAAG,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,SAAS;YAC3D,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,UAAU,EAAE,WAAW,CAAC,OAAO,CAAC,GAAG,CAAC,EAAE,SAAS,EAAE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,SAAS;YAC1F,WAAW,CAAC,SAAS,CAAC;YACtB,WAAW,CAAC,SAAS,CAAC;SACvB;aACE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC;aACxB,IAAI,EAAE;aACN,GAAG,CAAC,MAAM,CAAC;KACf,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,UAAe,EAAE,SAAc,EAAE,GAA0B;IACtF,MAAM,GAAG,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;IAC/C,MAAM,UAAU,GAAG,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACvD,IAAI,UAAU,IAAI,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC;IAExF,MAAM,MAAM,GAAG,UAAU,CAAC,kBAAkB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACvD,MAAM,UAAU,GAAG;QACjB,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC;QAChC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC;QAChC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,CAAC;KAClC,CAAC;IAEF,MAAM,gBAAgB,GAAG,UAAU,CAAC,UAAU,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IAE9D,OAAO;QACL,OAAO,EAAE,gBAAgB;QACzB,MAAM,EAAE,UAAU;QAClB,IAAI,EAAE;YACJ,CAAC,KAAK,EAAE,KAAK,CAAC;YACd,CAAC,KAAK,EAAE,sBAAsB,CAAC,EAAE,eAAe;YAChD,CAAC,UAAU,EAAE,kBAAkB,EAAE,UAAU,CAAC,EAAE,+CAA+C;YAC7F,CAAC,KAAK,EAAE,mBAAmB,CAAC,EAAE,kFAAkF;YAChH,CAAC,KAAK,EAAE,eAAe,CAAC;YACxB;;;eAGG;YACH,CAAC,KAAK,EAAE,2BAA2B,CAAC;YACpC,CAAC,KAAK,EAAE,aAAa,GAAG,CAAC,SAAS,EAAE,CAAC;YACrC,+BAA+B;YAC/B,CAAC,KAAK,EAAE,mBAAmB,GAAG,CAAC,cAAc,EAAE,CAAC;YAChD,CAAC,KAAK,EAAE,uBAAuB,GAAG,CAAC,kBAAkB,EAAE,CAAC;YACxD,CAAC,KAAK,EAAE,YAAY,GAAG,CAAC,WAAW,EAAE,CAAC;YACtC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,WAAW,GAAG,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,SAAS;YAC3D,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,eAAe,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,CAAC,SAAS;YACnE,GAAG,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,wBAAwB,GAAG,CAAC,iBAAiB,EAAE,CAAC,CAAC,CAAC,CAAC,SAAS;YAC5F,CAAC,KAAK,EAAE,eAAe,CAAC;YACxB,CAAC,KAAK,EAAE,cAAc,UAAU,CAAC,UAAU,CAAC,YAAY,EAAE,EAAE,CAAC;YAC7D,CAAC,KAAK,EAAE,UAAU,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC;YAC1C,CAAC,KAAK,EAAE,gBAAgB,EAAE,gBAAgB,CAAC;YAC3C,WAAW,CAAC,SAAS,CAAC;YACtB,WAAW,CAAC,UAAU,CAAC;SACxB;aACE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC;aACxB,IAAI,EAAE;aACN,GAAG,CAAC,MAAM,CAAC;KACf,CAAC;AACJ,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,UAAU,CAAC,UAAe,EAAE,KAAW,EAAE,GAA0B;IACjF,MAAM,GAAG,GAAG,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;IAE/C,MAAM,UAAU,GAAG,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IACvD,IAAI,UAAU,IAAI,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,6BAA6B,GAAG,GAAG,CAAC,UAAU,CAAC,CAAC;IAExF,MAAM,MAAM,GAAG,UAAU,CAAC,kBAAkB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACvD,MAAM,UAAU,GAAG,UAAU,CAAC,UAAU,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACxD,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,GAAG,UAAU,CAAC,CAAC;IAEnD,6CAA6C;IAC7C,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;IAEjF,OAAO;QACL,OAAO,EAAE,aAAa;QACtB,MAAM,EAAE,UAAU;QAClB,IAAI,EAAE;YACJ,CAAC,KAAK,EAAE,OAAO,CAAC;YAChB,CAAC,UAAU,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,8CAA8C;YACxE,CAAC,QAAQ,EAAE,GAAG,CAAC;YACf,CAAC,OAAO,EAAE,GAAG,KAAK,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,IAAI,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,0CAA0C;YACxG,CAAC,QAAQ,EAAE,UAAU,CAAC,UAAU,CAAC,YAAY,EAAE,CAAC;YAChD,CAAC,SAAS,EAAE,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAC;YAC5D,CAAC,KAAK,EAAE,cAAc,CAAC;YACvB,WAAW,CAAC,UAAU,CAAC;SACxB;aACE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC;aACxB,IAAI,EAAE;aACN,GAAG,CAAC,MAAM,CAAC;KACf,CAAC;AACJ,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,2BAA2B,CACzC,UAAe,EACf,SAAc,EACd,GAA0B,EAC1B,KAAa,EACb,MAAc;IAEd,MAAM,OAAO,GAAgB;QAC3B,OAAO,EAAE,gBAAgB;QACzB,MAAM,EAAE,UAAU;QAClB,IAAI,EAAE;YACJ,CAAC,IAAI,CAAC,EAAE,2BAA2B;YACnC,CAAC,QAAQ,CAAC,EAAE,8BAA8B;YAC1C,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,gBAAgB;YAChC,CAAC,SAAS,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,qBAAqB,EAAE,EAAE,GAAG,MAAM,EAAE,CAAC;YACtE,CAAC,QAAQ,EAAE,QAAQ,GAAG,CAAC,UAAU,EAAE,CAAC;YAEpC,sEAAsE;YACtE,CAAC,KAAK,EAAE,YAAY,CAAC;YACrB,CAAC,KAAK,EAAE,eAAe,CAAC;YACxB,CAAC,KAAK,EAAE,eAAe,CAAC;YACxB,CAAC,KAAK,EAAE,sBAAsB,CAAC,EAAE,eAAe;YAChD,CAAC,KAAK,EAAE,wBAAwB,CAAC;YACjC,CAAC,KAAK,EAAE,2BAA2B,CAAC;YACpC,CAAC,KAAK,EAAE,qBAAqB,CAAC;YAC9B,CAAC,KAAK,EAAE,6BAA6B,CAAC;YACtC,CAAC,KAAK,EAAE,aAAa,CAAC;YACtB,CAAC,KAAK,EAAE,gBAAgB,CAAC,EAAE,0BAA0B;YAErD,2FAA2F;YAC3F,CAAC,KAAK,EAAE,eAAe,CAAC;YACxB,WAAW,CAAC,SAAS,CAAC;YACtB,WAAW,CAAC,UAAU,CAAC;SACxB;aACE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,IAAI,CAAC;aACxB,IAAI,EAAE;aACN,GAAG,CAAC,MAAM,CAAC;KACf,CAAC;IAEF,OAAO,OAAO,CAAC;AACjB,CAAC"}

package/build/cogify/gdal/gdal.runner.d.ts

@@ -0,0 +1,42 @@
+import { LogType } from '@basemaps/shared';
+import { EventEmitter } from 'events';
+export interface GdalCommand {
+    /** Output file location */
+    output: URL;
+    /** GDAL command to use */
+    command: 'gdal_create' | 'gdalwarp' | 'gdalbuildvrt' | 'gdal_translate';
+    /** GDAL arguments to use */
+    args: string[];
+}
+export declare class GdalRunner {
+    parser: GdalProgressParser;
+    startTime: number;
+    cmd: GdalCommand;
+    constructor(cmd: GdalCommand);
+    /**
+     * Run a GDAL command
+     *
+     * @param logger logger to use
+     */
+    run(logger?: LogType): Promise<{
+        stdout: string;
+        stderr: string;
+        duration: number;
+    }>;
+}
+/**
+ * Emit a "progress" event every time a "." is recorded in the output
+ */
+export declare class GdalProgressParser extends EventEmitter<{
+    progress: [number];
+}> {
+    waitNewLine: boolean;
+    dotCount: number;
+    byteCount: number;
+    lastEmit: number;
+    debounceMs: number;
+    /** Reset the progress counter */
+    reset(): void;
+    get progress(): number;
+    data(data: Buffer): void;
+}