maplibre-gl-raster 0.2.0 → 0.2.2

package/dist/{RasterControl-Dsq7VIl7.js → RasterControl-BxGg08KC.js}

@@ -5,2132 +5,2184 @@ import { CompositeLayer, _GlobeViewport, assert } from "@deck.gl/core";
 import { SimpleMeshLayer } from "@deck.gl/mesh-layers";
 import { TileLayer, _Tileset2D } from "@deck.gl/geo-layers";
 import { Popup } from "maplibre-gl";
-//#region node_modules/@developmentseed/deck.gl-raster/dist/layer-utils.js
-function renderDebugTileOutline(id, tile, forwardTo4326) {
-	const { projectedCorners } = tile;
-	const { topLeft, topRight, bottomRight, bottomLeft } = projectedCorners;
-	const topLeftWgs84 = forwardTo4326(topLeft[0], topLeft[1]);
-	const topRightWgs84 = forwardTo4326(topRight[0], topRight[1]);
-	const bottomRightWgs84 = forwardTo4326(bottomRight[0], bottomRight[1]);
-	const path = [
-		topLeftWgs84,
-		topRightWgs84,
-		bottomRightWgs84,
-		forwardTo4326(bottomLeft[0], bottomLeft[1]),
-		topLeftWgs84
-	];
-	const center = [(topLeftWgs84[0] + bottomRightWgs84[0]) / 2, (topLeftWgs84[1] + bottomRightWgs84[1]) / 2];
-	const labelLayer = new TextLayer({
-		id: `${id}-label`,
-		data: [{
-			position: center,
-			text: `x=${tile.index.x} y=${tile.index.y} z=${tile.index.z}`
-		}],
-		getColor: [
-			255,
-			255,
-			255,
-			255
-		],
-		getSize: 24,
-		sizeUnits: "pixels",
-		outlineWidth: 3,
-		outlineColor: [
-			0,
-			0,
-			0,
-			255
-		],
-		fontSettings: { sdf: true }
-	});
-	return [new PathLayer({
-		id,
-		data: [path],
-		getPath: (d) => d,
-		getColor: [
-			255,
-			0,
-			0,
-			255
-		],
-		getWidth: 2,
-		widthUnits: "pixels",
-		pickable: false
-	}), labelLayer];
-}
-//#endregion
-//#region node_modules/@developmentseed/raster-reproject/dist/delatin.js
+//#region node_modules/@developmentseed/proj/dist/meters-per-unit.js
 /**
-* Define [**uv coordinates**](https://en.wikipedia.org/wiki/UV_mapping) as a float-valued image-local coordinate space where the top left is `(0, 0)` and the bottom right is `(1, 1)`.
-*
-* Define [**Barycentric coordinates**](https://en.wikipedia.org/wiki/Barycentric_coordinate_system) as float-valued triangle-local coordinates, represented as a 3-tuple of floats, where the tuple must add up to 1. The coordinate represents "how close to each vertex" a point in the interior of a triangle is. I.e. `(0, 0, 1)`, `(0, 1, 0)`, and `(1, 0, 0)`  are all valid barycentric coordinates that define one of the three vertices. `(1/3, 1/3, 1/3)` represents the centroid of a triangle. `(1/2, 1/2, 0)` represents a point that is halfway between vertices `a` and `b` and has "none" of vertex `c`.
-*
-*
-* ## Changes
+* Coefficient to convert the coordinate reference system (CRS)
+* units into meters (metersPerUnit).
 *
-* - Delatin coordinates are in terms of pixel space whereas here we use uv space.
+* From note g in http://docs.opengeospatial.org/is/17-083r2/17-083r2.html#table_2:
 *
-* Originally copied from https://github.com/mapbox/delatin under the ISC
-* license, then subject to further modifications.
-*/
-/**
-* Barycentric sample points in uv space for where to sample reprojection
-* errors.
-*/
-var SAMPLE_POINTS = [
-	[
-		1 / 3,
-		1 / 3,
-		1 / 3
-	],
-	[
-		.5,
-		.5,
-		0
-	],
-	[
-		.5,
-		0,
-		.5
-	],
-	[
-		0,
-		.5,
-		.5
-	]
-];
-var DEFAULT_MAX_ERROR$1 = .125;
-/**
-* RasterReprojector performs a Delaunay triangulation-based reprojection of a
-* raster image.
+* > If the CRS uses meters as units of measure for the horizontal dimensions,
+* > then metersPerUnit=1; if it has degrees, then metersPerUnit=2pa/360
+* > (a is the Earth maximum radius of the ellipsoid).
 *
-* It takes as input a set of functions to associate pixel positions with
-* coordinates in the input and output CRS, as well as the dimensions of the
-* output image, and it produces a triangulated mesh that can be used to
-* reproject the input raster onto the output raster with bounded error.
+* @param unit - The unit of the CRS.
+* @param semiMajorAxis - The semi-major axis of the ellipsoid, required if unit is 'degree'.
+* @returns The meters per unit conversion factor.
 */
-var RasterReprojector = class {
-	reprojectors;
-	/** Width of the image in pixels */
-	width;
-	/** Height of the image in pixels */
-	height;
-	/**
-	* UV vertex coordinates (x, y), i.e.
-	* [x0, y0, x1, y1, ...]
-	*
-	* These coordinates are floats that range from [0, 1] in both X and Y.
-	*/
-	uvs;
-	/**
-	* XY Positions in output CRS, computed via exact forward reprojection.
-	*/
-	exactOutputPositions;
-	/**
-	* triangle vertex indices
-	*/
-	triangles;
-	_halfedges;
-	/**
-	* The UV texture coordinates of candidates found from
-	* `findReprojectionCandidate`.
-	*
-	* Maybe in the future we'll want to store the barycentric coordinates instead
-	* of just the uv coordinates?
-	*/
-	_candidatesUV;
-	_queueIndices;
-	_queue;
-	_errors;
-	_pending;
-	_pendingLen;
-	constructor(reprojectors, width, height = width) {
-		this.reprojectors = reprojectors;
-		this.width = width;
-		this.height = height;
-		this.uvs = [];
-		this.exactOutputPositions = [];
-		this.triangles = [];
-		this._halfedges = [];
-		this._candidatesUV = [];
-		this._queueIndices = [];
-		this._queue = [];
-		this._errors = [];
-		this._pending = [];
-		this._pendingLen = 0;
-		const u1 = 1;
-		const v1 = 1;
-		const p0 = this._addPoint(0, 0);
-		const p1 = this._addPoint(u1, 0);
-		const p2 = this._addPoint(0, v1);
-		const p3 = this._addPoint(u1, v1);
-		const t0 = this._addTriangle(p3, p0, p2, -1, -1, -1);
-		this._addTriangle(p0, p3, p1, t0, -1, -1);
-		this._flush();
-	}
-	/**
-	* Refine the mesh until its maximum error gets below the given one
-	*
-	* @param maxError The maximum reprojection error in input pixels that the mesh should achieve.
-	* @param maxIterations Optional safeguard to prevent infinite loops in case of non-convergence. If the mesh fails to converge within this number of iterations, a warning will be logged and the function will return early.
-	*
-	* @return  {[type]}  [return description]
-	*/
-	run(maxError = DEFAULT_MAX_ERROR$1, { maxIterations = 1e4 } = {}) {
-		if (maxError <= 0) throw new Error("maxError must be positive");
-		let iterations = 0;
-		while (this.getMaxError() > maxError) {
-			this.refine();
-			if (++iterations > maxIterations) {
-				console.warn(`RasterReprojector: mesh refinement did not converge after ${iterations} iterations (maxError=${maxError}, currentError=${this.getMaxError()})`);
-				break;
-			}
-		}
+function metersPerUnit(unit, { semiMajorAxis } = {}) {
+	unit = unit.toLowerCase();
+	switch (unit) {
+		case "m":
+		case "metre":
+		case "meter":
+		case "meters": return 1;
+		case "foot": return .3048;
+		case "us survey foot": return 1200 / 3937;
 	}
-	refine() {
-		this._step();
-		this._flush();
+	if (unit === "degree") {
+		if (semiMajorAxis === void 0) throw new Error("CRS with degrees unit requires ellipsoid semi-major-axis");
+		return 2 * Math.PI * semiMajorAxis / 360;
 	}
-	getMaxError() {
-		return this._errors[0];
+	throw new Error(`Unsupported CRS units: ${unit} when computing metersPerUnit.`);
+}
+//#endregion
+//#region node_modules/wkt-parser/PROJJSONBuilderBase.js
+var PROJJSONBuilderBase = class {
+	static getId(node) {
+		const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
+		if (idNode && idNode.length >= 3) return {
+			authority: idNode[1],
+			code: parseInt(idNode[2], 10)
+		};
+		return null;
 	}
-	_flush() {
-		for (let i = 0; i < this._pendingLen; i++) {
-			const t = this._pending[i];
-			this._findReprojectionCandidate(t);
-		}
-		this._pendingLen = 0;
+	static convertUnit(node, type = "unit") {
+		if (!node || node.length < 3) return {
+			type,
+			name: "unknown",
+			conversion_factor: null
+		};
+		const name = node[1];
+		const conversionFactor = parseFloat(node[2]) || null;
+		const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
+		return {
+			type,
+			name,
+			conversion_factor: conversionFactor,
+			id: idNode ? {
+				authority: idNode[1],
+				code: parseInt(idNode[2], 10)
+			} : null
+		};
 	}
-	/**
-	* Conversion of upstream's `_findCandidate` for reprojection error handling.
-	*
-	* @param t The index (into `this.triangles`) of the pending triangle to process.
-	*
-	* @return Doesn't return; instead modifies internal state.
-	*/
-	_findReprojectionCandidate(t) {
-		const a = 2 * this.triangles[t * 3 + 0];
-		const b = 2 * this.triangles[t * 3 + 1];
-		const c = 2 * this.triangles[t * 3 + 2];
-		const p0u = this.uvs[a];
-		const p0v = this.uvs[a + 1];
-		const p1u = this.uvs[b];
-		const p1v = this.uvs[b + 1];
-		const p2u = this.uvs[c];
-		const p2v = this.uvs[c + 1];
-		const out0x = this.exactOutputPositions[a];
-		const out0y = this.exactOutputPositions[a + 1];
-		const out1x = this.exactOutputPositions[b];
-		const out1y = this.exactOutputPositions[b + 1];
-		const out2x = this.exactOutputPositions[c];
-		const out2y = this.exactOutputPositions[c + 1];
-		let maxError = 0;
-		let maxErrorU = 0;
-		let maxErrorV = 0;
-		for (const samplePoint of SAMPLE_POINTS) {
-			const uvSampleU = barycentricMix(p0u, p1u, p2u, samplePoint[0], samplePoint[1], samplePoint[2]);
-			const uvSampleV = barycentricMix(p0v, p1v, p2v, samplePoint[0], samplePoint[1], samplePoint[2]);
-			const outSampleX = barycentricMix(out0x, out1x, out2x, samplePoint[0], samplePoint[1], samplePoint[2]);
-			const outSampleY = barycentricMix(out0y, out1y, out2y, samplePoint[0], samplePoint[1], samplePoint[2]);
-			const pixelExactX = uvSampleU * (this.width - 1);
-			const pixelExactY = uvSampleV * (this.height - 1);
-			const inputCRSSampled = this.reprojectors.inverseReproject(outSampleX, outSampleY);
-			const pixelSampled = this.reprojectors.inverseTransform(inputCRSSampled[0], inputCRSSampled[1]);
-			const dx = pixelExactX - pixelSampled[0];
-			const dy = pixelExactY - pixelSampled[1];
-			const err = Math.hypot(dx, dy);
-			if (err > maxError) {
-				maxError = err;
-				maxErrorU = uvSampleU;
-				maxErrorV = uvSampleV;
-			}
-		}
-		if (maxErrorU === p0u && maxErrorV === p0v || maxErrorU === p1u && maxErrorV === p1v || maxErrorU === p2u && maxErrorV === p2v) maxError = 0;
-		this._candidatesUV[2 * t] = maxErrorU;
-		this._candidatesUV[2 * t + 1] = maxErrorV;
-		this._queuePush(t, maxError);
-	}
-	_step() {
-		const t = this._queuePop();
-		const e0 = t * 3 + 0;
-		const e1 = t * 3 + 1;
-		const e2 = t * 3 + 2;
-		const p0 = this.triangles[e0];
-		const p1 = this.triangles[e1];
-		const p2 = this.triangles[e2];
-		const au = this.uvs[2 * p0];
-		const av = this.uvs[2 * p0 + 1];
-		const bu = this.uvs[2 * p1];
-		const bv = this.uvs[2 * p1 + 1];
-		const cu = this.uvs[2 * p2];
-		const cv = this.uvs[2 * p2 + 1];
-		const pu = this._candidatesUV[2 * t];
-		const pv = this._candidatesUV[2 * t + 1];
-		const pn = this._addPoint(pu, pv);
-		if (orient(au, av, bu, bv, pu, pv) === 0) this._handleCollinear(pn, e0);
-		else if (orient(bu, bv, cu, cv, pu, pv) === 0) this._handleCollinear(pn, e1);
-		else if (orient(cu, cv, au, av, pu, pv) === 0) this._handleCollinear(pn, e2);
-		else {
-			const h0 = this._halfedges[e0];
-			const h1 = this._halfedges[e1];
-			const h2 = this._halfedges[e2];
-			const t0 = this._addTriangle(p0, p1, pn, h0, -1, -1, e0);
-			const t1 = this._addTriangle(p1, p2, pn, h1, -1, t0 + 1);
-			const t2 = this._addTriangle(p2, p0, pn, h2, t0 + 2, t1 + 1);
-			this._legalize(t0);
-			this._legalize(t1);
-			this._legalize(t2);
-		}
-	}
-	_addPoint(u, v) {
-		const i = this.uvs.length >> 1;
-		this.uvs.push(u, v);
-		const pixelX = u * (this.width - 1);
-		const pixelY = v * (this.height - 1);
-		const inputPosition = this.reprojectors.forwardTransform(pixelX, pixelY);
-		const exactOutputPosition = this.reprojectors.forwardReproject(inputPosition[0], inputPosition[1]);
-		this.exactOutputPositions.push(exactOutputPosition[0], exactOutputPosition[1]);
-		return i;
-	}
-	_addTriangle(a, b, c, ab, bc, ca, e = this.triangles.length) {
-		const t = e / 3;
-		this.triangles[e + 0] = a;
-		this.triangles[e + 1] = b;
-		this.triangles[e + 2] = c;
-		this._halfedges[e + 0] = ab;
-		this._halfedges[e + 1] = bc;
-		this._halfedges[e + 2] = ca;
-		if (ab >= 0) this._halfedges[ab] = e + 0;
-		if (bc >= 0) this._halfedges[bc] = e + 1;
-		if (ca >= 0) this._halfedges[ca] = e + 2;
-		this._candidatesUV[2 * t + 0] = 0;
-		this._candidatesUV[2 * t + 1] = 0;
-		this._queueIndices[t] = -1;
-		this._pending[this._pendingLen++] = t;
-		return e;
-	}
-	_legalize(a) {
-		const b = this._halfedges[a];
-		if (b < 0) return;
-		const a0 = a - a % 3;
-		const b0 = b - b % 3;
-		const al = a0 + (a + 1) % 3;
-		const ar = a0 + (a + 2) % 3;
-		const bl = b0 + (b + 2) % 3;
-		const br = b0 + (b + 1) % 3;
-		const p0 = this.triangles[ar];
-		const pr = this.triangles[a];
-		const pl = this.triangles[al];
-		const p1 = this.triangles[bl];
-		const uvs = this.uvs;
-		if (!inCircle(uvs[2 * p0], uvs[2 * p0 + 1], uvs[2 * pr], uvs[2 * pr + 1], uvs[2 * pl], uvs[2 * pl + 1], uvs[2 * p1], uvs[2 * p1 + 1])) return;
-		const hal = this._halfedges[al];
-		const har = this._halfedges[ar];
-		const hbl = this._halfedges[bl];
-		const hbr = this._halfedges[br];
-		this._queueRemove(a0 / 3);
-		this._queueRemove(b0 / 3);
-		const t0 = this._addTriangle(p0, p1, pl, -1, hbl, hal, a0);
-		const t1 = this._addTriangle(p1, p0, pr, t0, har, hbr, b0);
-		this._legalize(t0 + 1);
-		this._legalize(t1 + 2);
-	}
-	_handleCollinear(pn, a) {
-		const a0 = a - a % 3;
-		const al = a0 + (a + 1) % 3;
-		const ar = a0 + (a + 2) % 3;
-		const p0 = this.triangles[ar];
-		const pr = this.triangles[a];
-		const pl = this.triangles[al];
-		const hal = this._halfedges[al];
-		const har = this._halfedges[ar];
-		const b = this._halfedges[a];
-		if (b < 0) {
-			const t0 = this._addTriangle(pn, p0, pr, -1, har, -1, a0);
-			const t1 = this._addTriangle(p0, pn, pl, t0, -1, hal);
-			this._legalize(t0 + 1);
-			this._legalize(t1 + 2);
-			return;
-		}
-		const b0 = b - b % 3;
-		const bl = b0 + (b + 2) % 3;
-		const br = b0 + (b + 1) % 3;
-		const p1 = this.triangles[bl];
-		const hbl = this._halfedges[bl];
-		const hbr = this._halfedges[br];
-		this._queueRemove(b0 / 3);
-		const t0 = this._addTriangle(p0, pr, pn, har, -1, -1, a0);
-		const t1 = this._addTriangle(pr, p1, pn, hbr, -1, t0 + 1, b0);
-		const t2 = this._addTriangle(p1, pl, pn, hbl, -1, t1 + 1);
-		const t3 = this._addTriangle(pl, p0, pn, hal, t0 + 2, t2 + 1);
-		this._legalize(t0);
-		this._legalize(t1);
-		this._legalize(t2);
-		this._legalize(t3);
-	}
-	_queuePush(t, error) {
-		const i = this._queue.length;
-		this._queueIndices[t] = i;
-		this._queue.push(t);
-		this._errors.push(error);
-		this._queueUp(i);
-	}
-	_queuePop() {
-		const n = this._queue.length - 1;
-		this._queueSwap(0, n);
-		this._queueDown(0, n);
-		return this._queuePopBack();
-	}
-	_queuePopBack() {
-		const t = this._queue.pop();
-		this._errors.pop();
-		this._queueIndices[t] = -1;
-		return t;
-	}
-	_queueRemove(t) {
-		const i = this._queueIndices[t];
-		if (i < 0) {
-			const it = this._pending.indexOf(t);
-			if (it !== -1) this._pending[it] = this._pending[--this._pendingLen];
-			else throw new Error("Broken triangulation (something went wrong).");
-			return;
-		}
-		const n = this._queue.length - 1;
-		if (n !== i) {
-			this._queueSwap(i, n);
-			if (!this._queueDown(i, n)) this._queueUp(i);
-		}
-		this._queuePopBack();
-	}
-	_queueLess(i, j) {
-		return this._errors[i] > this._errors[j];
-	}
-	_queueSwap(i, j) {
-		const pi = this._queue[i];
-		const pj = this._queue[j];
-		this._queue[i] = pj;
-		this._queue[j] = pi;
-		this._queueIndices[pi] = j;
-		this._queueIndices[pj] = i;
-		const e = this._errors[i];
-		this._errors[i] = this._errors[j];
-		this._errors[j] = e;
-	}
-	_queueUp(j0) {
-		let j = j0;
-		while (true) {
-			const i = j - 1 >> 1;
-			if (i === j || !this._queueLess(j, i)) break;
-			this._queueSwap(i, j);
-			j = i;
-		}
-	}
-	_queueDown(i0, n) {
-		let i = i0;
-		while (true) {
-			const j1 = 2 * i + 1;
-			if (j1 >= n || j1 < 0) break;
-			const j2 = j1 + 1;
-			let j = j1;
-			if (j2 < n && this._queueLess(j2, j1)) j = j2;
-			if (!this._queueLess(j, i)) break;
-			this._queueSwap(i, j);
-			i = j;
-		}
-		return i > i0;
-	}
-};
-function orient(ax, ay, bx, by, cx, cy) {
-	return (bx - cx) * (ay - cy) - (by - cy) * (ax - cx);
-}
-function inCircle(ax, ay, bx, by, cx, cy, px, py) {
-	const dx = ax - px;
-	const dy = ay - py;
-	const ex = bx - px;
-	const ey = by - py;
-	const fx = cx - px;
-	const fy = cy - py;
-	const ap = dx * dx + dy * dy;
-	const bp = ex * ex + ey * ey;
-	const cp = fx * fx + fy * fy;
-	return dx * (ey * cp - bp * fy) - dy * (ex * cp - bp * fx) + ap * (ex * fy - ey * fx) < 0;
-}
-/**
-* Interpolate the value at a given barycentric coordinate within a triangle.
-*
-* I've seen the name "mix" used before in graphics programming to refer to
-* barycentric linear interpolation.
-*
-* Note: the caller must call this method twice: once for u and once again for
-* v. We do this because we want to avoid allocating an array for the return
-* value.
-*/
-function barycentricMix(a, b, c, t0, t1, t2) {
-	return t0 * a + t1 * b + t2 * c;
-}
-//#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/gpu-modules/create-texture.js
-/**
-* The base shader module for a render pipeline: samples a single input
-* texture into `color` so subsequent modules can transform it. Use this
-* when no decoding step (e.g. {@link CompositeBands}) is needed.
-*/
-var CreateTexture = {
-	name: "create-texture-unorm",
-	inject: {
-		"fs:#decl": `uniform sampler2D textureName;`,
-		"fs:DECKGL_FILTER_COLOR": `
-      color = texture(textureName, geometry.uv);
-    `
-	},
-	getUniforms: (props) => {
-		return { textureName: props.textureName };
-	}
-};
-//#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/mesh-layer/mesh-layer-fragment.glsl.js
-/**
-* This is a vendored copy of the SimpleMeshLayer's fragment shader:
-* https://github.com/visgl/deck.gl/blob/a15c8cea047993c8a861bf542835c1988f30165c/modules/mesh-layers/src/simple-mesh-layer/simple-mesh-layer-fragment.glsl.ts
-* under the MIT license.
-*
-* We edited this to remove the hard-coded texture uniform because we want to
-* support integer and signed integer textures, not only normalized unsigned
-* textures.
-*/
-var mesh_layer_fragment_glsl_default = `#version 300 es
-#define SHADER_NAME simple-mesh-layer-fs
-
-precision highp float;
-
-in vec2 vTexCoord;
-in vec3 cameraPosition;
-in vec3 normals_commonspace;
-in vec4 position_commonspace;
-in vec4 vColor;
-
-out vec4 fragColor;
-
-void main(void) {
-  geometry.uv = vTexCoord;
-
-  vec3 normal;
-  if (simpleMesh.flatShading) {
-
-  normal = normalize(cross(dFdx(position_commonspace.xyz), dFdy(position_commonspace.xyz)));
-  } else {
-    normal = normals_commonspace;
-  }
-
-  // We initialize color here before passing into DECKGL_FILTER_COLOR
-  vec4 color;
-  DECKGL_FILTER_COLOR(color, geometry);
-
-  vec3 lightColor = lighting_getLightColor(color.rgb, cameraPosition, position_commonspace.xyz, normal);
-  fragColor = vec4(lightColor, color.a * layer.opacity);
-}
-`;
-//#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/mesh-layer/mesh-layer.js
-var defaultProps$2 = {
-	...SimpleMeshLayer.defaultProps,
-	renderPipeline: {
-		type: "array",
-		value: [],
-		compare: true
-	},
-	material: {
-		ambient: 1,
-		diffuse: 0,
-		shininess: 0,
-		specularColor: [
-			0,
-			0,
-			0
-		]
-	}
-};
-/**
-* A small subclass of the SimpleMeshLayer to allow dynamic shader injections.
-*
-* In the future this may expand to diverge more from the SimpleMeshLayer, such
-* as allowing the texture to be a 2D _array_.
-*/
-var MeshTextureLayer = class extends SimpleMeshLayer {
-	static layerName = "mesh-texture-layer";
-	static defaultProps = defaultProps$2;
-	_resolveRenderPipeline() {
-		const { image, renderPipeline } = this.props;
-		return [...image ? [{
-			module: CreateTexture,
-			props: { textureName: image }
-		}] : [], ...renderPipeline ?? []];
-	}
-	updateState(params) {
-		if (this.hasRenderPipelineChanged(params)) params.changeFlags.extensionsChanged = true;
-		super.updateState(params);
-	}
-	/** Returns true if the render pipeline has changed between the old and new props. */
-	hasRenderPipelineChanged(params) {
-		const { oldProps, props: newProps } = params;
-		if (Boolean(oldProps.image) !== Boolean(newProps.image)) return true;
-		const oldPipeline = oldProps.renderPipeline ?? [];
-		const newPipeline = newProps.renderPipeline ?? [];
-		if (oldPipeline.length !== newPipeline.length) return true;
-		for (let i = 0; i < oldPipeline.length; i++) if (oldPipeline[i]?.module.name !== newPipeline[i]?.module.name) return true;
-		return false;
-	}
-	getShaders() {
-		const upstreamShaders = super.getShaders();
-		const modules = upstreamShaders.modules;
-		for (const m of this._resolveRenderPipeline()) modules.push(m.module);
-		return {
-			...upstreamShaders,
-			fs: mesh_layer_fragment_glsl_default,
-			modules
-		};
-	}
-	draw(opts) {
-		const shaderProps = {};
-		for (const m of this._resolveRenderPipeline()) shaderProps[m.module.name] = m.props || {};
-		for (const m of super.getModels()) m.shaderInputs.setProps(shaderProps);
-		super.draw(opts);
-	}
-};
-//#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-layer.js
-var DEFAULT_MAX_ERROR = .125;
-var DEBUG_COLORS = [
-	[
-		252,
-		73,
-		163
-	],
-	[
-		255,
-		51,
-		204
-	],
-	[
-		204,
-		102,
-		255
-	],
-	[
-		153,
-		51,
-		255
-	],
-	[
-		102,
-		204,
-		255
-	],
-	[
-		51,
-		153,
-		255
-	],
-	[
-		102,
-		255,
-		204
-	],
-	[
-		51,
-		255,
-		170
-	],
-	[
-		0,
-		255,
-		0
-	],
-	[
-		51,
-		204,
-		51
-	],
-	[
-		255,
-		204,
-		102
-	],
-	[
-		255,
-		179,
-		71
-	],
-	[
-		255,
-		102,
-		102
-	],
-	[
-		255,
-		80,
-		80
-	],
-	[
-		255,
-		0,
-		0
-	],
-	[
-		204,
-		0,
-		0
-	],
-	[
-		255,
-		128,
-		0
-	],
-	[
-		255,
-		153,
-		51
-	],
-	[
-		255,
-		255,
-		102
-	],
-	[
-		255,
-		255,
-		51
-	],
-	[
-		0,
-		255,
-		255
-	],
-	[
-		0,
-		204,
-		255
-	]
-];
-var defaultProps$1 = {
-	image: {
-		type: "image",
-		value: null,
-		async: true
-	},
-	renderPipeline: {
-		type: "array",
-		value: [],
-		compare: true
-	},
-	debug: false,
-	debugOpacity: .5
-};
-/**
-* Generic deck.gl layer for rendering geospatial raster data with client-side,
-* GPU-based reprojection and custom processing pipelines.
-*
-* This is a composite layer that uses {@link RasterReprojector} to generate an adaptive mesh
-* that accurately represents the reprojected raster, then renders it using
-* {@link MeshTextureLayer} (a small wrapper around a deck.gl
-* {@link SimpleMeshLayer}).
-*/
-var RasterLayer = class extends CompositeLayer {
-	static layerName = "RasterLayer";
-	static defaultProps = defaultProps$1;
-	initializeState() {
-		this.setState({});
+	static convertAxis(node) {
+		const name = node[1] || "Unknown";
+		let direction;
+		const abbreviationMatch = name.match(/^\((.)\)$/);
+		if (abbreviationMatch) {
+			const abbreviation = abbreviationMatch[1].toUpperCase();
+			if (abbreviation === "E") direction = "east";
+			else if (abbreviation === "N") direction = "north";
+			else if (abbreviation === "U") direction = "up";
+			else if (node[2]) direction = node[2];
+			else throw new Error(`Unknown axis abbreviation: ${abbreviation}`);
+		} else direction = node[2] || "unknown";
+		const orderNode = node.find((child) => Array.isArray(child) && child[0] === "ORDER");
+		const order = orderNode ? parseInt(orderNode[1], 10) : null;
+		const unitNode = node.find((child) => Array.isArray(child) && (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT"));
+		const unit = this.convertUnit(unitNode);
+		return {
+			name,
+			direction,
+			unit,
+			order
+		};
 	}
-	updateState(params) {
-		super.updateState(params);
-		const { props, oldProps, changeFlags } = params;
-		const reprojectionFnsChanged = props.reprojectionFns.forwardTransform !== oldProps.reprojectionFns?.forwardTransform || props.reprojectionFns.inverseTransform !== oldProps.reprojectionFns?.inverseTransform || props.reprojectionFns.forwardReproject !== oldProps.reprojectionFns?.forwardReproject || props.reprojectionFns.inverseReproject !== oldProps.reprojectionFns?.inverseReproject;
-		if (Boolean(changeFlags.dataChanged) || props.width !== oldProps.width || props.height !== oldProps.height || reprojectionFnsChanged || props.maxError !== oldProps.maxError) this._generateMesh();
+	static extractAxes(node) {
+		return node.filter((child) => Array.isArray(child) && child[0] === "AXIS").map((axis) => this.convertAxis(axis)).sort((a, b) => (a.order || 0) - (b.order || 0));
 	}
-	_generateMesh() {
-		const { width, height, reprojectionFns, maxError = DEFAULT_MAX_ERROR } = this.props;
-		const reprojector = new RasterReprojector(reprojectionFns, width + 1, height + 1);
-		reprojector.run(maxError);
-		const { indices, positions, texCoords } = reprojectorToMesh(reprojector);
-		this.setState({
-			reprojector,
-			mesh: {
-				indices: {
-					value: indices,
-					size: 1
-				},
-				attributes: {
-					POSITION: {
-						value: positions,
-						size: 3
-					},
-					TEXCOORD_0: {
-						value: texCoords,
-						size: 2
-					}
+	static convert(node, result = {}) {
+		switch (node[0]) {
+			case "PROJCRS":
+				result.type = "ProjectedCRS";
+				result.name = node[1];
+				result.base_crs = node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS")) : null;
+				result.conversion = node.find((child) => Array.isArray(child) && child[0] === "CONVERSION") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "CONVERSION")) : null;
+				const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
+				if (csNode) result.coordinate_system = {
+					type: csNode[1],
+					axis: this.extractAxes(node)
+				};
+				const lengthUnitNode = node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT");
+				if (lengthUnitNode) {
+					const unit = this.convertUnit(lengthUnitNode);
+					result.coordinate_system.unit = unit;
+				}
+				result.id = this.getId(node);
+				break;
+			case "BASEGEOGCRS":
+			case "GEOGCRS":
+			case "GEODCRS":
+				result.type = node[0] === "GEODCRS" ? "GeodeticCRS" : "GeographicCRS";
+				result.name = node[1];
+				const datumOrEnsembleNode = node.find((child) => Array.isArray(child) && (child[0] === "DATUM" || child[0] === "ENSEMBLE"));
+				if (datumOrEnsembleNode) {
+					const datumOrEnsemble = this.convert(datumOrEnsembleNode);
+					if (datumOrEnsembleNode[0] === "ENSEMBLE") result.datum_ensemble = datumOrEnsemble;
+					else result.datum = datumOrEnsemble;
+					const primem = node.find((child) => Array.isArray(child) && child[0] === "PRIMEM");
+					if (primem && primem[1] !== "Greenwich") datumOrEnsemble.prime_meridian = {
+						name: primem[1],
+						longitude: parseFloat(primem[2])
+					};
+				}
+				result.coordinate_system = {
+					type: "ellipsoidal",
+					axis: this.extractAxes(node)
+				};
+				result.id = this.getId(node);
+				break;
+			case "DATUM":
+				result.type = "GeodeticReferenceFrame";
+				result.name = node[1];
+				result.ellipsoid = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID")) : null;
+				break;
+			case "ENSEMBLE":
+				result.type = "DatumEnsemble";
+				result.name = node[1];
+				result.members = node.filter((child) => Array.isArray(child) && child[0] === "MEMBER").map((member) => ({
+					type: "DatumEnsembleMember",
+					name: member[1],
+					id: this.getId(member)
+				}));
+				const accuracyNode = node.find((child) => Array.isArray(child) && child[0] === "ENSEMBLEACCURACY");
+				if (accuracyNode) result.accuracy = parseFloat(accuracyNode[1]);
+				const ellipsoidNode = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID");
+				if (ellipsoidNode) result.ellipsoid = this.convert(ellipsoidNode);
+				result.id = this.getId(node);
+				break;
+			case "ELLIPSOID":
+				result.type = "Ellipsoid";
+				result.name = node[1];
+				result.semi_major_axis = parseFloat(node[2]);
+				result.inverse_flattening = parseFloat(node[3]);
+				node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT") && this.convert(node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT"), result);
+				break;
+			case "CONVERSION":
+				result.type = "Conversion";
+				result.name = node[1];
+				result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD")) : null;
+				result.parameters = node.filter((child) => Array.isArray(child) && child[0] === "PARAMETER").map((param) => this.convert(param));
+				break;
+			case "METHOD":
+				result.type = "Method";
+				result.name = node[1];
+				result.id = this.getId(node);
+				break;
+			case "PARAMETER":
+				result.type = "Parameter";
+				result.name = node[1];
+				result.value = parseFloat(node[2]);
+				result.unit = this.convertUnit(node.find((child) => Array.isArray(child) && (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT")));
+				result.id = this.getId(node);
+				break;
+			case "BOUNDCRS":
+				result.type = "BoundCRS";
+				const sourceCrsNode = node.find((child) => Array.isArray(child) && child[0] === "SOURCECRS");
+				if (sourceCrsNode) {
+					const sourceCrsContent = sourceCrsNode.find((child) => Array.isArray(child));
+					result.source_crs = sourceCrsContent ? this.convert(sourceCrsContent) : null;
+				}
+				const targetCrsNode = node.find((child) => Array.isArray(child) && child[0] === "TARGETCRS");
+				if (targetCrsNode) {
+					const targetCrsContent = targetCrsNode.find((child) => Array.isArray(child));
+					result.target_crs = targetCrsContent ? this.convert(targetCrsContent) : null;
+				}
+				const transformationNode = node.find((child) => Array.isArray(child) && child[0] === "ABRIDGEDTRANSFORMATION");
+				if (transformationNode) result.transformation = this.convert(transformationNode);
+				else result.transformation = null;
+				break;
+			case "ABRIDGEDTRANSFORMATION":
+				result.type = "Transformation";
+				result.name = node[1];
+				result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD")) : null;
+				result.parameters = node.filter((child) => Array.isArray(child) && (child[0] === "PARAMETER" || child[0] === "PARAMETERFILE")).map((param) => {
+					if (param[0] === "PARAMETER") return this.convert(param);
+					else if (param[0] === "PARAMETERFILE") return {
+						name: param[1],
+						value: param[2],
+						id: {
+							"authority": "EPSG",
+							"code": 8656
+						}
+					};
+				});
+				if (result.parameters.length === 7) {
+					const scaleDifference = result.parameters[6];
+					if (scaleDifference.name === "Scale difference") scaleDifference.value = Math.round((scaleDifference.value - 1) * 0xe8d4a51000) / 1e6;
 				}
-			}
-		});
-	}
-	renderDebugLayer() {
-		const { reprojector } = this.state;
-		const { debugOpacity } = this.props;
-		if (!reprojector) return null;
-		return new PolygonLayer(this.getSubLayerProps({
-			id: "polygon",
-			data: {
-				reprojector,
-				length: reprojector.triangles.length / 3
-			},
-			getPolygon: (_, { index, data }) => {
-				const triangles = data.reprojector.triangles;
-				const positions = reprojector.exactOutputPositions;
-				const a = triangles[index * 3];
-				const b = triangles[index * 3 + 1];
-				const c = triangles[index * 3 + 2];
-				return [
-					[positions[a * 2], positions[a * 2 + 1]],
-					[positions[b * 2], positions[b * 2 + 1]],
-					[positions[c * 2], positions[c * 2 + 1]],
-					[positions[a * 2], positions[a * 2 + 1]]
-				];
-			},
-			getFillColor: (_, { index, target }) => {
-				const color = DEBUG_COLORS[index % DEBUG_COLORS.length];
-				target[0] = color[0];
-				target[1] = color[1];
-				target[2] = color[2];
-				target[3] = 255;
-				return target;
-			},
-			getLineColor: [
-				0,
-				0,
-				0
-			],
-			getLineWidth: 1,
-			lineWidthUnits: "pixels",
-			opacity: debugOpacity !== void 0 && Number.isFinite(debugOpacity) ? Math.max(0, Math.min(1, debugOpacity)) : 1,
-			pickable: false
-		}));
-	}
-	renderLayers() {
-		const { mesh } = this.state;
-		const { debug, image, renderPipeline } = this.props;
-		if (!mesh || !image && (renderPipeline?.length ?? 0) === 0) return null;
-		const layers = [new MeshTextureLayer(this.getSubLayerProps({
-			id: "raster",
-			image,
-			renderPipeline,
-			data: [1],
-			mesh,
-			_instanced: false,
-			getPosition: [
-				0,
-				0,
-				0
-			],
-			getColor: [
-				255,
-				255,
-				255
-			]
-		}))];
-		if (debug) {
-			const debugLayer = this.renderDebugLayer();
-			if (debugLayer) layers.push(debugLayer);
+				result.id = this.getId(node);
+				break;
+			case "AXIS":
+				if (!result.coordinate_system) result.coordinate_system = {
+					type: "unspecified",
+					axis: []
+				};
+				result.coordinate_system.axis.push(this.convertAxis(node));
+				break;
+			case "LENGTHUNIT":
+				const unit = this.convertUnit(node, "LinearUnit");
+				if (result.coordinate_system && result.coordinate_system.axis) result.coordinate_system.axis.forEach((axis) => {
+					if (!axis.unit) axis.unit = unit;
+				});
+				if (unit.conversion_factor && unit.conversion_factor !== 1) {
+					if (result.semi_major_axis) result.semi_major_axis = {
+						value: result.semi_major_axis,
+						unit
+					};
+				}
+				break;
+			default:
+				result.keyword = node[0];
+				break;
 		}
-		return layers;
+		return result;
 	}
 };
-function reprojectorToMesh(reprojector) {
-	const numVertices = reprojector.uvs.length / 2;
-	const positions = new Float32Array(numVertices * 3);
-	const texCoords = new Float32Array(reprojector.uvs);
-	for (let i = 0; i < numVertices; i++) {
-		positions[i * 3] = reprojector.exactOutputPositions[i * 2];
-		positions[i * 3 + 1] = reprojector.exactOutputPositions[i * 2 + 1];
-		positions[i * 3 + 2] = 0;
+//#endregion
+//#region node_modules/wkt-parser/PROJJSONBuilder2015.js
+var PROJJSONBuilder2015 = class extends PROJJSONBuilderBase {
+	static convert(node, result = {}) {
+		super.convert(node, result);
+		if (result.coordinate_system && result.coordinate_system.subtype === "Cartesian") delete result.coordinate_system;
+		if (result.usage) delete result.usage;
+		return result;
 	}
-	return {
-		indices: new Uint32Array(reprojector.triangles),
-		positions,
-		texCoords
-	};
-}
+};
 //#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-tileset/affine-tileset.js
-/**
-* A {@link RasterTilesetDescriptor} backed by per-level affine transforms.
-*
-* Derives `projectedBounds` from the coarsest level's array. Everything else
-* is passed through from the constructor options.
-*/
-var AffineTileset = class {
-	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;
+//#region node_modules/wkt-parser/PROJJSONBuilder2019.js
+var PROJJSONBuilder2019 = class extends PROJJSONBuilderBase {
+	static convert(node, result = {}) {
+		super.convert(node, result);
+		const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
+		if (csNode) result.coordinate_system = {
+			subtype: csNode[1],
+			axis: this.extractAxes(node)
+		};
+		const usageNode = node.find((child) => Array.isArray(child) && child[0] === "USAGE");
+		if (usageNode) {
+			const scope = usageNode.find((child) => Array.isArray(child) && child[0] === "SCOPE");
+			const area = usageNode.find((child) => Array.isArray(child) && child[0] === "AREA");
+			const bbox = usageNode.find((child) => Array.isArray(child) && child[0] === "BBOX");
+			result.usage = {};
+			if (scope) result.usage.scope = scope[1];
+			if (area) result.usage.area = area[1];
+			if (bbox) result.usage.bbox = bbox.slice(1);
+		}
+		return result;
 	}
 };
 //#endregion
-//#region node_modules/@developmentseed/affine/dist/affine.js
-/**
-* Create a translation transform from an offset vector.
-*
-* @param xoff  Translation offset in x direction.
-* @param yoff  Translation offset in y direction.
-*
-* @return Transform that applies the given translation.
-*/
-function translation(xoff, yoff) {
-	return [
-		1,
-		0,
-		xoff,
-		0,
-		1,
-		yoff
-	];
-}
-/**
-* Create a scaling transform from a scalar or vector.
-*
-* You can pass either one or two scaling factors. Passing only a single scalar
-* value will scale in both dimensions equally. A vector scaling value scales
-* the dimensions independently.
-*
-* @param sx  Scaling factor in x direction.
-* @param sy  Scaling factor in y direction (defaults to sx if not provided).
-*
-* @return Transform that applies the given scaling.
-*/
-function scale$3(sx, sy = sx) {
-	return [
-		sx,
-		0,
-		0,
-		0,
-		sy,
-		0
-	];
-}
-/**
-* Apply a geotransform to a coordinate.
-*
-* That is, we apply this series of equations:
-*
-* ```
-*  x_out = a * x + b * y + c
-*  y_out = d * x + e * y + f
-* ```
-*
-* @param affine  The affine transform to apply.
-* @param x       The x coordinate.
-* @param y       The y coordinate.
-*
-* @return The transformed coordinates.
-*/
-function apply([a, b, c, d, e, f], x, y) {
-	return [a * x + b * y + c, d * x + e * y + f];
-}
+//#region node_modules/wkt-parser/buildPROJJSON.js
 /**
-* Compose two affine transforms: A×B (apply B **first**, then A).
-*
-* This is equivalent to `a @ b` in Python's `affine` library, and is equivalent
-* to multiplying the 3×3 matrices:
-* ```
-*   | a1 b1 c1 |   | a2 b2 c2 |
-*   | d1 e1 f1 | × | d2 e2 f2 |
-*   | 0  0  1  |   | 0  0  1  |
-* ```
-*
-* @param A The first affine transform to apply.
-* @param B The second affine transform to apply.
-*
-* @return The composed affine transform.
+* Detects the WKT2 version based on the structure of the WKT.
+* @param {Array} root The root WKT array node.
+* @returns {string} The detected version ("2015" or "2019").
 */
-function compose([a1, b1, c1, d1, e1, f1], [a2, b2, c2, d2, e2, f2]) {
-	return [
-		a1 * a2 + b1 * d2,
-		a1 * b2 + b1 * e2,
-		a1 * c2 + b1 * f2 + c1,
-		d1 * a2 + e1 * d2,
-		d1 * b2 + e1 * e2,
-		d1 * c2 + e1 * f2 + f1
-	];
+function detectWKT2Version(root) {
+	if (root.find((child) => Array.isArray(child) && child[0] === "USAGE")) return "2019";
+	if (root.find((child) => Array.isArray(child) && child[0] === "CS")) return "2015";
+	if (root[0] === "BOUNDCRS" || root[0] === "PROJCRS" || root[0] === "GEOGCRS") return "2015";
+	return "2015";
 }
 /**
-* Compute the inverse of an Affine.
-*
-* @param affine  The affine transform to invert.
-* @return The inverted affine transform.
-* @throws If the transform is degenerate and cannot be inverted.
+* Builds a PROJJSON object from a WKT array structure.
+* @param {Array} root The root WKT array node.
+* @returns {Object} The PROJJSON object.
 */
-function invert$2([sa, sb, sc, sd, se, sf]) {
-	const det = sa * se - sb * sd;
-	if (det === 0) throw new Error("Cannot invert degenerate transform");
-	const idet = 1 / det;
-	const ra = se * idet;
-	const rb = -sb * idet;
-	const rd = -sd * idet;
-	const re = sa * idet;
-	return [
-		ra,
-		rb,
-		-sc * ra - sf * rb,
-		rd,
-		re,
-		-sc * rd - sf * re
-	];
-}
-/** Get the 'a' component of an Affine transform. */
-function a(affine) {
-	return affine[0];
-}
-/** Get the 'e' component of an Affine transform. */
-function e(affine) {
-	return affine[4];
+function buildPROJJSON(root) {
+	return (detectWKT2Version(root) === "2019" ? PROJJSONBuilder2019 : PROJJSONBuilder2015).convert(root);
 }
 //#endregion
-//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-tileset/affine-tileset-level.js
+//#region node_modules/wkt-parser/detectWKTVersion.js
 /**
-* A {@link RasterTilesetLevel} 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 RasterTilesetLevel} manually.
+* Detects whether the WKT string is WKT1 or WKT2.
+* @param {string} wkt The WKT string.
+* @returns {string} The detected version ("WKT1" or "WKT2").
 */
-var AffineTilesetLevel = class {
-	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 = invert$2(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);
-		const a$1 = a(options.affine);
-		const e$1 = e(options.affine);
-		this.metersPerPixel = Math.sqrt(Math.abs(a$1 * e$1)) * options.mpu;
-		const corners = [
-			apply(options.affine, 0, 0),
-			apply(options.affine, options.arrayWidth, 0),
-			apply(options.affine, 0, options.arrayHeight),
-			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)
-		];
+function detectWKTVersion(wkt) {
+	const normalizedWKT = wkt.toUpperCase();
+	if (normalizedWKT.includes("PROJCRS") || normalizedWKT.includes("GEOGCRS") || normalizedWKT.includes("BOUNDCRS") || normalizedWKT.includes("VERTCRS") || normalizedWKT.includes("LENGTHUNIT") || normalizedWKT.includes("ANGLEUNIT") || normalizedWKT.includes("SCALEUNIT")) return "WKT2";
+	if (normalizedWKT.includes("PROJCS") || normalizedWKT.includes("GEOGCS") || normalizedWKT.includes("LOCAL_CS") || normalizedWKT.includes("VERT_CS") || normalizedWKT.includes("UNIT")) return "WKT1";
+	return "WKT1";
+}
+//#endregion
+//#region node_modules/wkt-parser/parser.js
+var parser_default = parseString;
+var NEUTRAL = 1;
+var KEYWORD = 2;
+var NUMBER = 3;
+var QUOTED = 4;
+var AFTERQUOTE = 5;
+var ENDED = -1;
+var whitespace = /\s/;
+var latin = /[A-Za-z]/;
+var keyword = /[A-Za-z84_]/;
+var endThings = /[,\]]/;
+var digets = /[\d\.E\-\+]/;
+function Parser(text) {
+	if (typeof text !== "string") throw new Error("not a string");
+	this.text = text.trim();
+	this.level = 0;
+	this.place = 0;
+	this.root = null;
+	this.stack = [];
+	this.currentObject = null;
+	this.state = NEUTRAL;
+}
+Parser.prototype.readCharicter = function() {
+	var char = this.text[this.place++];
+	if (this.state !== QUOTED) while (whitespace.test(char)) {
+		if (this.place >= this.text.length) return;
+		char = this.text[this.place++];
+	}
+	switch (this.state) {
+		case NEUTRAL: return this.neutral(char);
+		case KEYWORD: return this.keyword(char);
+		case QUOTED: return this.quoted(char);
+		case AFTERQUOTE: return this.afterquote(char);
+		case NUMBER: return this.number(char);
+		case ENDED: return;
+	}
+};
+Parser.prototype.afterquote = function(char) {
+	if (char === "\"") {
+		this.word += "\"";
+		this.state = QUOTED;
+		return;
+	}
+	if (endThings.test(char)) {
+		this.word = this.word.trim();
+		this.afterItem(char);
+		return;
+	}
+	throw new Error("havn't handled \"" + char + "\" in afterquote yet, index " + this.place);
+};
+Parser.prototype.afterItem = function(char) {
+	if (char === ",") {
+		if (this.word !== null) this.currentObject.push(this.word);
+		this.word = null;
+		this.state = NEUTRAL;
+		return;
+	}
+	if (char === "]") {
+		this.level--;
+		if (this.word !== null) {
+			this.currentObject.push(this.word);
+			this.word = null;
+		}
+		this.state = NEUTRAL;
+		this.currentObject = this.stack.pop();
+		if (!this.currentObject) this.state = ENDED;
+		return;
+	}
+};
+Parser.prototype.number = function(char) {
+	if (digets.test(char)) {
+		this.word += char;
+		return;
+	}
+	if (endThings.test(char)) {
+		this.word = parseFloat(this.word);
+		this.afterItem(char);
+		return;
+	}
+	throw new Error("havn't handled \"" + char + "\" in number yet, index " + this.place);
+};
+Parser.prototype.quoted = function(char) {
+	if (char === "\"") {
+		this.state = AFTERQUOTE;
+		return;
+	}
+	this.word += char;
+};
+Parser.prototype.keyword = function(char) {
+	if (keyword.test(char)) {
+		this.word += char;
+		return;
+	}
+	if (char === "[") {
+		var newObjects = [];
+		newObjects.push(this.word);
+		this.level++;
+		if (this.root === null) this.root = newObjects;
+		else this.currentObject.push(newObjects);
+		this.stack.push(this.currentObject);
+		this.currentObject = newObjects;
+		this.state = NEUTRAL;
+		return;
+	}
+	if (endThings.test(char)) {
+		this.afterItem(char);
+		return;
+	}
+	throw new Error("havn't handled \"" + char + "\" in keyword yet, index " + this.place);
+};
+Parser.prototype.neutral = function(char) {
+	if (latin.test(char)) {
+		this.word = char;
+		this.state = KEYWORD;
+		return;
+	}
+	if (char === "\"") {
+		this.word = "";
+		this.state = QUOTED;
+		return;
+	}
+	if (digets.test(char)) {
+		this.word = char;
+		this.state = NUMBER;
+		return;
+	}
+	if (endThings.test(char)) {
+		this.afterItem(char);
+		return;
+	}
+	throw new Error("havn't handled \"" + char + "\" in neutral yet, index " + this.place);
+};
+Parser.prototype.output = function() {
+	while (this.place < this.text.length) this.readCharicter();
+	if (this.state === ENDED) return this.root;
+	throw new Error("unable to parse string \"" + this.text + "\". State is " + this.state);
+};
+function parseString(txt) {
+	return new Parser(txt).output();
+}
+//#endregion
+//#region node_modules/wkt-parser/process.js
+function mapit(obj, key, value) {
+	if (Array.isArray(key)) {
+		value.unshift(key);
+		key = null;
+	}
+	var thing = key ? {} : obj;
+	var out = value.reduce(function(newObj, item) {
+		sExpr(item, newObj);
+		return newObj;
+	}, thing);
+	if (key) obj[key] = out;
+}
+function sExpr(v, obj) {
+	if (!Array.isArray(v)) {
+		obj[v] = true;
+		return;
 	}
-	projectedTileCorners(col, row) {
-		const tw = this.tileWidth;
-		const th = this.tileHeight;
-		const af = this._affine;
-		return {
-			topLeft: apply(af, col * tw, row * th),
-			topRight: apply(af, (col + 1) * tw, row * th),
-			bottomLeft: apply(af, col * tw, (row + 1) * th),
-			bottomRight: apply(af, (col + 1) * tw, (row + 1) * th)
-		};
+	var key = v.shift();
+	if (key === "PARAMETER") key = v.shift();
+	if (v.length === 1) {
+		if (Array.isArray(v[0])) {
+			obj[key] = {};
+			sExpr(v[0], obj[key]);
+			return;
+		}
+		obj[key] = v[0];
+		return;
 	}
-	tileTransform(col, row) {
-		const tileOffset = translation(col * this.tileWidth, row * this.tileHeight);
-		const tileAffine = compose(this._affine, tileOffset);
-		const invTileAffine = invert$2(tileAffine);
-		return {
-			forwardTransform: (x, y) => apply(tileAffine, x, y),
-			inverseTransform: (x, y) => apply(invTileAffine, x, y)
-		};
+	if (!v.length) {
+		obj[key] = true;
+		return;
 	}
-	crsBoundsToTileRange(projectedMinX, projectedMinY, projectedMaxX, projectedMaxY) {
-		const inv = this._invAffine;
-		const pixelCorners = [
-			apply(inv, projectedMinX, projectedMinY),
-			apply(inv, projectedMaxX, projectedMinY),
-			apply(inv, projectedMinX, projectedMaxY),
-			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;
-		return {
-			minCol: Math.max(0, Math.floor(pixMinX / tw)),
-			maxCol: Math.min(maxColIdx, Math.floor(pixMaxX / tw)),
-			minRow: Math.max(0, Math.floor(pixMinY / th)),
-			maxRow: Math.min(maxRowIdx, Math.floor(pixMaxY / th))
-		};
+	if (key === "TOWGS84") {
+		obj[key] = v;
+		return;
 	}
-};
-//#endregion
-//#region node_modules/@developmentseed/proj/dist/meters-per-unit.js
-/**
-* Coefficient to convert the coordinate reference system (CRS)
-* units into meters (metersPerUnit).
-*
-* From note g in http://docs.opengeospatial.org/is/17-083r2/17-083r2.html#table_2:
-*
-* > If the CRS uses meters as units of measure for the horizontal dimensions,
-* > then metersPerUnit=1; if it has degrees, then metersPerUnit=2pa/360
-* > (a is the Earth maximum radius of the ellipsoid).
-*
-* @param unit - The unit of the CRS.
-* @param semiMajorAxis - The semi-major axis of the ellipsoid, required if unit is 'degree'.
-* @returns The meters per unit conversion factor.
-*/
-function metersPerUnit(unit, { semiMajorAxis } = {}) {
-	unit = unit.toLowerCase();
-	switch (unit) {
-		case "m":
-		case "metre":
-		case "meter":
-		case "meters": return 1;
-		case "foot": return .3048;
-		case "us survey foot": return 1200 / 3937;
+	if (key === "AXIS") {
+		if (!(key in obj)) obj[key] = [];
+		obj[key].push(v);
+		return;
 	}
-	if (unit === "degree") {
-		if (semiMajorAxis === void 0) throw new Error("CRS with degrees unit requires ellipsoid semi-major-axis");
-		return 2 * Math.PI * semiMajorAxis / 360;
+	if (!Array.isArray(key)) obj[key] = {};
+	var i;
+	switch (key) {
+		case "UNIT":
+		case "PRIMEM":
+		case "VERT_DATUM":
+			obj[key] = {
+				name: v[0].toLowerCase(),
+				convert: v[1]
+			};
+			if (v.length === 3) sExpr(v[2], obj[key]);
+			return;
+		case "SPHEROID":
+		case "ELLIPSOID":
+			obj[key] = {
+				name: v[0],
+				a: v[1],
+				rf: v[2]
+			};
+			if (v.length === 4) sExpr(v[3], obj[key]);
+			return;
+		case "EDATUM":
+		case "ENGINEERINGDATUM":
+		case "LOCAL_DATUM":
+		case "DATUM":
+		case "VERT_CS":
+		case "VERTCRS":
+		case "VERTICALCRS":
+			v[0] = ["name", v[0]];
+			mapit(obj, key, v);
+			return;
+		case "COMPD_CS":
+		case "COMPOUNDCRS":
+		case "FITTED_CS":
+		case "PROJECTEDCRS":
+		case "PROJCRS":
+		case "GEOGCS":
+		case "GEOCCS":
+		case "PROJCS":
+		case "LOCAL_CS":
+		case "GEODCRS":
+		case "GEODETICCRS":
+		case "GEODETICDATUM":
+		case "ENGCRS":
+		case "ENGINEERINGCRS":
+			v[0] = ["name", v[0]];
+			mapit(obj, key, v);
+			obj[key].type = key;
+			return;
+		default:
+			i = -1;
+			while (++i < v.length) if (!Array.isArray(v[i])) return sExpr(v, obj[key]);
+			return mapit(obj, key, v);
 	}
-	throw new Error(`Unsupported CRS units: ${unit} when computing metersPerUnit.`);
 }
 //#endregion
-//#region node_modules/wkt-parser/PROJJSONBuilderBase.js
-var PROJJSONBuilderBase = class {
-	static getId(node) {
-		const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
-		if (idNode && idNode.length >= 3) return {
-			authority: idNode[1],
-			code: parseInt(idNode[2], 10)
-		};
-		return null;
+//#region node_modules/wkt-parser/util.js
+var D2R$1 = .017453292519943295;
+function d2r(input) {
+	return input * D2R$1;
+}
+function applyProjectionDefaults(wkt) {
+	const normalizedProjName = (wkt.projName || "").toLowerCase().replace(/_/g, " ");
+	if (wkt.long0 === void 0 && wkt.longc !== void 0) wkt.long0 = wkt.longc;
+	if (!wkt.lat_ts && wkt.lat1 && (normalizedProjName === "stereographic south pole" || normalizedProjName === "polar stereographic (variant b)")) {
+		wkt.lat0 = d2r(wkt.lat1 > 0 ? 90 : -90);
+		wkt.lat_ts = wkt.lat1;
+		delete wkt.lat1;
+	} else if (!wkt.lat_ts && wkt.lat0 && (normalizedProjName === "polar stereographic" || normalizedProjName === "polar stereographic (variant a)")) {
+		wkt.lat_ts = wkt.lat0;
+		wkt.lat0 = d2r(wkt.lat0 > 0 ? 90 : -90);
+		delete wkt.lat1;
 	}
-	static convertUnit(node, type = "unit") {
-		if (!node || node.length < 3) return {
-			type,
-			name: "unknown",
-			conversion_factor: null
-		};
-		const name = node[1];
-		const conversionFactor = parseFloat(node[2]) || null;
-		const idNode = node.find((child) => Array.isArray(child) && child[0] === "ID");
-		return {
-			type,
-			name,
-			conversion_factor: conversionFactor,
-			id: idNode ? {
-				authority: idNode[1],
-				code: parseInt(idNode[2], 10)
-			} : null
-		};
+}
+//#endregion
+//#region node_modules/wkt-parser/transformPROJJSON.js
+function processUnit(unit) {
+	let result = {
+		units: null,
+		to_meter: void 0
+	};
+	if (typeof unit === "string") {
+		result.units = unit.toLowerCase();
+		if (result.units === "metre") result.units = "meter";
+		if (result.units === "meter") result.to_meter = 1;
+	} else if (unit && unit.name) {
+		result.units = unit.name.toLowerCase();
+		if (result.units === "metre") result.units = "meter";
+		result.to_meter = unit.conversion_factor;
 	}
-	static convertAxis(node) {
-		const name = node[1] || "Unknown";
-		let direction;
-		const abbreviationMatch = name.match(/^\((.)\)$/);
-		if (abbreviationMatch) {
-			const abbreviation = abbreviationMatch[1].toUpperCase();
-			if (abbreviation === "E") direction = "east";
-			else if (abbreviation === "N") direction = "north";
-			else if (abbreviation === "U") direction = "up";
-			else if (node[2]) direction = node[2];
-			else throw new Error(`Unknown axis abbreviation: ${abbreviation}`);
-		} else direction = node[2] || "unknown";
-		const orderNode = node.find((child) => Array.isArray(child) && child[0] === "ORDER");
-		const order = orderNode ? parseInt(orderNode[1], 10) : null;
-		const unitNode = node.find((child) => Array.isArray(child) && (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT"));
-		const unit = this.convertUnit(unitNode);
-		return {
-			name,
-			direction,
-			unit,
-			order
-		};
+	return result;
+}
+function toValue(valueOrObject) {
+	if (typeof valueOrObject === "object") return valueOrObject.value * valueOrObject.unit.conversion_factor;
+	return valueOrObject;
+}
+function calculateEllipsoid(value, result) {
+	if (value.ellipsoid.radius) {
+		result.a = value.ellipsoid.radius;
+		result.rf = 0;
+	} else {
+		result.a = toValue(value.ellipsoid.semi_major_axis);
+		if (value.ellipsoid.inverse_flattening !== void 0) result.rf = value.ellipsoid.inverse_flattening;
+		else if (value.ellipsoid.semi_major_axis !== void 0 && value.ellipsoid.semi_minor_axis !== void 0) result.rf = result.a / (result.a - toValue(value.ellipsoid.semi_minor_axis));
 	}
-	static extractAxes(node) {
-		return node.filter((child) => Array.isArray(child) && child[0] === "AXIS").map((axis) => this.convertAxis(axis)).sort((a, b) => (a.order || 0) - (b.order || 0));
+}
+function transformPROJJSON(projjson, result = {}) {
+	if (!projjson || typeof projjson !== "object") return projjson;
+	if (projjson.type === "BoundCRS") {
+		transformPROJJSON(projjson.source_crs, result);
+		if (projjson.transformation) if (projjson.transformation.method && projjson.transformation.method.name === "NTv2") result.nadgrids = projjson.transformation.parameters[0].value;
+		else result.datum_params = projjson.transformation.parameters.map((param) => param.value);
+		return result;
 	}
-	static convert(node, result = {}) {
-		switch (node[0]) {
-			case "PROJCRS":
-				result.type = "ProjectedCRS";
-				result.name = node[1];
-				result.base_crs = node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "BASEGEOGCRS")) : null;
-				result.conversion = node.find((child) => Array.isArray(child) && child[0] === "CONVERSION") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "CONVERSION")) : null;
-				const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
-				if (csNode) result.coordinate_system = {
-					type: csNode[1],
-					axis: this.extractAxes(node)
-				};
-				const lengthUnitNode = node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT");
-				if (lengthUnitNode) {
-					const unit = this.convertUnit(lengthUnitNode);
-					result.coordinate_system.unit = unit;
-				}
-				result.id = this.getId(node);
-				break;
-			case "BASEGEOGCRS":
-			case "GEOGCRS":
-			case "GEODCRS":
-				result.type = node[0] === "GEODCRS" ? "GeodeticCRS" : "GeographicCRS";
-				result.name = node[1];
-				const datumOrEnsembleNode = node.find((child) => Array.isArray(child) && (child[0] === "DATUM" || child[0] === "ENSEMBLE"));
-				if (datumOrEnsembleNode) {
-					const datumOrEnsemble = this.convert(datumOrEnsembleNode);
-					if (datumOrEnsembleNode[0] === "ENSEMBLE") result.datum_ensemble = datumOrEnsemble;
-					else result.datum = datumOrEnsemble;
-					const primem = node.find((child) => Array.isArray(child) && child[0] === "PRIMEM");
-					if (primem && primem[1] !== "Greenwich") datumOrEnsemble.prime_meridian = {
-						name: primem[1],
-						longitude: parseFloat(primem[2])
-					};
-				}
-				result.coordinate_system = {
-					type: "ellipsoidal",
-					axis: this.extractAxes(node)
-				};
-				result.id = this.getId(node);
-				break;
-			case "DATUM":
-				result.type = "GeodeticReferenceFrame";
-				result.name = node[1];
-				result.ellipsoid = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID")) : null;
-				break;
-			case "ENSEMBLE":
-				result.type = "DatumEnsemble";
-				result.name = node[1];
-				result.members = node.filter((child) => Array.isArray(child) && child[0] === "MEMBER").map((member) => ({
-					type: "DatumEnsembleMember",
-					name: member[1],
-					id: this.getId(member)
-				}));
-				const accuracyNode = node.find((child) => Array.isArray(child) && child[0] === "ENSEMBLEACCURACY");
-				if (accuracyNode) result.accuracy = parseFloat(accuracyNode[1]);
-				const ellipsoidNode = node.find((child) => Array.isArray(child) && child[0] === "ELLIPSOID");
-				if (ellipsoidNode) result.ellipsoid = this.convert(ellipsoidNode);
-				result.id = this.getId(node);
-				break;
-			case "ELLIPSOID":
-				result.type = "Ellipsoid";
-				result.name = node[1];
-				result.semi_major_axis = parseFloat(node[2]);
-				result.inverse_flattening = parseFloat(node[3]);
-				node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT") && this.convert(node.find((child) => Array.isArray(child) && child[0] === "LENGTHUNIT"), result);
+	Object.keys(projjson).forEach((key) => {
+		const value = projjson[key];
+		if (value === null) return;
+		switch (key) {
+			case "name":
+				if (result.srsCode) break;
+				result.name = value;
+				result.srsCode = value;
 				break;
-			case "CONVERSION":
-				result.type = "Conversion";
-				result.name = node[1];
-				result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD")) : null;
-				result.parameters = node.filter((child) => Array.isArray(child) && child[0] === "PARAMETER").map((param) => this.convert(param));
+			case "type":
+				if (value === "GeographicCRS") result.projName = "longlat";
+				else if (value === "GeodeticCRS") if (projjson.coordinate_system && projjson.coordinate_system.subtype === "Cartesian") result.projName = "geocent";
+				else result.projName = "longlat";
+				else if (value === "ProjectedCRS" && projjson.conversion && projjson.conversion.method) result.projName = projjson.conversion.method.name;
 				break;
-			case "METHOD":
-				result.type = "Method";
-				result.name = node[1];
-				result.id = this.getId(node);
+			case "datum":
+			case "datum_ensemble":
+				if (value.ellipsoid) {
+					result.ellps = value.ellipsoid.name;
+					calculateEllipsoid(value, result);
+				}
+				if (value.prime_meridian) result.from_greenwich = value.prime_meridian.longitude * Math.PI / 180;
 				break;
-			case "PARAMETER":
-				result.type = "Parameter";
-				result.name = node[1];
-				result.value = parseFloat(node[2]);
-				result.unit = this.convertUnit(node.find((child) => Array.isArray(child) && (child[0] === "LENGTHUNIT" || child[0] === "ANGLEUNIT" || child[0] === "SCALEUNIT")));
-				result.id = this.getId(node);
+			case "ellipsoid":
+				result.ellps = value.name;
+				calculateEllipsoid(value, result);
 				break;
-			case "BOUNDCRS":
-				result.type = "BoundCRS";
-				const sourceCrsNode = node.find((child) => Array.isArray(child) && child[0] === "SOURCECRS");
-				if (sourceCrsNode) {
-					const sourceCrsContent = sourceCrsNode.find((child) => Array.isArray(child));
-					result.source_crs = sourceCrsContent ? this.convert(sourceCrsContent) : null;
-				}
-				const targetCrsNode = node.find((child) => Array.isArray(child) && child[0] === "TARGETCRS");
-				if (targetCrsNode) {
-					const targetCrsContent = targetCrsNode.find((child) => Array.isArray(child));
-					result.target_crs = targetCrsContent ? this.convert(targetCrsContent) : null;
-				}
-				const transformationNode = node.find((child) => Array.isArray(child) && child[0] === "ABRIDGEDTRANSFORMATION");
-				if (transformationNode) result.transformation = this.convert(transformationNode);
-				else result.transformation = null;
+			case "prime_meridian":
+				result.long0 = (value.longitude || 0) * Math.PI / 180;
 				break;
-			case "ABRIDGEDTRANSFORMATION":
-				result.type = "Transformation";
-				result.name = node[1];
-				result.method = node.find((child) => Array.isArray(child) && child[0] === "METHOD") ? this.convert(node.find((child) => Array.isArray(child) && child[0] === "METHOD")) : null;
-				result.parameters = node.filter((child) => Array.isArray(child) && (child[0] === "PARAMETER" || child[0] === "PARAMETERFILE")).map((param) => {
-					if (param[0] === "PARAMETER") return this.convert(param);
-					else if (param[0] === "PARAMETERFILE") return {
-						name: param[1],
-						value: param[2],
-						id: {
-							"authority": "EPSG",
-							"code": 8656
-						}
+			case "coordinate_system":
+				if (value.axis) {
+					const directionMap = {
+						"east": "e",
+						"north": "n",
+						"west": "w",
+						"south": "s",
+						"up": "u",
+						"down": "d",
+						"geocentricx": "e",
+						"geocentricy": "n",
+						"geocentricz": "u"
 					};
-				});
-				if (result.parameters.length === 7) {
-					const scaleDifference = result.parameters[6];
-					if (scaleDifference.name === "Scale difference") scaleDifference.value = Math.round((scaleDifference.value - 1) * 0xe8d4a51000) / 1e6;
+					const mapped = value.axis.map((axis) => directionMap[axis.direction.toLowerCase()]);
+					if (mapped.every(Boolean)) {
+						result.axis = mapped.join("");
+						if (result.axis.length === 2) result.axis += "u";
+					}
+					if (value.unit) {
+						const { units, to_meter } = processUnit(value.unit);
+						result.units = units;
+						result.to_meter = to_meter;
+					} else if (value.axis[0] && value.axis[0].unit) {
+						const { units, to_meter } = processUnit(value.axis[0].unit);
+						result.units = units;
+						result.to_meter = to_meter;
+					}
 				}
-				result.id = this.getId(node);
 				break;
-			case "AXIS":
-				if (!result.coordinate_system) result.coordinate_system = {
-					type: "unspecified",
-					axis: []
-				};
-				result.coordinate_system.axis.push(this.convertAxis(node));
+			case "id":
+				if (value.authority && value.code) result.title = value.authority + ":" + value.code;
 				break;
-			case "LENGTHUNIT":
-				const unit = this.convertUnit(node, "LinearUnit");
-				if (result.coordinate_system && result.coordinate_system.axis) result.coordinate_system.axis.forEach((axis) => {
-					if (!axis.unit) axis.unit = unit;
+			case "conversion":
+				if (value.method && value.method.name) result.projName = value.method.name;
+				if (value.parameters) value.parameters.forEach((param) => {
+					const paramName = param.name.toLowerCase().replace(/\s+/g, "_");
+					const paramValue = param.value;
+					if (param.unit && param.unit.conversion_factor) result[paramName] = paramValue * param.unit.conversion_factor;
+					else if (param.unit === "degree") result[paramName] = paramValue * Math.PI / 180;
+					else result[paramName] = paramValue;
 				});
-				if (unit.conversion_factor && unit.conversion_factor !== 1) {
-					if (result.semi_major_axis) result.semi_major_axis = {
-						value: result.semi_major_axis,
-						unit
-					};
+				break;
+			case "unit":
+				if (value.name) {
+					result.units = value.name.toLowerCase();
+					if (result.units === "metre") result.units = "meter";
 				}
+				if (value.conversion_factor) result.to_meter = value.conversion_factor;
 				break;
-			default:
-				result.keyword = node[0];
+			case "base_crs":
+				transformPROJJSON(value, result);
+				result.datumCode = value.id ? value.id.authority + "_" + value.id.code : value.name;
 				break;
+			default: break;
+		}
+	});
+	if (result.latitude_of_false_origin !== void 0) result.lat0 = result.latitude_of_false_origin;
+	if (result.longitude_of_false_origin !== void 0) result.long0 = result.longitude_of_false_origin;
+	if (result.latitude_of_standard_parallel !== void 0) {
+		result.lat0 = result.latitude_of_standard_parallel;
+		result.lat1 = result.latitude_of_standard_parallel;
+	}
+	if (result.latitude_of_1st_standard_parallel !== void 0) result.lat1 = result.latitude_of_1st_standard_parallel;
+	if (result.latitude_of_2nd_standard_parallel !== void 0) result.lat2 = result.latitude_of_2nd_standard_parallel;
+	if (result.latitude_of_projection_centre !== void 0) result.lat0 = result.latitude_of_projection_centre;
+	if (result.longitude_of_projection_centre !== void 0) result.longc = result.longitude_of_projection_centre;
+	if (result.easting_at_false_origin !== void 0) result.x0 = result.easting_at_false_origin;
+	if (result.northing_at_false_origin !== void 0) result.y0 = result.northing_at_false_origin;
+	if (result.latitude_of_natural_origin !== void 0) result.lat0 = result.latitude_of_natural_origin;
+	if (result.longitude_of_natural_origin !== void 0) result.long0 = result.longitude_of_natural_origin;
+	if (result.longitude_of_origin !== void 0) result.long0 = result.longitude_of_origin;
+	if (result.false_easting !== void 0) result.x0 = result.false_easting;
+	if (result.easting_at_projection_centre) result.x0 = result.easting_at_projection_centre;
+	if (result.false_northing !== void 0) result.y0 = result.false_northing;
+	if (result.northing_at_projection_centre) result.y0 = result.northing_at_projection_centre;
+	if (result.standard_parallel_1 !== void 0) result.lat1 = result.standard_parallel_1;
+	if (result.standard_parallel_2 !== void 0) result.lat2 = result.standard_parallel_2;
+	if (result.scale_factor_at_natural_origin !== void 0) result.k0 = result.scale_factor_at_natural_origin;
+	if (result.scale_factor_at_projection_centre !== void 0) result.k0 = result.scale_factor_at_projection_centre;
+	if (result.scale_factor_on_pseudo_standard_parallel !== void 0) result.k0 = result.scale_factor_on_pseudo_standard_parallel;
+	if (result.azimuth !== void 0) result.alpha = result.azimuth;
+	if (result.azimuth_at_projection_centre !== void 0) result.alpha = result.azimuth_at_projection_centre;
+	if (result.angle_from_rectified_to_skew_grid) result.rectified_grid_angle = result.angle_from_rectified_to_skew_grid;
+	applyProjectionDefaults(result);
+	return result;
+}
+//#endregion
+//#region node_modules/wkt-parser/index.js
+var knownTypes = [
+	"PROJECTEDCRS",
+	"PROJCRS",
+	"GEOGCS",
+	"GEOCCS",
+	"PROJCS",
+	"LOCAL_CS",
+	"GEODCRS",
+	"GEODETICCRS",
+	"GEODETICDATUM",
+	"ENGCRS",
+	"ENGINEERINGCRS"
+];
+function rename(obj, params) {
+	var outName = params[0];
+	var inName = params[1];
+	if (!(outName in obj) && inName in obj) {
+		obj[outName] = obj[inName];
+		if (params.length === 3) obj[outName] = params[2](obj[outName]);
+	}
+}
+function cleanWKT(wkt) {
+	var keys = Object.keys(wkt);
+	for (var i = 0, ii = keys.length; i < ii; ++i) {
+		var key = keys[i];
+		if (knownTypes.indexOf(key) !== -1) setPropertiesFromWkt(wkt[key]);
+		if (typeof wkt[key] === "object") cleanWKT(wkt[key]);
+	}
+}
+function setPropertiesFromWkt(wkt) {
+	if (wkt.AUTHORITY) {
+		var authority = Object.keys(wkt.AUTHORITY)[0];
+		if (authority && authority in wkt.AUTHORITY) wkt.title = authority + ":" + wkt.AUTHORITY[authority];
+	}
+	if (wkt.type === "GEOGCS") wkt.projName = "longlat";
+	else if (wkt.type === "LOCAL_CS") {
+		wkt.projName = "identity";
+		wkt.local = true;
+	} else if (typeof wkt.PROJECTION === "object") wkt.projName = Object.keys(wkt.PROJECTION)[0];
+	else wkt.projName = wkt.PROJECTION;
+	if (wkt.AXIS) {
+		var axisOrder = "";
+		for (var i = 0, ii = wkt.AXIS.length; i < ii; ++i) {
+			var axis = [wkt.AXIS[i][0].toLowerCase(), wkt.AXIS[i][1].toLowerCase()];
+			if (axis[0].indexOf("north") !== -1 || (axis[0] === "y" || axis[0] === "lat") && axis[1] === "north") axisOrder += "n";
+			else if (axis[0].indexOf("south") !== -1 || (axis[0] === "y" || axis[0] === "lat") && axis[1] === "south") axisOrder += "s";
+			else if (axis[0].indexOf("east") !== -1 || (axis[0] === "x" || axis[0] === "lon") && axis[1] === "east") axisOrder += "e";
+			else if (axis[0].indexOf("west") !== -1 || (axis[0] === "x" || axis[0] === "lon") && axis[1] === "west") axisOrder += "w";
 		}
-		return result;
+		if (axisOrder.length === 2) axisOrder += "u";
+		if (axisOrder.length === 3) wkt.axis = axisOrder;
 	}
-};
+	if (wkt.UNIT) {
+		wkt.units = wkt.UNIT.name.toLowerCase();
+		if (wkt.units === "metre") wkt.units = "meter";
+		if (wkt.UNIT.convert) if (wkt.type === "GEOGCS") {
+			if (wkt.DATUM && wkt.DATUM.SPHEROID) wkt.to_meter = wkt.UNIT.convert * wkt.DATUM.SPHEROID.a;
+		} else wkt.to_meter = wkt.UNIT.convert;
+	}
+	var geogcs = wkt.GEOGCS;
+	if (wkt.type === "GEOGCS") geogcs = wkt;
+	if (geogcs) {
+		if (geogcs.PRIMEM && geogcs.PRIMEM.convert) wkt.from_greenwich = d2r(geogcs.PRIMEM.convert);
+		if (geogcs.DATUM) wkt.datumCode = geogcs.DATUM.name.toLowerCase();
+		else wkt.datumCode = geogcs.name.toLowerCase();
+		if (wkt.datumCode.slice(0, 2) === "d_") wkt.datumCode = wkt.datumCode.slice(2);
+		if (wkt.datumCode === "new_zealand_1949") wkt.datumCode = "nzgd49";
+		if (wkt.datumCode === "wgs_1984" || wkt.datumCode === "world_geodetic_system_1984") {
+			if (wkt.PROJECTION === "Mercator_Auxiliary_Sphere") wkt.sphere = true;
+			wkt.datumCode = "wgs84";
+		}
+		if (wkt.datumCode === "belge_1972") wkt.datumCode = "rnb72";
+		if (geogcs.DATUM && geogcs.DATUM.SPHEROID) {
+			wkt.ellps = geogcs.DATUM.SPHEROID.name.replace("_19", "").replace(/[Cc]larke\_18/, "clrk");
+			if (wkt.ellps.toLowerCase().slice(0, 13) === "international") wkt.ellps = "intl";
+			wkt.a = geogcs.DATUM.SPHEROID.a;
+			wkt.rf = parseFloat(geogcs.DATUM.SPHEROID.rf);
+		}
+		if (geogcs.DATUM && geogcs.DATUM.TOWGS84) wkt.datum_params = geogcs.DATUM.TOWGS84;
+		if (~wkt.datumCode.indexOf("osgb_1936")) wkt.datumCode = "osgb36";
+		if (~wkt.datumCode.indexOf("osni_1952")) wkt.datumCode = "osni52";
+		if (~wkt.datumCode.indexOf("tm65") || ~wkt.datumCode.indexOf("geodetic_datum_of_1965")) wkt.datumCode = "ire65";
+		if (wkt.datumCode === "ch1903+") wkt.datumCode = "ch1903";
+		if (~wkt.datumCode.indexOf("israel")) wkt.datumCode = "isr93";
+	}
+	if (wkt.b && !isFinite(wkt.b)) wkt.b = wkt.a;
+	if (wkt.rectified_grid_angle) wkt.rectified_grid_angle = d2r(wkt.rectified_grid_angle);
+	function toMeter(input) {
+		return input * (wkt.to_meter || 1);
+	}
+	var renamer = function(a) {
+		return rename(wkt, a);
+	};
+	[
+		["standard_parallel_1", "Standard_Parallel_1"],
+		["standard_parallel_1", "Latitude of 1st standard parallel"],
+		["standard_parallel_2", "Standard_Parallel_2"],
+		["standard_parallel_2", "Latitude of 2nd standard parallel"],
+		["false_easting", "False_Easting"],
+		["false_easting", "False easting"],
+		["false-easting", "Easting at false origin"],
+		["false_northing", "False_Northing"],
+		["false_northing", "False northing"],
+		["false_northing", "Northing at false origin"],
+		["central_meridian", "Central_Meridian"],
+		["central_meridian", "Longitude of natural origin"],
+		["central_meridian", "Longitude of false origin"],
+		["latitude_of_origin", "Latitude_Of_Origin"],
+		["latitude_of_origin", "Central_Parallel"],
+		["latitude_of_origin", "Latitude of natural origin"],
+		["latitude_of_origin", "Latitude of false origin"],
+		["scale_factor", "Scale_Factor"],
+		["k0", "scale_factor"],
+		["latitude_of_center", "Latitude_Of_Center"],
+		["latitude_of_center", "Latitude_of_center"],
+		[
+			"lat0",
+			"latitude_of_center",
+			d2r
+		],
+		["longitude_of_center", "Longitude_Of_Center"],
+		["longitude_of_center", "Longitude_of_center"],
+		[
+			"longc",
+			"longitude_of_center",
+			d2r
+		],
+		[
+			"x0",
+			"false_easting",
+			toMeter
+		],
+		[
+			"y0",
+			"false_northing",
+			toMeter
+		],
+		[
+			"long0",
+			"central_meridian",
+			d2r
+		],
+		[
+			"lat0",
+			"latitude_of_origin",
+			d2r
+		],
+		[
+			"lat0",
+			"standard_parallel_1",
+			d2r
+		],
+		[
+			"lat1",
+			"standard_parallel_1",
+			d2r
+		],
+		[
+			"lat2",
+			"standard_parallel_2",
+			d2r
+		],
+		["azimuth", "Azimuth"],
+		[
+			"alpha",
+			"azimuth",
+			d2r
+		],
+		["srsCode", "name"]
+	].forEach(renamer);
+	applyProjectionDefaults(wkt);
+}
+function wkt_parser_default(wkt) {
+	if (typeof wkt === "object") return transformPROJJSON(wkt);
+	const version = detectWKTVersion(wkt);
+	var lisp = parser_default(wkt);
+	if (version === "WKT2") return transformPROJJSON(buildPROJJSON(lisp));
+	var type = lisp[0];
+	var obj = {};
+	sExpr(lisp, obj);
+	cleanWKT(obj);
+	return obj[type];
+}
 //#endregion
-//#region node_modules/wkt-parser/PROJJSONBuilder2015.js
-var PROJJSONBuilder2015 = class extends PROJJSONBuilderBase {
-	static convert(node, result = {}) {
-		super.convert(node, result);
-		if (result.coordinate_system && result.coordinate_system.subtype === "Cartesian") delete result.coordinate_system;
-		if (result.usage) delete result.usage;
-		return result;
+//#region node_modules/@developmentseed/proj/dist/parse-wkt.js
+/**
+* Parse a WKT string or PROJJSON object into a proj4-compatible projection
+* definition.
+*
+* This is a typed wrapper around the `wkt-parser` package.
+*/
+function parseWkt(input) {
+	const def = wkt_parser_default(input);
+	if (def.projName === "longlat" && (!def.units || def.units === "unknown")) {
+		def.units = "degree";
+		def.to_meter = void 0;
 	}
-};
+	return def;
+}
 //#endregion
-//#region node_modules/wkt-parser/PROJJSONBuilder2019.js
-var PROJJSONBuilder2019 = class extends PROJJSONBuilderBase {
-	static convert(node, result = {}) {
-		super.convert(node, result);
-		const csNode = node.find((child) => Array.isArray(child) && child[0] === "CS");
-		if (csNode) result.coordinate_system = {
-			subtype: csNode[1],
-			axis: this.extractAxes(node)
-		};
-		const usageNode = node.find((child) => Array.isArray(child) && child[0] === "USAGE");
-		if (usageNode) {
-			const scope = usageNode.find((child) => Array.isArray(child) && child[0] === "SCOPE");
-			const area = usageNode.find((child) => Array.isArray(child) && child[0] === "AREA");
-			const bbox = usageNode.find((child) => Array.isArray(child) && child[0] === "BBOX");
-			result.usage = {};
-			if (scope) result.usage.scope = scope[1];
-			if (area) result.usage.area = area[1];
-			if (bbox) result.usage.bbox = bbox.slice(1);
+//#region node_modules/@developmentseed/proj/dist/registry.js
+/**
+* A global registry holding parsed projection definitions.
+*/
+var PROJECTION_REGISTRY = /* @__PURE__ */ new Map();
+async function epsgResolver(epsg) {
+	const key = `EPSG:${epsg}`;
+	const cachedProj = PROJECTION_REGISTRY.get(key);
+	if (cachedProj !== void 0) return cachedProj;
+	const proj = parseWkt(await getProjjson(epsg));
+	PROJECTION_REGISTRY.set(key, proj);
+	return proj;
+}
+/** Query epsg.io for the PROJJSON corresponding to the given EPSG code. */
+async function getProjjson(epsg) {
+	const url = `https://epsg.io/${epsg}.json`;
+	const resp = await fetch(url);
+	if (!resp.ok) throw new Error(`Failed to fetch PROJJSON from ${url}`);
+	return await resp.json();
+}
+//#endregion
+//#region node_modules/@developmentseed/proj/dist/transform-bounds.js
+/**
+* Transform boundary densifying the edges to account for nonlinear
+* transformations along these edges and extracting the outermost bounds.
+*
+* @param project - function that maps (x, y) in source CRS to (x, y) in target CRS
+* @param left - min X in source CRS
+* @param bottom - min Y in source CRS
+* @param right - max X in source CRS
+* @param top - max Y in source CRS
+* @param options.densifyPts - number of intermediate points along each edge (default 21)
+* @returns [minX, minY, maxX, maxY] in the target CRS
+*/
+function transformBounds(project, left, bottom, right, top, options = {}) {
+	const { densifyPts = 21 } = options;
+	const cx = [
+		left,
+		right,
+		right,
+		left
+	];
+	const cy = [
+		bottom,
+		bottom,
+		top,
+		top
+	];
+	let outMinX = Infinity;
+	let outMinY = Infinity;
+	let outMaxX = -Infinity;
+	let outMaxY = -Infinity;
+	for (let i = 0; i < 4; i++) {
+		const fromX = cx[i];
+		const fromY = cy[i];
+		const toX = cx[(i + 1) % 4];
+		const toY = cy[(i + 1) % 4];
+		for (let j = 0; j <= densifyPts; j++) {
+			const t = j / (densifyPts + 1);
+			const [px, py] = project(fromX + (toX - fromX) * t, fromY + (toY - fromY) * t);
+			if (px < outMinX) outMinX = px;
+			if (py < outMinY) outMinY = py;
+			if (px > outMaxX) outMaxX = px;
+			if (py > outMaxY) outMaxY = py;
 		}
-		return result;
 	}
-};
+	return [
+		outMinX,
+		outMinY,
+		outMaxX,
+		outMaxY
+	];
+}
 //#endregion
-//#region node_modules/wkt-parser/buildPROJJSON.js
+//#region node_modules/@developmentseed/proj/dist/web-mercator.js
+var WGS84_ELLIPSOID_A$1 = 6378137;
+var MAX_WEB_MERCATOR_LAT$1 = 85.05112877980659;
 /**
-* Detects the WKT2 version based on the structure of the WKT.
-* @param {Array} root The root WKT array node.
-* @returns {string} The detected version ("2015" or "2019").
+* Convert a WGS84 longitude/latitude to EPSG:3857 meters analytically.
+* Valid for latitudes in [-MAX_WEB_MERCATOR_LAT, MAX_WEB_MERCATOR_LAT].
 */
-function detectWKT2Version(root) {
-	if (root.find((child) => Array.isArray(child) && child[0] === "USAGE")) return "2019";
-	if (root.find((child) => Array.isArray(child) && child[0] === "CS")) return "2015";
-	if (root[0] === "BOUNDCRS" || root[0] === "PROJCRS" || root[0] === "GEOGCRS") return "2015";
-	return "2015";
+function wgs84To3857(lon, lat) {
+	const x = lon * Math.PI * WGS84_ELLIPSOID_A$1 / 180;
+	const latRad = lat * Math.PI / 180;
+	return [x, Math.log(Math.tan(Math.PI / 4 + latRad / 2)) * WGS84_ELLIPSOID_A$1];
 }
 /**
-* Builds a PROJJSON object from a WKT array structure.
-* @param {Array} root The root WKT array node.
-* @returns {Object} The PROJJSON object.
+* Wrap a proj4 forward projection to EPSG:3857 so that it never returns NaN.
+*
+* proj4 returns [NaN, NaN] for points at the poles (lat = ±90°) because the
+* Mercator projection is undefined there. The wrapper falls back to:
+*   1. Project the input to WGS84 via `forwardTo4326`
+*   2. Clamp the latitude to the Web Mercator limit (±85.05°)
+*   3. Convert analytically from WGS84 to EPSG:3857
+*
+* This correctly handles any input CRS, not just EPSG:4326.
+*
+* NOTE: An identical copy of this function lives in `raster-tile-traversal.ts`.
+* The two packages cannot share code due to their dependency relationship
+* (deck.gl-geotiff depends on deck.gl-raster, not vice versa). If this logic
+* changes, update both copies.
+*
+* Perhaps in the future we'll make a `@developmentseed/projections` package to
+* hold shared projection utilities like this. *
 */
-function buildPROJJSON(root) {
-	return (detectWKT2Version(root) === "2019" ? PROJJSONBuilder2019 : PROJJSONBuilder2015).convert(root);
+function makeClampedForwardTo3857$1(forwardTo3857, forwardTo4326) {
+	return (x, y) => {
+		const [px, py] = forwardTo3857(x, y);
+		if (Number.isFinite(px) && Number.isFinite(py)) return [px, py];
+		const [lon, lat] = forwardTo4326(x, y);
+		return wgs84To3857(lon, Math.max(-85.05112877980659, Math.min(MAX_WEB_MERCATOR_LAT$1, lat)));
+	};
 }
 //#endregion
-//#region node_modules/wkt-parser/detectWKTVersion.js
+//#region src/lib/raster/epsg-resolver.ts
+var OFFLINE_DEFS = new Map([[4326, parseWkt({
+	type: "GeographicCRS",
+	name: "WGS 84",
+	datum: {
+		type: "GeodeticReferenceFrame",
+		name: "World Geodetic System 1984",
+		ellipsoid: {
+			name: "WGS 84",
+			semi_major_axis: 6378137,
+			inverse_flattening: 298.257223563
+		}
+	},
+	coordinate_system: {
+		subtype: "ellipsoidal",
+		axis: [{
+			name: "Geodetic latitude",
+			abbreviation: "Lat",
+			direction: "north",
+			unit: "degree"
+		}, {
+			name: "Geodetic longitude",
+			abbreviation: "Lon",
+			direction: "east",
+			unit: "degree"
+		}]
+	}
+})]]);
 /**
-* Detects whether the WKT string is WKT1 or WKT2.
-* @param {string} wkt The WKT string.
-* @returns {string} The detected version ("WKT1" or "WKT2").
+* Builds an {@link EpsgResolver} that resolves common CRS offline, delegates
+* the rest to {@link ResilientEpsgResolverOptions.fallback} (epsg.io by
+* default), and rethrows failures with a clear, actionable message.
+*
+* @param options - Optional fallback resolver override.
+* @returns A resolver suitable for the COGLayer `epsgResolver` prop.
 */
-function detectWKTVersion(wkt) {
-	const normalizedWKT = wkt.toUpperCase();
-	if (normalizedWKT.includes("PROJCRS") || normalizedWKT.includes("GEOGCRS") || normalizedWKT.includes("BOUNDCRS") || normalizedWKT.includes("VERTCRS") || normalizedWKT.includes("LENGTHUNIT") || normalizedWKT.includes("ANGLEUNIT") || normalizedWKT.includes("SCALEUNIT")) return "WKT2";
-	if (normalizedWKT.includes("PROJCS") || normalizedWKT.includes("GEOGCS") || normalizedWKT.includes("LOCAL_CS") || normalizedWKT.includes("VERT_CS") || normalizedWKT.includes("UNIT")) return "WKT1";
-	return "WKT1";
+function createResilientEpsgResolver(options = {}) {
+	const fallback = options.fallback ?? epsgResolver;
+	return async (epsg) => {
+		const offline = OFFLINE_DEFS.get(epsg);
+		if (offline) return offline;
+		try {
+			return await fallback(epsg);
+		} catch (cause) {
+			const detail = cause instanceof Error ? cause.message : String(cause);
+			const error = /* @__PURE__ */ new Error(`Could not resolve coordinate system EPSG:${epsg}. Coordinate systems are looked up from epsg.io; check your network connection or that the code is valid, then re-add the layer. (${detail})`);
+			if (cause instanceof Error) error.cause = cause;
+			throw error;
+		}
+	};
 }
 //#endregion
-//#region node_modules/wkt-parser/parser.js
-var parser_default = parseString;
-var NEUTRAL = 1;
-var KEYWORD = 2;
-var NUMBER = 3;
-var QUOTED = 4;
-var AFTERQUOTE = 5;
-var ENDED = -1;
-var whitespace = /\s/;
-var latin = /[A-Za-z]/;
-var keyword = /[A-Za-z84_]/;
-var endThings = /[,\]]/;
-var digets = /[\d\.E\-\+]/;
-function Parser(text) {
-	if (typeof text !== "string") throw new Error("not a string");
-	this.text = text.trim();
-	this.level = 0;
-	this.place = 0;
-	this.root = null;
-	this.stack = [];
-	this.currentObject = null;
-	this.state = NEUTRAL;
+//#region node_modules/@developmentseed/deck.gl-raster/dist/layer-utils.js
+function renderDebugTileOutline(id, tile, forwardTo4326) {
+	const { projectedCorners } = tile;
+	const { topLeft, topRight, bottomRight, bottomLeft } = projectedCorners;
+	const topLeftWgs84 = forwardTo4326(topLeft[0], topLeft[1]);
+	const topRightWgs84 = forwardTo4326(topRight[0], topRight[1]);
+	const bottomRightWgs84 = forwardTo4326(bottomRight[0], bottomRight[1]);
+	const path = [
+		topLeftWgs84,
+		topRightWgs84,
+		bottomRightWgs84,
+		forwardTo4326(bottomLeft[0], bottomLeft[1]),
+		topLeftWgs84
+	];
+	const center = [(topLeftWgs84[0] + bottomRightWgs84[0]) / 2, (topLeftWgs84[1] + bottomRightWgs84[1]) / 2];
+	const labelLayer = new TextLayer({
+		id: `${id}-label`,
+		data: [{
+			position: center,
+			text: `x=${tile.index.x} y=${tile.index.y} z=${tile.index.z}`
+		}],
+		getColor: [
+			255,
+			255,
+			255,
+			255
+		],
+		getSize: 24,
+		sizeUnits: "pixels",
+		outlineWidth: 3,
+		outlineColor: [
+			0,
+			0,
+			0,
+			255
+		],
+		fontSettings: { sdf: true }
+	});
+	return [new PathLayer({
+		id,
+		data: [path],
+		getPath: (d) => d,
+		getColor: [
+			255,
+			0,
+			0,
+			255
+		],
+		getWidth: 2,
+		widthUnits: "pixels",
+		pickable: false
+	}), labelLayer];
 }
-Parser.prototype.readCharicter = function() {
-	var char = this.text[this.place++];
-	if (this.state !== QUOTED) while (whitespace.test(char)) {
-		if (this.place >= this.text.length) return;
-		char = this.text[this.place++];
-	}
-	switch (this.state) {
-		case NEUTRAL: return this.neutral(char);
-		case KEYWORD: return this.keyword(char);
-		case QUOTED: return this.quoted(char);
-		case AFTERQUOTE: return this.afterquote(char);
-		case NUMBER: return this.number(char);
-		case ENDED: return;
-	}
-};
-Parser.prototype.afterquote = function(char) {
-	if (char === "\"") {
-		this.word += "\"";
-		this.state = QUOTED;
-		return;
-	}
-	if (endThings.test(char)) {
-		this.word = this.word.trim();
-		this.afterItem(char);
-		return;
-	}
-	throw new Error("havn't handled \"" + char + "\" in afterquote yet, index " + this.place);
-};
-Parser.prototype.afterItem = function(char) {
-	if (char === ",") {
-		if (this.word !== null) this.currentObject.push(this.word);
-		this.word = null;
-		this.state = NEUTRAL;
-		return;
+//#endregion
+//#region node_modules/@developmentseed/raster-reproject/dist/delatin.js
+/**
+* Define [**uv coordinates**](https://en.wikipedia.org/wiki/UV_mapping) as a float-valued image-local coordinate space where the top left is `(0, 0)` and the bottom right is `(1, 1)`.
+*
+* Define [**Barycentric coordinates**](https://en.wikipedia.org/wiki/Barycentric_coordinate_system) as float-valued triangle-local coordinates, represented as a 3-tuple of floats, where the tuple must add up to 1. The coordinate represents "how close to each vertex" a point in the interior of a triangle is. I.e. `(0, 0, 1)`, `(0, 1, 0)`, and `(1, 0, 0)`  are all valid barycentric coordinates that define one of the three vertices. `(1/3, 1/3, 1/3)` represents the centroid of a triangle. `(1/2, 1/2, 0)` represents a point that is halfway between vertices `a` and `b` and has "none" of vertex `c`.
+*
+*
+* ## Changes
+*
+* - Delatin coordinates are in terms of pixel space whereas here we use uv space.
+*
+* Originally copied from https://github.com/mapbox/delatin under the ISC
+* license, then subject to further modifications.
+*/
+/**
+* Barycentric sample points in uv space for where to sample reprojection
+* errors.
+*/
+var SAMPLE_POINTS = [
+	[
+		1 / 3,
+		1 / 3,
+		1 / 3
+	],
+	[
+		.5,
+		.5,
+		0
+	],
+	[
+		.5,
+		0,
+		.5
+	],
+	[
+		0,
+		.5,
+		.5
+	]
+];
+var DEFAULT_MAX_ERROR$1 = .125;
+/**
+* RasterReprojector performs a Delaunay triangulation-based reprojection of a
+* raster image.
+*
+* It takes as input a set of functions to associate pixel positions with
+* coordinates in the input and output CRS, as well as the dimensions of the
+* output image, and it produces a triangulated mesh that can be used to
+* reproject the input raster onto the output raster with bounded error.
+*/
+var RasterReprojector = class {
+	reprojectors;
+	/** Width of the image in pixels */
+	width;
+	/** Height of the image in pixels */
+	height;
+	/**
+	* UV vertex coordinates (x, y), i.e.
+	* [x0, y0, x1, y1, ...]
+	*
+	* These coordinates are floats that range from [0, 1] in both X and Y.
+	*/
+	uvs;
+	/**
+	* XY Positions in output CRS, computed via exact forward reprojection.
+	*/
+	exactOutputPositions;
+	/**
+	* triangle vertex indices
+	*/
+	triangles;
+	_halfedges;
+	/**
+	* The UV texture coordinates of candidates found from
+	* `findReprojectionCandidate`.
+	*
+	* Maybe in the future we'll want to store the barycentric coordinates instead
+	* of just the uv coordinates?
+	*/
+	_candidatesUV;
+	_queueIndices;
+	_queue;
+	_errors;
+	_pending;
+	_pendingLen;
+	constructor(reprojectors, width, height = width) {
+		this.reprojectors = reprojectors;
+		this.width = width;
+		this.height = height;
+		this.uvs = [];
+		this.exactOutputPositions = [];
+		this.triangles = [];
+		this._halfedges = [];
+		this._candidatesUV = [];
+		this._queueIndices = [];
+		this._queue = [];
+		this._errors = [];
+		this._pending = [];
+		this._pendingLen = 0;
+		const u1 = 1;
+		const v1 = 1;
+		const p0 = this._addPoint(0, 0);
+		const p1 = this._addPoint(u1, 0);
+		const p2 = this._addPoint(0, v1);
+		const p3 = this._addPoint(u1, v1);
+		const t0 = this._addTriangle(p3, p0, p2, -1, -1, -1);
+		this._addTriangle(p0, p3, p1, t0, -1, -1);
+		this._flush();
 	}
-	if (char === "]") {
-		this.level--;
-		if (this.word !== null) {
-			this.currentObject.push(this.word);
-			this.word = null;
+	/**
+	* Refine the mesh until its maximum error gets below the given one
+	*
+	* @param maxError The maximum reprojection error in input pixels that the mesh should achieve.
+	* @param maxIterations Optional safeguard to prevent infinite loops in case of non-convergence. If the mesh fails to converge within this number of iterations, a warning will be logged and the function will return early.
+	*
+	* @return  {[type]}  [return description]
+	*/
+	run(maxError = DEFAULT_MAX_ERROR$1, { maxIterations = 1e4 } = {}) {
+		if (maxError <= 0) throw new Error("maxError must be positive");
+		let iterations = 0;
+		while (this.getMaxError() > maxError) {
+			this.refine();
+			if (++iterations > maxIterations) {
+				console.warn(`RasterReprojector: mesh refinement did not converge after ${iterations} iterations (maxError=${maxError}, currentError=${this.getMaxError()})`);
+				break;
+			}
 		}
-		this.state = NEUTRAL;
-		this.currentObject = this.stack.pop();
-		if (!this.currentObject) this.state = ENDED;
-		return;
 	}
-};
-Parser.prototype.number = function(char) {
-	if (digets.test(char)) {
-		this.word += char;
-		return;
+	refine() {
+		this._step();
+		this._flush();
 	}
-	if (endThings.test(char)) {
-		this.word = parseFloat(this.word);
-		this.afterItem(char);
-		return;
+	getMaxError() {
+		return this._errors[0];
 	}
-	throw new Error("havn't handled \"" + char + "\" in number yet, index " + this.place);
-};
-Parser.prototype.quoted = function(char) {
-	if (char === "\"") {
-		this.state = AFTERQUOTE;
-		return;
+	_flush() {
+		for (let i = 0; i < this._pendingLen; i++) {
+			const t = this._pending[i];
+			this._findReprojectionCandidate(t);
+		}
+		this._pendingLen = 0;
 	}
-	this.word += char;
-};
-Parser.prototype.keyword = function(char) {
-	if (keyword.test(char)) {
-		this.word += char;
-		return;
+	/**
+	* Conversion of upstream's `_findCandidate` for reprojection error handling.
+	*
+	* @param t The index (into `this.triangles`) of the pending triangle to process.
+	*
+	* @return Doesn't return; instead modifies internal state.
+	*/
+	_findReprojectionCandidate(t) {
+		const a = 2 * this.triangles[t * 3 + 0];
+		const b = 2 * this.triangles[t * 3 + 1];
+		const c = 2 * this.triangles[t * 3 + 2];
+		const p0u = this.uvs[a];
+		const p0v = this.uvs[a + 1];
+		const p1u = this.uvs[b];
+		const p1v = this.uvs[b + 1];
+		const p2u = this.uvs[c];
+		const p2v = this.uvs[c + 1];
+		const out0x = this.exactOutputPositions[a];
+		const out0y = this.exactOutputPositions[a + 1];
+		const out1x = this.exactOutputPositions[b];
+		const out1y = this.exactOutputPositions[b + 1];
+		const out2x = this.exactOutputPositions[c];
+		const out2y = this.exactOutputPositions[c + 1];
+		let maxError = 0;
+		let maxErrorU = 0;
+		let maxErrorV = 0;
+		for (const samplePoint of SAMPLE_POINTS) {
+			const uvSampleU = barycentricMix(p0u, p1u, p2u, samplePoint[0], samplePoint[1], samplePoint[2]);
+			const uvSampleV = barycentricMix(p0v, p1v, p2v, samplePoint[0], samplePoint[1], samplePoint[2]);
+			const outSampleX = barycentricMix(out0x, out1x, out2x, samplePoint[0], samplePoint[1], samplePoint[2]);
+			const outSampleY = barycentricMix(out0y, out1y, out2y, samplePoint[0], samplePoint[1], samplePoint[2]);
+			const pixelExactX = uvSampleU * (this.width - 1);
+			const pixelExactY = uvSampleV * (this.height - 1);
+			const inputCRSSampled = this.reprojectors.inverseReproject(outSampleX, outSampleY);
+			const pixelSampled = this.reprojectors.inverseTransform(inputCRSSampled[0], inputCRSSampled[1]);
+			const dx = pixelExactX - pixelSampled[0];
+			const dy = pixelExactY - pixelSampled[1];
+			const err = Math.hypot(dx, dy);
+			if (err > maxError) {
+				maxError = err;
+				maxErrorU = uvSampleU;
+				maxErrorV = uvSampleV;
+			}
+		}
+		if (maxErrorU === p0u && maxErrorV === p0v || maxErrorU === p1u && maxErrorV === p1v || maxErrorU === p2u && maxErrorV === p2v) maxError = 0;
+		this._candidatesUV[2 * t] = maxErrorU;
+		this._candidatesUV[2 * t + 1] = maxErrorV;
+		this._queuePush(t, maxError);
 	}
-	if (char === "[") {
-		var newObjects = [];
-		newObjects.push(this.word);
-		this.level++;
-		if (this.root === null) this.root = newObjects;
-		else this.currentObject.push(newObjects);
-		this.stack.push(this.currentObject);
-		this.currentObject = newObjects;
-		this.state = NEUTRAL;
-		return;
+	_step() {
+		const t = this._queuePop();
+		const e0 = t * 3 + 0;
+		const e1 = t * 3 + 1;
+		const e2 = t * 3 + 2;
+		const p0 = this.triangles[e0];
+		const p1 = this.triangles[e1];
+		const p2 = this.triangles[e2];
+		const au = this.uvs[2 * p0];
+		const av = this.uvs[2 * p0 + 1];
+		const bu = this.uvs[2 * p1];
+		const bv = this.uvs[2 * p1 + 1];
+		const cu = this.uvs[2 * p2];
+		const cv = this.uvs[2 * p2 + 1];
+		const pu = this._candidatesUV[2 * t];
+		const pv = this._candidatesUV[2 * t + 1];
+		const pn = this._addPoint(pu, pv);
+		if (orient(au, av, bu, bv, pu, pv) === 0) this._handleCollinear(pn, e0);
+		else if (orient(bu, bv, cu, cv, pu, pv) === 0) this._handleCollinear(pn, e1);
+		else if (orient(cu, cv, au, av, pu, pv) === 0) this._handleCollinear(pn, e2);
+		else {
+			const h0 = this._halfedges[e0];
+			const h1 = this._halfedges[e1];
+			const h2 = this._halfedges[e2];
+			const t0 = this._addTriangle(p0, p1, pn, h0, -1, -1, e0);
+			const t1 = this._addTriangle(p1, p2, pn, h1, -1, t0 + 1);
+			const t2 = this._addTriangle(p2, p0, pn, h2, t0 + 2, t1 + 1);
+			this._legalize(t0);
+			this._legalize(t1);
+			this._legalize(t2);
+		}
 	}
-	if (endThings.test(char)) {
-		this.afterItem(char);
-		return;
+	_addPoint(u, v) {
+		const i = this.uvs.length >> 1;
+		this.uvs.push(u, v);
+		const pixelX = u * (this.width - 1);
+		const pixelY = v * (this.height - 1);
+		const inputPosition = this.reprojectors.forwardTransform(pixelX, pixelY);
+		const exactOutputPosition = this.reprojectors.forwardReproject(inputPosition[0], inputPosition[1]);
+		this.exactOutputPositions.push(exactOutputPosition[0], exactOutputPosition[1]);
+		return i;
 	}
-	throw new Error("havn't handled \"" + char + "\" in keyword yet, index " + this.place);
-};
-Parser.prototype.neutral = function(char) {
-	if (latin.test(char)) {
-		this.word = char;
-		this.state = KEYWORD;
-		return;
+	_addTriangle(a, b, c, ab, bc, ca, e = this.triangles.length) {
+		const t = e / 3;
+		this.triangles[e + 0] = a;
+		this.triangles[e + 1] = b;
+		this.triangles[e + 2] = c;
+		this._halfedges[e + 0] = ab;
+		this._halfedges[e + 1] = bc;
+		this._halfedges[e + 2] = ca;
+		if (ab >= 0) this._halfedges[ab] = e + 0;
+		if (bc >= 0) this._halfedges[bc] = e + 1;
+		if (ca >= 0) this._halfedges[ca] = e + 2;
+		this._candidatesUV[2 * t + 0] = 0;
+		this._candidatesUV[2 * t + 1] = 0;
+		this._queueIndices[t] = -1;
+		this._pending[this._pendingLen++] = t;
+		return e;
 	}
-	if (char === "\"") {
-		this.word = "";
-		this.state = QUOTED;
-		return;
+	_legalize(a) {
+		const b = this._halfedges[a];
+		if (b < 0) return;
+		const a0 = a - a % 3;
+		const b0 = b - b % 3;
+		const al = a0 + (a + 1) % 3;
+		const ar = a0 + (a + 2) % 3;
+		const bl = b0 + (b + 2) % 3;
+		const br = b0 + (b + 1) % 3;
+		const p0 = this.triangles[ar];
+		const pr = this.triangles[a];
+		const pl = this.triangles[al];
+		const p1 = this.triangles[bl];
+		const uvs = this.uvs;
+		if (!inCircle(uvs[2 * p0], uvs[2 * p0 + 1], uvs[2 * pr], uvs[2 * pr + 1], uvs[2 * pl], uvs[2 * pl + 1], uvs[2 * p1], uvs[2 * p1 + 1])) return;
+		const hal = this._halfedges[al];
+		const har = this._halfedges[ar];
+		const hbl = this._halfedges[bl];
+		const hbr = this._halfedges[br];
+		this._queueRemove(a0 / 3);
+		this._queueRemove(b0 / 3);
+		const t0 = this._addTriangle(p0, p1, pl, -1, hbl, hal, a0);
+		const t1 = this._addTriangle(p1, p0, pr, t0, har, hbr, b0);
+		this._legalize(t0 + 1);
+		this._legalize(t1 + 2);
 	}
-	if (digets.test(char)) {
-		this.word = char;
-		this.state = NUMBER;
-		return;
+	_handleCollinear(pn, a) {
+		const a0 = a - a % 3;
+		const al = a0 + (a + 1) % 3;
+		const ar = a0 + (a + 2) % 3;
+		const p0 = this.triangles[ar];
+		const pr = this.triangles[a];
+		const pl = this.triangles[al];
+		const hal = this._halfedges[al];
+		const har = this._halfedges[ar];
+		const b = this._halfedges[a];
+		if (b < 0) {
+			const t0 = this._addTriangle(pn, p0, pr, -1, har, -1, a0);
+			const t1 = this._addTriangle(p0, pn, pl, t0, -1, hal);
+			this._legalize(t0 + 1);
+			this._legalize(t1 + 2);
+			return;
+		}
+		const b0 = b - b % 3;
+		const bl = b0 + (b + 2) % 3;
+		const br = b0 + (b + 1) % 3;
+		const p1 = this.triangles[bl];
+		const hbl = this._halfedges[bl];
+		const hbr = this._halfedges[br];
+		this._queueRemove(b0 / 3);
+		const t0 = this._addTriangle(p0, pr, pn, har, -1, -1, a0);
+		const t1 = this._addTriangle(pr, p1, pn, hbr, -1, t0 + 1, b0);
+		const t2 = this._addTriangle(p1, pl, pn, hbl, -1, t1 + 1);
+		const t3 = this._addTriangle(pl, p0, pn, hal, t0 + 2, t2 + 1);
+		this._legalize(t0);
+		this._legalize(t1);
+		this._legalize(t2);
+		this._legalize(t3);
 	}
-	if (endThings.test(char)) {
-		this.afterItem(char);
-		return;
+	_queuePush(t, error) {
+		const i = this._queue.length;
+		this._queueIndices[t] = i;
+		this._queue.push(t);
+		this._errors.push(error);
+		this._queueUp(i);
 	}
-	throw new Error("havn't handled \"" + char + "\" in neutral yet, index " + this.place);
-};
-Parser.prototype.output = function() {
-	while (this.place < this.text.length) this.readCharicter();
-	if (this.state === ENDED) return this.root;
-	throw new Error("unable to parse string \"" + this.text + "\". State is " + this.state);
-};
-function parseString(txt) {
-	return new Parser(txt).output();
-}
-//#endregion
-//#region node_modules/wkt-parser/process.js
-function mapit(obj, key, value) {
-	if (Array.isArray(key)) {
-		value.unshift(key);
-		key = null;
+	_queuePop() {
+		const n = this._queue.length - 1;
+		this._queueSwap(0, n);
+		this._queueDown(0, n);
+		return this._queuePopBack();
 	}
-	var thing = key ? {} : obj;
-	var out = value.reduce(function(newObj, item) {
-		sExpr(item, newObj);
-		return newObj;
-	}, thing);
-	if (key) obj[key] = out;
-}
-function sExpr(v, obj) {
-	if (!Array.isArray(v)) {
-		obj[v] = true;
-		return;
+	_queuePopBack() {
+		const t = this._queue.pop();
+		this._errors.pop();
+		this._queueIndices[t] = -1;
+		return t;
 	}
-	var key = v.shift();
-	if (key === "PARAMETER") key = v.shift();
-	if (v.length === 1) {
-		if (Array.isArray(v[0])) {
-			obj[key] = {};
-			sExpr(v[0], obj[key]);
+	_queueRemove(t) {
+		const i = this._queueIndices[t];
+		if (i < 0) {
+			const it = this._pending.indexOf(t);
+			if (it !== -1) this._pending[it] = this._pending[--this._pendingLen];
+			else throw new Error("Broken triangulation (something went wrong).");
 			return;
 		}
-		obj[key] = v[0];
-		return;
+		const n = this._queue.length - 1;
+		if (n !== i) {
+			this._queueSwap(i, n);
+			if (!this._queueDown(i, n)) this._queueUp(i);
+		}
+		this._queuePopBack();
 	}
-	if (!v.length) {
-		obj[key] = true;
-		return;
+	_queueLess(i, j) {
+		return this._errors[i] > this._errors[j];
 	}
-	if (key === "TOWGS84") {
-		obj[key] = v;
-		return;
+	_queueSwap(i, j) {
+		const pi = this._queue[i];
+		const pj = this._queue[j];
+		this._queue[i] = pj;
+		this._queue[j] = pi;
+		this._queueIndices[pi] = j;
+		this._queueIndices[pj] = i;
+		const e = this._errors[i];
+		this._errors[i] = this._errors[j];
+		this._errors[j] = e;
 	}
-	if (key === "AXIS") {
-		if (!(key in obj)) obj[key] = [];
-		obj[key].push(v);
-		return;
+	_queueUp(j0) {
+		let j = j0;
+		while (true) {
+			const i = j - 1 >> 1;
+			if (i === j || !this._queueLess(j, i)) break;
+			this._queueSwap(i, j);
+			j = i;
+		}
 	}
-	if (!Array.isArray(key)) obj[key] = {};
-	var i;
-	switch (key) {
-		case "UNIT":
-		case "PRIMEM":
-		case "VERT_DATUM":
-			obj[key] = {
-				name: v[0].toLowerCase(),
-				convert: v[1]
-			};
-			if (v.length === 3) sExpr(v[2], obj[key]);
-			return;
-		case "SPHEROID":
-		case "ELLIPSOID":
-			obj[key] = {
-				name: v[0],
-				a: v[1],
-				rf: v[2]
-			};
-			if (v.length === 4) sExpr(v[3], obj[key]);
-			return;
-		case "EDATUM":
-		case "ENGINEERINGDATUM":
-		case "LOCAL_DATUM":
-		case "DATUM":
-		case "VERT_CS":
-		case "VERTCRS":
-		case "VERTICALCRS":
-			v[0] = ["name", v[0]];
-			mapit(obj, key, v);
-			return;
-		case "COMPD_CS":
-		case "COMPOUNDCRS":
-		case "FITTED_CS":
-		case "PROJECTEDCRS":
-		case "PROJCRS":
-		case "GEOGCS":
-		case "GEOCCS":
-		case "PROJCS":
-		case "LOCAL_CS":
-		case "GEODCRS":
-		case "GEODETICCRS":
-		case "GEODETICDATUM":
-		case "ENGCRS":
-		case "ENGINEERINGCRS":
-			v[0] = ["name", v[0]];
-			mapit(obj, key, v);
-			obj[key].type = key;
-			return;
-		default:
-			i = -1;
-			while (++i < v.length) if (!Array.isArray(v[i])) return sExpr(v, obj[key]);
-			return mapit(obj, key, v);
+	_queueDown(i0, n) {
+		let i = i0;
+		while (true) {
+			const j1 = 2 * i + 1;
+			if (j1 >= n || j1 < 0) break;
+			const j2 = j1 + 1;
+			let j = j1;
+			if (j2 < n && this._queueLess(j2, j1)) j = j2;
+			if (!this._queueLess(j, i)) break;
+			this._queueSwap(i, j);
+			i = j;
+		}
+		return i > i0;
 	}
+};
+function orient(ax, ay, bx, by, cx, cy) {
+	return (bx - cx) * (ay - cy) - (by - cy) * (ax - cx);
 }
-//#endregion
-//#region node_modules/wkt-parser/util.js
-var D2R$1 = .017453292519943295;
-function d2r(input) {
-	return input * D2R$1;
+function inCircle(ax, ay, bx, by, cx, cy, px, py) {
+	const dx = ax - px;
+	const dy = ay - py;
+	const ex = bx - px;
+	const ey = by - py;
+	const fx = cx - px;
+	const fy = cy - py;
+	const ap = dx * dx + dy * dy;
+	const bp = ex * ex + ey * ey;
+	const cp = fx * fx + fy * fy;
+	return dx * (ey * cp - bp * fy) - dy * (ex * cp - bp * fx) + ap * (ex * fy - ey * fx) < 0;
 }
-function applyProjectionDefaults(wkt) {
-	const normalizedProjName = (wkt.projName || "").toLowerCase().replace(/_/g, " ");
-	if (wkt.long0 === void 0 && wkt.longc !== void 0) wkt.long0 = wkt.longc;
-	if (!wkt.lat_ts && wkt.lat1 && (normalizedProjName === "stereographic south pole" || normalizedProjName === "polar stereographic (variant b)")) {
-		wkt.lat0 = d2r(wkt.lat1 > 0 ? 90 : -90);
-		wkt.lat_ts = wkt.lat1;
-		delete wkt.lat1;
-	} else if (!wkt.lat_ts && wkt.lat0 && (normalizedProjName === "polar stereographic" || normalizedProjName === "polar stereographic (variant a)")) {
-		wkt.lat_ts = wkt.lat0;
-		wkt.lat0 = d2r(wkt.lat0 > 0 ? 90 : -90);
-		delete wkt.lat1;
+/**
+* Interpolate the value at a given barycentric coordinate within a triangle.
+*
+* I've seen the name "mix" used before in graphics programming to refer to
+* barycentric linear interpolation.
+*
+* Note: the caller must call this method twice: once for u and once again for
+* v. We do this because we want to avoid allocating an array for the return
+* value.
+*/
+function barycentricMix(a, b, c, t0, t1, t2) {
+	return t0 * a + t1 * b + t2 * c;
+}
+//#endregion
+//#region node_modules/@developmentseed/deck.gl-raster/dist/gpu-modules/create-texture.js
+/**
+* The base shader module for a render pipeline: samples a single input
+* texture into `color` so subsequent modules can transform it. Use this
+* when no decoding step (e.g. {@link CompositeBands}) is needed.
+*/
+var CreateTexture = {
+	name: "create-texture-unorm",
+	inject: {
+		"fs:#decl": `uniform sampler2D textureName;`,
+		"fs:DECKGL_FILTER_COLOR": `
+      color = texture(textureName, geometry.uv);
+    `
+	},
+	getUniforms: (props) => {
+		return { textureName: props.textureName };
 	}
+};
+//#endregion
+//#region node_modules/@developmentseed/deck.gl-raster/dist/mesh-layer/mesh-layer-fragment.glsl.js
+/**
+* This is a vendored copy of the SimpleMeshLayer's fragment shader:
+* https://github.com/visgl/deck.gl/blob/a15c8cea047993c8a861bf542835c1988f30165c/modules/mesh-layers/src/simple-mesh-layer/simple-mesh-layer-fragment.glsl.ts
+* under the MIT license.
+*
+* We edited this to remove the hard-coded texture uniform because we want to
+* support integer and signed integer textures, not only normalized unsigned
+* textures.
+*/
+var mesh_layer_fragment_glsl_default = `#version 300 es
+#define SHADER_NAME simple-mesh-layer-fs
+
+precision highp float;
+
+in vec2 vTexCoord;
+in vec3 cameraPosition;
+in vec3 normals_commonspace;
+in vec4 position_commonspace;
+in vec4 vColor;
+
+out vec4 fragColor;
+
+void main(void) {
+  geometry.uv = vTexCoord;
+
+  vec3 normal;
+  if (simpleMesh.flatShading) {
+
+  normal = normalize(cross(dFdx(position_commonspace.xyz), dFdy(position_commonspace.xyz)));
+  } else {
+    normal = normals_commonspace;
+  }
+
+  // We initialize color here before passing into DECKGL_FILTER_COLOR
+  vec4 color;
+  DECKGL_FILTER_COLOR(color, geometry);
+
+  vec3 lightColor = lighting_getLightColor(color.rgb, cameraPosition, position_commonspace.xyz, normal);
+  fragColor = vec4(lightColor, color.a * layer.opacity);
 }
+`;
 //#endregion
-//#region node_modules/wkt-parser/transformPROJJSON.js
-function processUnit(unit) {
-	let result = {
-		units: null,
-		to_meter: void 0
-	};
-	if (typeof unit === "string") {
-		result.units = unit.toLowerCase();
-		if (result.units === "metre") result.units = "meter";
-		if (result.units === "meter") result.to_meter = 1;
-	} else if (unit && unit.name) {
-		result.units = unit.name.toLowerCase();
-		if (result.units === "metre") result.units = "meter";
-		result.to_meter = unit.conversion_factor;
+//#region node_modules/@developmentseed/deck.gl-raster/dist/mesh-layer/mesh-layer.js
+var defaultProps$2 = {
+	...SimpleMeshLayer.defaultProps,
+	renderPipeline: {
+		type: "array",
+		value: [],
+		compare: true
+	},
+	material: {
+		ambient: 1,
+		diffuse: 0,
+		shininess: 0,
+		specularColor: [
+			0,
+			0,
+			0
+		]
 	}
-	return result;
-}
-function toValue(valueOrObject) {
-	if (typeof valueOrObject === "object") return valueOrObject.value * valueOrObject.unit.conversion_factor;
-	return valueOrObject;
-}
-function calculateEllipsoid(value, result) {
-	if (value.ellipsoid.radius) {
-		result.a = value.ellipsoid.radius;
-		result.rf = 0;
-	} else {
-		result.a = toValue(value.ellipsoid.semi_major_axis);
-		if (value.ellipsoid.inverse_flattening !== void 0) result.rf = value.ellipsoid.inverse_flattening;
-		else if (value.ellipsoid.semi_major_axis !== void 0 && value.ellipsoid.semi_minor_axis !== void 0) result.rf = result.a / (result.a - toValue(value.ellipsoid.semi_minor_axis));
+};
+/**
+* A small subclass of the SimpleMeshLayer to allow dynamic shader injections.
+*
+* In the future this may expand to diverge more from the SimpleMeshLayer, such
+* as allowing the texture to be a 2D _array_.
+*/
+var MeshTextureLayer = class extends SimpleMeshLayer {
+	static layerName = "mesh-texture-layer";
+	static defaultProps = defaultProps$2;
+	_resolveRenderPipeline() {
+		const { image, renderPipeline } = this.props;
+		return [...image ? [{
+			module: CreateTexture,
+			props: { textureName: image }
+		}] : [], ...renderPipeline ?? []];
 	}
-}
-function transformPROJJSON(projjson, result = {}) {
-	if (!projjson || typeof projjson !== "object") return projjson;
-	if (projjson.type === "BoundCRS") {
-		transformPROJJSON(projjson.source_crs, result);
-		if (projjson.transformation) if (projjson.transformation.method && projjson.transformation.method.name === "NTv2") result.nadgrids = projjson.transformation.parameters[0].value;
-		else result.datum_params = projjson.transformation.parameters.map((param) => param.value);
-		return result;
+	updateState(params) {
+		if (this.hasRenderPipelineChanged(params)) params.changeFlags.extensionsChanged = true;
+		super.updateState(params);
 	}
-	Object.keys(projjson).forEach((key) => {
-		const value = projjson[key];
-		if (value === null) return;
-		switch (key) {
-			case "name":
-				if (result.srsCode) break;
-				result.name = value;
-				result.srsCode = value;
-				break;
-			case "type":
-				if (value === "GeographicCRS") result.projName = "longlat";
-				else if (value === "GeodeticCRS") if (projjson.coordinate_system && projjson.coordinate_system.subtype === "Cartesian") result.projName = "geocent";
-				else result.projName = "longlat";
-				else if (value === "ProjectedCRS" && projjson.conversion && projjson.conversion.method) result.projName = projjson.conversion.method.name;
-				break;
-			case "datum":
-			case "datum_ensemble":
-				if (value.ellipsoid) {
-					result.ellps = value.ellipsoid.name;
-					calculateEllipsoid(value, result);
-				}
-				if (value.prime_meridian) result.from_greenwich = value.prime_meridian.longitude * Math.PI / 180;
-				break;
-			case "ellipsoid":
-				result.ellps = value.name;
-				calculateEllipsoid(value, result);
-				break;
-			case "prime_meridian":
-				result.long0 = (value.longitude || 0) * Math.PI / 180;
-				break;
-			case "coordinate_system":
-				if (value.axis) {
-					const directionMap = {
-						"east": "e",
-						"north": "n",
-						"west": "w",
-						"south": "s",
-						"up": "u",
-						"down": "d",
-						"geocentricx": "e",
-						"geocentricy": "n",
-						"geocentricz": "u"
-					};
-					const mapped = value.axis.map((axis) => directionMap[axis.direction.toLowerCase()]);
-					if (mapped.every(Boolean)) {
-						result.axis = mapped.join("");
-						if (result.axis.length === 2) result.axis += "u";
-					}
-					if (value.unit) {
-						const { units, to_meter } = processUnit(value.unit);
-						result.units = units;
-						result.to_meter = to_meter;
-					} else if (value.axis[0] && value.axis[0].unit) {
-						const { units, to_meter } = processUnit(value.axis[0].unit);
-						result.units = units;
-						result.to_meter = to_meter;
-					}
-				}
-				break;
-			case "id":
-				if (value.authority && value.code) result.title = value.authority + ":" + value.code;
-				break;
-			case "conversion":
-				if (value.method && value.method.name) result.projName = value.method.name;
-				if (value.parameters) value.parameters.forEach((param) => {
-					const paramName = param.name.toLowerCase().replace(/\s+/g, "_");
-					const paramValue = param.value;
-					if (param.unit && param.unit.conversion_factor) result[paramName] = paramValue * param.unit.conversion_factor;
-					else if (param.unit === "degree") result[paramName] = paramValue * Math.PI / 180;
-					else result[paramName] = paramValue;
-				});
-				break;
-			case "unit":
-				if (value.name) {
-					result.units = value.name.toLowerCase();
-					if (result.units === "metre") result.units = "meter";
-				}
-				if (value.conversion_factor) result.to_meter = value.conversion_factor;
-				break;
-			case "base_crs":
-				transformPROJJSON(value, result);
-				result.datumCode = value.id ? value.id.authority + "_" + value.id.code : value.name;
-				break;
-			default: break;
-		}
-	});
-	if (result.latitude_of_false_origin !== void 0) result.lat0 = result.latitude_of_false_origin;
-	if (result.longitude_of_false_origin !== void 0) result.long0 = result.longitude_of_false_origin;
-	if (result.latitude_of_standard_parallel !== void 0) {
-		result.lat0 = result.latitude_of_standard_parallel;
-		result.lat1 = result.latitude_of_standard_parallel;
+	/** Returns true if the render pipeline has changed between the old and new props. */
+	hasRenderPipelineChanged(params) {
+		const { oldProps, props: newProps } = params;
+		if (Boolean(oldProps.image) !== Boolean(newProps.image)) return true;
+		const oldPipeline = oldProps.renderPipeline ?? [];
+		const newPipeline = newProps.renderPipeline ?? [];
+		if (oldPipeline.length !== newPipeline.length) return true;
+		for (let i = 0; i < oldPipeline.length; i++) if (oldPipeline[i]?.module.name !== newPipeline[i]?.module.name) return true;
+		return false;
 	}
-	if (result.latitude_of_1st_standard_parallel !== void 0) result.lat1 = result.latitude_of_1st_standard_parallel;
-	if (result.latitude_of_2nd_standard_parallel !== void 0) result.lat2 = result.latitude_of_2nd_standard_parallel;
-	if (result.latitude_of_projection_centre !== void 0) result.lat0 = result.latitude_of_projection_centre;
-	if (result.longitude_of_projection_centre !== void 0) result.longc = result.longitude_of_projection_centre;
-	if (result.easting_at_false_origin !== void 0) result.x0 = result.easting_at_false_origin;
-	if (result.northing_at_false_origin !== void 0) result.y0 = result.northing_at_false_origin;
-	if (result.latitude_of_natural_origin !== void 0) result.lat0 = result.latitude_of_natural_origin;
-	if (result.longitude_of_natural_origin !== void 0) result.long0 = result.longitude_of_natural_origin;
-	if (result.longitude_of_origin !== void 0) result.long0 = result.longitude_of_origin;
-	if (result.false_easting !== void 0) result.x0 = result.false_easting;
-	if (result.easting_at_projection_centre) result.x0 = result.easting_at_projection_centre;
-	if (result.false_northing !== void 0) result.y0 = result.false_northing;
-	if (result.northing_at_projection_centre) result.y0 = result.northing_at_projection_centre;
-	if (result.standard_parallel_1 !== void 0) result.lat1 = result.standard_parallel_1;
-	if (result.standard_parallel_2 !== void 0) result.lat2 = result.standard_parallel_2;
-	if (result.scale_factor_at_natural_origin !== void 0) result.k0 = result.scale_factor_at_natural_origin;
-	if (result.scale_factor_at_projection_centre !== void 0) result.k0 = result.scale_factor_at_projection_centre;
-	if (result.scale_factor_on_pseudo_standard_parallel !== void 0) result.k0 = result.scale_factor_on_pseudo_standard_parallel;
-	if (result.azimuth !== void 0) result.alpha = result.azimuth;
-	if (result.azimuth_at_projection_centre !== void 0) result.alpha = result.azimuth_at_projection_centre;
-	if (result.angle_from_rectified_to_skew_grid) result.rectified_grid_angle = result.angle_from_rectified_to_skew_grid;
-	applyProjectionDefaults(result);
-	return result;
-}
+	getShaders() {
+		const upstreamShaders = super.getShaders();
+		const modules = upstreamShaders.modules;
+		for (const m of this._resolveRenderPipeline()) modules.push(m.module);
+		return {
+			...upstreamShaders,
+			fs: mesh_layer_fragment_glsl_default,
+			modules
+		};
+	}
+	draw(opts) {
+		const shaderProps = {};
+		for (const m of this._resolveRenderPipeline()) shaderProps[m.module.name] = m.props || {};
+		for (const m of super.getModels()) m.shaderInputs.setProps(shaderProps);
+		super.draw(opts);
+	}
+};
 //#endregion
-//#region node_modules/wkt-parser/index.js
-var knownTypes = [
-	"PROJECTEDCRS",
-	"PROJCRS",
-	"GEOGCS",
-	"GEOCCS",
-	"PROJCS",
-	"LOCAL_CS",
-	"GEODCRS",
-	"GEODETICCRS",
-	"GEODETICDATUM",
-	"ENGCRS",
-	"ENGINEERINGCRS"
+//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-layer.js
+var DEFAULT_MAX_ERROR = .125;
+var DEBUG_COLORS = [
+	[
+		252,
+		73,
+		163
+	],
+	[
+		255,
+		51,
+		204
+	],
+	[
+		204,
+		102,
+		255
+	],
+	[
+		153,
+		51,
+		255
+	],
+	[
+		102,
+		204,
+		255
+	],
+	[
+		51,
+		153,
+		255
+	],
+	[
+		102,
+		255,
+		204
+	],
+	[
+		51,
+		255,
+		170
+	],
+	[
+		0,
+		255,
+		0
+	],
+	[
+		51,
+		204,
+		51
+	],
+	[
+		255,
+		204,
+		102
+	],
+	[
+		255,
+		179,
+		71
+	],
+	[
+		255,
+		102,
+		102
+	],
+	[
+		255,
+		80,
+		80
+	],
+	[
+		255,
+		0,
+		0
+	],
+	[
+		204,
+		0,
+		0
+	],
+	[
+		255,
+		128,
+		0
+	],
+	[
+		255,
+		153,
+		51
+	],
+	[
+		255,
+		255,
+		102
+	],
+	[
+		255,
+		255,
+		51
+	],
+	[
+		0,
+		255,
+		255
+	],
+	[
+		0,
+		204,
+		255
+	]
 ];
-function rename(obj, params) {
-	var outName = params[0];
-	var inName = params[1];
-	if (!(outName in obj) && inName in obj) {
-		obj[outName] = obj[inName];
-		if (params.length === 3) obj[outName] = params[2](obj[outName]);
-	}
-}
-function cleanWKT(wkt) {
-	var keys = Object.keys(wkt);
-	for (var i = 0, ii = keys.length; i < ii; ++i) {
-		var key = keys[i];
-		if (knownTypes.indexOf(key) !== -1) setPropertiesFromWkt(wkt[key]);
-		if (typeof wkt[key] === "object") cleanWKT(wkt[key]);
+var defaultProps$1 = {
+	image: {
+		type: "image",
+		value: null,
+		async: true
+	},
+	renderPipeline: {
+		type: "array",
+		value: [],
+		compare: true
+	},
+	debug: false,
+	debugOpacity: .5
+};
+/**
+* Generic deck.gl layer for rendering geospatial raster data with client-side,
+* GPU-based reprojection and custom processing pipelines.
+*
+* This is a composite layer that uses {@link RasterReprojector} to generate an adaptive mesh
+* that accurately represents the reprojected raster, then renders it using
+* {@link MeshTextureLayer} (a small wrapper around a deck.gl
+* {@link SimpleMeshLayer}).
+*/
+var RasterLayer = class extends CompositeLayer {
+	static layerName = "RasterLayer";
+	static defaultProps = defaultProps$1;
+	initializeState() {
+		this.setState({});
 	}
-}
-function setPropertiesFromWkt(wkt) {
-	if (wkt.AUTHORITY) {
-		var authority = Object.keys(wkt.AUTHORITY)[0];
-		if (authority && authority in wkt.AUTHORITY) wkt.title = authority + ":" + wkt.AUTHORITY[authority];
+	updateState(params) {
+		super.updateState(params);
+		const { props, oldProps, changeFlags } = params;
+		const reprojectionFnsChanged = props.reprojectionFns.forwardTransform !== oldProps.reprojectionFns?.forwardTransform || props.reprojectionFns.inverseTransform !== oldProps.reprojectionFns?.inverseTransform || props.reprojectionFns.forwardReproject !== oldProps.reprojectionFns?.forwardReproject || props.reprojectionFns.inverseReproject !== oldProps.reprojectionFns?.inverseReproject;
+		if (Boolean(changeFlags.dataChanged) || props.width !== oldProps.width || props.height !== oldProps.height || reprojectionFnsChanged || props.maxError !== oldProps.maxError) this._generateMesh();
 	}
-	if (wkt.type === "GEOGCS") wkt.projName = "longlat";
-	else if (wkt.type === "LOCAL_CS") {
-		wkt.projName = "identity";
-		wkt.local = true;
-	} else if (typeof wkt.PROJECTION === "object") wkt.projName = Object.keys(wkt.PROJECTION)[0];
-	else wkt.projName = wkt.PROJECTION;
-	if (wkt.AXIS) {
-		var axisOrder = "";
-		for (var i = 0, ii = wkt.AXIS.length; i < ii; ++i) {
-			var axis = [wkt.AXIS[i][0].toLowerCase(), wkt.AXIS[i][1].toLowerCase()];
-			if (axis[0].indexOf("north") !== -1 || (axis[0] === "y" || axis[0] === "lat") && axis[1] === "north") axisOrder += "n";
-			else if (axis[0].indexOf("south") !== -1 || (axis[0] === "y" || axis[0] === "lat") && axis[1] === "south") axisOrder += "s";
-			else if (axis[0].indexOf("east") !== -1 || (axis[0] === "x" || axis[0] === "lon") && axis[1] === "east") axisOrder += "e";
-			else if (axis[0].indexOf("west") !== -1 || (axis[0] === "x" || axis[0] === "lon") && axis[1] === "west") axisOrder += "w";
-		}
-		if (axisOrder.length === 2) axisOrder += "u";
-		if (axisOrder.length === 3) wkt.axis = axisOrder;
+	_generateMesh() {
+		const { width, height, reprojectionFns, maxError = DEFAULT_MAX_ERROR } = this.props;
+		const reprojector = new RasterReprojector(reprojectionFns, width + 1, height + 1);
+		reprojector.run(maxError);
+		const { indices, positions, texCoords } = reprojectorToMesh(reprojector);
+		this.setState({
+			reprojector,
+			mesh: {
+				indices: {
+					value: indices,
+					size: 1
+				},
+				attributes: {
+					POSITION: {
+						value: positions,
+						size: 3
+					},
+					TEXCOORD_0: {
+						value: texCoords,
+						size: 2
+					}
+				}
+			}
+		});
 	}
-	if (wkt.UNIT) {
-		wkt.units = wkt.UNIT.name.toLowerCase();
-		if (wkt.units === "metre") wkt.units = "meter";
-		if (wkt.UNIT.convert) if (wkt.type === "GEOGCS") {
-			if (wkt.DATUM && wkt.DATUM.SPHEROID) wkt.to_meter = wkt.UNIT.convert * wkt.DATUM.SPHEROID.a;
-		} else wkt.to_meter = wkt.UNIT.convert;
+	renderDebugLayer() {
+		const { reprojector } = this.state;
+		const { debugOpacity } = this.props;
+		if (!reprojector) return null;
+		return new PolygonLayer(this.getSubLayerProps({
+			id: "polygon",
+			data: {
+				reprojector,
+				length: reprojector.triangles.length / 3
+			},
+			getPolygon: (_, { index, data }) => {
+				const triangles = data.reprojector.triangles;
+				const positions = reprojector.exactOutputPositions;
+				const a = triangles[index * 3];
+				const b = triangles[index * 3 + 1];
+				const c = triangles[index * 3 + 2];
+				return [
+					[positions[a * 2], positions[a * 2 + 1]],
+					[positions[b * 2], positions[b * 2 + 1]],
+					[positions[c * 2], positions[c * 2 + 1]],
+					[positions[a * 2], positions[a * 2 + 1]]
+				];
+			},
+			getFillColor: (_, { index, target }) => {
+				const color = DEBUG_COLORS[index % DEBUG_COLORS.length];
+				target[0] = color[0];
+				target[1] = color[1];
+				target[2] = color[2];
+				target[3] = 255;
+				return target;
+			},
+			getLineColor: [
+				0,
+				0,
+				0
+			],
+			getLineWidth: 1,
+			lineWidthUnits: "pixels",
+			opacity: debugOpacity !== void 0 && Number.isFinite(debugOpacity) ? Math.max(0, Math.min(1, debugOpacity)) : 1,
+			pickable: false
+		}));
 	}
-	var geogcs = wkt.GEOGCS;
-	if (wkt.type === "GEOGCS") geogcs = wkt;
-	if (geogcs) {
-		if (geogcs.PRIMEM && geogcs.PRIMEM.convert) wkt.from_greenwich = d2r(geogcs.PRIMEM.convert);
-		if (geogcs.DATUM) wkt.datumCode = geogcs.DATUM.name.toLowerCase();
-		else wkt.datumCode = geogcs.name.toLowerCase();
-		if (wkt.datumCode.slice(0, 2) === "d_") wkt.datumCode = wkt.datumCode.slice(2);
-		if (wkt.datumCode === "new_zealand_1949") wkt.datumCode = "nzgd49";
-		if (wkt.datumCode === "wgs_1984" || wkt.datumCode === "world_geodetic_system_1984") {
-			if (wkt.PROJECTION === "Mercator_Auxiliary_Sphere") wkt.sphere = true;
-			wkt.datumCode = "wgs84";
-		}
-		if (wkt.datumCode === "belge_1972") wkt.datumCode = "rnb72";
-		if (geogcs.DATUM && geogcs.DATUM.SPHEROID) {
-			wkt.ellps = geogcs.DATUM.SPHEROID.name.replace("_19", "").replace(/[Cc]larke\_18/, "clrk");
-			if (wkt.ellps.toLowerCase().slice(0, 13) === "international") wkt.ellps = "intl";
-			wkt.a = geogcs.DATUM.SPHEROID.a;
-			wkt.rf = parseFloat(geogcs.DATUM.SPHEROID.rf);
+	renderLayers() {
+		const { mesh } = this.state;
+		const { debug, image, renderPipeline } = this.props;
+		if (!mesh || !image && (renderPipeline?.length ?? 0) === 0) return null;
+		const layers = [new MeshTextureLayer(this.getSubLayerProps({
+			id: "raster",
+			image,
+			renderPipeline,
+			data: [1],
+			mesh,
+			_instanced: false,
+			getPosition: [
+				0,
+				0,
+				0
+			],
+			getColor: [
+				255,
+				255,
+				255
+			]
+		}))];
+		if (debug) {
+			const debugLayer = this.renderDebugLayer();
+			if (debugLayer) layers.push(debugLayer);
 		}
-		if (geogcs.DATUM && geogcs.DATUM.TOWGS84) wkt.datum_params = geogcs.DATUM.TOWGS84;
-		if (~wkt.datumCode.indexOf("osgb_1936")) wkt.datumCode = "osgb36";
-		if (~wkt.datumCode.indexOf("osni_1952")) wkt.datumCode = "osni52";
-		if (~wkt.datumCode.indexOf("tm65") || ~wkt.datumCode.indexOf("geodetic_datum_of_1965")) wkt.datumCode = "ire65";
-		if (wkt.datumCode === "ch1903+") wkt.datumCode = "ch1903";
-		if (~wkt.datumCode.indexOf("israel")) wkt.datumCode = "isr93";
+		return layers;
 	}
-	if (wkt.b && !isFinite(wkt.b)) wkt.b = wkt.a;
-	if (wkt.rectified_grid_angle) wkt.rectified_grid_angle = d2r(wkt.rectified_grid_angle);
-	function toMeter(input) {
-		return input * (wkt.to_meter || 1);
+};
+function reprojectorToMesh(reprojector) {
+	const numVertices = reprojector.uvs.length / 2;
+	const positions = new Float32Array(numVertices * 3);
+	const texCoords = new Float32Array(reprojector.uvs);
+	for (let i = 0; i < numVertices; i++) {
+		positions[i * 3] = reprojector.exactOutputPositions[i * 2];
+		positions[i * 3 + 1] = reprojector.exactOutputPositions[i * 2 + 1];
+		positions[i * 3 + 2] = 0;
 	}
-	var renamer = function(a) {
-		return rename(wkt, a);
+	return {
+		indices: new Uint32Array(reprojector.triangles),
+		positions,
+		texCoords
 	};
-	[
-		["standard_parallel_1", "Standard_Parallel_1"],
-		["standard_parallel_1", "Latitude of 1st standard parallel"],
-		["standard_parallel_2", "Standard_Parallel_2"],
-		["standard_parallel_2", "Latitude of 2nd standard parallel"],
-		["false_easting", "False_Easting"],
-		["false_easting", "False easting"],
-		["false-easting", "Easting at false origin"],
-		["false_northing", "False_Northing"],
-		["false_northing", "False northing"],
-		["false_northing", "Northing at false origin"],
-		["central_meridian", "Central_Meridian"],
-		["central_meridian", "Longitude of natural origin"],
-		["central_meridian", "Longitude of false origin"],
-		["latitude_of_origin", "Latitude_Of_Origin"],
-		["latitude_of_origin", "Central_Parallel"],
-		["latitude_of_origin", "Latitude of natural origin"],
-		["latitude_of_origin", "Latitude of false origin"],
-		["scale_factor", "Scale_Factor"],
-		["k0", "scale_factor"],
-		["latitude_of_center", "Latitude_Of_Center"],
-		["latitude_of_center", "Latitude_of_center"],
-		[
-			"lat0",
-			"latitude_of_center",
-			d2r
-		],
-		["longitude_of_center", "Longitude_Of_Center"],
-		["longitude_of_center", "Longitude_of_center"],
-		[
-			"longc",
-			"longitude_of_center",
-			d2r
-		],
-		[
-			"x0",
-			"false_easting",
-			toMeter
-		],
-		[
-			"y0",
-			"false_northing",
-			toMeter
-		],
-		[
-			"long0",
-			"central_meridian",
-			d2r
-		],
-		[
-			"lat0",
-			"latitude_of_origin",
-			d2r
-		],
-		[
-			"lat0",
-			"standard_parallel_1",
-			d2r
-		],
-		[
-			"lat1",
-			"standard_parallel_1",
-			d2r
-		],
-		[
-			"lat2",
-			"standard_parallel_2",
-			d2r
-		],
-		["azimuth", "Azimuth"],
-		[
-			"alpha",
-			"azimuth",
-			d2r
-		],
-		["srsCode", "name"]
-	].forEach(renamer);
-	applyProjectionDefaults(wkt);
-}
-function wkt_parser_default(wkt) {
-	if (typeof wkt === "object") return transformPROJJSON(wkt);
-	const version = detectWKTVersion(wkt);
-	var lisp = parser_default(wkt);
-	if (version === "WKT2") return transformPROJJSON(buildPROJJSON(lisp));
-	var type = lisp[0];
-	var obj = {};
-	sExpr(lisp, obj);
-	cleanWKT(obj);
-	return obj[type];
 }
 //#endregion
-//#region node_modules/@developmentseed/proj/dist/parse-wkt.js
+//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-tileset/affine-tileset.js
 /**
-* Parse a WKT string or PROJJSON object into a proj4-compatible projection
-* definition.
+* A {@link RasterTilesetDescriptor} backed by per-level affine transforms.
 *
-* This is a typed wrapper around the `wkt-parser` package.
+* Derives `projectedBounds` from the coarsest level's array. Everything else
+* is passed through from the constructor options.
 */
-function parseWkt(input) {
-	const def = wkt_parser_default(input);
-	if (def.projName === "longlat" && (!def.units || def.units === "unknown")) {
-		def.units = "degree";
-		def.to_meter = void 0;
+var AffineTileset = class {
+	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;
 	}
-	return def;
-}
+};
 //#endregion
-//#region node_modules/@developmentseed/proj/dist/registry.js
+//#region node_modules/@developmentseed/affine/dist/affine.js
 /**
-* A global registry holding parsed projection definitions.
+* Create a translation transform from an offset vector.
+*
+* @param xoff  Translation offset in x direction.
+* @param yoff  Translation offset in y direction.
+*
+* @return Transform that applies the given translation.
 */
-var PROJECTION_REGISTRY = /* @__PURE__ */ new Map();
-async function epsgResolver(epsg) {
-	const key = `EPSG:${epsg}`;
-	const cachedProj = PROJECTION_REGISTRY.get(key);
-	if (cachedProj !== void 0) return cachedProj;
-	const proj = parseWkt(await getProjjson(epsg));
-	PROJECTION_REGISTRY.set(key, proj);
-	return proj;
-}
-/** Query epsg.io for the PROJJSON corresponding to the given EPSG code. */
-async function getProjjson(epsg) {
-	const url = `https://epsg.io/${epsg}.json`;
-	const resp = await fetch(url);
-	if (!resp.ok) throw new Error(`Failed to fetch PROJJSON from ${url}`);
-	return await resp.json();
+function translation(xoff, yoff) {
+	return [
+		1,
+		0,
+		xoff,
+		0,
+		1,
+		yoff
+	];
 }
-//#endregion
-//#region node_modules/@developmentseed/proj/dist/transform-bounds.js
 /**
-* Transform boundary densifying the edges to account for nonlinear
-* transformations along these edges and extracting the outermost bounds.
+* Create a scaling transform from a scalar or vector.
+*
+* You can pass either one or two scaling factors. Passing only a single scalar
+* value will scale in both dimensions equally. A vector scaling value scales
+* the dimensions independently.
+*
+* @param sx  Scaling factor in x direction.
+* @param sy  Scaling factor in y direction (defaults to sx if not provided).
 *
-* @param project - function that maps (x, y) in source CRS to (x, y) in target CRS
-* @param left - min X in source CRS
-* @param bottom - min Y in source CRS
-* @param right - max X in source CRS
-* @param top - max Y in source CRS
-* @param options.densifyPts - number of intermediate points along each edge (default 21)
-* @returns [minX, minY, maxX, maxY] in the target CRS
+* @return Transform that applies the given scaling.
 */
-function transformBounds(project, left, bottom, right, top, options = {}) {
-	const { densifyPts = 21 } = options;
-	const cx = [
-		left,
-		right,
-		right,
-		left
-	];
-	const cy = [
-		bottom,
-		bottom,
-		top,
-		top
-	];
-	let outMinX = Infinity;
-	let outMinY = Infinity;
-	let outMaxX = -Infinity;
-	let outMaxY = -Infinity;
-	for (let i = 0; i < 4; i++) {
-		const fromX = cx[i];
-		const fromY = cy[i];
-		const toX = cx[(i + 1) % 4];
-		const toY = cy[(i + 1) % 4];
-		for (let j = 0; j <= densifyPts; j++) {
-			const t = j / (densifyPts + 1);
-			const [px, py] = project(fromX + (toX - fromX) * t, fromY + (toY - fromY) * t);
-			if (px < outMinX) outMinX = px;
-			if (py < outMinY) outMinY = py;
-			if (px > outMaxX) outMaxX = px;
-			if (py > outMaxY) outMaxY = py;
-		}
-	}
+function scale$3(sx, sy = sx) {
 	return [
-		outMinX,
-		outMinY,
-		outMaxX,
-		outMaxY
+		sx,
+		0,
+		0,
+		0,
+		sy,
+		0
 	];
 }
-//#endregion
-//#region node_modules/@developmentseed/proj/dist/web-mercator.js
-var WGS84_ELLIPSOID_A$1 = 6378137;
-var MAX_WEB_MERCATOR_LAT$1 = 85.05112877980659;
 /**
-* Convert a WGS84 longitude/latitude to EPSG:3857 meters analytically.
-* Valid for latitudes in [-MAX_WEB_MERCATOR_LAT, MAX_WEB_MERCATOR_LAT].
+* Apply a geotransform to a coordinate.
+*
+* That is, we apply this series of equations:
+*
+* ```
+*  x_out = a * x + b * y + c
+*  y_out = d * x + e * y + f
+* ```
+*
+* @param affine  The affine transform to apply.
+* @param x       The x coordinate.
+* @param y       The y coordinate.
+*
+* @return The transformed coordinates.
 */
-function wgs84To3857(lon, lat) {
-	const x = lon * Math.PI * WGS84_ELLIPSOID_A$1 / 180;
-	const latRad = lat * Math.PI / 180;
-	return [x, Math.log(Math.tan(Math.PI / 4 + latRad / 2)) * WGS84_ELLIPSOID_A$1];
+function apply([a, b, c, d, e, f], x, y) {
+	return [a * x + b * y + c, d * x + e * y + f];
 }
 /**
-* Wrap a proj4 forward projection to EPSG:3857 so that it never returns NaN.
+* Compose two affine transforms: A×B (apply B **first**, then A).
 *
-* proj4 returns [NaN, NaN] for points at the poles (lat = ±90°) because the
-* Mercator projection is undefined there. The wrapper falls back to:
-*   1. Project the input to WGS84 via `forwardTo4326`
-*   2. Clamp the latitude to the Web Mercator limit (±85.05°)
-*   3. Convert analytically from WGS84 to EPSG:3857
+* This is equivalent to `a @ b` in Python's `affine` library, and is equivalent
+* to multiplying the 3×3 matrices:
+* ```
+*   | a1 b1 c1 |   | a2 b2 c2 |
+*   | d1 e1 f1 | × | d2 e2 f2 |
+*   | 0  0  1  |   | 0  0  1  |
+* ```
 *
-* This correctly handles any input CRS, not just EPSG:4326.
+* @param A The first affine transform to apply.
+* @param B The second affine transform to apply.
 *
-* NOTE: An identical copy of this function lives in `raster-tile-traversal.ts`.
-* The two packages cannot share code due to their dependency relationship
-* (deck.gl-geotiff depends on deck.gl-raster, not vice versa). If this logic
-* changes, update both copies.
+* @return The composed affine transform.
+*/
+function compose([a1, b1, c1, d1, e1, f1], [a2, b2, c2, d2, e2, f2]) {
+	return [
+		a1 * a2 + b1 * d2,
+		a1 * b2 + b1 * e2,
+		a1 * c2 + b1 * f2 + c1,
+		d1 * a2 + e1 * d2,
+		d1 * b2 + e1 * e2,
+		d1 * c2 + e1 * f2 + f1
+	];
+}
+/**
+* Compute the inverse of an Affine.
 *
-* Perhaps in the future we'll make a `@developmentseed/projections` package to
-* hold shared projection utilities like this. *
+* @param affine  The affine transform to invert.
+* @return The inverted affine transform.
+* @throws If the transform is degenerate and cannot be inverted.
 */
-function makeClampedForwardTo3857$1(forwardTo3857, forwardTo4326) {
-	return (x, y) => {
-		const [px, py] = forwardTo3857(x, y);
-		if (Number.isFinite(px) && Number.isFinite(py)) return [px, py];
-		const [lon, lat] = forwardTo4326(x, y);
-		return wgs84To3857(lon, Math.max(-85.05112877980659, Math.min(MAX_WEB_MERCATOR_LAT$1, lat)));
-	};
+function invert$2([sa, sb, sc, sd, se, sf]) {
+	const det = sa * se - sb * sd;
+	if (det === 0) throw new Error("Cannot invert degenerate transform");
+	const idet = 1 / det;
+	const ra = se * idet;
+	const rb = -sb * idet;
+	const rd = -sd * idet;
+	const re = sa * idet;
+	return [
+		ra,
+		rb,
+		-sc * ra - sf * rb,
+		rd,
+		re,
+		-sc * rd - sf * re
+	];
+}
+/** Get the 'a' component of an Affine transform. */
+function a(affine) {
+	return affine[0];
+}
+/** Get the 'e' component of an Affine transform. */
+function e(affine) {
+	return affine[4];
 }
 //#endregion
+//#region node_modules/@developmentseed/deck.gl-raster/dist/raster-tileset/affine-tileset-level.js
+/**
+* A {@link RasterTilesetLevel} 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 RasterTilesetLevel} manually.
+*/
+var AffineTilesetLevel = class {
+	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 = invert$2(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);
+		const a$1 = a(options.affine);
+		const e$1 = e(options.affine);
+		this.metersPerPixel = Math.sqrt(Math.abs(a$1 * e$1)) * options.mpu;
+		const corners = [
+			apply(options.affine, 0, 0),
+			apply(options.affine, options.arrayWidth, 0),
+			apply(options.affine, 0, options.arrayHeight),
+			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: apply(af, col * tw, row * th),
+			topRight: apply(af, (col + 1) * tw, row * th),
+			bottomLeft: apply(af, col * tw, (row + 1) * th),
+			bottomRight: apply(af, (col + 1) * tw, (row + 1) * th)
+		};
+	}
+	tileTransform(col, row) {
+		const tileOffset = translation(col * this.tileWidth, row * this.tileHeight);
+		const tileAffine = compose(this._affine, tileOffset);
+		const invTileAffine = invert$2(tileAffine);
+		return {
+			forwardTransform: (x, y) => apply(tileAffine, x, y),
+			inverseTransform: (x, y) => apply(invTileAffine, x, y)
+		};
+	}
+	crsBoundsToTileRange(projectedMinX, projectedMinY, projectedMaxX, projectedMaxY) {
+		const inv = this._invAffine;
+		const pixelCorners = [
+			apply(inv, projectedMinX, projectedMinY),
+			apply(inv, projectedMaxX, projectedMinY),
+			apply(inv, projectedMinX, projectedMaxY),
+			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;
+		return {
+			minCol: Math.max(0, Math.floor(pixMinX / tw)),
+			maxCol: Math.min(maxColIdx, Math.floor(pixMaxX / tw)),
+			minRow: Math.max(0, Math.floor(pixMinY / th)),
+			maxRow: Math.min(maxRowIdx, Math.floor(pixMaxY / th))
+		};
+	}
+};
+//#endregion
 //#region node_modules/@developmentseed/deck.gl-raster/dist/raster-tileset/bounding-volume-cache.js
 var DEFAULT_MAX_ENTRIES = 65536;
 /**
@@ -16284,6 +16336,44 @@ var COLORMAP_OPTIONS = COLORMAP_NAMES.map((name) => ({
 	rowIndex: COLORMAP_INDEX[name]
 }));
 //#endregion
+//#region src/lib/raster/repair-geokeys.ts
+/**
+* Repair GeoTIFFs whose projected CRS is tagged with a "user-defined" model
+* type so {@link import('@developmentseed/geotiff').Overview.crs} can build it.
+*
+* `crsFromGeoKeys` only accepts `GTModelTypeGeoKey` of 1 (projected) or 2
+* (geographic) and throws `Unsupported GeoTIFF model type: 32767` for anything
+* else. Some exporters (notably ArcGIS, which writes an ESRI PE string for the
+* `WGS_1984_Web_Mercator` / "Popular Visualisation CRS" auxiliary-sphere
+* Mercator) emit a fully specified projected CRS via the projection geo keys
+* but leave `GTModelTypeGeoKey` and `ProjectedCSTypeGeoKey` as user-defined
+* (32767). COGLayer then fails to resolve the CRS and the raster never draws,
+* even though the projection is completely described by the keys already
+* present. See GeoLibre issue #393.
+*
+* When a projection coordinate-transformation method is present
+* (`ProjMethodGeoKey`, exposed as `gkd.projMethod`), the CRS is unambiguously
+* projected, so we set `modelType` to projected. `crsFromGeoKeys` then takes
+* the projected path and builds a PROJJSON CRS from the projection keys (or
+* returns the EPSG code when `ProjectedCSTypeGeoKey` carries one). Files with a
+* valid model type, or user-defined geographic CRSes with no projection method,
+* are left untouched.
+*
+* Must run before any `Overview.crs` access, which caches its result; callers
+* invoke it immediately after opening the GeoTIFF and before handing it to
+* COGLayer.
+*
+* @param tiff - The opened GeoTIFF to repair in place.
+*/
+function repairUserDefinedProjectedCrs(tiff) {
+	const MODEL_TYPE_PROJECTED = 1;
+	const MODEL_TYPE_USER_DEFINED = 32767;
+	for (const overview of tiff.overviews) {
+		const gkd = overview.gkd;
+		if ((gkd.modelType === null || gkd.modelType === MODEL_TYPE_USER_DEFINED) && gkd.projMethod !== null) gkd.modelType = MODEL_TYPE_PROJECTED;
+	}
+}
+//#endregion
 //#region src/lib/raster/load-geotiff.ts
 /**
 * COG access pattern is hundreds of distinct Range requests against the
@@ -16362,10 +16452,12 @@ function loadGeoTIFF(url) {
 	const promise = (async () => {
 		const source = new CorsSafeSourceHttp(url, {});
 		const view = new SourceView(source, [new SourceChunk({ size: CHUNK_SIZE }), new SourceCache({ size: CACHE_SIZE })]);
-		return await GeoTIFF.open({
+		const tiff = await GeoTIFF.open({
 			dataSource: source,
 			headerSource: view
 		});
+		repairUserDefinedProjectedCrs(tiff);
+		return tiff;
 	})();
 	inflight.set(url, promise);
 	promise.finally(() => inflight.delete(url)).catch(() => {});
@@ -17318,6 +17410,7 @@ function extractPalette(tiff) {
 var DEFAULT_DEPS$1 = {
 	loadGeoTIFF,
 	computeAutoStats,
+	epsgResolver: createResilientEpsgResolver(),
 	createOverlay: (map, options) => {
 		const overlay = new MapboxOverlay({
 			interleaved: options.interleaved,
@@ -17351,6 +17444,7 @@ var LayerManager = class {
 	_device = null;
 	_colormapTexture = null;
 	_handlers = new globalThis.Map();
+	_crsFailed = /* @__PURE__ */ new Set();
 	_destroyed = false;
 	/**
 	* Creates a LayerManager bound to a map.
@@ -17486,6 +17580,7 @@ var LayerManager = class {
 		const [layer] = this._layers.splice(index, 1);
 		layer.abort.abort();
 		if (layer.source.kind === "file") URL.revokeObjectURL(layer.source.objectUrl);
+		this._crsFailed.delete(id);
 		this._destroyPaletteTexture(layer);
 		if (this._selectedId === id) this.select(this._layers[this._layers.length - 1]?.id ?? null);
 		this._rebuild();
@@ -17679,7 +17774,7 @@ var LayerManager = class {
 	* layer's tile cache across rebuilds. */
 	_rebuild() {
 		if (!this._overlay) return;
-		const layers = this._layers.filter((l) => l.geotiff && l.state.visible).map((l) => this._buildCogLayer(l));
+		const layers = this._layers.filter((l) => l.geotiff && l.state.visible && !this._crsFailed.has(l.id)).map((l) => this._buildCogLayer(l));
 		this._overlay.setProps({ layers });
 	}
 	_buildCogLayer(layer) {
@@ -17691,6 +17786,10 @@ var LayerManager = class {
 			getTileData,
 			renderTile,
 			beforeId: this._resolveBeforeId(layer.beforeId),
+			epsgResolver: (epsg) => this._deps.epsgResolver(epsg).catch((err) => {
+				this._failLayerCrs(layer, err);
+				throw err;
+			}),
 			onGeoTIFFLoad: (_tiff, options) => {
 				const boundsArrived = !layer.bounds;
 				layer.bounds = options.geographicBounds;
@@ -17705,6 +17804,26 @@ var LayerManager = class {
 			}
 		});
 	}
+	/** Records a CRS-resolution failure once: marks the layer errored, drops it
+	* from future rebuilds (so the resolver is not retried in a loop), and emits
+	* the error to consumers. Runs after the resolver promise settles, so it is
+	* never synchronous with a deck.gl render. */
+	_failLayerCrs(layer, err) {
+		if (this._destroyed || !this.getLayer(layer.id) || this._crsFailed.has(layer.id)) return;
+		this._crsFailed.add(layer.id);
+		const error = err instanceof Error ? err : new Error(String(err));
+		layer.loading = false;
+		layer.error = error;
+		this._emit({
+			type: "error",
+			layerId: layer.id,
+			error
+		});
+		this._emit({
+			type: "rasterchange",
+			layerId: layer.id
+		});
+	}
 	/** Returns the beforeId only when that layer exists in the map's current
 	* style (warns otherwise). */
 	_resolveBeforeId(beforeId) {
@@ -19151,7 +19270,8 @@ var DEFAULT_OPTIONS = {
 	className: "",
 	interleaved: true,
 	defaultUrl: "",
-	autoLoad: false
+	autoLoad: false,
+	epsgResolver: createResilientEpsgResolver()
 };
 /**
 * A MapLibre GL control for visualizing local and remote raster datasets
@@ -19212,7 +19332,7 @@ var RasterControl = class {
 		this._container = this._createContainer();
 		this._panel = this._createPanel();
 		this._mapContainer.appendChild(this._panel);
-		this._layerManager = new LayerManager(map, { interleaved: this._options.interleaved });
+		this._layerManager = new LayerManager(map, { interleaved: this._options.interleaved }, { epsgResolver: this._options.epsgResolver });
 		this._forwardLayerManagerEvents(this._layerManager);
 		const manager = this._layerManager;
 		this._inspector = new PixelInspector(map, () => {
@@ -19603,6 +19723,6 @@ var RasterControl = class {
 	}
 };
 //#endregion
-export { colormaps_default as _, debounce as a, throttle as c, percentileFromHistogram as d, readBandNames as f, COLORMAP_ROW_COUNT as g, COLORMAP_OPTIONS as h, classNames as i, MAX_SAMPLE_TILES as l, COLORMAP_NAMES as m, PALETTE_COLORMAP as n, formatNumericValue as o, loadGeoTIFF as p, clamp as r, generateId as s, RasterControl as t, computeAutoStats as u };
+export { colormaps_default as _, debounce as a, throttle as c, percentileFromHistogram as d, readBandNames as f, COLORMAP_ROW_COUNT as g, COLORMAP_OPTIONS as h, classNames as i, MAX_SAMPLE_TILES as l, COLORMAP_NAMES as m, PALETTE_COLORMAP as n, formatNumericValue as o, loadGeoTIFF as p, clamp as r, generateId as s, RasterControl as t, computeAutoStats as u, createResilientEpsgResolver as v };
 
-//# sourceMappingURL=RasterControl-Dsq7VIl7.js.map
+//# sourceMappingURL=RasterControl-BxGg08KC.js.map