feat: Overhaul to use our geotiff package & generic TileMatrixSet support (#208)

This is a really big PR, but unfortunately it really had to be a big PR
because we effectively had to change our geotiff reader _and_ our
TileMatrixSet traversal at the same time.

### Change list

- Fully removes `geotiff.js` as a dependency. We now only use
`@developmentseed/geotiff` and the underlying `@cogeotiff/core` for
reading data.
- The `TileMatrixSetTileset` now accepts an _arbitrary_,
_spec-compliant_, OGC TileMatrixSet instance. This will make it much
more interoperable. Closes
#216. This
should make it easier for Lonboard to generate a TMS from async-geotiff,
pass the TMS definition to JS, and then have it auto-fetch tiles from
Python according to that tiling system.

- Removes `GeoTIFFLayer` from our public API for now, since our geotiff
implementation doesn't support fetching the entire full-resolution image
and stitching tiles together. (And it's probably a candidate to fully
remove in the future in favor of `COGLayer`)
- Define new projection resolver `epsgResolver`. Our default
implementation caches all fetched projections in our
`PROJECTION_REGISTRY` (closes
#185). But users
can easily override this behavior by passing a custom `epsgResolver`
- Removes multiple copies of affine handling, ensuring we put all affine
transformation handling into `@developmentseed/affine`
- Update `cog-basic` example for latest changes. (Other examples haven't
been updated yet)

`@cogeotiff/core` is _very_ pleasant to use. It's nice to use Typescript
Enums instead of hard-coding integer matching.

### Todo

- [ ] Ensure request cancellation (needs
blacha/cogeotiff#1400 upstream). Should make the
deck.gl visualization quite a bit faster.
- [ ] Move image decoding onto a thread pool in workers.

Closes #216

Author: kylebarron

examples/cog-basic/package.json

@@ -14,12 +14,11 @@
     "@deck.gl/layers": "^9.2.7",
     "@deck.gl/mapbox": "^9.2.7",
     "@deck.gl/mesh-layers": "^9.2.7",
+    "@developmentseed/geotiff": "workspace:^",
     "@developmentseed/deck.gl-geotiff": "workspace:^",
     "@developmentseed/deck.gl-raster": "workspace:^",
     "@luma.gl/core": "9.2.6",
     "@luma.gl/shadertools": "9.2.6",
-    "geotiff": "^2.1.3",
-    "geotiff-geokeys-to-proj4": "^2024.4.13",
     "maplibre-gl": "^5.17.0",
     "proj4": "^2.20.2",
     "react": "^19.2.4",

examples/cog-basic/src/App.tsx

@@ -1,7 +1,6 @@
 import type { DeckProps } from "@deck.gl/core";
 import { MapboxOverlay } from "@deck.gl/mapbox";
-import { COGLayer, proj } from "@developmentseed/deck.gl-geotiff";
-import { toProj4 } from "geotiff-geokeys-to-proj4";
+import { COGLayer } from "@developmentseed/deck.gl-geotiff";
 import "maplibre-gl/dist/maplibre-gl.css";
 import { useRef, useState } from "react";
 import type { MapRef } from "react-map-gl/maplibre";
@@ -13,26 +12,14 @@ function DeckGLOverlay(props: DeckProps) {
   return null;
 }
 
-async function geoKeysParser(
-  geoKeys: Record<string, any>,
-): Promise<proj.ProjectionInfo> {
-  const projDefinition = toProj4(geoKeys as any);
-
-  return {
-    def: projDefinition.proj4,
-    parsed: proj.parseCrs(projDefinition.proj4),
-    coordinatesUnits: projDefinition.coordinatesUnits as proj.SupportedCrsUnit,
-  };
-}
-
-// const COG_URL =
-//   "https://nz-imagery.s3-ap-southeast-2.amazonaws.com/new-zealand/new-zealand_2024-2025_10m/rgb/2193/CC11.tiff";
+const COG_URL =
+  "https://nz-imagery.s3-ap-southeast-2.amazonaws.com/new-zealand/new-zealand_2024-2025_10m/rgb/2193/CC11.tiff";
 
 // const COG_URL =
 //   "https://ds-wheels.s3.us-east-1.amazonaws.com/m_4007307_sw_18_060_20220803.tif";
 
-const COG_URL =
-  "https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/18/T/WL/2026/1/S2B_18TWL_20260101_0_L2A/TCI.tif";
+// const COG_URL =
+//   "https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/18/T/WL/2026/1/S2B_18TWL_20260101_0_L2A/TCI.tif";
 // const COG_URL =
 //   "https://ds-wheels.s3.us-east-1.amazonaws.com/Annual_NLCD_LndCov_2023_CU_C1V0.tif";
 
@@ -62,7 +49,6 @@ export default function App() {
     geotiff: COG_URL,
     debug,
     debugOpacity,
-    geoKeysParser,
     onGeoTIFFLoad: (tiff, options) => {
       (window as any).tiff = tiff;
       const { west, south, east, north } = options.geographicBounds;

packages/affine/src/affine.ts

@@ -81,7 +81,11 @@ export function invert([sa, sb, sc, sd, se, sf]: Affine): Affine {
   const rd = -sd * idet;
   const re = sa * idet;
 
-  return [ra, rb, -sc * ra - sf * rb, rd, re, -sc * rd - sf * re];
+  // biome-ignore format: array
+  return [
+    ra, rb, -sc * ra - sf * rb,
+    rd, re, -sc * rd - sf * re
+  ];
 }
 
 /** Get the 'a' component of an Affine transform. */

packages/deck.gl-geotiff/package.json

@@ -46,11 +46,15 @@
     "vitest": "^4.0.18"
   },
   "dependencies": {
+    "@cogeotiff/core": "^9.1.2",
+    "@developmentseed/affine": "workspace:^",
     "@developmentseed/deck.gl-raster": "workspace:^",
+    "@developmentseed/geotiff": "workspace:^",
+    "@developmentseed/morecantile": "workspace:^",
     "@developmentseed/raster-reproject": "workspace:^",
     "flatbush": "^4.5.0",
-    "geotiff": "2.1.3",
-    "proj4": "^2.20.2"
+    "proj4": "^2.20.2",
+    "wkt-parser": "^1.5.2"
   },
   "peerDependencies": {
     "@deck.gl/core": "^9.2.7",