@loaders.gl/tiles 5.0.0-alpha.8 → 5.0.0-alpha.9

package/dist/dist.dev.js

@@ -3162,6 +3162,88 @@ var __exports__ = (() => {
     }
   };
 
+  // ../../node_modules/@math.gl/core/dist/spheroid.js
+  function spheroidToCartesian(point, geometry, result = point) {
+    const longitude = point.x, latitude = point.y, height = point.z;
+    const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
+    const sin = Math.sin(latitude), cos = Math.cos(latitude);
+    const ratio = b / a;
+    const radius = es < 0.9 ? a / Math.sqrt(1 - es * sin * sin) : a / Math.hypot(cos, ratio * sin);
+    const x = (radius + height) * cos * Math.cos(longitude);
+    const y = (radius + height) * cos * Math.sin(longitude);
+    const z = (radius * (es < 0.9 ? 1 - es : ratio * ratio) + height) * sin;
+    return commit(result, x, y, z);
+  }
+  function cartesianToSpheroid(point, geometry, result = point) {
+    const x = point.x, y = point.y, z = point.z;
+    const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = geometry;
+    const p = Math.hypot(x, y), rr = Math.hypot(x, y, z);
+    if (!(rr > 0) || !Number.isFinite(rr))
+      return false;
+    if (es === 0)
+      return commit(result, p === 0 ? 0 : Math.atan2(y, x), Math.atan2(z, p), rr - a);
+    if (p === 0)
+      return commit(result, 0, Math.sign(z) * Math.PI / 2, Math.abs(z) - b);
+    if (es >= 0.9)
+      return exteriorInverse(x, y, z, p, a, b, result);
+    const ct = z / rr, st = p / rr;
+    let rx = 1 / Math.sqrt(1 - es * (2 - es) * st * st);
+    let cos = st * (1 - es) * rx, sin = ct * rx;
+    for (let i = 0; i < 30; i++) {
+      const rn = a / Math.sqrt(1 - es * sin * sin);
+      const height = p * cos + z * sin - rn * (1 - es * sin * sin);
+      const rk = es * rn / (rn + height);
+      rx = 1 / Math.sqrt(1 - rk * (2 - rk) * st * st);
+      const nextCos = st * (1 - rk) * rx, nextSin = ct * rx;
+      if (Math.abs(nextSin * cos - nextCos * sin) <= 1e-12) {
+        if (!Number.isFinite(height) || !Number.isFinite(nextSin) || !Number.isFinite(nextCos))
+          break;
+        return commit(result, Math.atan2(y, x), Math.atan2(nextSin, Math.abs(nextCos)), height);
+      }
+      cos = nextCos;
+      sin = nextSin;
+    }
+    return exteriorInverse(x, y, z, p, a, b, result);
+  }
+  function exteriorInverse(x, y, z, p, a, b, result) {
+    const ratio = b / a, ratio2 = ratio * ratio, horizontal = p / a, vertical = z / a;
+    const norm = Math.hypot(horizontal, z / b);
+    if (!Number.isFinite(norm) || norm < 1 - 8 * Number.EPSILON || !(ratio2 > 0))
+      return false;
+    let lower = 0, upper = Math.hypot(horizontal, vertical * ratio);
+    if (!Number.isFinite(upper))
+      return false;
+    let multiplier = norm <= 1 + 8 * Number.EPSILON ? 0 : upper / 2;
+    for (let i = 0; i < 96; i++) {
+      const equatorial = 1 + multiplier, polar = ratio2 + multiplier;
+      const u = horizontal / equatorial, v = vertical * ratio / polar;
+      const length4 = Math.hypot(u, v);
+      const residual = length4 - 1;
+      if (Math.abs(residual) <= 4 * Number.EPSILON) {
+        const latitude = Math.atan2(vertical / polar, horizontal / equatorial);
+        const sin = Math.sin(latitude), cos = Math.cos(latitude);
+        const height = p * cos + z * sin - a * Math.hypot(cos, ratio * sin);
+        return commit(result, Math.atan2(y, x), latitude, height);
+      }
+      if (residual > 0)
+        lower = multiplier;
+      else
+        upper = multiplier;
+      const derivative = -(u * u / equatorial + v * v / polar) / length4;
+      const next = multiplier - residual / derivative;
+      multiplier = Number.isFinite(next) && next > lower && next < upper ? next : (lower + upper) / 2;
+    }
+    return false;
+  }
+  function commit(result, x, y, z) {
+    if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))
+      return false;
+    result.x = x;
+    result.y = y;
+    result.z = z;
+    return true;
+  }
+
   // ../../node_modules/@math.gl/geospatial/dist/constants.js
   var WGS84_RADIUS_X = 6378137;
   var WGS84_RADIUS_Y = 6378137;
@@ -3208,25 +3290,24 @@ var __exports__ = (() => {
   function fromCartographicToRadians(cartographic, vector = []) {
     return fromCartographic(cartographic, vector, config._cartographicRadians ? identity3 : toRadians);
   }
-  function toCartographic(vector, cartographic, map2 = identity3) {
+  function toCartographicFromRadiansComponents(longitude, latitude, height, cartographic) {
+    const map2 = config._cartographicRadians ? identity3 : toDegrees;
+    const mappedLongitude = map2(longitude), mappedLatitude = map2(latitude);
     if ("longitude" in cartographic) {
-      cartographic.longitude = map2(vector[0]);
-      cartographic.latitude = map2(vector[1]);
-      cartographic.height = vector[2];
+      cartographic.longitude = mappedLongitude;
+      cartographic.latitude = mappedLatitude;
+      cartographic.height = height;
     } else if ("x" in cartographic) {
-      cartographic.x = map2(vector[0]);
-      cartographic.y = map2(vector[1]);
-      cartographic.z = vector[2];
+      cartographic.x = mappedLongitude;
+      cartographic.y = mappedLatitude;
+      cartographic.z = height;
     } else {
-      cartographic[0] = map2(vector[0]);
-      cartographic[1] = map2(vector[1]);
-      cartographic[2] = vector[2];
+      cartographic[0] = mappedLongitude;
+      cartographic[1] = mappedLatitude;
+      cartographic[2] = height;
     }
     return cartographic;
   }
-  function toCartographicFromRadians(vector, cartographic) {
-    return toCartographic(vector, cartographic, config._cartographicRadians ? identity3 : toDegrees);
-  }
 
   // ../../node_modules/@math.gl/geospatial/dist/ellipsoid-helpers/ellipsoid-transform.js
   var EPSILON142 = 1e-14;
@@ -3348,10 +3429,14 @@ var __exports__ = (() => {
   var scaleToGeodeticSurfaceGradient = new Vector3();
   function scaleToGeodeticSurface(cartesian, ellipsoid, result = []) {
     const { oneOverRadii, oneOverRadiiSquared, centerToleranceSquared } = ellipsoid;
-    scratchVector4.from(cartesian);
-    const positionX = scratchVector4.x;
-    const positionY = scratchVector4.y;
-    const positionZ = scratchVector4.z;
+    const object2 = cartesian;
+    const positionX = "x" in cartesian ? object2.x : cartesian[0];
+    const positionY = "x" in cartesian ? object2.y : cartesian[1];
+    const positionZ = "x" in cartesian ? object2.z : cartesian[2];
+    if (!Number.isFinite(positionX) || !Number.isFinite(positionY) || !Number.isFinite(positionZ)) {
+      return void 0;
+    }
+    scratchVector4.set(positionX, positionY, positionZ);
     const oneOverRadiiX = oneOverRadii.x;
     const oneOverRadiiY = oneOverRadii.y;
     const oneOverRadiiZ = oneOverRadii.z;
@@ -3360,13 +3445,13 @@ var __exports__ = (() => {
     const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;
     const squaredNorm = x2 + y2 + z2;
     const ratio = Math.sqrt(1 / squaredNorm);
-    if (!Number.isFinite(ratio)) {
+    if (!Number.isFinite(ratio) || !Number.isFinite(squaredNorm)) {
       return void 0;
     }
     const intersection = scaleToGeodeticSurfaceIntersection;
-    intersection.copy(cartesian).scale(ratio);
+    intersection.set(positionX, positionY, positionZ).scale(ratio);
     if (squaredNorm < centerToleranceSquared) {
-      return intersection.to(result);
+      return writeResult(result, intersection.x, intersection.y, intersection.z);
     }
     const oneOverRadiiSquaredX = oneOverRadiiSquared.x;
     const oneOverRadiiSquaredY = oneOverRadiiSquared.y;
@@ -3379,7 +3464,7 @@ var __exports__ = (() => {
     let yMultiplier;
     let zMultiplier;
     let func;
-    do {
+    for (let iteration = 0; iteration < 64; iteration++) {
       lambda -= correction;
       xMultiplier = 1 / (1 + lambda * oneOverRadiiSquaredX);
       yMultiplier = 1 / (1 + lambda * oneOverRadiiSquaredY);
@@ -3391,21 +3476,52 @@ var __exports__ = (() => {
       const yMultiplier3 = yMultiplier2 * yMultiplier;
       const zMultiplier3 = zMultiplier2 * zMultiplier;
       func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1;
+      if (!Number.isFinite(func))
+        return void 0;
+      if (Math.abs(func) <= math_utils_exports.EPSILON12) {
+        return writeResult(result, positionX * xMultiplier, positionY * yMultiplier, positionZ * zMultiplier);
+      }
       const denominator = x2 * xMultiplier3 * oneOverRadiiSquaredX + y2 * yMultiplier3 * oneOverRadiiSquaredY + z2 * zMultiplier3 * oneOverRadiiSquaredZ;
       const derivative = -2 * denominator;
+      if (!Number.isFinite(derivative) || derivative === 0)
+        return void 0;
       correction = func / derivative;
-    } while (Math.abs(func) > math_utils_exports.EPSILON12);
-    return scratchVector4.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);
+      if (!Number.isFinite(correction))
+        return void 0;
+    }
+    return void 0;
+  }
+  function writeResult(result, x, y, z) {
+    if (isArray(result)) {
+      result[0] = x;
+      result[1] = y;
+      result[2] = z;
+    } else {
+      const point = result;
+      point.x = x;
+      point.y = y;
+      point.z = z;
+    }
+    return result;
   }
 
   // ../../node_modules/@math.gl/geospatial/dist/ellipsoid.js
+  var spheroidScratch = { x: 0, y: 0, z: 0 };
   var scratchVector5 = new Vector3();
   var scratchNormal = new Vector3();
   var scratchK = new Vector3();
   var scratchPosition = new Vector3();
   var scratchHeight = new Vector3();
   var scratchCartesian = new Vector3();
-  var Ellipsoid = class {
+  var Ellipsoid = class _Ellipsoid {
+    /** Construct from two-axis geometry, including normalized projection ellipsoids.
+     * The equatorial X/Y radii are equal; no CRS or datum information is inferred.
+     */
+    static fromSpheroid(parameters) {
+      const { semiMajorAxis, semiMinorAxis } = parameters;
+      validateSpheroid(semiMajorAxis, semiMinorAxis);
+      return new _Ellipsoid(semiMajorAxis, semiMajorAxis, semiMinorAxis);
+    }
     constructor(x = 0, y = 0, z = 0) {
       this.centerToleranceSquared = math_utils_exports.EPSILON1;
       assert(x >= 0);
@@ -3421,8 +3537,25 @@ var __exports__ = (() => {
       if (this.radiiSquared.z !== 0) {
         this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;
       }
+      if (x === y && z > 0 && z <= x && Number.isFinite(this.radiiSquared.x) && this.radiiSquared.z > 0)
+        this.spheroidGeometry = Object.freeze({
+          semiMajorAxis: x,
+          semiMinorAxis: z,
+          eccentricitySquared: 1 - this.radiiSquared.z / this.radiiSquared.x
+        });
       Object.freeze(this);
     }
+    /** Return an owned geometry snapshot for a sphere or oblate spheroid.
+     * Triaxial, prolate and degenerate ellipsoids cannot be represented by this contract.
+     * Call during setup, outside coordinate loops.
+     */
+    toSpheroid() {
+      const { x, y, z } = this.radii;
+      if (x !== y)
+        throw new Error("Spheroid requires equal X/Y radii");
+      validateSpheroid(x, z);
+      return Object.freeze({ semiMajorAxis: x, semiMinorAxis: z });
+    }
     /** Compares this Ellipsoid against the provided Ellipsoid componentwise */
     equals(right) {
       return this === right || Boolean(right && this.radii.equals(right.radii));
@@ -3432,6 +3565,15 @@ var __exports__ = (() => {
       return this.radii.toString();
     }
     cartographicToCartesian(cartographic, result = [0, 0, 0]) {
+      if (this.spheroidGeometry) {
+        const llh = fromCartographicToRadians(cartographic, scratchCartesian);
+        spheroidScratch.x = llh[0];
+        spheroidScratch.y = llh[1];
+        spheroidScratch.z = llh[2];
+        if (!spheroidToCartesian(spheroidScratch, this.spheroidGeometry))
+          return void 0;
+        return writeResult(result, spheroidScratch.x, spheroidScratch.y, spheroidScratch.z);
+      }
       const normal = scratchNormal;
       const k = scratchK;
       const [, , height] = cartographic;
@@ -3444,7 +3586,24 @@ var __exports__ = (() => {
       return k.to(result);
     }
     cartesianToCartographic(cartesian, result = [0, 0, 0]) {
-      scratchCartesian.from(cartesian);
+      const object2 = cartesian;
+      const x = "x" in cartesian ? object2.x : cartesian[0];
+      const y = "x" in cartesian ? object2.y : cartesian[1];
+      const z = "x" in cartesian ? object2.z : cartesian[2];
+      if (!Number.isFinite(x) || !Number.isFinite(y) || !Number.isFinite(z))
+        return void 0;
+      const a = this.radii.x, b = this.radii.z;
+      const normSquared = (x / a) ** 2 + (y / a) ** 2 + (z / b) ** 2;
+      if (this.spheroidGeometry && Number.isFinite(normSquared) && normSquared >= 1 - 16 * Number.EPSILON && Number.isFinite(x * x + y * y + z * z)) {
+        spheroidScratch.x = x;
+        spheroidScratch.y = y;
+        spheroidScratch.z = z;
+        if (!cartesianToSpheroid(spheroidScratch, this.spheroidGeometry))
+          return void 0;
+        const longitude2 = x === 0 && y === 0 ? Math.atan2(y, x) : spheroidScratch.x;
+        return toCartographicFromRadiansComponents(longitude2, spheroidScratch.y, spheroidScratch.z, result);
+      }
+      scratchCartesian.set(x, y, z);
       const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);
       if (!point) {
         return void 0;
@@ -3453,9 +3612,9 @@ var __exports__ = (() => {
       const h = scratchHeight;
       h.copy(scratchCartesian).subtract(point);
       const longitude = Math.atan2(normal.y, normal.x);
-      const latitude = Math.asin(normal.z);
+      const latitude = Math.atan2(normal.z, Math.hypot(normal.x, normal.y));
       const height = Math.sign(dot(h, scratchCartesian)) * length(h);
-      return toCartographicFromRadians([longitude, latitude, height], result);
+      return toCartographicFromRadiansComponents(longitude, latitude, height, result);
     }
     eastNorthUpToFixedFrame(origin, result = new Matrix4()) {
       return localFrameToFixedFrame(this, "east", "north", "up", origin, result);
@@ -3519,6 +3678,10 @@ var __exports__ = (() => {
     }
   };
   Ellipsoid.WGS84 = new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);
+  function validateSpheroid(a, b) {
+    if (!Number.isFinite(a) || !Number.isFinite(b) || a <= 0 || b <= 0 || b > a)
+      throw new Error("Spheroid requires finite positive axes with semiMinorAxis <= semiMajorAxis");
+  }
 
   // ../../node_modules/@math.gl/culling/dist/constants.js
   var INTERSECTION;
@@ -11088,8 +11251,8 @@ var __exports__ = (() => {
     return void 0;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projection-descriptor.js
-  var CACHE_KEY = /* @__PURE__ */ Symbol.for("@math.gl/proj4/projection-implementations");
+  // ../../node_modules/@math.gl/projection/dist/experimental/projection-descriptor.js
+  var CACHE_KEY = /* @__PURE__ */ Symbol.for("@math.gl/projection/projection-implementations");
   function cache() {
     const scope = globalThis;
     scope[CACHE_KEY] || (scope[CACHE_KEY] = /* @__PURE__ */ new WeakMap());
@@ -11123,18 +11286,45 @@ var __exports__ = (() => {
     return entry.pending;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/mutable-projection.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projection-scratch.js
+  var ProjectionScratch = class {
+    constructor() {
+      this.point = { x: 0, y: 0, z: 0 };
+      this.busy = false;
+    }
+    acquire() {
+      if (this.busy)
+        return { x: 0, y: 0, z: 0 };
+      this.busy = true;
+      return this.point;
+    }
+    release(point) {
+      if (point === this.point)
+        this.busy = false;
+    }
+  };
+
+  // ../../node_modules/@math.gl/projection/dist/experimental/mutable-projection.js
   function createProjection(forward33, inverse33) {
+    const scratch = new ProjectionScratch();
+    const scalar = (x, y, operation2) => {
+      const point = scratch.acquire();
+      try {
+        point.x = x;
+        point.y = y;
+        point.z = 0;
+        operation2(point);
+        return [point.x, point.y];
+      } finally {
+        scratch.release(point);
+      }
+    };
     return {
       forward(x, y) {
-        const point = { x, y, z: 0 };
-        forward33(point);
-        return [point.x, point.y];
+        return scalar(x, y, forward33);
       },
       inverse(x, y) {
-        const point = { x, y, z: 0 };
-        inverse33(point);
-        return [point.x, point.y];
+        return scalar(x, y, inverse33);
       },
       forwardInPlace: forward33,
       inverseInPlace: inverse33
@@ -11161,9 +11351,43 @@ var __exports__ = (() => {
     };
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/types.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/scalar-output.js
+  var FLOAT32_MAX = 34028234663852886e22;
+  function validateScalarOutput(coordinate, output, length4) {
+    if (!(Array.isArray(coordinate) || coordinate instanceof Float32Array || coordinate instanceof Float64Array))
+      throw new Error("Scalar input requires a number array, Float32Array or Float64Array");
+    if (!(Array.isArray(output) || output instanceof Float32Array || output instanceof Float64Array) || output.length < length4)
+      throw new Error("Scalar output requires a preallocated number array, Float32Array or Float64Array with sufficient capacity");
+    if (coordinate === output || Array.isArray(coordinate) || Array.isArray(output))
+      return;
+    const input = coordinate;
+    const result = output;
+    const shared = typeof SharedArrayBuffer !== "undefined" && input.buffer instanceof SharedArrayBuffer && result.buffer instanceof SharedArrayBuffer;
+    if ((input.buffer === result.buffer || shared) && input.byteOffset < result.byteOffset + length4 * result.BYTES_PER_ELEMENT && result.byteOffset < input.byteOffset + input.byteLength)
+      throw new Error("Scalar input/output views overlap; use the same view for in-place projection");
+  }
+  function writeScalarOutput(coordinate, output, point, length4) {
+    if (output instanceof Float32Array) {
+      if (Math.abs(point.x) > FLOAT32_MAX || Math.abs(point.y) > FLOAT32_MAX || length4 >= 3 && Math.abs(point.z) > FLOAT32_MAX)
+        throw new Error("Projected coordinate exceeds Float32 range");
+      for (let i = 3; i < coordinate.length; i++) {
+        const value2 = coordinate[i];
+        if (Number.isFinite(value2) && Math.abs(value2) > FLOAT32_MAX)
+          throw new Error("Trailing ordinate exceeds Float32 range");
+      }
+    }
+    for (let i = 3; i < coordinate.length; i++)
+      output[i] = coordinate[i];
+    output[0] = point.x;
+    output[1] = point.y;
+    if (length4 >= 3)
+      output[2] = point.z;
+    return output;
+  }
+
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/types.js
   var _a;
-  var CRS_ERROR = /* @__PURE__ */ Symbol.for("@math.gl/proj4/TypeScriptCRSError");
+  var CRS_ERROR = /* @__PURE__ */ Symbol.for("@math.gl/projection/TypeScriptCRSError");
   var TypeScriptCRSError = class _TypeScriptCRSError extends Error {
     static [(_a = CRS_ERROR, Symbol.hasInstance)](value2) {
       if (this !== _TypeScriptCRSError)
@@ -11181,7 +11405,7 @@ var __exports__ = (() => {
     throw new TypeScriptCRSError("missing-transform-stage", message);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/spatial-reference.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/spatial-reference.js
   function resolveCRSInput(input, options) {
     if (typeof input === "string" || "type" in input)
       return { definition: input, lossy: false };
@@ -11255,7 +11479,7 @@ var __exports__ = (() => {
     return Object.freeze({ ...crs, storedAxis });
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/parameters.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/parameters.js
   var DEGREES_TO_RADIANS2 = Math.PI / 180;
   function numberParameter(parameters, name, fallback) {
     if (!Object.prototype.hasOwnProperty.call(parameters, name))
@@ -11282,7 +11506,7 @@ var __exports__ = (() => {
     return ((longitude + Math.PI) % circle + circle) % circle - Math.PI;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/datum-table.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/datum-table.js
   var datums = {
     wgs84: {
       towgs84: "0,0,0",
@@ -12644,7 +12868,7 @@ var __exports__ = (() => {
   }
   var datum_table_default = datums;
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/ellipsoid-table.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/ellipsoid-table.js
   var ellipsoids = {
     MERIT: {
       a: 6378137,
@@ -12870,7 +13094,7 @@ var __exports__ = (() => {
   };
   var ellipsoid_table_default = ellipsoids;
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/primemeridian-table.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/primemeridian-table.js
   var primeMeridian = {};
   primeMeridian["greenwich"] = 0;
   primeMeridian["lisbon"] = -9.131906111111;
@@ -12887,7 +13111,7 @@ var __exports__ = (() => {
   primeMeridian["oslo"] = 10.722916666667;
   var primemeridian_table_default = primeMeridian;
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/units-table.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/units-table.js
   var units_table_default = {
     mm: { to_meter: 1e-3 },
     cm: { to_meter: 0.01 },
@@ -12911,7 +13135,7 @@ var __exports__ = (() => {
     "us-yd": { to_meter: 0.914401828803658 }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/normalize.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/normalize.js
   var ALIASES = {
     WGS84: "+proj=longlat +datum=WGS84",
     "EPSG:4326": "+proj=longlat +datum=WGS84",
@@ -13142,7 +13366,7 @@ var __exports__ = (() => {
     });
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/datum.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/datum.js
   var ARC_SECOND = Math.PI / (180 * 3600);
   var WGS84 = {
     ellipsoid: {
@@ -13153,47 +13377,29 @@ var __exports__ = (() => {
     towgs84: [0, 0, 0]
   };
   function geodeticToGeocentricInPlace(point, ellipsoid) {
-    const lon = point.x, lat = point.y, height = point.z;
-    const { semiMajorAxis: a, eccentricitySquared: es } = ellipsoid;
-    if (Math.abs(lat) > Math.PI / 2)
+    if (Math.abs(point.y) > Math.PI / 2)
       throw new Error("Latitude outside geocentric domain");
-    const sin = Math.sin(lat), cos = Math.cos(lat);
-    const radius = a / Math.sqrt(1 - es * sin * sin);
-    point.x = (radius + height) * cos * Math.cos(lon);
-    point.y = (radius + height) * cos * Math.sin(lon);
-    point.z = (radius * (1 - es) + height) * sin;
+    if (!spheroidToCartesian(point, ellipsoid))
+      throw new Error("Geocentric coordinate must be finite");
   }
   function geocentricToGeodeticInPlace(point, ellipsoid) {
     const x = point.x, y = point.y, z = point.z;
     const { semiMajorAxis: a, semiMinorAxis: b, eccentricitySquared: es } = ellipsoid;
-    const p = Math.hypot(x, y), rr = Math.hypot(x, y, z);
-    if (rr === 0)
+    if (x === 0 && y === 0 && z === 0)
       throw new Error("Geodetic coordinates are undefined at the Earth center");
-    if (p < 1e-12 * a) {
+    if (es !== 0 && Math.abs(x) < 1e-12 * a && Math.abs(y) < 1e-12 * a && Math.hypot(x, y) < 1e-12 * a) {
+      if (!Number.isFinite(Math.hypot(x, y, z)))
+        throw new Error("Geocentric radius must be finite");
       point.x = 0;
       point.y = Math.sign(z) * Math.PI / 2;
       point.z = Math.abs(z) - b;
       return;
     }
-    const ct = z / rr, st = p / rr;
-    let rx = 1 / Math.sqrt(1 - es * (2 - es) * st * st);
-    let cos = st * (1 - es) * rx, sin = ct * rx;
-    for (let i = 0; i < 30; i++) {
-      const rn = a / Math.sqrt(1 - es * sin * sin);
-      const height = p * cos + z * sin - rn * (1 - es * sin * sin);
-      const rk = es * rn / (rn + height);
-      rx = 1 / Math.sqrt(1 - rk * (2 - rk) * st * st);
-      const nextCos = st * (1 - rk) * rx, nextSin = ct * rx;
-      if (Math.abs(nextSin * cos - nextCos * sin) <= 1e-12) {
-        point.x = Math.atan2(y, x);
-        point.y = Math.atan2(nextSin, Math.abs(nextCos));
-        point.z = height;
-        return;
-      }
-      cos = nextCos;
-      sin = nextSin;
+    if (!cartesianToSpheroid(point, ellipsoid)) {
+      if (!Number.isFinite(Math.hypot(x, y, z)))
+        throw new Error("Geocentric radius must be finite");
+      throw new Error("Geocentric inverse did not converge");
     }
-    throw new Error("Geocentric inverse did not converge");
   }
   function createHelmert(values, inverse33) {
     const dx = values[0], dy = values[1], dz = values[2];
@@ -13295,7 +13501,7 @@ var __exports__ = (() => {
     return [point.x, point.y, point.z];
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/vertical-datum.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/vertical-datum.js
   function compileVerticalGrid(names, grids2 = {}) {
     if (names === void 0)
       return void 0;
@@ -13332,7 +13538,7 @@ var __exports__ = (() => {
     };
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/typescript-projection.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/typescript-projection.js
   var ProjectionEngine = class _ProjectionEngine {
     /** Resolve required descriptors once, then return an ordinary synchronous instance. */
     static async create(options = {}) {
@@ -13340,6 +13546,7 @@ var __exports__ = (() => {
       return new _ProjectionEngine({ ...options, projections });
     }
     constructor(options = {}) {
+      this.coordinateScratch = new ProjectionScratch();
       const registrations = options.projections || [];
       const plugins = registry(registrations);
       if (registrations.some((projection) => !("create" in projection))) {
@@ -13362,6 +13569,10 @@ var __exports__ = (() => {
       }
       this.project = this.project.bind(this);
       this.unproject = this.unproject.bind(this);
+      this.projectTo = this.projectTo.bind(this);
+      this.unprojectTo = this.unprojectTo.bind(this);
+      this.projectToSync = this.projectToSync.bind(this);
+      this.unprojectToSync = this.unprojectToSync.bind(this);
       this.projectFlat = this.projectFlat.bind(this);
       this.unprojectFlat = this.unprojectFlat.bind(this);
       this.projectSync = this.projectSync.bind(this);
@@ -13401,16 +13612,16 @@ var __exports__ = (() => {
     }
     /** Synchronous variants never import algorithms; descriptors must be preloaded. */
     projectSync(coordinate) {
-      return this.deferred ? this.loadSync().project(coordinate) : transformScalar(coordinate, this.forwardTransform);
+      return this.deferred ? this.loadSync().project(coordinate) : transformScalar(coordinate, this.forwardTransform, this.coordinateScratch);
     }
     unprojectSync(coordinate) {
-      return this.deferred ? this.loadSync().unproject(coordinate) : transformScalar(coordinate, this.inverseTransform);
+      return this.deferred ? this.loadSync().unproject(coordinate) : transformScalar(coordinate, this.inverseTransform, this.coordinateScratch);
     }
     projectFlatSync(coordinates, dimension = 2) {
-      return this.deferred ? this.loadSync().projectFlat(coordinates, dimension) : transformInPlace(coordinates, dimension, this.forwardTransform);
+      return this.deferred ? this.loadSync().projectFlat(coordinates, dimension) : transformInPlace(coordinates, dimension, this.forwardTransform, this.coordinateScratch);
     }
     unprojectFlatSync(coordinates, dimension = 2) {
-      return this.deferred ? this.loadSync().unprojectFlat(coordinates, dimension) : transformInPlace(coordinates, dimension, this.inverseTransform);
+      return this.deferred ? this.loadSync().unprojectFlat(coordinates, dimension) : transformInPlace(coordinates, dimension, this.inverseTransform, this.coordinateScratch);
     }
     /** Optional warm-up. Coordinate methods also load automatically on first use. */
     async preload() {
@@ -13418,24 +13629,48 @@ var __exports__ = (() => {
     }
     /** Eager plugins return arrays; descriptor-backed instances return promises. */
     project(coordinate) {
-      const result = this.deferred ? this.projectLoaded(coordinate, false) : transformScalar(coordinate, this.forwardTransform);
+      const result = this.deferred ? this.projectLoaded(coordinate, false) : transformScalar(coordinate, this.forwardTransform, this.coordinateScratch);
       return result;
     }
     unproject(coordinate) {
-      const result = this.deferred ? this.projectLoaded(coordinate, true) : transformScalar(coordinate, this.inverseTransform);
+      const result = this.deferred ? this.projectLoaded(coordinate, true) : transformScalar(coordinate, this.inverseTransform, this.coordinateScratch);
       return result;
     }
+    /** Write one coordinate into caller-owned storage; return that same output. */
+    projectTo(coordinate, output) {
+      return this.deferred ? this.projectToLoaded(coordinate, output, false) : this.projectToSync(coordinate, output);
+    }
+    unprojectTo(coordinate, output) {
+      return this.deferred ? this.projectToLoaded(coordinate, output, true) : this.unprojectToSync(coordinate, output);
+    }
+    projectToSync(coordinate, output) {
+      if (this.deferred)
+        return this.loadSync().projectTo(coordinate, output);
+      validateScalarOutput(coordinate, output, Math.max(coordinate.length, this.forwardTransform.geocentricOutput ? 3 : 2));
+      return transformScalar(coordinate, this.forwardTransform, this.coordinateScratch, output);
+    }
+    unprojectToSync(coordinate, output) {
+      if (this.deferred)
+        return this.loadSync().unprojectTo(coordinate, output);
+      validateScalarOutput(coordinate, output, Math.max(coordinate.length, this.inverseTransform.geocentricOutput ? 3 : 2));
+      return transformScalar(coordinate, this.inverseTransform, this.coordinateScratch, output);
+    }
+    projectToLoaded(coordinate, output, inverse33) {
+      validateScalarOutput(coordinate, output, coordinate.length);
+      const input = coordinate.slice();
+      return this.load().then((projection) => inverse33 ? projection.unprojectTo(input, output) : projection.projectTo(input, output));
+    }
     projectLoaded(coordinate, inverse33) {
       const input = coordinate.slice();
       return this.load().then((projection) => inverse33 ? projection.unproject(input) : projection.project(input));
     }
     /** In-place; descriptor callers must await completion before reusing the buffer. */
     projectFlat(coordinates, dimension = 2) {
-      const result = this.deferred ? this.load().then((projection) => projection.projectFlat(coordinates, dimension)) : transformInPlace(coordinates, dimension, this.forwardTransform);
+      const result = this.deferred ? this.load().then((projection) => projection.projectFlat(coordinates, dimension)) : transformInPlace(coordinates, dimension, this.forwardTransform, this.coordinateScratch);
       return result;
     }
     unprojectFlat(coordinates, dimension = 2) {
-      const result = this.deferred ? this.load().then((projection) => projection.unprojectFlat(coordinates, dimension)) : transformInPlace(coordinates, dimension, this.inverseTransform);
+      const result = this.deferred ? this.load().then((projection) => projection.unprojectFlat(coordinates, dimension)) : transformInPlace(coordinates, dimension, this.inverseTransform, this.coordinateScratch);
       return result;
     }
   };
@@ -13593,21 +13828,31 @@ var __exports__ = (() => {
       }
     };
   }
-  function transformScalar(coordinate, operation2) {
+  function transformScalar(coordinate, operation2, scratch, output) {
     if (coordinate.length < 2)
       throw new Error("Coordinates must contain finite x, y and optional z values");
     if (coordinate.length < 3 && operation2.requiresInputZ)
       throw new Error("This transform requires three ordinates");
-    const point = { x: coordinate[0], y: coordinate[1], z: coordinate.length >= 3 ? coordinate[2] : 0 };
-    operation2.run(point);
-    const output = coordinate.slice();
-    output[0] = point.x;
-    output[1] = point.y;
-    if (coordinate.length >= 3 || operation2.geocentricOutput)
-      output[2] = point.z;
-    return output;
+    const length4 = Math.max(coordinate.length, operation2.geocentricOutput ? 3 : 2);
+    const point = scratch.acquire();
+    try {
+      point.x = coordinate[0];
+      point.y = coordinate[1];
+      point.z = coordinate.length >= 3 ? coordinate[2] : 0;
+      operation2.run(point);
+      if (output)
+        return writeScalarOutput(coordinate, output, point, length4);
+      const result = coordinate.slice();
+      result[0] = point.x;
+      result[1] = point.y;
+      if (coordinate.length >= 3 || operation2.geocentricOutput)
+        result[2] = point.z;
+      return result;
+    } finally {
+      scratch.release(point);
+    }
   }
-  function transformInPlace(coordinates, dimension, operation2) {
+  function transformInPlace(coordinates, dimension, operation2, scratch) {
     if (!(coordinates instanceof Float32Array || coordinates instanceof Float64Array))
       throw new Error("In-place projection requires a Float32Array or Float64Array");
     if (!Number.isSafeInteger(dimension) || dimension < 2 || coordinates.length % dimension !== 0)
@@ -13618,24 +13863,28 @@ var __exports__ = (() => {
       operation2.flat(coordinates, dimension);
       return coordinates;
     }
-    const point = { x: 0, y: 0, z: 0 };
-    const float32 = coordinates instanceof Float32Array;
-    for (let offset = 0; offset < coordinates.length; offset += dimension) {
-      point.x = coordinates[offset];
-      point.y = coordinates[offset + 1];
-      point.z = dimension >= 3 ? coordinates[offset + 2] : 0;
-      operation2.run(point);
-      if (float32 && (Math.abs(point.x) > 34028234663852886e22 || Math.abs(point.y) > 34028234663852886e22 || dimension >= 3 && Math.abs(point.z) > 34028234663852886e22))
-        throw new Error("Projected coordinate exceeds Float32 range");
-      coordinates[offset] = point.x;
-      coordinates[offset + 1] = point.y;
-      if (dimension >= 3)
-        coordinates[offset + 2] = point.z;
+    const point = scratch.acquire();
+    try {
+      const float32 = coordinates instanceof Float32Array;
+      for (let offset = 0; offset < coordinates.length; offset += dimension) {
+        point.x = coordinates[offset];
+        point.y = coordinates[offset + 1];
+        point.z = dimension >= 3 ? coordinates[offset + 2] : 0;
+        operation2.run(point);
+        if (float32 && (Math.abs(point.x) > 34028234663852886e22 || Math.abs(point.y) > 34028234663852886e22 || dimension >= 3 && Math.abs(point.z) > 34028234663852886e22))
+          throw new Error("Projected coordinate exceeds Float32 range");
+        coordinates[offset] = point.x;
+        coordinates[offset + 1] = point.y;
+        if (dimension >= 3)
+          coordinates[offset + 2] = point.z;
+      }
+      return coordinates;
+    } finally {
+      scratch.release(point);
     }
-    return coordinates;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/flat-projection.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/flat-projection.js
   function createFlatProjection(forward33, inverse33) {
     return {
       ...createProjection(forward33, inverse33),
@@ -13683,7 +13932,7 @@ var __exports__ = (() => {
     };
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/mercator.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/mercator.js
   var mercator = {
     name: "merc",
     aliases: [
@@ -13739,7 +13988,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/equidistant-cylindrical.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/equidistant-cylindrical.js
   var equidistantCylindrical = {
     name: "eqc",
     aliases: ["Equirectangular", "Equidistant_Cylindrical", "Equidistant_Cylindrical_Spherical"],
@@ -13761,7 +14010,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernel.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernel.js
   var ORIGIN_PARAMETERS = ["lon_0", "lat_0", "x_0", "y_0"];
   var SCALE_PARAMETERS = ["k", "k_0"];
   function angleParameter(parameters, name, fallback = 0) {
@@ -13829,13 +14078,13 @@ var __exports__ = (() => {
     }, projectionOperation(implementation, true));
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/msfnz.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/msfnz.js
   function msfnz_default(eccent, sinphi, cosphi) {
     var con = eccent * sinphi;
     return cosphi / Math.sqrt(1 - con * con);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/constants.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/constants.js
   var EPSLN = 1e-10;
   var HALF_PI = Math.PI / 2;
   var FORTPI = Math.PI / 4;
@@ -13845,7 +14094,7 @@ var __exports__ = (() => {
   var R2D = 180 / Math.PI;
   var SEC_TO_RAD = Math.PI / (180 * 3600);
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/tsfnz.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/tsfnz.js
   function tsfnz_default(eccent, phi, sinphi) {
     var con = eccent * sinphi;
     var com = 0.5 * eccent;
@@ -13853,12 +14102,12 @@ var __exports__ = (() => {
     return Math.tan(0.5 * (HALF_PI - phi)) / con;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/sign.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/sign.js
   function sign_default(x) {
     return x < 0 ? -1 : 1;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/adjust_lon.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/adjust_lon.js
   function adjust_lon_default(x, skipAdjust) {
     if (skipAdjust) {
       return x;
@@ -13866,7 +14115,7 @@ var __exports__ = (() => {
     return Math.abs(x) <= SPI ? x : x - sign_default(x) * TWO_PI2;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/phi2z.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/phi2z.js
   function phi2z_default(eccent, ts) {
     var eccnth = 0.5 * eccent;
     var con, dphi;
@@ -13882,7 +14131,7 @@ var __exports__ = (() => {
     return -9999;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/lcc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/lcc.js
   function initialize(state) {
     if (!state.k0) {
       state.k0 = 1;
@@ -13980,7 +14229,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/lcc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/lcc.js
   var lambertConformalConic = {
     name: "lcc",
     aliases: [
@@ -14007,7 +14256,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/qsfnz.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/qsfnz.js
   function qsfnz_default(eccent, sinphi) {
     var con;
     if (eccent > 1e-7) {
@@ -14018,7 +14267,7 @@ var __exports__ = (() => {
     }
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/asinz.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/asinz.js
   function asinz_default(x) {
     if (Math.abs(x) > 1) {
       x = x > 1 ? 1 : -1;
@@ -14026,7 +14275,7 @@ var __exports__ = (() => {
     return Math.asin(x);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/aea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/aea.js
   function initialize2(state) {
     if (Math.abs(state.lat1 + state.lat2) < EPSLN) {
       return;
@@ -14144,7 +14393,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/aea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/aea.js
   var albersEqualArea = {
     name: "aea",
     aliases: ["Albers_Conic_Equal_Area", "Albers_Equal_Area", "Albers"],
@@ -14156,7 +14405,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/pj_enfn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/pj_enfn.js
   var C00 = 1;
   var C02 = 0.25;
   var C04 = 0.046875;
@@ -14181,14 +14430,14 @@ var __exports__ = (() => {
     return en;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/pj_mlfn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/pj_mlfn.js
   function pj_mlfn_default(phi, sphi, cphi, en) {
     cphi *= sphi;
     sphi *= sphi;
     return en[0] * phi - cphi * (en[1] + sphi * (en[2] + sphi * (en[3] + sphi * en[4])));
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/pj_inv_mlfn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/pj_inv_mlfn.js
   var MAX_ITER = 20;
   function pj_inv_mlfn_default(arg, es, en) {
     var k = 1 / (1 - es);
@@ -14205,12 +14454,12 @@ var __exports__ = (() => {
     return phi;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/adjust_lat.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/adjust_lat.js
   function adjust_lat_default(x) {
     return Math.abs(x) < HALF_PI ? x : x - sign_default(x) * Math.PI;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/eqdc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/eqdc.js
   function initialize3(state) {
     if (Math.abs(state.lat1 + state.lat2) < EPSLN) {
       return;
@@ -14303,7 +14552,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/eqdc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/eqdc.js
   var equidistantConic = {
     name: "eqdc",
     aliases: ["Equidistant_Conic"],
@@ -14315,7 +14564,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/authset.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/authset.js
   var P00 = 0.3333333333333333;
   var P01 = 0.17222222222222222;
   var P02 = 0.10257936507936508;
@@ -14336,13 +14585,13 @@ var __exports__ = (() => {
     return APA;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/authlat.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/authlat.js
   function authlat(beta, APA) {
     var t = beta + beta;
     return beta + APA[0] * Math.sin(t) + APA[1] * Math.sin(t + t) + APA[2] * Math.sin(t + t + t);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/laea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/laea.js
   var S_POLE = 1;
   var N_POLE = 2;
   var EQUIT = 3;
@@ -14583,7 +14832,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/laea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/laea.js
   var lambertAzimuthalEqualArea = {
     name: "laea",
     aliases: ["Lambert Azimuthal Equal Area", "Lambert_Azimuthal_Equal_Area"],
@@ -14594,7 +14843,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/stere.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/stere.js
   function ssfn_(phit, sinphi, eccen) {
     sinphi *= eccen;
     return Math.tan(0.5 * (HALF_PI + phit)) * Math.pow((1 - sinphi) / (1 + sinphi), 0.5 * eccen);
@@ -14739,7 +14988,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/stere.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/stere.js
   var stereographic = {
     name: "stere",
     aliases: [
@@ -14755,12 +15004,12 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/srat.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/srat.js
   function srat_default(esinp, exp) {
     return Math.pow((1 - esinp) / (1 + esinp), exp);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/gauss.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/gauss.js
   var MAX_ITER2 = 20;
   function initialize6(state) {
     var sphi = Math.sin(state.lat0);
@@ -14811,7 +15060,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/hypot.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/hypot.js
   function hypot_default(x, y) {
     x = Math.abs(x);
     y = Math.abs(y);
@@ -14820,7 +15069,7 @@ var __exports__ = (() => {
     return a * Math.sqrt(1 + Math.pow(b, 2));
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/sterea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/sterea.js
   function initialize7(state) {
     if (!state.rc) {
       return;
@@ -14878,7 +15127,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/sterea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/sterea.js
   var obliqueStereographic = {
     name: "sterea",
     aliases: [
@@ -14896,32 +15145,32 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/mlfn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/mlfn.js
   function mlfn_default(e0, e1, e2, e3, phi) {
     return e0 * phi - e1 * Math.sin(2 * phi) + e2 * Math.sin(4 * phi) - e3 * Math.sin(6 * phi);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/e0fn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/e0fn.js
   function e0fn_default(x) {
     return 1 - 0.25 * x * (1 + x / 16 * (3 + 1.25 * x));
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/e1fn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/e1fn.js
   function e1fn_default(x) {
     return 0.375 * x * (1 + 0.25 * x * (1 + 0.46875 * x));
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/e2fn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/e2fn.js
   function e2fn_default(x) {
     return 0.05859375 * x * x * (1 + 0.75 * x);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/e3fn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/e3fn.js
   function e3fn_default(x) {
     return x * x * x * (35 / 3072);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/imlfn.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/imlfn.js
   function imlfn_default(ml, e0, e1, e2, e3) {
     var phi;
     var dphi;
@@ -14936,8 +15185,8 @@ var __exports__ = (() => {
     return NaN;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/vincenty.js
-  function vincentyInverse(lat1, lon1, lat2, lon2, a, f) {
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/vincenty.js
+  function vincentyInverse(lat1, lon1, lat2, lon2, a, f, output) {
     const L = lon2 - lon1;
     const U1 = Math.atan((1 - f) * Math.tan(lat1));
     const U2 = Math.atan((1 - f) * Math.tan(lat2));
@@ -14951,7 +15200,9 @@ var __exports__ = (() => {
       cosLambda = Math.cos(lambda);
       sinSigma = Math.sqrt(cosU2 * sinLambda * (cosU2 * sinLambda) + (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda) * (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda));
       if (sinSigma === 0) {
-        return { azi1: 0, s12: 0 };
+        output.x = 0;
+        output.y = 0;
+        return;
       }
       cosSigma = sinU1 * sinU2 + cosU1 * cosU2 * cosLambda;
       sigma = Math.atan2(sinSigma, cosSigma);
@@ -14963,7 +15214,9 @@ var __exports__ = (() => {
       lambda = L + (1 - C) * f * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));
     } while (Math.abs(lambda - lambdaP) > 1e-12 && --iterLimit > 0);
     if (iterLimit === 0) {
-      return { azi1: NaN, s12: NaN };
+      output.x = NaN;
+      output.y = NaN;
+      return;
     }
     uSq = cos2Alpha * (a * a - a * (1 - f) * (a * (1 - f))) / (a * (1 - f) * (a * (1 - f)));
     A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
@@ -14971,9 +15224,10 @@ var __exports__ = (() => {
     deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) - B / 6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM)));
     s = a * (1 - f) * A * (sigma - deltaSigma);
     const azi1 = Math.atan2(cosU2 * sinLambda, cosU1 * sinU2 - sinU1 * cosU2 * cosLambda);
-    return { azi1, s12: s };
+    output.x = azi1;
+    output.y = s;
   }
-  function vincentyDirect(lat1, lon1, azi1, s12, a, f) {
+  function vincentyDirect(lat1, lon1, azi1, s12, a, f, output) {
     const U1 = Math.atan((1 - f) * Math.tan(lat1));
     const sinU1 = Math.sin(U1), cosU1 = Math.cos(U1);
     const sinAlpha1 = Math.sin(azi1), cosAlpha1 = Math.cos(azi1);
@@ -14994,7 +15248,9 @@ var __exports__ = (() => {
       sigma = s12 / (a * (1 - f) * A) + deltaSigma;
     } while (Math.abs(sigma - sigmaP) > 1e-12 && --iterLimit > 0);
     if (iterLimit === 0) {
-      return { lat2: NaN, lon2: NaN };
+      output.x = NaN;
+      output.y = NaN;
+      return;
     }
     const tmp = sinU1 * sinSigma - cosU1 * cosSigma * cosAlpha1;
     const lat2 = Math.atan2(sinU1 * cosSigma + cosU1 * sinSigma * cosAlpha1, (1 - f) * Math.sqrt(sinAlpha * sinAlpha + tmp * tmp));
@@ -15002,10 +15258,11 @@ var __exports__ = (() => {
     const C = f / 16 * cos2Alpha * (4 + f * (4 - 3 * cos2Alpha));
     const L = lambda - (1 - C) * f * sinAlpha * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));
     const lon2 = lon1 + L;
-    return { lat2, lon2 };
+    output.x = lon2;
+    output.y = lat2;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/aeqd.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/aeqd.js
   function initialize8(state) {
     state.sin_p12 = Math.sin(state.lat0);
     state.cos_p12 = Math.cos(state.lat0);
@@ -15020,7 +15277,7 @@ var __exports__ = (() => {
     var sinphi = Math.sin(p.y);
     var cosphi = Math.cos(p.y);
     var dlon = adjust_lon_default(lon - state.long0, state.over);
-    var e0, e1, e2, e3, Mlp, Ml, c, kp, cos_c, vars, azi1;
+    var e0, e1, e2, e3, Mlp, Ml, c, kp, cos_c, azi1;
     if (state.sphere) {
       if (Math.abs(state.sin_p12 - 1) <= EPSLN) {
         p.x = state.x0 + state.a * (HALF_PI - lat) * Math.sin(dlon);
@@ -15061,10 +15318,11 @@ var __exports__ = (() => {
           p.y = state.y0;
           return p;
         }
-        vars = vincentyInverse(state.lat0, state.long0, lat, lon, state.a, state.f);
-        azi1 = vars.azi1;
-        p.x = state.x0 + vars.s12 * Math.sin(azi1);
-        p.y = state.y0 + vars.s12 * Math.cos(azi1);
+        vincentyInverse(state.lat0, state.long0, lat, lon, state.a, state.f, p);
+        azi1 = p.x;
+        const distance = p.y;
+        p.x = state.x0 + distance * Math.sin(azi1);
+        p.y = state.y0 + distance * Math.cos(azi1);
         return p;
       }
     }
@@ -15072,7 +15330,7 @@ var __exports__ = (() => {
   function inverse8(state, p) {
     p.x -= state.x0;
     p.y -= state.y0;
-    var rh, z, sinz, cosz, lon, lat, con, e0, e1, e2, e3, Mlp, M2, azi1, s12, vars;
+    var rh, z, sinz, cosz, lon, lat, con, e0, e1, e2, e3, Mlp, M2, azi1, s12;
     if (state.sphere) {
       rh = Math.sqrt(p.x * p.x + p.y * p.y);
       if (rh > 2 * HALF_PI * state.a) {
@@ -15126,9 +15384,7 @@ var __exports__ = (() => {
       } else {
         azi1 = Math.atan2(p.x, p.y);
         s12 = Math.sqrt(p.x * p.x + p.y * p.y);
-        vars = vincentyDirect(state.lat0, state.long0, azi1, s12, state.a, state.f);
-        p.x = vars.lon2;
-        p.y = vars.lat2;
+        vincentyDirect(state.lat0, state.long0, azi1, s12, state.a, state.f, p);
         return p;
       }
     }
@@ -15144,7 +15400,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/aeqd.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/aeqd.js
   var azimuthalEquidistant = {
     name: "aeqd",
     aliases: ["Azimuthal_Equidistant"],
@@ -15155,28 +15411,28 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/sinh.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/sinh.js
   function sinh_default(x) {
     var r = Math.exp(x);
     r = (r - 1 / r) / 2;
     return r;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/log1py.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/log1py.js
   function log1py_default(x) {
     var y = 1 + x;
     var z = y - 1;
     return z === 0 ? x : x * Math.log(y) / z;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/asinhy.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/asinhy.js
   function asinhy_default(x) {
     var y = Math.abs(x);
     y = log1py_default(y * (1 + y / (hypot_default(1, y) + 1)));
     return x < 0 ? -y : y;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/gatg.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/gatg.js
   function gatg_default(pp, B) {
     var cos_2B = 2 * Math.cos(2 * B);
     var i = pp.length - 1;
@@ -15191,7 +15447,7 @@ var __exports__ = (() => {
     return B + h * Math.sin(2 * B);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/clens.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/clens.js
   function clens_default(pp, arg_r) {
     var r = 2 * Math.cos(arg_r);
     var i = pp.length - 1;
@@ -15206,14 +15462,14 @@ var __exports__ = (() => {
     return Math.sin(arg_r) * hr;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/cosh.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/cosh.js
   function cosh_default(x) {
     var r = Math.exp(x);
     r = (r + 1 / r) / 2;
     return r;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/clens_cmplx.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/clens_cmplx.js
   function clens_cmplx_default(pp, arg_r, arg_i, output) {
     var sin_arg_r = Math.sin(arg_r);
     var cos_arg_r = Math.cos(arg_r);
@@ -15242,7 +15498,7 @@ var __exports__ = (() => {
     output.y = r * hi + i * hr;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/etmerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/etmerc.js
   function initialize9(state) {
     state.x0 = state.x0 !== void 0 ? state.x0 : 0;
     state.y0 = state.y0 !== void 0 ? state.y0 : 0;
@@ -15363,7 +15619,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/tmerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/tmerc.js
   function initialize10(state) {
     state.x0 = state.x0 !== void 0 ? state.x0 : 0;
     state.y0 = state.y0 !== void 0 ? state.y0 : 0;
@@ -15478,7 +15734,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/transverse-mercator.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/transverse-mercator.js
   function createTransverseMercator(context) {
     const base = kernelParameters(context);
     if (flagParameter(context.parameters, "approx")) {
@@ -15510,7 +15766,7 @@ var __exports__ = (() => {
     create: createTransverseMercator
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/utm.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/utm.js
   var universalTransverseMercator = {
     name: "utm",
     aliases: ["Universal Transverse Mercator System"],
@@ -15536,16 +15792,16 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/geocentric.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/geocentric.js
   var geocentric = {
     name: "geocent",
     aliases: ["Geocentric", "geocentric", "Geocent"],
     parameters: [],
-    create({ semiMajorAxis, eccentricitySquared }) {
+    create({ semiMajorAxis, semiMinorAxis, eccentricitySquared }) {
       const ellipsoid = {
         semiMajorAxis,
         eccentricitySquared,
-        semiMinorAxis: semiMajorAxis * Math.sqrt(1 - eccentricitySquared)
+        semiMinorAxis: semiMinorAxis ?? semiMajorAxis * Math.sqrt(1 - eccentricitySquared)
       };
       return {
         forward() {
@@ -15562,7 +15818,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/catalogue-kernel.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/catalogue-kernel.js
   function catalogueParameters(context) {
     const base = kernelParameters(context);
     return { ...base, x0: 0, y0: 0, lat_ts: Number.isNaN(base.lat_ts) ? 0 : base.lat_ts };
@@ -15589,7 +15845,7 @@ var __exports__ = (() => {
     };
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/bonne.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/bonne.js
   var EPS10 = 1e-10;
   function initialize11(state) {
     var c;
@@ -15696,7 +15952,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/bonne.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/bonne.js
   var bonne = {
     name: "bonne",
     aliases: ["Bonne (Werner lat_1=90)"],
@@ -15710,13 +15966,13 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/gN.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/gN.js
   function gN_default(a, e, sinphi) {
     var temp = e * sinphi;
     return a / Math.sqrt(1 - temp * temp);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/cass.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/cass.js
   function initialize12(state) {
     if (!state.sphere) {
       state.en = pj_enfn_default(state.es);
@@ -15816,7 +16072,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/cass.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/cass.js
   var cassiniSoldner = {
     name: "cass",
     aliases: ["Cassini", "Cassini_Soldner"],
@@ -15828,7 +16084,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/iqsfnz.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/iqsfnz.js
   function iqsfnz_default(eccent, q) {
     var temp = 1 - (1 - eccent * eccent) / (2 * eccent) * Math.log((1 - eccent) / (1 + eccent));
     const roundoff = 4 * Number.EPSILON * temp;
@@ -15859,7 +16115,7 @@ var __exports__ = (() => {
     return NaN;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/cea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/cea.js
   function initialize13(state) {
     if (!state.sphere) {
       state.k0 = msfnz_default(state.e, Math.sin(state.lat_ts), Math.cos(state.lat_ts));
@@ -15906,7 +16162,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/cea.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/cea.js
   var cylindricalEqualArea = {
     name: "cea",
     aliases: [],
@@ -15920,7 +16176,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/sinu.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/sinu.js
   var MAX_ITER3 = 20;
   function initialize14(state) {
     state.long0 = state.long0 || 0;
@@ -16001,7 +16257,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/eck6.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/eck6.js
   function initialize15(state) {
     state.sphere = true;
     state.b = state.a;
@@ -16019,7 +16275,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/eck6.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/eck6.js
   var eckertVI = {
     name: "eck6",
     aliases: ["Eckert_VI"],
@@ -16031,7 +16287,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/eqearth.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/eqearth.js
   var A1 = 1.340264;
   var A2 = -0.081106;
   var A3 = 893e-6;
@@ -16099,7 +16355,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/eqearth.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/eqearth.js
   var equalEarth = {
     name: "eqearth",
     aliases: ["Equal Earth", "Equal_Earth"],
@@ -16111,7 +16367,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/equi.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/equi.js
   function initialize17(state) {
     state.x0 = state.x0 || 0;
     state.y0 = state.y0 || 0;
@@ -16145,7 +16401,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/equi.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/equi.js
   var equirectangular = {
     name: "equi",
     aliases: [],
@@ -16159,7 +16415,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/geos.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/geos.js
   function initialize18(state) {
     state.flip_axis = state.sweep === "x" ? 1 : 0;
     state.h = Number(state.h);
@@ -16301,7 +16557,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/geos.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/geos.js
   var geostationary = {
     name: "geos",
     aliases: ["Geostationary Satellite View", "Geostationary_Satellite"],
@@ -16319,7 +16575,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/gnom.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/gnom.js
   function initialize19(state) {
     state.sin_p14 = Math.sin(state.lat0);
     state.cos_p14 = Math.cos(state.lat0);
@@ -16391,7 +16647,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/gnom.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/gnom.js
   var gnomonic = {
     name: "gnom",
     aliases: [],
@@ -16403,7 +16659,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/latiso.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/latiso.js
   function latiso_default(eccent, phi, sinphi) {
     if (Math.abs(phi) > HALF_PI) {
       return Number.NaN;
@@ -16418,12 +16674,12 @@ var __exports__ = (() => {
     return Math.log(Math.tan((HALF_PI + phi) / 2)) + eccent * Math.log((1 - con) / (1 + con)) / 2;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/fL.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/fL.js
   function fL_default(x, L) {
     return 2 * Math.atan(x * Math.exp(L)) - HALF_PI;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/common/invlatiso.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/common/invlatiso.js
   function invlatiso_default(eccent, ts) {
     var phi = fL_default(1, ts);
     var Iphi = 0;
@@ -16439,7 +16695,7 @@ var __exports__ = (() => {
     return phi;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/gstmerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/gstmerc.js
   function initialize20(state) {
     if (!state.k0) {
       state.k0 = 1;
@@ -16483,7 +16739,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/gstmerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/gstmerc.js
   var gaussSchreiberTransverseMercator = {
     name: "gstmerc",
     aliases: ["gstmerg"],
@@ -16495,7 +16751,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/krovak.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/krovak.js
   function initialize21(state) {
     state.a = 6377397155e-3;
     state.es = 0.006674372230614;
@@ -16595,7 +16851,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/krovak.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/krovak.js
   var krovak = {
     name: "krovak",
     aliases: [
@@ -16622,7 +16878,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/mill.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/mill.js
   function initialize22(state) {
   }
   function forward22(state, p) {
@@ -16653,7 +16909,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/mill.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/mill.js
   var millerCylindrical = {
     name: "mill",
     aliases: ["Miller_Cylindrical"],
@@ -16665,7 +16921,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/moll.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/moll.js
   function initialize23(state) {
     state.x0 = state.x0 !== void 0 ? state.x0 : 0;
     state.y0 = state.y0 !== void 0 ? state.y0 : 0;
@@ -16735,7 +16991,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/moll.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/moll.js
   var mollweide = {
     name: "moll",
     aliases: ["Mollweide"],
@@ -16747,7 +17003,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/nzmg.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/nzmg.js
   function initialize24(state) {
     state.A = [];
     state.A[1] = 0.6399175073;
@@ -16916,7 +17172,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/nzmg.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/nzmg.js
   var newZealandMapGrid = {
     name: "nzmg",
     aliases: ["New_Zealand_Map_Grid"],
@@ -16933,7 +17189,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/omerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/omerc.js
   var TOL = 1e-7;
   function isTypeA(P) {
     return P.no_uoff || P.no_off || [
@@ -17058,7 +17314,7 @@ var __exports__ = (() => {
     state.v_pole_s = state.ArB * Math.log(Math.tan(FORTPI + F));
   }
   function forward25(state, p) {
-    var coords = { x: 0, y: 0 };
+    var coords = p;
     var S, T, U, V, W, temp, u, v;
     p.x = p.x - state.lam0;
     if (Math.abs(Math.abs(p.y) - HALF_PI) > EPSLN) {
@@ -17096,7 +17352,7 @@ var __exports__ = (() => {
   }
   function inverse25(state, p) {
     var u, v, Qp, Sp, Tp, Vp, Up;
-    var coords = { x: 0, y: 0 };
+    var coords = p;
     p.x = (p.x - state.x0) * (1 / state.a);
     p.y = (p.y - state.y0) * (1 / state.a);
     if (state.no_rot) {
@@ -17157,7 +17413,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/omerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/omerc.js
   var obliqueMercator = {
     name: "omerc",
     aliases: [
@@ -17213,7 +17469,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/ortho.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/ortho.js
   function initialize26(state) {
     state.sin_p14 = Math.sin(state.lat0 || 0);
     state.cos_p14 = Math.cos(state.lat0 || 0);
@@ -17287,7 +17543,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/ortho.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/ortho.js
   var orthographic = {
     name: "ortho",
     aliases: [],
@@ -17299,7 +17555,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/poly.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/poly.js
   var MAX_ITER4 = 20;
   function initialize27(state) {
     state.temp = state.b / state.a;
@@ -17402,7 +17658,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/poly.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/poly.js
   var polyconic = {
     name: "poly",
     aliases: ["Polyconic", "American_Polyconic"],
@@ -17414,7 +17670,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/qsc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/qsc.js
   var FACE_ENUM = {
     FRONT: 1,
     RIGHT: 2,
@@ -17452,11 +17708,11 @@ var __exports__ = (() => {
     }
   }
   function forward28(state, p) {
-    var xy = { x: 0, y: 0 };
+    var xy = p;
     var lat, lon;
     var theta, phi;
     var t, mu;
-    var area = { value: 0 };
+    var area = 0;
     p.x -= state.long0;
     if (state.es !== 0) {
       lat = Math.atan(state.one_minus_f_squared * Math.tan(p.y));
@@ -17467,31 +17723,31 @@ var __exports__ = (() => {
     if (state.face === FACE_ENUM.TOP) {
       phi = HALF_PI - lat;
       if (lon >= FORTPI && lon <= HALF_PI + FORTPI) {
-        area.value = AREA_ENUM.AREA_0;
+        area = AREA_ENUM.AREA_0;
         theta = lon - HALF_PI;
       } else if (lon > HALF_PI + FORTPI || lon <= -(HALF_PI + FORTPI)) {
-        area.value = AREA_ENUM.AREA_1;
+        area = AREA_ENUM.AREA_1;
         theta = lon > 0 ? lon - SPI : lon + SPI;
       } else if (lon > -(HALF_PI + FORTPI) && lon <= -FORTPI) {
-        area.value = AREA_ENUM.AREA_2;
+        area = AREA_ENUM.AREA_2;
         theta = lon + HALF_PI;
       } else {
-        area.value = AREA_ENUM.AREA_3;
+        area = AREA_ENUM.AREA_3;
         theta = lon;
       }
     } else if (state.face === FACE_ENUM.BOTTOM) {
       phi = HALF_PI + lat;
       if (lon >= FORTPI && lon <= HALF_PI + FORTPI) {
-        area.value = AREA_ENUM.AREA_0;
+        area = AREA_ENUM.AREA_0;
         theta = -lon + HALF_PI;
       } else if (lon < FORTPI && lon >= -FORTPI) {
-        area.value = AREA_ENUM.AREA_1;
+        area = AREA_ENUM.AREA_1;
         theta = -lon;
       } else if (lon < -FORTPI && lon >= -(HALF_PI + FORTPI)) {
-        area.value = AREA_ENUM.AREA_2;
+        area = AREA_ENUM.AREA_2;
         theta = -lon - HALF_PI;
       } else {
-        area.value = AREA_ENUM.AREA_3;
+        area = AREA_ENUM.AREA_3;
         theta = lon > 0 ? -lon + SPI : -lon - SPI;
       }
     } else {
@@ -17514,58 +17770,67 @@ var __exports__ = (() => {
       s = sinlat;
       if (state.face === FACE_ENUM.FRONT) {
         phi = Math.acos(q);
-        theta = qsc_fwd_equat_face_theta(phi, s, r, area);
+        theta = qsc_fwd_equat_face_theta(phi, s, r);
       } else if (state.face === FACE_ENUM.RIGHT) {
         phi = Math.acos(r);
-        theta = qsc_fwd_equat_face_theta(phi, s, -q, area);
+        theta = qsc_fwd_equat_face_theta(phi, s, -q);
       } else if (state.face === FACE_ENUM.BACK) {
         phi = Math.acos(-q);
-        theta = qsc_fwd_equat_face_theta(phi, s, -r, area);
+        theta = qsc_fwd_equat_face_theta(phi, s, -r);
       } else if (state.face === FACE_ENUM.LEFT) {
         phi = Math.acos(-r);
-        theta = qsc_fwd_equat_face_theta(phi, s, q, area);
+        theta = qsc_fwd_equat_face_theta(phi, s, q);
       } else {
         phi = theta = 0;
-        area.value = AREA_ENUM.AREA_0;
+      }
+      if (Math.abs(theta) <= FORTPI) {
+        area = AREA_ENUM.AREA_0;
+      } else if (theta > FORTPI && theta <= HALF_PI + FORTPI) {
+        area = AREA_ENUM.AREA_1;
+        theta -= HALF_PI;
+      } else if (theta > HALF_PI + FORTPI || theta <= -(HALF_PI + FORTPI)) {
+        area = AREA_ENUM.AREA_2;
+        theta = theta >= 0 ? theta - SPI : theta + SPI;
+      } else {
+        area = AREA_ENUM.AREA_3;
+        theta += HALF_PI;
       }
     }
     mu = Math.atan(12 / SPI * (theta + Math.acos(Math.sin(theta) * Math.cos(FORTPI)) - HALF_PI));
     t = Math.sqrt((1 - Math.cos(phi)) / (Math.cos(mu) * Math.cos(mu)) / (1 - Math.cos(Math.atan(1 / Math.cos(theta)))));
-    if (area.value === AREA_ENUM.AREA_1) {
+    if (area === AREA_ENUM.AREA_1) {
       mu += HALF_PI;
-    } else if (area.value === AREA_ENUM.AREA_2) {
+    } else if (area === AREA_ENUM.AREA_2) {
       mu += SPI;
-    } else if (area.value === AREA_ENUM.AREA_3) {
+    } else if (area === AREA_ENUM.AREA_3) {
       mu += 1.5 * SPI;
     }
     xy.x = t * Math.cos(mu);
     xy.y = t * Math.sin(mu);
     xy.x = xy.x * state.a + state.x0;
     xy.y = xy.y * state.a + state.y0;
-    p.x = xy.x;
-    p.y = xy.y;
     return p;
   }
   function inverse28(state, p) {
-    var lp = { lam: 0, phi: 0 };
+    var lam = 0, latitude = 0;
     var mu, nu, cosmu, tannu;
     var tantheta, theta, cosphi, phi;
     var t;
-    var area = { value: 0 };
+    var area = 0;
     p.x = (p.x - state.x0) / state.a;
     p.y = (p.y - state.y0) / state.a;
     nu = Math.atan(Math.sqrt(p.x * p.x + p.y * p.y));
     mu = Math.atan2(p.y, p.x);
     if (p.x >= 0 && p.x >= Math.abs(p.y)) {
-      area.value = AREA_ENUM.AREA_0;
+      area = AREA_ENUM.AREA_0;
     } else if (p.y >= 0 && p.y >= Math.abs(p.x)) {
-      area.value = AREA_ENUM.AREA_1;
+      area = AREA_ENUM.AREA_1;
       mu -= HALF_PI;
     } else if (p.x < 0 && -p.x >= Math.abs(p.y)) {
-      area.value = AREA_ENUM.AREA_2;
+      area = AREA_ENUM.AREA_2;
       mu = mu < 0 ? mu + SPI : mu - SPI;
     } else {
-      area.value = AREA_ENUM.AREA_3;
+      area = AREA_ENUM.AREA_3;
       mu += HALF_PI;
     }
     t = SPI / 12 * Math.tan(mu);
@@ -17581,27 +17846,27 @@ var __exports__ = (() => {
     }
     if (state.face === FACE_ENUM.TOP) {
       phi = Math.acos(cosphi);
-      lp.phi = HALF_PI - phi;
-      if (area.value === AREA_ENUM.AREA_0) {
-        lp.lam = theta + HALF_PI;
-      } else if (area.value === AREA_ENUM.AREA_1) {
-        lp.lam = theta < 0 ? theta + SPI : theta - SPI;
-      } else if (area.value === AREA_ENUM.AREA_2) {
-        lp.lam = theta - HALF_PI;
+      latitude = HALF_PI - phi;
+      if (area === AREA_ENUM.AREA_0) {
+        lam = theta + HALF_PI;
+      } else if (area === AREA_ENUM.AREA_1) {
+        lam = theta < 0 ? theta + SPI : theta - SPI;
+      } else if (area === AREA_ENUM.AREA_2) {
+        lam = theta - HALF_PI;
       } else {
-        lp.lam = theta;
+        lam = theta;
       }
     } else if (state.face === FACE_ENUM.BOTTOM) {
       phi = Math.acos(cosphi);
-      lp.phi = phi - HALF_PI;
-      if (area.value === AREA_ENUM.AREA_0) {
-        lp.lam = -theta + HALF_PI;
-      } else if (area.value === AREA_ENUM.AREA_1) {
-        lp.lam = -theta;
-      } else if (area.value === AREA_ENUM.AREA_2) {
-        lp.lam = -theta - HALF_PI;
+      latitude = phi - HALF_PI;
+      if (area === AREA_ENUM.AREA_0) {
+        lam = -theta + HALF_PI;
+      } else if (area === AREA_ENUM.AREA_1) {
+        lam = -theta;
+      } else if (area === AREA_ENUM.AREA_2) {
+        lam = -theta - HALF_PI;
       } else {
-        lp.lam = theta < 0 ? -theta - SPI : -theta + SPI;
+        lam = theta < 0 ? -theta - SPI : -theta + SPI;
       }
     } else {
       var q, r, s;
@@ -17618,14 +17883,14 @@ var __exports__ = (() => {
       } else {
         r = Math.sqrt(1 - t);
       }
-      if (area.value === AREA_ENUM.AREA_1) {
+      if (area === AREA_ENUM.AREA_1) {
         t = r;
         r = -s;
         s = t;
-      } else if (area.value === AREA_ENUM.AREA_2) {
+      } else if (area === AREA_ENUM.AREA_2) {
         r = -r;
         s = -s;
-      } else if (area.value === AREA_ENUM.AREA_3) {
+      } else if (area === AREA_ENUM.AREA_3) {
         t = r;
         r = s;
         s = -t;
@@ -17642,53 +17907,34 @@ var __exports__ = (() => {
         q = r;
         r = -t;
       }
-      lp.phi = Math.acos(-s) - HALF_PI;
-      lp.lam = Math.atan2(r, q);
+      latitude = Math.acos(-s) - HALF_PI;
+      lam = Math.atan2(r, q);
       if (state.face === FACE_ENUM.RIGHT) {
-        lp.lam = qsc_shift_lon_origin(lp.lam, -HALF_PI);
+        lam = qsc_shift_lon_origin(lam, -HALF_PI);
       } else if (state.face === FACE_ENUM.BACK) {
-        lp.lam = qsc_shift_lon_origin(lp.lam, -SPI);
+        lam = qsc_shift_lon_origin(lam, -SPI);
       } else if (state.face === FACE_ENUM.LEFT) {
-        lp.lam = qsc_shift_lon_origin(lp.lam, +HALF_PI);
+        lam = qsc_shift_lon_origin(lam, +HALF_PI);
       }
     }
     if (state.es !== 0) {
       var invert_sign;
       var tanphi, xa;
-      invert_sign = lp.phi < 0 ? 1 : 0;
-      tanphi = Math.tan(lp.phi);
+      invert_sign = latitude < 0 ? 1 : 0;
+      tanphi = Math.tan(latitude);
       xa = state.b / Math.sqrt(tanphi * tanphi + state.one_minus_f_squared);
-      lp.phi = Math.atan(Math.sqrt(state.a * state.a - xa * xa) / (state.one_minus_f * xa));
+      latitude = Math.atan(Math.sqrt(state.a * state.a - xa * xa) / (state.one_minus_f * xa));
       if (invert_sign) {
-        lp.phi = -lp.phi;
+        latitude = -latitude;
       }
     }
-    lp.lam += state.long0;
-    p.x = lp.lam;
-    p.y = lp.phi;
+    lam += state.long0;
+    p.x = lam;
+    p.y = latitude;
     return p;
   }
-  function qsc_fwd_equat_face_theta(phi, y, x, area) {
-    var theta;
-    if (phi < EPSLN) {
-      area.value = AREA_ENUM.AREA_0;
-      theta = 0;
-    } else {
-      theta = Math.atan2(y, x);
-      if (Math.abs(theta) <= FORTPI) {
-        area.value = AREA_ENUM.AREA_0;
-      } else if (theta > FORTPI && theta <= HALF_PI + FORTPI) {
-        area.value = AREA_ENUM.AREA_1;
-        theta -= HALF_PI;
-      } else if (theta > HALF_PI + FORTPI || theta <= -(HALF_PI + FORTPI)) {
-        area.value = AREA_ENUM.AREA_2;
-        theta = theta >= 0 ? theta - SPI : theta + SPI;
-      } else {
-        area.value = AREA_ENUM.AREA_3;
-        theta += HALF_PI;
-      }
-    }
-    return theta;
+  function qsc_fwd_equat_face_theta(phi, y, x) {
+    return phi < EPSLN ? 0 : Math.atan2(y, x);
   }
   function qsc_shift_lon_origin(lon, offset) {
     var slon = lon + offset;
@@ -17705,7 +17951,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/qsc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/qsc.js
   var quadrilateralizedSphericalCube = {
     name: "qsc",
     aliases: ["Quadrilateralized Spherical Cube", "Quadrilateralized_Spherical_Cube"],
@@ -17717,7 +17963,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/robin.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/robin.js
   var COEFS_X = [
     [1, 22199e-21, -715515e-10, 31103e-10],
     [0.9986, -482243e-9, -24897e-9, -13309e-10],
@@ -17773,10 +18019,10 @@ var __exports__ = (() => {
   var poly3_der = function(coefs, x) {
     return coefs[1] + x * (2 * coefs[2] + x * 3 * coefs[3]);
   };
-  function newton_rapshon(f_df, start, max_err, iters) {
+  function newton_rapshon(coefs, target, start, max_err, iters) {
     var x = start;
     for (; iters; --iters) {
-      var upd = f_df(x);
+      var upd = (poly3_val(coefs, x) - target) / poly3_der(coefs, x);
       x -= upd;
       if (Math.abs(upd) < max_err) {
         break;
@@ -17800,32 +18046,30 @@ var __exports__ = (() => {
       i = NODES;
     }
     dphi = R2D * (dphi - RC1 * i);
-    var xy = {
-      x: poly3_val(COEFS_X[i], dphi) * lon,
-      y: poly3_val(COEFS_Y[i], dphi)
-    };
+    var x = poly3_val(COEFS_X[i], dphi) * lon;
+    var y = poly3_val(COEFS_Y[i], dphi);
     if (ll.y < 0) {
-      xy.y = -xy.y;
+      y = -y;
     }
-    xy.x = xy.x * state.a * FXC + state.x0;
-    xy.y = xy.y * state.a * FYC + state.y0;
-    return xy;
+    ll.x = x * state.a * FXC + state.x0;
+    ll.y = y * state.a * FYC + state.y0;
+    return ll;
   }
   function inverse29(state, xy) {
-    var ll = {
-      x: (xy.x - state.x0) / (state.a * FXC),
-      y: Math.abs(xy.y - state.y0) / (state.a * FYC)
-    };
+    var originalY = xy.y;
+    var ll = xy;
+    ll.x = (ll.x - state.x0) / (state.a * FXC);
+    ll.y = Math.abs(originalY - state.y0) / (state.a * FYC);
     for (let knot = 0; knot <= NODES; knot++) {
       if (Math.abs(ll.y - COEFS_Y[knot][0]) <= 8 * Number.EPSILON) {
         ll.x = adjust_lon_default(ll.x / COEFS_X[knot][0] + state.long0, state.over);
-        ll.y = (xy.y < state.y0 ? -1 : 1) * knot * 5 * D2R;
+        ll.y = (originalY < state.y0 ? -1 : 1) * knot * 5 * D2R;
         return ll;
       }
     }
     if (ll.y >= 1) {
       ll.x /= COEFS_X[NODES][0];
-      ll.y = xy.y < 0 ? -HALF_PI : HALF_PI;
+      ll.y = originalY < 0 ? -HALF_PI : HALF_PI;
     } else {
       var i = Math.floor(ll.y * NODES);
       if (i < 0) {
@@ -17844,12 +18088,10 @@ var __exports__ = (() => {
       }
       var coefs = COEFS_Y[i];
       var t = 5 * (ll.y - coefs[0]) / (COEFS_Y[i + 1][0] - coefs[0]);
-      t = newton_rapshon(function(x) {
-        return (poly3_val(coefs, x) - ll.y) / poly3_der(coefs, x);
-      }, t, EPSLN, 100);
+      t = newton_rapshon(coefs, ll.y, t, EPSLN, 100);
       ll.x /= poly3_val(COEFS_X[i], t);
       ll.y = (5 * i + t) * D2R;
-      if (xy.y < 0) {
+      if (originalY < 0) {
         ll.y = -ll.y;
       }
     }
@@ -17865,7 +18107,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/robin.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/robin.js
   var robinson = {
     name: "robin",
     aliases: ["Robinson"],
@@ -17877,7 +18119,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/sinu.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/sinu.js
   var sinusoidal = {
     name: "sinu",
     aliases: ["Sinusoidal"],
@@ -17889,7 +18131,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/somerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/somerc.js
   function initialize30(state) {
     if (!state.k0) {
       state.k0 = 1;
@@ -17953,7 +18195,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/somerc.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/somerc.js
   var swissObliqueMercator = {
     name: "somerc",
     aliases: [],
@@ -17965,7 +18207,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/tpers.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/tpers.js
   var mode = {
     N_POLE: 0,
     S_POLE: 1,
@@ -18046,7 +18288,7 @@ var __exports__ = (() => {
   function inverse31(state, p) {
     p.x /= state.a;
     p.y /= state.a;
-    var r = { x: p.x, y: p.y };
+    var rx = p.x, ry = p.y;
     var bm, bq, yt;
     yt = 1 / (state.pn1 - p.y * state.sw);
     bm = state.pn1 * p.x * yt;
@@ -18055,8 +18297,8 @@ var __exports__ = (() => {
     p.y = bq * state.cg - bm * state.sg;
     var rh = hypot_default(p.x, p.y);
     if (Math.abs(rh) < EPSLN) {
-      r.x = 0;
-      r.y = state.lat0;
+      rx = 0;
+      ry = state.lat0;
     } else {
       var cosz, sinz;
       sinz = 1 - rh * rh * state.pfact;
@@ -18064,27 +18306,27 @@ var __exports__ = (() => {
       cosz = Math.sqrt(1 - sinz * sinz);
       switch (state.mode) {
         case mode.OBLIQ:
-          r.y = Math.asin(cosz * state.sinph0 + p.y * sinz * state.cosph0 / rh);
-          p.y = (cosz - state.sinph0 * Math.sin(r.y)) * rh;
+          ry = Math.asin(cosz * state.sinph0 + p.y * sinz * state.cosph0 / rh);
+          p.y = (cosz - state.sinph0 * Math.sin(ry)) * rh;
           p.x *= sinz * state.cosph0;
           break;
         case mode.EQUIT:
-          r.y = Math.asin(p.y * sinz / rh);
+          ry = Math.asin(p.y * sinz / rh);
           p.y = cosz * rh;
           p.x *= sinz;
           break;
         case mode.N_POLE:
-          r.y = Math.asin(cosz);
+          ry = Math.asin(cosz);
           p.y = -p.y;
           break;
         case mode.S_POLE:
-          r.y = -Math.asin(cosz);
+          ry = -Math.asin(cosz);
           break;
       }
-      r.x = Math.atan2(p.x, p.y);
+      rx = Math.atan2(p.x, p.y);
     }
-    p.x = r.x + state.long0;
-    p.y = r.y;
+    p.x = rx + state.long0;
+    p.y = ry;
     return p;
   }
   function createState31(base, options = {}) {
@@ -18111,7 +18353,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/tpers.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/tpers.js
   var tiltedPerspective = {
     name: "tpers",
     aliases: ["Tilted_Perspective"],
@@ -18128,7 +18370,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/kernels/vandg.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/kernels/vandg.js
   function initialize32(state) {
     state.R = state.a;
   }
@@ -18234,7 +18476,7 @@ var __exports__ = (() => {
     return state;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/vandg.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/vandg.js
   var vanDerGrinten = {
     name: "vandg",
     aliases: ["Van_der_Grinten_I", "VanDerGrinten", "Van_der_Grinten"],
@@ -18246,7 +18488,7 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/projections/ob-tran.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/projections/ob-tran.js
   var rotationParameters = [
     "o_proj",
     "o_alpha",
@@ -18351,7 +18593,7 @@ var __exports__ = (() => {
     };
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/structured.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/structured.js
   function object(value2) {
     if (!value2 || typeof value2 !== "object" || Array.isArray(value2))
       throw new Error("Expected a CRS object");
@@ -18761,7 +19003,7 @@ var __exports__ = (() => {
     return result;
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/wkt.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/wkt.js
   function children(node, ...names) {
     return node.values.filter((value2) => value2.type === "node" && names.includes(value2.keyword.toUpperCase()));
   }
@@ -18922,14 +19164,14 @@ var __exports__ = (() => {
     }
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/crs/projjson.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/crs/projjson.js
   var projJSONCRSParser = {
     name: "projjson",
     canParse: (definition) => inferCRSRepresentation(definition) === "projjson",
     parse: (definition, options) => readStructuredCRS(object(definition), options)
   };
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/grids/grid.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/grids/grid.js
   var TOLERANCE = 1e-12;
   function interpolate(grid, x, y, output, clamp2 = false) {
     let col = (x - grid.origin[0]) / grid.step[0], row = (y - grid.origin[1]) / grid.step[1];
@@ -18946,11 +19188,12 @@ var __exports__ = (() => {
       const weight = (right ? u : 1 - u) * (upper ? v : 1 - v);
       if (weight === 0)
         continue;
-      const node = grid.shifts[(j + (upper ? 1 : 0)) * width + i + right];
-      if (!Number.isFinite(node[0]) || !Number.isFinite(node[1]))
+      const index = ((j + (upper ? 1 : 0)) * width + i + right) * 2;
+      const xShift = grid.shifts[index], yShift = grid.shifts[index + 1];
+      if (!Number.isFinite(xShift) || !Number.isFinite(yShift))
         return false;
-      dx += weight * node[0];
-      dy += weight * node[1];
+      dx += weight * xShift;
+      dy += weight * yShift;
     }
     output.x = dx;
     output.y = dy;
@@ -18993,13 +19236,14 @@ var __exports__ = (() => {
       throw new Error("Datum grid requires at least one subgrid");
     const prepared = subgrids.map((grid) => {
       const [width, height] = grid.size;
-      if (![width, height].every((n) => Number.isSafeInteger(n) && n >= 2) || width * height !== grid.shifts.length || !grid.origin.every(Number.isFinite) || !grid.step.every((n) => Number.isFinite(n) && n > 0))
+      if (![width, height].every((n) => Number.isSafeInteger(n) && n >= 2) || width * height * 2 !== grid.shifts.length || !grid.origin.every(Number.isFinite) || !grid.step.every((n) => Number.isFinite(n) && n > 0))
         throw new Error("Invalid datum subgrid geometry or node count");
       const maxShift = [0, 0];
-      for (const node of grid.shifts)
-        for (let i = 0; i < 2; i++)
-          if (Number.isFinite(node[i]))
-            maxShift[i] = Math.max(maxShift[i], Math.abs(node[i]));
+      for (let index = 0; index < grid.shifts.length; index++) {
+        const value2 = grid.shifts[index];
+        if (Number.isFinite(value2))
+          maxShift[index % 2] = Math.max(maxShift[index % 2], Math.abs(value2));
+      }
       return { ...grid, maxShift };
     });
     const shiftInPlace = (point, inverse33) => {
@@ -19016,17 +19260,19 @@ var __exports__ = (() => {
       point.y = latitude;
       return false;
     };
+    const tupleScratch = { x: NaN, y: NaN };
     return Object.freeze({
       subgridCount: prepared.length,
       shiftInPlace,
       shift(longitude, latitude, inverse33) {
-        const point = { x: longitude, y: latitude };
-        return shiftInPlace(point, inverse33) ? [point.x, point.y] : void 0;
+        tupleScratch.x = longitude;
+        tupleScratch.y = latitude;
+        return shiftInPlace(tupleScratch, inverse33) ? [tupleScratch.x, tupleScratch.y] : void 0;
       }
     });
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/experimental/grids/ntv2.js
+  // ../../node_modules/@math.gl/projection/dist/experimental/grids/ntv2.js
   var ARC_SECOND2 = Math.PI / 648e3;
   function parseNTv2Grid(data, options = {}) {
     const view = new DataView(data);
@@ -19056,12 +19302,11 @@ var __exports__ = (() => {
       const start = offset + 176, end = start + nodes * stride;
       if (end > view.byteLength)
         throw new Error("Truncated NTv2 node records");
-      const shifts = [];
-      for (let i = 0; i < nodes; i++)
-        shifts.push([
-          view.getFloat32(start + i * stride + 4, little) * ARC_SECOND2,
-          view.getFloat32(start + i * stride, little) * ARC_SECOND2
-        ]);
+      const shifts = new Float64Array(nodes * 2);
+      for (let i = 0; i < nodes; i++) {
+        shifts[i * 2] = view.getFloat32(start + i * stride + 4, little) * ARC_SECOND2;
+        shifts[i * 2 + 1] = view.getFloat32(start + i * stride, little) * ARC_SECOND2;
+      }
       subgrids.push({
         origin: [lowerLon * ARC_SECOND2, lowerLat * ARC_SECOND2],
         step: [lonStep * ARC_SECOND2, latStep * ARC_SECOND2],
@@ -19073,7 +19318,7 @@ var __exports__ = (() => {
     return createDatumGrid(subgrids);
   }
 
-  // ../../node_modules/@math.gl/proj4/dist/lib/typescript-proj4-projection.js
+  // ../../node_modules/@math.gl/projection/dist/lib/projection.js
   var aliases = /* @__PURE__ */ Object.create(null);
   var grids = /* @__PURE__ */ Object.create(null);
   var preset;
@@ -19143,14 +19388,14 @@ var __exports__ = (() => {
     return preset;
   }
   var Projection = class extends ProjectionEngine {
-    /** Register aliases for subsequently constructed TypeScript wrappers. */
+    /** Register aliases for subsequently constructed projections. */
     static defineProjectionAliases(definitions) {
       for (const name of Object.keys(definitions)) {
         const definition = definitions[name];
         aliases[name] = typeof definition === "string" ? definition : JSON.parse(JSON.stringify(definition));
       }
     }
-    /** Register an owned NTv2 grid snapshot for subsequent TypeScript wrappers. */
+    /** Register an owned NTv2 grid snapshot for subsequent projections. */
     static registerDatumGrid(name, grid, options) {
       grids[name] = parseNTv2Grid(grid, options);
     }
@@ -19892,15 +20137,15 @@ var __exports__ = (() => {
       this.outputIsGeographic = isWgs84Geographic(outputCrs);
       if (spatialReference.sourceCrs && spatialReference.targetCrs && !this.sourceIsGeocentric && !this.outputIsGeocentric) {
         this.horizontalProjection = new Projection({
-          from: getHorizontalProj4Definition(spatialReference.sourceCrs),
-          to: getHorizontalProj4Definition(spatialReference.targetCrs),
+          from: getHorizontalProjectionDefinition(spatialReference.sourceCrs),
+          to: getHorizontalProjectionDefinition(spatialReference.targetCrs),
           enforceAxis: false
         });
       }
       if (requiresHeightTransformation(spatialReference) || this.sourceIsGeocentric || this.outputIsGeocentric) {
         if (!this.sourceIsGeographic && !this.sourceIsGeocentric) {
           this.geographicProjection = new Projection({
-            from: getHorizontalProj4Definition(spatialReference.sourceCrs),
+            from: getHorizontalProjectionDefinition(spatialReference.sourceCrs),
             to: WGS84_GEOGRAPHIC_CRS,
             enforceAxis: false
           });
@@ -19908,7 +20153,7 @@ var __exports__ = (() => {
         if (!this.outputIsGeographic && !this.outputIsGeocentric) {
           this.heightOutputProjection = new Projection({
             from: WGS84_GEOGRAPHIC_CRS,
-            to: getHorizontalProj4Definition(outputCrs),
+            to: getHorizontalProjectionDefinition(outputCrs),
             enforceAxis: false
           });
         }
@@ -20062,7 +20307,7 @@ var __exports__ = (() => {
       (row) => (row[0] * vector[0] + row[1] * vector[1] + row[2] * vector[2]) / determinant3
     );
   }
-  function getHorizontalProj4Definition(definition) {
+  function getHorizontalProjectionDefinition(definition) {
     if (!definition) {
       throw new Error("Cannot construct a projection because the CRS is unknown");
     }
@@ -20073,7 +20318,7 @@ var __exports__ = (() => {
       if (horizontalComponents.length !== 1) {
         throw new Error("CompoundCRS requires exactly one horizontal component");
       }
-      return getHorizontalProj4Definition(horizontalComponents[0]);
+      return getHorizontalProjectionDefinition(horizontalComponents[0]);
     }
     return definition;
   }
@@ -21393,6 +21638,8 @@ var __exports__ = (() => {
     tileUrlCache = /* @__PURE__ */ new Map();
     /** Original content descriptors retained across {@link Tile3D.unloadContent} calls. */
     tileContentHeaders = /* @__PURE__ */ new WeakMap();
+    /** Installed nested roots by parent placement and ordered nested-content slot, across reloads. */
+    nestedTilesetRoots = /* @__PURE__ */ new WeakMap();
     /** Source-local namespaces assigned to distinct embedded glTF package file collections. */
     packageNamespaces = /* @__PURE__ */ new WeakMap();
     /** Next source-local package namespace. */
@@ -21823,7 +22070,9 @@ var __exports__ = (() => {
       return this.extensionsUsed.includes(extensionName);
     }
     /**
-     * Updates content-format flags and installs nested tileset subtrees.
+     * Updates content-format flags and installs nested tileset subtrees once per content slot.
+     * Reloading payloads retains attached roots and their descendants, including duplicate URLs
+     * declared in separate slots or parent placements. Detached roots are reinstalled in declaration order.
      */
     onTileLoaded(tileset, tile, loadResult) {
       const contents = tile.contents.length ? tile.contents : [tile.content];
@@ -21838,8 +22087,23 @@ var __exports__ = (() => {
         this.contentFormats.ktx2 = true;
       }
       const nestedTilesets = loadResult.nestedTilesets || (loadResult.nestedTileset ? [loadResult.nestedTileset] : []);
-      for (const nestedTileset of nestedTilesets) {
-        tileset._initializeTileHeaders(this.prepareNestedTileset(tileset, nestedTileset), tile);
+      const installedRoots = this.nestedTilesetRoots.get(tile) || [];
+      for (const [nestedIndex, nestedTileset] of nestedTilesets.entries()) {
+        const installedRoot = installedRoots[nestedIndex];
+        if (installedRoot && tile.children.includes(installedRoot)) {
+          continue;
+        }
+        const nextRoot = installedRoots.slice(nestedIndex + 1).find((root) => tile.children.includes(root));
+        const nestedRoot = tileset._initializeTileHeaders(
+          this.prepareNestedTileset(tileset, nestedTileset),
+          tile
+        );
+        if (nextRoot) {
+          tile.children.splice(tile.children.indexOf(nestedRoot), 1);
+          tile.children.splice(tile.children.indexOf(nextRoot), 0, nestedRoot);
+        }
+        installedRoots[nestedIndex] = nestedRoot;
+        this.nestedTilesetRoots.set(tile, installedRoots);
       }
     }
     /**