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

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;YACI,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;SACtH;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":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n// package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { Ellipsoid } from \"../Ellipsoid\";\r\nimport { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class LambertConical2SP defines a Lambert Conical Conformal map projection with 2 standard parallels.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised May 2005\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.1.1 -> 1.4.1.4\r\n *\r\n * @version 1.0 July 2005\r\n */\r\n/** @internal */\r\nexport class LambertConical2SP extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9802;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n    /** The half value of PI */\r\n    private static readonly hPI: float64 = 0.5 * Math.PI;\r\n\r\n    /** Latitude of false origin */\r\n    private _latF: float64;\r\n    /** Longitude of false origin (the central meridian) */\r\n    private _lonF: float64;\r\n    /** Latitude of first standard parallel */\r\n    private _lat1: float64;\r\n    /** Latitude of second standard parallel */\r\n    private _lat2: float64;\r\n    /** Easting at false origin */\r\n    private _eF: float64;\r\n    /** Northing at false origin */\r\n    private _nF: float64;\r\n\r\n    /* Calculation parameters */\r\n    private _e: float64;\r\n    private _a: float64;\r\n    private _m1: float64;\r\n    private _m2: float64;\r\n    private _t1: float64;\r\n    private _t2: float64;\r\n    private _tF: float64;\r\n    private _n: float64;\r\n    private _F: float64;\r\n    private _rF: float64;\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param parameters the values of the parameters.\r\n       */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(LambertConical2SP.METHOD_CODE, \"Lambert Conic Conformal (2SP)\", parameters);\r\n        /* Store the parameters */\r\n        this._latF = parameters.getValue(8821);\r\n        this._lonF = parameters.getValue(8822);\r\n        this._lat1 = parameters.getValue(8823);\r\n        this._lat2 = parameters.getValue(8824);\r\n        this._eF = parameters.getValue(8826);\r\n        this._nF = parameters.getValue(8827);\r\n    }\r\n\r\n    /**\r\n       * Get the sign of a number.\r\n       */\r\n    private static sign(v: float64): float64 {\r\n        return (v < 0.0) ? (-1.0) : (1.0);\r\n    }\r\n\r\n    /**\r\n       * Get the square of a number.\r\n       */\r\n    private static square(v: float64): float64 {\r\n        return (v * v);\r\n    }\r\n\r\n    /**\r\n       * Calculate M.\r\n       */\r\n    private static calcM(e: float64, lat: float64): float64 {\r\n        return Math.cos(lat) / Math.sqrt(1.0 - LambertConical2SP.square(e * Math.sin(lat)));\r\n    }\r\n\r\n    /**\r\n       * Calculate T.\r\n       */\r\n    private static calcT(e: float64, lat: float64): float64 {\r\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);\r\n    }\r\n\r\n    /**\r\n       * Calculate R.\r\n       */\r\n    private static calcR(a: float64, F: float64, t: float64, n: float64): float64 {\r\n        if (Math.abs(t) < 1.0e-6) return 0.0;\r\n        return a * F * Math.pow(t, n);\r\n    }\r\n\r\n    /**\r\n       * Initialize the projection.\r\n       * @param ellipsoid the ellipsoid to use.\r\n       * @return this projection (for convenience).\r\n       */\r\n    public initializeProjection(ellipsoid: Ellipsoid): LambertConical2SP {\r\n        /* Get the ellipsoid parameters */\r\n        this._e = ellipsoid.getE();\r\n        this._a = ellipsoid.getA();\r\n        /* Make the calculation */\r\n        this._m1 = LambertConical2SP.calcM(this._e, this._lat1);\r\n        this._m2 = LambertConical2SP.calcM(this._e, this._lat2);\r\n        this._t1 = LambertConical2SP.calcT(this._e, this._lat1);\r\n        this._t2 = LambertConical2SP.calcT(this._e, this._lat2);\r\n        this._tF = LambertConical2SP.calcT(this._e, this._latF);\r\n        this._n = (Math.log(this._m1) - Math.log(this._m2)) / (Math.log(this._t1) - Math.log(this._t2));\r\n        this._F = this._m1 / (this._n * Math.pow(this._t1, this._n));\r\n        this._rF = LambertConical2SP.calcR(this._a, this._F, this._tF, this._n);\r\n        /* Return the projection */\r\n        return this;\r\n    }\r\n\r\n    /**\r\n       * Do the projection.\r\n       * @param lon the longitude (radians).\r\n       * @param lat the latitude (radians).\r\n       * @param projected the target projected coordinate (X and Y will be set).\r\n       */\r\n    public toProjection(lon: float64, lat: float64, projected: Coordinate): void {\r\n        /* Make the calculation */\r\n        const t: float64 = LambertConical2SP.calcT(this._e, lat);\r\n        const r: float64 = LambertConical2SP.calcR(this._a, this._F, t, this._n);\r\n        const theta: float64 = this._n * (lon - this._lonF);\r\n        const E: float64 = this._eF + r * Math.sin(theta);\r\n        const N: float64 = this._nF + this._rF - r * Math.cos(theta);\r\n        /* Save the position */\r\n        projected.setX(E);\r\n        projected.setY(N);\r\n    }\r\n\r\n    /**\r\n       * Do the inverse projection.\r\n       * @param x the easting.\r\n       * @param y the northing.\r\n       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).\r\n       */\r\n    public toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = x;\r\n        const N: float64 = y;\r\n        /* Make the calculation */\r\n        const r_: float64 = LambertConical2SP.sign(this._n) * Math.sqrt(LambertConical2SP.square(E - this._eF) + LambertConical2SP.square(this._rF - (N - this._nF)));\r\n        const t_: float64 = Math.pow(r_ / (this._a * this._F), 1.0 / this._n);\r\n        const theta_: float64 = Math.atan((E - this._eF) / (this._rF - (N - this._nF)));\r\n        let lat: float64 = LambertConical2SP.PI * 0.5 - 2.0 * Math.atan(t_);\r\n        const he: float64 = 0.5 * this._e;\r\n        for (let i: number = 0; i < 7; i++) // double-checked iteration count. LER, 24/11/2011\r\n        {\r\n            const eSin: float64 = this._e * Math.sin(lat);\r\n            lat = LambertConical2SP.hPI - 2.0 * Math.atan(t_ * Math.pow((1.0 - eSin) / (1.0 + eSin), he)); // recursive formula\r\n        }\r\n        const lon: float64 = (theta_) / this._n + this._lonF;\r\n        /* Save the position */\r\n        geographic.setX(lon);\r\n        geographic.setY(lat);\r\n    }\r\n\r\n    /**\r\n       * OperationMethod method.\r\n       * @see OperationMethod#initialize\r\n       */\r\n    public override initialize(operation: Operation): void {\r\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#forward\r\n       */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const lon: float64 = source.getX();\r\n        const lat: float64 = source.getY();\r\n        /* Do the projection */\r\n        this.toProjection(lon, lat, target);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#reverse\r\n       */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Do the inverse projection */\r\n        this.toGeoGraphic(E, N, source);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\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,MAAa,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;YACI,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;SACtH;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;SAP5C,iBAAiB","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n// package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { Ellipsoid } from \"../Ellipsoid\";\r\nimport { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class LambertConical2SP defines a Lambert Conical Conformal map projection with 2 standard parallels.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised May 2005\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.1.1 -> 1.4.1.4\r\n *\r\n * @version 1.0 July 2005\r\n */\r\n/** @internal */\r\nexport class LambertConical2SP extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9802;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n    /** The half value of PI */\r\n    private static readonly hPI: float64 = 0.5 * Math.PI;\r\n\r\n    /** Latitude of false origin */\r\n    private _latF: float64;\r\n    /** Longitude of false origin (the central meridian) */\r\n    private _lonF: float64;\r\n    /** Latitude of first standard parallel */\r\n    private _lat1: float64;\r\n    /** Latitude of second standard parallel */\r\n    private _lat2: float64;\r\n    /** Easting at false origin */\r\n    private _eF: float64;\r\n    /** Northing at false origin */\r\n    private _nF: float64;\r\n\r\n    /* Calculation parameters */\r\n    private _e: float64;\r\n    private _a: float64;\r\n    private _m1: float64;\r\n    private _m2: float64;\r\n    private _t1: float64;\r\n    private _t2: float64;\r\n    private _tF: float64;\r\n    private _n: float64;\r\n    private _F: float64;\r\n    private _rF: float64;\r\n\r\n    /**\r\n       * Create a new projection.\r\n       * @param parameters the values of the parameters.\r\n       */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(LambertConical2SP.METHOD_CODE, \"Lambert Conic Conformal (2SP)\", parameters);\r\n        /* Store the parameters */\r\n        this._latF = parameters.getValue(8821);\r\n        this._lonF = parameters.getValue(8822);\r\n        this._lat1 = parameters.getValue(8823);\r\n        this._lat2 = parameters.getValue(8824);\r\n        this._eF = parameters.getValue(8826);\r\n        this._nF = parameters.getValue(8827);\r\n    }\r\n\r\n    /**\r\n       * Get the sign of a number.\r\n       */\r\n    private static sign(v: float64): float64 {\r\n        return (v < 0.0) ? (-1.0) : (1.0);\r\n    }\r\n\r\n    /**\r\n       * Get the square of a number.\r\n       */\r\n    private static square(v: float64): float64 {\r\n        return (v * v);\r\n    }\r\n\r\n    /**\r\n       * Calculate M.\r\n       */\r\n    private static calcM(e: float64, lat: float64): float64 {\r\n        return Math.cos(lat) / Math.sqrt(1.0 - LambertConical2SP.square(e * Math.sin(lat)));\r\n    }\r\n\r\n    /**\r\n       * Calculate T.\r\n       */\r\n    private static calcT(e: float64, lat: float64): float64 {\r\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);\r\n    }\r\n\r\n    /**\r\n       * Calculate R.\r\n       */\r\n    private static calcR(a: float64, F: float64, t: float64, n: float64): float64 {\r\n        if (Math.abs(t) < 1.0e-6) return 0.0;\r\n        return a * F * Math.pow(t, n);\r\n    }\r\n\r\n    /**\r\n       * Initialize the projection.\r\n       * @param ellipsoid the ellipsoid to use.\r\n       * @return this projection (for convenience).\r\n       */\r\n    public initializeProjection(ellipsoid: Ellipsoid): LambertConical2SP {\r\n        /* Get the ellipsoid parameters */\r\n        this._e = ellipsoid.getE();\r\n        this._a = ellipsoid.getA();\r\n        /* Make the calculation */\r\n        this._m1 = LambertConical2SP.calcM(this._e, this._lat1);\r\n        this._m2 = LambertConical2SP.calcM(this._e, this._lat2);\r\n        this._t1 = LambertConical2SP.calcT(this._e, this._lat1);\r\n        this._t2 = LambertConical2SP.calcT(this._e, this._lat2);\r\n        this._tF = LambertConical2SP.calcT(this._e, this._latF);\r\n        this._n = (Math.log(this._m1) - Math.log(this._m2)) / (Math.log(this._t1) - Math.log(this._t2));\r\n        this._F = this._m1 / (this._n * Math.pow(this._t1, this._n));\r\n        this._rF = LambertConical2SP.calcR(this._a, this._F, this._tF, this._n);\r\n        /* Return the projection */\r\n        return this;\r\n    }\r\n\r\n    /**\r\n       * Do the projection.\r\n       * @param lon the longitude (radians).\r\n       * @param lat the latitude (radians).\r\n       * @param projected the target projected coordinate (X and Y will be set).\r\n       */\r\n    public toProjection(lon: float64, lat: float64, projected: Coordinate): void {\r\n        /* Make the calculation */\r\n        const t: float64 = LambertConical2SP.calcT(this._e, lat);\r\n        const r: float64 = LambertConical2SP.calcR(this._a, this._F, t, this._n);\r\n        const theta: float64 = this._n * (lon - this._lonF);\r\n        const E: float64 = this._eF + r * Math.sin(theta);\r\n        const N: float64 = this._nF + this._rF - r * Math.cos(theta);\r\n        /* Save the position */\r\n        projected.setX(E);\r\n        projected.setY(N);\r\n    }\r\n\r\n    /**\r\n       * Do the inverse projection.\r\n       * @param x the easting.\r\n       * @param y the northing.\r\n       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).\r\n       */\r\n    public toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = x;\r\n        const N: float64 = y;\r\n        /* Make the calculation */\r\n        const r_: float64 = LambertConical2SP.sign(this._n) * Math.sqrt(LambertConical2SP.square(E - this._eF) + LambertConical2SP.square(this._rF - (N - this._nF)));\r\n        const t_: float64 = Math.pow(r_ / (this._a * this._F), 1.0 / this._n);\r\n        const theta_: float64 = Math.atan((E - this._eF) / (this._rF - (N - this._nF)));\r\n        let lat: float64 = LambertConical2SP.PI * 0.5 - 2.0 * Math.atan(t_);\r\n        const he: float64 = 0.5 * this._e;\r\n        for (let i: number = 0; i < 7; i++) // double-checked iteration count. LER, 24/11/2011\r\n        {\r\n            const eSin: float64 = this._e * Math.sin(lat);\r\n            lat = LambertConical2SP.hPI - 2.0 * Math.atan(t_ * Math.pow((1.0 - eSin) / (1.0 + eSin), he)); // recursive formula\r\n        }\r\n        const lon: float64 = (theta_) / this._n + this._lonF;\r\n        /* Save the position */\r\n        geographic.setX(lon);\r\n        geographic.setY(lat);\r\n    }\r\n\r\n    /**\r\n       * OperationMethod method.\r\n       * @see OperationMethod#initialize\r\n       */\r\n    public override initialize(operation: Operation): void {\r\n        this.initializeProjection(operation.getSourceCRS().getEllipsoid());\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#forward\r\n       */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const lon: float64 = source.getX();\r\n        const lat: float64 = source.getY();\r\n        /* Do the projection */\r\n        this.toProjection(lon, lat, target);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n       * OperationMethod interface method.\r\n       * @see OperationMethod#reverse\r\n       */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        const E: float64 = target.getX();\r\n        const N: float64 = target.getY();\r\n        /* Do the inverse projection */\r\n        this.toGeoGraphic(E, N, source);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}

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

@@ -1,74 +1,74 @@
-/** @packageDocumentation
- * @module OrbitGT
- */
-declare type int32 = number;
-declare type float64 = number;
-import { Coordinate } from "../../geom/Coordinate";
-import { CRS } from "../CRS";
-import { OperationMethod } from "../OperationMethod";
-import { ParameterValueList } from "../ParameterValueList";
-/**
- * Class Mercator1SP defines the Mercator method with one standard parallel.
- *
- * The 'source' CRS is the geographic CRS.
- * The 'target' CRS is the projected CRS.
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised January 2009
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.3.3.1 -> 1.3.3.2
- *
- * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)
- * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must
- * have a value of zero, for completeness in CRS labelling.
-
- * @version 1.0 April 2009
- */
-/** @internal */
-export declare class Mercator1SP extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** The value of PI/4 */
-    private static readonly QPI;
-    /** Latitude of natural origin */
-    private _latN;
-    /** Longitude of natural origin (the central meridian) */
-    private _lonN;
-    /** Scale factor at natural origin */
-    private _k0;
-    /** False easting */
-    private _fE;
-    /** False northing */
-    private _fN;
-    /**
-     * Create a new projection.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters: ParameterValueList);
-    /**
-     * Create a new projection.
-     * @param latN latitude of natural origin.
-     * @param lonN longitude of natural origin (the central meridian).
-     * @param k0 scale factor at natural origin.
-     * @param fE false easting.
-     * @param fN false northing.
-     * @return the projection.
-     */
-    static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP;
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#forward
-     */
-    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#reverse
-     */
-    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-}
-export {};
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+type int32 = number;
+type float64 = number;
+import { Coordinate } from "../../geom/Coordinate";
+import { CRS } from "../CRS";
+import { OperationMethod } from "../OperationMethod";
+import { ParameterValueList } from "../ParameterValueList";
+/**
+ * Class Mercator1SP defines the Mercator method with one standard parallel.
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised January 2009
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.3.3.1 -> 1.3.3.2
+ *
+ * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)
+ * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must
+ * have a value of zero, for completeness in CRS labelling.
+
+ * @version 1.0 April 2009
+ */
+/** @internal */
+export declare class Mercator1SP extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The value of PI/4 */
+    private static readonly QPI;
+    /** Latitude of natural origin */
+    private _latN;
+    /** Longitude of natural origin (the central meridian) */
+    private _lonN;
+    /** Scale factor at natural origin */
+    private _k0;
+    /** False easting */
+    private _fE;
+    /** False northing */
+    private _fN;
+    /**
+     * Create a new projection.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters: ParameterValueList);
+    /**
+     * Create a new projection.
+     * @param latN latitude of natural origin.
+     * @param lonN longitude of natural origin (the central meridian).
+     * @param k0 scale factor at natural origin.
+     * @param fE false easting.
+     * @param fN false northing.
+     * @return the projection.
+     */
+    static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP;
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#forward
+     */
+    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#reverse
+     */
+    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+}
+export {};
 //# sourceMappingURL=Mercator1SP.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"Mercator1SP.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/Mercator1SP.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,aAAK,KAAK,GAAG,MAAM,CAAC;AAEpB,aAAK,OAAO,GAAG,MAAM,CAAC;AAEtB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,gBAAgB;AAChB,qBAAa,WAAY,SAAQ,eAAe;IAC5C,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,wBAAwB;IACxB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAA0B;IAErD,iCAAiC;IACjC,OAAO,CAAC,KAAK,CAAU;IACvB,yDAAyD;IACzD,OAAO,CAAC,KAAK,CAAU;IACvB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAUjD;;;;;;;;OAQG;WACW,MAAM,CAAC,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,WAAW;IAUtG;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAkB5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAgB/F"}
+{"version":3,"file":"Mercator1SP.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/Mercator1SP.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,KAAK,KAAK,GAAG,MAAM,CAAC;AAEpB,KAAK,OAAO,GAAG,MAAM,CAAC;AAEtB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;;;;;GAoBG;AACH,gBAAgB;AAChB,qBAAa,WAAY,SAAQ,eAAe;IAC5C,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,wBAAwB;IACxB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAA0B;IAErD,iCAAiC;IACjC,OAAO,CAAC,KAAK,CAAU;IACvB,yDAAyD;IACzD,OAAO,CAAC,KAAK,CAAU;IACvB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAUjD;;;;;;;;OAQG;WACW,MAAM,CAAC,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,WAAW;IAUtG;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAkB5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAgB/F"}

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

@@ -1,109 +1,110 @@
-/*---------------------------------------------------------------------------------------------
-* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
-* See LICENSE.md in the project root for license terms and full copyright notice.
-*--------------------------------------------------------------------------------------------*/
-/** @packageDocumentation
- * @module OrbitGT
- */
-import { OperationMethod } from "../OperationMethod";
-/**
- * Class Mercator1SP defines the Mercator method with one standard parallel.
- *
- * The 'source' CRS is the geographic CRS.
- * The 'target' CRS is the projected CRS.
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised January 2009
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.3.3.1 -> 1.3.3.2
- *
- * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)
- * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must
- * have a value of zero, for completeness in CRS labelling.
-
- * @version 1.0 April 2009
- */
-/** @internal */
-export class Mercator1SP extends OperationMethod {
-    /**
-     * Create a new projection.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters) {
-        super(Mercator1SP.METHOD_CODE, "Mercator (1SP)", parameters);
-        /* Store the parameters */
-        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;
-        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;
-        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;
-        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;
-        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;
-    }
-    /**
-     * Create a new projection.
-     * @param latN latitude of natural origin.
-     * @param lonN longitude of natural origin (the central meridian).
-     * @param k0 scale factor at natural origin.
-     * @param fE false easting.
-     * @param fN false northing.
-     * @return the projection.
-     */
-    static create(latN, lonN, k0, fE, fN) {
-        let projection = new Mercator1SP(null);
-        projection._latN = latN;
-        projection._lonN = lonN;
-        projection._k0 = k0;
-        projection._fE = fE;
-        projection._fN = fN;
-        return projection;
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#forward
-     */
-    forward(sourceCRS, source, targetCRS, target) {
-        /* Get the geographic parameters (radians) */
-        let lon = source.getX();
-        let lat = source.getY();
-        /* Get the radius of the sphere */
-        let R = sourceCRS.getEllipsoid().getA();
-        /* Calculate easting */
-        let E = this._fE + R * (lon - this._lonN);
-        /* Calculate northing */
-        let sin = Math.sin(lat);
-        let N = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));
-        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));
-        /* Save the position */
-        target.setX(E);
-        target.setY(N);
-        target.setZ(source.getZ()); // Keep the Z value
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#reverse
-     */
-    reverse(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        let E = target.getX();
-        let N = target.getY();
-        /* Get the radius of the sphere */
-        let R = sourceCRS.getEllipsoid().getA();
-        /* Calculate longitude */
-        let lon = ((E - this._fE) / R) + this._lonN;
-        /* Calculate latitude */
-        let D = (this._fN - N) / R;
-        let lat = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));
-        /* Save the position */
-        source.setX(lon);
-        source.setY(lat);
-        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;
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+import { OperationMethod } from "../OperationMethod";
+/**
+ * Class Mercator1SP defines the Mercator method with one standard parallel.
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised January 2009
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.3.3.1 -> 1.3.3.2
+ *
+ * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)
+ * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must
+ * have a value of zero, for completeness in CRS labelling.
+
+ * @version 1.0 April 2009
+ */
+/** @internal */
+class Mercator1SP extends OperationMethod {
+    /**
+     * Create a new projection.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters) {
+        super(Mercator1SP.METHOD_CODE, "Mercator (1SP)", parameters);
+        /* Store the parameters */
+        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;
+        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;
+        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;
+        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;
+        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;
+    }
+    /**
+     * Create a new projection.
+     * @param latN latitude of natural origin.
+     * @param lonN longitude of natural origin (the central meridian).
+     * @param k0 scale factor at natural origin.
+     * @param fE false easting.
+     * @param fN false northing.
+     * @return the projection.
+     */
+    static create(latN, lonN, k0, fE, fN) {
+        let projection = new Mercator1SP(null);
+        projection._latN = latN;
+        projection._lonN = lonN;
+        projection._k0 = k0;
+        projection._fE = fE;
+        projection._fN = fN;
+        return projection;
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#forward
+     */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the geographic parameters (radians) */
+        let lon = source.getX();
+        let lat = source.getY();
+        /* Get the radius of the sphere */
+        let R = sourceCRS.getEllipsoid().getA();
+        /* Calculate easting */
+        let E = this._fE + R * (lon - this._lonN);
+        /* Calculate northing */
+        let sin = Math.sin(lat);
+        let N = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));
+        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));
+        /* Save the position */
+        target.setX(E);
+        target.setY(N);
+        target.setZ(source.getZ()); // Keep the Z value
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#reverse
+     */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        let E = target.getX();
+        let N = target.getY();
+        /* Get the radius of the sphere */
+        let R = sourceCRS.getEllipsoid().getA();
+        /* Calculate longitude */
+        let lon = ((E - this._fE) / R) + this._lonN;
+        /* Calculate latitude */
+        let D = (this._fN - N) / R;
+        let lat = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));
+        /* Save the position */
+        source.setX(lon);
+        source.setY(lat);
+        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;
+export { Mercator1SP };
 //# 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":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n//package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class Mercator1SP defines the Mercator method with one standard parallel.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised January 2009\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.3.3.1 -> 1.3.3.2\r\n *\r\n * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)\r\n * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must\r\n * have a value of zero, for completeness in CRS labelling.\r\n\r\n * @version 1.0 April 2009\r\n */\r\n/** @internal */\r\nexport class Mercator1SP extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9804;\r\n\r\n    /** The value of PI/4 */\r\n    private static readonly QPI: float64 = Math.PI / 4.0;\r\n\r\n    /** Latitude of natural origin */\r\n    private _latN: float64;\r\n    /** Longitude of natural origin (the central meridian) */\r\n    private _lonN: float64;\r\n    /** Scale factor at natural origin */\r\n    private _k0: float64;\r\n    /** False easting */\r\n    private _fE: float64;\r\n    /** False northing */\r\n    private _fN: float64;\r\n\r\n    /**\r\n     * Create a new projection.\r\n     * @param parameters the values of the parameters.\r\n     */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(Mercator1SP.METHOD_CODE, \"Mercator (1SP)\", parameters);\r\n        /* Store the parameters */\r\n        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;\r\n        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;\r\n        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;\r\n        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;\r\n        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;\r\n    }\r\n\r\n    /**\r\n     * Create a new projection.\r\n     * @param latN latitude of natural origin.\r\n     * @param lonN longitude of natural origin (the central meridian).\r\n     * @param k0 scale factor at natural origin.\r\n     * @param fE false easting.\r\n     * @param fN false northing.\r\n     * @return the projection.\r\n     */\r\n    public static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP {\r\n        let projection: Mercator1SP = new Mercator1SP(null);\r\n        projection._latN = latN;\r\n        projection._lonN = lonN;\r\n        projection._k0 = k0;\r\n        projection._fE = fE;\r\n        projection._fN = fN;\r\n        return projection;\r\n    }\r\n\r\n    /**\r\n     * OperationMethod interface method.\r\n     * @see OperationMethod#forward\r\n     */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the geographic parameters (radians) */\r\n        let lon: float64 = source.getX();\r\n        let lat: float64 = source.getY();\r\n        /* Get the radius of the sphere */\r\n        let R: float64 = sourceCRS.getEllipsoid().getA();\r\n        /* Calculate easting */\r\n        let E: float64 = this._fE + R * (lon - this._lonN);\r\n        /* Calculate northing */\r\n        let sin: float64 = Math.sin(lat);\r\n        let N: float64 = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));\r\n        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));\r\n        /* Save the position */\r\n        target.setX(E);\r\n        target.setY(N);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n     * OperationMethod interface method.\r\n     * @see OperationMethod#reverse\r\n     */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        let E: float64 = target.getX();\r\n        let N: float64 = target.getY();\r\n        /* Get the radius of the sphere */\r\n        let R: float64 = sourceCRS.getEllipsoid().getA();\r\n        /* Calculate longitude */\r\n        let lon: float64 = ((E - this._fE) / R) + this._lonN;\r\n        /* Calculate latitude */\r\n        let D: float64 = (this._fN - N) / R;\r\n        let lat: float64 = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(lat);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"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,MAAa,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;SAL5C,WAAW","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n//package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class Mercator1SP defines the Mercator method with one standard parallel.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised January 2009\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.3.3.1 -> 1.3.3.2\r\n *\r\n * Note that in these formulas common to both 1SP and 2SP cases, the parameter latitude of natural origin (latN)\r\n * is not used. However for the Merctor (1SP) method the EPSG dataset includes this parameter, which must\r\n * have a value of zero, for completeness in CRS labelling.\r\n\r\n * @version 1.0 April 2009\r\n */\r\n/** @internal */\r\nexport class Mercator1SP extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9804;\r\n\r\n    /** The value of PI/4 */\r\n    private static readonly QPI: float64 = Math.PI / 4.0;\r\n\r\n    /** Latitude of natural origin */\r\n    private _latN: float64;\r\n    /** Longitude of natural origin (the central meridian) */\r\n    private _lonN: float64;\r\n    /** Scale factor at natural origin */\r\n    private _k0: float64;\r\n    /** False easting */\r\n    private _fE: float64;\r\n    /** False northing */\r\n    private _fN: float64;\r\n\r\n    /**\r\n     * Create a new projection.\r\n     * @param parameters the values of the parameters.\r\n     */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(Mercator1SP.METHOD_CODE, \"Mercator (1SP)\", parameters);\r\n        /* Store the parameters */\r\n        this._latN = (parameters != null) ? parameters.getValue2(8801, 0.0) : 0.0;\r\n        this._lonN = (parameters != null) ? parameters.getValue2(8802, 0.0) : 0.0;\r\n        this._k0 = (parameters != null) ? parameters.getValue2(8805, 1.0) : 1.0;\r\n        this._fE = (parameters != null) ? parameters.getValue2(8806, 0.0) : 0.0;\r\n        this._fN = (parameters != null) ? parameters.getValue2(8807, 0.0) : 0.0;\r\n    }\r\n\r\n    /**\r\n     * Create a new projection.\r\n     * @param latN latitude of natural origin.\r\n     * @param lonN longitude of natural origin (the central meridian).\r\n     * @param k0 scale factor at natural origin.\r\n     * @param fE false easting.\r\n     * @param fN false northing.\r\n     * @return the projection.\r\n     */\r\n    public static create(latN: float64, lonN: float64, k0: float64, fE: float64, fN: float64): Mercator1SP {\r\n        let projection: Mercator1SP = new Mercator1SP(null);\r\n        projection._latN = latN;\r\n        projection._lonN = lonN;\r\n        projection._k0 = k0;\r\n        projection._fE = fE;\r\n        projection._fN = fN;\r\n        return projection;\r\n    }\r\n\r\n    /**\r\n     * OperationMethod interface method.\r\n     * @see OperationMethod#forward\r\n     */\r\n    public forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the geographic parameters (radians) */\r\n        let lon: float64 = source.getX();\r\n        let lat: float64 = source.getY();\r\n        /* Get the radius of the sphere */\r\n        let R: float64 = sourceCRS.getEllipsoid().getA();\r\n        /* Calculate easting */\r\n        let E: float64 = this._fE + R * (lon - this._lonN);\r\n        /* Calculate northing */\r\n        let sin: float64 = Math.sin(lat);\r\n        let N: float64 = this._fN + 0.5 * R * Math.log((1.0 + sin) / (1.0 - sin));\r\n        //        double N = this._fN+R*Math.log(Math.tan(QPI+0.5*lat));\r\n        /* Save the position */\r\n        target.setX(E);\r\n        target.setY(N);\r\n        target.setZ(source.getZ()); // Keep the Z value\r\n    }\r\n\r\n    /**\r\n     * OperationMethod interface method.\r\n     * @see OperationMethod#reverse\r\n     */\r\n    public reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void {\r\n        /* Get the parameters */\r\n        let E: float64 = target.getX();\r\n        let N: float64 = target.getY();\r\n        /* Get the radius of the sphere */\r\n        let R: float64 = sourceCRS.getEllipsoid().getA();\r\n        /* Calculate longitude */\r\n        let lon: float64 = ((E - this._fE) / R) + this._lonN;\r\n        /* Calculate latitude */\r\n        let D: float64 = (this._fN - N) / R;\r\n        let lat: float64 = 2.0 * (Mercator1SP.QPI - Math.atan(Math.exp(D)));\r\n        /* Save the position */\r\n        source.setX(lon);\r\n        source.setY(lat);\r\n        source.setZ(target.getZ()); // Keep the Z value\r\n    }\r\n}\r\n"]}