@itwin/core-orbitgt 3.0.0-dev.77 → 3.0.0-dev.81

package/lib/esm/spatial/ecrs/projection/ObliqueMercator.d.ts

@@ -0,0 +1,97 @@
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+declare type int32 = number;
+import { Coordinate } from "../../geom/Coordinate";
+import { CRS } from "../CRS";
+import { Operation } from "../Operation";
+import { OperationMethod } from "../OperationMethod";
+import { ParameterValueList } from "../ParameterValueList";
+/**
+ * Class ObliqueMercator defines an Oblique Mercator projection (Hotine Variant B).
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.6 (spec v6.13): Oblique Mercator and Hotine Oblique Mercator
+ *
+ * @version 1.0 May 2008
+ */
+/** @internal */
+export declare class ObliqueMercator extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The value of PI */
+    private static readonly PI;
+    /** Latitude of the projection center */
+    private _latC;
+    /** Longitude of the projection center */
+    private _lonC;
+    /** Azimuth of the initial line passing through the projection center */
+    private _aziC;
+    /** Angle from the recified grid to the skew (oblique) grid */
+    private _gamC;
+    /** Scale factor on the initial line of projection */
+    private _kC;
+    /** False easting at the center of projection */
+    private _eC;
+    /** False northing at the center of projection */
+    private _nC;
+    /** Forward calculation parameters */
+    private _a;
+    private _e;
+    private _e2;
+    private _B;
+    private _A;
+    private _tO;
+    private _D;
+    private _F;
+    private _H;
+    private _G;
+    private _gamO;
+    private _lonO;
+    private _vC;
+    private _uC;
+    /** Reverse calculation parameters */
+    private _rev1;
+    private _rev2;
+    private _rev3;
+    private _rev4;
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters: ParameterValueList);
+    /**
+       * Get the sign of a number.
+       */
+    private static sign;
+    /**
+       * Get the power of a number (must be able to handle negative 'n' values).
+       */
+    private static pow;
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation: Operation): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+}
+export {};
+//# sourceMappingURL=ObliqueMercator.d.ts.map

package/lib/esm/spatial/ecrs/projection/ObliqueMercator.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObliqueMercator.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/ObliqueMercator.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,aAAK,KAAK,GAAG,MAAM,CAAC;AAIpB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,eAAgB,SAAQ,eAAe;IAChD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,sBAAsB;IACtB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAoB;IAE9C,wCAAwC;IACxC,OAAO,CAAC,KAAK,CAAU;IACvB,yCAAyC;IACzC,OAAO,CAAC,KAAK,CAAU;IACvB,wEAAwE;IACxE,OAAO,CAAC,KAAK,CAAU;IACvB,8DAA8D;IAC9D,OAAO,CAAC,KAAK,CAAU;IACvB,qDAAqD;IACrD,OAAO,CAAC,GAAG,CAAU;IACrB,gDAAgD;IAChD,OAAO,CAAC,GAAG,CAAU;IACrB,iDAAiD;IACjD,OAAO,CAAC,GAAG,CAAU;IAErB,qCAAqC;IACrC,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,qCAAqC;IACrC,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,KAAK,CAAU;IAEvB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAYjD;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,IAAI;IAInB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAMlB;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAoCtD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IA+B5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAwB/F"}

package/lib/esm/spatial/ecrs/projection/ObliqueMercator.js

@@ -0,0 +1,170 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+import { OperationMethod } from "../OperationMethod";
+/**
+ * Class ObliqueMercator defines an Oblique Mercator projection (Hotine Variant B).
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.6 (spec v6.13): Oblique Mercator and Hotine Oblique Mercator
+ *
+ * @version 1.0 May 2008
+ */
+/** @internal */
+export class ObliqueMercator extends OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(ObliqueMercator.METHOD_CODE, "Oblique Mercator", parameters);
+        /* Store the parameters */
+        this._latC = parameters.getValue(8811);
+        this._lonC = parameters.getValue(8812);
+        this._aziC = parameters.getValue(8813);
+        this._gamC = parameters.getValue(8814);
+        this._kC = parameters.getValue(8815);
+        this._eC = parameters.getValue(8816);
+        this._nC = parameters.getValue(8817);
+    }
+    /**
+       * Get the sign of a number.
+       */
+    static sign(v) {
+        return (v < 0.0) ? (-1.0) : (1.0);
+    }
+    /**
+       * Get the power of a number (must be able to handle negative 'n' values).
+       */
+    static pow(n, e) {
+        let p = 1.0;
+        for (let i = 0; i < e; i++)
+            p *= n;
+        return p;
+    }
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation) {
+        /* Get the ellipsoid */
+        const sourceCRS = operation.getSourceCRS();
+        /* Prepare the forward parameters */
+        this._a = sourceCRS.getEllipsoid().getA();
+        this._e = sourceCRS.getEllipsoid().getE();
+        this._e2 = this._e * this._e;
+        this._B = Math.sqrt(1.0 + (this._e2 * ObliqueMercator.pow(Math.cos(this._latC), 4) / (1.0 - this._e2)));
+        const esinLatC = this._e * Math.sin(this._latC);
+        this._A = this._a * this._B * this._kC * Math.sqrt(1.0 - this._e2) / (1.0 - ObliqueMercator.pow(esinLatC, 2));
+        this._tO = Math.tan(ObliqueMercator.PI / 4.0 - this._latC / 2.0) / Math.pow((1.0 - esinLatC) / (1.0 + esinLatC), this._e / 2.0);
+        this._D = this._B * Math.sqrt(1.0 - this._e2) / (Math.cos(this._latC) * Math.sqrt(1.0 - ObliqueMercator.pow(esinLatC, 2)));
+        const D2 = (this._D < 1.0) ? (1.0) : (this._D * this._D);
+        this._F = this._D + Math.sqrt(D2 - 1.0) * ObliqueMercator.sign(this._latC);
+        this._H = this._F * Math.pow(this._tO, this._B);
+        this._G = (this._F - 1.0 / this._F) / 2.0;
+        this._gamO = Math.asin(Math.sin(this._aziC) / this._D);
+        // fix for Switzerland CRS 2056 projection error. LER, 19/06/2018
+        let w = this._G * Math.tan(this._gamO);
+        if (Math.abs(1.0 - w) < 1.0e-13)
+            w = 1.0; // Math.asin(0.9999999999999999) gives a 1.49e-8 difference from the asin(1.0) value, which leads to a 19 cm projection error
+        //        this._lonO = lonC-Math.asin(G*Math.tan(gamO))/B;
+        this._lonO = this._lonC - Math.asin(w) / this._B;
+        // end of fix
+        this._vC = 0.0;
+        if (Math.abs(this._aziC - 0.5 * ObliqueMercator.PI) < 1.0e-12)
+            this._uC = this._A * (this._lonC - this._lonO); // special case aziC==90 (Hungary, Switzerland)
+        else
+            this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * ObliqueMercator.sign(this._latC);
+        /* Prepare the reverse parameters */
+        const e4 = this._e2 * this._e2;
+        const e6 = e4 * this._e2;
+        const e8 = e4 * e4;
+        this._rev1 = (this._e2 / 2.0 + 5.0 * e4 / 24.0 + e6 / 12.0 + 13.0 * e8 / 360.0);
+        this._rev2 = (7.0 * e4 / 48.0 + 29.0 * e6 / 240.0 + 811.0 * e8 / 11520.0);
+        this._rev3 = (7.0 * e6 / 120.0 + 81.0 * e8 / 1120.0);
+        this._rev4 = (4279.0 * e8 / 161280.0);
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters (radians) */
+        const lon = source.getX();
+        const lat = source.getY();
+        /* Make the calculation */
+        const esinLat = this._e * Math.sin(lat);
+        const t = Math.tan(ObliqueMercator.PI / 4.0 - lat / 2.0) / Math.pow((1.0 - esinLat) / (1.0 + esinLat), this._e / 2.0);
+        const Q = this._H / Math.pow(t, this._B);
+        const S = (Q - 1.0 / Q) / 2.0;
+        const T = (Q + 1.0 / Q) / 2.0;
+        const V = Math.sin(this._B * (lon - this._lonO));
+        const U = (-V * Math.cos(this._gamO) + S * Math.sin(this._gamO)) / T;
+        const v = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);
+        let u;
+        if (Math.abs(this._aziC - 0.5 * ObliqueMercator.PI) < 1.0e-12) // aziC==90? new in guidance note 2 version 8.5. LER, 21/09/2016
+         {
+            if (lon == this._lonC)
+                u = 0.0;
+            else
+                u = (this._A * Math.atan((S * Math.cos(this._gamO) + V * Math.sin(this._gamO)) / Math.cos(this._B * (lon - this._lonO))) / this._B) - (Math.abs(this._uC) * ObliqueMercator.sign(this._latC) * ObliqueMercator.sign(this._lonC - lon));
+        }
+        else {
+            u = (this._A * Math.atan((S * Math.cos(this._gamO) + V * Math.sin(this._gamO)) / Math.cos(this._B * (lon - this._lonO))) / this._B) - (Math.abs(this._uC) * ObliqueMercator.sign(this._latC));
+        }
+        const sinGamC = Math.sin(this._gamC);
+        const cosGamC = Math.cos(this._gamC);
+        const E = v * cosGamC + u * sinGamC + this._eC;
+        const N = u * cosGamC - v * sinGamC + this._nC;
+        /* Save the position */
+        target.setX(E);
+        target.setY(N);
+        target.setZ(source.getZ()); // Keep the Z value
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const E = target.getX();
+        const N = target.getY();
+        /* Make the calculation */
+        const sinGamC = Math.sin(this._gamC);
+        const cosGamC = Math.cos(this._gamC);
+        const v = (E - this._eC) * cosGamC - (N - this._nC) * sinGamC;
+        const u = (N - this._nC) * cosGamC + (E - this._eC) * sinGamC + (Math.abs(this._uC) * ObliqueMercator.sign(this._latC));
+        const Q = Math.exp(-this._B * v / this._A);
+        const S = (Q - 1.0 / Q) / 2.0;
+        const T = (Q + 1.0 / Q) / 2.0;
+        const V = Math.sin(this._B * u / this._A);
+        const U = (V * Math.cos(this._gamO) + S * Math.sin(this._gamO)) / T;
+        const t = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);
+        const chi = ObliqueMercator.PI / 2.0 - 2.0 * Math.atan(t);
+        const lat = chi + Math.sin(2.0 * chi) * this._rev1 + Math.sin(4.0 * chi) * this._rev2 + Math.sin(6.0 * chi) * this._rev3 + Math.sin(8.0 * chi) * this._rev4;
+        //        double lon = lonO-Math.atan((S*Math.cos(gamO)-V*Math.sin(gamO))/Math.cos(B*u/A))/B;
+        const lon = this._lonO - Math.atan2((S * Math.cos(this._gamO) - V * Math.sin(this._gamO)), Math.cos(this._B * u / this._A)) / this._B; // replaced atan(y/x) by atan2(y,x). LER, 23/05/2017
+        /* Save the position */
+        source.setX(lon);
+        source.setY(lat);
+        source.setZ(target.getZ()); // Keep the Z value
+    }
+}
+/** The code of this method */
+ObliqueMercator.METHOD_CODE = 9815; // Hotine Variant B
+/** The value of PI */
+ObliqueMercator.PI = Math.PI;
+//# sourceMappingURL=ObliqueMercator.js.map

package/lib/esm/spatial/ecrs/projection/ObliqueMercator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObliqueMercator.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/ObliqueMercator.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAaH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,eAAgB,SAAQ,eAAe;IA2ChD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,eAAe,CAAC,WAAW,EAAE,kBAAkB,EAAE,UAAU,CAAC,CAAC;QACnE,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU,EAAE,CAAQ;QACnC,IAAI,CAAC,GAAY,GAAG,CAAC;QACrB,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,IAAI,CAAC,CAAC;QAC3C,OAAO,CAAC,CAAC;IACb,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,uBAAuB;QACvB,MAAM,SAAS,GAAQ,SAAS,CAAC,YAAY,EAAE,CAAC;QAChD,oCAAoC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QAC1C,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QAC1C,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC7B,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,eAAe,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACxG,MAAM,QAAQ,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,eAAe,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9G,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;QAChI,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,eAAe,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3H,MAAM,EAAE,GAAY,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAClE,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC3E,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAChD,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;QAC1C,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACvD,iEAAiE;QACjE,IAAI,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAChD,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,OAAO;YAAE,CAAC,GAAG,GAAG,CAAC,CAAE,6HAA6H;QACxK,0DAA0D;QAC1D,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACjD,aAAa;QACb,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,eAAe,CAAC,EAAE,CAAC,GAAG,OAAO;YAAE,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,+CAA+C;;YACzJ,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC/H,oCAAoC;QACpC,MAAM,EAAE,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QACxC,MAAM,EAAE,GAAY,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC;QAClC,MAAM,EAAE,GAAY,EAAE,GAAG,EAAE,CAAC;QAC5B,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;QAChF,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,GAAG,KAAK,GAAG,EAAE,GAAG,OAAO,CAAC,CAAC;QAC1E,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,KAAK,GAAG,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC,CAAC;QACrD,IAAI,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,EAAE,GAAG,QAAQ,CAAC,CAAC;IAC1C,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,kCAAkC;QAClC,MAAM,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACnC,MAAM,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACnC,0BAA0B;QAC1B,MAAM,OAAO,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,eAAe,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;QAC/H,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAClD,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACvC,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACvC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAY,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QAC9E,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAC/E,IAAI,CAAU,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,eAAe,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,gEAAgE;SAC/H;YACI,IAAI,GAAG,IAAI,IAAI,CAAC,KAAK;gBAAE,CAAC,GAAG,GAAG,CAAC;;gBAC1B,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;SAC/O;aAAM;YACH,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;SACjM;QACD,MAAM,OAAO,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9C,MAAM,OAAO,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAY,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC;QACxD,MAAM,CAAC,GAAY,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC;QACxD,uBAAuB;QACvB,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,0BAA0B;QAC1B,MAAM,OAAO,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9C,MAAM,OAAO,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC;QACvE,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QACjI,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACpD,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACvC,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACvC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACnD,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC;QAC7E,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACvF,MAAM,GAAG,GAAY,eAAe,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnE,MAAM,GAAG,GAAY,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrK,6FAA6F;QAC7F,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,oDAAoD;QACpM,uBAAuB;QACvB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;;AAhLD,8BAA8B;AACP,2BAAW,GAAU,IAAI,CAAC,CAAC,mBAAmB;AAErE,sBAAsB;AACE,kBAAE,GAAY,IAAI,CAAC,EAAE,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n// package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { Operation } from \"../Operation\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class ObliqueMercator defines an Oblique Mercator projection (Hotine Variant B).\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised August 2006\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.6 (spec v6.13): Oblique Mercator and Hotine Oblique Mercator\n *\n * @version 1.0 May 2008\n */\n/** @internal */\nexport class ObliqueMercator extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9815; // Hotine Variant B\n\n    /** The value of PI */\n    private static readonly PI: float64 = Math.PI;\n\n    /** Latitude of the projection center */\n    private _latC: float64;\n    /** Longitude of the projection center */\n    private _lonC: float64;\n    /** Azimuth of the initial line passing through the projection center */\n    private _aziC: float64;\n    /** Angle from the recified grid to the skew (oblique) grid */\n    private _gamC: float64;\n    /** Scale factor on the initial line of projection */\n    private _kC: float64;\n    /** False easting at the center of projection */\n    private _eC: float64;\n    /** False northing at the center of projection */\n    private _nC: float64;\n\n    /** Forward calculation parameters */\n    private _a: float64;\n    private _e: float64;\n    private _e2: float64;\n    private _B: float64;\n    private _A: float64;\n    private _tO: float64;\n    private _D: float64;\n    private _F: float64;\n    private _H: float64;\n    private _G: float64;\n    private _gamO: float64;\n    private _lonO: float64;\n    private _vC: float64;\n    private _uC: float64;\n    /** Reverse calculation parameters */\n    private _rev1: float64;\n    private _rev2: float64;\n    private _rev3: float64;\n    private _rev4: float64;\n\n    /**\n       * Create a new projection.\n       * @param parameters the values of the parameters.\n       */\n    public constructor(parameters: ParameterValueList) {\n        super(ObliqueMercator.METHOD_CODE, \"Oblique Mercator\", parameters);\n        /* Store the parameters */\n        this._latC = parameters.getValue(8811);\n        this._lonC = parameters.getValue(8812);\n        this._aziC = parameters.getValue(8813);\n        this._gamC = parameters.getValue(8814);\n        this._kC = parameters.getValue(8815);\n        this._eC = parameters.getValue(8816);\n        this._nC = parameters.getValue(8817);\n    }\n\n    /**\n       * Get the sign of a number.\n       */\n    private static sign(v: float64): float64 {\n        return (v < 0.0) ? (-1.0) : (1.0);\n    }\n\n    /**\n       * Get the power of a number (must be able to handle negative 'n' values).\n       */\n    private static pow(n: float64, e: int32): float64 {\n        let p: float64 = 1.0;\n        for (let i: number = 0; i < e; i++) p *= n;\n        return p;\n    }\n\n    /**\n       * OperationMethod method.\n       * @see OperationMethod#initialize\n       */\n    public override initialize(operation: Operation): void {\n        /* Get the ellipsoid */\n        const sourceCRS: CRS = operation.getSourceCRS();\n        /* Prepare the forward parameters */\n        this._a = sourceCRS.getEllipsoid().getA();\n        this._e = sourceCRS.getEllipsoid().getE();\n        this._e2 = this._e * this._e;\n        this._B = Math.sqrt(1.0 + (this._e2 * ObliqueMercator.pow(Math.cos(this._latC), 4) / (1.0 - this._e2)));\n        const esinLatC: float64 = this._e * Math.sin(this._latC);\n        this._A = this._a * this._B * this._kC * Math.sqrt(1.0 - this._e2) / (1.0 - ObliqueMercator.pow(esinLatC, 2));\n        this._tO = Math.tan(ObliqueMercator.PI / 4.0 - this._latC / 2.0) / Math.pow((1.0 - esinLatC) / (1.0 + esinLatC), this._e / 2.0);\n        this._D = this._B * Math.sqrt(1.0 - this._e2) / (Math.cos(this._latC) * Math.sqrt(1.0 - ObliqueMercator.pow(esinLatC, 2)));\n        const D2: float64 = (this._D < 1.0) ? (1.0) : (this._D * this._D);\n        this._F = this._D + Math.sqrt(D2 - 1.0) * ObliqueMercator.sign(this._latC);\n        this._H = this._F * Math.pow(this._tO, this._B);\n        this._G = (this._F - 1.0 / this._F) / 2.0;\n        this._gamO = Math.asin(Math.sin(this._aziC) / this._D);\n        // fix for Switzerland CRS 2056 projection error. LER, 19/06/2018\n        let w: float64 = this._G * Math.tan(this._gamO);\n        if (Math.abs(1.0 - w) < 1.0e-13) w = 1.0;  // Math.asin(0.9999999999999999) gives a 1.49e-8 difference from the asin(1.0) value, which leads to a 19 cm projection error\n        //        this._lonO = lonC-Math.asin(G*Math.tan(gamO))/B;\n        this._lonO = this._lonC - Math.asin(w) / this._B;\n        // end of fix\n        this._vC = 0.0;\n        if (Math.abs(this._aziC - 0.5 * ObliqueMercator.PI) < 1.0e-12) this._uC = this._A * (this._lonC - this._lonO); // special case aziC==90 (Hungary, Switzerland)\n        else this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * ObliqueMercator.sign(this._latC);\n        /* Prepare the reverse parameters */\n        const e4: float64 = this._e2 * this._e2;\n        const e6: float64 = e4 * this._e2;\n        const e8: float64 = e4 * e4;\n        this._rev1 = (this._e2 / 2.0 + 5.0 * e4 / 24.0 + e6 / 12.0 + 13.0 * e8 / 360.0);\n        this._rev2 = (7.0 * e4 / 48.0 + 29.0 * e6 / 240.0 + 811.0 * e8 / 11520.0);\n        this._rev3 = (7.0 * e6 / 120.0 + 81.0 * e8 / 1120.0);\n        this._rev4 = (4279.0 * e8 / 161280.0);\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#forward\n       */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters (radians) */\n        const lon: float64 = source.getX();\n        const lat: float64 = source.getY();\n        /* Make the calculation */\n        const esinLat: float64 = this._e * Math.sin(lat);\n        const t: float64 = Math.tan(ObliqueMercator.PI / 4.0 - lat / 2.0) / Math.pow((1.0 - esinLat) / (1.0 + esinLat), this._e / 2.0);\n        const Q: float64 = this._H / Math.pow(t, this._B);\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\n        const V: float64 = Math.sin(this._B * (lon - this._lonO));\n        const U: float64 = (-V * Math.cos(this._gamO) + S * Math.sin(this._gamO)) / T;\n        const v: float64 = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);\n        let u: float64;\n        if (Math.abs(this._aziC - 0.5 * ObliqueMercator.PI) < 1.0e-12) // aziC==90? new in guidance note 2 version 8.5. LER, 21/09/2016\n        {\n            if (lon == this._lonC) u = 0.0;\n            else u = (this._A * Math.atan((S * Math.cos(this._gamO) + V * Math.sin(this._gamO)) / Math.cos(this._B * (lon - this._lonO))) / this._B) - (Math.abs(this._uC) * ObliqueMercator.sign(this._latC) * ObliqueMercator.sign(this._lonC - lon));\n        } else {\n            u = (this._A * Math.atan((S * Math.cos(this._gamO) + V * Math.sin(this._gamO)) / Math.cos(this._B * (lon - this._lonO))) / this._B) - (Math.abs(this._uC) * ObliqueMercator.sign(this._latC));\n        }\n        const sinGamC: float64 = Math.sin(this._gamC);\n        const cosGamC: float64 = Math.cos(this._gamC);\n        const E: float64 = v * cosGamC + u * sinGamC + this._eC;\n        const N: float64 = u * cosGamC - v * sinGamC + this._nC;\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#reverse\n       */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = target.getX();\n        const N: float64 = target.getY();\n        /* Make the calculation */\n        const sinGamC: float64 = Math.sin(this._gamC);\n        const cosGamC: float64 = Math.cos(this._gamC);\n        const v: float64 = (E - this._eC) * cosGamC - (N - this._nC) * sinGamC;\n        const u: float64 = (N - this._nC) * cosGamC + (E - this._eC) * sinGamC + (Math.abs(this._uC) * ObliqueMercator.sign(this._latC));\n        const Q: float64 = Math.exp(-this._B * v / this._A);\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\n        const V: float64 = Math.sin(this._B * u / this._A);\n        const U: float64 = (V * Math.cos(this._gamO) + S * Math.sin(this._gamO)) / T;\n        const t: float64 = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);\n        const chi: float64 = ObliqueMercator.PI / 2.0 - 2.0 * Math.atan(t);\n        const lat: float64 = chi + Math.sin(2.0 * chi) * this._rev1 + Math.sin(4.0 * chi) * this._rev2 + Math.sin(6.0 * chi) * this._rev3 + Math.sin(8.0 * chi) * this._rev4;\n        //        double lon = lonO-Math.atan((S*Math.cos(gamO)-V*Math.sin(gamO))/Math.cos(B*u/A))/B;\n        const lon: float64 = this._lonO - Math.atan2((S * Math.cos(this._gamO) - V * Math.sin(this._gamO)), Math.cos(this._B * u / this._A)) / this._B; // replaced atan(y/x) by atan2(y,x). LER, 23/05/2017\n        /* Save the position */\n        source.setX(lon);\n        source.setY(lat);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}

package/lib/esm/spatial/ecrs/projection/ObliqueStereographic.d.ts

@@ -0,0 +1,124 @@
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+declare type int32 = number;
+declare type float64 = number;
+import { Coordinate } from "../../geom/Coordinate";
+import { CRS } from "../CRS";
+import { Ellipsoid } from "../Ellipsoid";
+import { Operation } from "../Operation";
+import { OperationMethod } from "../OperationMethod";
+import { ParameterValueList } from "../ParameterValueList";
+/**
+ * Class ObliqueStereographic defines an Oblique Stereographic projection.
+ *
+ * The 'source' CRS is the geodetic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.7.1 'Oblique and Equatorial Stereographic cases'
+ *
+ * @version 1.0 December 2006
+ */
+/** @internal */
+export declare class ObliqueStereographic extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The value of PI */
+    private static readonly PI;
+    /** Latitude of false origin */
+    private _lat0;
+    /** Longitude of false origin */
+    private _lon0;
+    /** Scale factor at natural origin */
+    private _k0;
+    /** False easting */
+    private _fe;
+    /** False northing */
+    private _fn;
+    /** The parameters for the conformal sphere */
+    private _e;
+    private _e2;
+    private _R;
+    private _n;
+    private _c;
+    private _chi0;
+    private _lambda0;
+    private _g;
+    private _h;
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters: ParameterValueList);
+    /**
+       * Create a new projection.
+       * @param lat0 latitude of natural origin (radians).
+       * @param lon0 longitude of natural origin (radians).
+       * @param k0 scale factor at natural origin.
+       * @param fe false easting.
+       * @param fn false northing.
+       */
+    static create(lat0: float64, lon0: float64, k0: float64, fe: float64, fn: float64): ObliqueStereographic;
+    /**
+       * Get the power of a number.
+       */
+    private static pow;
+    /**
+       * Calculate sinus.
+       */
+    private static sin;
+    /**
+       * Calculate cosinus.
+       */
+    private static cos;
+    /**
+       * Calculate tangent.
+       */
+    private static tan;
+    /**
+       * Calculate arcsinus.
+       */
+    private static asin;
+    /**
+       * Calculate arctangent.
+       */
+    private static atan;
+    /**
+       * Calculate log.
+       */
+    private static log;
+    /**
+       * Calculate exp.
+       */
+    private static exp;
+    /**
+       * Initialize the projection.
+       * @param ellipsoid the ellipsoid to use.
+       * @return this projection (for convenience).
+       */
+    initializeProjection(ellipsoid: Ellipsoid): ObliqueStereographic;
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation: Operation): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+}
+export {};
+//# sourceMappingURL=ObliqueStereographic.d.ts.map

package/lib/esm/spatial/ecrs/projection/ObliqueStereographic.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObliqueStereographic.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/ObliqueStereographic.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,aAAK,KAAK,GAAG,MAAM,CAAC;AAEpB,aAAK,OAAO,GAAG,MAAM,CAAC;AAEtB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,oBAAqB,SAAQ,eAAe;IACrD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,sBAAsB;IACtB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAoB;IAE9C,+BAA+B;IAC/B,OAAO,CAAC,KAAK,CAAU;IACvB,gCAAgC;IAChC,OAAO,CAAC,KAAK,CAAU;IACvB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB,8CAA8C;IAC9C,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,KAAK,CAAU;IACvB,OAAO,CAAC,QAAQ,CAAU;IAC1B,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IAEpB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAUjD;;;;;;;SAOK;WACS,MAAM,CAAC,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,oBAAoB;IAU/G;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,IAAI;IAInB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,IAAI;IAInB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;;;SAIK;IACE,oBAAoB,CAAC,SAAS,EAAE,SAAS,GAAG,oBAAoB;IA2BvE;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAItD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAqB5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CA2B/F"}

package/lib/esm/spatial/ecrs/projection/ObliqueStereographic.js

@@ -0,0 +1,205 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+import { OperationMethod } from "../OperationMethod";
+/**
+ * Class ObliqueStereographic defines an Oblique Stereographic projection.
+ *
+ * The 'source' CRS is the geodetic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.7.1 'Oblique and Equatorial Stereographic cases'
+ *
+ * @version 1.0 December 2006
+ */
+/** @internal */
+export class ObliqueStereographic extends OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(ObliqueStereographic.METHOD_CODE, "Oblique Stereographic", parameters);
+        /* Store the parameters */
+        this._lat0 = (parameters != null) ? parameters.getValue(8801) : 0.0;
+        this._lon0 = (parameters != null) ? parameters.getValue(8802) : 0.0;
+        this._k0 = (parameters != null) ? parameters.getValue(8805) : 0.0;
+        this._fe = (parameters != null) ? parameters.getValue(8806) : 0.0;
+        this._fn = (parameters != null) ? parameters.getValue(8807) : 0.0;
+    }
+    /**
+       * Create a new projection.
+       * @param lat0 latitude of natural origin (radians).
+       * @param lon0 longitude of natural origin (radians).
+       * @param k0 scale factor at natural origin.
+       * @param fe false easting.
+       * @param fn false northing.
+       */
+    static create(lat0, lon0, k0, fe, fn) {
+        const projection = new ObliqueStereographic(null);
+        projection._lat0 = lat0;
+        projection._lon0 = lon0;
+        projection._k0 = k0;
+        projection._fe = fe;
+        projection._fn = fn;
+        return projection;
+    }
+    /**
+       * Get the power of a number.
+       */
+    static pow(n, e) {
+        return Math.exp(e * Math.log(n));
+    }
+    /**
+       * Calculate sinus.
+       */
+    static sin(a) {
+        return Math.sin(a);
+    }
+    /**
+       * Calculate cosinus.
+       */
+    static cos(a) {
+        return Math.cos(a);
+    }
+    /**
+       * Calculate tangent.
+       */
+    static tan(a) {
+        return Math.tan(a);
+    }
+    /**
+       * Calculate arcsinus.
+       */
+    static asin(v) {
+        return Math.asin(v);
+    }
+    /**
+       * Calculate arctangent.
+       */
+    static atan(v) {
+        return Math.atan(v);
+    }
+    /**
+       * Calculate log.
+       */
+    static log(v) {
+        return Math.log(v);
+    }
+    /**
+       * Calculate exp.
+       */
+    static exp(v) {
+        return Math.exp(v);
+    }
+    /**
+       * Initialize the projection.
+       * @param ellipsoid the ellipsoid to use.
+       * @return this projection (for convenience).
+       */
+    initializeProjection(ellipsoid) {
+        /* Get the ellipsoid parameters */
+        const a = ellipsoid.getA();
+        this._e = ellipsoid.getE();
+        /* Calculate conformat latitude (chi0) and longitude (lambda0) */
+        this._e2 = this._e * this._e;
+        const sinLat0 = ObliqueStereographic.sin(this._lat0);
+        const t0 = 1.0 - this._e2 * sinLat0 * sinLat0;
+        const rho0 = a * (1 - this._e2) / ObliqueStereographic.pow(t0, 1.5);
+        const nu0 = a / ObliqueStereographic.pow(t0, 0.5);
+        this._R = ObliqueStereographic.pow(rho0 * nu0, 0.5);
+        this._n = ObliqueStereographic.pow(1.0 + (this._e2 * ObliqueStereographic.pow(ObliqueStereographic.cos(this._lat0), 4)) / (1.0 - this._e2), 0.5);
+        const S1 = (1.0 + sinLat0) / (1.0 - sinLat0);
+        const S2 = (1.0 - this._e * sinLat0) / (1.0 + this._e * sinLat0);
+        const w1 = ObliqueStereographic.pow(S1 * ObliqueStereographic.pow(S2, this._e), this._n);
+        const sinChi0 = (w1 - 1.0) / (w1 + 1.0);
+        this._c = (this._n + sinLat0) * (1.0 - sinChi0) / ((this._n - sinLat0) * (1.0 + sinChi0));
+        const w2 = this._c * w1;
+        this._chi0 = ObliqueStereographic.asin((w2 - 1.0) / (w2 + 1.0));
+        this._lambda0 = this._lon0;
+        /* For reverse calculation */
+        this._g = 2.0 * this._R * this._k0 * ObliqueStereographic.tan(ObliqueStereographic.PI / 4.0 - this._chi0 / 2.0);
+        this._h = 4.0 * this._R * this._k0 * ObliqueStereographic.tan(this._chi0) + this._g;
+        /* Return the projection */
+        return this;
+    }
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation) {
+        this.initializeProjection(operation.getSourceCRS().getEllipsoid());
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const lon = source.getX();
+        const lat = source.getY();
+        const sinLat = ObliqueStereographic.sin(lat);
+        /* Calculate conformal latitude and longitude */
+        const lambda = this._n * (lon - this._lambda0) + this._lambda0;
+        const Sa = (1.0 + sinLat) / (1.0 - sinLat);
+        const Sb = (1.0 - this._e * sinLat) / (1.0 + this._e * sinLat);
+        const w = this._c * ObliqueStereographic.pow(Sa * ObliqueStereographic.pow(Sb, this._e), this._n);
+        const chi = ObliqueStereographic.asin((w - 1.0) / (w + 1.0));
+        /* Calculate easting and northing */
+        const B = 1.0 + ObliqueStereographic.sin(chi) * ObliqueStereographic.sin(this._chi0) + ObliqueStereographic.cos(chi) * ObliqueStereographic.cos(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0);
+        const E = this._fe + 2.0 * this._R * this._k0 * ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(lambda - this._lambda0) / B;
+        const N = this._fn + 2.0 * this._R * this._k0 * (ObliqueStereographic.sin(chi) * ObliqueStereographic.cos(this._chi0) - ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0)) / B;
+        /* Save the position */
+        target.setX(E);
+        target.setY(N);
+        target.setZ(source.getZ()); // Keep the Z value
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const E = target.getX() - this._fe;
+        const N = target.getY() - this._fn;
+        /* Calculate conformal latitude and longitude */
+        const i = ObliqueStereographic.atan(E / (this._h + N));
+        const j = ObliqueStereographic.atan(E / (this._g - N)) - i;
+        const chi = this._chi0 + 2.0 * ObliqueStereographic.atan((N - E * ObliqueStereographic.tan(j / 2.0)) / (2.0 * this._R * this._k0));
+        const lambda = j + 2.0 * i + this._lambda0;
+        const lon = (lambda - this._lambda0) / this._n + this._lambda0;
+        /* Iterate to get isometric (psi) and geodetic (phi) latitude */
+        const sinChi = ObliqueStereographic.sin(chi);
+        const psi = 0.5 * ObliqueStereographic.log((1.0 + sinChi) / (this._c * (1.0 - sinChi))) / this._n;
+        let phiK = 2.0 * ObliqueStereographic.atan(ObliqueStereographic.exp(psi)) - ObliqueStereographic.PI / 2.0;
+        for (let k = 1; k < 8; k++) {
+            const esinPhi = this._e * ObliqueStereographic.sin(phiK);
+            const psiK = ObliqueStereographic.log(ObliqueStereographic.tan(phiK / 2.0 + ObliqueStereographic.PI / 4.0) * ObliqueStereographic.pow((1.0 - esinPhi) / (1.0 + esinPhi), this._e / 2.0));
+            const phiNextK = phiK - (psiK - psi) * ObliqueStereographic.cos(phiK) * (1.0 - esinPhi * esinPhi) / (1.0 - this._e2);
+            const change = Math.abs(phiNextK - phiK);
+            phiK = phiNextK;
+            if (change < 1.0e-12)
+                break; // normally after 4 iterations
+        }
+        /* Save the position */
+        source.setX(lon);
+        source.setY(phiK);
+        source.setZ(target.getZ()); // Keep the Z value
+    }
+}
+/** The code of this method */
+ObliqueStereographic.METHOD_CODE = 9809;
+/** The value of PI */
+ObliqueStereographic.PI = Math.PI;
+//# sourceMappingURL=ObliqueStereographic.js.map