@itwin/core-orbitgt 4.0.0-dev.8 → 4.0.0-dev.80

package/lib/cjs/spatial/ecrs/projection/HotineObliqueMercator.js

@@ -1,163 +1,163 @@
-"use strict";
-/*---------------------------------------------------------------------------------------------
-* 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
- */
-Object.defineProperty(exports, "__esModule", { value: true });
-exports.HotineObliqueMercator = void 0;
-const OperationMethod_1 = require("../OperationMethod");
-/**
- * Class HotineObliqueMercator defines a Hotine Oblique Mercator projection.
- *
- * 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: Oblique Mercator and Hotine Oblique Mercator
- *
- * @version 1.0 May 2008
- */
-/** @internal */
-class HotineObliqueMercator extends OperationMethod_1.OperationMethod {
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters) {
-        super(HotineObliqueMercator.METHOD_CODE, "Hotine 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._fe = parameters.getValue(8806);
-        this._fn = parameters.getValue(8807);
-    }
-    /**
-       * 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;
-    }
-    /**
-       * Initialize the projection.
-       * @param ellipsoid the ellipsoid to use.
-       * @return this projection (for convenience).
-       */
-    initializeProjection(ellipsoid) {
-        /* Prepare the forward parameters */
-        this._a = ellipsoid.getA();
-        this._e = ellipsoid.getE();
-        this._e2 = this._e * this._e;
-        this._B = Math.sqrt(1.0 + (this._e2 * HotineObliqueMercator.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 - HotineObliqueMercator.pow(esinLatC, 2));
-        this._tO = Math.tan(Math.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 - HotineObliqueMercator.pow(esinLatC, 2)));
-        const D2 = (this._D < 1.0) ? (1.0) : (this._D * this._D);
-        this._F = this._D + Math.sqrt(D2 - 1.0) * HotineObliqueMercator.sign(this._latC);
-        this._H = this._F * Math.pow(this._tO, this._B);
-        this._G = (this._F - 1.0 / this._F) / 2.0;
-        this._latO = Math.asin(Math.sin(this._aziC) / this._D);
-        this._lonO = this._lonC - Math.asin(this._G * Math.tan(this._latO)) / this._B;
-        this._vC = 0.0;
-        if (Math.abs(this._aziC - 0.5 * Math.PI) < 0.00001)
-            this._uC = this._A * (this._lonC - this._lonO); // special case Hungary, Switzerland
-        else
-            this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * HotineObliqueMercator.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 + 5 * e4 / 24 + e6 / 12 + 13 * e8 / 360);
-        this._rev2 = (7 * e4 / 48 + 29 * e6 / 240 + 811 * e8 / 11520);
-        this._rev3 = (7 * e6 / 120 + 81 * e8 / 1120);
-        this._rev4 = (4279 * e8 / 161280);
-        /* 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();
-        /* Make the calculation */
-        const esinLat = this._e * Math.sin(lat);
-        const t = Math.tan(Math.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._latO) + S * Math.sin(this._latO)) / T;
-        const v = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);
-        const u = (this._A * Math.atan((S * Math.cos(this._latO) + V * Math.sin(this._latO)) / Math.cos(this._B * (lon - this._lonO))) / this._B); // possibly related to method 9815 atan2 problem? (LER 15/06/2018)
-        const sinGamC = Math.sin(this._gamC);
-        const cosGamC = Math.cos(this._gamC);
-        const E = v * cosGamC + u * sinGamC + this._fe;
-        const N = u * cosGamC - v * sinGamC + this._fn;
-        /* 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._fe) * cosGamC - (N - this._fn) * sinGamC;
-        const u = (N - this._fn) * cosGamC + (E - this._fe) * sinGamC;
-        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._latO) + S * Math.sin(this._latO)) / T;
-        const t = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);
-        const chi = Math.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;
-        const lon = this._lonO - Math.atan((S * cosGamC - V * sinGamC) / Math.cos(this._B * u / this._A)) / this._B; // possibly related to method 9815 atan2 problem? (LER 15/06/2018)
-        /* Save the position */
-        source.setX(lon);
-        source.setY(lat);
-        source.setZ(target.getZ()); // Keep the Z value
-    }
-}
-exports.HotineObliqueMercator = HotineObliqueMercator;
-/** The code of this method */
-HotineObliqueMercator.METHOD_CODE = 9812;
+"use strict";
+/*---------------------------------------------------------------------------------------------
+* 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
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.HotineObliqueMercator = void 0;
+const OperationMethod_1 = require("../OperationMethod");
+/**
+ * Class HotineObliqueMercator defines a Hotine Oblique Mercator projection.
+ *
+ * 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: Oblique Mercator and Hotine Oblique Mercator
+ *
+ * @version 1.0 May 2008
+ */
+/** @internal */
+class HotineObliqueMercator extends OperationMethod_1.OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(HotineObliqueMercator.METHOD_CODE, "Hotine 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._fe = parameters.getValue(8806);
+        this._fn = parameters.getValue(8807);
+    }
+    /**
+       * 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;
+    }
+    /**
+       * Initialize the projection.
+       * @param ellipsoid the ellipsoid to use.
+       * @return this projection (for convenience).
+       */
+    initializeProjection(ellipsoid) {
+        /* Prepare the forward parameters */
+        this._a = ellipsoid.getA();
+        this._e = ellipsoid.getE();
+        this._e2 = this._e * this._e;
+        this._B = Math.sqrt(1.0 + (this._e2 * HotineObliqueMercator.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 - HotineObliqueMercator.pow(esinLatC, 2));
+        this._tO = Math.tan(Math.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 - HotineObliqueMercator.pow(esinLatC, 2)));
+        const D2 = (this._D < 1.0) ? (1.0) : (this._D * this._D);
+        this._F = this._D + Math.sqrt(D2 - 1.0) * HotineObliqueMercator.sign(this._latC);
+        this._H = this._F * Math.pow(this._tO, this._B);
+        this._G = (this._F - 1.0 / this._F) / 2.0;
+        this._latO = Math.asin(Math.sin(this._aziC) / this._D);
+        this._lonO = this._lonC - Math.asin(this._G * Math.tan(this._latO)) / this._B;
+        this._vC = 0.0;
+        if (Math.abs(this._aziC - 0.5 * Math.PI) < 0.00001)
+            this._uC = this._A * (this._lonC - this._lonO); // special case Hungary, Switzerland
+        else
+            this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * HotineObliqueMercator.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 + 5 * e4 / 24 + e6 / 12 + 13 * e8 / 360);
+        this._rev2 = (7 * e4 / 48 + 29 * e6 / 240 + 811 * e8 / 11520);
+        this._rev3 = (7 * e6 / 120 + 81 * e8 / 1120);
+        this._rev4 = (4279 * e8 / 161280);
+        /* 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();
+        /* Make the calculation */
+        const esinLat = this._e * Math.sin(lat);
+        const t = Math.tan(Math.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._latO) + S * Math.sin(this._latO)) / T;
+        const v = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);
+        const u = (this._A * Math.atan((S * Math.cos(this._latO) + V * Math.sin(this._latO)) / Math.cos(this._B * (lon - this._lonO))) / this._B); // possibly related to method 9815 atan2 problem? (LER 15/06/2018)
+        const sinGamC = Math.sin(this._gamC);
+        const cosGamC = Math.cos(this._gamC);
+        const E = v * cosGamC + u * sinGamC + this._fe;
+        const N = u * cosGamC - v * sinGamC + this._fn;
+        /* 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._fe) * cosGamC - (N - this._fn) * sinGamC;
+        const u = (N - this._fn) * cosGamC + (E - this._fe) * sinGamC;
+        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._latO) + S * Math.sin(this._latO)) / T;
+        const t = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);
+        const chi = Math.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;
+        const lon = this._lonO - Math.atan((S * cosGamC - V * sinGamC) / Math.cos(this._B * u / this._A)) / this._B; // possibly related to method 9815 atan2 problem? (LER 15/06/2018)
+        /* Save the position */
+        source.setX(lon);
+        source.setY(lat);
+        source.setZ(target.getZ()); // Keep the Z value
+    }
+}
+/** The code of this method */
+HotineObliqueMercator.METHOD_CODE = 9812;
+exports.HotineObliqueMercator = HotineObliqueMercator;
 //# sourceMappingURL=HotineObliqueMercator.js.map

package/lib/cjs/spatial/ecrs/projection/HotineObliqueMercator.js.map

@@ -1 +1 @@
-{"version":3,"file":"HotineObliqueMercator.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/HotineObliqueMercator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAcH,wDAAqD;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAa,qBAAsB,SAAQ,iCAAe;IAwCtD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,qBAAqB,CAAC,WAAW,EAAE,yBAAyB,EAAE,UAAU,CAAC,CAAC;QAChF,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;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,oCAAoC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,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,qBAAqB,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;QAC9G,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,qBAAqB,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC;QACpH,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,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;QACrH,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,qBAAqB,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACjI,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,qBAAqB,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjF,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,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC9E,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,IAAI,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,oCAAoC;;YACnI,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,qBAAqB,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACrI,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,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;QACpE,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;QAC9D,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC;QAC7C,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC,CAAC;QAClC,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,IAAI,CAAC,oBAAoB,CAAC,SAAS,CAAC,YAAY,EAAE,CAAC,YAAY,EAAE,CAAC,CAAC;IACvE,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,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,IAAI,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;QACpH,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,MAAM,CAAC,GAAY,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,CAAC,CAAC,kEAAkE;QACtN,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,CAAC;QACvE,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,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,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,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,kEAAkE;QACxL,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;;AA1KL,sDA2KC;AA1KG,8BAA8B;AACP,iCAAW,GAAU,IAAI,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n// package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { Ellipsoid } from \"../Ellipsoid\";\r\nimport { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class HotineObliqueMercator defines a Hotine Oblique Mercator projection.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised August 2006\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.6: Oblique Mercator and Hotine Oblique Mercator\r\n *\r\n * @version 1.0 May 2008\r\n */\r\n/** @internal */\r\nexport class HotineObliqueMercator extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9812;\r\n\r\n    /** Latitude of the projection center */\r\n    private _latC: float64;\r\n    /** Longitude of the projection center */\r\n    private _lonC: float64;\r\n    /** Azimuth of the initial line passing through the projection center */\r\n    private _aziC: float64;\r\n    /** Angle from the recified grid to the skew (oblique) grid */\r\n    private _gamC: float64;\r\n    /** Scale factor on the initial line of projection */\r\n    private _kC: float64;\r\n    /** False easting at the natural origin */\r\n    private _fe: float64;\r\n    /** False northing at the natural origin */\r\n    private _fn: float64;\r\n\r\n    /** Forward calculation parameters */\r\n    private _a: float64;\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _B: float64;\r\n    private _A: float64;\r\n    private _tO: float64;\r\n    private _D: float64;\r\n    private _F: float64;\r\n    private _H: float64;\r\n    private _G: float64;\r\n    private _latO: float64;\r\n    private _lonO: float64;\r\n    private _vC: float64;\r\n    private _uC: float64;\r\n    /** Reverse calculation parameters */\r\n    private _rev1: float64;\r\n    private _rev2: float64;\r\n    private _rev3: float64;\r\n    private _rev4: float64;\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param parameters the values of the parameters.\r\n       */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(HotineObliqueMercator.METHOD_CODE, \"Hotine Oblique Mercator\", parameters);\r\n        /* Store the parameters */\r\n        this._latC = parameters.getValue(8811);\r\n        this._lonC = parameters.getValue(8812);\r\n        this._aziC = parameters.getValue(8813);\r\n        this._gamC = parameters.getValue(8814);\r\n        this._kC = parameters.getValue(8815);\r\n        this._fe = parameters.getValue(8806);\r\n        this._fn = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the sign of a number.\r\n       */\r\n    private static sign(v: float64): float64 {\r\n        return (v < 0.0) ? (-1.0) : (1.0);\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number (must be able to handle negative 'n' values).\r\n       */\r\n    private static pow(n: float64, e: int32): float64 {\r\n        let p: float64 = 1.0;\r\n        for (let i: number = 0; i < e; i++) p *= n;\r\n        return p;\r\n    }\r\n\r\n    /**\r\n       * Initialize the projection.\r\n       * @param ellipsoid the ellipsoid to use.\r\n       * @return this projection (for convenience).\r\n       */\r\n    public initializeProjection(ellipsoid: Ellipsoid): HotineObliqueMercator {\r\n        /* Prepare the forward parameters */\r\n        this._a = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        this._e2 = this._e * this._e;\r\n        this._B = Math.sqrt(1.0 + (this._e2 * HotineObliqueMercator.pow(Math.cos(this._latC), 4) / (1.0 - this._e2)));\r\n        const esinLatC: float64 = this._e * Math.sin(this._latC);\r\n        this._A = this._a * this._B * this._kC * Math.sqrt(1.0 - this._e2) / (1.0 - HotineObliqueMercator.pow(esinLatC, 2));\r\n        this._tO = Math.tan(Math.PI / 4.0 - this._latC / 2.0) / Math.pow((1.0 - esinLatC) / (1.0 + esinLatC), this._e / 2.0);\r\n        this._D = this._B * Math.sqrt(1.0 - this._e2) / (Math.cos(this._latC) * Math.sqrt(1.0 - HotineObliqueMercator.pow(esinLatC, 2)));\r\n        const D2: float64 = (this._D < 1.0) ? (1.0) : (this._D * this._D);\r\n        this._F = this._D + Math.sqrt(D2 - 1.0) * HotineObliqueMercator.sign(this._latC);\r\n        this._H = this._F * Math.pow(this._tO, this._B);\r\n        this._G = (this._F - 1.0 / this._F) / 2.0;\r\n        this._latO = Math.asin(Math.sin(this._aziC) / this._D);\r\n        this._lonO = this._lonC - Math.asin(this._G * Math.tan(this._latO)) / this._B;\r\n        this._vC = 0.0;\r\n        if (Math.abs(this._aziC - 0.5 * Math.PI) < 0.00001) this._uC = this._A * (this._lonC - this._lonO); // special case Hungary, Switzerland\r\n        else this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * HotineObliqueMercator.sign(this._latC);\r\n        /* Prepare the reverse parameters */\r\n        const e4: float64 = this._e2 * this._e2;\r\n        const e6: float64 = e4 * this._e2;\r\n        const e8: float64 = e4 * e4;\r\n        this._rev1 = (this._e2 / 2 + 5 * e4 / 24 + e6 / 12 + 13 * e8 / 360);\r\n        this._rev2 = (7 * e4 / 48 + 29 * e6 / 240 + 811 * e8 / 11520);\r\n        this._rev3 = (7 * e6 / 120 + 81 * e8 / 1120);\r\n        this._rev4 = (4279 * e8 / 161280);\r\n        /* Return the projection */\r\n        return this;\r\n    }\r\n\r\n    /**\r\n       * OperationMethod method.\r\n       * @see OperationMethod#initialize\r\n       */\r\n    public override initialize(operation: Operation): void {\r\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#forward\r\n       */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const lon: float64 = source.getX();\r\n        const lat: float64 = source.getY();\r\n        /* Make the calculation */\r\n        const esinLat: float64 = this._e * Math.sin(lat);\r\n        const t: float64 = Math.tan(Math.PI / 4.0 - lat / 2.0) / Math.pow((1.0 - esinLat) / (1.0 + esinLat), this._e / 2.0);\r\n        const Q: float64 = this._H / Math.pow(t, this._B);\r\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\r\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\r\n        const V: float64 = Math.sin(this._B * (lon - this._lonO));\r\n        const U: float64 = (-V * Math.cos(this._latO) + S * Math.sin(this._latO)) / T;\r\n        const v: float64 = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);\r\n        const u: float64 = (this._A * Math.atan((S * Math.cos(this._latO) + V * Math.sin(this._latO)) / Math.cos(this._B * (lon - this._lonO))) / this._B); // possibly related to method 9815 atan2 problem? (LER 15/06/2018)\r\n        const sinGamC: float64 = Math.sin(this._gamC);\r\n        const cosGamC: float64 = Math.cos(this._gamC);\r\n        const E: float64 = v * cosGamC + u * sinGamC + this._fe;\r\n        const N: float64 = u * cosGamC - v * sinGamC + this._fn;\r\n        /* Save the position */\r\n        target.setX(E);\r\n        target.setY(N);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#reverse\r\n       */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Make the calculation */\r\n        const sinGamC: float64 = Math.sin(this._gamC);\r\n        const cosGamC: float64 = Math.cos(this._gamC);\r\n        const v: float64 = (E - this._fe) * cosGamC - (N - this._fn) * sinGamC;\r\n        const u: float64 = (N - this._fn) * cosGamC + (E - this._fe) * sinGamC;\r\n        const Q: float64 = Math.exp(-this._B * v / this._A);\r\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\r\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\r\n        const V: float64 = Math.sin(this._B * u / this._A);\r\n        const U: float64 = (V * Math.cos(this._latO) + S * Math.sin(this._latO)) / T;\r\n        const t: float64 = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);\r\n        const chi: float64 = Math.PI / 2.0 - 2.0 * Math.atan(t);\r\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;\r\n        const lon: float64 = this._lonO - Math.atan((S * cosGamC - V * sinGamC) / Math.cos(this._B * u / this._A)) / this._B; // possibly related to method 9815 atan2 problem? (LER 15/06/2018)\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(lat);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"HotineObliqueMercator.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/HotineObliqueMercator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAcH,wDAAqD;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAa,qBAAsB,SAAQ,iCAAe;IAwCtD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,qBAAqB,CAAC,WAAW,EAAE,yBAAyB,EAAE,UAAU,CAAC,CAAC;QAChF,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;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,oCAAoC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,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,qBAAqB,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;QAC9G,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,qBAAqB,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC;QACpH,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,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;QACrH,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,qBAAqB,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACjI,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,qBAAqB,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACjF,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,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC9E,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,IAAI,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,oCAAoC;;YACnI,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,qBAAqB,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACrI,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,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;QACpE,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;QAC9D,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC;QAC7C,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC,CAAC;QAClC,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,IAAI,CAAC,oBAAoB,CAAC,SAAS,CAAC,YAAY,EAAE,CAAC,YAAY,EAAE,CAAC,CAAC;IACvE,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,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,IAAI,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;QACpH,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,MAAM,CAAC,GAAY,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,CAAC,CAAC,kEAAkE;QACtN,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,CAAC;QACvE,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,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxD,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,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,GAAG,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,kEAAkE;QACxL,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;;AAzKD,8BAA8B;AACP,iCAAW,GAAU,IAAI,CAAC;AAFxC,sDAAqB","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n// package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { Ellipsoid } from \"../Ellipsoid\";\r\nimport { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class HotineObliqueMercator defines a Hotine Oblique Mercator projection.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised August 2006\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.6: Oblique Mercator and Hotine Oblique Mercator\r\n *\r\n * @version 1.0 May 2008\r\n */\r\n/** @internal */\r\nexport class HotineObliqueMercator extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9812;\r\n\r\n    /** Latitude of the projection center */\r\n    private _latC: float64;\r\n    /** Longitude of the projection center */\r\n    private _lonC: float64;\r\n    /** Azimuth of the initial line passing through the projection center */\r\n    private _aziC: float64;\r\n    /** Angle from the recified grid to the skew (oblique) grid */\r\n    private _gamC: float64;\r\n    /** Scale factor on the initial line of projection */\r\n    private _kC: float64;\r\n    /** False easting at the natural origin */\r\n    private _fe: float64;\r\n    /** False northing at the natural origin */\r\n    private _fn: float64;\r\n\r\n    /** Forward calculation parameters */\r\n    private _a: float64;\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _B: float64;\r\n    private _A: float64;\r\n    private _tO: float64;\r\n    private _D: float64;\r\n    private _F: float64;\r\n    private _H: float64;\r\n    private _G: float64;\r\n    private _latO: float64;\r\n    private _lonO: float64;\r\n    private _vC: float64;\r\n    private _uC: float64;\r\n    /** Reverse calculation parameters */\r\n    private _rev1: float64;\r\n    private _rev2: float64;\r\n    private _rev3: float64;\r\n    private _rev4: float64;\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param parameters the values of the parameters.\r\n       */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(HotineObliqueMercator.METHOD_CODE, \"Hotine Oblique Mercator\", parameters);\r\n        /* Store the parameters */\r\n        this._latC = parameters.getValue(8811);\r\n        this._lonC = parameters.getValue(8812);\r\n        this._aziC = parameters.getValue(8813);\r\n        this._gamC = parameters.getValue(8814);\r\n        this._kC = parameters.getValue(8815);\r\n        this._fe = parameters.getValue(8806);\r\n        this._fn = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the sign of a number.\r\n       */\r\n    private static sign(v: float64): float64 {\r\n        return (v < 0.0) ? (-1.0) : (1.0);\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number (must be able to handle negative 'n' values).\r\n       */\r\n    private static pow(n: float64, e: int32): float64 {\r\n        let p: float64 = 1.0;\r\n        for (let i: number = 0; i < e; i++) p *= n;\r\n        return p;\r\n    }\r\n\r\n    /**\r\n       * Initialize the projection.\r\n       * @param ellipsoid the ellipsoid to use.\r\n       * @return this projection (for convenience).\r\n       */\r\n    public initializeProjection(ellipsoid: Ellipsoid): HotineObliqueMercator {\r\n        /* Prepare the forward parameters */\r\n        this._a = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        this._e2 = this._e * this._e;\r\n        this._B = Math.sqrt(1.0 + (this._e2 * HotineObliqueMercator.pow(Math.cos(this._latC), 4) / (1.0 - this._e2)));\r\n        const esinLatC: float64 = this._e * Math.sin(this._latC);\r\n        this._A = this._a * this._B * this._kC * Math.sqrt(1.0 - this._e2) / (1.0 - HotineObliqueMercator.pow(esinLatC, 2));\r\n        this._tO = Math.tan(Math.PI / 4.0 - this._latC / 2.0) / Math.pow((1.0 - esinLatC) / (1.0 + esinLatC), this._e / 2.0);\r\n        this._D = this._B * Math.sqrt(1.0 - this._e2) / (Math.cos(this._latC) * Math.sqrt(1.0 - HotineObliqueMercator.pow(esinLatC, 2)));\r\n        const D2: float64 = (this._D < 1.0) ? (1.0) : (this._D * this._D);\r\n        this._F = this._D + Math.sqrt(D2 - 1.0) * HotineObliqueMercator.sign(this._latC);\r\n        this._H = this._F * Math.pow(this._tO, this._B);\r\n        this._G = (this._F - 1.0 / this._F) / 2.0;\r\n        this._latO = Math.asin(Math.sin(this._aziC) / this._D);\r\n        this._lonO = this._lonC - Math.asin(this._G * Math.tan(this._latO)) / this._B;\r\n        this._vC = 0.0;\r\n        if (Math.abs(this._aziC - 0.5 * Math.PI) < 0.00001) this._uC = this._A * (this._lonC - this._lonO); // special case Hungary, Switzerland\r\n        else this._uC = (this._A / this._B) * Math.atan(Math.sqrt(D2 - 1.0) / Math.cos(this._aziC)) * HotineObliqueMercator.sign(this._latC);\r\n        /* Prepare the reverse parameters */\r\n        const e4: float64 = this._e2 * this._e2;\r\n        const e6: float64 = e4 * this._e2;\r\n        const e8: float64 = e4 * e4;\r\n        this._rev1 = (this._e2 / 2 + 5 * e4 / 24 + e6 / 12 + 13 * e8 / 360);\r\n        this._rev2 = (7 * e4 / 48 + 29 * e6 / 240 + 811 * e8 / 11520);\r\n        this._rev3 = (7 * e6 / 120 + 81 * e8 / 1120);\r\n        this._rev4 = (4279 * e8 / 161280);\r\n        /* Return the projection */\r\n        return this;\r\n    }\r\n\r\n    /**\r\n       * OperationMethod method.\r\n       * @see OperationMethod#initialize\r\n       */\r\n    public override initialize(operation: Operation): void {\r\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#forward\r\n       */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const lon: float64 = source.getX();\r\n        const lat: float64 = source.getY();\r\n        /* Make the calculation */\r\n        const esinLat: float64 = this._e * Math.sin(lat);\r\n        const t: float64 = Math.tan(Math.PI / 4.0 - lat / 2.0) / Math.pow((1.0 - esinLat) / (1.0 + esinLat), this._e / 2.0);\r\n        const Q: float64 = this._H / Math.pow(t, this._B);\r\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\r\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\r\n        const V: float64 = Math.sin(this._B * (lon - this._lonO));\r\n        const U: float64 = (-V * Math.cos(this._latO) + S * Math.sin(this._latO)) / T;\r\n        const v: float64 = this._A * Math.log((1.0 - U) / (1.0 + U)) / (2.0 * this._B);\r\n        const u: float64 = (this._A * Math.atan((S * Math.cos(this._latO) + V * Math.sin(this._latO)) / Math.cos(this._B * (lon - this._lonO))) / this._B); // possibly related to method 9815 atan2 problem? (LER 15/06/2018)\r\n        const sinGamC: float64 = Math.sin(this._gamC);\r\n        const cosGamC: float64 = Math.cos(this._gamC);\r\n        const E: float64 = v * cosGamC + u * sinGamC + this._fe;\r\n        const N: float64 = u * cosGamC - v * sinGamC + this._fn;\r\n        /* Save the position */\r\n        target.setX(E);\r\n        target.setY(N);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#reverse\r\n       */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Make the calculation */\r\n        const sinGamC: float64 = Math.sin(this._gamC);\r\n        const cosGamC: float64 = Math.cos(this._gamC);\r\n        const v: float64 = (E - this._fe) * cosGamC - (N - this._fn) * sinGamC;\r\n        const u: float64 = (N - this._fn) * cosGamC + (E - this._fe) * sinGamC;\r\n        const Q: float64 = Math.exp(-this._B * v / this._A);\r\n        const S: float64 = (Q - 1.0 / Q) / 2.0;\r\n        const T: float64 = (Q + 1.0 / Q) / 2.0;\r\n        const V: float64 = Math.sin(this._B * u / this._A);\r\n        const U: float64 = (V * Math.cos(this._latO) + S * Math.sin(this._latO)) / T;\r\n        const t: float64 = Math.pow(this._H / Math.sqrt((1.0 + U) / (1.0 - U)), 1.0 / this._B);\r\n        const chi: float64 = Math.PI / 2.0 - 2.0 * Math.atan(t);\r\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;\r\n        const lon: float64 = this._lonO - Math.atan((S * cosGamC - V * sinGamC) / Math.cos(this._B * u / this._A)) / this._B; // possibly related to method 9815 atan2 problem? (LER 15/06/2018)\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(lat);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}

package/lib/cjs/spatial/ecrs/projection/KrovakObliqueConformalConic.d.ts

@@ -1,80 +1,80 @@
-/** @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 KrovakObliqueConformalConic defines the Krovak Oblique Conformal Conic projection as used in the Czech Republic and Slovakia.
- *
- * 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 (version 6.18.3)
- * Revised January 2009
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.3.2 "Krovak Oblique Conformal Conic"
- *
- * NOTE: the X axis runs north to south and the Y axis runs east to west. This differs from most other projections.
- *
- * @version 1.0 October 2009
- */
-/** @internal */
-export declare class KrovakObliqueConformalConic extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** Latitude of projection center */
-    private _latC;
-    /** Longitude of origin */
-    private _lonO;
-    /** Azimuth of initial line through the projection center */
-    private _aziC;
-    /** Latitude of pseudo standard parallel */
-    private _latP;
-    /** Scale factor on pseudo standard parallel */
-    private _kP;
-    /** Easting at grid origin */
-    private _FE;
-    /** Northing at grid origin */
-    private _FN;
-    /** The constants */
-    private _A;
-    private _B;
-    private _gO;
-    private _tO;
-    private _n;
-    private _rO;
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters: ParameterValueList);
-    /**
-       * Get the square of a number.
-       */
-    private static sqr;
-    /**
-       * 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 {};
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+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 KrovakObliqueConformalConic defines the Krovak Oblique Conformal Conic projection as used in the Czech Republic and Slovakia.
+ *
+ * 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 (version 6.18.3)
+ * Revised January 2009
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.3.2 "Krovak Oblique Conformal Conic"
+ *
+ * NOTE: the X axis runs north to south and the Y axis runs east to west. This differs from most other projections.
+ *
+ * @version 1.0 October 2009
+ */
+/** @internal */
+export declare class KrovakObliqueConformalConic extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** Latitude of projection center */
+    private _latC;
+    /** Longitude of origin */
+    private _lonO;
+    /** Azimuth of initial line through the projection center */
+    private _aziC;
+    /** Latitude of pseudo standard parallel */
+    private _latP;
+    /** Scale factor on pseudo standard parallel */
+    private _kP;
+    /** Easting at grid origin */
+    private _FE;
+    /** Northing at grid origin */
+    private _FN;
+    /** The constants */
+    private _A;
+    private _B;
+    private _gO;
+    private _tO;
+    private _n;
+    private _rO;
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters: ParameterValueList);
+    /**
+       * Get the square of a number.
+       */
+    private static sqr;
+    /**
+       * 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=KrovakObliqueConformalConic.d.ts.map

package/lib/cjs/spatial/ecrs/projection/KrovakObliqueConformalConic.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"KrovakObliqueConformalConic.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConic.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;AAE7B,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;;;GAkBG;AACH,gBAAgB;AAChB,qBAAa,2BAA4B,SAAQ,eAAe;IAC5D,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,oCAAoC;IACpC,OAAO,CAAC,KAAK,CAAU;IACvB,0BAA0B;IAC1B,OAAO,CAAC,KAAK,CAAU;IACvB,4DAA4D;IAC5D,OAAO,CAAC,KAAK,CAAU;IACvB,2CAA2C;IAC3C,OAAO,CAAC,KAAK,CAAU;IACvB,+CAA+C;IAC/C,OAAO,CAAC,GAAG,CAAU;IACrB,6BAA6B;IAC7B,OAAO,CAAC,GAAG,CAAU;IACrB,8BAA8B;IAC9B,OAAO,CAAC,GAAG,CAAU;IAErB,oBAAoB;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IAErB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAYjD;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IActD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAoB5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAwB/F"}
+{"version":3,"file":"KrovakObliqueConformalConic.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConic.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,KAAK,KAAK,GAAG,MAAM,CAAC;AAIpB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAE7B,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;;;GAkBG;AACH,gBAAgB;AAChB,qBAAa,2BAA4B,SAAQ,eAAe;IAC5D,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,oCAAoC;IACpC,OAAO,CAAC,KAAK,CAAU;IACvB,0BAA0B;IAC1B,OAAO,CAAC,KAAK,CAAU;IACvB,4DAA4D;IAC5D,OAAO,CAAC,KAAK,CAAU;IACvB,2CAA2C;IAC3C,OAAO,CAAC,KAAK,CAAU;IACvB,+CAA+C;IAC/C,OAAO,CAAC,GAAG,CAAU;IACrB,6BAA6B;IAC7B,OAAO,CAAC,GAAG,CAAU;IACrB,8BAA8B;IAC9B,OAAO,CAAC,GAAG,CAAU;IAErB,oBAAoB;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IAErB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAYjD;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IActD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAoB5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAwB/F"}