@itwin/core-orbitgt 5.0.0-dev.8 → 5.0.0-dev.83

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

@@ -25,6 +25,18 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class LambertConical1SP extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9801;
+    /** Latitude of natural origin */
+    _latN;
+    /** Longitude of natural origin (the central meridian) */
+    _lonN;
+    /** Scale factor at natural origin */
+    _k0;
+    /** False easting */
+    _fE;
+    /** False northing */
+    _fN;
     /**
      * Create a new projection.
      * @param parameters the values of the parameters.
@@ -121,6 +133,4 @@ export class LambertConical1SP extends OperationMethod {
         source.setZ(target.getZ()); // Keep the Z value
     }
 }
-/** The code of this method */
-LambertConical1SP.METHOD_CODE = 9801;
 //# sourceMappingURL=LambertConical1SP.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"LambertConical1SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical1SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAYH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,iBAAkB,SAAQ,eAAe;IAelD;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,iBAAiB,CAAC,WAAW,EAAE,+BAA+B,EAAE,UAAU,CAAC,CAAC;QAClF,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,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;;OAEG;IACK,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;IACjF,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;IAC3H,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,EAAW;QAC5E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM;YAAE,OAAO,GAAG,CAAC;QACrC,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IACvC,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,0BAA0B;QAC1B,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACtC,IAAI,CAAC,GAAY,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACjE,IAAI,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACjD,IAAI,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QAC/D,IAAI,KAAK,GAAY,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QAC5C,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAChD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACrD,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,0BAA0B;QAC1B,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACtC,IAAI,CAAC,GAAY,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACjE,IAAI,EAAE,GAAY,iBAAiB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;QAC9H,IAAI,EAAE,GAAY,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;QAC7D,IAAI,MAAM,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACxE,IAAI,GAAG,GAAY,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvD,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,oBAAoB;QAC9K,IAAI,GAAG,GAAY,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QAC7C,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;;AAhHD,8BAA8B;AACP,6BAAW,GAAU,IAAI,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n//package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised August 2006\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.1.1 -> 1.4.1.4\n *\n * @version 1.0 January 2007\n */\n/** @internal */\nexport class LambertConical1SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9801;\n\n    /** Latitude of natural origin */\n    private _latN: float64;\n    /** Longitude of natural origin (the central meridian) */\n    private _lonN: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fE: float64;\n    /** False northing */\n    private _fN: float64;\n\n    /**\n     * Create a new projection.\n     * @param parameters the values of the parameters.\n     */\n    public constructor(parameters: ParameterValueList) {\n        super(LambertConical1SP.METHOD_CODE, \"Lambert Conic Conformal (1SP)\", parameters);\n        /* Store the parameters */\n        this._latN = parameters.getValue(8801);\n        this._lonN = parameters.getValue(8802);\n        this._k0 = parameters.getValue(8805);\n        this._fE = parameters.getValue(8806);\n        this._fN = parameters.getValue(8807);\n    }\n\n    /**\n     * Get the sign of a number.\n     */\n    private static sign(v: float64): float64 {\n        return (v < 0.0) ? (-1.0) : (1.0);\n    }\n\n    /**\n     * Calculate M.\n     */\n    private static calcM(e: float64, lat: float64): float64 {\n        return Math.cos(lat) / Math.pow(1.0 - Math.pow(e * Math.sin(lat), 2.0), 0.5);\n    }\n\n    /**\n     * Calculate T.\n     */\n    private static calcT(e: float64, lat: float64): float64 {\n        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);\n    }\n\n    /**\n     * Calculate R.\n     */\n    private static calcR(a: float64, F: float64, t: float64, n: float64, k0: float64): float64 {\n        if (Math.abs(t) < 1.0e-6) return 0.0;\n        return a * F * Math.pow(t, n) * k0;\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#forward\n     */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let lon: float64 = source.getX();\n        let lat: float64 = source.getY();\n        /* Get the ellipsoid parameters */\n        let e: float64 = sourceCRS.getEllipsoid().getE();\n        let a: float64 = sourceCRS.getEllipsoid().getA();\n        /* Make the calculation */\n        let m0: float64 = LambertConical1SP.calcM(e, this._latN);\n        let t0: float64 = LambertConical1SP.calcT(e, this._latN);\n        let n: float64 = Math.sin(this._latN);\n        let F: float64 = m0 / (n * Math.pow(t0, n));\n        let r0: float64 = LambertConical1SP.calcR(a, F, t0, n, this._k0);\n        let t: float64 = LambertConical1SP.calcT(e, lat);\n        let r: float64 = LambertConical1SP.calcR(a, F, t, n, this._k0);\n        let theta: float64 = n * (lon - this._lonN);\n        let E: float64 = this._fE + r * Math.sin(theta);\n        let N: float64 = this._fN + r0 - r * Math.cos(theta);\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#reverse\n     */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let E: float64 = target.getX();\n        let N: float64 = target.getY();\n        /* Get the ellipsoid parameters */\n        let e: float64 = sourceCRS.getEllipsoid().getE();\n        let a: float64 = sourceCRS.getEllipsoid().getA();\n        /* Make the calculation */\n        let m0: float64 = LambertConical1SP.calcM(e, this._latN);\n        let t0: float64 = LambertConical1SP.calcT(e, this._latN);\n        let n: float64 = Math.sin(this._latN);\n        let F: float64 = m0 / (n * Math.pow(t0, n));\n        let r0: float64 = LambertConical1SP.calcR(a, F, t0, n, this._k0);\n        let r_: float64 = LambertConical1SP.sign(n) * Math.pow(Math.pow(E - this._fE, 2.0) + Math.pow(r0 - (N - this._fN), 2.0), 0.5);\n        let t_: float64 = Math.pow(r_ / (a * this._k0 * F), 1.0 / n);\n        let theta_: float64 = Math.atan((E - this._fE) / (r0 - (N - this._fN)));\n        let lat: float64 = Math.PI * 0.5 - 2.0 * Math.atan(t_);\n        for (let i: number = 0; i < 7; i++) lat = Math.PI * 0.5 - 2.0 * Math.atan(t_ * Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5)); // recursive formula\n        let lon: float64 = (theta_) / n + this._lonN;\n        /* Save the position */\n        source.setX(lon);\n        source.setY(lat);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}
+{"version":3,"file":"LambertConical1SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical1SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAYH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,iBAAkB,SAAQ,eAAe;IAClD,8BAA8B;IACvB,MAAM,CAAU,WAAW,GAAU,IAAI,CAAC;IAEjD,iCAAiC;IACzB,KAAK,CAAU;IACvB,yDAAyD;IACjD,KAAK,CAAU;IACvB,qCAAqC;IAC7B,GAAG,CAAU;IACrB,oBAAoB;IACZ,GAAG,CAAU;IACrB,qBAAqB;IACb,GAAG,CAAU;IAErB;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,iBAAiB,CAAC,WAAW,EAAE,+BAA+B,EAAE,UAAU,CAAC,CAAC;QAClF,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,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;;OAEG;IACK,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;IACjF,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;IAC3H,CAAC;IAED;;OAEG;IACK,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,EAAW;QAC5E,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM;YAAE,OAAO,GAAG,CAAC;QACrC,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;IACvC,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,0BAA0B;QAC1B,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACtC,IAAI,CAAC,GAAY,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACjE,IAAI,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACjD,IAAI,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QAC/D,IAAI,KAAK,GAAY,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QAC5C,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAChD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACrD,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,0BAA0B;QAC1B,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACtC,IAAI,CAAC,GAAY,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAC5C,IAAI,EAAE,GAAY,iBAAiB,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACjE,IAAI,EAAE,GAAY,iBAAiB,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;QAC9H,IAAI,EAAE,GAAY,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;QAC7D,IAAI,MAAM,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACxE,IAAI,GAAG,GAAY,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvD,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,oBAAoB;QAC9K,IAAI,GAAG,GAAY,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QAC7C,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","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n//package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised August 2006\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.1.1 -> 1.4.1.4\n *\n * @version 1.0 January 2007\n */\n/** @internal */\nexport class LambertConical1SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9801;\n\n    /** Latitude of natural origin */\n    private _latN: float64;\n    /** Longitude of natural origin (the central meridian) */\n    private _lonN: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fE: float64;\n    /** False northing */\n    private _fN: float64;\n\n    /**\n     * Create a new projection.\n     * @param parameters the values of the parameters.\n     */\n    public constructor(parameters: ParameterValueList) {\n        super(LambertConical1SP.METHOD_CODE, \"Lambert Conic Conformal (1SP)\", parameters);\n        /* Store the parameters */\n        this._latN = parameters.getValue(8801);\n        this._lonN = parameters.getValue(8802);\n        this._k0 = parameters.getValue(8805);\n        this._fE = parameters.getValue(8806);\n        this._fN = parameters.getValue(8807);\n    }\n\n    /**\n     * Get the sign of a number.\n     */\n    private static sign(v: float64): float64 {\n        return (v < 0.0) ? (-1.0) : (1.0);\n    }\n\n    /**\n     * Calculate M.\n     */\n    private static calcM(e: float64, lat: float64): float64 {\n        return Math.cos(lat) / Math.pow(1.0 - Math.pow(e * Math.sin(lat), 2.0), 0.5);\n    }\n\n    /**\n     * Calculate T.\n     */\n    private static calcT(e: float64, lat: float64): float64 {\n        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);\n    }\n\n    /**\n     * Calculate R.\n     */\n    private static calcR(a: float64, F: float64, t: float64, n: float64, k0: float64): float64 {\n        if (Math.abs(t) < 1.0e-6) return 0.0;\n        return a * F * Math.pow(t, n) * k0;\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#forward\n     */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let lon: float64 = source.getX();\n        let lat: float64 = source.getY();\n        /* Get the ellipsoid parameters */\n        let e: float64 = sourceCRS.getEllipsoid().getE();\n        let a: float64 = sourceCRS.getEllipsoid().getA();\n        /* Make the calculation */\n        let m0: float64 = LambertConical1SP.calcM(e, this._latN);\n        let t0: float64 = LambertConical1SP.calcT(e, this._latN);\n        let n: float64 = Math.sin(this._latN);\n        let F: float64 = m0 / (n * Math.pow(t0, n));\n        let r0: float64 = LambertConical1SP.calcR(a, F, t0, n, this._k0);\n        let t: float64 = LambertConical1SP.calcT(e, lat);\n        let r: float64 = LambertConical1SP.calcR(a, F, t, n, this._k0);\n        let theta: float64 = n * (lon - this._lonN);\n        let E: float64 = this._fE + r * Math.sin(theta);\n        let N: float64 = this._fN + r0 - r * Math.cos(theta);\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#reverse\n     */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let E: float64 = target.getX();\n        let N: float64 = target.getY();\n        /* Get the ellipsoid parameters */\n        let e: float64 = sourceCRS.getEllipsoid().getE();\n        let a: float64 = sourceCRS.getEllipsoid().getA();\n        /* Make the calculation */\n        let m0: float64 = LambertConical1SP.calcM(e, this._latN);\n        let t0: float64 = LambertConical1SP.calcT(e, this._latN);\n        let n: float64 = Math.sin(this._latN);\n        let F: float64 = m0 / (n * Math.pow(t0, n));\n        let r0: float64 = LambertConical1SP.calcR(a, F, t0, n, this._k0);\n        let r_: float64 = LambertConical1SP.sign(n) * Math.pow(Math.pow(E - this._fE, 2.0) + Math.pow(r0 - (N - this._fN), 2.0), 0.5);\n        let t_: float64 = Math.pow(r_ / (a * this._k0 * F), 1.0 / n);\n        let theta_: float64 = Math.atan((E - this._fE) / (r0 - (N - this._fN)));\n        let lat: float64 = Math.PI * 0.5 - 2.0 * Math.atan(t_);\n        for (let i: number = 0; i < 7; i++) lat = Math.PI * 0.5 - 2.0 * Math.atan(t_ * Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5)); // recursive formula\n        let lon: float64 = (theta_) / n + this._lonN;\n        /* Save the position */\n        source.setX(lon);\n        source.setY(lat);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}

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

@@ -25,6 +25,35 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class LambertConical2SP extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9802;
+    /** The value of PI */
+    static PI = Math.PI;
+    /** The half value of PI */
+    static hPI = 0.5 * Math.PI;
+    /** Latitude of false origin */
+    _latF;
+    /** Longitude of false origin (the central meridian) */
+    _lonF;
+    /** Latitude of first standard parallel */
+    _lat1;
+    /** Latitude of second standard parallel */
+    _lat2;
+    /** Easting at false origin */
+    _eF;
+    /** Northing at false origin */
+    _nF;
+    /* Calculation parameters */
+    _e;
+    _a;
+    _m1;
+    _m2;
+    _t1;
+    _t2;
+    _tF;
+    _n;
+    _F;
+    _rF;
     /**
        * Create a new projection.
        * @param parameters the values of the parameters.
@@ -167,10 +196,4 @@ export class LambertConical2SP extends OperationMethod {
         source.setZ(target.getZ()); // Keep the Z value
     }
 }
-/** The code of this method */
-LambertConical2SP.METHOD_CODE = 9802;
-/** The value of PI */
-LambertConical2SP.PI = Math.PI;
-/** The half value of PI */
-LambertConical2SP.hPI = 0.5 * Math.PI;
 //# sourceMappingURL=LambertConical2SP.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"LambertConical2SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical2SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,iBAAkB,SAAQ,eAAe;IAkClD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,iBAAiB,CAAC,WAAW,EAAE,+BAA+B,EAAE,UAAU,CAAC,CAAC;QAClF,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,MAAM,CAAC,CAAU;QAC5B,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxF,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;IAC3H,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,CAAU;QAC/D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM;YAAE,OAAO,GAAG,CAAC;QACrC,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClC,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,kCAAkC;QAClC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,0BAA0B;QAC1B,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAChG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACxE,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,GAAY,EAAE,GAAY,EAAE,SAAqB;QACjE,0BAA0B;QAC1B,MAAM,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QACzD,MAAM,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACzE,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACpD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7D,uBAAuB;QACvB,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,CAAU,EAAE,CAAU,EAAE,UAAsB;QAC9D,wBAAwB;QACxB,MAAM,CAAC,GAAY,CAAC,CAAC;QACrB,MAAM,CAAC,GAAY,CAAC,CAAC;QACrB,0BAA0B;QAC1B,MAAM,EAAE,GAAY,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,iBAAiB,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC9J,MAAM,EAAE,GAAY,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACtE,MAAM,MAAM,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAChF,IAAI,GAAG,GAAY,iBAAiB,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACpE,MAAM,EAAE,GAAY,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;QAClC,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,kDAAkD;SACtF,CAAC;YACG,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC9C,GAAG,GAAG,iBAAiB,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,oBAAoB;QACvH,CAAC;QACD,MAAM,GAAG,GAAY,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QACrD,uBAAuB;QACvB,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACzB,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,uBAAuB;QACvB,IAAI,CAAC,YAAY,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;QACpC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,+BAA+B;QAC/B,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAChC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;;AAvLD,8BAA8B;AACP,6BAAW,GAAU,IAAI,CAAC;AAEjD,sBAAsB;AACE,oBAAE,GAAY,IAAI,CAAC,EAAE,CAAC;AAC9C,2BAA2B;AACH,qBAAG,GAAY,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n// package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { Ellipsoid } from \"../Ellipsoid\";\nimport { Operation } from \"../Operation\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class LambertConical2SP defines a Lambert Conical Conformal map projection with 2 standard parallels.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised May 2005\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.1.1 -> 1.4.1.4\n *\n * @version 1.0 July 2005\n */\n/** @internal */\nexport class LambertConical2SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9802;\n\n    /** The value of PI */\n    private static readonly PI: float64 = Math.PI;\n    /** The half value of PI */\n    private static readonly hPI: float64 = 0.5 * Math.PI;\n\n    /** Latitude of false origin */\n    private _latF: float64;\n    /** Longitude of false origin (the central meridian) */\n    private _lonF: float64;\n    /** Latitude of first standard parallel */\n    private _lat1: float64;\n    /** Latitude of second standard parallel */\n    private _lat2: float64;\n    /** Easting at false origin */\n    private _eF: float64;\n    /** Northing at false origin */\n    private _nF: float64;\n\n    /* Calculation parameters */\n    private _e: float64;\n    private _a: float64;\n    private _m1: float64;\n    private _m2: float64;\n    private _t1: float64;\n    private _t2: float64;\n    private _tF: float64;\n    private _n: float64;\n    private _F: float64;\n    private _rF: float64;\n\n    /**\n       * Create a new projection.\n       * @param parameters the values of the parameters.\n       */\n    public constructor(parameters: ParameterValueList) {\n        super(LambertConical2SP.METHOD_CODE, \"Lambert Conic Conformal (2SP)\", parameters);\n        /* Store the parameters */\n        this._latF = parameters.getValue(8821);\n        this._lonF = parameters.getValue(8822);\n        this._lat1 = parameters.getValue(8823);\n        this._lat2 = parameters.getValue(8824);\n        this._eF = parameters.getValue(8826);\n        this._nF = parameters.getValue(8827);\n    }\n\n    /**\n       * Get the sign of a number.\n       */\n    private static sign(v: float64): float64 {\n        return (v < 0.0) ? (-1.0) : (1.0);\n    }\n\n    /**\n       * Get the square of a number.\n       */\n    private static square(v: float64): float64 {\n        return (v * v);\n    }\n\n    /**\n       * Calculate M.\n       */\n    private static calcM(e: float64, lat: float64): float64 {\n        return Math.cos(lat) / Math.sqrt(1.0 - LambertConical2SP.square(e * Math.sin(lat)));\n    }\n\n    /**\n       * Calculate T.\n       */\n    private static calcT(e: float64, lat: float64): float64 {\n        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);\n    }\n\n    /**\n       * Calculate R.\n       */\n    private static calcR(a: float64, F: float64, t: float64, n: float64): float64 {\n        if (Math.abs(t) < 1.0e-6) return 0.0;\n        return a * F * Math.pow(t, n);\n    }\n\n    /**\n       * Initialize the projection.\n       * @param ellipsoid the ellipsoid to use.\n       * @return this projection (for convenience).\n       */\n    public initializeProjection(ellipsoid: Ellipsoid): LambertConical2SP {\n        /* Get the ellipsoid parameters */\n        this._e = ellipsoid.getE();\n        this._a = ellipsoid.getA();\n        /* Make the calculation */\n        this._m1 = LambertConical2SP.calcM(this._e, this._lat1);\n        this._m2 = LambertConical2SP.calcM(this._e, this._lat2);\n        this._t1 = LambertConical2SP.calcT(this._e, this._lat1);\n        this._t2 = LambertConical2SP.calcT(this._e, this._lat2);\n        this._tF = LambertConical2SP.calcT(this._e, this._latF);\n        this._n = (Math.log(this._m1) - Math.log(this._m2)) / (Math.log(this._t1) - Math.log(this._t2));\n        this._F = this._m1 / (this._n * Math.pow(this._t1, this._n));\n        this._rF = LambertConical2SP.calcR(this._a, this._F, this._tF, this._n);\n        /* Return the projection */\n        return this;\n    }\n\n    /**\n       * Do the projection.\n       * @param lon the longitude (radians).\n       * @param lat the latitude (radians).\n       * @param projected the target projected coordinate (X and Y will be set).\n       */\n    public toProjection(lon: float64, lat: float64, projected: Coordinate): void {\n        /* Make the calculation */\n        const t: float64 = LambertConical2SP.calcT(this._e, lat);\n        const r: float64 = LambertConical2SP.calcR(this._a, this._F, t, this._n);\n        const theta: float64 = this._n * (lon - this._lonF);\n        const E: float64 = this._eF + r * Math.sin(theta);\n        const N: float64 = this._nF + this._rF - r * Math.cos(theta);\n        /* Save the position */\n        projected.setX(E);\n        projected.setY(N);\n    }\n\n    /**\n       * Do the inverse projection.\n       * @param x the easting.\n       * @param y the northing.\n       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).\n       */\n    public toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = x;\n        const N: float64 = y;\n        /* Make the calculation */\n        const r_: float64 = LambertConical2SP.sign(this._n) * Math.sqrt(LambertConical2SP.square(E - this._eF) + LambertConical2SP.square(this._rF - (N - this._nF)));\n        const t_: float64 = Math.pow(r_ / (this._a * this._F), 1.0 / this._n);\n        const theta_: float64 = Math.atan((E - this._eF) / (this._rF - (N - this._nF)));\n        let lat: float64 = LambertConical2SP.PI * 0.5 - 2.0 * Math.atan(t_);\n        const he: float64 = 0.5 * this._e;\n        for (let i: number = 0; i < 7; i++) // double-checked iteration count. LER, 24/11/2011\n        {\n            const eSin: float64 = this._e * Math.sin(lat);\n            lat = LambertConical2SP.hPI - 2.0 * Math.atan(t_ * Math.pow((1.0 - eSin) / (1.0 + eSin), he)); // recursive formula\n        }\n        const lon: float64 = (theta_) / this._n + this._lonF;\n        /* Save the position */\n        geographic.setX(lon);\n        geographic.setY(lat);\n    }\n\n    /**\n       * OperationMethod method.\n       * @see OperationMethod#initialize\n       */\n    public override initialize(operation: Operation): void {\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#forward\n       */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const lon: float64 = source.getX();\n        const lat: float64 = source.getY();\n        /* Do the projection */\n        this.toProjection(lon, lat, target);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#reverse\n       */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = target.getX();\n        const N: float64 = target.getY();\n        /* Do the inverse projection */\n        this.toGeoGraphic(E, N, source);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}
+{"version":3,"file":"LambertConical2SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical2SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAcH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAM,OAAO,iBAAkB,SAAQ,eAAe;IAClD,8BAA8B;IACvB,MAAM,CAAU,WAAW,GAAU,IAAI,CAAC;IAEjD,sBAAsB;IACd,MAAM,CAAU,EAAE,GAAY,IAAI,CAAC,EAAE,CAAC;IAC9C,2BAA2B;IACnB,MAAM,CAAU,GAAG,GAAY,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;IAErD,+BAA+B;IACvB,KAAK,CAAU;IACvB,uDAAuD;IAC/C,KAAK,CAAU;IACvB,0CAA0C;IAClC,KAAK,CAAU;IACvB,2CAA2C;IACnC,KAAK,CAAU;IACvB,8BAA8B;IACtB,GAAG,CAAU;IACrB,+BAA+B;IACvB,GAAG,CAAU;IAErB,4BAA4B;IACpB,EAAE,CAAU;IACZ,EAAE,CAAU;IACZ,GAAG,CAAU;IACb,GAAG,CAAU;IACb,GAAG,CAAU;IACb,GAAG,CAAU;IACb,GAAG,CAAU;IACb,EAAE,CAAU;IACZ,EAAE,CAAU;IACZ,GAAG,CAAU;IAErB;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,iBAAiB,CAAC,WAAW,EAAE,+BAA+B,EAAE,UAAU,CAAC,CAAC;QAClF,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,KAAK,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACvC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,IAAI,CAAC,CAAU;QAC1B,OAAO,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,MAAM,CAAC,CAAU;QAC5B,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnB,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACxF,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,GAAY;QACzC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;IAC3H,CAAC;IAED;;SAEK;IACG,MAAM,CAAC,KAAK,CAAC,CAAU,EAAE,CAAU,EAAE,CAAU,EAAE,CAAU;QAC/D,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM;YAAE,OAAO,GAAG,CAAC;QACrC,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAClC,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,kCAAkC;QAClC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,0BAA0B;QAC1B,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACxD,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAChG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7D,IAAI,CAAC,GAAG,GAAG,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACxE,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,GAAY,EAAE,GAAY,EAAE,SAAqB;QACjE,0BAA0B;QAC1B,MAAM,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QACzD,MAAM,CAAC,GAAY,iBAAiB,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;QACzE,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACpD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAClD,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC7D,uBAAuB;QACvB,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAClB,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,CAAU,EAAE,CAAU,EAAE,UAAsB;QAC9D,wBAAwB;QACxB,MAAM,CAAC,GAAY,CAAC,CAAC;QACrB,MAAM,CAAC,GAAY,CAAC,CAAC;QACrB,0BAA0B;QAC1B,MAAM,EAAE,GAAY,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,iBAAiB,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC9J,MAAM,EAAE,GAAY,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACtE,MAAM,MAAM,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAChF,IAAI,GAAG,GAAY,iBAAiB,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACpE,MAAM,EAAE,GAAY,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;QAClC,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,kDAAkD;SACtF,CAAC;YACG,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC9C,GAAG,GAAG,iBAAiB,CAAC,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,oBAAoB;QACvH,CAAC;QACD,MAAM,GAAG,GAAY,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC;QACrD,uBAAuB;QACvB,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACzB,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,uBAAuB;QACvB,IAAI,CAAC,YAAY,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;QACpC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC;IAED;;;SAGK;IACE,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,MAAM,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,+BAA+B;QAC/B,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAChC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,mBAAmB;IACnD,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n// package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { Ellipsoid } from \"../Ellipsoid\";\nimport { Operation } from \"../Operation\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class LambertConical2SP defines a Lambert Conical Conformal map projection with 2 standard parallels.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised May 2005\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.1.1 -> 1.4.1.4\n *\n * @version 1.0 July 2005\n */\n/** @internal */\nexport class LambertConical2SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9802;\n\n    /** The value of PI */\n    private static readonly PI: float64 = Math.PI;\n    /** The half value of PI */\n    private static readonly hPI: float64 = 0.5 * Math.PI;\n\n    /** Latitude of false origin */\n    private _latF: float64;\n    /** Longitude of false origin (the central meridian) */\n    private _lonF: float64;\n    /** Latitude of first standard parallel */\n    private _lat1: float64;\n    /** Latitude of second standard parallel */\n    private _lat2: float64;\n    /** Easting at false origin */\n    private _eF: float64;\n    /** Northing at false origin */\n    private _nF: float64;\n\n    /* Calculation parameters */\n    private _e: float64;\n    private _a: float64;\n    private _m1: float64;\n    private _m2: float64;\n    private _t1: float64;\n    private _t2: float64;\n    private _tF: float64;\n    private _n: float64;\n    private _F: float64;\n    private _rF: float64;\n\n    /**\n       * Create a new projection.\n       * @param parameters the values of the parameters.\n       */\n    public constructor(parameters: ParameterValueList) {\n        super(LambertConical2SP.METHOD_CODE, \"Lambert Conic Conformal (2SP)\", parameters);\n        /* Store the parameters */\n        this._latF = parameters.getValue(8821);\n        this._lonF = parameters.getValue(8822);\n        this._lat1 = parameters.getValue(8823);\n        this._lat2 = parameters.getValue(8824);\n        this._eF = parameters.getValue(8826);\n        this._nF = parameters.getValue(8827);\n    }\n\n    /**\n       * Get the sign of a number.\n       */\n    private static sign(v: float64): float64 {\n        return (v < 0.0) ? (-1.0) : (1.0);\n    }\n\n    /**\n       * Get the square of a number.\n       */\n    private static square(v: float64): float64 {\n        return (v * v);\n    }\n\n    /**\n       * Calculate M.\n       */\n    private static calcM(e: float64, lat: float64): float64 {\n        return Math.cos(lat) / Math.sqrt(1.0 - LambertConical2SP.square(e * Math.sin(lat)));\n    }\n\n    /**\n       * Calculate T.\n       */\n    private static calcT(e: float64, lat: float64): float64 {\n        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);\n    }\n\n    /**\n       * Calculate R.\n       */\n    private static calcR(a: float64, F: float64, t: float64, n: float64): float64 {\n        if (Math.abs(t) < 1.0e-6) return 0.0;\n        return a * F * Math.pow(t, n);\n    }\n\n    /**\n       * Initialize the projection.\n       * @param ellipsoid the ellipsoid to use.\n       * @return this projection (for convenience).\n       */\n    public initializeProjection(ellipsoid: Ellipsoid): LambertConical2SP {\n        /* Get the ellipsoid parameters */\n        this._e = ellipsoid.getE();\n        this._a = ellipsoid.getA();\n        /* Make the calculation */\n        this._m1 = LambertConical2SP.calcM(this._e, this._lat1);\n        this._m2 = LambertConical2SP.calcM(this._e, this._lat2);\n        this._t1 = LambertConical2SP.calcT(this._e, this._lat1);\n        this._t2 = LambertConical2SP.calcT(this._e, this._lat2);\n        this._tF = LambertConical2SP.calcT(this._e, this._latF);\n        this._n = (Math.log(this._m1) - Math.log(this._m2)) / (Math.log(this._t1) - Math.log(this._t2));\n        this._F = this._m1 / (this._n * Math.pow(this._t1, this._n));\n        this._rF = LambertConical2SP.calcR(this._a, this._F, this._tF, this._n);\n        /* Return the projection */\n        return this;\n    }\n\n    /**\n       * Do the projection.\n       * @param lon the longitude (radians).\n       * @param lat the latitude (radians).\n       * @param projected the target projected coordinate (X and Y will be set).\n       */\n    public toProjection(lon: float64, lat: float64, projected: Coordinate): void {\n        /* Make the calculation */\n        const t: float64 = LambertConical2SP.calcT(this._e, lat);\n        const r: float64 = LambertConical2SP.calcR(this._a, this._F, t, this._n);\n        const theta: float64 = this._n * (lon - this._lonF);\n        const E: float64 = this._eF + r * Math.sin(theta);\n        const N: float64 = this._nF + this._rF - r * Math.cos(theta);\n        /* Save the position */\n        projected.setX(E);\n        projected.setY(N);\n    }\n\n    /**\n       * Do the inverse projection.\n       * @param x the easting.\n       * @param y the northing.\n       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).\n       */\n    public toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = x;\n        const N: float64 = y;\n        /* Make the calculation */\n        const r_: float64 = LambertConical2SP.sign(this._n) * Math.sqrt(LambertConical2SP.square(E - this._eF) + LambertConical2SP.square(this._rF - (N - this._nF)));\n        const t_: float64 = Math.pow(r_ / (this._a * this._F), 1.0 / this._n);\n        const theta_: float64 = Math.atan((E - this._eF) / (this._rF - (N - this._nF)));\n        let lat: float64 = LambertConical2SP.PI * 0.5 - 2.0 * Math.atan(t_);\n        const he: float64 = 0.5 * this._e;\n        for (let i: number = 0; i < 7; i++) // double-checked iteration count. LER, 24/11/2011\n        {\n            const eSin: float64 = this._e * Math.sin(lat);\n            lat = LambertConical2SP.hPI - 2.0 * Math.atan(t_ * Math.pow((1.0 - eSin) / (1.0 + eSin), he)); // recursive formula\n        }\n        const lon: float64 = (theta_) / this._n + this._lonF;\n        /* Save the position */\n        geographic.setX(lon);\n        geographic.setY(lat);\n    }\n\n    /**\n       * OperationMethod method.\n       * @see OperationMethod#initialize\n       */\n    public override initialize(operation: Operation): void {\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#forward\n       */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const lon: float64 = source.getX();\n        const lat: float64 = source.getY();\n        /* Do the projection */\n        this.toProjection(lon, lat, target);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#reverse\n       */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = target.getX();\n        const N: float64 = target.getY();\n        /* Do the inverse projection */\n        this.toGeoGraphic(E, N, source);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}

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

@@ -29,6 +29,20 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class Mercator1SP extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9804;
+    /** The value of PI/4 */
+    static QPI = Math.PI / 4.0;
+    /** Latitude of natural origin */
+    _latN;
+    /** Longitude of natural origin (the central meridian) */
+    _lonN;
+    /** Scale factor at natural origin */
+    _k0;
+    /** False easting */
+    _fE;
+    /** False northing */
+    _fN;
     /**
      * Create a new projection.
      * @param parameters the values of the parameters.
@@ -102,8 +116,4 @@ export class Mercator1SP extends OperationMethod {
         source.setZ(target.getZ()); // Keep the Z value
     }
 }
-/** The code of this method */
-Mercator1SP.METHOD_CODE = 9804;
-/** The value of PI/4 */
-Mercator1SP.QPI = Math.PI / 4.0;
 //# sourceMappingURL=Mercator1SP.js.map

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

@@ -1 +1 @@
-{"version":3,"file":"Mercator1SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/Mercator1SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAYH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,gBAAgB;AAChB,MAAM,OAAO,WAAY,SAAQ,eAAe;IAkB5C;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,WAAW,CAAC,WAAW,EAAE,gBAAgB,EAAE,UAAU,CAAC,CAAC;QAC7D,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1E,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1E,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACxE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACxE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IAC5E,CAAC;IAED;;;;;;;;OAQG;IACI,MAAM,CAAC,MAAM,CAAC,IAAa,EAAE,IAAa,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QACpF,IAAI,UAAU,GAAgB,IAAI,WAAW,CAAC,IAAI,CAAC,CAAC;QACpD,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,6CAA6C;QAC7C,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,uBAAuB;QACvB,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACnD,wBAAwB;QACxB,IAAI,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QAC1E,gEAAgE;QAChE,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,yBAAyB;QACzB,IAAI,GAAG,GAAY,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrD,wBAAwB;QACxB,IAAI,CAAC,GAAY,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACpC,IAAI,GAAG,GAAY,GAAG,GAAG,CAAC,WAAW,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,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;;AA3FD,8BAA8B;AACP,uBAAW,GAAU,IAAI,CAAC;AAEjD,wBAAwB;AACA,eAAG,GAAY,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n//package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class Mercator1SP defines the Mercator method with one standard parallel.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised January 2009\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.3.3.1 -> 1.3.3.2\n *\n * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)\n * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must\n * have a value of zero, for completeness in CRS labelling.\n\n * @version 1.0 April 2009\n */\n/** @internal */\nexport class Mercator1SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9804;\n\n    /** The value of PI/4 */\n    private static readonly QPI: float64 = Math.PI / 4.0;\n\n    /** Latitude of natural origin */\n    private _latN: float64;\n    /** Longitude of natural origin (the central meridian) */\n    private _lonN: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fE: float64;\n    /** False northing */\n    private _fN: float64;\n\n    /**\n     * Create a new projection.\n     * @param parameters the values of the parameters.\n     */\n    public constructor(parameters: ParameterValueList) {\n        super(Mercator1SP.METHOD_CODE, \"Mercator (1SP)\", parameters);\n        /* Store the parameters */\n        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;\n        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;\n        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;\n        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;\n        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;\n    }\n\n    /**\n     * Create a new projection.\n     * @param latN latitude of natural origin.\n     * @param lonN longitude of natural origin (the central meridian).\n     * @param k0 scale factor at natural origin.\n     * @param fE false easting.\n     * @param fN false northing.\n     * @return the projection.\n     */\n    public static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP {\n        let projection: Mercator1SP = new Mercator1SP(null);\n        projection._latN = latN;\n        projection._lonN = lonN;\n        projection._k0 = k0;\n        projection._fE = fE;\n        projection._fN = fN;\n        return projection;\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#forward\n     */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the geographic parameters (radians) */\n        let lon: float64 = source.getX();\n        let lat: float64 = source.getY();\n        /* Get the radius of the sphere */\n        let R: float64 = sourceCRS.getEllipsoid().getA();\n        /* Calculate easting */\n        let E: float64 = this._fE + R * (lon - this._lonN);\n        /* Calculate northing */\n        let sin: float64 = Math.sin(lat);\n        let N: float64 = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));\n        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#reverse\n     */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let E: float64 = target.getX();\n        let N: float64 = target.getY();\n        /* Get the radius of the sphere */\n        let R: float64 = sourceCRS.getEllipsoid().getA();\n        /* Calculate longitude */\n        let lon: float64 = ((E - this._fE) / R) + this._lonN;\n        /* Calculate latitude */\n        let D: float64 = (this._fN - N) / R;\n        let lat: float64 = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));\n        /* Save the position */\n        source.setX(lon);\n        source.setY(lat);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}
+{"version":3,"file":"Mercator1SP.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/Mercator1SP.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAYH,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AAGrD;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,gBAAgB;AAChB,MAAM,OAAO,WAAY,SAAQ,eAAe;IAC5C,8BAA8B;IACvB,MAAM,CAAU,WAAW,GAAU,IAAI,CAAC;IAEjD,wBAAwB;IAChB,MAAM,CAAU,GAAG,GAAY,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IAErD,iCAAiC;IACzB,KAAK,CAAU;IACvB,yDAAyD;IACjD,KAAK,CAAU;IACvB,qCAAqC;IAC7B,GAAG,CAAU;IACrB,oBAAoB;IACZ,GAAG,CAAU;IACrB,qBAAqB;IACb,GAAG,CAAU;IAErB;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,WAAW,CAAC,WAAW,EAAE,gBAAgB,EAAE,UAAU,CAAC,CAAC;QAC7D,0BAA0B;QAC1B,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1E,IAAI,CAAC,KAAK,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1E,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACxE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACxE,IAAI,CAAC,GAAG,GAAG,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,SAAS,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IAC5E,CAAC;IAED;;;;;;;;OAQG;IACI,MAAM,CAAC,MAAM,CAAC,IAAa,EAAE,IAAa,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QACpF,IAAI,UAAU,GAAgB,IAAI,WAAW,CAAC,IAAI,CAAC,CAAC;QACpD,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,6CAA6C;QAC7C,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,IAAI,GAAG,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QACjC,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,uBAAuB;QACvB,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACnD,wBAAwB;QACxB,IAAI,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QAC1E,gEAAgE;QAChE,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;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,wBAAwB;QACxB,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,kCAAkC;QAClC,IAAI,CAAC,GAAY,SAAS,CAAC,YAAY,EAAE,CAAC,IAAI,EAAE,CAAC;QACjD,yBAAyB;QACzB,IAAI,GAAG,GAAY,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;QACrD,wBAAwB;QACxB,IAAI,CAAC,GAAY,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACpC,IAAI,GAAG,GAAY,GAAG,GAAG,CAAC,WAAW,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,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","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n//package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class Mercator1SP defines the Mercator method with one standard parallel.\n *\n * The 'source' CRS is the geographic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised January 2009\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.3.3.1 -> 1.3.3.2\n *\n * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)\n * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must\n * have a value of zero, for completeness in CRS labelling.\n\n * @version 1.0 April 2009\n */\n/** @internal */\nexport class Mercator1SP extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9804;\n\n    /** The value of PI/4 */\n    private static readonly QPI: float64 = Math.PI / 4.0;\n\n    /** Latitude of natural origin */\n    private _latN: float64;\n    /** Longitude of natural origin (the central meridian) */\n    private _lonN: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fE: float64;\n    /** False northing */\n    private _fN: float64;\n\n    /**\n     * Create a new projection.\n     * @param parameters the values of the parameters.\n     */\n    public constructor(parameters: ParameterValueList) {\n        super(Mercator1SP.METHOD_CODE, \"Mercator (1SP)\", parameters);\n        /* Store the parameters */\n        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;\n        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;\n        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;\n        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;\n        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;\n    }\n\n    /**\n     * Create a new projection.\n     * @param latN latitude of natural origin.\n     * @param lonN longitude of natural origin (the central meridian).\n     * @param k0 scale factor at natural origin.\n     * @param fE false easting.\n     * @param fN false northing.\n     * @return the projection.\n     */\n    public static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP {\n        let projection: Mercator1SP = new Mercator1SP(null);\n        projection._latN = latN;\n        projection._lonN = lonN;\n        projection._k0 = k0;\n        projection._fE = fE;\n        projection._fN = fN;\n        return projection;\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#forward\n     */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the geographic parameters (radians) */\n        let lon: float64 = source.getX();\n        let lat: float64 = source.getY();\n        /* Get the radius of the sphere */\n        let R: float64 = sourceCRS.getEllipsoid().getA();\n        /* Calculate easting */\n        let E: float64 = this._fE + R * (lon - this._lonN);\n        /* Calculate northing */\n        let sin: float64 = Math.sin(lat);\n        let N: float64 = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));\n        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n     * OperationMethod interface method.\n     * @see OperationMethod#reverse\n     */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        let E: float64 = target.getX();\n        let N: float64 = target.getY();\n        /* Get the radius of the sphere */\n        let R: float64 = sourceCRS.getEllipsoid().getA();\n        /* Calculate longitude */\n        let lon: float64 = ((E - this._fE) / R) + this._lonN;\n        /* Calculate latitude */\n        let D: float64 = (this._fN - N) / R;\n        let lat: float64 = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));\n        /* Save the position */\n        source.setX(lon);\n        source.setY(lat);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}

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

@@ -25,6 +25,44 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class ObliqueMercator extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9815; // Hotine Variant B
+    /** The value of PI */
+    static PI = Math.PI;
+    /** Latitude of the projection center */
+    _latC;
+    /** Longitude of the projection center */
+    _lonC;
+    /** Azimuth of the initial line passing through the projection center */
+    _aziC;
+    /** Angle from the recified grid to the skew (oblique) grid */
+    _gamC;
+    /** Scale factor on the initial line of projection */
+    _kC;
+    /** False easting at the center of projection */
+    _eC;
+    /** False northing at the center of projection */
+    _nC;
+    /** Forward calculation parameters */
+    _a;
+    _e;
+    _e2;
+    _B;
+    _A;
+    _tO;
+    _D;
+    _F;
+    _H;
+    _G;
+    _gamO;
+    _lonO;
+    _vC;
+    _uC;
+    /** Reverse calculation parameters */
+    _rev1;
+    _rev2;
+    _rev3;
+    _rev4;
     /**
        * Create a new projection.
        * @param parameters the values of the parameters.
@@ -163,8 +201,4 @@ export class ObliqueMercator extends OperationMethod {
         source.setZ(target.getZ()); // Keep the Z value
     }
 }
-/** The code of this method */
-ObliqueMercator.METHOD_CODE = 9815; // Hotine Variant B
-/** The value of PI */
-ObliqueMercator.PI = Math.PI;
 //# sourceMappingURL=ObliqueMercator.js.map

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

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

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

@@ -25,6 +25,30 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class ObliqueStereographic extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9809;
+    /** The value of PI */
+    static PI = Math.PI;
+    /** Latitude of false origin */
+    _lat0;
+    /** Longitude of false origin */
+    _lon0;
+    /** Scale factor at natural origin */
+    _k0;
+    /** False easting */
+    _fe;
+    /** False northing */
+    _fn;
+    /** The parameters for the conformal sphere */
+    _e;
+    _e2;
+    _R;
+    _n;
+    _c;
+    _chi0;
+    _lambda0;
+    _g;
+    _h;
     /**
        * Create a new projection.
        * @param parameters the values of the parameters.
@@ -198,8 +222,4 @@ export class ObliqueStereographic extends OperationMethod {
         source.setZ(target.getZ()); // Keep the Z value
     }
 }
-/** The code of this method */
-ObliqueStereographic.METHOD_CODE = 9809;
-/** The value of PI */
-ObliqueStereographic.PI = Math.PI;
 //# sourceMappingURL=ObliqueStereographic.js.map

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,CAAC;YACjC,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;QAC/D,CAAC;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":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n// package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { Ellipsoid } from \"../Ellipsoid\";\nimport { Operation } from \"../Operation\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class ObliqueStereographic defines an Oblique Stereographic projection.\n *\n * The 'source' CRS is the geodetic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised August 2006\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.7.1 'Oblique and Equatorial Stereographic cases'\n *\n * @version 1.0 December 2006\n */\n/** @internal */\nexport class ObliqueStereographic extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9809;\n\n    /** The value of PI */\n    private static readonly PI: float64 = Math.PI;\n\n    /** Latitude of false origin */\n    private _lat0: float64;\n    /** Longitude of false origin */\n    private _lon0: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fe: float64;\n    /** False northing */\n    private _fn: float64;\n\n    /** The parameters for the conformal sphere */\n    private _e: float64;\n    private _e2: float64;\n    private _R: float64;\n    private _n: float64;\n    private _c: float64;\n    private _chi0: float64;\n    private _lambda0: float64;\n    private _g: float64;\n    private _h: float64;\n\n    /**\n       * Create a new projection.\n       * @param parameters the values of the parameters.\n       */\n    public constructor(parameters: ParameterValueList) {\n        super(ObliqueStereographic.METHOD_CODE, \"Oblique Stereographic\", parameters);\n        /* Store the parameters */\n        this._lat0 = (parameters != null) ? parameters.getValue(8801) : 0.0;\n        this._lon0 = (parameters != null) ? parameters.getValue(8802) : 0.0;\n        this._k0 = (parameters != null) ? parameters.getValue(8805) : 0.0;\n        this._fe = (parameters != null) ? parameters.getValue(8806) : 0.0;\n        this._fn = (parameters != null) ? parameters.getValue(8807) : 0.0;\n    }\n\n    /**\n       * Create a new projection.\n       * @param lat0 latitude of natural origin (radians).\n       * @param lon0 longitude of natural origin (radians).\n       * @param k0 scale factor at natural origin.\n       * @param fe false easting.\n       * @param fn false northing.\n       */\n    public static create(lat0: float64, lon0: float64, k0: float64, fe: float64, fn: float64): ObliqueStereographic {\n        const projection: ObliqueStereographic = new ObliqueStereographic(null);\n        projection._lat0 = lat0;\n        projection._lon0 = lon0;\n        projection._k0 = k0;\n        projection._fe = fe;\n        projection._fn = fn;\n        return projection;\n    }\n\n    /**\n       * Get the power of a number.\n       */\n    private static pow(n: float64, e: float64): float64 {\n        return Math.exp(e * Math.log(n));\n    }\n\n    /**\n       * Calculate sinus.\n       */\n    private static sin(a: float64): float64 {\n        return Math.sin(a);\n    }\n\n    /**\n       * Calculate cosinus.\n       */\n    private static cos(a: float64): float64 {\n        return Math.cos(a);\n    }\n\n    /**\n       * Calculate tangent.\n       */\n    private static tan(a: float64): float64 {\n        return Math.tan(a);\n    }\n\n    /**\n       * Calculate arcsinus.\n       */\n    private static asin(v: float64): float64 {\n        return Math.asin(v);\n    }\n\n    /**\n       * Calculate arctangent.\n       */\n    private static atan(v: float64): float64 {\n        return Math.atan(v);\n    }\n\n    /**\n       * Calculate log.\n       */\n    private static log(v: float64): float64 {\n        return Math.log(v);\n    }\n\n    /**\n       * Calculate exp.\n       */\n    private static exp(v: float64): float64 {\n        return Math.exp(v);\n    }\n\n    /**\n       * Initialize the projection.\n       * @param ellipsoid the ellipsoid to use.\n       * @return this projection (for convenience).\n       */\n    public initializeProjection(ellipsoid: Ellipsoid): ObliqueStereographic {\n        /* Get the ellipsoid parameters */\n        const a: float64 = ellipsoid.getA();\n        this._e = ellipsoid.getE();\n        /* Calculate conformat latitude (chi0) and longitude (lambda0) */\n        this._e2 = this._e * this._e;\n        const sinLat0: float64 = ObliqueStereographic.sin(this._lat0);\n        const t0: float64 = 1.0 - this._e2 * sinLat0 * sinLat0;\n        const rho0: float64 = a * (1 - this._e2) / ObliqueStereographic.pow(t0, 1.5);\n        const nu0: float64 = a / ObliqueStereographic.pow(t0, 0.5);\n        this._R = ObliqueStereographic.pow(rho0 * nu0, 0.5);\n        this._n = ObliqueStereographic.pow(1.0 + (this._e2 * ObliqueStereographic.pow(ObliqueStereographic.cos(this._lat0), 4)) / (1.0 - this._e2), 0.5);\n        const S1: float64 = (1.0 + sinLat0) / (1.0 - sinLat0);\n        const S2: float64 = (1.0 - this._e * sinLat0) / (1.0 + this._e * sinLat0);\n        const w1: float64 = ObliqueStereographic.pow(S1 * ObliqueStereographic.pow(S2, this._e), this._n);\n        const sinChi0: float64 = (w1 - 1.0) / (w1 + 1.0);\n        this._c = (this._n + sinLat0) * (1.0 - sinChi0) / ((this._n - sinLat0) * (1.0 + sinChi0));\n        const w2: float64 = this._c * w1;\n        this._chi0 = ObliqueStereographic.asin((w2 - 1.0) / (w2 + 1.0));\n        this._lambda0 = this._lon0;\n        /* For reverse calculation */\n        this._g = 2.0 * this._R * this._k0 * ObliqueStereographic.tan(ObliqueStereographic.PI / 4.0 - this._chi0 / 2.0);\n        this._h = 4.0 * this._R * this._k0 * ObliqueStereographic.tan(this._chi0) + this._g;\n        /* Return the projection */\n        return this;\n    }\n\n    /**\n       * OperationMethod method.\n       * @see OperationMethod#initialize\n       */\n    public override initialize(operation: Operation): void {\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#forward\n       */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const lon: float64 = source.getX();\n        const lat: float64 = source.getY();\n        const sinLat: float64 = ObliqueStereographic.sin(lat);\n        /* Calculate conformal latitude and longitude */\n        const lambda: float64 = this._n * (lon - this._lambda0) + this._lambda0;\n        const Sa: float64 = (1.0 + sinLat) / (1.0 - sinLat);\n        const Sb: float64 = (1.0 - this._e * sinLat) / (1.0 + this._e * sinLat);\n        const w: float64 = this._c * ObliqueStereographic.pow(Sa * ObliqueStereographic.pow(Sb, this._e), this._n);\n        const chi: float64 = ObliqueStereographic.asin((w - 1.0) / (w + 1.0));\n        /* Calculate easting and northing */\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);\n        const E: float64 = this._fe + 2.0 * this._R * this._k0 * ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(lambda - this._lambda0) / B;\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;\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#reverse\n       */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = target.getX() - this._fe;\n        const N: float64 = target.getY() - this._fn;\n        /* Calculate conformal latitude and longitude */\n        const i: float64 = ObliqueStereographic.atan(E / (this._h + N));\n        const j: float64 = ObliqueStereographic.atan(E / (this._g - N)) - i;\n        const chi: float64 = this._chi0 + 2.0 * ObliqueStereographic.atan((N - E * ObliqueStereographic.tan(j / 2.0)) / (2.0 * this._R * this._k0));\n        const lambda: float64 = j + 2.0 * i + this._lambda0;\n        const lon: float64 = (lambda - this._lambda0) / this._n + this._lambda0;\n        /* Iterate to get isometric (psi) and geodetic (phi) latitude */\n        const sinChi: float64 = ObliqueStereographic.sin(chi);\n        const psi: float64 = 0.5 * ObliqueStereographic.log((1.0 + sinChi) / (this._c * (1.0 - sinChi))) / this._n;\n        let phiK: float64 = 2.0 * ObliqueStereographic.atan(ObliqueStereographic.exp(psi)) - ObliqueStereographic.PI / 2.0;\n        for (let k: number = 1; k < 8; k++) {\n            const esinPhi: float64 = this._e * ObliqueStereographic.sin(phiK);\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));\n            const phiNextK: float64 = phiK - (psiK - psi) * ObliqueStereographic.cos(phiK) * (1.0 - esinPhi * esinPhi) / (1.0 - this._e2);\n            const change: float64 = Math.abs(phiNextK - phiK);\n            phiK = phiNextK;\n            if (change < 1.0e-12) break; // normally after 4 iterations\n        }\n        /* Save the position */\n        source.setX(lon);\n        source.setY(phiK);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\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,MAAM,OAAO,oBAAqB,SAAQ,eAAe;IACrD,8BAA8B;IACvB,MAAM,CAAU,WAAW,GAAU,IAAI,CAAC;IAEjD,sBAAsB;IACd,MAAM,CAAU,EAAE,GAAY,IAAI,CAAC,EAAE,CAAC;IAE9C,+BAA+B;IACvB,KAAK,CAAU;IACvB,gCAAgC;IACxB,KAAK,CAAU;IACvB,qCAAqC;IAC7B,GAAG,CAAU;IACrB,oBAAoB;IACZ,GAAG,CAAU;IACrB,qBAAqB;IACb,GAAG,CAAU;IAErB,8CAA8C;IACtC,EAAE,CAAU;IACZ,GAAG,CAAU;IACb,EAAE,CAAU;IACZ,EAAE,CAAU;IACZ,EAAE,CAAU;IACZ,KAAK,CAAU;IACf,QAAQ,CAAU;IAClB,EAAE,CAAU;IACZ,EAAE,CAAU;IAEpB;;;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,CAAC;YACjC,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;QAC/D,CAAC;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","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module OrbitGT\n */\n\n// package orbitgt.spatial.ecrs.projection;\n\ntype int8 = number;\ntype int16 = number;\ntype int32 = number;\ntype float32 = number;\ntype float64 = number;\n\nimport { Coordinate } from \"../../geom/Coordinate\";\nimport { CRS } from \"../CRS\";\nimport { Ellipsoid } from \"../Ellipsoid\";\nimport { Operation } from \"../Operation\";\nimport { OperationMethod } from \"../OperationMethod\";\nimport { ParameterValueList } from \"../ParameterValueList\";\n\n/**\n * Class ObliqueStereographic defines an Oblique Stereographic projection.\n *\n * The 'source' CRS is the geodetic CRS.\n * The 'target' CRS is the projected CRS.\n *\n * Based on the following document:\n *\n * Coordinate Conversions and Transformations including Formulas\n * Guidance Note Number 7, part 2\n * Revised August 2006\n * Available at: http://www.epsg.org/\n *\n * Formulas: see 1.4.7.1 'Oblique and Equatorial Stereographic cases'\n *\n * @version 1.0 December 2006\n */\n/** @internal */\nexport class ObliqueStereographic extends OperationMethod {\n    /** The code of this method */\n    public static readonly METHOD_CODE: int32 = 9809;\n\n    /** The value of PI */\n    private static readonly PI: float64 = Math.PI;\n\n    /** Latitude of false origin */\n    private _lat0: float64;\n    /** Longitude of false origin */\n    private _lon0: float64;\n    /** Scale factor at natural origin */\n    private _k0: float64;\n    /** False easting */\n    private _fe: float64;\n    /** False northing */\n    private _fn: float64;\n\n    /** The parameters for the conformal sphere */\n    private _e: float64;\n    private _e2: float64;\n    private _R: float64;\n    private _n: float64;\n    private _c: float64;\n    private _chi0: float64;\n    private _lambda0: float64;\n    private _g: float64;\n    private _h: float64;\n\n    /**\n       * Create a new projection.\n       * @param parameters the values of the parameters.\n       */\n    public constructor(parameters: ParameterValueList) {\n        super(ObliqueStereographic.METHOD_CODE, \"Oblique Stereographic\", parameters);\n        /* Store the parameters */\n        this._lat0 = (parameters != null) ? parameters.getValue(8801) : 0.0;\n        this._lon0 = (parameters != null) ? parameters.getValue(8802) : 0.0;\n        this._k0 = (parameters != null) ? parameters.getValue(8805) : 0.0;\n        this._fe = (parameters != null) ? parameters.getValue(8806) : 0.0;\n        this._fn = (parameters != null) ? parameters.getValue(8807) : 0.0;\n    }\n\n    /**\n       * Create a new projection.\n       * @param lat0 latitude of natural origin (radians).\n       * @param lon0 longitude of natural origin (radians).\n       * @param k0 scale factor at natural origin.\n       * @param fe false easting.\n       * @param fn false northing.\n       */\n    public static create(lat0: float64, lon0: float64, k0: float64, fe: float64, fn: float64): ObliqueStereographic {\n        const projection: ObliqueStereographic = new ObliqueStereographic(null);\n        projection._lat0 = lat0;\n        projection._lon0 = lon0;\n        projection._k0 = k0;\n        projection._fe = fe;\n        projection._fn = fn;\n        return projection;\n    }\n\n    /**\n       * Get the power of a number.\n       */\n    private static pow(n: float64, e: float64): float64 {\n        return Math.exp(e * Math.log(n));\n    }\n\n    /**\n       * Calculate sinus.\n       */\n    private static sin(a: float64): float64 {\n        return Math.sin(a);\n    }\n\n    /**\n       * Calculate cosinus.\n       */\n    private static cos(a: float64): float64 {\n        return Math.cos(a);\n    }\n\n    /**\n       * Calculate tangent.\n       */\n    private static tan(a: float64): float64 {\n        return Math.tan(a);\n    }\n\n    /**\n       * Calculate arcsinus.\n       */\n    private static asin(v: float64): float64 {\n        return Math.asin(v);\n    }\n\n    /**\n       * Calculate arctangent.\n       */\n    private static atan(v: float64): float64 {\n        return Math.atan(v);\n    }\n\n    /**\n       * Calculate log.\n       */\n    private static log(v: float64): float64 {\n        return Math.log(v);\n    }\n\n    /**\n       * Calculate exp.\n       */\n    private static exp(v: float64): float64 {\n        return Math.exp(v);\n    }\n\n    /**\n       * Initialize the projection.\n       * @param ellipsoid the ellipsoid to use.\n       * @return this projection (for convenience).\n       */\n    public initializeProjection(ellipsoid: Ellipsoid): ObliqueStereographic {\n        /* Get the ellipsoid parameters */\n        const a: float64 = ellipsoid.getA();\n        this._e = ellipsoid.getE();\n        /* Calculate conformat latitude (chi0) and longitude (lambda0) */\n        this._e2 = this._e * this._e;\n        const sinLat0: float64 = ObliqueStereographic.sin(this._lat0);\n        const t0: float64 = 1.0 - this._e2 * sinLat0 * sinLat0;\n        const rho0: float64 = a * (1 - this._e2) / ObliqueStereographic.pow(t0, 1.5);\n        const nu0: float64 = a / ObliqueStereographic.pow(t0, 0.5);\n        this._R = ObliqueStereographic.pow(rho0 * nu0, 0.5);\n        this._n = ObliqueStereographic.pow(1.0 + (this._e2 * ObliqueStereographic.pow(ObliqueStereographic.cos(this._lat0), 4)) / (1.0 - this._e2), 0.5);\n        const S1: float64 = (1.0 + sinLat0) / (1.0 - sinLat0);\n        const S2: float64 = (1.0 - this._e * sinLat0) / (1.0 + this._e * sinLat0);\n        const w1: float64 = ObliqueStereographic.pow(S1 * ObliqueStereographic.pow(S2, this._e), this._n);\n        const sinChi0: float64 = (w1 - 1.0) / (w1 + 1.0);\n        this._c = (this._n + sinLat0) * (1.0 - sinChi0) / ((this._n - sinLat0) * (1.0 + sinChi0));\n        const w2: float64 = this._c * w1;\n        this._chi0 = ObliqueStereographic.asin((w2 - 1.0) / (w2 + 1.0));\n        this._lambda0 = this._lon0;\n        /* For reverse calculation */\n        this._g = 2.0 * this._R * this._k0 * ObliqueStereographic.tan(ObliqueStereographic.PI / 4.0 - this._chi0 / 2.0);\n        this._h = 4.0 * this._R * this._k0 * ObliqueStereographic.tan(this._chi0) + this._g;\n        /* Return the projection */\n        return this;\n    }\n\n    /**\n       * OperationMethod method.\n       * @see OperationMethod#initialize\n       */\n    public override initialize(operation: Operation): void {\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#forward\n       */\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const lon: float64 = source.getX();\n        const lat: float64 = source.getY();\n        const sinLat: float64 = ObliqueStereographic.sin(lat);\n        /* Calculate conformal latitude and longitude */\n        const lambda: float64 = this._n * (lon - this._lambda0) + this._lambda0;\n        const Sa: float64 = (1.0 + sinLat) / (1.0 - sinLat);\n        const Sb: float64 = (1.0 - this._e * sinLat) / (1.0 + this._e * sinLat);\n        const w: float64 = this._c * ObliqueStereographic.pow(Sa * ObliqueStereographic.pow(Sb, this._e), this._n);\n        const chi: float64 = ObliqueStereographic.asin((w - 1.0) / (w + 1.0));\n        /* Calculate easting and northing */\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);\n        const E: float64 = this._fe + 2.0 * this._R * this._k0 * ObliqueStereographic.cos(chi) * ObliqueStereographic.sin(lambda - this._lambda0) / B;\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;\n        /* Save the position */\n        target.setX(E);\n        target.setY(N);\n        target.setZ(source.getZ()); // Keep the Z value\n    }\n\n    /**\n       * OperationMethod interface method.\n       * @see OperationMethod#reverse\n       */\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\n        /* Get the parameters */\n        const E: float64 = target.getX() - this._fe;\n        const N: float64 = target.getY() - this._fn;\n        /* Calculate conformal latitude and longitude */\n        const i: float64 = ObliqueStereographic.atan(E / (this._h + N));\n        const j: float64 = ObliqueStereographic.atan(E / (this._g - N)) - i;\n        const chi: float64 = this._chi0 + 2.0 * ObliqueStereographic.atan((N - E * ObliqueStereographic.tan(j / 2.0)) / (2.0 * this._R * this._k0));\n        const lambda: float64 = j + 2.0 * i + this._lambda0;\n        const lon: float64 = (lambda - this._lambda0) / this._n + this._lambda0;\n        /* Iterate to get isometric (psi) and geodetic (phi) latitude */\n        const sinChi: float64 = ObliqueStereographic.sin(chi);\n        const psi: float64 = 0.5 * ObliqueStereographic.log((1.0 + sinChi) / (this._c * (1.0 - sinChi))) / this._n;\n        let phiK: float64 = 2.0 * ObliqueStereographic.atan(ObliqueStereographic.exp(psi)) - ObliqueStereographic.PI / 2.0;\n        for (let k: number = 1; k < 8; k++) {\n            const esinPhi: float64 = this._e * ObliqueStereographic.sin(phiK);\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));\n            const phiNextK: float64 = phiK - (psiK - psi) * ObliqueStereographic.cos(phiK) * (1.0 - esinPhi * esinPhi) / (1.0 - this._e2);\n            const change: float64 = Math.abs(phiNextK - phiK);\n            phiK = phiNextK;\n            if (change < 1.0e-12) break; // normally after 4 iterations\n        }\n        /* Save the position */\n        source.setX(lon);\n        source.setY(phiK);\n        source.setZ(target.getZ()); // Keep the Z value\n    }\n}\n"]}

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

@@ -25,6 +25,34 @@ import { OperationMethod } from "../OperationMethod";
  */
 /** @internal */
 export class TransverseMercator extends OperationMethod {
+    /** The code of this method */
+    static METHOD_CODE = 9807;
+    /** The value of PI */
+    static PI = Math.PI;
+    /** Latitude of natural origin (radians) */
+    _latN;
+    /** Longitude of natural origin (the central meridian) (radians) */
+    _lonN;
+    /** The scale factor at the natural origin */
+    _k0;
+    /** False easting */
+    _fE;
+    /** False northing */
+    _fN;
+    /** Forward calculation parameters */
+    _a;
+    _e;
+    _e2;
+    _e4;
+    _e6;
+    _e_2;
+    _m0;
+    /** Reverse calculation parameters */
+    _e1;
+    _R1;
+    _R2;
+    _R3;
+    _R4;
     /**
        * Create a new projection.
        * @param parameters the values of the parameters.
@@ -173,8 +201,4 @@ export class TransverseMercator extends OperationMethod {
         source.setZ(target.getZ()); // Copy the Z value
     }
 }
-/** The code of this method */
-TransverseMercator.METHOD_CODE = 9807;
-/** The value of PI */
-TransverseMercator.PI = Math.PI;
 //# sourceMappingURL=TransverseMercator.js.map