@math.gl/geospatial 5.0.0-alpha.9 → 5.0.0-beta.2

package/dist/ellipsoid-helpers/tangent-plane-numerics.js

@@ -0,0 +1,88 @@
+import { ratio } from "./exact-nearest-surface.js";
+const SPLITTER = 134217729;
+function productError(a, b, product) {
+    const ca = SPLITTER * a, cb = SPLITTER * b;
+    const ah = ca - (ca - a), al = a - ah;
+    const bh = cb - (cb - b), bl = b - bh;
+    return ah * bh - product + ah * bl + al * bh + al * bl;
+}
+function sumError(a, b, sum) {
+    const v = sum - a;
+    return a - (sum - v) + (b - v);
+}
+/** Exact power-of-two scaling keeps products finite without rounding ratios. */
+function binaryScale(value) {
+    return 2 ** Math.min(1023, Math.floor(Math.log2(value)));
+}
+export function compensatedDot3(x, y, z, a, b, c) {
+    const magnitude = Math.max(Math.abs(x), Math.abs(y), Math.abs(z));
+    if (magnitude === 0)
+        return 0;
+    const scale = binaryScale(magnitude);
+    const sx = x / scale, sy = y / scale, sz = z / scale;
+    const p = sx * a, q = sy * b, r = sz * c;
+    const s = p + q, t = s + r;
+    return ((t +
+        (sumError(p, q, s) +
+            sumError(s, r, t) +
+            productError(sx, a, p) +
+            productError(sy, b, q) +
+            productError(sz, c, r))) *
+        scale);
+}
+/** Rare exact-sign fallback for a normal line near an ellipsoid silhouette.
+ * The stored binary64 origin/basis and local coordinates define the affine line.
+ */
+export function exactNormalIntersection(radii, origin, east, north, normal, point, sphereRadius) {
+    const values = sphereRadius
+        ? [
+            sphereRadius,
+            sphereRadius,
+            sphereRadius,
+            sphereRadius,
+            0,
+            0,
+            0,
+            1,
+            0,
+            0,
+            0,
+            1,
+            1,
+            0,
+            0,
+            point[0],
+            point[1]
+        ]
+        : [...radii, ...origin, ...east, ...north, ...normal, ...point];
+    const v = values.map(value => {
+        const [whole, fraction = ''] = Math.abs(value).toString(2).split('.');
+        return (((value < 0 ? -1n : 1n) * BigInt('0b' + whole + fraction)) << BigInt(1074 - fraction.length));
+    });
+    const unit = 1n << 1074n, q = 1n << 64n;
+    const factors = [v[1] * v[2], v[0] * v[2], v[0] * v[1]];
+    const p = [0, 1, 2].map(i => (v[3 + i] * unit + v[6 + i] * v[15] + v[9 + i] * v[16]) * factors[i]);
+    const n = v.slice(12, 15).map((x, i) => x * factors[i]);
+    const dot = (x, y) => x.reduce((s, v, i) => s + v * y[i], 0n);
+    const a = dot(n, n), b = dot(n, p), c = dot(p, p) - (v[0] * v[1] * v[2] * unit) ** 2n;
+    const discriminant = b * b - a * c;
+    if (a === 0n || discriminant < 0n)
+        return undefined;
+    const value = discriminant << 128n;
+    let root = value === 0n ? 0n : 1n << BigInt(Math.ceil(value.toString(2).length / 2));
+    if (root)
+        for (;;) {
+            const next = (root + value / root) >> 1n;
+            if (next >= root)
+                break;
+            root = next;
+        }
+    const output = p.map((x, i) => signedRatio(x * a * q + n[i] * (root - b * q), a * q * factors[i]));
+    return sphereRadius
+        ? output.map((_, i) => east[i] * output[1] + north[i] * output[2] + normal[i] * output[0])
+        : output;
+}
+function signedRatio(n, d) {
+    return Math.sign(Number(n)) * ratio(n < 0n ? -n : n, d, -2148);
+}
+//# sourceMappingURL=tangent-plane-numerics.js.map

package/dist/ellipsoid-helpers/tangent-plane-numerics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"tangent-plane-numerics.js","sourceRoot":"","sources":["../../src/ellipsoid-helpers/tangent-plane-numerics.ts"],"names":[],"mappings":"AAKA,OAAO,EAAC,KAAK,EAAC,mCAAgC;AAC9C,MAAM,QAAQ,GAAG,SAAS,CAAC;AAC3B,SAAS,YAAY,CAAC,CAAS,EAAE,CAAS,EAAE,OAAe;IACzD,MAAM,EAAE,GAAG,QAAQ,GAAG,CAAC,EACrB,EAAE,GAAG,QAAQ,GAAG,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EACtB,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EACtB,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC;IACd,OAAO,EAAE,GAAG,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AACzD,CAAC;AACD,SAAS,QAAQ,CAAC,CAAS,EAAE,CAAS,EAAE,GAAW;IACjD,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IAClB,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;AACjC,CAAC;AACD,gFAAgF;AAChF,SAAS,WAAW,CAAC,KAAa;IAChC,OAAO,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAC3D,CAAC;AACD,MAAM,UAAU,eAAe,CAC7B,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS,EACT,CAAS;IAET,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,IAAI,SAAS,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,WAAW,CAAC,SAAS,CAAC,CAAC;IACrC,MAAM,EAAE,GAAG,CAAC,GAAG,KAAK,EAClB,EAAE,GAAG,CAAC,GAAG,KAAK,EACd,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC;IACjB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,EACd,CAAC,GAAG,EAAE,GAAG,CAAC,EACV,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACb,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,EACb,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACZ,OAAO,CACL,CAAC,CAAC;QACA,CAAC,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YAChB,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACjB,YAAY,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;YACtB,YAAY,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC;YACtB,YAAY,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5B,KAAK,CACN,CAAC;AACJ,CAAC;AACD;;GAEG;AACH,MAAM,UAAU,uBAAuB,CACrC,KAAwB,EACxB,MAAyB,EACzB,IAAuB,EACvB,KAAwB,EACxB,MAAyB,EACzB,KAA6B,EAC7B,YAAqB;IAErB,MAAM,MAAM,GAAG,YAAY;QACzB,CAAC,CAAC;YACE,YAAY;YACZ,YAAY;YACZ,YAAY;YACZ,YAAY;YACZ,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,CAAC;YACD,KAAK,CAAC,CAAC,CAAC;YACR,KAAK,CAAC,CAAC,CAAC;SACT;QACH,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,GAAG,MAAM,EAAE,GAAG,IAAI,EAAE,GAAG,KAAK,EAAE,GAAG,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC;IAClE,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE;QAC3B,MAAM,CAAC,KAAK,EAAE,QAAQ,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACtE,OAAO,CACL,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,IAAI,GAAG,KAAK,GAAG,QAAQ,CAAC,CAAC,IAAI,MAAM,CAAC,IAAI,GAAG,QAAQ,CAAC,MAAM,CAAC,CAC7F,CAAC;IACJ,CAAC,CAAC,CAAC;IACH,MAAM,IAAI,GAAG,EAAE,IAAI,KAAK,EACtB,CAAC,GAAG,EAAE,IAAI,GAAG,CAAC;IAChB,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CACrB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAC1E,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,MAAM,GAAG,GAAG,CAAC,CAAW,EAAE,CAAW,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IAClF,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EACjB,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,EACb,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;IACpD,MAAM,YAAY,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACnC,IAAI,CAAC,KAAK,EAAE,IAAI,YAAY,GAAG,EAAE;QAAE,OAAO,SAAS,CAAC;IACpD,MAAM,KAAK,GAAG,YAAY,IAAI,IAAI,CAAC;IACnC,IAAI,IAAI,GAAG,KAAK,KAAK,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IACrF,IAAI,IAAI;QACN,SAAS,CAAC;YACR,MAAM,IAAI,GAAG,CAAC,IAAI,GAAG,KAAK,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;YACzC,IAAI,IAAI,IAAI,IAAI;gBAAE,MAAM;YACxB,IAAI,GAAG,IAAI,CAAC;QACd,CAAC;IACH,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAC5B,WAAW,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CACnE,CAAC;IACF,OAAO,YAAY;QACjB,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAC1F,CAAC,CAAC,MAAM,CAAC;AACb,CAAC;AAED,SAAS,WAAW,CAAC,CAAS,EAAE,CAAS;IACvC,OAAO,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC;AACjE,CAAC"}

package/dist/ellipsoid-metrics.d.ts

@@ -0,0 +1,31 @@
+import type { GeodesicParameters } from "./ellipsoid-geodesic.js";
+/** All latitude values are degrees, including isometric latitude (radians * 180/pi). */
+export type AuxiliaryLatitude = 'geocentric' | 'parametric' | 'rectifying' | 'authalic' | 'conformal' | 'isometric';
+/** Metrics on a revolution ellipsoid. Qualified for 0.5 <= b/a <= 2.
+ * Axes use caller length units; areas/volumes use their square/cube.
+ * Separate from Ellipsoid so Cartesian and projection users pay no bundle cost.
+ */
+export declare class EllipsoidMetrics {
+    readonly semiMajorAxis: number;
+    readonly semiMinorAxis: number;
+    readonly surfaceArea: number;
+    readonly volume: number;
+    readonly quarterMeridian: number;
+    private readonly ratio;
+    private readonly e2;
+    private readonly qPole;
+    private readonly quarter;
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }: GeodesicParameters);
+    meridianDistance(latitude: number): number;
+    meridionalCurvatureRadius(latitude: number): number;
+    transverseCurvatureRadius(latitude: number): number;
+    normalCurvatureRadius(latitude: number, azimuth: number): number;
+    latitudeToAuxiliary(latitude: number, kind: AuxiliaryLatitude): number;
+    latitudeFromAuxiliary(latitude: number, kind: AuxiliaryLatitude): number;
+    private eccentricIntegral;
+    private meridianIntegrand;
+    private meridian;
+}
+/** @internal Shared bounded quadrature for original revolution integrals. */
+export declare function integrate(f: (t: number) => number, lo: number, hi: number): number;
+//# sourceMappingURL=ellipsoid-metrics.d.ts.map

package/dist/ellipsoid-metrics.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-metrics.d.ts","sourceRoot":"","sources":["../src/ellipsoid-metrics.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,kBAAkB,EAAC,gCAA6B;AAE7D,wFAAwF;AACxF,MAAM,MAAM,iBAAiB,GACzB,YAAY,GACZ,YAAY,GACZ,YAAY,GACZ,UAAU,GACV,WAAW,GACX,WAAW,CAAC;AAqBhB;;;GAGG;AACH,qBAAa,gBAAgB;IAC3B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAC;IAC/B,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,QAAQ,CAAC,eAAe,EAAE,MAAM,CAAC;IACjC,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAS;IAC/B,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAS;IAC5B,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAS;IAC/B,OAAO,CAAC,QAAQ,CAAC,OAAO,CAAS;gBAErB,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAC,EAAE,kBAAkB;IAsBpE,gBAAgB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAK1C,yBAAyB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAKnD,yBAAyB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;IAMnD,qBAAqB,CAAC,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAG,MAAM;IAUhE,mBAAmB,CAAC,QAAQ,EAAE,MAAM,EAAE,IAAI,EAAE,iBAAiB,GAAG,MAAM;IAyDtE,qBAAqB,CAAC,QAAQ,EAAE,MAAM,EAAE,IAAI,EAAE,iBAAiB,GAAG,MAAM;IA2BxE,OAAO,CAAC,iBAAiB;IAKzB,OAAO,CAAC,iBAAiB;IAIzB,OAAO,CAAC,QAAQ;CAKjB;AAiBD,6EAA6E;AAC7E,wBAAgB,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAWlF"}

package/dist/ellipsoid-metrics.js

@@ -0,0 +1,186 @@
+const DEG = Math.PI / 180;
+// Positive half of the 16-point Gauss-Legendre rule.
+const NODES = [
+    0.09501250983763744, 0.2816035507792589, 0.4580167776572274, 0.6178762444026438,
+    0.755404408355003, 0.8656312023878318, 0.9445750230732326, 0.9894009349916499
+];
+const WEIGHTS = [
+    0.1894506104550685, 0.1826034150449236, 0.16915651939500254, 0.14959598881657673,
+    0.12462897125553388, 0.09515851168249279, 0.06225352393864789, 0.027152459411754095
+];
+const KINDS = [
+    'geocentric',
+    'parametric',
+    'rectifying',
+    'authalic',
+    'conformal',
+    'isometric'
+];
+/** Metrics on a revolution ellipsoid. Qualified for 0.5 <= b/a <= 2.
+ * Axes use caller length units; areas/volumes use their square/cube.
+ * Separate from Ellipsoid so Cartesian and projection users pay no bundle cost.
+ */
+export class EllipsoidMetrics {
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }) {
+        const r = b / a;
+        if (!(Number.isFinite(a) && Number.isFinite(b) && a > 0 && b > 0 && r >= 0.5 && r <= 2)) {
+            throw new Error('Ellipsoid metrics require finite positive axes with 0.5 <= b/a <= 2');
+        }
+        this.semiMajorAxis = a;
+        this.semiMinorAxis = b;
+        this.ratio = r;
+        this.e2 = (1 - r) * (1 + r);
+        this.qPole = 1 + r * r * this.eccentricIntegral(1);
+        this.quarter = integrate(t => this.meridianIntegrand(t), 0, Math.PI / 2);
+        this.quarterMeridian = a * this.quarter;
+        this.surfaceArea = 2 * Math.PI * a * a * this.qPole;
+        this.volume = ((4 * Math.PI) / 3) * a * b * a;
+        if (![this.quarterMeridian, this.surfaceArea, this.volume].every(v => Number.isFinite(v) && v > 0)) {
+            throw new Error('Ellipsoid metrics exceed representable length/area/volume range');
+        }
+        Object.freeze(this);
+    }
+    meridianDistance(latitude) {
+        validateLatitude(latitude);
+        return Math.sign(latitude) * this.semiMajorAxis * this.meridian(Math.abs(latitude));
+    }
+    meridionalCurvatureRadius(latitude) {
+        validateLatitude(latitude);
+        return this.semiMajorAxis * this.meridianIntegrand(latitude * DEG);
+    }
+    transverseCurvatureRadius(latitude) {
+        validateLatitude(latitude);
+        const [s, c] = sincos(latitude);
+        return this.semiMajorAxis / Math.hypot(c, this.ratio * s);
+    }
+    normalCurvatureRadius(latitude, azimuth) {
+        if (!Number.isFinite(azimuth))
+            throw new Error('Azimuth must be finite');
+        const [s, c] = sincos(azimuth % 180);
+        return (1 /
+            ((c * c) / this.meridionalCurvatureRadius(latitude) +
+                (s * s) / this.transverseCurvatureRadius(latitude)));
+    }
+    latitudeToAuxiliary(latitude, kind) {
+        validateKind(kind);
+        validateLatitude(latitude);
+        if (latitude === 0)
+            return latitude;
+        if (Math.abs(latitude) === 90)
+            return kind === 'isometric' ? Math.sign(latitude) * Infinity : latitude;
+        if (this.ratio === 1 && kind !== 'isometric')
+            return latitude;
+        const sign = Math.sign(latitude);
+        const phi = Math.abs(latitude);
+        const [s, c] = sincos(phi);
+        switch (kind) {
+            case 'geocentric':
+                return sign * angle(this.ratio * this.ratio * s, c);
+            case 'parametric':
+                return sign * angle(this.ratio * s, c);
+            case 'rectifying':
+                return (sign *
+                    (phi <= 45
+                        ? (this.meridian(phi) / this.quarter) * 90
+                        : 90 -
+                            (integrate(t => this.meridianIntegrand(Math.PI / 2 - t), 0, (90 - phi) * DEG) /
+                                this.quarter) *
+                                90));
+            case 'authalic': {
+                if (phi <= 45) {
+                    const q = this.ratio ** 2 * (s / (c * c + (this.ratio * s) ** 2) + this.eccentricIntegral(s));
+                    return (sign * Math.asin(q / this.qPole)) / DEG;
+                }
+                const cap = integrate(t => {
+                    const st = Math.sin(t), ct = Math.cos(t);
+                    const d = st * st + (this.ratio * ct) ** 2;
+                    return (2 * this.ratio ** 2 * st) / (d * d);
+                }, 0, (90 - phi) * DEG);
+                return sign * (90 - (2 * Math.asin(Math.sqrt(cap / (2 * this.qPole)))) / DEG);
+            }
+            default: {
+                const psi = Math.asinh(s / c) - this.e2 * this.eccentricIntegral(s);
+                return (sign *
+                    (kind === 'isometric'
+                        ? psi / DEG
+                        : psi < 1
+                            ? Math.atan(Math.sinh(psi)) / DEG
+                            : 90 - (2 * Math.atan(Math.exp(-psi))) / DEG));
+            }
+        }
+    }
+    latitudeFromAuxiliary(latitude, kind) {
+        validateKind(kind);
+        if (kind === 'isometric' && Math.abs(latitude) === Infinity)
+            return Math.sign(latitude) * 90;
+        if (kind !== 'isometric')
+            validateLatitude(latitude);
+        else if (!Number.isFinite(latitude))
+            throw new Error('Isometric latitude must be finite or a pole infinity');
+        if (latitude === 0)
+            return latitude;
+        if (kind !== 'isometric' && Math.abs(latitude) === 90)
+            return latitude;
+        if (this.ratio === 1 && kind !== 'isometric')
+            return latitude;
+        const sign = Math.sign(latitude), target = Math.abs(latitude);
+        if (kind === 'geocentric' || kind === 'parametric') {
+            const [s, c] = sincos(target);
+            return sign * angle(s, c * (kind === 'geocentric' ? this.ratio ** 2 : this.ratio));
+        }
+        // Monotonic, bounded solve in degrees; preserves input complements near poles.
+        let lo = 0, hi = 90;
+        for (let i = 0; i < 64; ++i) {
+            const mid = (lo + hi) / 2;
+            if (mid === lo || mid === hi)
+                break;
+            if (this.latitudeToAuxiliary(mid, kind) < target)
+                lo = mid;
+            else
+                hi = mid;
+        }
+        return sign * ((lo + hi) / 2);
+    }
+    eccentricIntegral(s) {
+        const e = Math.sqrt(Math.abs(this.e2));
+        return e === 0 ? s : (this.e2 > 0 ? Math.atanh(e * s) : Math.atan(e * s)) / e;
+    }
+    meridianIntegrand(t) {
+        return this.ratio ** 2 / Math.hypot(Math.cos(t), this.ratio * Math.sin(t)) ** 3;
+    }
+    meridian(phi) {
+        return phi <= 45
+            ? integrate(t => this.meridianIntegrand(t), 0, phi * DEG)
+            : this.quarter - integrate(t => this.meridianIntegrand(Math.PI / 2 - t), 0, (90 - phi) * DEG);
+    }
+}
+function validateLatitude(value) {
+    if (!(Number.isFinite(value) && Math.abs(value) <= 90))
+        throw new Error('Latitude must be in [-90, 90] degrees');
+}
+function validateKind(kind) {
+    if (!KINDS.includes(kind))
+        throw new Error('Unknown auxiliary latitude');
+}
+function sincos(degrees) {
+    const a = Math.abs(degrees);
+    if (a <= 45)
+        return [Math.sin(degrees * DEG), Math.cos(degrees * DEG)];
+    return [Math.sign(degrees) * Math.cos((90 - a) * DEG), Math.sin((90 - a) * DEG)];
+}
+function angle(s, c) {
+    return s > c ? 90 - Math.atan2(c, s) / DEG : Math.atan2(s, c) / DEG;
+}
+/** @internal Shared bounded quadrature for original revolution integrals. */
+export function integrate(f, lo, hi) {
+    let sum = 0;
+    for (let j = 0; j < 8; ++j) {
+        const h = (hi - lo) / 16, mid = lo + (2 * j + 1) * h;
+        let part = 0;
+        for (let i = 0; i < 8; ++i)
+            part += WEIGHTS[i] * (f(mid - h * NODES[i]) + f(mid + h * NODES[i]));
+        sum += h * part;
+    }
+    return sum;
+}
+//# sourceMappingURL=ellipsoid-metrics.js.map

package/dist/ellipsoid-metrics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-metrics.js","sourceRoot":"","sources":["../src/ellipsoid-metrics.ts"],"names":[],"mappings":"AAeA,MAAM,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;AAC1B,qDAAqD;AACrD,MAAM,KAAK,GAAG;IACZ,mBAAmB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;IAC/E,iBAAiB,EAAE,kBAAkB,EAAE,kBAAkB,EAAE,kBAAkB;CAC9E,CAAC;AACF,MAAM,OAAO,GAAG;IACd,kBAAkB,EAAE,kBAAkB,EAAE,mBAAmB,EAAE,mBAAmB;IAChF,mBAAmB,EAAE,mBAAmB,EAAE,mBAAmB,EAAE,oBAAoB;CACpF,CAAC;AACF,MAAM,KAAK,GAAsB;IAC/B,YAAY;IACZ,YAAY;IACZ,YAAY;IACZ,UAAU;IACV,WAAW;IACX,WAAW;CACZ,CAAC;AAEF;;;GAGG;AACH,MAAM,OAAO,gBAAgB;IAW3B,YAAY,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAqB;QAClE,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC;YACxF,MAAM,IAAI,KAAK,CAAC,qEAAqE,CAAC,CAAC;QACzF,CAAC;QACD,IAAI,CAAC,aAAa,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,aAAa,GAAG,CAAC,CAAC;QACvB,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5B,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;QACnD,IAAI,CAAC,OAAO,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACzE,IAAI,CAAC,eAAe,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC;QACxC,IAAI,CAAC,WAAW,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACpD,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9C,IACE,CAAC,CAAC,IAAI,CAAC,eAAe,EAAE,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAC9F,CAAC;YACD,MAAM,IAAI,KAAK,CAAC,iEAAiE,CAAC,CAAC;QACrF,CAAC;QACD,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IACtB,CAAC;IAED,gBAAgB,CAAC,QAAgB;QAC/B,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAC3B,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;IACtF,CAAC;IAED,yBAAyB,CAAC,QAAgB;QACxC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAC3B,OAAO,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,iBAAiB,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;IACrE,CAAC;IAED,yBAAyB,CAAC,QAAgB;QACxC,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAC3B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC;QAChC,OAAO,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,qBAAqB,CAAC,QAAgB,EAAE,OAAe;QACrD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,wBAAwB,CAAC,CAAC;QACzE,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC;QACrC,OAAO,CACL,CAAC;YACD,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,yBAAyB,CAAC,QAAQ,CAAC;gBACjD,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,yBAAyB,CAAC,QAAQ,CAAC,CAAC,CACtD,CAAC;IACJ,CAAC;IAED,mBAAmB,CAAC,QAAgB,EAAE,IAAuB;QAC3D,YAAY,CAAC,IAAI,CAAC,CAAC;QACnB,gBAAgB,CAAC,QAAQ,CAAC,CAAC;QAC3B,IAAI,QAAQ,KAAK,CAAC;YAAE,OAAO,QAAQ,CAAC;QACpC,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE;YAC3B,OAAO,IAAI,KAAK,WAAW,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC;QAC1E,IAAI,IAAI,CAAC,KAAK,KAAK,CAAC,IAAI,IAAI,KAAK,WAAW;YAAE,OAAO,QAAQ,CAAC;QAC9D,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACjC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;QAC/B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC3B,QAAQ,IAAI,EAAE,CAAC;YACb,KAAK,YAAY;gBACf,OAAO,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YACtD,KAAK,YAAY;gBACf,OAAO,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YACzC,KAAK,YAAY;gBACf,OAAO,CACL,IAAI;oBACJ,CAAC,GAAG,IAAI,EAAE;wBACR,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE;wBAC1C,CAAC,CAAC,EAAE;4BACF,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,iBAAiB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;gCAC3E,IAAI,CAAC,OAAO,CAAC;gCACb,EAAE,CAAC,CACV,CAAC;YACJ,KAAK,UAAU,CAAC,CAAC,CAAC;gBAChB,IAAI,GAAG,IAAI,EAAE,EAAE,CAAC;oBACd,MAAM,CAAC,GACL,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC;oBACtF,OAAO,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC;gBAClD,CAAC;gBACD,MAAM,GAAG,GAAG,SAAS,CACnB,CAAC,CAAC,EAAE;oBACF,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EACpB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBACnB,MAAM,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC;oBAC3C,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9C,CAAC,EACD,CAAC,EACD,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,CACjB,CAAC;gBACF,OAAO,IAAI,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;YAChF,CAAC;YACD,OAAO,CAAC,CAAC,CAAC;gBACR,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;gBACpE,OAAO,CACL,IAAI;oBACJ,CAAC,IAAI,KAAK,WAAW;wBACnB,CAAC,CAAC,GAAG,GAAG,GAAG;wBACX,CAAC,CAAC,GAAG,GAAG,CAAC;4BACP,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG;4BACjC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAClD,CAAC;YACJ,CAAC;QACH,CAAC;IACH,CAAC;IAED,qBAAqB,CAAC,QAAgB,EAAE,IAAuB;QAC7D,YAAY,CAAC,IAAI,CAAC,CAAC;QACnB,IAAI,IAAI,KAAK,WAAW,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,QAAQ;YAAE,OAAO,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,CAAC;QAC7F,IAAI,IAAI,KAAK,WAAW;YAAE,gBAAgB,CAAC,QAAQ,CAAC,CAAC;aAChD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YACjC,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;QAC1E,IAAI,QAAQ,KAAK,CAAC;YAAE,OAAO,QAAQ,CAAC;QACpC,IAAI,IAAI,KAAK,WAAW,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,EAAE;YAAE,OAAO,QAAQ,CAAC;QACvE,IAAI,IAAI,CAAC,KAAK,KAAK,CAAC,IAAI,IAAI,KAAK,WAAW;YAAE,OAAO,QAAQ,CAAC;QAC9D,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,EAC9B,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;QAC9B,IAAI,IAAI,KAAK,YAAY,IAAI,IAAI,KAAK,YAAY,EAAE,CAAC;YACnD,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;YAC9B,OAAO,IAAI,GAAG,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,KAAK,YAAY,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QACrF,CAAC;QACD,+EAA+E;QAC/E,IAAI,EAAE,GAAG,CAAC,EACR,EAAE,GAAG,EAAE,CAAC;QACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;YAC5B,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YAC1B,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE;gBAAE,MAAM;YACpC,IAAI,IAAI,CAAC,mBAAmB,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,MAAM;gBAAE,EAAE,GAAG,GAAG,CAAC;;gBACtD,EAAE,GAAG,GAAG,CAAC;QAChB,CAAC;QACD,OAAO,IAAI,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;IAChC,CAAC;IAEO,iBAAiB,CAAC,CAAS;QACjC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACvC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAChF,CAAC;IAEO,iBAAiB,CAAC,CAAS;QACjC,OAAO,IAAI,CAAC,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IAClF,CAAC;IAEO,QAAQ,CAAC,GAAW;QAC1B,OAAO,GAAG,IAAI,EAAE;YACd,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,iBAAiB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC;YACzD,CAAC,CAAC,IAAI,CAAC,OAAO,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,iBAAiB,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;IAClG,CAAC;CACF;AAED,SAAS,gBAAgB,CAAC,KAAa;IACrC,IAAI,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;QACpD,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;AAC7D,CAAC;AACD,SAAS,YAAY,CAAC,IAAuB;IAC3C,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,CAAC;AAC3E,CAAC;AACD,SAAS,MAAM,CAAC,OAAe;IAC7B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;IAC5B,IAAI,CAAC,IAAI,EAAE;QAAE,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC,CAAC;IACvE,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AACnF,CAAC;AACD,SAAS,KAAK,CAAC,CAAS,EAAE,CAAS;IACjC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;AACtE,CAAC;AACD,6EAA6E;AAC7E,MAAM,UAAU,SAAS,CAAC,CAAwB,EAAE,EAAU,EAAE,EAAU;IACxE,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,EACtB,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC;YACxB,IAAI,IAAI,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACvE,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC;IAClB,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/ellipsoid-occluder.d.ts

@@ -0,0 +1,28 @@
+/** The limb is center + axis1*cos(t) + axis2*sin(t), in input Cartesian units. */
+export type EllipsoidHorizon = {
+    center: number[];
+    axis1: number[];
+    axis2: number[];
+};
+/** Axis-aligned ellipsoid queries, independent of viewport and vector classes. */
+export declare class EllipsoidOccluder {
+    readonly radii: readonly number[];
+    readonly center: readonly number[];
+    constructor(radii: Readonly<ArrayLike<number>>, center?: Readonly<ArrayLike<number>>);
+    /** Nonnegative ray parameter interval, including tangency and inside-origin rays.
+     * Direction need not be normalized. Undefined is a true miss; never clamps a miss
+     * to the limb. Result is untouched on miss. Distances equal parameters only for
+     * unit directions. Zero direction is a miss.
+     */
+    intersectRay(origin: Readonly<ArrayLike<number>>, direction: Readonly<ArrayLike<number>>, result?: number[]): number[] | undefined;
+    /** True if the finite camera-to-point segment enters the solid interior.
+     * Surface contact/tangency alone is not occlusion. Works for elevated points;
+     * unlike a facing-normal test it does not discard visible points beyond the
+     * surface horizon. tolerance is a nonnegative scaled-space squared-radius margin.
+     */
+    isPointOccluded(camera: Readonly<ArrayLike<number>>, point: Readonly<ArrayLike<number>>, tolerance?: number): boolean;
+    /** Exact apparent limb ellipse for an exterior camera; undefined on/inside surface. */
+    getHorizon(camera: Readonly<ArrayLike<number>>): EllipsoidHorizon | undefined;
+    private scalePoint;
+}
+//# sourceMappingURL=ellipsoid-occluder.d.ts.map

package/dist/ellipsoid-occluder.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-occluder.d.ts","sourceRoot":"","sources":["../src/ellipsoid-occluder.ts"],"names":[],"mappings":"AAKA,kFAAkF;AAClF,MAAM,MAAM,gBAAgB,GAAG;IAC7B,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,KAAK,EAAE,MAAM,EAAE,CAAC;IAChB,KAAK,EAAE,MAAM,EAAE,CAAC;CACjB,CAAC;AAEF,kFAAkF;AAClF,qBAAa,iBAAiB;IAC5B,QAAQ,CAAC,KAAK,EAAE,SAAS,MAAM,EAAE,CAAC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,MAAM,EAAE,CAAC;gBAEvB,KAAK,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAAE,MAAM,GAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAa;IAQ/F;;;;OAIG;IACH,YAAY,CACV,MAAM,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EACnC,SAAS,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EACtC,MAAM,GAAE,MAAM,EAAW,GACxB,MAAM,EAAE,GAAG,SAAS;IA8BvB;;;;OAIG;IACH,eAAe,CACb,MAAM,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EACnC,KAAK,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,EAClC,SAAS,SAAQ,GAChB,OAAO;IAuBV,uFAAuF;IACvF,UAAU,CAAC,MAAM,EAAE,QAAQ,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC,GAAG,gBAAgB,GAAG,SAAS;IAuB7E,OAAO,CAAC,UAAU;CAKnB"}

package/dist/ellipsoid-occluder.js

@@ -0,0 +1,112 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original scaled-space analytic geometry. Consumer requirements: https://github.com/keplergl/kepler.gl/pull/3557 and https://github.com/keplergl/kepler.gl/pull/3593. No implementation code copied.
+/** Axis-aligned ellipsoid queries, independent of viewport and vector classes. */
+export class EllipsoidOccluder {
+    constructor(radii, center = [0, 0, 0]) {
+        validateVector(radii);
+        validateVector(center);
+        if (Array.from(radii).some(value => value <= 0))
+            throw new RangeError('Radii must be positive');
+        this.radii = Object.freeze(Array.from(radii));
+        this.center = Object.freeze(Array.from(center));
+    }
+    /** Nonnegative ray parameter interval, including tangency and inside-origin rays.
+     * Direction need not be normalized. Undefined is a true miss; never clamps a miss
+     * to the limb. Result is untouched on miss. Distances equal parameters only for
+     * unit directions. Zero direction is a miss.
+     */
+    intersectRay(origin, direction, result = [0, 0]) {
+        validateVector(origin);
+        validateVector(direction);
+        const [x, y, z] = this.scalePoint(origin);
+        const dx = direction[0] / this.radii[0];
+        const dy = direction[1] / this.radii[1];
+        const dz = direction[2] / this.radii[2];
+        // Normalize in scaled space: avoid a tiny quadratic coefficient for ECEF rays.
+        const length = Math.hypot(dx, dy, dz);
+        if (length === 0 || !Number.isFinite(length))
+            return undefined;
+        const ux = dx / length, uy = dy / length, uz = dz / length;
+        const along = -(x * ux + y * uy + z * uz);
+        const perpendicular = Math.hypot(y * uz - z * uy, z * ux - x * uz, x * uy - y * ux);
+        if (perpendicular > 1)
+            return undefined;
+        const half = Math.sqrt((1 - perpendicular) * (1 + perpendicular));
+        let near = along - half;
+        let far = along + half;
+        // Recover the small root by the root product rather than subtracting large terms.
+        const radius = Math.hypot(x, y, z);
+        if (!Number.isFinite(radius))
+            throw new RangeError('Scaled origin magnitude must be finite');
+        if (half > 0 && along > 0 && far !== 0)
+            near = ((radius - 1) / far) * (radius + 1);
+        if (half > 0 && along < 0 && near !== 0)
+            far = ((radius - 1) / near) * (radius + 1);
+        if (far < 0)
+            return undefined;
+        result[0] = Math.max(0, Math.min(near, far)) / length;
+        result[1] = Math.max(0, far) / length;
+        return result;
+    }
+    /** True if the finite camera-to-point segment enters the solid interior.
+     * Surface contact/tangency alone is not occlusion. Works for elevated points;
+     * unlike a facing-normal test it does not discard visible points beyond the
+     * surface horizon. tolerance is a nonnegative scaled-space squared-radius margin.
+     */
+    isPointOccluded(camera, point, tolerance = 1e-12) {
+        validateVector(camera);
+        validateVector(point);
+        if (!Number.isFinite(tolerance) || tolerance < 0 || tolerance >= 1)
+            throw new RangeError('Tolerance must be in [0, 1)');
+        const a = this.scalePoint(camera);
+        const b = this.scalePoint(point);
+        const dx = b[0] - a[0], dy = b[1] - a[1], dz = b[2] - a[2];
+        const length = Math.hypot(dx, dy, dz);
+        if (!Number.isFinite(length))
+            throw new RangeError('Segment magnitude must be finite');
+        if (length === 0)
+            return Math.hypot(...a) ** 2 < 1 - tolerance;
+        const ux = dx / length, uy = dy / length, uz = dz / length;
+        const distance = Math.max(0, Math.min(length, -(a[0] * ux + a[1] * uy + a[2] * uz)));
+        return (Math.hypot(a[0] + distance * ux, a[1] + distance * uy, a[2] + distance * uz) ** 2 <
+            1 - tolerance);
+    }
+    /** Exact apparent limb ellipse for an exterior camera; undefined on/inside surface. */
+    getHorizon(camera) {
+        validateVector(camera);
+        const scaled = this.scalePoint(camera);
+        const distance = Math.hypot(...scaled);
+        if (!Number.isFinite(distance))
+            throw new RangeError('Camera magnitude must be finite');
+        if (distance <= 1)
+            return undefined;
+        const normal = scaled.map(value => value / distance);
+        // Choose a basis away from parallel to the normal, including polar views.
+        const reference = Math.abs(normal[2]) < 0.9 ? [0, 0, 1] : [0, 1, 0];
+        const u = cross(normal, reference);
+        const uLength = Math.hypot(...u);
+        const v = cross(normal, u.map(value => value / uLength));
+        const radius = Math.sqrt((1 - 1 / distance) * (1 + 1 / distance));
+        return {
+            center: normal.map((value, i) => this.center[i] + (this.radii[i] * value) / distance),
+            axis1: u.map((value, i) => (this.radii[i] * radius * value) / uLength),
+            axis2: v.map((value, i) => this.radii[i] * radius * value)
+        };
+    }
+    scalePoint(point) {
+        const result = this.radii.map((radius, i) => (point[i] - this.center[i]) / radius);
+        if (!result.every(Number.isFinite))
+            throw new RangeError('Scaled coordinates must be finite');
+        return result;
+    }
+}
+function validateVector(vector) {
+    if (vector.length !== 3 || !Array.from(vector).every(Number.isFinite))
+        throw new RangeError('Expected three finite Cartesian components');
+}
+function cross(a, b) {
+    return [a[1] * b[2] - a[2] * b[1], a[2] * b[0] - a[0] * b[2], a[0] * b[1] - a[1] * b[0]];
+}
+//# sourceMappingURL=ellipsoid-occluder.js.map

package/dist/ellipsoid-occluder.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-occluder.js","sourceRoot":"","sources":["../src/ellipsoid-occluder.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,4DAA4D;AAC5D,wNAAwN;AASxN,kFAAkF;AAClF,MAAM,OAAO,iBAAiB;IAI5B,YAAY,KAAkC,EAAE,SAAsC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;QAC7F,cAAc,CAAC,KAAK,CAAC,CAAC;QACtB,cAAc,CAAC,MAAM,CAAC,CAAC;QACvB,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,IAAI,CAAC,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,wBAAwB,CAAC,CAAC;QAChG,IAAI,CAAC,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QAC9C,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;IAClD,CAAC;IAED;;;;OAIG;IACH,YAAY,CACV,MAAmC,EACnC,SAAsC,EACtC,SAAmB,CAAC,CAAC,EAAE,CAAC,CAAC;QAEzB,cAAc,CAAC,MAAM,CAAC,CAAC;QACvB,cAAc,CAAC,SAAS,CAAC,CAAC;QAC1B,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACxC,+EAA+E;QAC/E,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,IAAI,MAAM,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YAAE,OAAO,SAAS,CAAC;QAC/D,MAAM,EAAE,GAAG,EAAE,GAAG,MAAM,EACpB,EAAE,GAAG,EAAE,GAAG,MAAM,EAChB,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC;QACnB,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAC1C,MAAM,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QACpF,IAAI,aAAa,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QACxC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC;QAClE,IAAI,IAAI,GAAG,KAAK,GAAG,IAAI,CAAC;QACxB,IAAI,GAAG,GAAG,KAAK,GAAG,IAAI,CAAC;QACvB,kFAAkF;QAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACnC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,wCAAwC,CAAC,CAAC;QAC7F,IAAI,IAAI,GAAG,CAAC,IAAI,KAAK,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC;YAAE,IAAI,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACnF,IAAI,IAAI,GAAG,CAAC,IAAI,KAAK,GAAG,CAAC,IAAI,IAAI,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QACpF,IAAI,GAAG,GAAG,CAAC;YAAE,OAAO,SAAS,CAAC;QAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QACtD,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC;QACtC,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;;OAIG;IACH,eAAe,CACb,MAAmC,EACnC,KAAkC,EAClC,SAAS,GAAG,KAAK;QAEjB,cAAc,CAAC,MAAM,CAAC,CAAC;QACvB,cAAc,CAAC,KAAK,CAAC,CAAC;QACtB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,SAAS,GAAG,CAAC,IAAI,SAAS,IAAI,CAAC;YAChE,MAAM,IAAI,UAAU,CAAC,6BAA6B,CAAC,CAAC;QACtD,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QAClC,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EACpB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAChB,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACtC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,kCAAkC,CAAC,CAAC;QACvF,IAAI,MAAM,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC;QAC/D,MAAM,EAAE,GAAG,EAAE,GAAG,MAAM,EACpB,EAAE,GAAG,EAAE,GAAG,MAAM,EAChB,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC;QACnB,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACrF,OAAO,CACL,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,CAAC,IAAI,CAAC;YACjF,CAAC,GAAG,SAAS,CACd,CAAC;IACJ,CAAC;IAED,uFAAuF;IACvF,UAAU,CAAC,MAAmC;QAC5C,cAAc,CAAC,MAAM,CAAC,CAAC;QACvB,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QACvC,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;QACvC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;QACxF,IAAI,QAAQ,IAAI,CAAC;YAAE,OAAO,SAAS,CAAC;QACpC,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,QAAQ,CAAC,CAAC;QACrD,0EAA0E;QAC1E,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACpE,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACnC,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,KAAK,CACb,MAAM,EACN,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,KAAK,GAAG,OAAO,CAAC,CAChC,CAAC;QACF,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC;QAClE,OAAO;YACL,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,QAAQ,CAAC;YACrF,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC,GAAG,OAAO,CAAC;YACtE,KAAK,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;SAC3D,CAAC;IACJ,CAAC;IAEO,UAAU,CAAC,KAAkC;QACnD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;QACnF,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;YAAE,MAAM,IAAI,UAAU,CAAC,mCAAmC,CAAC,CAAC;QAC9F,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAED,SAAS,cAAc,CAAC,MAAmC;IACzD,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QACnE,MAAM,IAAI,UAAU,CAAC,4CAA4C,CAAC,CAAC;AACvE,CAAC;AAED,SAAS,KAAK,CAAC,CAAW,EAAE,CAAW;IACrC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC3F,CAAC"}

package/dist/ellipsoid-rhumb.d.ts

@@ -0,0 +1,52 @@
+import type { GeodesicParameters } from "./ellipsoid-geodesic.js";
+export type RhumbOptions = {
+    area?: boolean;
+    unrollLongitude?: boolean;
+    /** Extra whole turns added to the shortest inverse longitude difference. */
+    longitudeTurns?: number;
+};
+export type RhumbResult = {
+    lat1: number;
+    lon1: number;
+    lat2: number;
+    lon2: number;
+    azi12: number;
+    s12: number;
+    S12?: number;
+};
+export interface RhumbLine {
+    position(distance: number): RhumbResult;
+}
+/** A physical rhumb endpoint is undefined here, as distinct from invalid input or failed convergence. */
+export declare class RhumbDomainError extends RangeError {
+    readonly reason: 'pole' | 'pole-crossing';
+    constructor(reason: 'pole' | 'pole-crossing');
+}
+/** Constant-bearing paths on a revolution ellipsoid; explicit degrees and axis units.
+ * Original integral algorithm, qualified for 0.5 <= b/a <= 2. Exact poles are
+ * rejected rather than replaced by an arbitrary epsilon latitude.
+ */
+export declare class EllipsoidRhumb {
+    readonly surfaceArea: number;
+    private readonly a;
+    private readonly r;
+    private readonly b;
+    private readonly metrics;
+    constructor({ semiMajorAxis: a, semiMinorAxis: b }: GeodesicParameters);
+    inverse(lat1: number, lon1: number, lat2: number, lon2: number, options?: RhumbOptions): RhumbResult;
+    direct(latitude: number, longitude: number, azimuth: number, distance: number, options?: RhumbOptions): RhumbResult;
+    line(latitude: number, longitude: number, azimuth: number, options?: RhumbOptions): RhumbLine;
+    private means;
+    private psi;
+    private polarEndpoint;
+    private polarPsi;
+    private meridianMean;
+    private meridianDelta;
+    private inverseMeridian;
+    private inversePsi;
+    private potential;
+    private area;
+}
+/** @internal Shared with RhumbPolygon so even its first vertex rejects exact poles. */
+export declare function validateRhumbPoint(latitude: number, longitude: number): void;
+//# sourceMappingURL=ellipsoid-rhumb.d.ts.map

package/dist/ellipsoid-rhumb.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ellipsoid-rhumb.d.ts","sourceRoot":"","sources":["../src/ellipsoid-rhumb.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAC,kBAAkB,EAAC,gCAA6B;AAE7D,MAAM,MAAM,YAAY,GAAG;IACzB,IAAI,CAAC,EAAE,OAAO,CAAC;IACf,eAAe,CAAC,EAAE,OAAO,CAAC;IAC1B,4EAA4E;IAC5E,cAAc,CAAC,EAAE,MAAM,CAAC;CACzB,CAAC;AACF,MAAM,MAAM,WAAW,GAAG;IACxB,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,MAAM,CAAC;IACd,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,CAAC,EAAE,MAAM,CAAC;CACd,CAAC;AACF,MAAM,WAAW,SAAS;IACxB,QAAQ,CAAC,QAAQ,EAAE,MAAM,GAAG,WAAW,CAAC;CACzC;AAED,yGAAyG;AACzG,qBAAa,gBAAiB,SAAQ,UAAU;IAClC,QAAQ,CAAC,MAAM,EAAE,MAAM,GAAG,eAAe;gBAAhC,MAAM,EAAE,MAAM,GAAG,eAAe;CAQtD;AAGD;;;GAGG;AACH,qBAAa,cAAc;IACzB,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAS;IAC3B,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAS;IAC3B,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAS;IAC3B,OAAO,CAAC,QAAQ,CAAC,OAAO,CAAmB;gBAC/B,EAAC,aAAa,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAC,EAAE,kBAAkB;IAapE,OAAO,CACL,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,MAAM,EACZ,OAAO,GAAE,YAAiB,GACzB,WAAW;IA0Bd,MAAM,CACJ,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,EAAE,MAAM,EACf,QAAQ,EAAE,MAAM,EAChB,OAAO,GAAE,YAAiB,GACzB,WAAW;IAId,IAAI,CACF,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,EACjB,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,YAAiB,GACzB,SAAS;IA0DZ,OAAO,CAAC,KAAK;IAoBb,OAAO,CAAC,GAAG;IAIX,OAAO,CAAC,aAAa;IAyDrB,OAAO,CAAC,QAAQ;IAShB,OAAO,CAAC,YAAY;IAcpB,OAAO,CAAC,aAAa;IAIrB,OAAO,CAAC,eAAe;IA2BvB,OAAO,CAAC,UAAU;IAqBlB,OAAO,CAAC,SAAS;IAQjB,OAAO,CAAC,IAAI;CAwBb;AAED,uFAAuF;AACvF,wBAAgB,kBAAkB,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,IAAI,CAI5E"}