maplibre-gl-raster 0.2.1 → 0.2.3

package/dist/{RasterControl-DNAVxbPt.js → RasterControl-COXL7MI-.js}

@@ -6578,9 +6578,9 @@ var SourceMemory = class SourceMemory {
 //#endregion
 //#region node_modules/@developmentseed/geotiff/dist/crs.js
 var PROJJSON_SCHEMA = "https://proj.org/schemas/v0.7/projjson.schema.json";
-var MODEL_TYPE_PROJECTED = 1;
+var MODEL_TYPE_PROJECTED$1 = 1;
 var MODEL_TYPE_GEOGRAPHIC = 2;
-var USER_DEFINED = 32767;
+var USER_DEFINED$1 = 32767;
 var CT_TRANSVERSE_MERCATOR = 1;
 var CT_TRANSVERSE_MERCATOR_SOUTH = 2;
 var CT_OBLIQUE_MERCATOR = 3;
@@ -6633,31 +6633,31 @@ var LINEAR_UNIT = {
 */
 function crsFromGeoKeys(gkd) {
 	const modelType = gkd.modelType;
-	if (modelType === MODEL_TYPE_PROJECTED) return _projectedCrs(gkd);
+	if (modelType === MODEL_TYPE_PROJECTED$1) return _projectedCrs(gkd);
 	if (modelType === MODEL_TYPE_GEOGRAPHIC) return _geographicCrs(gkd);
 	throw new Error(`Unsupported GeoTIFF model type: ${modelType}`);
 }
 function _geographicCrs(gkd) {
 	const epsg = gkd.geodeticCRS;
-	if (epsg !== null && epsg !== USER_DEFINED) return epsg;
+	if (epsg !== null && epsg !== USER_DEFINED$1) return epsg;
 	return _buildGeographicCrs(gkd, PROJJSON_SCHEMA);
 }
 function _projectedCrs(gkd) {
 	const epsg = gkd.projectedCRS;
-	if (epsg !== null && epsg !== USER_DEFINED) return epsg;
+	if (epsg !== null && epsg !== USER_DEFINED$1) return epsg;
 	return _buildProjectedCrs(gkd);
 }
 function _buildGeographicCrs(gkd, schema) {
 	const ellipsoid = _buildEllipsoid(gkd);
 	let pmName = "Greenwich";
 	let pmLongitude = 0;
-	if (gkd.primeMeridian !== null && gkd.primeMeridian !== USER_DEFINED) pmName = `EPSG:${gkd.primeMeridian}`;
+	if (gkd.primeMeridian !== null && gkd.primeMeridian !== USER_DEFINED$1) pmName = `EPSG:${gkd.primeMeridian}`;
 	else if (gkd.primeMeridianLongitude !== null) {
 		pmLongitude = gkd.primeMeridianLongitude;
 		pmName = "User-defined";
 	}
 	let datum;
-	if (gkd.geodeticDatum !== null && gkd.geodeticDatum !== USER_DEFINED) datum = {
+	if (gkd.geodeticDatum !== null && gkd.geodeticDatum !== USER_DEFINED$1) datum = {
 		type: "GeodeticReferenceFrame",
 		name: `Unknown datum based upon EPSG ${gkd.geodeticDatum} ellipsoid`
 	};
@@ -6693,7 +6693,7 @@ function _buildProjectedCrs(gkd) {
 	};
 }
 function _buildEllipsoid(gkd) {
-	if (gkd.ellipsoid !== null && gkd.ellipsoid !== USER_DEFINED) {
+	if (gkd.ellipsoid !== null && gkd.ellipsoid !== USER_DEFINED$1) {
 		const ellipsoid = { name: `EPSG ellipsoid ${gkd.ellipsoid}` };
 		if (gkd.ellipsoidSemiMajorAxis !== null) ellipsoid.semi_major_axis = gkd.ellipsoidSemiMajorAxis;
 		if (gkd.ellipsoidInvFlattening !== null) ellipsoid.inverse_flattening = gkd.ellipsoidInvFlattening;
@@ -16336,6 +16336,63 @@ var COLORMAP_OPTIONS = COLORMAP_NAMES.map((name) => ({
 	rowIndex: COLORMAP_INDEX[name]
 }));
 //#endregion
+//#region src/lib/raster/repair-geokeys.ts
+var MODEL_TYPE_PROJECTED = 1;
+var USER_DEFINED = 32767;
+var EPSG_WEB_MERCATOR = 3857;
+var WEB_MERCATOR_ALIASES = new Set([
+	3785,
+	900913,
+	102100,
+	102113
+]);
+var WEB_MERCATOR_CITATION = /web[ _]?mercator|pseudo[ _-]?mercator|popular[ _]?vis|auxiliary[ _]?sphere/i;
+function citationText(gkd) {
+	return [
+		gkd.citation,
+		gkd.projectedCitation,
+		gkd.geodeticCitation
+	].filter((value) => typeof value === "string").join(" ");
+}
+/**
+* Repair GeoTIFF CRS geo keys so {@link import('@developmentseed/geotiff').Overview.crs}
+* resolves them to the right projection, in place. Two problems are fixed:
+*
+* 1. **User-defined model type for a projected CRS.** `crsFromGeoKeys` only
+*    accepts `GTModelTypeGeoKey` 1 (projected) or 2 (geographic) and throws
+*    `Unsupported GeoTIFF model type: 32767` otherwise, so the raster never
+*    draws. ArcGIS exports the ESRI Web Mercator with a fully specified
+*    projection but leaves `GTModelTypeGeoKey` (and `ProjectedCSTypeGeoKey`)
+*    user-defined. When a projection method is present (`ProjMethodGeoKey`, i.e.
+*    `gkd.projMethod`) the CRS is unambiguously projected, so the model type is
+*    set to projected and the CRS is built from the keys (or the EPSG code in
+*    `ProjectedCSTypeGeoKey`). This is projection-agnostic.
+*
+* 2. **Spherical Web Mercator written as an ellipsoidal CRS.** ESRI's "Web
+*    Mercator Auxiliary Sphere" (a.k.a. Popular Visualisation CRS, EPSG:3857)
+*    projects on a sphere of radius 6378137, but the geo keys can only declare
+*    the WGS84 ellipsoid (the sphere lives only in the citation). Building the
+*    CRS from those keys yields an *ellipsoidal* Mercator, which misplaces the
+*    raster by kilometres in latitude. When the projected CRS is a Web Mercator
+*    alias code, or is user-defined with a Web-Mercator citation, it is replaced
+*    with EPSG:3857 so proj4's correct spherical definition is used.
+*
+* Files with a valid model type and a non-Web-Mercator EPSG/projection are left
+* untouched. Must run before any `Overview.crs` access (which caches its
+* result), so callers invoke it immediately after opening the GeoTIFF.
+*
+* @param tiff - The opened GeoTIFF to repair in place.
+*/
+function repairUserDefinedProjectedCrs(tiff) {
+	for (const overview of tiff.overviews) {
+		const gkd = overview.gkd;
+		const projectedIsUserDefined = gkd.projectedCRS === null || gkd.projectedCRS === USER_DEFINED;
+		if (gkd.projectedCRS !== null && WEB_MERCATOR_ALIASES.has(gkd.projectedCRS)) gkd.projectedCRS = EPSG_WEB_MERCATOR;
+		else if (projectedIsUserDefined && gkd.projMethod !== null && WEB_MERCATOR_CITATION.test(citationText(gkd))) gkd.projectedCRS = EPSG_WEB_MERCATOR;
+		if ((gkd.modelType === null || gkd.modelType === 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
@@ -16414,10 +16471,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(() => {});
@@ -19685,4 +19744,4 @@ 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, createResilientEpsgResolver as v };
 
-//# sourceMappingURL=RasterControl-DNAVxbPt.js.map
+//# sourceMappingURL=RasterControl-COXL7MI-.js.map