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

package/src/projection/experimental/catalogue-kernel.ts

@@ -0,0 +1,65 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original adapter for kernels directly ported from proj4js 2.22.0.
+import {createProjection} from './mutable-projection';
+import {bindKernel, kernelParameters} from './kernel';
+import type {KernelParameters, Point} from './kernel';
+import {numberParameter} from './parameters';
+import type {ProjectionContext, ProjectionImplementation} from './types';
+
+export function catalogueParameters(context: ProjectionContext): KernelParameters {
+  const base = kernelParameters(context);
+  // False offsets are applied once by the adapter, including kernels that omit them upstream.
+  return {...base, x0: 0, y0: 0, lat_ts: Number.isNaN(base.lat_ts) ? 0 : base.lat_ts};
+}
+/** Bounds include false offsets that the numerical-kernel adapter applies externally. */
+export function catalogueBoundsParameters<State extends KernelParameters>(
+  state: State,
+  context: ProjectionContext
+): State {
+  return {
+    ...state,
+    x0: numberParameter(context.parameters, 'x_0', 0),
+    y0: numberParameter(context.parameters, 'y_0', 0)
+  };
+}
+export function bindCatalogueKernel<State>(
+  name: string,
+  context: ProjectionContext,
+  state: State,
+  forward: (state: State, point: Point) => Point | null | undefined | number,
+  inverse: (state: State, point: Point) => Point | null | undefined | number,
+  guard?: (longitude: number, latitude: number) => void
+): ProjectionImplementation {
+  const kernel = bindKernel(name, state, forward, inverse);
+  const x0 = numberParameter(context.parameters, 'x_0', 0),
+    y0 = numberParameter(context.parameters, 'y_0', 0);
+  return createProjection(
+    point => {
+      guard?.(point.x, point.y);
+      kernel.forwardInPlace(point);
+      point.x += x0;
+      point.y += y0;
+    },
+    point => {
+      point.x -= x0;
+      point.y -= y0;
+      kernel.inverseInPlace(point);
+    }
+  );
+}
+/** Reject invisible points instead of returning a finite horizon approximation. */
+export function horizonGuard(
+  base: KernelParameters,
+  minimumCosine = 0,
+  includeBoundary = false
+): (lon: number, lat: number) => void {
+  return (longitude, latitude) => {
+    const cosine =
+      Math.sin(base.lat0) * Math.sin(latitude) +
+      Math.cos(base.lat0) * Math.cos(latitude) * Math.cos(longitude - base.long0);
+    if (includeBoundary ? cosine < minimumCosine - 1e-14 : cosine <= minimumCosine + 1e-14)
+      throw new Error('Coordinate is outside the visible projection horizon');
+  };
+}

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

@@ -0,0 +1,11 @@
+// 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';
+import sign from './sign';
+
+export default function (x: number) {
+  return Math.abs(x) < HALF_PI ? x : x - sign(x) * Math.PI;
+}

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

@@ -0,0 +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 {TWO_PI, SPI} from './constants';
+import sign from './sign';
+
+export default function (x: number, skipAdjust: boolean) {
+  if (skipAdjust) {
+    return x;
+  }
+  return Math.abs(x) <= SPI ? x : x - sign(x) * TWO_PI;
+}

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

@@ -0,0 +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 hypot from './hypot';
+import log1py from './log1py';
+
+export default function (x: number) {
+  var y = Math.abs(x);
+  y = log1py(y * (1 + y / (hypot(1, y) + 1)));
+
+  return x < 0 ? -y : y;
+}

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

@@ -0,0 +1,11 @@
+// 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.
+
+export default function (x: number) {
+  if (Math.abs(x) > 1) {
+    x = x > 1 ? 1 : -1;
+  }
+  return Math.asin(x);
+}

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

@@ -0,0 +1,9 @@
+// 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.
+
+export default function authlat(beta: number, APA: readonly number[]) {
+  var t = beta + beta;
+  return beta + APA[0] * Math.sin(t) + APA[1] * Math.sin(t + t) + APA[2] * Math.sin(t + t + t);
+}

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

@@ -0,0 +1,25 @@
+// 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.
+
+var P00 = 0.33333333333333333333;
+var P01 = 0.17222222222222222222;
+var P02 = 0.10257936507936507936;
+var P10 = 0.06388888888888888888;
+var P11 = 0.06640211640211640211;
+var P20 = 0.01641501294219154443;
+
+export default function authset(es: number) {
+  var t;
+  var APA = [];
+  APA[0] = es * P00;
+  t = es * es;
+  APA[0] += t * P01;
+  APA[1] = t * P10;
+  t *= es;
+  APA[0] += t * P02;
+  APA[1] += t * P11;
+  APA[2] = t * P20;
+  return APA;
+}

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

@@ -0,0 +1,20 @@
+// 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.
+
+export default function (pp: readonly number[], arg_r: number) {
+  var r = 2 * Math.cos(arg_r);
+  var i = pp.length - 1;
+  var hr1 = pp[i];
+  var hr2 = 0;
+  var hr;
+
+  while (--i >= 0) {
+    hr = -hr2 + r * hr1 + pp[i];
+    hr2 = hr1;
+    hr1 = hr;
+  }
+
+  return Math.sin(arg_r) * hr;
+}

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

@@ -0,0 +1,43 @@
+// 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. Modified to write the complex result into caller-owned scratch storage.
+
+import sinh from './sinh';
+import cosh from './cosh';
+
+export default function (
+  pp: readonly number[],
+  arg_r: number,
+  arg_i: number,
+  output: {x: number; y: number}
+) {
+  var sin_arg_r = Math.sin(arg_r);
+  var cos_arg_r = Math.cos(arg_r);
+  var sinh_arg_i = sinh(arg_i);
+  var cosh_arg_i = cosh(arg_i);
+  var r = 2 * cos_arg_r * cosh_arg_i;
+  var i = -2 * sin_arg_r * sinh_arg_i;
+  var j = pp.length - 1;
+  var hr = pp[j];
+  var hi1 = 0;
+  var hr1 = 0;
+  var hi = 0;
+  var hr2;
+  var hi2;
+
+  while (--j >= 0) {
+    hr2 = hr1;
+    hi2 = hi1;
+    hr1 = hr;
+    hi1 = hi;
+    hr = -hr2 + r * hr1 - i * hi1 + pp[j];
+    hi = -hi2 + i * hr1 + r * hi1;
+  }
+
+  r = sin_arg_r * cosh_arg_i;
+  i = cos_arg_r * sinh_arg_i;
+
+  output.x = r * hr - i * hi;
+  output.y = r * hi + i * hr;
+}

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

@@ -0,0 +1,15 @@
+// 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.
+
+export const EPSLN = 1e-10;
+export const HALF_PI = Math.PI / 2;
+export const FORTPI = Math.PI / 4;
+export const TWO_PI = Math.PI * 2;
+// Preserve upstream tolerance at the antimeridian.
+export const SPI = 3.14159265359;
+
+export const D2R = Math.PI / 180;
+export const R2D = 180 / Math.PI;
+export const SEC_TO_RAD = Math.PI / (180 * 3600);

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

@@ -0,0 +1,10 @@
+// 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.
+
+export default function (x: number) {
+  var r = Math.exp(x);
+  r = (r + 1 / r) / 2;
+  return r;
+}

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

@@ -0,0 +1,8 @@
+// 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.
+
+export default function (x: number) {
+  return 1 - 0.25 * x * (1 + (x / 16) * (3 + 1.25 * x));
+}

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

@@ -0,0 +1,8 @@
+// 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.
+
+export default function (x: number) {
+  return 0.375 * x * (1 + 0.25 * x * (1 + 0.46875 * x));
+}

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

@@ -0,0 +1,8 @@
+// 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.
+
+export default function (x: number) {
+  return 0.05859375 * x * x * (1 + 0.75 * x);
+}

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

@@ -0,0 +1,8 @@
+// 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.
+
+export default function (x: number) {
+  return x * x * x * (35 / 3072);
+}

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

@@ -0,0 +1,10 @@
+// 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';
+
+export default function (x: number, L: number) {
+  return 2 * Math.atan(x * Math.exp(L)) - HALF_PI;
+}

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

@@ -0,0 +1,9 @@
+// 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.
+
+export default function (a: number, e: number, sinphi: number) {
+  var temp = e * sinphi;
+  return a / Math.sqrt(1 - temp * temp);
+}

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

@@ -0,0 +1,20 @@
+// 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.
+
+export default function (pp: readonly number[], B: number) {
+  var cos_2B = 2 * Math.cos(2 * B);
+  var i = pp.length - 1;
+  var h1 = pp[i];
+  var h2 = 0;
+  var h;
+
+  while (--i >= 0) {
+    h = -h2 + cos_2B * h1 + pp[i];
+    h2 = h1;
+    h1 = h;
+  }
+
+  return B + h * Math.sin(2 * B);
+}

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

@@ -0,0 +1,13 @@
+// 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.
+
+export default function (x: number, y: number) {
+  x = Math.abs(x);
+  y = Math.abs(y);
+  var a = Math.max(x, y);
+  var b = Math.min(x, y) / (a ? a : 1);
+
+  return a * Math.sqrt(1 + Math.pow(b, 2));
+}

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

@@ -0,0 +1,23 @@
+// 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.
+
+export default function (ml: number, e0: number, e1: number, e2: number, e3: number) {
+  var phi;
+  var dphi;
+
+  phi = ml / e0;
+  for (var i = 0; i < 15; i++) {
+    dphi =
+      (ml - (e0 * phi - e1 * Math.sin(2 * phi) + e2 * Math.sin(4 * phi) - e3 * Math.sin(6 * phi))) /
+      (e0 - 2 * e1 * Math.cos(2 * phi) + 4 * e2 * Math.cos(4 * phi) - 6 * e3 * Math.cos(6 * phi));
+    phi += dphi;
+    if (Math.abs(dphi) <= 0.0000000001) {
+      return phi;
+    }
+  }
+
+  // ..reportError("IMLFN-CONV:Latitude failed to converge after 15 iterations");
+  return NaN;
+}

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

@@ -0,0 +1,20 @@
+// 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 fL from './fL';
+
+export default function (eccent: number, ts: number) {
+  var phi = fL(1, ts);
+  var Iphi = 0;
+  var con = 0;
+  let iterations = 0;
+  do {
+    if (++iterations > 50) throw new Error('Isometric latitude inverse did not converge');
+    Iphi = phi;
+    con = eccent * Math.sin(Iphi);
+    phi = fL(Math.exp((eccent * Math.log((1 + con) / (1 - con))) / 2), ts);
+  } while (Math.abs(phi - Iphi) > 1.0e-12);
+  return phi;
+}

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

@@ -0,0 +1,45 @@
+// 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.
+
+import {HALF_PI} from './constants';
+
+// Original math.gl correction: the pole cutoff is roundoff-sized, not 1e-6
+// in authalic q (which collapsed latitudes as far as 0.05 degrees from a pole).
+
+export default function (eccent: number, q: number) {
+  var temp = 1 - ((1 - eccent * eccent) / (2 * eccent)) * Math.log((1 - eccent) / (1 + eccent));
+  const roundoff = 4 * Number.EPSILON * temp;
+  if (Math.abs(q) > temp + roundoff) return NaN;
+  if (Math.abs(Math.abs(q) - temp) <= roundoff) {
+    if (q < 0) {
+      return -1 * HALF_PI;
+    } else {
+      return HALF_PI;
+    }
+  }
+  // var phi = 0.5* q/(1-eccent*eccent);
+  var phi = Math.asin(0.5 * q);
+  var dphi;
+  var sin_phi;
+  var cos_phi;
+  var con;
+  for (var i = 0; i < 30; i++) {
+    sin_phi = Math.sin(phi);
+    cos_phi = Math.cos(phi);
+    con = eccent * sin_phi;
+    dphi =
+      (Math.pow(1 - con * con, 2) / (2 * cos_phi)) *
+      (q / (1 - eccent * eccent) -
+        sin_phi / (1 - con * con) +
+        (0.5 / eccent) * Math.log((1 - con) / (1 + con)));
+    phi += dphi;
+    if (Math.abs(dphi) <= 0.0000000001) {
+      return phi;
+    }
+  }
+
+  // console.log("IQSFN-CONV:Latitude failed to converge after 30 iterations");
+  return NaN;
+}

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

@@ -0,0 +1,21 @@
+// 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';
+
+export default function (eccent: number, phi: number, sinphi: number) {
+  if (Math.abs(phi) > HALF_PI) {
+    return Number.NaN;
+  }
+  if (phi === HALF_PI) {
+    return Number.POSITIVE_INFINITY;
+  }
+  if (phi === -1 * HALF_PI) {
+    return Number.NEGATIVE_INFINITY;
+  }
+
+  var con = eccent * sinphi;
+  return Math.log(Math.tan((HALF_PI + phi) / 2)) + (eccent * Math.log((1 - con) / (1 + con))) / 2;
+}

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

@@ -0,0 +1,11 @@
+// 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.
+
+export default function (x: number) {
+  var y = 1 + x;
+  var z = y - 1;
+
+  return z === 0 ? x : (x * Math.log(y)) / z;
+}

package/src/projection/experimental/common/meridian-distance.ts

@@ -0,0 +1,45 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original math.gl Gauss-Legendre integration of the meridional radius.
+
+function radius(es: number, latitude: number): number {
+  const d = 1 - es * Math.sin(latitude) ** 2;
+  return (1 - es) / (d * Math.sqrt(d));
+}
+function pair(es: number, midpoint: number, half: number, node: number): number {
+  return radius(es, midpoint - half * node) + radius(es, midpoint + half * node);
+}
+/** Signed meridian distance. Four eight-point panels avoid the legacy e^6 truncation. */
+export function meridianDistance(a: number, es: number, start: number, end: number): number {
+  const half = (end - start) / 8;
+  let distance = 0;
+  for (let i = 0; i < 4; i++) {
+    const midpoint = start + (2 * i + 1) * half;
+    distance +=
+      half *
+      (0.362683783378362 * pair(es, midpoint, half, 0.1834346424956498) +
+        0.3137066458778873 * pair(es, midpoint, half, 0.525532409916329) +
+        0.2223810344533745 * pair(es, midpoint, half, 0.7966664774136267) +
+        0.1012285362903763 * pair(es, midpoint, half, 0.9602898564975363));
+  }
+  return a * distance;
+}
+/** Bounded inverse of the same integral; NaN marks an invalid radius or failed solve. */
+export function inverseMeridianDistance(a: number, es: number, distance: number): number {
+  const pole = Math.PI / 2;
+  const limit = meridianDistance(a, es, 0, pole);
+  if (!Number.isFinite(distance) || Math.abs(distance) > limit) return NaN;
+  let lower = -pole,
+    upper = pole,
+    latitude = (distance / limit) * pole;
+  for (let i = 0; i < 48; i++) {
+    const residual = meridianDistance(a, es, 0, latitude) - distance;
+    if (Math.abs(residual) <= 8 * Number.EPSILON * a) return latitude;
+    if (residual > 0) upper = latitude;
+    else lower = latitude;
+    const next = latitude - residual / (a * radius(es, latitude));
+    latitude = next > lower && next < upper ? next : (lower + upper) / 2;
+  }
+  return NaN;
+}

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

@@ -0,0 +1,8 @@
+// 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.
+
+export default function (e0: number, e1: number, e2: number, e3: number, phi: number) {
+  return e0 * phi - e1 * Math.sin(2 * phi) + e2 * Math.sin(4 * phi) - e3 * Math.sin(6 * phi);
+}

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

@@ -0,0 +1,9 @@
+// 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.
+
+export default function (eccent: number, sinphi: number, cosphi: number) {
+  var con = eccent * sinphi;
+  return cosphi / Math.sqrt(1 - con * con);
+}

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

@@ -0,0 +1,22 @@
+// 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';
+
+export default function (eccent: number, ts: number) {
+  var eccnth = 0.5 * eccent;
+  var con, dphi;
+  var phi = HALF_PI - 2 * Math.atan(ts);
+  for (var i = 0; i <= 15; i++) {
+    con = eccent * Math.sin(phi);
+    dphi = HALF_PI - 2 * Math.atan(ts * Math.pow((1 - con) / (1 + con), eccnth)) - phi;
+    phi += dphi;
+    if (Math.abs(dphi) <= 0.0000000001) {
+      return phi;
+    }
+  }
+  // console.log("phi2z has NoConvergence");
+  return -9999;
+}

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

@@ -0,0 +1,29 @@
+// 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.
+
+var C00 = 1;
+var C02 = 0.25;
+var C04 = 0.046875;
+var C06 = 0.01953125;
+var C08 = 0.01068115234375;
+var C22 = 0.75;
+var C44 = 0.46875;
+var C46 = 0.01302083333333333333;
+var C48 = 0.00712076822916666666;
+var C66 = 0.36458333333333333333;
+var C68 = 0.00569661458333333333;
+var C88 = 0.3076171875;
+
+export default function (es: number) {
+  var en = [];
+  en[0] = C00 - es * (C02 + es * (C04 + es * (C06 + es * C08)));
+  en[1] = es * (C22 - es * (C04 + es * (C06 + es * C08)));
+  var t = es * es;
+  en[2] = t * (C44 - es * (C46 + es * C48));
+  t *= es;
+  en[3] = t * (C66 - es * C68);
+  en[4] = t * es * C88;
+  return en;
+}

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

@@ -0,0 +1,28 @@
+// 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 pj_mlfn from './pj_mlfn';
+import {EPSLN} from './constants';
+
+var MAX_ITER = 20;
+
+export default function (arg: number, es: number, en: readonly number[]) {
+  var k = 1 / (1 - es);
+  var phi = arg;
+  for (var i = MAX_ITER; i; --i) {
+    /* rarely goes over 2 iterations */
+    var s = Math.sin(phi);
+    var t = 1 - es * s * s;
+    // t = this.pj_mlfn(phi, s, Math.cos(phi), en) - arg;
+    // phi -= t * (t * Math.sqrt(t)) * k;
+    t = (pj_mlfn(phi, s, Math.cos(phi), en) - arg) * (t * Math.sqrt(t)) * k;
+    phi -= t;
+    if (Math.abs(t) < EPSLN) {
+      return phi;
+    }
+  }
+  // ..reportError("cass:pj_inv_mlfn: Convergence error");
+  return phi;
+}

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

@@ -0,0 +1,10 @@
+// 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.
+
+export default function (phi: number, sphi: number, cphi: number, en: readonly number[]) {
+  cphi *= sphi;
+  sphi *= sphi;
+  return en[0] * phi - cphi * (en[1] + sphi * (en[2] + sphi * (en[3] + sphi * en[4])));
+}

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

@@ -0,0 +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.
+
+export default function (eccent: number, sinphi: number) {
+  var con;
+  if (eccent > 1.0e-7) {
+    con = eccent * sinphi;
+    return (
+      (1 - eccent * eccent) *
+      (sinphi / (1 - con * con) - (0.5 / eccent) * Math.log((1 - con) / (1 + con)))
+    );
+  } else {
+    return 2 * sinphi;
+  }
+}