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

package/dist/projection/experimental/kernels/robin.js

@@ -0,0 +1,166 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileComment: Direct TypeScript port of proj4js 2.22.0, modified to reuse the input point and avoid per-coordinate result objects and Newton callbacks. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution. Table precision and interval selection follow PROJ 9.5.1 robin.cpp. See ../../../../projection/PROJ-LICENSE.txt for PROJ attribution.
+// Robinson projection
+// Based on https://github.com/OSGeo/proj.4/blob/master/src/PJ_robin.c
+// Polynomial coeficients from http://article.gmane.org/gmane.comp.gis.proj-4.devel/6039
+import { HALF_PI, D2R, R2D, EPSLN } from "../common/constants.js";
+import adjust_lon from "../common/adjust_lon.js";
+var COEFS_X = [
+    [1.0, 2.2199e-17, -7.15515e-5, 3.1103e-6],
+    [0.9986, -0.000482243, -2.4897e-5, -1.3309e-6],
+    [0.9954, -0.00083103, -4.48605e-5, -9.86701e-7],
+    [0.99, -0.00135364, -5.9661e-5, 3.6777e-6],
+    [0.9822, -0.00167442, -4.49547e-6, -5.72411e-6],
+    [0.973, -0.00214868, -9.03571e-5, 1.8736e-8],
+    [0.96, -0.00305085, -9.00761e-5, 1.64917e-6],
+    [0.9427, -0.00382792, -6.53386e-5, -2.6154e-6],
+    [0.9216, -0.00467746, -0.00010457, 4.81243e-6],
+    [0.8962, -0.00536223, -3.23831e-5, -5.43432e-6],
+    [0.8679, -0.00609363, -0.000113898, 3.32484e-6],
+    [0.835, -0.00698325, -6.40253e-5, 9.34959e-7],
+    [0.7986, -0.00755338, -5.00009e-5, 9.35324e-7],
+    [0.7597, -0.00798324, -3.5971e-5, -2.27626e-6],
+    [0.7186, -0.00851367, -7.01149e-5, -8.6303e-6],
+    [0.6732, -0.00986209, -0.000199569, 1.91974e-5],
+    [0.6213, -0.010418, 8.83923e-5, 6.24051e-6],
+    [0.5722, -0.00906601, 0.000182, 6.24051e-6],
+    [0.5322, -0.00677797, 0.000275608, 6.24051e-6]
+];
+var COEFS_Y = [
+    [-5.20417e-18, 0.0124, 1.21431e-18, -8.45284e-11],
+    [0.062, 0.0124, -1.26793e-9, 4.22642e-10],
+    [0.124, 0.0124, 5.07171e-9, -1.60604e-9],
+    [0.186, 0.0123999, -1.90189e-8, 6.00152e-9],
+    [0.248, 0.0124002, 7.10039e-8, -2.24e-8],
+    [0.31, 0.0123992, -2.64997e-7, 8.35986e-8],
+    [0.372, 0.0124029, 9.88983e-7, -3.11994e-7],
+    [0.434, 0.0123893, -3.69093e-6, -4.35621e-7],
+    [0.4958, 0.0123198, -1.02252e-5, -3.45523e-7],
+    [0.5571, 0.0121916, -1.54081e-5, -5.82288e-7],
+    [0.6176, 0.0119938, -2.41424e-5, -5.25327e-7],
+    [0.6769, 0.011713, -3.20223e-5, -5.16405e-7],
+    [0.7346, 0.0113541, -3.97684e-5, -6.09052e-7],
+    [0.7903, 0.0109107, -4.89042e-5, -1.04739e-6],
+    [0.8435, 0.0103431, -6.4615e-5, -1.40374e-9],
+    [0.8936, 0.00969686, -6.4636e-5, -8.547e-6],
+    [0.9394, 0.00840947, -0.000192841, -4.2106e-6],
+    [0.9761, 0.00616527, -0.000256, -4.2106e-6],
+    [1.0, 0.00328947, -0.000319159, -4.2106e-6]
+];
+// PROJ defines the published polynomial coefficients as float constants. Preserve
+// those values instead of silently changing the interpolation table to doubles.
+COEFS_X = COEFS_X.map(row => row.map(Math.fround));
+COEFS_Y = COEFS_Y.map(row => row.map(Math.fround));
+var FXC = 0.8487;
+var FYC = 1.3523;
+var C1 = R2D / 5; // rad to 5-degree interval
+var RC1 = 1 / C1;
+var NODES = 18;
+var poly3_val = function (coefs, x) {
+    return coefs[0] + x * (coefs[1] + x * (coefs[2] + x * coefs[3]));
+};
+var poly3_der = function (coefs, x) {
+    return coefs[1] + x * (2 * coefs[2] + x * 3 * coefs[3]);
+};
+function newton_rapshon(coefs, target, start, max_err, iters) {
+    var x = start;
+    for (; iters; --iters) {
+        var upd = (poly3_val(coefs, x) - target) / poly3_der(coefs, x);
+        x -= upd;
+        if (Math.abs(upd) < max_err) {
+            break;
+        }
+    }
+    return x;
+}
+function initialize(state) {
+    state.x0 = state.x0 || 0;
+    state.y0 = state.y0 || 0;
+    state.long0 = state.long0 || 0;
+    state.es = 0;
+}
+export function forward(state, ll) {
+    var lon = adjust_lon(ll.x - state.long0, state.over);
+    var dphi = Math.abs(ll.y);
+    var i = Math.floor(dphi * C1 + 1e-15);
+    if (i < 0) {
+        i = 0;
+    }
+    else if (i >= NODES) {
+        // The final table row represents the exact pole, as in PROJ.
+        i = NODES;
+    }
+    dphi = R2D * (dphi - RC1 * i);
+    var x = poly3_val(COEFS_X[i], dphi) * lon;
+    var y = poly3_val(COEFS_Y[i], dphi);
+    if (ll.y < 0) {
+        y = -y;
+    }
+    ll.x = x * state.a * FXC + state.x0;
+    ll.y = y * state.a * FYC + state.y0;
+    return ll;
+}
+export function inverse(state, xy) {
+    var originalY = xy.y;
+    var ll = xy;
+    ll.x = (ll.x - state.x0) / (state.a * FXC);
+    ll.y = Math.abs(originalY - state.y0) / (state.a * FYC);
+    // Polynomial pieces have small discontinuities at table knots. Roundoff in
+    // de-scaling must not select a neighbouring piece for an exact knot output.
+    for (let knot = 0; knot <= NODES; knot++) {
+        if (Math.abs(ll.y - COEFS_Y[knot][0]) <= 8 * Number.EPSILON) {
+            ll.x = adjust_lon(ll.x / COEFS_X[knot][0] + state.long0, state.over);
+            ll.y = (originalY < state.y0 ? -1 : 1) * knot * 5 * D2R;
+            return ll;
+        }
+    }
+    if (ll.y >= 1) {
+        // pathologic case
+        ll.x /= COEFS_X[NODES][0];
+        ll.y = originalY < 0 ? -HALF_PI : HALF_PI;
+    }
+    else {
+        // find table interval
+        var i = Math.floor(ll.y * NODES);
+        if (i < 0) {
+            i = 0;
+        }
+        else if (i >= NODES) {
+            i = NODES - 1;
+        }
+        for (;;) {
+            if (COEFS_Y[i][0] > ll.y) {
+                --i;
+            }
+            else if (COEFS_Y[i + 1][0] <= ll.y) {
+                ++i;
+            }
+            else {
+                break;
+            }
+        }
+        // linear interpolation in 5 degree interval
+        var coefs = COEFS_Y[i];
+        var t = (5 * (ll.y - coefs[0])) / (COEFS_Y[i + 1][0] - coefs[0]);
+        // find t so that poly3_val(coefs, t) = ll.y
+        t = newton_rapshon(coefs, ll.y, t, EPSLN, 100);
+        ll.x /= poly3_val(COEFS_X[i], t);
+        ll.y = (5 * i + t) * D2R;
+        if (originalY < 0) {
+            ll.y = -ll.y;
+        }
+    }
+    ll.x = adjust_lon(ll.x + state.long0, state.over);
+    return ll;
+}
+export function createState(base, options = {}) {
+    const state = {
+        ...base,
+        ...options
+    };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=robin.js.map

package/dist/projection/experimental/kernels/robin.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"robin.js","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/robin.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,0XAA0X;AAI1X,sBAAsB;AACtB,sEAAsE;AACtE,wFAAwF;AAExF,OAAO,EAAC,OAAO,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAC,+BAA4B;AAC7D,OAAO,UAAU,gCAA6B;AAE9C,IAAI,OAAO,GAAG;IACZ,CAAC,GAAG,EAAE,UAAU,EAAE,CAAC,UAAU,EAAE,SAAS,CAAC;IACzC,CAAC,MAAM,EAAE,CAAC,WAAW,EAAE,CAAC,SAAS,EAAE,CAAC,SAAS,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC/C,CAAC,IAAI,EAAE,CAAC,UAAU,EAAE,CAAC,SAAS,EAAE,SAAS,CAAC;IAC1C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC/C,CAAC,KAAK,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,SAAS,CAAC;IAC5C,CAAC,IAAI,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC5C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,SAAS,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC/C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,WAAW,EAAE,UAAU,CAAC;IAC/C,CAAC,KAAK,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,SAAS,EAAE,CAAC,UAAU,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC,SAAS,CAAC;IAC9C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,CAAC,WAAW,EAAE,UAAU,CAAC;IAC/C,CAAC,MAAM,EAAE,CAAC,QAAQ,EAAE,UAAU,EAAE,UAAU,CAAC;IAC3C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,QAAQ,EAAE,UAAU,CAAC;IAC3C,CAAC,MAAM,EAAE,CAAC,UAAU,EAAE,WAAW,EAAE,UAAU,CAAC;CAC/C,CAAC;AAEF,IAAI,OAAO,GAAG;IACZ,CAAC,CAAC,WAAW,EAAE,MAAM,EAAE,WAAW,EAAE,CAAC,WAAW,CAAC;IACjD,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,UAAU,EAAE,WAAW,CAAC;IACzC,CAAC,KAAK,EAAE,MAAM,EAAE,UAAU,EAAE,CAAC,UAAU,CAAC;IACxC,CAAC,KAAK,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC3C,CAAC,KAAK,EAAE,SAAS,EAAE,UAAU,EAAE,CAAC,OAAO,CAAC;IACxC,CAAC,IAAI,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,UAAU,CAAC;IAC1C,CAAC,KAAK,EAAE,SAAS,EAAE,UAAU,EAAE,CAAC,UAAU,CAAC;IAC3C,CAAC,KAAK,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC5C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,QAAQ,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC5C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,UAAU,EAAE,CAAC,UAAU,CAAC;IAC7C,CAAC,MAAM,EAAE,SAAS,EAAE,CAAC,SAAS,EAAE,CAAC,UAAU,CAAC;IAC5C,CAAC,MAAM,EAAE,UAAU,EAAE,CAAC,SAAS,EAAE,CAAC,QAAQ,CAAC;IAC3C,CAAC,MAAM,EAAE,UAAU,EAAE,CAAC,WAAW,EAAE,CAAC,SAAS,CAAC;IAC9C,CAAC,MAAM,EAAE,UAAU,EAAE,CAAC,QAAQ,EAAE,CAAC,SAAS,CAAC;IAC3C,CAAC,GAAG,EAAE,UAAU,EAAE,CAAC,WAAW,EAAE,CAAC,SAAS,CAAC;CAC5C,CAAC;AAEF,kFAAkF;AAClF,gFAAgF;AAChF,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;AACnD,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC;AAEnD,IAAI,GAAG,GAAG,MAAM,CAAC;AACjB,IAAI,GAAG,GAAG,MAAM,CAAC;AACjB,IAAI,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,2BAA2B;AAC7C,IAAI,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;AACjB,IAAI,KAAK,GAAG,EAAE,CAAC;AAEf,IAAI,SAAS,GAAG,UAAU,KAAwB,EAAE,CAAS;IAC3D,OAAO,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC,CAAC;AAEF,IAAI,SAAS,GAAG,UAAU,KAAwB,EAAE,CAAS;IAC3D,OAAO,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1D,CAAC,CAAC;AAEF,SAAS,cAAc,CACrB,KAAwB,EACxB,MAAc,EACd,KAAa,EACb,OAAe,EACf,KAAa;IAEb,IAAI,CAAC,GAAG,KAAK,CAAC;IACd,OAAO,KAAK,EAAE,EAAE,KAAK,EAAE,CAAC;QACtB,IAAI,GAAG,GAAG,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;QAC/D,CAAC,IAAI,GAAG,CAAC;QACT,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,OAAO,EAAE,CAAC;YAC5B,MAAM;QACR,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,UAAU,CAAC,KAAY;IAC9B,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC;IACzB,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC;IACzB,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,IAAI,CAAC,CAAC;IAC/B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;AACf,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,EAAS;IAC7C,IAAI,GAAG,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IAErD,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC1B,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;IACtC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QACV,CAAC,GAAG,CAAC,CAAC;IACR,CAAC;SAAM,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;QACtB,6DAA6D;QAC7D,CAAC,GAAG,KAAK,CAAC;IACZ,CAAC;IACD,IAAI,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;IAC9B,IAAI,CAAC,GAAG,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;IAC1C,IAAI,CAAC,GAAG,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACpC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,CAAC,GAAG,CAAC,CAAC,CAAC;IACT,CAAC;IAED,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC;IACpC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC;IAEpC,OAAO,EAAE,CAAC;AACZ,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,EAAS;IAC7C,IAAI,SAAS,GAAG,EAAE,CAAC,CAAC,CAAC;IACrB,IAAI,EAAE,GAAG,EAAE,CAAC;IACZ,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAC3C,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAExD,2EAA2E;IAC3E,4EAA4E;IAC5E,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,IAAI,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC;QACzC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC;YAC5D,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;YACrE,EAAE,CAAC,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC;YACxD,OAAO,EAAE,CAAC;QACZ,CAAC;IACH,CAAC;IACD,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC;QACd,kBAAkB;QAClB,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1B,EAAE,CAAC,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC;IAC5C,CAAC;SAAM,CAAC;QACN,sBAAsB;QACtB,IAAI,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACjC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,CAAC,GAAG,CAAC,CAAC;QACR,CAAC;aAAM,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;YACtB,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;QAChB,CAAC;QACD,SAAS,CAAC;YACR,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC;gBACzB,EAAE,CAAC,CAAC;YACN,CAAC;iBAAM,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC;gBACrC,EAAE,CAAC,CAAC;YACN,CAAC;iBAAM,CAAC;gBACN,MAAM;YACR,CAAC;QACH,CAAC;QACD,4CAA4C;QAC5C,IAAI,KAAK,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QACjE,4CAA4C;QAC5C,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAE/C,EAAE,CAAC,CAAC,IAAI,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACjC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACzB,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;YAClB,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QACf,CAAC;IACH,CAAC;IAED,EAAE,CAAC,CAAC,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IAElD,OAAO,EAAE,CAAC;AACZ,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU;QACnB,GAAG,IAAI;QAEP,GAAG,OAAO;KACX,CAAC;IACF,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/projection/experimental/kernels/sinu.d.ts

@@ -0,0 +1,12 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    en: number[];
+    n: number;
+    m: number;
+    C_y: number;
+    C_x: number;
+};
+export declare function forward(state: State, p: Point): Point | null | undefined | number;
+export declare function inverse(state: State, p: Point): Point | null | undefined | number;
+export declare function createState(base: KernelParameters, options?: Partial<State>): State;
+//# sourceMappingURL=sinu.d.ts.map

package/dist/projection/experimental/kernels/sinu.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sinu.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/sinu.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,EAAE,EAAE,MAAM,EAAE,CAAC;IACb,CAAC,EAAE,MAAM,CAAC;IACV,CAAC,EAAE,MAAM,CAAC;IACV,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;CACb,CAAC;AA6BF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAiCjF;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAmCjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAIvF"}

package/dist/projection/experimental/kernels/sinu.js

@@ -0,0 +1,107 @@
+// 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 with bounded Eckert VI inverse pole-line validation. See ../../../../projection/PROJ4-LICENSE.md for upstream attribution. Independent corrections acknowledge PROJ 9.5.1 src/projections/gn_sinu.cpp; see ../../../../projection/PROJ-LICENSE.txt.
+import adjust_lon from "../common/adjust_lon.js";
+import adjust_lat from "../common/adjust_lat.js";
+import pj_enfn from "../common/pj_enfn.js";
+var MAX_ITER = 20;
+import pj_mlfn from "../common/pj_mlfn.js";
+import pj_inv_mlfn from "../common/pj_inv_mlfn.js";
+import { EPSLN, HALF_PI } from "../common/constants.js";
+import asinz from "../common/asinz.js";
+function initialize(state) {
+    /* Place parameters in static storage for common use
+      ------------------------------------------------- */
+    state.long0 = state.long0 || 0;
+    if (!state.sphere) {
+        state.en = pj_enfn(state.es);
+    }
+    else {
+        state.n = 1;
+        state.m = 0;
+        state.es = 0;
+        state.C_y = Math.sqrt((state.m + 1) / state.n);
+        state.C_x = state.C_y / (state.m + 1);
+    }
+}
+/* Sinusoidal forward equations--mapping lat,long to x,y
+  ----------------------------------------------------- */
+export function forward(state, p) {
+    var x, y;
+    var lon = p.x;
+    var lat = p.y;
+    /* Forward equations
+      ----------------- */
+    lon = adjust_lon(lon - state.long0, state.over);
+    if (state.sphere) {
+        if (!state.m) {
+            lat = state.n !== 1 ? Math.asin(state.n * Math.sin(lat)) : lat;
+        }
+        else {
+            var k = state.n * Math.sin(lat);
+            for (var i = MAX_ITER; i; --i) {
+                var V = (state.m * lat + Math.sin(lat) - k) / (state.m + Math.cos(lat));
+                lat -= V;
+                if (Math.abs(V) < EPSLN) {
+                    break;
+                }
+            }
+        }
+        x = state.a * state.C_x * lon * (state.m + Math.cos(lat));
+        y = state.a * state.C_y * lat;
+    }
+    else {
+        var s = Math.sin(lat);
+        var c = Math.cos(lat);
+        y = state.a * pj_mlfn(lat, s, c, state.en);
+        x = (state.a * lon * c) / Math.sqrt(1 - state.es * s * s);
+    }
+    p.x = x;
+    p.y = y;
+    return p;
+}
+export function inverse(state, p) {
+    var lat, temp, lon, s;
+    p.x -= state.x0;
+    lon = p.x / state.a;
+    p.y -= state.y0;
+    lat = p.y / state.a;
+    if (state.sphere) {
+        lat /= state.C_y;
+        lon = lon / (state.C_x * (state.m + Math.cos(lat)));
+        if (state.m) {
+            const sinLatitude = (state.m * lat + Math.sin(lat)) / state.n;
+            if (Math.abs(sinLatitude) > 1 + 8 * Number.EPSILON)
+                return null;
+            lat = asinz(sinLatitude);
+        }
+        else if (state.n !== 1) {
+            lat = asinz(Math.sin(lat) / state.n);
+        }
+        lon = adjust_lon(lon + state.long0, state.over);
+        lat = adjust_lat(lat);
+    }
+    else {
+        lat = pj_inv_mlfn(p.y / state.a, state.es, state.en);
+        s = Math.abs(lat);
+        if (s < HALF_PI) {
+            s = Math.sin(lat);
+            temp = state.long0 + (p.x * Math.sqrt(1 - state.es * s * s)) / (state.a * Math.cos(lat));
+            // temp = state.long0 + p.x / (state.a * Math.cos(lat));
+            lon = adjust_lon(temp, state.over);
+        }
+        else if (s - EPSLN < HALF_PI) {
+            lon = state.long0;
+        }
+    }
+    p.x = lon;
+    p.y = lat;
+    return p;
+}
+export function createState(base, options = {}) {
+    const state = { ...base, en: [], n: 0, m: 0, C_y: 0, C_x: 0, ...options };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=sinu.js.map

package/dist/projection/experimental/kernels/sinu.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sinu.js","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/sinu.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,oSAAoS;AAUpS,OAAO,UAAU,gCAA6B;AAC9C,OAAO,UAAU,gCAA6B;AAC9C,OAAO,OAAO,6BAA0B;AACxC,IAAI,QAAQ,GAAG,EAAE,CAAC;AAClB,OAAO,OAAO,6BAA0B;AACxC,OAAO,WAAW,iCAA8B;AAChD,OAAO,EAAC,KAAK,EAAE,OAAO,EAAC,+BAA4B;AAEnD,OAAO,KAAK,2BAAwB;AAEpC,SAAS,UAAU,CAAC,KAAY;IAC9B;0DACsD;IACtD,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,IAAI,CAAC,CAAC;IAE/B,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QAClB,KAAK,CAAC,EAAE,GAAG,OAAO,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;IAC/B,CAAC;SAAM,CAAC;QACN,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QACZ,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QACZ,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;QACb,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QAC/C,KAAK,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxC,CAAC;AACH,CAAC;AAED;0DAC0D;AAC1D,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,CAAC,EAAE,CAAC,CAAC;IACT,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;IACd;0BACsB;IACtB,GAAG,GAAG,UAAU,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;IAEhD,IAAI,KAAK,CAAC,MAAM,EAAE,CAAC;QACjB,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;YACb,GAAG,GAAG,KAAK,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACjE,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChC,KAAK,IAAI,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC;gBAC9B,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBACxE,GAAG,IAAI,CAAC,CAAC;gBACT,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE,CAAC;oBACxB,MAAM;gBACR,CAAC;YACH,CAAC;QACH,CAAC;QACD,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1D,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,GAAG,GAAG,CAAC;IAChC,CAAC;SAAM,CAAC;QACN,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtB,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtB,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,OAAO,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC,CAAC;QAC3C,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACR,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACR,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,IAAI,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC;IAEtB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;IAChB,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;IAChB,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAEpB,IAAI,KAAK,CAAC,MAAM,EAAE,CAAC;QACjB,GAAG,IAAI,KAAK,CAAC,GAAG,CAAC;QACjB,GAAG,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpD,IAAI,KAAK,CAAC,CAAC,EAAE,CAAC;YACZ,MAAM,WAAW,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;YAC9D,IAAI,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,OAAO;gBAAE,OAAO,IAAI,CAAC;YAChE,GAAG,GAAG,KAAK,CAAC,WAAW,CAAC,CAAC;QAC3B,CAAC;aAAM,IAAI,KAAK,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;YACzB,GAAG,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACvC,CAAC;QACD,GAAG,GAAG,UAAU,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QAChD,GAAG,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;IACxB,CAAC;SAAM,CAAC;QACN,GAAG,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,EAAE,KAAK,CAAC,EAAE,CAAC,CAAC;QACrD,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAClB,IAAI,CAAC,GAAG,OAAO,EAAE,CAAC;YAChB,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAClB,IAAI,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YACzF,wDAAwD;YACxD,GAAG,GAAG,UAAU,CAAC,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC;QACrC,CAAC;aAAM,IAAI,CAAC,GAAG,KAAK,GAAG,OAAO,EAAE,CAAC;YAC/B,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC;QACpB,CAAC;IACH,CAAC;IACD,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACV,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU,EAAC,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,OAAO,EAAC,CAAC;IAC/E,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/projection/experimental/kernels/somerc.d.ts

@@ -0,0 +1,13 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    lambda0: number;
+    rf: number;
+    R: number;
+    alpha: number;
+    b0: number;
+    K: number;
+};
+export declare function forward(state: State, p: Point): Point | null;
+export declare function inverse(state: State, p: Point): Point | null;
+export declare function createState(base: KernelParameters, options?: Partial<State>): State;
+//# sourceMappingURL=somerc.d.ts.map

package/dist/projection/experimental/kernels/somerc.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"somerc.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/somerc.ts"],"names":[],"mappings":"AAMA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AACvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,OAAO,EAAE,MAAM,CAAC;IAChB,EAAE,EAAE,MAAM,CAAC;IACX,CAAC,EAAE,MAAM,CAAC;IACV,KAAK,EAAE,MAAM,CAAC;IACd,EAAE,EAAE,MAAM,CAAC;IACX,CAAC,EAAE,MAAM,CAAC;CACX,CAAC;AA+BF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,CAsB5D;AAED,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,CA0B5D;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE,OAAO,CAAC,KAAK,CAAM,GAAG,KAAK,CAIvF"}

package/dist/projection/experimental/kernels/somerc.js

@@ -0,0 +1,79 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileCopyrightText: Copyright (c) PROJ contributors
+// SPDX-FileComment: Initialization adapted from proj4js 2.22.0; forward branch and inverse Newton equations adapted from PROJ 9.5.1 src/projections/somerc.cpp. See ../../../../projection/PROJ-LICENSE.txt. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
+/*
+  references:
+    Formules et constantes pour le Calcul pour la
+    projection cylindrique conforme à axe oblique et pour la transformation entre
+    des systèmes de référence.
+    http://www.swisstopo.admin.ch/internet/swisstopo/fr/home/topics/survey/sys/refsys/switzerland.parsysrelated1.31216.downloadList.77004.DownloadFile.tmp/swissprojectionfr.pdf
+  */
+function initialize(state) {
+    if (!state.k0) {
+        state.k0 = 1;
+    }
+    var phy0 = state.lat0;
+    state.lambda0 = state.long0;
+    var sinPhy0 = Math.sin(phy0);
+    var semiMajorAxis = state.a;
+    var invF = state.rf;
+    var flattening = 1 / invF;
+    var e2 = 2 * flattening - Math.pow(flattening, 2);
+    state.e = Math.sqrt(e2);
+    var e = state.e;
+    state.R = (state.k0 * semiMajorAxis * Math.sqrt(1 - e2)) / (1 - e2 * Math.pow(sinPhy0, 2));
+    state.alpha = Math.sqrt(1 + (e2 / (1 - e2)) * Math.pow(Math.cos(phy0), 4));
+    state.b0 = Math.asin(sinPhy0 / state.alpha);
+    var k1 = Math.log(Math.tan(Math.PI / 4 + state.b0 / 2));
+    var k2 = Math.log(Math.tan(Math.PI / 4 + phy0 / 2));
+    var k3 = Math.log((1 + e * sinPhy0) / (1 - e * sinPhy0));
+    state.K = k1 - state.alpha * k2 + ((state.alpha * e) / 2) * k3;
+}
+export function forward(state, p) {
+    const sp = state.e * Math.sin(p.y);
+    const latitude = 2 *
+        Math.atan(Math.exp(state.alpha *
+            (Math.log(Math.tan(Math.PI / 4 + p.y / 2)) -
+                (state.e / 2) * Math.log((1 + sp) / (1 - sp))) +
+            state.K)) -
+        Math.PI / 2;
+    const longitude = state.alpha * (p.x - state.lambda0);
+    const cos = Math.cos(latitude);
+    const rotatedLatitude = Math.asin(Math.cos(state.b0) * Math.sin(latitude) - Math.sin(state.b0) * cos * Math.cos(longitude));
+    const rotatedLongitude = Math.asin((cos * Math.sin(longitude)) / Math.cos(rotatedLatitude));
+    p.x = state.x0 + state.R * rotatedLongitude;
+    p.y = state.y0 + state.R * Math.log(Math.tan(Math.PI / 4 + rotatedLatitude / 2));
+    return p;
+}
+export function inverse(state, p) {
+    const longitude = (p.x - state.x0) / state.R;
+    const latitude = 2 * (Math.atan(Math.exp((p.y - state.y0) / state.R)) - Math.PI / 4);
+    let phi = Math.asin(Math.cos(state.b0) * Math.sin(latitude) +
+        Math.sin(state.b0) * Math.cos(latitude) * Math.cos(longitude));
+    const lambda = Math.asin((Math.cos(latitude) * Math.sin(longitude)) / Math.cos(phi));
+    const con = (state.K - Math.log(Math.tan(Math.PI / 4 + phi / 2))) / state.alpha;
+    for (let iteration = 0; iteration < 6; iteration++) {
+        const esp = state.e * Math.sin(phi);
+        const step = ((con +
+            Math.log(Math.tan(Math.PI / 4 + phi / 2)) -
+            (state.e / 2) * Math.log((1 + esp) / (1 - esp))) *
+            (1 - esp * esp) *
+            Math.cos(phi)) /
+            (1 - state.es);
+        phi -= step;
+        if (Math.abs(step) < 1e-10) {
+            p.x = state.lambda0 + lambda / state.alpha;
+            p.y = phi;
+            return p;
+        }
+    }
+    return null;
+}
+export function createState(base, options = {}) {
+    const state = { ...base, lambda0: 0, rf: 0, R: 0, alpha: 0, b0: 0, K: 0, ...options };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=somerc.js.map

package/dist/projection/experimental/kernels/somerc.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"somerc.js","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/somerc.ts"],"names":[],"mappings":"AAAA,UAAU;AACV,+BAA+B;AAC/B,yIAAyI;AACzI,0DAA0D;AAC1D,mSAAmS;AAWnS;;;;;;IAMI;AAEJ,SAAS,UAAU,CAAC,KAAY;IAC9B,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC;QACd,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;IACf,CAAC;IACD,IAAI,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC;IACtB,KAAK,CAAC,OAAO,GAAG,KAAK,CAAC,KAAK,CAAC;IAC5B,IAAI,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAC7B,IAAI,aAAa,GAAG,KAAK,CAAC,CAAC,CAAC;IAC5B,IAAI,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;IACpB,IAAI,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC;IAC1B,IAAI,EAAE,GAAG,CAAC,GAAG,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC;IAClD,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxB,IAAI,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAChB,KAAK,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,aAAa,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC;IAC3F,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC3E,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC;IAC5C,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;IACxD,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IACpD,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC;IACzD,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,KAAK,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;AACjE,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,QAAQ,GACZ,CAAC;QACC,IAAI,CAAC,IAAI,CACP,IAAI,CAAC,GAAG,CACN,KAAK,CAAC,KAAK;YACT,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACxC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YAChD,KAAK,CAAC,CAAC,CACV,CACF;QACH,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACd,MAAM,SAAS,GAAG,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC;IACtD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC;IAC/B,MAAM,eAAe,GAAG,IAAI,CAAC,IAAI,CAC/B,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CACzF,CAAC;IACF,MAAM,gBAAgB,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC,CAAC;IAC5F,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,gBAAgB,CAAC;IAC5C,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,eAAe,GAAG,CAAC,CAAC,CAAC,CAAC;IACjF,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,KAAY,EAAE,CAAQ;IAC5C,MAAM,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC7C,MAAM,QAAQ,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACrF,IAAI,GAAG,GAAG,IAAI,CAAC,IAAI,CACjB,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;QACrC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAChE,CAAC;IACF,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACrF,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC;IAChF,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,EAAE,SAAS,EAAE,EAAE,CAAC;QACnD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACpC,MAAM,IAAI,GACR,CAAC,CAAC,GAAG;YACH,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;YACzC,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;YAChD,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YACf,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;QACjB,GAAG,IAAI,IAAI,CAAC;QACZ,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC;YAC3B,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,OAAO,GAAG,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC;YAC3C,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACV,OAAO,CAAC,CAAC;QACX,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,IAAsB,EAAE,UAA0B,EAAE;IAC9E,MAAM,KAAK,GAAU,EAAC,GAAG,IAAI,EAAE,OAAO,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,OAAO,EAAC,CAAC;IAC3F,UAAU,CAAC,KAAK,CAAC,CAAC;IAClB,OAAO,KAAK,CAAC;AACf,CAAC"}

package/dist/projection/experimental/kernels/stere.d.ts

@@ -0,0 +1,15 @@
+import type { KernelParameters, Point } from "../kernel.js";
+export type State = KernelParameters & {
+    coslat0: number;
+    sinlat0: number;
+    con: number;
+    cons: number;
+    ms1: number;
+    X0: number;
+    cosX0: number;
+    sinX0: number;
+};
+export declare function forward(state: State, p: Point): Point | null | undefined | number;
+export declare function inverse(state: State, p: Point): Point | null | undefined | number;
+export declare function createState(base: KernelParameters): State;
+//# sourceMappingURL=stere.d.ts.map

package/dist/projection/experimental/kernels/stere.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"stere.d.ts","sourceRoot":"","sources":["../../../../src/projection/experimental/kernels/stere.ts"],"names":[],"mappings":"AAOA,OAAO,KAAK,EAAC,gBAAgB,EAAE,KAAK,EAAC,qBAAkB;AAEvD,MAAM,MAAM,KAAK,GAAG,gBAAgB,GAAG;IACrC,OAAO,EAAE,MAAM,CAAC;IAChB,OAAO,EAAE,MAAM,CAAC;IAChB,GAAG,EAAE,MAAM,CAAC;IACZ,IAAI,EAAE,MAAM,CAAC;IACb,GAAG,EAAE,MAAM,CAAC;IACZ,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;CACf,CAAC;AA2DF,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CA+DjF;AAGD,wBAAgB,OAAO,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,GAAG,KAAK,GAAG,IAAI,GAAG,SAAS,GAAG,MAAM,CAyFjF;AAED,wBAAgB,WAAW,CAAC,IAAI,EAAE,gBAAgB,GAAG,KAAK,CAczD"}

package/dist/projection/experimental/kernels/stere.js

@@ -0,0 +1,207 @@
+// math.gl
+// SPDX-License-Identifier: MIT
+// SPDX-FileCopyrightText: 2014 Mike Adair, Richard Greenwood, Didier Richard, Stephen Irons, Olivier Terral and Calvin Metcalf (proj4js)
+// SPDX-FileCopyrightText: Copyright (c) PROJ contributors
+// SPDX-FileCopyrightText: Copyright (c) vis.gl contributors
+// SPDX-FileComment: Original math.gl inverse Cartesian rotation correction. Adapted from proj4js 2.22.0, with polar half-angle equations from PROJ 9.5.1 src/projections/stere.cpp. See ../../../../projection/PROJ-LICENSE.txt. See ../../../../projection/PROJ4-LICENSE.md for the upstream license and attribution.
+import { inverseSphericalLatitude } from "../inverse-spherical-latitude.js";
+import { EPSLN, HALF_PI } from "../common/constants.js";
+import sign from "../common/sign.js";
+import msfnz from "../common/msfnz.js";
+import tsfnz from "../common/tsfnz.js";
+import phi2z from "../common/phi2z.js";
+import adjust_lon from "../common/adjust_lon.js";
+function ssfn_(phit, sinphi, eccen) {
+    sinphi *= eccen;
+    return Math.tan(0.5 * (HALF_PI + phit)) * Math.pow((1 - sinphi) / (1 + sinphi), 0.5 * eccen);
+}
+function initialize(state) {
+    // setting default parameters
+    state.x0 = state.x0 || 0;
+    state.y0 = state.y0 || 0;
+    state.lat0 = state.lat0 || 0;
+    state.long0 = state.long0 || 0;
+    state.coslat0 = Math.cos(state.lat0);
+    state.sinlat0 = Math.sin(state.lat0);
+    if (state.sphere) {
+        if (!Number.isNaN(state.lat_ts) && Math.abs(state.coslat0) <= EPSLN) {
+            state.k0 = 0.5 * (1 + sign(state.lat0) * Math.sin(state.lat_ts));
+        }
+    }
+    else {
+        if (Math.abs(state.coslat0) <= EPSLN) {
+            if (state.lat0 > 0) {
+                // North pole
+                // trace('stere:north pole');
+                state.con = 1;
+            }
+            else {
+                // South pole
+                // trace('stere:south pole');
+                state.con = -1;
+            }
+        }
+        state.cons = Math.sqrt(Math.pow(1 + state.e, 1 + state.e) * Math.pow(1 - state.e, 1 - state.e));
+        if (!Number.isNaN(state.lat_ts) &&
+            Math.abs(state.coslat0) <= EPSLN &&
+            Math.abs(Math.cos(state.lat_ts)) > EPSLN) {
+            // Polar Stereographic variant B: k0 is derived from lat_ts (EPSG guidance note 7-2, method 9829)
+            // https://epsg.org/coord-operation-method_9829/Polar-Stereographic-variant-B.html
+            state.k0 =
+                (0.5 * state.cons * msfnz(state.e, Math.sin(state.lat_ts), Math.cos(state.lat_ts))) /
+                    tsfnz(state.e, state.con * state.lat_ts, state.con * Math.sin(state.lat_ts));
+        }
+        state.ms1 = msfnz(state.e, state.sinlat0, state.coslat0);
+        state.X0 = 2 * Math.atan(ssfn_(state.lat0, state.sinlat0, state.e)) - HALF_PI;
+        state.cosX0 = Math.cos(state.X0);
+        state.sinX0 = Math.sin(state.X0);
+    }
+}
+// Stereographic forward equations--mapping lat,long to x,y
+export function forward(state, p) {
+    var lon = p.x;
+    var lat = p.y;
+    var sinlat = Math.sin(lat);
+    var coslat = Math.cos(lat);
+    var A, X, sinX, cosX, ts, rh;
+    var dlon = adjust_lon(lon - state.long0, state.over);
+    if (Math.abs(Math.abs(lon - state.long0) - Math.PI) <= EPSLN &&
+        Math.abs(lat + state.lat0) <= EPSLN) {
+        // case of the origine point
+        // trace('stere:this is the origin point');
+        p.x = NaN;
+        p.y = NaN;
+        return p;
+    }
+    if (state.sphere && Math.abs(state.coslat0) <= EPSLN) {
+        // PROJ's polar half-angle form avoids subtracting nearly equal values in
+        // 1 +/- sin(latitude) close to the opposite pole.
+        const sign = state.lat0 > 0 ? 1 : -1;
+        const radius = 2 * state.a * state.k0 * Math.tan(Math.PI / 4 - (sign * lat) / 2);
+        p.x = state.x0 + radius * Math.sin(dlon);
+        p.y = state.y0 - sign * radius * Math.cos(dlon);
+        return p;
+    }
+    if (state.sphere) {
+        // trace('stere:sphere case');
+        A = (2 * state.k0) / (1 + state.sinlat0 * sinlat + state.coslat0 * coslat * Math.cos(dlon));
+        p.x = state.a * A * coslat * Math.sin(dlon) + state.x0;
+        p.y =
+            state.a * A * (state.coslat0 * sinlat - state.sinlat0 * coslat * Math.cos(dlon)) + state.y0;
+        return p;
+    }
+    else {
+        X = 2 * Math.atan(ssfn_(lat, sinlat, state.e)) - HALF_PI;
+        cosX = Math.cos(X);
+        sinX = Math.sin(X);
+        if (Math.abs(state.coslat0) <= EPSLN) {
+            ts = tsfnz(state.e, lat * state.con, state.con * sinlat);
+            rh = (2 * state.a * state.k0 * ts) / state.cons;
+            p.x = state.x0 + rh * Math.sin(lon - state.long0);
+            p.y = state.y0 - state.con * rh * Math.cos(lon - state.long0);
+            // trace(p.toString());
+            return p;
+        }
+        else if (Math.abs(state.sinlat0) < EPSLN) {
+            // Eq
+            // trace('stere:equateur');
+            A = (2 * state.a * state.k0) / (1 + cosX * Math.cos(dlon));
+            // Correct the missing false northing in proj4js 2.22.0 equatorial stereographic.
+            p.y = A * sinX + state.y0;
+        }
+        else {
+            // other case
+            // trace('stere:normal case');
+            A =
+                (2 * state.a * state.k0 * state.ms1) /
+                    (state.cosX0 * (1 + state.sinX0 * sinX + state.cosX0 * cosX * Math.cos(dlon)));
+            p.y = A * (state.cosX0 * sinX - state.sinX0 * cosX * Math.cos(dlon)) + state.y0;
+        }
+        p.x = A * cosX * Math.sin(dlon) + state.x0;
+    }
+    // trace(p.toString());
+    return p;
+}
+//* Stereographic inverse equations--mapping x,y to lat/long
+export function inverse(state, p) {
+    p.x -= state.x0;
+    p.y -= state.y0;
+    var lon, lat, ts, ce, Chi;
+    var rh = Math.sqrt(p.x * p.x + p.y * p.y);
+    if (state.sphere) {
+        var c = 2 * Math.atan(rh / (2 * state.a * state.k0));
+        lon = state.long0;
+        lat = state.lat0;
+        if (rh <= EPSLN) {
+            p.x = lon;
+            p.y = lat;
+            return p;
+        }
+        lat = inverseSphericalLatitude(state.sinlat0, state.coslat0, Math.sin(c), Math.cos(c), p.x, p.y, rh);
+        if (Math.abs(state.coslat0) < EPSLN) {
+            if (state.lat0 > 0) {
+                lon = adjust_lon(state.long0 + Math.atan2(p.x, -1 * p.y), state.over);
+            }
+            else {
+                lon = adjust_lon(state.long0 + Math.atan2(p.x, p.y), state.over);
+            }
+        }
+        else {
+            lon = adjust_lon(state.long0 +
+                Math.atan2(p.x * Math.sin(c), rh * state.coslat0 * Math.cos(c) - p.y * state.sinlat0 * Math.sin(c)), state.over);
+        }
+        p.x = lon;
+        p.y = lat;
+        return p;
+    }
+    else {
+        if (Math.abs(state.coslat0) <= EPSLN) {
+            if (rh <= EPSLN) {
+                lat = state.lat0;
+                lon = state.long0;
+                p.x = lon;
+                p.y = lat;
+                // trace(p.toString());
+                return p;
+            }
+            p.x *= state.con;
+            p.y *= state.con;
+            ts = (rh * state.cons) / (2 * state.a * state.k0);
+            lat = state.con * phi2z(state.e, ts);
+            lon = state.con * adjust_lon(state.con * state.long0 + Math.atan2(p.x, -1 * p.y), state.over);
+        }
+        else {
+            ce = 2 * Math.atan((rh * state.cosX0) / (2 * state.a * state.k0 * state.ms1));
+            lon = state.long0;
+            if (rh <= EPSLN) {
+                Chi = state.X0;
+            }
+            else {
+                Chi = inverseSphericalLatitude(state.sinX0, state.cosX0, Math.sin(ce), Math.cos(ce), p.x, p.y, rh);
+                lon = adjust_lon(state.long0 +
+                    Math.atan2(p.x * Math.sin(ce), rh * state.cosX0 * Math.cos(ce) - p.y * state.sinX0 * Math.sin(ce)), state.over);
+            }
+            lat = -1 * phi2z(state.e, Math.tan(0.5 * (HALF_PI + Chi)));
+        }
+    }
+    p.x = lon;
+    p.y = lat;
+    // trace(p.toString());
+    return p;
+}
+export function createState(base) {
+    const state = {
+        ...base,
+        coslat0: 0,
+        sinlat0: 0,
+        con: 0,
+        cons: 0,
+        ms1: 0,
+        X0: 0,
+        cosX0: 0,
+        sinX0: 0
+    };
+    initialize(state);
+    return state;
+}
+//# sourceMappingURL=stere.js.map