@math.gl/geospatial 5.0.0-alpha.19 → 5.0.0-alpha.20

package/src/projection/experimental/common/phi2z.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import {HALF_PI} from './constants';
 

package/src/projection/experimental/common/pj_enfn.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 var C00 = 1;
 var C02 = 0.25;

package/src/projection/experimental/common/pj_inv_mlfn.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import pj_mlfn from './pj_mlfn';
 import {EPSLN} from './constants';

package/src/projection/experimental/common/pj_mlfn.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (phi: number, sphi: number, cphi: number, en: readonly number[]) {
   cphi *= sphi;

package/src/projection/experimental/common/qsfnz.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (eccent: number, sinphi: number) {
   var con;

package/src/projection/experimental/common/sign.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (x: number) {
   return x < 0 ? -1 : 1;

package/src/projection/experimental/common/sinh.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (x: number) {
   var r = Math.exp(x);

package/src/projection/experimental/common/srat.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (esinp: number, exp: number) {
   return Math.pow((1 - esinp) / (1 + esinp), exp);

package/src/projection/experimental/common/tsfnz.ts

@@ -1,13 +1,14 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
-
-import {HALF_PI} from './constants';
+// SPDX-FileCopyrightText: Copyright (c) PROJ contributors
+// SPDX-FileComment: Adapted from proj4js 2.22.0; stable half-angle evaluation adapted from PROJ 9.5.1 src/tsfn.cpp (MIT). See ../../../../projection/PROJ-LICENSE.txt. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 export default function (eccent: number, phi: number, sinphi: number) {
-  var con = eccent * sinphi;
-  var com = 0.5 * eccent;
-  con = Math.pow((1 - con) / (1 + con), com);
-  return Math.tan(0.5 * (HALF_PI - phi)) / con;
+  // PROJ's sign-dependent half-angle identities avoid cancellation near poles.
+  const cosphi = Math.cos(phi);
+  return (
+    Math.exp(eccent * Math.atanh(eccent * sinphi)) *
+    (sinphi > 0 ? cosphi / (1 + sinphi) : (1 - sinphi) / cosphi)
+  );
 }

package/src/projection/experimental/common/vincenty.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0, modified to write into caller-owned storage without temporary result objects. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0, modified to write into caller-owned storage without temporary result objects. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {Point} from '../kernel';
 

package/src/projection/experimental/crs/datum-table.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md.
 
 const datums: Record<
   string,

package/src/projection/experimental/crs/ellipsoid-table.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md.
 
 const ellipsoids: Record<string, {a: number; b?: number; rf?: number; ellipseName: string}> = {
   MERIT: {

package/src/projection/experimental/crs/primemeridian-table.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md.
 
 const primeMeridian: Record<string, number> = {};
 

package/src/projection/experimental/crs/structured.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Original adapter; datum lookup and method/parameter normalization rules are adapted from proj4js 2.22.0 and its MIT-licensed wkt-parser dependency. See ../../../PROJ4-LICENSE.md.
+// SPDX-FileComment: Original adapter; datum lookup and method/parameter normalization rules are adapted from proj4js 2.22.0 and its MIT-licensed wkt-parser dependency. See ../../../../projection/PROJ4-LICENSE.md.
 import {datumNameKey, getDatumDefinitions} from './datum-catalog';
 import {DEGREES_TO_RADIANS} from '../parameters';
 import {unsupportedStage} from './types';

package/src/projection/experimental/crs/wkt.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Original adapter; named-axis fallback inspired by proj4js 2.22.0 / wkt-parser (MIT). See ../../../PROJ4-LICENSE.md for upstream attribution.
+// SPDX-FileComment: Original adapter; named-axis fallback inspired by proj4js 2.22.0 / wkt-parser (MIT). See ../../../../projection/PROJ4-LICENSE.md for upstream attribution.
 import {inferCRSRepresentation} from '../../../crs/spatial-reference';
 import {parseWKTCRS} from '../../../crs/wkt-crs';
 import type {WKTCRSNode} from '../../../crs/index';

package/src/projection/experimental/datum.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Helmert equations directly adapted from proj4js 2.22.0 datumUtils.js. See ../../PROJ4-LICENSE.md. Geocentric equations now delegate to the attributed @math.gl/core/spheroid leaf.
+// SPDX-FileComment: Helmert equations directly adapted from proj4js 2.22.0 datumUtils.js. See ../../../projection/PROJ4-LICENSE.md. Geocentric equations now delegate to the attributed @math.gl/core/spheroid leaf.
 import {spheroidToCartesian, cartesianToSpheroid} from '@math.gl/core/spheroid';
 import type {Datum, Ellipsoid} from './crs/types';
 import type {ProjectionPoint} from './types';

package/src/projection/experimental/ellipsoidal-gnomonic.ts

@@ -0,0 +1,83 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) PROJ contributors
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Ellipsoidal gnomonic dispatch and inverse Newton equations adapted from PROJ 9.5.1 src/projections/gnom.cpp (Karney geodesic formulation). See ../../../projection/PROJ-LICENSE.txt. Geodesic evaluation uses the shared ellipsoid engine; public optional polygon/line facade code stays out of projection bundles.
+import {createGeodesicEngine, nonAreaGeodesicMask} from '../../ellipsoid-geodesic';
+import type {KernelParameters, Point} from './kernel';
+import adjust_lon from './common/adjust_lon';
+
+const degrees = 180 / Math.PI;
+export type EllipsoidalGnomonic = {
+  engine: ReturnType<typeof createGeodesicEngine>;
+  latitude: number;
+};
+
+export function prepareEllipsoidalGnomonic(base: KernelParameters): EllipsoidalGnomonic {
+  return {
+    engine: createGeodesicEngine({semiMajorAxis: 1, semiMinorAxis: Math.sqrt(1 - base.es)}),
+    latitude: base.lat0 * degrees
+  };
+}
+
+export function forwardEllipsoidalGnomonic(
+  base: KernelParameters,
+  state: EllipsoidalGnomonic,
+  p: Point
+): Point {
+  const result = state.engine.Inverse(
+    state.latitude,
+    0,
+    p.y * degrees,
+    adjust_lon(p.x - base.long0, base.over) * degrees,
+    nonAreaGeodesicMask
+  );
+  if (!(result.M12 > 0)) throw new Error('Coordinate is outside the ellipsoidal gnomonic horizon');
+  const radius = (base.a * result.m12) / result.M12;
+  const azimuth = result.azi1 / degrees;
+  p.x = base.x0 + radius * Math.sin(azimuth);
+  p.y = base.y0 + radius * Math.cos(azimuth);
+  return p;
+}
+
+export function inverseEllipsoidalGnomonic(
+  base: KernelParameters,
+  state: EllipsoidalGnomonic,
+  p: Point
+): Point {
+  const x = (p.x - base.x0) / base.a;
+  const y = (p.y - base.y0) / base.a;
+  const radius = Math.hypot(x, y);
+  if (radius === 0) {
+    p.x = base.long0;
+    p.y = base.lat0;
+    return p;
+  }
+  const line = state.engine.Line(
+    state.latitude,
+    0,
+    Math.atan2(x, y) * degrees,
+    nonAreaGeodesicMask
+  );
+  const little = radius <= 1;
+  const target = little ? radius : 1 / radius;
+  let distance = Math.atan(radius);
+  // d(m/M)/ds = 1/M². Use its reciprocal near the horizon to avoid
+  // overflow, and evaluate once more at the converged distance.
+  for (let iteration = 0; iteration < 20; iteration++) {
+    const result = line.Position(distance, nonAreaGeodesicMask);
+    const step = little
+      ? (result.m12 - target * result.M12) * result.M12
+      : (target * result.m12 - result.M12) * result.m12;
+    if (!Number.isFinite(step)) break;
+    distance -= step;
+    if (Math.abs(step) < 1e-14) {
+      const final = line.Position(distance, nonAreaGeodesicMask);
+      if (!(final.M12 > 0)) break;
+      p.x = adjust_lon(base.long0 + final.lon2 / degrees, base.over);
+      p.y = final.lat2 / degrees;
+      return p;
+    }
+  }
+  throw new Error('Ellipsoidal gnomonic inverse did not converge');
+}

package/src/projection/experimental/exact-helmert.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: Copyright (c) 2016, Thomas Knudsen / SDFE (PROJ)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Exact rotation matrix adapted from PROJ 9.5.1 src/transformations/helmert.cpp. See ../../PROJ-LICENSE.txt. Modified into prepared, typed, allocation-free point operations. Existing small-angle Helmert equations remain the proj4js adaptation in datum.ts.
+// SPDX-FileComment: Exact rotation matrix adapted from PROJ 9.5.1 src/transformations/helmert.cpp. See ../../../projection/PROJ-LICENSE.txt. Modified into prepared, typed, allocation-free point operations. Existing small-angle Helmert equations remain the proj4js adaptation in datum.ts.
 import type {ProjectionPoint} from './types';
 
 type Operation = (point: ProjectionPoint) => void;

package/src/projection/experimental/grids/geotiff.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: GeoTIFF node orientation directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../PROJ4-LICENSE.md. Modified for explicit async loading, validation and nodata handling; no geotiff runtime import.
+// SPDX-FileComment: GeoTIFF node orientation directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../../projection/PROJ4-LICENSE.md. Modified for explicit async loading, validation and nodata handling; no geotiff runtime import.
 import {createDatumGrid} from './grid';
 import type {DatumGrid, Subgrid} from './types';
 

package/src/projection/experimental/grids/grid.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Bilinear interpolation and inverse grid equations directly adapted from proj4js 2.22.0 lib/datum_transform.js. See ../../../PROJ4-LICENSE.md. Modified for immutable explicit state, inclusive edges, nodata and bounded two-axis convergence.
+// SPDX-FileComment: Bilinear interpolation and inverse grid equations directly adapted from proj4js 2.22.0 lib/datum_transform.js. See ../../../../projection/PROJ4-LICENSE.md. Modified for immutable explicit state, inclusive edges, nodata and bounded two-axis convergence.
 import type {DatumGrid, Subgrid} from './types';
 import {wrapLongitude} from '../parameters';
 

package/src/projection/experimental/grids/ntv2.ts

@@ -2,7 +2,7 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: NTv2 decoding directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../PROJ4-LICENSE.md. Modified for explicit ownership, bounds/geometry validation and no global registration.
+// SPDX-FileComment: NTv2 decoding directly adapted from proj4js 2.22.0 lib/nadgrid.js. See ../../../../projection/PROJ4-LICENSE.md. Modified for explicit ownership, bounds/geometry validation and no global registration.
 import {createDatumGrid} from './grid';
 import type {DatumGrid, Subgrid} from './types';
 

package/src/projection/experimental/helmert-flat.ts

@@ -3,7 +3,7 @@
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) 2016, Thomas Knudsen / SDFE (PROJ)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Original whole-buffer dispatch. Small-angle equations follow the proj4js adaptation in datum.ts; see ../../PROJ4-LICENSE.md. Exact equations follow exact-helmert.ts, adapted from PROJ 9.5.1 helmert.cpp; see ../../PROJ-LICENSE.txt. Consume its prepared matrix without recomputing rotations.
+// SPDX-FileComment: Original whole-buffer dispatch. Small-angle equations follow the proj4js adaptation in datum.ts; see ../../../projection/PROJ4-LICENSE.md. Exact equations follow exact-helmert.ts, adapted from PROJ 9.5.1 helmert.cpp; see ../../../projection/PROJ-LICENSE.txt. Consume its prepared matrix without recomputing rotations.
 import type {NumericFlatOperation} from './numeric-flat';
 
 /** Setup-owned coefficients: small-angle rx/ry/rz or row-major exact matrix, then dx/dy/dz/scale. */

package/src/projection/experimental/inverse-spherical-latitude.ts

@@ -0,0 +1,19 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original Cartesian rotation identity; preserves the horizontal components near a pole instead of subtracting a rounded sine from one.
+
+export function inverseSphericalLatitude(
+  sin0: number,
+  cos0: number,
+  sinDistance: number,
+  cosDistance: number,
+  x: number,
+  y: number,
+  radius: number
+): number {
+  const east = (x * sinDistance) / radius;
+  const meridian = cosDistance * cos0 - (y * sinDistance * sin0) / radius;
+  const vertical = cosDistance * sin0 + (y * sinDistance * cos0) / radius;
+  return Math.atan2(vertical, Math.hypot(east, meridian));
+}

package/src/projection/experimental/kernels/aea.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Adapted from proj4js 2.22.0 with a roundoff-bounded authalic inverse at poles. See ../../../../projection/PROJ4-LICENSE.md for upstream attribution. Independent corrections acknowledge PROJ 9.5.1 src/projections/aea.cpp; see ../../../../projection/PROJ-LICENSE.txt.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {
@@ -27,7 +27,7 @@ export type State = KernelParameters & {
 import msfnz from '../common/msfnz';
 import qsfnz from '../common/qsfnz';
 import adjust_lon from '../common/adjust_lon';
-import asinz from '../common/asinz';
+import iqsfnz from '../common/iqsfnz';
 import {EPSLN} from '../common/constants';
 
 function initialize(state: State): void {
@@ -104,10 +104,16 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
   }
   con = (rh1 * state.ns0) / state.a;
   if (state.sphere) {
-    lat = Math.asin((state.c - con * con) / (2 * state.ns0));
+    const sinLatitude = (state.c - con * con) / (2 * state.ns0);
+    if (Math.abs(sinLatitude) > 1 + 8 * Number.EPSILON) return null;
+    lat = Math.asin(Math.max(-1, Math.min(1, sinLatitude)));
   } else {
     qs = (state.c - con * con) / state.ns0;
-    lat = phi1z(state.e3, qs);
+    const poleQ = qsfnz(state.e3, 1);
+    lat =
+      Math.abs(Math.abs(qs) - poleQ) <= 16 * Number.EPSILON * poleQ
+        ? (Math.sign(qs) * Math.PI) / 2
+        : iqsfnz(state.e3, qs);
   }
 
   lon = adjust_lon(theta / state.ns0 + state.long0, state.over);
@@ -116,33 +122,6 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
   return p;
 }
 
-/* Function to compute phi1, the latitude for the inverse of the
-   Albers Conical Equal-Area projection.
-------------------------------------------- */
-function phi1z(eccent: number, qs: number): number | null {
-  var sinphi, cosphi, con, com, dphi;
-  var phi = asinz(0.5 * qs);
-  if (eccent < EPSLN) {
-    return phi;
-  }
-
-  var eccnts = eccent * eccent;
-  for (var i = 1; i <= 25; i++) {
-    sinphi = Math.sin(phi);
-    cosphi = Math.cos(phi);
-    con = eccent * sinphi;
-    com = 1 - con * con;
-    dphi =
-      ((0.5 * com * com) / cosphi) *
-      (qs / (1 - eccnts) - sinphi / com + (0.5 / eccent) * Math.log((1 - con) / (1 + con)));
-    phi = phi + dphi;
-    if (Math.abs(dphi) <= 1e-7) {
-      return phi;
-    }
-  }
-  return null;
-}
-
 export function createState(base: KernelParameters): State {
   const state: State = {
     ...base,

package/src/projection/experimental/kernels/aeqd.ts

@@ -2,9 +2,10 @@
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
 // SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
-// SPDX-FileComment: Adapted from proj4js 2.22.0; polar meridian integration/inversion is original math.gl code. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Original math.gl inverse Cartesian rotation correction. Adapted from proj4js 2.22.0; polar meridian integration/inversion is original math.gl code. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
+import {inverseSphericalLatitude} from '../inverse-spherical-latitude';
 export type State = KernelParameters & {
   sin_p12: number;
   cos_p12: number;
@@ -12,7 +13,6 @@ export type State = KernelParameters & {
 };
 import adjust_lon from '../common/adjust_lon';
 import {HALF_PI, EPSLN} from '../common/constants';
-import asinz from '../common/asinz';
 import {meridianDistance, inverseMeridianDistance} from '../common/meridian-distance';
 import {vincentyDirect, vincentyInverse} from '../common/vincenty';
 
@@ -104,7 +104,7 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
     if (Math.abs(rh) <= EPSLN) {
       lat = state.lat0;
     } else {
-      lat = asinz(cosz * state.sin_p12 + (p.y * sinz * state.cos_p12) / rh);
+      lat = inverseSphericalLatitude(state.sin_p12, state.cos_p12, sinz, cosz, p.x, p.y, rh);
       con = Math.abs(state.lat0) - HALF_PI;
       if (Math.abs(con) <= EPSLN) {
         if (state.lat0 >= 0) {

package/src/projection/experimental/kernels/bonne.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/cass.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution. Corrected series signs and inverse refinement informed by PROJ 9.5.1 cass.cpp. See ../../../PROJ-LICENSE.txt. Newton solver implementation is original math.gl code.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution. Corrected series signs and inverse refinement informed by PROJ 9.5.1 cass.cpp. See ../../../../projection/PROJ-LICENSE.txt. Newton solver implementation is original math.gl code.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/cea.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {};

package/src/projection/experimental/kernels/eck6.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/eqdc.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Adapted from proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution. math.gl uses the higher-order meridional series shared with Cassini, rather than the upstream truncated e0/e1/e2/e3 series. Qualified against PROJ 9.5.1.
+// SPDX-FileComment: Adapted from proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution. math.gl uses the higher-order meridional series shared with Cassini, rather than the upstream truncated e0/e1/e2/e3 series. Qualified against PROJ 9.5.1.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/eqearth.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: Apache-2.0
 // SPDX-FileCopyrightText: 2018 Bernie Jenny, Monash University, Melbourne, Australia (Equal Earth, ported through proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0 lib/projections/eqearth.js. Modified to use explicit TypeScript state. Original Apache-2.0 notice retained below. See ../../../APACHE-2.0-LICENSE.txt and ../../../THIRD-PARTY-NOTICES.md.
+// SPDX-FileComment: Adapted from proj4js 2.22.0 lib/projections/eqearth.js with explicit TypeScript state and bounded inverse pole roundoff. Original Apache-2.0 notice retained below. See ../../../../projection/APACHE-2.0-LICENSE.txt and ../../../../projection/THIRD-PARTY-NOTICES.md. Independent corrections acknowledge PROJ 9.5.1 src/projections/eqearth.cpp and OSGeo/PROJ #1097; see ../../../../projection/PROJ-LICENSE.txt.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {
@@ -94,6 +94,11 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
     p.y /= state.rqda;
   }
 
+  const poleParam = Math.PI / 3;
+  const poleSq = poleParam * poleParam;
+  const poleY = poleParam * (A1 + A2 * poleSq + poleSq ** 3 * (A3 + A4 * poleSq));
+  if (Math.abs(p.y) > poleY + 8 * Number.EPSILON * poleY) return null;
+
   var EPS = 1e-9,
     NITER = 12,
     paramLat = p.y,
@@ -119,7 +124,12 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
   p.x =
     (M * p.x * (A1 + 3 * A2 * paramLatSq + paramLatPow6 * (7 * A3 + 9 * A4 * paramLatSq))) /
     Math.cos(paramLat);
-  p.y = Math.asin(Math.sin(paramLat) / M);
+  // Clamp only roundoff at the pole; points beyond the map remain invalid.
+  const sinBeta = Math.sin(paramLat) / M;
+  const poleRoundoff = 8 * Number.EPSILON;
+  if (Math.abs(sinBeta) > 1 + poleRoundoff) return null;
+  p.y =
+    Math.abs(sinBeta) >= 1 - poleRoundoff ? (Math.sign(sinBeta) * Math.PI) / 2 : Math.asin(sinBeta);
 
   if (state.es !== 0) {
     p.y = authlat(p.y, state.apa);

package/src/projection/experimental/kernels/equi.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Adapted from proj4js 2.22.0 with bounded inverse pole roundoff. See ../../../../projection/PROJ4-LICENSE.md for upstream attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {
@@ -39,6 +39,8 @@ export function inverse(state: State, p: Point): Point | null | undefined | numb
   p.x -= state.x0;
   p.y -= state.y0;
   var lat = p.y / state.a;
+  if (Math.abs(lat) > Math.PI / 2 && Math.abs(lat) <= Math.PI / 2 + 8 * Number.EPSILON)
+    lat = (Math.sign(lat) * Math.PI) / 2;
 
   var lon = adjust_lon(state.long0 + p.x / (state.a * Math.cos(state.lat0)), state.over);
   p.x = lon;

package/src/projection/experimental/kernels/etmerc.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/gauss.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/geos.ts

@@ -1,7 +1,7 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
 export type State = KernelParameters & {

package/src/projection/experimental/kernels/gnom.ts

@@ -1,10 +1,17 @@
 // math.gl
 // SPDX-License-Identifier: MIT
 // SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
-// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../PROJ4-LICENSE.md for the upstream license and attribution.
+// SPDX-FileComment: Spherical path adapted from proj4js 2.22.0; ellipsoidal path delegates to the attributed PROJ/GeographicLib adapter. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
 
 import type {KernelParameters, Point} from '../kernel';
+import {
+  prepareEllipsoidalGnomonic,
+  forwardEllipsoidalGnomonic,
+  inverseEllipsoidalGnomonic
+} from '../ellipsoidal-gnomonic';
+import type {EllipsoidalGnomonic} from '../ellipsoidal-gnomonic';
 export type State = KernelParameters & {
+  ellipsoid?: EllipsoidalGnomonic;
   sin_p14: number;
   cos_p14: number;
   infinity_dist: number;
@@ -16,6 +23,7 @@ import asinz from '../common/asinz';
 import {EPSLN} from '../common/constants';
 
 function initialize(state: State): void {
+  if (!state.sphere) state.ellipsoid = prepareEllipsoidalGnomonic(state);
   /* Place parameters in static storage for common use
       ------------------------------------------------- */
   state.sin_p14 = Math.sin(state.lat0);
@@ -28,6 +36,7 @@ function initialize(state: State): void {
 /* Gnomonic forward equations--mapping lat,long to x,y
     --------------------------------------------------- */
 export function forward(state: State, p: Point): Point | null | undefined | number {
+  if (state.ellipsoid) return forwardEllipsoidalGnomonic(state, state.ellipsoid, p);
   var sinphi, cosphi; /* sin and cos value        */
   var dlon; /* delta longitude value      */
   var coslon; /* cos of longitude        */
@@ -66,6 +75,7 @@ export function forward(state: State, p: Point): Point | null | undefined | numb
 }
 
 export function inverse(state: State, p: Point): Point | null | undefined | number {
+  if (state.ellipsoid) return inverseEllipsoidalGnomonic(state, state.ellipsoid, p);
   var rh; /* Rho */
   var sinc, cosc;
   var c;