@itwin/core-orbitgt 4.0.0-dev.6 → 4.0.0-dev.61

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

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

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

@@ -1 +1 @@
-{"version":3,"file":"ObliqueStereographic.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/ObliqueStereographic.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,oBAAqB,SAAQ,eAAe;IA6BrD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,oBAAoB,CAAC,WAAW,EAAE,uBAAuB,EAAE,UAAU,CAAC,CAAC;QAC7E,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACpE,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACpE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACtE,CAAC;IAED;;;;;;;SAOK;IACE,MAAM,CAAC,MAAM,CAAC,IAAa,EAAE,IAAa,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QACpF,MAAM,UAAU,GAAyB,IAAI,oBAAoB,CAAC,IAAI,CAAC,CAAC;QACxE,UAAU,CAAC,KAAK,GAAG,IAAI,CAAC;QACxB,UAAU,CAAC,KAAK,GAAG,IAAI,CAAC;QACxB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,OAAO,UAAU,CAAC;IACtB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU,EAAE,CAAU;QACrC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,kCAAkC;QAClC,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,iEAAiE;QACjE,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAY,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9D,MAAM,EAAE,GAAY,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC;QACvD,MAAM,IAAI,GAAY,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAC7E,MAAM,GAAG,GAAY,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAC3D,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC;QACpD,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;QACjJ,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC;QACtD,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC;QAC1E,MAAM,EAAE,GAAY,oBAAoB,CAAC,GAAG,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAClG,MAAM,OAAO,GAAY,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;QACjD,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC;QAC1F,MAAM,EAAE,GAAY,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACjC,IAAI,CAAC,KAAK,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;QAChE,IAAI,CAAC,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC;QAC3B,6BAA6B;QAC7B,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;QAChH,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACpF,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,MAAM,MAAM,GAAY,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtD,gDAAgD;QAChD,MAAM,MAAM,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;QACxE,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpD,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QACxE,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC3G,MAAM,GAAG,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACtE,oCAAoC;QACpC,MAAM,CAAC,GAAY,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxN,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;QAC9I,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC;QAC9P,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,GAAG,IAAI,CAAC,GAAG,CAAC;QAC5C,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC;QAC5C,gDAAgD;QAChD,MAAM,CAAC,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QAChE,MAAM,CAAC,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5I,MAAM,MAAM,GAAY,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;QACpD,MAAM,GAAG,GAAY,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC;QACxE,gEAAgE;QAChE,MAAM,MAAM,GAAY,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtD,MAAM,GAAG,GAAY,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC3G,IAAI,IAAI,GAAY,GAAG,GAAG,oBAAoB,CAAC,IAAI,CAAC,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,oBAAoB,CAAC,EAAE,GAAG,GAAG,CAAC;QACnH,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,MAAM,OAAO,GAAY,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YAClE,MAAM,IAAI,GAAY,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,oBAAoB,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;YAClM,MAAM,QAAQ,GAAY,IAAI,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YAC9H,MAAM,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC;YAClD,IAAI,GAAG,QAAQ,CAAC;YAChB,IAAI,MAAM,GAAG,OAAO;gBAAE,MAAM,CAAC,8BAA8B;SAC9D;QACD,uBAAuB;QACvB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;;AAnND,8BAA8B;AACP,gCAAW,GAAU,IAAI,CAAC;AAEjD,sBAAsB;AACE,uBAAE,GAAY,IAAI,CAAC,EAAE,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 ObliqueStereographic defines an Oblique Stereographic projection.\r\n *\r\n * The 'source' CRS is the geodetic 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.7.1 'Oblique and Equatorial Stereographic cases'\r\n *\r\n * @version 1.0 December 2006\r\n */\r\n/** @internal */\r\nexport class ObliqueStereographic extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9809;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n\r\n    /** Latitude of false origin */\r\n    private _lat0: float64;\r\n    /** Longitude of false origin */\r\n    private _lon0: float64;\r\n    /** Scale factor at natural origin */\r\n    private _k0: float64;\r\n    /** False easting */\r\n    private _fe: float64;\r\n    /** False northing */\r\n    private _fn: float64;\r\n\r\n    /** The parameters for the conformal sphere */\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _R: float64;\r\n    private _n: float64;\r\n    private _c: float64;\r\n    private _chi0: float64;\r\n    private _lambda0: float64;\r\n    private _g: float64;\r\n    private _h: 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(ObliqueStereographic.METHOD_CODE, \"Oblique Stereographic\", parameters);\r\n        /* Store the parameters */\r\n        this._lat0 = (parameters != null) ? parameters.getValue(8801) : 0.0;\r\n        this._lon0 = (parameters != null) ? parameters.getValue(8802) : 0.0;\r\n        this._k0 = (parameters != null) ? parameters.getValue(8805) : 0.0;\r\n        this._fe = (parameters != null) ? parameters.getValue(8806) : 0.0;\r\n        this._fn = (parameters != null) ? parameters.getValue(8807) : 0.0;\r\n    }\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param lat0 latitude of natural origin (radians).\r\n       * @param lon0 longitude of natural origin (radians).\r\n       * @param k0 scale factor at natural origin.\r\n       * @param fe false easting.\r\n       * @param fn false northing.\r\n       */\r\n    public static create(lat0: float64, lon0: float64, k0: float64, fe: float64, fn: float64): ObliqueStereographic {\r\n        const projection: ObliqueStereographic = new ObliqueStereographic(null);\r\n        projection._lat0 = lat0;\r\n        projection._lon0 = lon0;\r\n        projection._k0 = k0;\r\n        projection._fe = fe;\r\n        projection._fn = fn;\r\n        return projection;\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number.\r\n       */\r\n    private static pow(n: float64, e: float64): float64 {\r\n        return Math.exp(e * Math.log(n));\r\n    }\r\n\r\n    /**\r\n       * Calculate sinus.\r\n       */\r\n    private static sin(a: float64): float64 {\r\n        return Math.sin(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate cosinus.\r\n       */\r\n    private static cos(a: float64): float64 {\r\n        return Math.cos(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate tangent.\r\n       */\r\n    private static tan(a: float64): float64 {\r\n        return Math.tan(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate arcsinus.\r\n       */\r\n    private static asin(v: float64): float64 {\r\n        return Math.asin(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate arctangent.\r\n       */\r\n    private static atan(v: float64): float64 {\r\n        return Math.atan(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate log.\r\n       */\r\n    private static log(v: float64): float64 {\r\n        return Math.log(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate exp.\r\n       */\r\n    private static exp(v: float64): float64 {\r\n        return Math.exp(v);\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): ObliqueStereographic {\r\n        /* Get the ellipsoid parameters */\r\n        const a: float64 = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        /* Calculate conformat latitude (chi0) and longitude (lambda0) */\r\n        this._e2 = this._e * this._e;\r\n        const sinLat0: float64 = ObliqueStereographic.sin(this._lat0);\r\n        const t0: float64 = 1.0 - this._e2 * sinLat0 * sinLat0;\r\n        const rho0: float64 = a * (1 - this._e2) / ObliqueStereographic.pow(t0, 1.5);\r\n        const nu0: float64 = a / ObliqueStereographic.pow(t0, 0.5);\r\n        this._R = ObliqueStereographic.pow(rho0 * nu0, 0.5);\r\n        this._n = ObliqueStereographic.pow(1.0 + (this._e2 * ObliqueStereographic.pow(ObliqueStereographic.cos(this._lat0), 4)) / (1.0 - this._e2), 0.5);\r\n        const S1: float64 = (1.0 + sinLat0) / (1.0 - sinLat0);\r\n        const S2: float64 = (1.0 - this._e * sinLat0) / (1.0 + this._e * sinLat0);\r\n        const w1: float64 = ObliqueStereographic.pow(S1 * ObliqueStereographic.pow(S2, this._e), this._n);\r\n        const sinChi0: float64 = (w1 - 1.0) / (w1 + 1.0);\r\n        this._c = (this._n + sinLat0) * (1.0 - sinChi0) / ((this._n - sinLat0) * (1.0 + sinChi0));\r\n        const w2: float64 = this._c * w1;\r\n        this._chi0 = ObliqueStereographic.asin((w2 - 1.0) / (w2 + 1.0));\r\n        this._lambda0 = this._lon0;\r\n        /* For reverse calculation */\r\n        this._g = 2.0 * this._R * this._k0 * ObliqueStereographic.tan(ObliqueStereographic.PI / 4.0 - this._chi0 / 2.0);\r\n        this._h = 4.0 * this._R * this._k0 * ObliqueStereographic.tan(this._chi0) + this._g;\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        const sinLat: float64 = ObliqueStereographic.sin(lat);\r\n        /* Calculate conformal latitude and longitude */\r\n        const lambda: float64 = this._n * (lon - this._lambda0) + this._lambda0;\r\n        const Sa: float64 = (1.0 + sinLat) / (1.0 - sinLat);\r\n        const Sb: float64 = (1.0 - this._e * sinLat) / (1.0 + this._e * sinLat);\r\n        const w: float64 = this._c * ObliqueStereographic.pow(Sa * ObliqueStereographic.pow(Sb, this._e), this._n);\r\n        const chi: float64 = ObliqueStereographic.asin((w - 1.0) / (w + 1.0));\r\n        /* Calculate easting and northing */\r\n        const B: float64 = 1.0 + ObliqueStereographic.sin(chi) * ObliqueStereographic.sin(this._chi0) + ObliqueStereographic.cos(chi) * ObliqueStereographic.cos(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0);\r\n        const E: float64 = this._fe + 2.0 * this._R * this._k0 * ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(lambda - this._lambda0) / B;\r\n        const N: float64 = this._fn + 2.0 * this._R * this._k0 * (ObliqueStereographic.sin(chi) * ObliqueStereographic.cos(this._chi0) - ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0)) / B;\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() - this._fe;\r\n        const N: float64 = target.getY() - this._fn;\r\n        /* Calculate conformal latitude and longitude */\r\n        const i: float64 = ObliqueStereographic.atan(E / (this._h + N));\r\n        const j: float64 = ObliqueStereographic.atan(E / (this._g - N)) - i;\r\n        const chi: float64 = this._chi0 + 2.0 * ObliqueStereographic.atan((N - E * ObliqueStereographic.tan(j / 2.0)) / (2.0 * this._R * this._k0));\r\n        const lambda: float64 = j + 2.0 * i + this._lambda0;\r\n        const lon: float64 = (lambda - this._lambda0) / this._n + this._lambda0;\r\n        /* Iterate to get isometric (psi) and geodetic (phi) latitude */\r\n        const sinChi: float64 = ObliqueStereographic.sin(chi);\r\n        const psi: float64 = 0.5 * ObliqueStereographic.log((1.0 + sinChi) / (this._c * (1.0 - sinChi))) / this._n;\r\n        let phiK: float64 = 2.0 * ObliqueStereographic.atan(ObliqueStereographic.exp(psi)) - ObliqueStereographic.PI / 2.0;\r\n        for (let k: number = 1; k < 8; k++) {\r\n            const esinPhi: float64 = this._e * ObliqueStereographic.sin(phiK);\r\n            const psiK: float64 = ObliqueStereographic.log(ObliqueStereographic.tan(phiK / 2.0 + ObliqueStereographic.PI / 4.0) * ObliqueStereographic.pow((1.0 - esinPhi) / (1.0 + esinPhi), this._e / 2.0));\r\n            const phiNextK: float64 = phiK - (psiK - psi) * ObliqueStereographic.cos(phiK) * (1.0 - esinPhi * esinPhi) / (1.0 - this._e2);\r\n            const change: float64 = Math.abs(phiNextK - phiK);\r\n            phiK = phiNextK;\r\n            if (change < 1.0e-12) break; // normally after 4 iterations\r\n        }\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(phiK);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"ObliqueStereographic.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/ObliqueStereographic.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAa,oBAAqB,SAAQ,eAAe;IA6BrD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,oBAAoB,CAAC,WAAW,EAAE,uBAAuB,EAAE,UAAU,CAAC,CAAC;QAC7E,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACpE,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACpE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAClE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACtE,CAAC;IAED;;;;;;;SAOK;IACE,MAAM,CAAC,MAAM,CAAC,IAAa,EAAE,IAAa,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QACpF,MAAM,UAAU,GAAyB,IAAI,oBAAoB,CAAC,IAAI,CAAC,CAAC;QACxE,UAAU,CAAC,KAAK,GAAG,IAAI,CAAC;QACxB,UAAU,CAAC,KAAK,GAAG,IAAI,CAAC;QACxB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,UAAU,CAAC,GAAG,GAAG,EAAE,CAAC;QACpB,OAAO,UAAU,CAAC;IACtB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU,EAAE,CAAU;QACrC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,kCAAkC;QAClC,MAAM,CAAC,GAAY,SAAS,CAAC,IAAI,EAAE,CAAC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,iEAAiE;QACjE,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAY,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC9D,MAAM,EAAE,GAAY,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC;QACvD,MAAM,IAAI,GAAY,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAC7E,MAAM,GAAG,GAAY,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QAC3D,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC;QACpD,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;QACjJ,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC;QACtD,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC;QAC1E,MAAM,EAAE,GAAY,oBAAoB,CAAC,GAAG,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAClG,MAAM,OAAO,GAAY,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;QACjD,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC,CAAC;QAC1F,MAAM,EAAE,GAAY,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACjC,IAAI,CAAC,KAAK,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;QAChE,IAAI,CAAC,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC;QAC3B,6BAA6B;QAC7B,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;QAChH,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACpF,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,MAAM,MAAM,GAAY,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtD,gDAAgD;QAChD,MAAM,MAAM,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;QACxE,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpD,MAAM,EAAE,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QACxE,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QAC3G,MAAM,GAAG,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACtE,oCAAoC;QACpC,MAAM,CAAC,GAAY,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxN,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;QAC9I,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC;QAC9P,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,GAAG,IAAI,CAAC,GAAG,CAAC;QAC5C,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC;QAC5C,gDAAgD;QAChD,MAAM,CAAC,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QAChE,MAAM,CAAC,GAAY,oBAAoB,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,GAAG,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5I,MAAM,MAAM,GAAY,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC;QACpD,MAAM,GAAG,GAAY,CAAC,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC;QACxE,gEAAgE;QAChE,MAAM,MAAM,GAAY,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACtD,MAAM,GAAG,GAAY,GAAG,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC3G,IAAI,IAAI,GAAY,GAAG,GAAG,oBAAoB,CAAC,IAAI,CAAC,oBAAoB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,oBAAoB,CAAC,EAAE,GAAG,GAAG,CAAC;QACnH,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,MAAM,OAAO,GAAY,IAAI,CAAC,EAAE,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YAClE,MAAM,IAAI,GAAY,oBAAoB,CAAC,GAAG,CAAC,oBAAoB,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,oBAAoB,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;YAClM,MAAM,QAAQ,GAAY,IAAI,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,oBAAoB,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YAC9H,MAAM,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,CAAC;YAClD,IAAI,GAAG,QAAQ,CAAC;YAChB,IAAI,MAAM,GAAG,OAAO;gBAAE,MAAM,CAAC,8BAA8B;SAC9D;QACD,uBAAuB;QACvB,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;;AAnND,8BAA8B;AACP,gCAAW,GAAU,IAAI,CAAC;AAEjD,sBAAsB;AACE,uBAAE,GAAY,IAAI,CAAC,EAAE,CAAC;SALrC,oBAAoB","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 ObliqueStereographic defines an Oblique Stereographic projection.\r\n *\r\n * The 'source' CRS is the geodetic 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.7.1 'Oblique and Equatorial Stereographic cases'\r\n *\r\n * @version 1.0 December 2006\r\n */\r\n/** @internal */\r\nexport class ObliqueStereographic extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9809;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n\r\n    /** Latitude of false origin */\r\n    private _lat0: float64;\r\n    /** Longitude of false origin */\r\n    private _lon0: float64;\r\n    /** Scale factor at natural origin */\r\n    private _k0: float64;\r\n    /** False easting */\r\n    private _fe: float64;\r\n    /** False northing */\r\n    private _fn: float64;\r\n\r\n    /** The parameters for the conformal sphere */\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _R: float64;\r\n    private _n: float64;\r\n    private _c: float64;\r\n    private _chi0: float64;\r\n    private _lambda0: float64;\r\n    private _g: float64;\r\n    private _h: 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(ObliqueStereographic.METHOD_CODE, \"Oblique Stereographic\", parameters);\r\n        /* Store the parameters */\r\n        this._lat0 = (parameters != null) ? parameters.getValue(8801) : 0.0;\r\n        this._lon0 = (parameters != null) ? parameters.getValue(8802) : 0.0;\r\n        this._k0 = (parameters != null) ? parameters.getValue(8805) : 0.0;\r\n        this._fe = (parameters != null) ? parameters.getValue(8806) : 0.0;\r\n        this._fn = (parameters != null) ? parameters.getValue(8807) : 0.0;\r\n    }\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param lat0 latitude of natural origin (radians).\r\n       * @param lon0 longitude of natural origin (radians).\r\n       * @param k0 scale factor at natural origin.\r\n       * @param fe false easting.\r\n       * @param fn false northing.\r\n       */\r\n    public static create(lat0: float64, lon0: float64, k0: float64, fe: float64, fn: float64): ObliqueStereographic {\r\n        const projection: ObliqueStereographic = new ObliqueStereographic(null);\r\n        projection._lat0 = lat0;\r\n        projection._lon0 = lon0;\r\n        projection._k0 = k0;\r\n        projection._fe = fe;\r\n        projection._fn = fn;\r\n        return projection;\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number.\r\n       */\r\n    private static pow(n: float64, e: float64): float64 {\r\n        return Math.exp(e * Math.log(n));\r\n    }\r\n\r\n    /**\r\n       * Calculate sinus.\r\n       */\r\n    private static sin(a: float64): float64 {\r\n        return Math.sin(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate cosinus.\r\n       */\r\n    private static cos(a: float64): float64 {\r\n        return Math.cos(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate tangent.\r\n       */\r\n    private static tan(a: float64): float64 {\r\n        return Math.tan(a);\r\n    }\r\n\r\n    /**\r\n       * Calculate arcsinus.\r\n       */\r\n    private static asin(v: float64): float64 {\r\n        return Math.asin(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate arctangent.\r\n       */\r\n    private static atan(v: float64): float64 {\r\n        return Math.atan(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate log.\r\n       */\r\n    private static log(v: float64): float64 {\r\n        return Math.log(v);\r\n    }\r\n\r\n    /**\r\n       * Calculate exp.\r\n       */\r\n    private static exp(v: float64): float64 {\r\n        return Math.exp(v);\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): ObliqueStereographic {\r\n        /* Get the ellipsoid parameters */\r\n        const a: float64 = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        /* Calculate conformat latitude (chi0) and longitude (lambda0) */\r\n        this._e2 = this._e * this._e;\r\n        const sinLat0: float64 = ObliqueStereographic.sin(this._lat0);\r\n        const t0: float64 = 1.0 - this._e2 * sinLat0 * sinLat0;\r\n        const rho0: float64 = a * (1 - this._e2) / ObliqueStereographic.pow(t0, 1.5);\r\n        const nu0: float64 = a / ObliqueStereographic.pow(t0, 0.5);\r\n        this._R = ObliqueStereographic.pow(rho0 * nu0, 0.5);\r\n        this._n = ObliqueStereographic.pow(1.0 + (this._e2 * ObliqueStereographic.pow(ObliqueStereographic.cos(this._lat0), 4)) / (1.0 - this._e2), 0.5);\r\n        const S1: float64 = (1.0 + sinLat0) / (1.0 - sinLat0);\r\n        const S2: float64 = (1.0 - this._e * sinLat0) / (1.0 + this._e * sinLat0);\r\n        const w1: float64 = ObliqueStereographic.pow(S1 * ObliqueStereographic.pow(S2, this._e), this._n);\r\n        const sinChi0: float64 = (w1 - 1.0) / (w1 + 1.0);\r\n        this._c = (this._n + sinLat0) * (1.0 - sinChi0) / ((this._n - sinLat0) * (1.0 + sinChi0));\r\n        const w2: float64 = this._c * w1;\r\n        this._chi0 = ObliqueStereographic.asin((w2 - 1.0) / (w2 + 1.0));\r\n        this._lambda0 = this._lon0;\r\n        /* For reverse calculation */\r\n        this._g = 2.0 * this._R * this._k0 * ObliqueStereographic.tan(ObliqueStereographic.PI / 4.0 - this._chi0 / 2.0);\r\n        this._h = 4.0 * this._R * this._k0 * ObliqueStereographic.tan(this._chi0) + this._g;\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        const sinLat: float64 = ObliqueStereographic.sin(lat);\r\n        /* Calculate conformal latitude and longitude */\r\n        const lambda: float64 = this._n * (lon - this._lambda0) + this._lambda0;\r\n        const Sa: float64 = (1.0 + sinLat) / (1.0 - sinLat);\r\n        const Sb: float64 = (1.0 - this._e * sinLat) / (1.0 + this._e * sinLat);\r\n        const w: float64 = this._c * ObliqueStereographic.pow(Sa * ObliqueStereographic.pow(Sb, this._e), this._n);\r\n        const chi: float64 = ObliqueStereographic.asin((w - 1.0) / (w + 1.0));\r\n        /* Calculate easting and northing */\r\n        const B: float64 = 1.0 + ObliqueStereographic.sin(chi) * ObliqueStereographic.sin(this._chi0) + ObliqueStereographic.cos(chi) * ObliqueStereographic.cos(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0);\r\n        const E: float64 = this._fe + 2.0 * this._R * this._k0 * ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(lambda - this._lambda0) / B;\r\n        const N: float64 = this._fn + 2.0 * this._R * this._k0 * (ObliqueStereographic.sin(chi) * ObliqueStereographic.cos(this._chi0) - ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(this._chi0) * ObliqueStereographic.cos(lambda - this._lambda0)) / B;\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() - this._fe;\r\n        const N: float64 = target.getY() - this._fn;\r\n        /* Calculate conformal latitude and longitude */\r\n        const i: float64 = ObliqueStereographic.atan(E / (this._h + N));\r\n        const j: float64 = ObliqueStereographic.atan(E / (this._g - N)) - i;\r\n        const chi: float64 = this._chi0 + 2.0 * ObliqueStereographic.atan((N - E * ObliqueStereographic.tan(j / 2.0)) / (2.0 * this._R * this._k0));\r\n        const lambda: float64 = j + 2.0 * i + this._lambda0;\r\n        const lon: float64 = (lambda - this._lambda0) / this._n + this._lambda0;\r\n        /* Iterate to get isometric (psi) and geodetic (phi) latitude */\r\n        const sinChi: float64 = ObliqueStereographic.sin(chi);\r\n        const psi: float64 = 0.5 * ObliqueStereographic.log((1.0 + sinChi) / (this._c * (1.0 - sinChi))) / this._n;\r\n        let phiK: float64 = 2.0 * ObliqueStereographic.atan(ObliqueStereographic.exp(psi)) - ObliqueStereographic.PI / 2.0;\r\n        for (let k: number = 1; k < 8; k++) {\r\n            const esinPhi: float64 = this._e * ObliqueStereographic.sin(phiK);\r\n            const psiK: float64 = ObliqueStereographic.log(ObliqueStereographic.tan(phiK / 2.0 + ObliqueStereographic.PI / 4.0) * ObliqueStereographic.pow((1.0 - esinPhi) / (1.0 + esinPhi), this._e / 2.0));\r\n            const phiNextK: float64 = phiK - (psiK - psi) * ObliqueStereographic.cos(phiK) * (1.0 - esinPhi * esinPhi) / (1.0 - this._e2);\r\n            const change: float64 = Math.abs(phiNextK - phiK);\r\n            phiK = phiNextK;\r\n            if (change < 1.0e-12) break; // normally after 4 iterations\r\n        }\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(phiK);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}