@developmentseed/deck.gl-geotiff 0.1.0-beta.4

package/dist/ellipsoids.js

@@ -0,0 +1,153 @@
+/**
+ * Vendored and edited from proj4.js.
+ *
+ * In the implementation of metersPerUnit while generated a COG TileMatrixSet,
+ * we need to know the size of the semi major axis in the case the CRS is in
+ * degrees.
+ *
+ * https://github.com/proj4js/proj4js/blob/e90e5fa6872a1ffc40edb161cbeb4bd5e3bd9db5/lib/constants/Ellipsoid.js
+ */
+const ellipsoids = {
+    MERIT: {
+        a: 6378137,
+    },
+    SGS85: {
+        a: 6378136,
+    },
+    GRS80: {
+        a: 6378137,
+    },
+    IAU76: {
+        a: 6378140,
+    },
+    airy: {
+        a: 6377563.396,
+        b: 6356256.91,
+    },
+    APL4: {
+        a: 6378137,
+    },
+    NWL9D: {
+        a: 6378145,
+    },
+    mod_airy: {
+        a: 6377340.189,
+        b: 6356034.446,
+    },
+    andrae: {
+        a: 6377104.43,
+    },
+    aust_SA: {
+        a: 6378160,
+    },
+    GRS67: {
+        a: 6378160,
+    },
+    bessel: {
+        a: 6377397.155,
+    },
+    bess_nam: {
+        a: 6377483.865,
+    },
+    clrk66: {
+        a: 6378206.4,
+        b: 6356583.8,
+    },
+    clrk80: {
+        a: 6378249.145,
+    },
+    clrk80ign: {
+        a: 6378249.2,
+        b: 6356515,
+    },
+    clrk58: {
+        a: 6378293.645208759,
+    },
+    CPM: {
+        a: 6375738.7,
+    },
+    delmbr: {
+        a: 6376428,
+    },
+    engelis: {
+        a: 6378136.05,
+    },
+    evrst30: {
+        a: 6377276.345,
+    },
+    evrst48: {
+        a: 6377304.063,
+    },
+    evrst56: {
+        a: 6377301.243,
+    },
+    evrst69: {
+        a: 6377295.664,
+    },
+    evrstSS: {
+        a: 6377298.556,
+    },
+    fschr60: {
+        a: 6378166,
+    },
+    fschr60m: {
+        a: 6378155,
+    },
+    fschr68: {
+        a: 6378150,
+    },
+    helmert: {
+        a: 6378200,
+    },
+    hough: {
+        a: 6378270,
+    },
+    intl: {
+        a: 6378388,
+    },
+    kaula: {
+        a: 6378163,
+    },
+    lerch: {
+        a: 6378139,
+    },
+    mprts: {
+        a: 6397300,
+    },
+    new_intl: {
+        a: 6378157.5,
+        b: 6356772.2,
+    },
+    plessis: {
+        a: 6376523,
+    },
+    krass: {
+        a: 6378245,
+    },
+    SEasia: {
+        a: 6378155,
+        b: 6356773.3205,
+    },
+    walbeck: {
+        a: 6376896,
+        b: 6355834.8467,
+    },
+    WGS60: {
+        a: 6378165,
+    },
+    WGS66: {
+        a: 6378145,
+    },
+    WGS7: {
+        a: 6378135,
+    },
+    WGS84: {
+        a: 6378137,
+    },
+    sphere: {
+        a: 6370997,
+        b: 6370997,
+    },
+};
+export default ellipsoids;
+//# sourceMappingURL=ellipsoids.js.map

package/dist/ellipsoids.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoids.js","sourceRoot":"","sources":["../src/ellipsoids.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,MAAM,UAAU,GAAG;IACjB,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,IAAI,EAAE;QACJ,CAAC,EAAE,WAAW;QACd,CAAC,EAAE,UAAU;KACd;IACD,IAAI,EAAE;QACJ,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,QAAQ,EAAE;QACR,CAAC,EAAE,WAAW;QACd,CAAC,EAAE,WAAW;KACf;IACD,MAAM,EAAE;QACN,CAAC,EAAE,UAAU;KACd;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,MAAM,EAAE;QACN,CAAC,EAAE,WAAW;KACf;IACD,QAAQ,EAAE;QACR,CAAC,EAAE,WAAW;KACf;IACD,MAAM,EAAE;QACN,CAAC,EAAE,SAAS;QACZ,CAAC,EAAE,SAAS;KACb;IACD,MAAM,EAAE;QACN,CAAC,EAAE,WAAW;KACf;IACD,SAAS,EAAE;QACT,CAAC,EAAE,SAAS;QACZ,CAAC,EAAE,OAAO;KACX;IACD,MAAM,EAAE;QACN,CAAC,EAAE,iBAAiB;KACrB;IACD,GAAG,EAAE;QACH,CAAC,EAAE,SAAS;KACb;IACD,MAAM,EAAE;QACN,CAAC,EAAE,OAAO;KACX;IACD,OAAO,EAAE;QACP,CAAC,EAAE,UAAU;KACd;IACD,OAAO,EAAE;QACP,CAAC,EAAE,WAAW;KACf;IACD,OAAO,EAAE;QACP,CAAC,EAAE,WAAW;KACf;IACD,OAAO,EAAE;QACP,CAAC,EAAE,WAAW;KACf;IACD,OAAO,EAAE;QACP,CAAC,EAAE,WAAW;KACf;IACD,OAAO,EAAE;QACP,CAAC,EAAE,WAAW;KACf;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;KACX;IACD,QAAQ,EAAE;QACR,CAAC,EAAE,OAAO;KACX;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;KACX;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,IAAI,EAAE;QACJ,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,QAAQ,EAAE;QACR,CAAC,EAAE,SAAS;QACZ,CAAC,EAAE,SAAS;KACb;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,MAAM,EAAE;QACN,CAAC,EAAE,OAAO;QACV,CAAC,EAAE,YAAY;KAChB;IACD,OAAO,EAAE;QACP,CAAC,EAAE,OAAO;QACV,CAAC,EAAE,YAAY;KAChB;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,IAAI,EAAE;QACJ,CAAC,EAAE,OAAO;KACX;IACD,KAAK,EAAE;QACL,CAAC,EAAE,OAAO;KACX;IACD,MAAM,EAAE;QACN,CAAC,EAAE,OAAO;QACV,CAAC,EAAE,OAAO;KACX;CACF,CAAC;AAEF,eAAe,UAAU,CAAC"}

package/dist/geotiff-layer.d.ts

@@ -0,0 +1,55 @@
+import type { CompositeLayerProps, UpdateParameters } from "@deck.gl/core";
+import { CompositeLayer } from "@deck.gl/core";
+import { RasterLayer } from "@developmentseed/deck.gl-raster";
+import type { ReprojectionFns } from "@developmentseed/raster-reproject";
+import type { GeoTIFF, Pool } from "geotiff";
+export interface GeoTIFFLayerProps extends CompositeLayerProps {
+    geotiff: GeoTIFF;
+    /**
+     * GeoTIFF.js Pool for decoding image chunks.
+     *
+     * If none is provided, a default Pool will be created and shared between all
+     * COGLayer and GeoTIFFLayer instances.
+     */
+    pool?: Pool;
+    /**
+     * Maximum reprojection error in pixels for mesh refinement.
+     * Lower values create denser meshes with higher accuracy.
+     * @default 0.125
+     */
+    maxError?: number;
+    /**
+     * Enable debug visualization showing the triangulation mesh
+     * @default false
+     */
+    debug?: boolean;
+    /**
+     * Opacity of the debug mesh overlay (0-1)
+     * @default 0.5
+     */
+    debugOpacity?: number;
+}
+/**
+ * GeoTIFFLayer renders a GeoTIFF file from an arbitrary projection.
+ *
+ * The GeoTIFFLayer differs from the COGLayer in that it doesn't assume any
+ * internal tiling. Rather, it fetches the entire full-resolution image and
+ * displays it directly.
+ */
+export declare class GeoTIFFLayer extends CompositeLayer<GeoTIFFLayerProps> {
+    static layerName: string;
+    static defaultProps: {
+        maxError: number;
+    };
+    state: {
+        reprojectionFns?: ReprojectionFns;
+        imageData?: ImageData;
+        height?: number;
+        width?: number;
+    };
+    initializeState(): void;
+    updateState(params: UpdateParameters<this>): void;
+    _parseGeoTIFF(): Promise<void>;
+    renderLayers(): RasterLayer | null;
+}
+//# sourceMappingURL=geotiff-layer.d.ts.map

package/dist/geotiff-layer.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff-layer.d.ts","sourceRoot":"","sources":["../src/geotiff-layer.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,mBAAmB,EAAE,gBAAgB,EAAE,MAAM,eAAe,CAAC;AAC3E,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAC/C,OAAO,EAAE,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAC9D,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,mCAAmC,CAAC;AACzE,OAAO,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,SAAS,CAAC;AAM7C,MAAM,WAAW,iBAAkB,SAAQ,mBAAmB;IAC5D,OAAO,EAAE,OAAO,CAAC;IAEjB;;;;;OAKG;IACH,IAAI,CAAC,EAAE,IAAI,CAAC;IAEZ;;;;OAIG;IACH,QAAQ,CAAC,EAAE,MAAM,CAAC;IAElB;;;OAGG;IACH,KAAK,CAAC,EAAE,OAAO,CAAC;IAEhB;;;OAGG;IACH,YAAY,CAAC,EAAE,MAAM,CAAC;CACvB;AAMD;;;;;;GAMG;AACH,qBAAa,YAAa,SAAQ,cAAc,CAAC,iBAAiB,CAAC;IACjE,OAAgB,SAAS,SAAkB;IAC3C,OAAgB,YAAY;;MAAgB;IAEpC,KAAK,EAAE;QACb,eAAe,CAAC,EAAE,eAAe,CAAC;QAClC,SAAS,CAAC,EAAE,SAAS,CAAC;QACtB,MAAM,CAAC,EAAE,MAAM,CAAC;QAChB,KAAK,CAAC,EAAE,MAAM,CAAC;KAChB,CAAC;IAEO,eAAe,IAAI,IAAI;IAIvB,WAAW,CAAC,MAAM,EAAE,gBAAgB,CAAC,IAAI,CAAC;IAe7C,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC;IAiBpC,YAAY;CAsBb"}

package/dist/geotiff-layer.js

@@ -0,0 +1,64 @@
+import { CompositeLayer } from "@deck.gl/core";
+import { RasterLayer } from "@developmentseed/deck.gl-raster";
+import { extractGeotiffReprojectors } from "./geotiff-reprojection.js";
+import { defaultPool, loadRgbImage } from "./geotiff.js";
+const DEFAULT_MAX_ERROR = 0.125;
+const defaultProps = {
+    maxError: DEFAULT_MAX_ERROR,
+};
+/**
+ * GeoTIFFLayer renders a GeoTIFF file from an arbitrary projection.
+ *
+ * The GeoTIFFLayer differs from the COGLayer in that it doesn't assume any
+ * internal tiling. Rather, it fetches the entire full-resolution image and
+ * displays it directly.
+ */
+export class GeoTIFFLayer extends CompositeLayer {
+    static layerName = "GeoTIFFLayer";
+    static defaultProps = defaultProps;
+    initializeState() {
+        this.setState({});
+    }
+    updateState(params) {
+        super.updateState(params);
+        const { props, oldProps, changeFlags } = params;
+        const needsUpdate = Boolean(changeFlags.dataChanged) ||
+            props.geotiff !== oldProps.geotiff ||
+            props.maxError !== oldProps.maxError;
+        if (needsUpdate) {
+            this._parseGeoTIFF();
+        }
+    }
+    async _parseGeoTIFF() {
+        const { geotiff } = this.props;
+        const reprojectionFns = await extractGeotiffReprojectors(geotiff);
+        const image = await geotiff.getImage();
+        const { imageData, height, width } = await loadRgbImage(image, {
+            pool: this.props.pool || defaultPool(),
+        });
+        this.setState({
+            reprojectionFns,
+            imageData,
+            height,
+            width,
+        });
+    }
+    renderLayers() {
+        const { reprojectionFns, imageData, height, width } = this.state;
+        if (!reprojectionFns || !imageData || !height || !width) {
+            return null;
+        }
+        const { maxError, debug, debugOpacity } = this.props;
+        return new RasterLayer(this.getSubLayerProps({
+            id: "raster",
+            width,
+            height,
+            reprojectionFns,
+            maxError,
+            texture: imageData,
+            debug,
+            debugOpacity,
+        }));
+    }
+}
+//# sourceMappingURL=geotiff-layer.js.map

package/dist/geotiff-layer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff-layer.js","sourceRoot":"","sources":["../src/geotiff-layer.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAC/C,OAAO,EAAE,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAG9D,OAAO,EAAE,0BAA0B,EAAE,MAAM,2BAA2B,CAAC;AACvE,OAAO,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,cAAc,CAAC;AAEzD,MAAM,iBAAiB,GAAG,KAAK,CAAC;AAiChC,MAAM,YAAY,GAAG;IACnB,QAAQ,EAAE,iBAAiB;CAC5B,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,OAAO,YAAa,SAAQ,cAAiC;IACjE,MAAM,CAAU,SAAS,GAAG,cAAc,CAAC;IAC3C,MAAM,CAAU,YAAY,GAAG,YAAY,CAAC;IASnC,eAAe;QACtB,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IACpB,CAAC;IAEQ,WAAW,CAAC,MAA8B;QACjD,KAAK,CAAC,WAAW,CAAC,MAAM,CAAC,CAAC;QAE1B,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,MAAM,CAAC;QAEhD,MAAM,WAAW,GACf,OAAO,CAAC,WAAW,CAAC,WAAW,CAAC;YAChC,KAAK,CAAC,OAAO,KAAK,QAAQ,CAAC,OAAO;YAClC,KAAK,CAAC,QAAQ,KAAK,QAAQ,CAAC,QAAQ,CAAC;QAEvC,IAAI,WAAW,EAAE,CAAC;YAChB,IAAI,CAAC,aAAa,EAAE,CAAC;QACvB,CAAC;IACH,CAAC;IAED,KAAK,CAAC,aAAa;QACjB,MAAM,EAAE,OAAO,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QAE/B,MAAM,eAAe,GAAG,MAAM,0BAA0B,CAAC,OAAO,CAAC,CAAC;QAClE,MAAM,KAAK,GAAG,MAAM,OAAO,CAAC,QAAQ,EAAE,CAAC;QACvC,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,MAAM,YAAY,CAAC,KAAK,EAAE;YAC7D,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,IAAI,WAAW,EAAE;SACvC,CAAC,CAAC;QAEH,IAAI,CAAC,QAAQ,CAAC;YACZ,eAAe;YACf,SAAS;YACT,MAAM;YACN,KAAK;SACN,CAAC,CAAC;IACL,CAAC;IAED,YAAY;QACV,MAAM,EAAE,eAAe,EAAE,SAAS,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QAEjE,IAAI,CAAC,eAAe,IAAI,CAAC,SAAS,IAAI,CAAC,MAAM,IAAI,CAAC,KAAK,EAAE,CAAC;YACxD,OAAO,IAAI,CAAC;QACd,CAAC;QAED,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,YAAY,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QAErD,OAAO,IAAI,WAAW,CACpB,IAAI,CAAC,gBAAgB,CAAC;YACpB,EAAE,EAAE,QAAQ;YACZ,KAAK;YACL,MAAM;YACN,eAAe;YACf,QAAQ;YACR,OAAO,EAAE,SAAS;YAClB,KAAK;YACL,YAAY;SACb,CAAC,CACH,CAAC;IACJ,CAAC"}

package/dist/geotiff-reprojection.d.ts

@@ -0,0 +1,31 @@
+import type { GeoTIFF, GeoTIFFImage } from "geotiff";
+import type { ReprojectionFns } from "@developmentseed/raster-reproject";
+import type { PROJJSONDefinition } from "proj4/dist/lib/core";
+import type Projection from "proj4/dist/lib/Proj";
+export declare function extractGeotiffReprojectors(tiff: GeoTIFF, outputCrs?: string | PROJJSONDefinition | Projection): Promise<ReprojectionFns>;
+export declare function fromGeoTransform(geotransform: [number, number, number, number, number, number]): {
+    pixelToInputCRS: (x: number, y: number) => [number, number];
+    inputCRSToPixel: (x: number, y: number) => [number, number];
+};
+/**
+ * Get the Projection of a GeoTIFF
+ *
+ * The first `image` must be passed in, as only the top-level IFD contains geo
+ * keys.
+ */
+export declare function getGeoTIFFProjection(image: GeoTIFFImage): Promise<PROJJSONDefinition | null>;
+/**
+ * Extract affine geotransform from a GeoTIFF image.
+ *
+ * The first `image` must be passed in, as only the top-level IFD contains geo
+ * keys.
+ *
+ * Returns a 6-element array in Python `affine` package ordering:
+ * [a, b, c, d, e, f] where:
+ * - x_geo = a * col + b * row + c
+ * - y_geo = d * col + e * row + f
+ *
+ * This is NOT GDAL ordering, which is [c, a, b, f, d, e].
+ */
+export declare function extractGeotransform(image: GeoTIFFImage): [number, number, number, number, number, number];
+//# sourceMappingURL=geotiff-reprojection.d.ts.map

package/dist/geotiff-reprojection.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff-reprojection.d.ts","sourceRoot":"","sources":["../src/geotiff-reprojection.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,OAAO,EAAE,YAAY,EAAE,MAAM,SAAS,CAAC;AACrD,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,mCAAmC,CAAC;AAMzE,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AAC9D,OAAO,KAAK,UAAU,MAAM,qBAAqB,CAAC;AA+ElD,wBAAsB,0BAA0B,CAC9C,IAAI,EAAE,OAAO,EACb,SAAS,GAAE,MAAM,GAAG,kBAAkB,GAAG,UAAmB,GAC3D,OAAO,CAAC,eAAe,CAAC,CAyB1B;AAED,wBAAgB,gBAAgB,CAC9B,YAAY,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,GAC7D;IACD,eAAe,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC5D,eAAe,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAC7D,CAOA;AAED;;;;;GAKG;AACH,wBAAsB,oBAAoB,CACxC,KAAK,EAAE,YAAY,GAClB,OAAO,CAAC,kBAAkB,GAAG,IAAI,CAAC,CAOpC;AAiBD;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CACjC,KAAK,EAAE,YAAY,GAClB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAiClD"}

package/dist/geotiff-reprojection.js

@@ -0,0 +1,172 @@
+/* eslint-env browser */
+import { applyAffine, invertGeoTransform, } from "@developmentseed/raster-reproject";
+import proj4 from "proj4";
+const OGC_84 = {
+    $schema: "https://proj.org/schemas/v0.7/projjson.schema.json",
+    type: "GeographicCRS",
+    name: "WGS 84 (CRS84)",
+    datum_ensemble: {
+        name: "World Geodetic System 1984 ensemble",
+        members: [
+            {
+                name: "World Geodetic System 1984 (Transit)",
+                id: { authority: "EPSG", code: 1166 },
+            },
+            {
+                name: "World Geodetic System 1984 (G730)",
+                id: { authority: "EPSG", code: 1152 },
+            },
+            {
+                name: "World Geodetic System 1984 (G873)",
+                id: { authority: "EPSG", code: 1153 },
+            },
+            {
+                name: "World Geodetic System 1984 (G1150)",
+                id: { authority: "EPSG", code: 1154 },
+            },
+            {
+                name: "World Geodetic System 1984 (G1674)",
+                id: { authority: "EPSG", code: 1155 },
+            },
+            {
+                name: "World Geodetic System 1984 (G1762)",
+                id: { authority: "EPSG", code: 1156 },
+            },
+            {
+                name: "World Geodetic System 1984 (G2139)",
+                id: { authority: "EPSG", code: 1309 },
+            },
+        ],
+        ellipsoid: {
+            name: "WGS 84",
+            semi_major_axis: 6378137,
+            inverse_flattening: 298.257223563,
+        },
+        accuracy: "2.0",
+        id: { authority: "EPSG", code: 6326 },
+    },
+    coordinate_system: {
+        subtype: "ellipsoidal",
+        axis: [
+            {
+                name: "Geodetic longitude",
+                abbreviation: "Lon",
+                direction: "east",
+                unit: "degree",
+            },
+            {
+                name: "Geodetic latitude",
+                abbreviation: "Lat",
+                direction: "north",
+                unit: "degree",
+            },
+        ],
+    },
+    scope: "Not known.",
+    area: "World.",
+    bbox: {
+        south_latitude: -90,
+        west_longitude: -180,
+        north_latitude: 90,
+        east_longitude: 180,
+    },
+    // @ts-expect-error - proj4 types are incomplete
+    id: { authority: "OGC", code: "CRS84" },
+};
+// Derived from existing work here:
+// https://github.com/developmentseed/lonboard/blob/35a1f3d691604ad9e083bf10a4bfde4158171486/src/cog-tileset/claude-tileset-2d-improved.ts#L141
+//
+// TODO: return a RasterReprojector instance, given the IFD and tile of interest?
+export async function extractGeotiffReprojectors(tiff, outputCrs = OGC_84) {
+    const image = await tiff.getImage();
+    // Extract geotransform from full-resolution image
+    // Only the top-level IFD has geo keys, so we'll derive overviews from this
+    const baseGeotransform = extractGeotransform(image);
+    const sourceProjection = await getGeoTIFFProjection(image);
+    if (sourceProjection === null) {
+        throw new Error("Could not determine source projection from GeoTIFF geo keys");
+    }
+    const converter = proj4(sourceProjection, outputCrs);
+    const { pixelToInputCRS, inputCRSToPixel } = fromGeoTransform(baseGeotransform);
+    return {
+        pixelToInputCRS,
+        inputCRSToPixel,
+        forwardReproject: (x, y) => converter.forward([x, y], false),
+        inverseReproject: (x, y) => converter.inverse([x, y], false),
+    };
+}
+export function fromGeoTransform(geotransform) {
+    const inverseGeotransform = invertGeoTransform(geotransform);
+    return {
+        pixelToInputCRS: (x, y) => applyAffine(x, y, geotransform),
+        inputCRSToPixel: (x, y) => applyAffine(x, y, inverseGeotransform),
+    };
+}
+/**
+ * Get the Projection of a GeoTIFF
+ *
+ * The first `image` must be passed in, as only the top-level IFD contains geo
+ * keys.
+ */
+export async function getGeoTIFFProjection(image) {
+    const geoKeys = image.getGeoKeys();
+    const projectionCode = geoKeys.ProjectedCSTypeGeoKey || geoKeys.GeographicTypeGeoKey || null;
+    const sourceProjection = await getProjjson(projectionCode);
+    return sourceProjection;
+}
+/** Query epsg.io for the PROJJSON corresponding to the given EPSG code. */
+async function getProjjson(projectionCode) {
+    if (projectionCode === null) {
+        return null;
+    }
+    const url = `https://epsg.io/${projectionCode}.json`;
+    const response = await fetch(url);
+    if (!response.ok) {
+        throw new Error(`Failed to fetch projection data from ${url}`);
+    }
+    const data = await response.json();
+    return data;
+}
+/**
+ * Extract affine geotransform from a GeoTIFF image.
+ *
+ * The first `image` must be passed in, as only the top-level IFD contains geo
+ * keys.
+ *
+ * Returns a 6-element array in Python `affine` package ordering:
+ * [a, b, c, d, e, f] where:
+ * - x_geo = a * col + b * row + c
+ * - y_geo = d * col + e * row + f
+ *
+ * This is NOT GDAL ordering, which is [c, a, b, f, d, e].
+ */
+export function extractGeotransform(image) {
+    const origin = image.getOrigin();
+    const resolution = image.getResolution();
+    // origin: [x, y, z]
+    // resolution: [x_res, y_res, z_res]
+    // Check for rotation/skew in the file directory
+    const fileDirectory = image.getFileDirectory();
+    const modelTransformation = fileDirectory.ModelTransformation;
+    let b = 0; // row rotation
+    let d = 0; // column rotation
+    if (modelTransformation && modelTransformation.length >= 16) {
+        // ModelTransformation is a 4x4 matrix in row-major order
+        // [0  1  2  3 ]   [a  b  0  c]
+        // [4  5  6  7 ] = [d  e  0  f]
+        // [8  9  10 11]   [0  0  1  0]
+        // [12 13 14 15]   [0  0  0  1]
+        b = modelTransformation[1];
+        d = modelTransformation[4];
+    }
+    // Return in affine package ordering: [a, b, c, d, e, f]
+    return [
+        resolution[0], // a: pixel width
+        b, // b: row rotation
+        origin[0], // c: x origin
+        d, // d: column rotation
+        resolution[1], // e: pixel height (often negative)
+        origin[1],
+    ];
+}
+//# sourceMappingURL=geotiff-reprojection.js.map

package/dist/geotiff-reprojection.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff-reprojection.js","sourceRoot":"","sources":["../src/geotiff-reprojection.ts"],"names":[],"mappings":"AAAA,wBAAwB;AAIxB,OAAO,EACL,WAAW,EACX,kBAAkB,GACnB,MAAM,mCAAmC,CAAC;AAC3C,OAAO,KAAK,MAAM,OAAO,CAAC;AAI1B,MAAM,MAAM,GAAuB;IACjC,OAAO,EAAE,oDAAoD;IAC7D,IAAI,EAAE,eAAe;IACrB,IAAI,EAAE,gBAAgB;IACtB,cAAc,EAAE;QACd,IAAI,EAAE,qCAAqC;QAC3C,OAAO,EAAE;YACP;gBACE,IAAI,EAAE,sCAAsC;gBAC5C,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,mCAAmC;gBACzC,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,mCAAmC;gBACzC,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,oCAAoC;gBAC1C,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,oCAAoC;gBAC1C,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,oCAAoC;gBAC1C,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;YACD;gBACE,IAAI,EAAE,oCAAoC;gBAC1C,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;aACtC;SACF;QACD,SAAS,EAAE;YACT,IAAI,EAAE,QAAQ;YACd,eAAe,EAAE,OAAO;YACxB,kBAAkB,EAAE,aAAa;SAClC;QACD,QAAQ,EAAE,KAAK;QACf,EAAE,EAAE,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE;KACtC;IACD,iBAAiB,EAAE;QACjB,OAAO,EAAE,aAAa;QACtB,IAAI,EAAE;YACJ;gBACE,IAAI,EAAE,oBAAoB;gBAC1B,YAAY,EAAE,KAAK;gBACnB,SAAS,EAAE,MAAM;gBACjB,IAAI,EAAE,QAAQ;aACf;YACD;gBACE,IAAI,EAAE,mBAAmB;gBACzB,YAAY,EAAE,KAAK;gBACnB,SAAS,EAAE,OAAO;gBAClB,IAAI,EAAE,QAAQ;aACf;SACF;KACF;IACD,KAAK,EAAE,YAAY;IACnB,IAAI,EAAE,QAAQ;IACd,IAAI,EAAE;QACJ,cAAc,EAAE,CAAC,EAAE;QACnB,cAAc,EAAE,CAAC,GAAG;QACpB,cAAc,EAAE,EAAE;QAClB,cAAc,EAAE,GAAG;KACpB;IACD,gDAAgD;IAChD,EAAE,EAAE,EAAE,SAAS,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,EAAE;CACxC,CAAC;AAEF,mCAAmC;AACnC,+IAA+I;AAC/I,EAAE;AACF,iFAAiF;AACjF,MAAM,CAAC,KAAK,UAAU,0BAA0B,CAC9C,IAAa,EACb,YAAsD,MAAM;IAE5D,MAAM,KAAK,GAAG,MAAM,IAAI,CAAC,QAAQ,EAAE,CAAC;IAEpC,kDAAkD;IAClD,2EAA2E;IAC3E,MAAM,gBAAgB,GAAG,mBAAmB,CAAC,KAAK,CAAC,CAAC;IAEpD,MAAM,gBAAgB,GAAG,MAAM,oBAAoB,CAAC,KAAK,CAAC,CAAC;IAC3D,IAAI,gBAAgB,KAAK,IAAI,EAAE,CAAC;QAC9B,MAAM,IAAI,KAAK,CACb,6DAA6D,CAC9D,CAAC;IACJ,CAAC;IACD,MAAM,SAAS,GAAG,KAAK,CAAC,gBAAgB,EAAE,SAAS,CAAC,CAAC;IACrD,MAAM,EAAE,eAAe,EAAE,eAAe,EAAE,GACxC,gBAAgB,CAAC,gBAAgB,CAAC,CAAC;IAErC,OAAO;QACL,eAAe;QACf,eAAe;QACf,gBAAgB,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CACzC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC;QAClC,gBAAgB,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CACzC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,KAAK,CAAC;KACnC,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,gBAAgB,CAC9B,YAA8D;IAK9D,MAAM,mBAAmB,GAAG,kBAAkB,CAAC,YAAY,CAAC,CAAC;IAC7D,OAAO;QACL,eAAe,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC;QAC1E,eAAe,EAAE,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CACxC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,mBAAmB,CAAC;KACzC,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,MAAM,CAAC,KAAK,UAAU,oBAAoB,CACxC,KAAmB;IAEnB,MAAM,OAAO,GAAG,KAAK,CAAC,UAAU,EAAE,CAAC;IACnC,MAAM,cAAc,GAClB,OAAO,CAAC,qBAAqB,IAAI,OAAO,CAAC,oBAAoB,IAAI,IAAI,CAAC;IAExE,MAAM,gBAAgB,GAAG,MAAM,WAAW,CAAC,cAAc,CAAC,CAAC;IAC3D,OAAO,gBAAgB,CAAC;AAC1B,CAAC;AAED,2EAA2E;AAC3E,KAAK,UAAU,WAAW,CAAC,cAA6B;IACtD,IAAI,cAAc,KAAK,IAAI,EAAE,CAAC;QAC5B,OAAO,IAAI,CAAC;IACd,CAAC;IAED,MAAM,GAAG,GAAG,mBAAmB,cAAc,OAAO,CAAC;IACrD,MAAM,QAAQ,GAAG,MAAM,KAAK,CAAC,GAAG,CAAC,CAAC;IAClC,IAAI,CAAC,QAAQ,CAAC,EAAE,EAAE,CAAC;QACjB,MAAM,IAAI,KAAK,CAAC,wCAAwC,GAAG,EAAE,CAAC,CAAC;IACjE,CAAC;IACD,MAAM,IAAI,GAAG,MAAM,QAAQ,CAAC,IAAI,EAAE,CAAC;IACnC,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CACjC,KAAmB;IAEnB,MAAM,MAAM,GAAG,KAAK,CAAC,SAAS,EAAE,CAAC;IACjC,MAAM,UAAU,GAAG,KAAK,CAAC,aAAa,EAAE,CAAC;IAEzC,oBAAoB;IACpB,oCAAoC;IAEpC,gDAAgD;IAChD,MAAM,aAAa,GAAG,KAAK,CAAC,gBAAgB,EAAE,CAAC;IAC/C,MAAM,mBAAmB,GAAG,aAAa,CAAC,mBAAmB,CAAC;IAE9D,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,eAAe;IAC1B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,kBAAkB;IAE7B,IAAI,mBAAmB,IAAI,mBAAmB,CAAC,MAAM,IAAI,EAAE,EAAE,CAAC;QAC5D,yDAAyD;QACzD,+BAA+B;QAC/B,+BAA+B;QAC/B,+BAA+B;QAC/B,+BAA+B;QAC/B,CAAC,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;QAC3B,CAAC,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAC7B,CAAC;IAED,wDAAwD;IACxD,OAAO;QACL,UAAU,CAAC,CAAC,CAAE,EAAE,iBAAiB;QACjC,CAAC,EAAE,kBAAkB;QACrB,MAAM,CAAC,CAAC,CAAE,EAAE,cAAc;QAC1B,CAAC,EAAE,qBAAqB;QACxB,UAAU,CAAC,CAAC,CAAE,EAAE,mCAAmC;QACnD,MAAM,CAAC,CAAC,CAAE;KACX,CAAC;AACJ,CAAC"}

package/dist/geotiff.d.ts

@@ -0,0 +1,31 @@
+import type { GeoTIFFImage } from "geotiff";
+import { Pool } from "geotiff";
+/**
+ * Options that may be passed when reading image data from geotiff.js
+ */
+type ReadRasterOptions = {
+    /** the subset to read data from in pixels. */
+    window?: [number, number, number, number];
+    /** The optional decoder pool to use. */
+    pool?: Pool;
+    /** An AbortSignal that may be signalled if the request is to be aborted */
+    signal?: AbortSignal;
+};
+/**
+ * Retrieve the default geotiff.js decoder Pool.
+ *
+ * If a Pool has not yet been created, it will be created on first call.
+ *
+ * The Pool will be shared between all COGLayer and GeoTIFFLayer instances.
+ */
+export declare function defaultPool(): Pool;
+/**
+ * Load an RGBA image from a GeoTIFFImage.
+ */
+export declare function loadRgbImage(image: GeoTIFFImage, options?: ReadRasterOptions): Promise<{
+    imageData: ImageData;
+    height: number;
+    width: number;
+}>;
+export {};
+//# sourceMappingURL=geotiff.d.ts.map

package/dist/geotiff.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff.d.ts","sourceRoot":"","sources":["../src/geotiff.ts"],"names":[],"mappings":"AAEA,OAAO,KAAK,EAAE,YAAY,EAA4B,MAAM,SAAS,CAAC;AACtE,OAAO,EAAE,IAAI,EAAE,MAAM,SAAS,CAAC;AAE/B;;GAEG;AACH,KAAK,iBAAiB,GAAG;IACvB,8CAA8C;IAC9C,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE1C,wCAAwC;IACxC,IAAI,CAAC,EAAE,IAAI,CAAC;IAEZ,2EAA2E;IAC3E,MAAM,CAAC,EAAE,WAAW,CAAC;CACtB,CAAC;AASF;;;;;;GAMG;AACH,wBAAgB,WAAW,IAAI,IAAI,CAMlC;AAED;;GAEG;AACH,wBAAsB,YAAY,CAChC,KAAK,EAAE,YAAY,EACnB,OAAO,CAAC,EAAE,iBAAiB,GAC1B,OAAO,CAAC;IAAE,SAAS,EAAE,SAAS,CAAC;IAAC,MAAM,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC,CAmBlE"}

package/dist/geotiff.js

@@ -0,0 +1,70 @@
+// Utilities for interacting with geotiff.js.
+import { Pool } from "geotiff";
+/**
+ * A default geotiff.js decoder pool instance.
+ *
+ * It will be created on first call of `defaultPool`.
+ */
+let DEFAULT_POOL = null;
+/**
+ * Retrieve the default geotiff.js decoder Pool.
+ *
+ * If a Pool has not yet been created, it will be created on first call.
+ *
+ * The Pool will be shared between all COGLayer and GeoTIFFLayer instances.
+ */
+export function defaultPool() {
+    if (DEFAULT_POOL === null) {
+        DEFAULT_POOL = new Pool();
+    }
+    return DEFAULT_POOL;
+}
+/**
+ * Load an RGBA image from a GeoTIFFImage.
+ */
+export async function loadRgbImage(image, options) {
+    const mergedOptions = {
+        ...options,
+        interleave: true,
+        enableAlpha: true,
+    };
+    // Since we set interleave: true, the result is a single array with all
+    // samples, so we cast to TypedArrayWithDimensions
+    // https://github.com/geotiffjs/geotiff.js/issues/486
+    const rgbImage = (await image.readRGB(mergedOptions));
+    const imageData = addAlphaChannel(rgbImage);
+    return {
+        imageData,
+        height: rgbImage.height,
+        width: rgbImage.width,
+    };
+}
+/**
+ * Add an alpha channel to an RGB image array.
+ *
+ * Only supports input arrays with 3 (RGB) or 4 (RGBA) channels. If the input is
+ * already RGBA, it is returned unchanged.
+ */
+function addAlphaChannel(rgbImage) {
+    const { height, width } = rgbImage;
+    if (rgbImage.length === height * width * 4) {
+        // Already has alpha channel
+        return new ImageData(new Uint8ClampedArray(rgbImage), width, height);
+    }
+    else if (rgbImage.length === height * width * 3) {
+        // Need to add alpha channel
+        const rgbaLength = (rgbImage.length / 3) * 4;
+        const rgbaArray = new Uint8ClampedArray(rgbaLength);
+        for (let i = 0; i < rgbImage.length / 3; ++i) {
+            rgbaArray[i * 4] = rgbImage[i * 3];
+            rgbaArray[i * 4 + 1] = rgbImage[i * 3 + 1];
+            rgbaArray[i * 4 + 2] = rgbImage[i * 3 + 2];
+            rgbaArray[i * 4 + 3] = 255;
+        }
+        return new ImageData(rgbaArray, width, height);
+    }
+    else {
+        throw new Error(`Unexpected number of channels in raster data: ${rgbImage.length / (height * width)}`);
+    }
+}
+//# sourceMappingURL=geotiff.js.map

package/dist/geotiff.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geotiff.js","sourceRoot":"","sources":["../src/geotiff.ts"],"names":[],"mappings":"AAAA,6CAA6C;AAG7C,OAAO,EAAE,IAAI,EAAE,MAAM,SAAS,CAAC;AAgB/B;;;;GAIG;AACH,IAAI,YAAY,GAAgB,IAAI,CAAC;AAErC;;;;;;GAMG;AACH,MAAM,UAAU,WAAW;IACzB,IAAI,YAAY,KAAK,IAAI,EAAE,CAAC;QAC1B,YAAY,GAAG,IAAI,IAAI,EAAE,CAAC;IAC5B,CAAC;IAED,OAAO,YAAY,CAAC;AACtB,CAAC;AAED;;GAEG;AACH,MAAM,CAAC,KAAK,UAAU,YAAY,CAChC,KAAmB,EACnB,OAA2B;IAE3B,MAAM,aAAa,GAAG;QACpB,GAAG,OAAO;QACV,UAAU,EAAE,IAAI;QAChB,WAAW,EAAE,IAAI;KAClB,CAAC;IACF,uEAAuE;IACvE,kDAAkD;IAClD,qDAAqD;IACrD,MAAM,QAAQ,GAAG,CAAC,MAAM,KAAK,CAAC,OAAO,CACnC,aAAa,CACd,CAA6B,CAAC;IAC/B,MAAM,SAAS,GAAG,eAAe,CAAC,QAAQ,CAAC,CAAC;IAE5C,OAAO;QACL,SAAS;QACT,MAAM,EAAE,QAAQ,CAAC,MAAM;QACvB,KAAK,EAAE,QAAQ,CAAC,KAAK;KACtB,CAAC;AACJ,CAAC;AAED;;;;;GAKG;AACH,SAAS,eAAe,CAAC,QAAkC;IACzD,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,QAAQ,CAAC;IAEnC,IAAI,QAAQ,CAAC,MAAM,KAAK,MAAM,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC;QAC3C,4BAA4B;QAC5B,OAAO,IAAI,SAAS,CAAC,IAAI,iBAAiB,CAAC,QAAQ,CAAC,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACvE,CAAC;SAAM,IAAI,QAAQ,CAAC,MAAM,KAAK,MAAM,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC;QAClD,4BAA4B;QAE5B,MAAM,UAAU,GAAG,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC7C,MAAM,SAAS,GAAG,IAAI,iBAAiB,CAAC,UAAU,CAAC,CAAC;QACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;YAC7C,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;YACpC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;YAC5C,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;YAC5C,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7B,CAAC;QAED,OAAO,IAAI,SAAS,CAAC,SAAS,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACjD,CAAC;SAAM,CAAC;QACN,MAAM,IAAI,KAAK,CACb,iDAAiD,QAAQ,CAAC,MAAM,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,EAAE,CACtF,CAAC;IACJ,CAAC;AACH,CAAC"}