@itwin/core-orbitgt 4.0.0-dev.52 → 4.0.0-dev.55

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

@@ -1,122 +1,123 @@
-/*---------------------------------------------------------------------------------------------
-* 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 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 class KrovakObliqueConformalConic extends OperationMethod {
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters) {
-        super(KrovakObliqueConformalConic.METHOD_CODE, "Krovak Oblique Conic Conformal", parameters);
-        /* Store the parameters */
-        this._latC = parameters.getValue(8811);
-        this._lonO = parameters.getValue(8833);
-        this._aziC = parameters.getValue(1036); // Parameter 8813 has been changed to 1036 for method 9819 in EPSG version 7.7. updated on 31/08/2017.
-        this._latP = parameters.getValue(8818);
-        this._kP = parameters.getValue(8819);
-        this._FE = parameters.getValue(8806);
-        this._FN = parameters.getValue(8807);
-    }
-    /**
-       * Get the square of a number.
-       */
-    static sqr(v) {
-        return (v * v);
-    }
-    /**
-       * OperationMethod method.
-       * @see OperationMethod#initialize
-       */
-    initialize(operation) {
-        /* Prepare the parameters */
-        const ellipsoid = operation.getSourceCRS().getEllipsoid();
-        const a = ellipsoid.getA();
-        const e = ellipsoid.getE();
-        const e2 = e * e;
-        this._A = a * Math.sqrt(1.0 - e2) / (1.0 - KrovakObliqueConformalConic.sqr(e * Math.sin(this._latC)));
-        this._B = Math.sqrt(1.0 + e2 * Math.pow(Math.cos(this._latC), 4) / (1.0 - e2));
-        this._gO = Math.asin(Math.sin(this._latC) / this._B);
-        this._tO = Math.tan(Math.PI / 4 + this._gO / 2) * Math.pow((1 + e * Math.sin(this._latC)) / (1 - e * Math.sin(this._latC)), e * this._B / 2) / Math.pow(Math.tan(Math.PI / 4 + this._latC / 2), this._B);
-        this._n = Math.sin(this._latP);
-        this._rO = this._kP * this._A / Math.tan(this._latP);
-    }
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#forward
-       */
-    forward(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        const lon = source.getX();
-        const lat = source.getY();
-        /* Calculate easting and northing */
-        const e = sourceCRS.getEllipsoid().getE();
-        const U = 2 * (Math.atan(this._tO * Math.pow(Math.tan(lat / 2 + Math.PI / 4), this._B) / Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e * this._B / 2)) - Math.PI / 4);
-        const V = this._B * (this._lonO - lon);
-        const S = Math.asin(Math.cos(this._aziC) * Math.sin(U) + Math.sin(this._aziC) * Math.cos(U) * Math.cos(V));
-        const D = Math.asin(Math.cos(U) * Math.sin(V) / Math.cos(S));
-        const theta = this._n * D;
-        const r = this._rO * Math.pow(Math.tan(Math.PI / 4 + this._latP / 2), this._n) / Math.pow(Math.tan(S / 2 + Math.PI / 4), this._n);
-        const E = this._FE + r * Math.cos(theta);
-        const N = this._FN + r * Math.sin(theta);
-        /* 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 easting and northing */
-        const E = target.getX();
-        const N = target.getY();
-        /* Get the ellipsoid parameters */
-        const e = sourceCRS.getEllipsoid().getE();
-        /* Calculate longitude and latitude */
-        const r = Math.sqrt(KrovakObliqueConformalConic.sqr(E - this._FE) + KrovakObliqueConformalConic.sqr(N - this._FN));
-        const theta = Math.atan2((N - this._FN), (E - this._FE));
-        const D = theta / Math.sin(this._latP);
-        const S = 2 * (Math.atan(Math.pow(this._rO / r, 1 / this._n) * Math.tan(Math.PI / 4 + this._latP / 2)) - Math.PI / 4);
-        const U = Math.asin(Math.cos(this._aziC) * Math.sin(S) - Math.sin(this._aziC) * Math.cos(S) * Math.cos(D));
-        const V = Math.asin(Math.cos(S) * Math.sin(D) / Math.cos(U));
-        /* Iterate */
-        let lat = U;
-        for (let i = 1; i < 5; i++) {
-            lat = 2 * (Math.atan(Math.pow(this._tO, -1 / this._B) * Math.pow(Math.tan(U / 2 + Math.PI / 4), 1 / this._B) * Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e / 2)) - Math.PI / 4);
-        }
-        const lon = this._lonO - V / this._B;
-        /* Save the position */
-        source.setX(lon);
-        source.setY(lat);
-        source.setZ(target.getZ()); // Keep the Z value
-    }
-}
-/** The code of this method */
-KrovakObliqueConformalConic.METHOD_CODE = 9819;
+/*---------------------------------------------------------------------------------------------
+* 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 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 */
+class KrovakObliqueConformalConic extends OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(KrovakObliqueConformalConic.METHOD_CODE, "Krovak Oblique Conic Conformal", parameters);
+        /* Store the parameters */
+        this._latC = parameters.getValue(8811);
+        this._lonO = parameters.getValue(8833);
+        this._aziC = parameters.getValue(1036); // Parameter 8813 has been changed to 1036 for method 9819 in EPSG version 7.7. updated on 31/08/2017.
+        this._latP = parameters.getValue(8818);
+        this._kP = parameters.getValue(8819);
+        this._FE = parameters.getValue(8806);
+        this._FN = parameters.getValue(8807);
+    }
+    /**
+       * Get the square of a number.
+       */
+    static sqr(v) {
+        return (v * v);
+    }
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation) {
+        /* Prepare the parameters */
+        const ellipsoid = operation.getSourceCRS().getEllipsoid();
+        const a = ellipsoid.getA();
+        const e = ellipsoid.getE();
+        const e2 = e * e;
+        this._A = a * Math.sqrt(1.0 - e2) / (1.0 - KrovakObliqueConformalConic.sqr(e * Math.sin(this._latC)));
+        this._B = Math.sqrt(1.0 + e2 * Math.pow(Math.cos(this._latC), 4) / (1.0 - e2));
+        this._gO = Math.asin(Math.sin(this._latC) / this._B);
+        this._tO = Math.tan(Math.PI / 4 + this._gO / 2) * Math.pow((1 + e * Math.sin(this._latC)) / (1 - e * Math.sin(this._latC)), e * this._B / 2) / Math.pow(Math.tan(Math.PI / 4 + this._latC / 2), this._B);
+        this._n = Math.sin(this._latP);
+        this._rO = this._kP * this._A / Math.tan(this._latP);
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const lon = source.getX();
+        const lat = source.getY();
+        /* Calculate easting and northing */
+        const e = sourceCRS.getEllipsoid().getE();
+        const U = 2 * (Math.atan(this._tO * Math.pow(Math.tan(lat / 2 + Math.PI / 4), this._B) / Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e * this._B / 2)) - Math.PI / 4);
+        const V = this._B * (this._lonO - lon);
+        const S = Math.asin(Math.cos(this._aziC) * Math.sin(U) + Math.sin(this._aziC) * Math.cos(U) * Math.cos(V));
+        const D = Math.asin(Math.cos(U) * Math.sin(V) / Math.cos(S));
+        const theta = this._n * D;
+        const r = this._rO * Math.pow(Math.tan(Math.PI / 4 + this._latP / 2), this._n) / Math.pow(Math.tan(S / 2 + Math.PI / 4), this._n);
+        const E = this._FE + r * Math.cos(theta);
+        const N = this._FN + r * Math.sin(theta);
+        /* 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 easting and northing */
+        const E = target.getX();
+        const N = target.getY();
+        /* Get the ellipsoid parameters */
+        const e = sourceCRS.getEllipsoid().getE();
+        /* Calculate longitude and latitude */
+        const r = Math.sqrt(KrovakObliqueConformalConic.sqr(E - this._FE) + KrovakObliqueConformalConic.sqr(N - this._FN));
+        const theta = Math.atan2((N - this._FN), (E - this._FE));
+        const D = theta / Math.sin(this._latP);
+        const S = 2 * (Math.atan(Math.pow(this._rO / r, 1 / this._n) * Math.tan(Math.PI / 4 + this._latP / 2)) - Math.PI / 4);
+        const U = Math.asin(Math.cos(this._aziC) * Math.sin(S) - Math.sin(this._aziC) * Math.cos(S) * Math.cos(D));
+        const V = Math.asin(Math.cos(S) * Math.sin(D) / Math.cos(U));
+        /* Iterate */
+        let lat = U;
+        for (let i = 1; i < 5; i++) {
+            lat = 2 * (Math.atan(Math.pow(this._tO, -1 / this._B) * Math.pow(Math.tan(U / 2 + Math.PI / 4), 1 / this._B) * Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e / 2)) - Math.PI / 4);
+        }
+        const lon = this._lonO - V / this._B;
+        /* Save the position */
+        source.setX(lon);
+        source.setY(lat);
+        source.setZ(target.getZ()); // Keep the Z value
+    }
+}
+/** The code of this method */
+KrovakObliqueConformalConic.METHOD_CODE = 9819;
+export { KrovakObliqueConformalConic };
 //# sourceMappingURL=KrovakObliqueConformalConic.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"KrovakObliqueConformalConic.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConic.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;;;GAkBG;AACH,gBAAgB;AAChB,MAAM,OAAO,2BAA4B,SAAQ,eAAe;IA2B5D;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,2BAA2B,CAAC,WAAW,EAAE,gCAAgC,EAAE,UAAU,CAAC,CAAC;QAC7F,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,CAAC,sGAAsG;QAC9I,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,GAAG,CAAC,CAAU;QACzB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnB,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,4BAA4B;QAC5B,MAAM,SAAS,GAAc,SAAS,CAAC,YAAY,EAAE,CAAC,YAAY,EAAE,CAAC;QACrE,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,MAAM,EAAE,GAAY,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACtG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;QAC/E,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACrD,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACzM,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACzD,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,oCAAoC;QACpC,MAAM,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACnD,MAAM,CAAC,GAAY,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/L,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;QAChD,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACpH,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC3I,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,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,kCAAkC;QAClC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,MAAM,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACnD,sCAAsC;QACtC,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5H,MAAM,KAAK,GAAY,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAClE,MAAM,CAAC,GAAY,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAChD,MAAM,CAAC,GAAY,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/H,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACpH,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,aAAa;QACb,IAAI,GAAG,GAAY,CAAC,CAAC;QACrB,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;SACrM;QACD,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC9C,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;;AAtHD,8BAA8B;AACP,uCAAW,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 KrovakObliqueConformalConic defines the Krovak Oblique Conformal Conic projection as used in the Czech Republic and Slovakia.\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 (version 6.18.3)\r\n * Revised January 2009\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.3.2 \"Krovak Oblique Conformal Conic\"\r\n *\r\n * NOTE: the X axis runs north to south and the Y axis runs east to west. This differs from most other projections.\r\n *\r\n * @version 1.0 October 2009\r\n */\r\n/** @internal */\r\nexport class KrovakObliqueConformalConic extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9819;\r\n\r\n    /** Latitude of projection center */\r\n    private _latC: float64;\r\n    /** Longitude of origin */\r\n    private _lonO: float64;\r\n    /** Azimuth of initial line through the projection center */\r\n    private _aziC: float64;\r\n    /** Latitude of pseudo standard parallel */\r\n    private _latP: float64;\r\n    /** Scale factor on pseudo standard parallel */\r\n    private _kP: float64;\r\n    /** Easting at grid origin */\r\n    private _FE: float64;\r\n    /** Northing at grid origin */\r\n    private _FN: float64;\r\n\r\n    /** The constants */\r\n    private _A: float64;\r\n    private _B: float64;\r\n    private _gO: float64;\r\n    private _tO: float64;\r\n    private _n: float64;\r\n    private _rO: 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(KrovakObliqueConformalConic.METHOD_CODE, \"Krovak Oblique Conic Conformal\", parameters);\r\n        /* Store the parameters */\r\n        this._latC = parameters.getValue(8811);\r\n        this._lonO = parameters.getValue(8833);\r\n        this._aziC = parameters.getValue(1036); // Parameter 8813 has been changed to 1036 for method 9819 in EPSG version 7.7. updated on 31/08/2017.\r\n        this._latP = parameters.getValue(8818);\r\n        this._kP = parameters.getValue(8819);\r\n        this._FE = parameters.getValue(8806);\r\n        this._FN = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the square of a number.\r\n       */\r\n    private static sqr(v: float64): float64 {\r\n        return (v * v);\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        /* Prepare the parameters */\r\n        const ellipsoid: Ellipsoid = operation.getSourceCRS().getEllipsoid();\r\n        const a: float64 = ellipsoid.getA();\r\n        const e: float64 = ellipsoid.getE();\r\n        const e2: float64 = e * e;\r\n        this._A = a * Math.sqrt(1.0 - e2) / (1.0 - KrovakObliqueConformalConic.sqr(e * Math.sin(this._latC)));\r\n        this._B = Math.sqrt(1.0 + e2 * Math.pow(Math.cos(this._latC), 4) / (1.0 - e2));\r\n        this._gO = Math.asin(Math.sin(this._latC) / this._B);\r\n        this._tO = Math.tan(Math.PI / 4 + this._gO / 2) * Math.pow((1 + e * Math.sin(this._latC)) / (1 - e * Math.sin(this._latC)), e * this._B / 2) / Math.pow(Math.tan(Math.PI / 4 + this._latC / 2), this._B);\r\n        this._n = Math.sin(this._latP);\r\n        this._rO = this._kP * this._A / Math.tan(this._latP);\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        /* Calculate easting and northing */\r\n        const e: float64 = sourceCRS.getEllipsoid().getE();\r\n        const U: float64 = 2 * (Math.atan(this._tO * Math.pow(Math.tan(lat / 2 + Math.PI / 4), this._B) / Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e * this._B / 2)) - Math.PI / 4);\r\n        const V: float64 = this._B * (this._lonO - lon);\r\n        const S: float64 = Math.asin(Math.cos(this._aziC) * Math.sin(U) + Math.sin(this._aziC) * Math.cos(U) * Math.cos(V));\r\n        const D: float64 = Math.asin(Math.cos(U) * Math.sin(V) / Math.cos(S));\r\n        const theta: float64 = this._n * D;\r\n        const r: float64 = this._rO * Math.pow(Math.tan(Math.PI / 4 + this._latP / 2), this._n) / Math.pow(Math.tan(S / 2 + Math.PI / 4), this._n);\r\n        const E: float64 = this._FE + r * Math.cos(theta);\r\n        const N: float64 = this._FN + r * Math.sin(theta);\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 easting and northing */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Get the ellipsoid parameters */\r\n        const e: float64 = sourceCRS.getEllipsoid().getE();\r\n        /* Calculate longitude and latitude */\r\n        const r: float64 = Math.sqrt(KrovakObliqueConformalConic.sqr(E - this._FE) + KrovakObliqueConformalConic.sqr(N - this._FN));\r\n        const theta: float64 = Math.atan2((N - this._FN), (E - this._FE));\r\n        const D: float64 = theta / Math.sin(this._latP);\r\n        const S: float64 = 2 * (Math.atan(Math.pow(this._rO / r, 1 / this._n) * Math.tan(Math.PI / 4 + this._latP / 2)) - Math.PI / 4);\r\n        const U: float64 = Math.asin(Math.cos(this._aziC) * Math.sin(S) - Math.sin(this._aziC) * Math.cos(S) * Math.cos(D));\r\n        const V: float64 = Math.asin(Math.cos(S) * Math.sin(D) / Math.cos(U));\r\n        /* Iterate */\r\n        let lat: float64 = U;\r\n        for (let i: number = 1; i < 5; i++) {\r\n            lat = 2 * (Math.atan(Math.pow(this._tO, -1 / this._B) * Math.pow(Math.tan(U / 2 + Math.PI / 4), 1 / this._B) * Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e / 2)) - Math.PI / 4);\r\n        }\r\n        const lon: float64 = this._lonO - V / this._B;\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":"KrovakObliqueConformalConic.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConic.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;;;GAkBG;AACH,gBAAgB;AAChB,MAAa,2BAA4B,SAAQ,eAAe;IA2B5D;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,2BAA2B,CAAC,WAAW,EAAE,gCAAgC,EAAE,UAAU,CAAC,CAAC;QAC7F,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,CAAC,sGAAsG;QAC9I,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,GAAG,CAAC,CAAU;QACzB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnB,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,4BAA4B;QAC5B,MAAM,SAAS,GAAc,SAAS,CAAC,YAAY,EAAE,CAAC,YAAY,EAAE,CAAC;QACrE,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,MAAM,EAAE,GAAY,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACtG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC;QAC/E,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACrD,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACzM,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC/B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACzD,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,oCAAoC;QACpC,MAAM,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACnD,MAAM,CAAC,GAAY,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/L,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;QAChD,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACpH,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC3I,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,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,kCAAkC;QAClC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,MAAM,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACnD,sCAAsC;QACtC,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,2BAA2B,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5H,MAAM,KAAK,GAAY,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAClE,MAAM,CAAC,GAAY,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAChD,MAAM,CAAC,GAAY,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/H,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACpH,MAAM,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,aAAa;QACb,IAAI,GAAG,GAAY,CAAC,CAAC;QACrB,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;SACrM;QACD,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC9C,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;;AAtHD,8BAA8B;AACP,uCAAW,GAAU,IAAI,CAAC;SAFxC,2BAA2B","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 KrovakObliqueConformalConic defines the Krovak Oblique Conformal Conic projection as used in the Czech Republic and Slovakia.\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 (version 6.18.3)\r\n * Revised January 2009\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.3.2 \"Krovak Oblique Conformal Conic\"\r\n *\r\n * NOTE: the X axis runs north to south and the Y axis runs east to west. This differs from most other projections.\r\n *\r\n * @version 1.0 October 2009\r\n */\r\n/** @internal */\r\nexport class KrovakObliqueConformalConic extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9819;\r\n\r\n    /** Latitude of projection center */\r\n    private _latC: float64;\r\n    /** Longitude of origin */\r\n    private _lonO: float64;\r\n    /** Azimuth of initial line through the projection center */\r\n    private _aziC: float64;\r\n    /** Latitude of pseudo standard parallel */\r\n    private _latP: float64;\r\n    /** Scale factor on pseudo standard parallel */\r\n    private _kP: float64;\r\n    /** Easting at grid origin */\r\n    private _FE: float64;\r\n    /** Northing at grid origin */\r\n    private _FN: float64;\r\n\r\n    /** The constants */\r\n    private _A: float64;\r\n    private _B: float64;\r\n    private _gO: float64;\r\n    private _tO: float64;\r\n    private _n: float64;\r\n    private _rO: 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(KrovakObliqueConformalConic.METHOD_CODE, \"Krovak Oblique Conic Conformal\", parameters);\r\n        /* Store the parameters */\r\n        this._latC = parameters.getValue(8811);\r\n        this._lonO = parameters.getValue(8833);\r\n        this._aziC = parameters.getValue(1036); // Parameter 8813 has been changed to 1036 for method 9819 in EPSG version 7.7. updated on 31/08/2017.\r\n        this._latP = parameters.getValue(8818);\r\n        this._kP = parameters.getValue(8819);\r\n        this._FE = parameters.getValue(8806);\r\n        this._FN = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the square of a number.\r\n       */\r\n    private static sqr(v: float64): float64 {\r\n        return (v * v);\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        /* Prepare the parameters */\r\n        const ellipsoid: Ellipsoid = operation.getSourceCRS().getEllipsoid();\r\n        const a: float64 = ellipsoid.getA();\r\n        const e: float64 = ellipsoid.getE();\r\n        const e2: float64 = e * e;\r\n        this._A = a * Math.sqrt(1.0 - e2) / (1.0 - KrovakObliqueConformalConic.sqr(e * Math.sin(this._latC)));\r\n        this._B = Math.sqrt(1.0 + e2 * Math.pow(Math.cos(this._latC), 4) / (1.0 - e2));\r\n        this._gO = Math.asin(Math.sin(this._latC) / this._B);\r\n        this._tO = Math.tan(Math.PI / 4 + this._gO / 2) * Math.pow((1 + e * Math.sin(this._latC)) / (1 - e * Math.sin(this._latC)), e * this._B / 2) / Math.pow(Math.tan(Math.PI / 4 + this._latC / 2), this._B);\r\n        this._n = Math.sin(this._latP);\r\n        this._rO = this._kP * this._A / Math.tan(this._latP);\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        /* Calculate easting and northing */\r\n        const e: float64 = sourceCRS.getEllipsoid().getE();\r\n        const U: float64 = 2 * (Math.atan(this._tO * Math.pow(Math.tan(lat / 2 + Math.PI / 4), this._B) / Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e * this._B / 2)) - Math.PI / 4);\r\n        const V: float64 = this._B * (this._lonO - lon);\r\n        const S: float64 = Math.asin(Math.cos(this._aziC) * Math.sin(U) + Math.sin(this._aziC) * Math.cos(U) * Math.cos(V));\r\n        const D: float64 = Math.asin(Math.cos(U) * Math.sin(V) / Math.cos(S));\r\n        const theta: float64 = this._n * D;\r\n        const r: float64 = this._rO * Math.pow(Math.tan(Math.PI / 4 + this._latP / 2), this._n) / Math.pow(Math.tan(S / 2 + Math.PI / 4), this._n);\r\n        const E: float64 = this._FE + r * Math.cos(theta);\r\n        const N: float64 = this._FN + r * Math.sin(theta);\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 easting and northing */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Get the ellipsoid parameters */\r\n        const e: float64 = sourceCRS.getEllipsoid().getE();\r\n        /* Calculate longitude and latitude */\r\n        const r: float64 = Math.sqrt(KrovakObliqueConformalConic.sqr(E - this._FE) + KrovakObliqueConformalConic.sqr(N - this._FN));\r\n        const theta: float64 = Math.atan2((N - this._FN), (E - this._FE));\r\n        const D: float64 = theta / Math.sin(this._latP);\r\n        const S: float64 = 2 * (Math.atan(Math.pow(this._rO / r, 1 / this._n) * Math.tan(Math.PI / 4 + this._latP / 2)) - Math.PI / 4);\r\n        const U: float64 = Math.asin(Math.cos(this._aziC) * Math.sin(S) - Math.sin(this._aziC) * Math.cos(S) * Math.cos(D));\r\n        const V: float64 = Math.asin(Math.cos(S) * Math.sin(D) / Math.cos(U));\r\n        /* Iterate */\r\n        let lat: float64 = U;\r\n        for (let i: number = 1; i < 5; i++) {\r\n            lat = 2 * (Math.atan(Math.pow(this._tO, -1 / this._B) * Math.pow(Math.tan(U / 2 + Math.PI / 4), 1 / this._B) * Math.pow((1 + e * Math.sin(lat)) / (1 - e * Math.sin(lat)), e / 2)) - Math.PI / 4);\r\n        }\r\n        const lon: float64 = this._lonO - V / this._B;\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/esm/spatial/ecrs/projection/KrovakObliqueConformalConicEN.d.ts

@@ -1,48 +1,48 @@
-/** @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 KrovakObliqueConformalConicEN defines the "East-North" variant of the Krovak Oblique Conformal Conic projection.
- *
- * In the standard projection X runs from north to south and Y runs from east to west.
- * This seems awkward to most users. They prefer X from west to east and Y from south to north.
- * Therefore an alternative is provided in which both axes are inversed and swapped.
- * This leads to a more traditional map.
- *
- * @version 1.0 October 2009
- */
-/** @internal */
-export declare class KrovakObliqueConformalConicEN extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** The original projection */
-    private projection;
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters: ParameterValueList);
-    /**
-       * 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 KrovakObliqueConformalConicEN defines the "East-North" variant of the Krovak Oblique Conformal Conic projection.
+ *
+ * In the standard projection X runs from north to south and Y runs from east to west.
+ * This seems awkward to most users. They prefer X from west to east and Y from south to north.
+ * Therefore an alternative is provided in which both axes are inversed and swapped.
+ * This leads to a more traditional map.
+ *
+ * @version 1.0 October 2009
+ */
+/** @internal */
+export declare class KrovakObliqueConformalConicEN extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The original projection */
+    private projection;
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters: ParameterValueList);
+    /**
+       * 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=KrovakObliqueConformalConicEN.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"KrovakObliqueConformalConicEN.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConicEN.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;AAG3D;;;;;;;;;GASG;AACH,gBAAgB;AAChB,qBAAa,6BAA8B,SAAQ,eAAe;IAC9D,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,8BAA8B;IAC9B,OAAO,CAAC,UAAU,CAA8B;IAEhD;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAMjD;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAKtD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAW5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAU/F"}
+{"version":3,"file":"KrovakObliqueConformalConicEN.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConicEN.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;AAC7B,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAG3D;;;;;;;;;GASG;AACH,gBAAgB;AAChB,qBAAa,6BAA8B,SAAQ,eAAe;IAC9D,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,8BAA8B;IAC9B,OAAO,CAAC,UAAU,CAA8B;IAEhD;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAMjD;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAKtD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAW5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAU/F"}

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

@@ -1,70 +1,71 @@
-/*---------------------------------------------------------------------------------------------
-* 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";
-import { KrovakObliqueConformalConic } from "./KrovakObliqueConformalConic";
-/**
- * Class KrovakObliqueConformalConicEN defines the "East-North" variant of the Krovak Oblique Conformal Conic projection.
- *
- * In the standard projection X runs from north to south and Y runs from east to west.
- * This seems awkward to most users. They prefer X from west to east and Y from south to north.
- * Therefore an alternative is provided in which both axes are inversed and swapped.
- * This leads to a more traditional map.
- *
- * @version 1.0 October 2009
- */
-/** @internal */
-export class KrovakObliqueConformalConicEN extends OperationMethod {
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters) {
-        super(KrovakObliqueConformalConicEN.METHOD_CODE, "Krovak Oblique Conic Conformal East-North", parameters);
-        /* Create the projection */
-        this.projection = new KrovakObliqueConformalConic(parameters);
-    }
-    /**
-       * OperationMethod method.
-       * @see OperationMethod#initialize
-       */
-    initialize(operation) {
-        /* Prepare the projection */
-        this.projection.initialize(operation);
-    }
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#forward
-       */
-    forward(sourceCRS, source, targetCRS, target) {
-        /* Do the forward projection */
-        this.projection.forward(sourceCRS, source, targetCRS, target);
-        /* Get the original position */
-        const E = target.getX();
-        const N = target.getY();
-        /* Swap */
-        target.setX(-N);
-        target.setY(-E);
-    }
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#reverse
-       */
-    reverse(sourceCRS, source, targetCRS, target) {
-        /* Get the swapped position */
-        const N = -target.getX();
-        const E = -target.getY();
-        /* Swap (leave the target coordinate untouched) */
-        source.setX(E);
-        source.setY(N);
-        /* Do the reverse projection */
-        this.projection.reverse(sourceCRS, source, targetCRS, source /* target*/);
-    }
-}
-/** The code of this method */
-KrovakObliqueConformalConicEN.METHOD_CODE = 1041;
+/*---------------------------------------------------------------------------------------------
+* 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";
+import { KrovakObliqueConformalConic } from "./KrovakObliqueConformalConic";
+/**
+ * Class KrovakObliqueConformalConicEN defines the "East-North" variant of the Krovak Oblique Conformal Conic projection.
+ *
+ * In the standard projection X runs from north to south and Y runs from east to west.
+ * This seems awkward to most users. They prefer X from west to east and Y from south to north.
+ * Therefore an alternative is provided in which both axes are inversed and swapped.
+ * This leads to a more traditional map.
+ *
+ * @version 1.0 October 2009
+ */
+/** @internal */
+class KrovakObliqueConformalConicEN extends OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(KrovakObliqueConformalConicEN.METHOD_CODE, "Krovak Oblique Conic Conformal East-North", parameters);
+        /* Create the projection */
+        this.projection = new KrovakObliqueConformalConic(parameters);
+    }
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation) {
+        /* Prepare the projection */
+        this.projection.initialize(operation);
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Do the forward projection */
+        this.projection.forward(sourceCRS, source, targetCRS, target);
+        /* Get the original position */
+        const E = target.getX();
+        const N = target.getY();
+        /* Swap */
+        target.setX(-N);
+        target.setY(-E);
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the swapped position */
+        const N = -target.getX();
+        const E = -target.getY();
+        /* Swap (leave the target coordinate untouched) */
+        source.setX(E);
+        source.setY(N);
+        /* Do the reverse projection */
+        this.projection.reverse(sourceCRS, source, targetCRS, source /* target*/);
+    }
+}
+/** The code of this method */
+KrovakObliqueConformalConicEN.METHOD_CODE = 1041;
+export { KrovakObliqueConformalConicEN };
 //# sourceMappingURL=KrovakObliqueConformalConicEN.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"KrovakObliqueConformalConicEN.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConicEN.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAaH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAErD,OAAO,EAAE,2BAA2B,EAAE,MAAM,+BAA+B,CAAC;AAE5E;;;;;;;;;GASG;AACH,gBAAgB;AAChB,MAAM,OAAO,6BAA8B,SAAQ,eAAe;IAO9D;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,6BAA6B,CAAC,WAAW,EAAE,2CAA2C,EAAE,UAAU,CAAC,CAAC;QAC1G,2BAA2B;QAC3B,IAAI,CAAC,UAAU,GAAG,IAAI,2BAA2B,CAAC,UAAU,CAAC,CAAC;IAClE,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,4BAA4B;QAC5B,IAAI,CAAC,UAAU,CAAC,UAAU,CAAC,SAAS,CAAC,CAAC;IAC1C,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,+BAA+B;QAC/B,IAAI,CAAC,UAAU,CAAC,OAAO,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,CAAC,CAAC;QAC9D,+BAA+B;QAC/B,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,UAAU;QACV,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,8BAA8B;QAC9B,MAAM,CAAC,GAAY,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;QAClC,MAAM,CAAC,GAAY,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;QAClC,kDAAkD;QAClD,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,+BAA+B;QAC/B,IAAI,CAAC,UAAU,CAAC,OAAO,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,CAAA,WAAW,CAAC,CAAC;IAC7E,CAAC;;AArDD,8BAA8B;AACP,yCAAW,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 { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\nimport { KrovakObliqueConformalConic } from \"./KrovakObliqueConformalConic\";\r\n\r\n/**\r\n * Class KrovakObliqueConformalConicEN defines the \"East-North\" variant of the Krovak Oblique Conformal Conic projection.\r\n *\r\n * In the standard projection X runs from north to south and Y runs from east to west.\r\n * This seems awkward to most users. They prefer X from west to east and Y from south to north.\r\n * Therefore an alternative is provided in which both axes are inversed and swapped.\r\n * This leads to a more traditional map.\r\n *\r\n * @version 1.0 October 2009\r\n */\r\n/** @internal */\r\nexport class KrovakObliqueConformalConicEN extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 1041;\r\n\r\n    /** The original projection */\r\n    private projection: KrovakObliqueConformalConic;\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(KrovakObliqueConformalConicEN.METHOD_CODE, \"Krovak Oblique Conic Conformal East-North\", parameters);\r\n        /* Create the projection */\r\n        this.projection = new KrovakObliqueConformalConic(parameters);\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        /* Prepare the projection */\r\n        this.projection.initialize(operation);\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        /* Do the forward projection */\r\n        this.projection.forward(sourceCRS, source, targetCRS, target);\r\n        /* Get the original position */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Swap */\r\n        target.setX(-N);\r\n        target.setY(-E);\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 swapped position */\r\n        const N: float64 = -target.getX();\r\n        const E: float64 = -target.getY();\r\n        /* Swap (leave the target coordinate untouched) */\r\n        source.setX(E);\r\n        source.setY(N);\r\n        /* Do the reverse projection */\r\n        this.projection.reverse(sourceCRS, source, targetCRS, source/* target*/);\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"KrovakObliqueConformalConicEN.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/KrovakObliqueConformalConicEN.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAaH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAErD,OAAO,EAAE,2BAA2B,EAAE,MAAM,+BAA+B,CAAC;AAE5E;;;;;;;;;GASG;AACH,gBAAgB;AAChB,MAAa,6BAA8B,SAAQ,eAAe;IAO9D;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,6BAA6B,CAAC,WAAW,EAAE,2CAA2C,EAAE,UAAU,CAAC,CAAC;QAC1G,2BAA2B;QAC3B,IAAI,CAAC,UAAU,GAAG,IAAI,2BAA2B,CAAC,UAAU,CAAC,CAAC;IAClE,CAAC;IAED;;;SAGK;IACW,UAAU,CAAC,SAAoB;QAC3C,4BAA4B;QAC5B,IAAI,CAAC,UAAU,CAAC,UAAU,CAAC,SAAS,CAAC,CAAC;IAC1C,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,+BAA+B;QAC/B,IAAI,CAAC,UAAU,CAAC,OAAO,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,CAAC,CAAC;QAC9D,+BAA+B;QAC/B,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,UAAU;QACV,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,8BAA8B;QAC9B,MAAM,CAAC,GAAY,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;QAClC,MAAM,CAAC,GAAY,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC;QAClC,kDAAkD;QAClD,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACf,+BAA+B;QAC/B,IAAI,CAAC,UAAU,CAAC,OAAO,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,CAAA,WAAW,CAAC,CAAC;IAC7E,CAAC;;AArDD,8BAA8B;AACP,yCAAW,GAAU,IAAI,CAAC;SAFxC,6BAA6B","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 { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\nimport { KrovakObliqueConformalConic } from \"./KrovakObliqueConformalConic\";\r\n\r\n/**\r\n * Class KrovakObliqueConformalConicEN defines the \"East-North\" variant of the Krovak Oblique Conformal Conic projection.\r\n *\r\n * In the standard projection X runs from north to south and Y runs from east to west.\r\n * This seems awkward to most users. They prefer X from west to east and Y from south to north.\r\n * Therefore an alternative is provided in which both axes are inversed and swapped.\r\n * This leads to a more traditional map.\r\n *\r\n * @version 1.0 October 2009\r\n */\r\n/** @internal */\r\nexport class KrovakObliqueConformalConicEN extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 1041;\r\n\r\n    /** The original projection */\r\n    private projection: KrovakObliqueConformalConic;\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(KrovakObliqueConformalConicEN.METHOD_CODE, \"Krovak Oblique Conic Conformal East-North\", parameters);\r\n        /* Create the projection */\r\n        this.projection = new KrovakObliqueConformalConic(parameters);\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        /* Prepare the projection */\r\n        this.projection.initialize(operation);\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        /* Do the forward projection */\r\n        this.projection.forward(sourceCRS, source, targetCRS, target);\r\n        /* Get the original position */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Swap */\r\n        target.setX(-N);\r\n        target.setY(-E);\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 swapped position */\r\n        const N: float64 = -target.getX();\r\n        const E: float64 = -target.getY();\r\n        /* Swap (leave the target coordinate untouched) */\r\n        source.setX(E);\r\n        source.setY(N);\r\n        /* Do the reverse projection */\r\n        this.projection.reverse(sourceCRS, source, targetCRS, source/* target*/);\r\n    }\r\n}\r\n"]}