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

package/src/local-reference-frame.ts

@@ -0,0 +1,320 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original geodetic frame and pose composition over math.gl local-frame kernels.
+
+import {Matrix4, Quaternion, Vector3, config} from '@math.gl/core';
+import {
+  createLocalFrameBasis,
+  eastNorthUpBasis,
+  localFrameToMatrix
+} from '@math.gl/core/local-frame';
+import type {LocalFrameAxis} from '@math.gl/core/local-frame';
+import type {NumericArray} from '@math.gl/core/types';
+import {Ellipsoid} from './ellipsoid';
+
+export type LocalReferenceFrameOptions = {
+  ellipsoid?: Ellipsoid;
+  /** Right-handed signed geographic axes. Default east/north/up. */
+  axes?: readonly [LocalFrameAxis, LocalFrameAxis, LocalFrameAxis];
+};
+/** Radians: Rz(-heading) Ry(-pitch) Rx(roll) in the selected local axes.
+ * Zero attitude aligns model XYZ with local XYZ (ENU by default).
+ */
+export type HeadingPitchRoll = {heading: number; pitch: number; roll: number};
+
+type Coordinate = Readonly<NumericArray>;
+
+/** A snapshot of a rigid local geodetic frame on a sphere or oblate spheroid.
+ * Origins/positions use fixed Cartesian axis units; vectors retain their input units.
+ * The basis follows geodetic longitude/latitude, including elevated origins.
+ */
+export class LocalReferenceFrame {
+  private readonly matrix: readonly number[];
+  private readonly up: readonly number[];
+
+  constructor(origin: Coordinate, options: LocalReferenceFrameOptions = {}) {
+    finiteCoordinate(origin);
+    const ellipsoid = options.ellipsoid || Ellipsoid.WGS84;
+    const {x: a, y, z: b} = ellipsoid.radii;
+    if (!(a > 0 && a === y && b > 0 && b <= a && Number.isFinite(a))) {
+      throw new RangeError('LocalReferenceFrame requires a positive sphere or oblate spheroid');
+    }
+    const llh = ellipsoid.cartesianToCartographic(origin);
+    if (!llh) throw new RangeError('Origin has no unique geodetic frame');
+    const basis = createLocalFrameBasis();
+    const radians = config._cartographicRadians ? 1 : Math.PI / 180;
+    if (!eastNorthUpBasis(llh[0] * radians, llh[1] * radians, basis))
+      throw new RangeError('Cannot resolve local geodetic basis');
+    const axes = options.axes || ['east', 'north', 'up'];
+    const matrix: number[] = [];
+    if (!localFrameToMatrix(basis, {x: origin[0], y: origin[1], z: origin[2]}, ...axes, matrix)) {
+      throw new RangeError('Expected three right-handed geographic axes');
+    }
+    this.matrix = Object.freeze(matrix);
+    this.up = Object.freeze([basis.upX, basis.upY, basis.upZ]);
+  }
+
+  /** Independent copies: callers cannot mutate the cached frame. */
+  get localToFixedMatrix(): Matrix4 {
+    return new Matrix4(this.matrix);
+  }
+
+  get fixedToLocalMatrix(): Matrix4 {
+    return new Matrix4(this.matrix).invert();
+  }
+
+  localToFixedPoint(point: Coordinate, result = new Vector3()): Vector3 {
+    return this.transform(point, result, false, true);
+  }
+
+  fixedToLocalPoint(point: Coordinate, result = new Vector3()): Vector3 {
+    return this.transform(point, result, true, true);
+  }
+
+  localToFixedVector(vector: Coordinate, result = new Vector3()): Vector3 {
+    return this.transform(vector, result, false, false);
+  }
+
+  fixedToLocalVector(vector: Coordinate, result = new Vector3()): Vector3 {
+    return this.transform(vector, result, true, false);
+  }
+
+  /** Packed XYZ points. Exact in-place operation is supported; other overlap is rejected.
+   * All validation and representability checks finish before writing caller storage.
+   */
+  localToFixedPointsFlat(input: Coordinate, result?: Float64Array): Float64Array {
+    return this.transformFlat(input, result, false, true);
+  }
+
+  fixedToLocalPointsFlat(input: Coordinate, result?: Float64Array): Float64Array {
+    return this.transformFlat(input, result, true, true);
+  }
+
+  localToFixedVectorsFlat(input: Coordinate, result?: Float64Array): Float64Array {
+    return this.transformFlat(input, result, false, false);
+  }
+
+  fixedToLocalVectorsFlat(input: Coordinate, result?: Float64Array): Float64Array {
+    return this.transformFlat(input, result, true, false);
+  }
+
+  /** Model-to-fixed affine pose, with unit scale. Angles are radians. */
+  headingPitchRollToFixedFrame(attitude: HeadingPitchRoll, result = new Matrix4()): Matrix4 {
+    const {heading, pitch, roll} = attitude;
+    if (![heading, pitch, roll].every(Number.isFinite)) throw new RangeError('Nonfinite attitude');
+    return result.copy(this.matrix).rotateZ(-heading).rotateY(-pitch).rotateX(roll);
+  }
+
+  headingPitchRollQuaternion(attitude: HeadingPitchRoll, result = new Quaternion()): Quaternion {
+    const m = this.headingPitchRollToFixedFrame(attitude);
+    return result.fromMatrix3([m[0], m[1], m[2], m[4], m[5], m[6], m[8], m[9], m[10]]).normalize();
+  }
+
+  /** Extract attitude from a rigid model-to-fixed pose. Translation is ignored.
+   * Pitch is canonicalized to [-pi/2, pi/2]. At gimbal lock, roll is zero.
+   */
+  fixedFrameToHeadingPitchRoll(pose: Coordinate): HeadingPitchRoll {
+    validateRigidPose(pose);
+    // Rotation extraction must not multiply irrelevant (possibly huge) translations.
+    const local = new Matrix4();
+    const m = this.matrix;
+    for (let column = 0; column < 3; column++) {
+      for (let row = 0; row < 3; row++) {
+        local[column * 4 + row] =
+          m[row * 4] * pose[column * 4] +
+          m[row * 4 + 1] * pose[column * 4 + 1] +
+          m[row * 4 + 2] * pose[column * 4 + 2];
+      }
+    }
+    const pitchY = Math.atan2(-local[2], Math.hypot(local[0], local[1]));
+    const locked = Math.hypot(local[0], local[1]) < 1e-12;
+    const roll = locked ? 0 : Math.atan2(local[6], local[10]);
+    const yaw = locked ? Math.atan2(-local[4], local[5]) : Math.atan2(local[1], local[0]);
+    return {heading: -yaw, pitch: -pitchY, roll};
+  }
+
+  /** Fixed-space velocity defines model +X. +Z is as close as possible to geodetic up.
+   * Returns undefined without writing result for zero or nearly vertical velocity.
+   * Stateless: callers own temporal continuity across these singularities.
+   */
+  velocityQuaternion(velocity: Coordinate, result = new Quaternion()): Quaternion | undefined {
+    finiteCoordinate(velocity);
+    const scale = Math.max(Math.abs(velocity[0]), Math.abs(velocity[1]), Math.abs(velocity[2]));
+    if (scale === 0) return undefined;
+    const x = new Vector3(velocity[0] / scale, velocity[1] / scale, velocity[2] / scale);
+    x.scale(1 / Math.hypot(x[0], x[1], x[2]));
+    const y = new Vector3(this.up).cross(x);
+    if (y.len() < 1e-12) return undefined;
+    y.normalize();
+    const z = new Vector3(x).cross(y).normalize();
+    return result.fromMatrix3([...x, ...y, ...z]).normalize();
+  }
+
+  private transformFlat(
+    input: Coordinate,
+    result: Float64Array | undefined,
+    inverse: boolean,
+    position: boolean
+  ): Float64Array {
+    if (
+      !Number.isSafeInteger(input.length) ||
+      input.length % 3 ||
+      (result && (!(result instanceof Float64Array) || result.length !== input.length))
+    )
+      throw new RangeError('Expected packed XYZ and matching Float64 output');
+    if (
+      result &&
+      input !== result &&
+      ArrayBuffer.isView(input) &&
+      input.buffer === result.buffer &&
+      input.byteOffset < result.byteOffset + result.byteLength &&
+      result.byteOffset < input.byteOffset + input.byteLength
+    )
+      throw new RangeError('Packed frame buffers overlap');
+    let bounded =
+      !position ||
+      Math.max(Math.abs(this.matrix[12]), Math.abs(this.matrix[13]), Math.abs(this.matrix[14])) <=
+        Number.MAX_VALUE / 8;
+    for (let i = 0; i < input.length; i++) {
+      if (!Number.isFinite(input[i])) throw new RangeError('Expected finite packed coordinates');
+      bounded &&= Math.abs(input[i]) <= Number.MAX_VALUE / 8;
+    }
+    // Ordinary bounded coordinates cannot overflow a rigid 3D transform. Only
+    // extreme batches need a preflight pass to preserve atomic failure.
+    const scratch = bounded ? undefined : new Vector3();
+    if (scratch)
+      for (let i = 0; i < input.length; i += 3)
+        this.transform(
+          scratch.set(input[i], input[i + 1], input[i + 2]),
+          scratch,
+          inverse,
+          position
+        );
+    const output = result || new Float64Array(input.length),
+      [xx, xy, xz, , yx, yy, yz, , zx, zy, zz, , ox, oy, oz] = this.matrix;
+    for (let i = 0; i < input.length; i += 3) {
+      const x = input[i] - (inverse && position ? ox : 0),
+        y = input[i + 1] - (inverse && position ? oy : 0),
+        z = input[i + 2] - (inverse && position ? oz : 0);
+      if (scratch) {
+        this.transform(
+          scratch.set(input[i], input[i + 1], input[i + 2]),
+          scratch,
+          inverse,
+          position
+        );
+        output.set(scratch, i);
+      } else {
+        output[i] = inverse
+          ? xx * x + xy * y + xz * z
+          : xx * x + yx * y + zx * z + (position ? ox : 0);
+        output[i + 1] = inverse
+          ? yx * x + yy * y + yz * z
+          : xy * x + yy * y + zy * z + (position ? oy : 0);
+        output[i + 2] = inverse
+          ? zx * x + zy * y + zz * z
+          : xz * x + yz * y + zz * z + (position ? oz : 0);
+      }
+    }
+    return output;
+  }
+
+  private transform(
+    point: Coordinate,
+    result: Vector3,
+    inverse: boolean,
+    position: boolean
+  ): Vector3 {
+    finiteCoordinate(point);
+    const m = this.matrix;
+    const x = point[0] - (inverse && position ? m[12] : 0);
+    const y = point[1] - (inverse && position ? m[13] : 0);
+    const z = point[2] - (inverse && position ? m[14] : 0);
+    let ox = inverse
+      ? m[0] * x + m[1] * y + m[2] * z
+      : m[0] * x + m[4] * y + m[8] * z + (position ? m[12] : 0);
+    let oy = inverse
+      ? m[4] * x + m[5] * y + m[6] * z
+      : m[1] * x + m[5] * y + m[9] * z + (position ? m[13] : 0);
+    let oz = inverse
+      ? m[8] * x + m[9] * y + m[10] * z
+      : m[2] * x + m[6] * y + m[10] * z + (position ? m[14] : 0);
+    if (!Number.isFinite(ox) || !Number.isFinite(oy) || !Number.isFinite(oz)) {
+      // Preserve the ordinary subtract-first path (small offsets at large
+      // origins). Rescale only when an intermediate overflowed before cancellation.
+      const scale = Math.max(
+        Math.abs(point[0]),
+        Math.abs(point[1]),
+        Math.abs(point[2]),
+        position ? Math.abs(m[12]) : 0,
+        position ? Math.abs(m[13]) : 0,
+        position ? Math.abs(m[14]) : 0
+      );
+      const sx = point[0] / scale - (inverse && position ? m[12] / scale : 0);
+      const sy = point[1] / scale - (inverse && position ? m[13] / scale : 0);
+      const sz = point[2] / scale - (inverse && position ? m[14] / scale : 0);
+      ox =
+        (inverse
+          ? m[0] * sx + m[1] * sy + m[2] * sz
+          : m[0] * sx + m[4] * sy + m[8] * sz + (position ? m[12] / scale : 0)) * scale;
+      oy =
+        (inverse
+          ? m[4] * sx + m[5] * sy + m[6] * sz
+          : m[1] * sx + m[5] * sy + m[9] * sz + (position ? m[13] / scale : 0)) * scale;
+      oz =
+        (inverse
+          ? m[8] * sx + m[9] * sy + m[10] * sz
+          : m[2] * sx + m[6] * sy + m[10] * sz + (position ? m[14] / scale : 0)) * scale;
+    }
+    if (!Number.isFinite(ox) || !Number.isFinite(oy) || !Number.isFinite(oz))
+      throw new RangeError('Frame transformation overflow');
+    return result.set(ox, oy, oz);
+  }
+}
+
+function finiteCoordinate(point: Coordinate): void {
+  if (
+    point.length !== 3 ||
+    !Number.isFinite(point[0]) ||
+    !Number.isFinite(point[1]) ||
+    !Number.isFinite(point[2])
+  ) {
+    throw new RangeError('Expected three finite coordinates');
+  }
+}
+
+function validateRigidPose(pose: Coordinate): void {
+  if (
+    pose.length !== 16 ||
+    !Array.from(pose).every(Number.isFinite) ||
+    pose[3] !== 0 ||
+    pose[7] !== 0 ||
+    pose[11] !== 0 ||
+    pose[15] !== 1
+  ) {
+    throw new RangeError('Expected a finite affine pose');
+  }
+  const x = new Vector3(pose[0], pose[1], pose[2]);
+  const y = new Vector3(pose[4], pose[5], pose[6]);
+  const z = new Vector3(pose[8], pose[9], pose[10]);
+  // Reject oversized/scaled columns before vector cross products can overflow,
+  // including when math.gl debug validation is enabled.
+  for (const axis of [x, y, z]) {
+    if (Math.abs(Math.hypot(axis[0], axis[1], axis[2]) - 1) > 1e-6) {
+      throw new RangeError('Pose rotation must be orthonormal and right-handed');
+    }
+  }
+  const errors = [
+    Math.abs(x.len() - 1),
+    Math.abs(y.len() - 1),
+    Math.abs(z.len() - 1),
+    Math.abs(x.dot(y)),
+    Math.abs(x.dot(z)),
+    Math.abs(y.dot(z)),
+    Math.abs(new Vector3(x).cross(y).dot(z) - 1)
+  ];
+  if (!errors.every(Number.isFinite) || Math.max(...errors) > 1e-6) {
+    throw new RangeError('Pose rotation must be orthonormal and right-handed');
+  }
+}

package/src/make-obb-from-region.ts

@@ -5,7 +5,7 @@
 // SPDX-FileComment: Derived from Cesium. See the repository LICENSE for upstream attribution and Apache-2.0 terms.
 
 // This file is derived from the Cesium math library under Apache 2 license
-// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md
+// See ../LICENSE for the retained Cesium notice and Apache-2.0 terms.
 
 import {Matrix3, Matrix4, Vector2, Vector3, degrees, radians, _MathUtils} from '@math.gl/core';
 import {OrientedBoundingBox, Plane} from '@math.gl/geometry/culling';

package/src/projection/experimental/authalic-latitude.ts

@@ -0,0 +1,49 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original stable evaluation and inversion of the standard authalic integral. No upstream implementation copied.
+
+/** qp - |q(phi)|, evaluated without subtracting almost equal pole integrals. */
+export function authalicPoleDifference(e: number, latitude: number): number {
+  const sin = Math.abs(Math.sin(latitude));
+  const cos = Math.cos(latitude);
+  const delta = (cos * cos) / (1 + sin);
+  if (e === 0) return 2 * delta;
+  const es = e * e;
+  return (
+    (delta * (1 + es * sin)) / (1 - es * sin * sin) +
+    ((1 - es) / e) * Math.atanh((e * delta) / (1 - es * sin))
+  );
+}
+
+/** Stable authalic latitude, including the cosine retained close to a pole. */
+export function authalicLatitude(e: number, qp: number, latitude: number): number {
+  if (latitude === 0 || e === 0) return latitude;
+  const fraction = authalicPoleDifference(e, latitude) / qp;
+  return Math.atan2(Math.sign(latitude) * (1 - fraction), Math.sqrt(fraction * (2 - fraction)));
+}
+
+/** Recover geodetic latitude from authalic latitude, preserving pole complements. */
+export function inverseAuthalicLatitude(e: number, qp: number, beta: number): number {
+  if (e === 0 || beta === 0 || Math.abs(beta) === Math.PI / 2) return beta;
+  const sinBeta = Math.abs(Math.sin(beta));
+  const cosBeta = Math.cos(beta);
+  const target = (qp * cosBeta * cosBeta) / (1 + sinBeta);
+  let latitude = Math.abs(beta);
+  const es = e * e;
+  let lower = 0;
+  let upper = Math.PI / 2;
+  for (let i = 0; i < 48; i++) {
+    const sin = Math.sin(latitude);
+    const denominator = 1 - es * sin * sin;
+    const derivative = (2 * (1 - es) * Math.cos(latitude)) / (denominator * denominator);
+    const residual = authalicPoleDifference(e, latitude) - target;
+    const step = residual / derivative;
+    if (Math.abs(step) < 2e-15) return Math.sign(beta) * latitude;
+    if (residual > 0) lower = latitude;
+    else upper = latitude;
+    const candidate = latitude + step;
+    latitude = candidate > lower && candidate < upper ? candidate : (lower + upper) / 2;
+  }
+  throw new Error('Authalic latitude inverse did not converge');
+}

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

@@ -1,11 +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.
+// SPDX-FileComment: Adapted from proj4js 2.22.0; corrected the inclusive pole boundary. See ../../../../projection/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;
+  return Math.abs(x) <= HALF_PI ? x : x - sign(x) * Math.PI;
 }

package/src/projection/experimental/common/adjust_lon.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 {TWO_PI, SPI} from './constants';
 import sign from './sign';

package/src/projection/experimental/common/asinhy.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 hypot from './hypot';
 import log1py from './log1py';

package/src/projection/experimental/common/asinz.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) {
   if (Math.abs(x) > 1) {

package/src/projection/experimental/common/authlat.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 authlat(beta: number, APA: readonly number[]) {
   var t = beta + beta;

package/src/projection/experimental/common/authset.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 P00 = 0.33333333333333333333;
 var P01 = 0.17222222222222222222;

package/src/projection/experimental/common/clens.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 (pp: readonly number[], arg_r: number) {
   var r = 2 * Math.cos(arg_r);

package/src/projection/experimental/common/clens_cmplx.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. Modified to write the complex result into caller-owned scratch storage.
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0. See ../../../../projection/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';

package/src/projection/experimental/common/constants.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 const EPSLN = 1e-10;
 export const HALF_PI = Math.PI / 2;

package/src/projection/experimental/common/cosh.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/e0fn.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 1 - 0.25 * x * (1 + (x / 16) * (3 + 1.25 * x));

package/src/projection/experimental/common/e1fn.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 0.375 * x * (1 + 0.25 * x * (1 + 0.46875 * x));

package/src/projection/experimental/common/e2fn.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 0.05859375 * x * x * (1 + 0.75 * x);

package/src/projection/experimental/common/e3fn.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 * x * x * (35 / 3072);

package/src/projection/experimental/common/fL.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/gN.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 (a: number, e: number, sinphi: number) {
   var temp = e * sinphi;

package/src/projection/experimental/common/gatg.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 (pp: readonly number[], B: number) {
   var cos_2B = 2 * Math.cos(2 * B);

package/src/projection/experimental/common/hypot.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, y: number) {
   x = Math.abs(x);

package/src/projection/experimental/common/imlfn.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 (ml: number, e0: number, e1: number, e2: number, e3: number) {
   var phi;

package/src/projection/experimental/common/invlatiso.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 fL from './fL';
 

package/src/projection/experimental/common/iqsfnz.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.
+// SPDX-FileComment: Adapted from 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/latiso.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/log1py.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 y = 1 + x;

package/src/projection/experimental/common/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 (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

@@ -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, cosphi: number) {
   var con = eccent * sinphi;