@itwin/core-orbitgt 4.0.0-dev.7 → 4.0.0-dev.72

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

@@ -1 +1 @@
-{"version":3,"file":"TransverseMercator.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/TransverseMercator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAcH,wDAAqD;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAa,kBAAmB,SAAQ,iCAAe;IAiCnD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,kBAAkB,CAAC,WAAW,EAAE,qBAAqB,EAAE,UAAU,CAAC,CAAC;QACzE,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;;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;;SAEK;IACG,KAAK,CAAC,GAAY;QACtB,MAAM,CAAC,GAAY,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACjS,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACvB,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,oCAAoC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC7B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC/B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC/B,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACxC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAClC,oCAAoC;QACpC,MAAM,IAAI,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC9C,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACvC,MAAM,IAAI,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC1C,MAAM,IAAI,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC;QACtC,MAAM,IAAI,GAAY,IAAI,GAAG,IAAI,CAAC;QAClC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC/C,IAAI,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC7B,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC;QAC/B,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,GAAY,EAAE,GAAY,EAAE,SAAqB;QACjE,0BAA0B;QAC1B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,GAAG,GAAG,GAAG,CAAC;QAC7B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QAC5C,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAY,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACnC,0BAA0B;QAC1B,MAAM,KAAK,GAAY,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;QACpK,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACrD,MAAM,KAAK,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QAC3N,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC;QACjD,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,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzD,MAAM,EAAE,GAAY,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QACnG,MAAM,IAAI,GAAY,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACjJ,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAChD,MAAM,IAAI,GAAY,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;QACrD,MAAM,EAAE,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC;QACnC,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACnF,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACpC,MAAM,EAAE,GAAY,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,IAAI,GAAY,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACpC,MAAM,EAAE,GAAY,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;QAC1C,MAAM,IAAI,GAAY,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpD,MAAM,GAAG,GAAY,CAAC,GAAG,CAAC,CAAC;QAC3B,0BAA0B;QAC1B,MAAM,KAAK,GAAY,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;QACrO,MAAM,GAAG,GAAY,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACxD,MAAM,KAAK,GAAY,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;QACrL,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QAChD,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,CAAC,2DAA2D;IACnI,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;;AA1LL,gDA2LC;AA1LG,8BAA8B;AACP,8BAAW,GAAU,IAAI,CAAC;AAEjD,sBAAsB;AACE,qBAAE,GAAY,IAAI,CAAC,EAAE,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n// package orbitgt.spatial.ecrs.projection;\r\n\r\ntype int8 = number;\r\ntype int16 = number;\r\ntype int32 = number;\r\ntype float32 = number;\r\ntype float64 = number;\r\n\r\nimport { Coordinate } from \"../../geom/Coordinate\";\r\nimport { CRS } from \"../CRS\";\r\nimport { Ellipsoid } from \"../Ellipsoid\";\r\nimport { Operation } from \"../Operation\";\r\nimport { OperationMethod } from \"../OperationMethod\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class TransverseMercator defines a Transverse Mercator projection.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised August 2006\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.5.1: General Case\r\n *\r\n * @version 1.0 January 2007\r\n */\r\n/** @internal */\r\nexport class TransverseMercator extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9807;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n\r\n    /** Latitude of natural origin (radians) */\r\n    private _latN: float64;\r\n    /** Longitude of natural origin (the central meridian) (radians) */\r\n    private _lonN: float64;\r\n    /** The scale factor at the 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    /** Forward calculation parameters */\r\n    private _a: float64;\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _e4: float64;\r\n    private _e6: float64;\r\n    private _e_2: float64;\r\n    private _m0: float64;\r\n    /** Reverse calculation parameters */\r\n    private _e1: float64;\r\n    private _R1: float64;\r\n    private _R2: float64;\r\n    private _R3: float64;\r\n    private _R4: 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(TransverseMercator.METHOD_CODE, \"Transverse Mercator\", parameters);\r\n        /* Store the parameters */\r\n        this._latN = parameters.getValue(8801);\r\n        this._lonN = parameters.getValue(8802);\r\n        this._k0 = parameters.getValue(8805);\r\n        this._fE = parameters.getValue(8806);\r\n        this._fN = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number (must be able to handle negative 'n' values).\r\n       */\r\n    private static pow(n: float64, e: int32): float64 {\r\n        let p: float64 = 1.0;\r\n        for (let i: number = 0; i < e; i++) p *= n;\r\n        return p;\r\n    }\r\n\r\n    /**\r\n       * Calculate M.\r\n       */\r\n    private calcM(lat: float64): float64 {\r\n        const M: float64 = lat * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256) - Math.sin(2.0 * lat) * (3 * this._e2 / 8 + 3 * this._e4 / 32 + 45 * this._e6 / 1024) + Math.sin(4.0 * lat) * (15 * this._e4 / 256 + 45 * this._e6 / 1024) - Math.sin(6.0 * lat) * (35 * this._e6 / 3072);\r\n        return this._a * M;\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): TransverseMercator {\r\n        /* Prepare the forward parameters */\r\n        this._a = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        this._e2 = this._e * this._e;\r\n        this._e4 = this._e2 * this._e2;\r\n        this._e6 = this._e2 * this._e4;\r\n        this._e_2 = this._e2 / (1.0 - this._e2);\r\n        this._m0 = this.calcM(this._latN);\r\n        /* Prepare the reverse parameters */\r\n        const temp: float64 = Math.sqrt(1 - this._e2);\r\n        this._e1 = (1.0 - temp) / (1.0 + temp);\r\n        const e1_2: float64 = this._e1 * this._e1;\r\n        const e1_3: float64 = this._e1 * e1_2;\r\n        const e1_4: float64 = e1_2 * e1_2;\r\n        this._R1 = (3 * this._e1 / 2 - 27 * e1_3 / 32);\r\n        this._R2 = (21 * e1_2 / 16 - 55 * e1_4 / 32);\r\n        this._R3 = (151 * e1_3 / 96);\r\n        this._R4 = (1097 * e1_4 / 512);\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        /* Calculate parameters */\r\n        const tan: float64 = Math.tan(lat);\r\n        const T: float64 = tan * tan;\r\n        const cos: float64 = Math.cos(lat);\r\n        const C: float64 = this._e2 * cos * cos / (1.0 - this._e2);\r\n        const A: float64 = (lon - this._lonN) * cos;\r\n        const eSin: float64 = this._e * Math.sin(lat);\r\n        const v: float64 = this._a / Math.sqrt(1.0 - eSin * eSin);\r\n        const m: float64 = this.calcM(lat);\r\n        /* Make the calculation */\r\n        const temp1: float64 = A + (1.0 - T + C) * TransverseMercator.pow(A, 3) / 6 + (5 - 18 * T + T * T + 72 * C - 58 * this._e_2) * TransverseMercator.pow(A, 5) / 120.0;\r\n        const E: float64 = this._fE + this._k0 * v * (temp1);\r\n        const temp2: float64 = m - this._m0 + v * tan * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * TransverseMercator.pow(A, 4) / 24.0 + (61 - 58 * T + T * T + 600 * C - 330 * this._e_2) * TransverseMercator.pow(A, 6) / 720.0);\r\n        const N: float64 = this._fN + this._k0 * (temp2);\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        /* Calculate parameters */\r\n        const m1: float64 = this._m0 + (N - this._fN) / this._k0;\r\n        const u1: float64 = m1 / (this._a * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256));\r\n        const lat1: float64 = u1 + this._R1 * Math.sin(2 * u1) + this._R2 * Math.sin(4 * u1) + this._R3 * Math.sin(6 * u1) + this._R4 * Math.sin(8 * u1);\r\n        const eSin1: float64 = this._e * Math.sin(lat1);\r\n        const temp: float64 = Math.sqrt(1.0 - eSin1 * eSin1);\r\n        const v1: float64 = this._a / temp;\r\n        const rho1: float64 = this._a * (1.0 - this._e2) / TransverseMercator.pow(temp, 3);\r\n        const tan: float64 = Math.tan(lat1);\r\n        const T1: float64 = tan * tan;\r\n        const T1_2: float64 = T1 * T1;\r\n        const cos: float64 = Math.cos(lat1);\r\n        const C1: float64 = this._e_2 * cos * cos;\r\n        const C1_2: float64 = C1 * C1;\r\n        const D: float64 = (E - this._fE) / (v1 * this._k0);\r\n        const D_2: float64 = D * D;\r\n        /* Make the calculation */\r\n        const temp1: float64 = D_2 / 2 - (5 + 3 * T1 + 10 * C1 - 4 * C1_2 - 9 * this._e_2) * TransverseMercator.pow(D, 4) / 24.0 + (61 + 90 * T1 + 298 * C1 + 45 * T1_2 - 252 * this._e_2 - 3 * C1_2) * TransverseMercator.pow(D, 6) / 720.0;\r\n        const lat: float64 = lat1 - (v1 * tan / rho1) * (temp1);\r\n        const temp2: float64 = D - (1 + 2 * T1 + C1) * TransverseMercator.pow(D, 3) / 6 + (5 - 2 * C1 + 28 * T1 - 3 * C1_2 + 8 * this._e_2 + 24 * T1_2) * TransverseMercator.pow(D, 5) / 120;\r\n        const lon: float64 = this._lonN + (temp2) / cos;\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()); // should this use getTargetCRS? (while testing OSTN15 RoI)\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()); // Copy 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()); // Copy the Z value\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"TransverseMercator.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/TransverseMercator.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAcH,wDAAqD;AAGrD;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,MAAa,kBAAmB,SAAQ,iCAAe;IAiCnD;;;SAGK;IACL,YAAmB,UAA8B;QAC7C,KAAK,CAAC,kBAAkB,CAAC,WAAW,EAAE,qBAAqB,EAAE,UAAU,CAAC,CAAC;QACzE,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;;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;;SAEK;IACG,KAAK,CAAC,GAAY;QACtB,MAAM,CAAC,GAAY,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACjS,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACvB,CAAC;IAED;;;;SAIK;IACE,oBAAoB,CAAC,SAAoB;QAC5C,oCAAoC;QACpC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,IAAI,EAAE,CAAC;QAC3B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC7B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC/B,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC/B,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACxC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAClC,oCAAoC;QACpC,MAAM,IAAI,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC9C,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;QACvC,MAAM,IAAI,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC;QAC1C,MAAM,IAAI,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC;QACtC,MAAM,IAAI,GAAY,IAAI,GAAG,IAAI,CAAC;QAClC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC/C,IAAI,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;QAC7B,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC;QAC/B,2BAA2B;QAC3B,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;;;SAKK;IACE,YAAY,CAAC,GAAY,EAAE,GAAY,EAAE,SAAqB;QACjE,0BAA0B;QAC1B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,GAAG,GAAG,GAAG,CAAC;QAC7B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QAC5C,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9C,MAAM,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAY,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QACnC,0BAA0B;QAC1B,MAAM,KAAK,GAAY,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;QACpK,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACrD,MAAM,KAAK,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QAC3N,MAAM,CAAC,GAAY,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC;QACjD,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,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzD,MAAM,EAAE,GAAY,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;QACnG,MAAM,IAAI,GAAY,EAAE,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACjJ,MAAM,KAAK,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAChD,MAAM,IAAI,GAAY,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;QACrD,MAAM,EAAE,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC;QACnC,MAAM,IAAI,GAAY,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACnF,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACpC,MAAM,EAAE,GAAY,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,IAAI,GAAY,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAY,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACpC,MAAM,EAAE,GAAY,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;QAC1C,MAAM,IAAI,GAAY,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,CAAC,GAAY,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpD,MAAM,GAAG,GAAY,CAAC,GAAG,CAAC,CAAC;QAC3B,0BAA0B;QAC1B,MAAM,KAAK,GAAY,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;QACrO,MAAM,GAAG,GAAY,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACxD,MAAM,KAAK,GAAY,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;QACrL,MAAM,GAAG,GAAY,IAAI,CAAC,KAAK,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;QAChD,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,CAAC,2DAA2D;IACnI,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;;AAzLD,8BAA8B;AACP,8BAAW,GAAU,IAAI,CAAC;AAEjD,sBAAsB;AACE,qBAAE,GAAY,IAAI,CAAC,EAAE,CAAC;AALrC,gDAAkB","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 TransverseMercator defines a Transverse Mercator projection.\r\n *\r\n * The 'source' CRS is the geographic CRS.\r\n * The 'target' CRS is the projected CRS.\r\n *\r\n * Based on the following document:\r\n *\r\n * Coordinate Conversions and Transformations including Formulas\r\n * Guidance Note Number 7, part 2\r\n * Revised August 2006\r\n * Available at: http://www.epsg.org/\r\n *\r\n * Formulas: see 1.4.5.1: General Case\r\n *\r\n * @version 1.0 January 2007\r\n */\r\n/** @internal */\r\nexport class TransverseMercator extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9807;\r\n\r\n    /** The value of PI */\r\n    private static readonly PI: float64 = Math.PI;\r\n\r\n    /** Latitude of natural origin (radians) */\r\n    private _latN: float64;\r\n    /** Longitude of natural origin (the central meridian) (radians) */\r\n    private _lonN: float64;\r\n    /** The scale factor at the 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    /** Forward calculation parameters */\r\n    private _a: float64;\r\n    private _e: float64;\r\n    private _e2: float64;\r\n    private _e4: float64;\r\n    private _e6: float64;\r\n    private _e_2: float64;\r\n    private _m0: float64;\r\n    /** Reverse calculation parameters */\r\n    private _e1: float64;\r\n    private _R1: float64;\r\n    private _R2: float64;\r\n    private _R3: float64;\r\n    private _R4: 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(TransverseMercator.METHOD_CODE, \"Transverse Mercator\", parameters);\r\n        /* Store the parameters */\r\n        this._latN = parameters.getValue(8801);\r\n        this._lonN = parameters.getValue(8802);\r\n        this._k0 = parameters.getValue(8805);\r\n        this._fE = parameters.getValue(8806);\r\n        this._fN = parameters.getValue(8807);\r\n    }\r\n\r\n    /**\r\n       * Get the power of a number (must be able to handle negative 'n' values).\r\n       */\r\n    private static pow(n: float64, e: int32): float64 {\r\n        let p: float64 = 1.0;\r\n        for (let i: number = 0; i < e; i++) p *= n;\r\n        return p;\r\n    }\r\n\r\n    /**\r\n       * Calculate M.\r\n       */\r\n    private calcM(lat: float64): float64 {\r\n        const M: float64 = lat * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256) - Math.sin(2.0 * lat) * (3 * this._e2 / 8 + 3 * this._e4 / 32 + 45 * this._e6 / 1024) + Math.sin(4.0 * lat) * (15 * this._e4 / 256 + 45 * this._e6 / 1024) - Math.sin(6.0 * lat) * (35 * this._e6 / 3072);\r\n        return this._a * M;\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): TransverseMercator {\r\n        /* Prepare the forward parameters */\r\n        this._a = ellipsoid.getA();\r\n        this._e = ellipsoid.getE();\r\n        this._e2 = this._e * this._e;\r\n        this._e4 = this._e2 * this._e2;\r\n        this._e6 = this._e2 * this._e4;\r\n        this._e_2 = this._e2 / (1.0 - this._e2);\r\n        this._m0 = this.calcM(this._latN);\r\n        /* Prepare the reverse parameters */\r\n        const temp: float64 = Math.sqrt(1 - this._e2);\r\n        this._e1 = (1.0 - temp) / (1.0 + temp);\r\n        const e1_2: float64 = this._e1 * this._e1;\r\n        const e1_3: float64 = this._e1 * e1_2;\r\n        const e1_4: float64 = e1_2 * e1_2;\r\n        this._R1 = (3 * this._e1 / 2 - 27 * e1_3 / 32);\r\n        this._R2 = (21 * e1_2 / 16 - 55 * e1_4 / 32);\r\n        this._R3 = (151 * e1_3 / 96);\r\n        this._R4 = (1097 * e1_4 / 512);\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        /* Calculate parameters */\r\n        const tan: float64 = Math.tan(lat);\r\n        const T: float64 = tan * tan;\r\n        const cos: float64 = Math.cos(lat);\r\n        const C: float64 = this._e2 * cos * cos / (1.0 - this._e2);\r\n        const A: float64 = (lon - this._lonN) * cos;\r\n        const eSin: float64 = this._e * Math.sin(lat);\r\n        const v: float64 = this._a / Math.sqrt(1.0 - eSin * eSin);\r\n        const m: float64 = this.calcM(lat);\r\n        /* Make the calculation */\r\n        const temp1: float64 = A + (1.0 - T + C) * TransverseMercator.pow(A, 3) / 6 + (5 - 18 * T + T * T + 72 * C - 58 * this._e_2) * TransverseMercator.pow(A, 5) / 120.0;\r\n        const E: float64 = this._fE + this._k0 * v * (temp1);\r\n        const temp2: float64 = m - this._m0 + v * tan * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * TransverseMercator.pow(A, 4) / 24.0 + (61 - 58 * T + T * T + 600 * C - 330 * this._e_2) * TransverseMercator.pow(A, 6) / 720.0);\r\n        const N: float64 = this._fN + this._k0 * (temp2);\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        /* Calculate parameters */\r\n        const m1: float64 = this._m0 + (N - this._fN) / this._k0;\r\n        const u1: float64 = m1 / (this._a * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256));\r\n        const lat1: float64 = u1 + this._R1 * Math.sin(2 * u1) + this._R2 * Math.sin(4 * u1) + this._R3 * Math.sin(6 * u1) + this._R4 * Math.sin(8 * u1);\r\n        const eSin1: float64 = this._e * Math.sin(lat1);\r\n        const temp: float64 = Math.sqrt(1.0 - eSin1 * eSin1);\r\n        const v1: float64 = this._a / temp;\r\n        const rho1: float64 = this._a * (1.0 - this._e2) / TransverseMercator.pow(temp, 3);\r\n        const tan: float64 = Math.tan(lat1);\r\n        const T1: float64 = tan * tan;\r\n        const T1_2: float64 = T1 * T1;\r\n        const cos: float64 = Math.cos(lat1);\r\n        const C1: float64 = this._e_2 * cos * cos;\r\n        const C1_2: float64 = C1 * C1;\r\n        const D: float64 = (E - this._fE) / (v1 * this._k0);\r\n        const D_2: float64 = D * D;\r\n        /* Make the calculation */\r\n        const temp1: float64 = D_2 / 2 - (5 + 3 * T1 + 10 * C1 - 4 * C1_2 - 9 * this._e_2) * TransverseMercator.pow(D, 4) / 24.0 + (61 + 90 * T1 + 298 * C1 + 45 * T1_2 - 252 * this._e_2 - 3 * C1_2) * TransverseMercator.pow(D, 6) / 720.0;\r\n        const lat: float64 = lat1 - (v1 * tan / rho1) * (temp1);\r\n        const temp2: float64 = D - (1 + 2 * T1 + C1) * TransverseMercator.pow(D, 3) / 6 + (5 - 2 * C1 + 28 * T1 - 3 * C1_2 + 8 * this._e_2 + 24 * T1_2) * TransverseMercator.pow(D, 5) / 120;\r\n        const lon: float64 = this._lonN + (temp2) / cos;\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()); // should this use getTargetCRS? (while testing OSTN15 RoI)\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()); // Copy 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()); // Copy the Z value\r\n    }\r\n}\r\n"]}

package/lib/cjs/spatial/ecrs/transformation/PositionVector.d.ts

@@ -1,91 +1,91 @@
-/** @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 PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).
- *
- * NOTE: this method is also known as "Bursa-Wolf".
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised May 2005
- * Available at: http://www.epsg.org/
- *
- * Used in Europe. To convert to "Coordinate Frame Rotation" inverse the sign of the three rotations.
- *
- * @version 1.0 July 2005
- */
-/** @internal */
-export declare class PositionVector extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** The translation vector */
-    private _dX;
-    private _dY;
-    private _dZ;
-    /** The rotations (radians) */
-    private _rX;
-    private _rY;
-    private _rZ;
-    /** The scale factor (offset from 1.0) */
-    private _dS;
-    /**
-     * Create a new transformation.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters: ParameterValueList);
-    /**
-     * Create a new transformation.
-     * @param dX the translation vector x.
-     * @param dY the translation vector y.
-     * @param dZ the translation vector z.
-     * @param rX the first rotation (radians).
-     * @param rY the second rotation (radians).
-     * @param rZ the third rotation (radians).
-     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).
-     * @return the new transformation.
-     */
-    static create(dX: float64, dY: float64, dZ: float64, rX: float64, rY: float64, rZ: float64, dS: float64): PositionVector;
-    /**
-     * Get the 7 parameters of the method.
-     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.
-     */
-    getParameters(): Float64Array;
-    /**
-     * Make a forward transformation.
-     * @param source the source coordinate (the coordinate to transform).
-     * @param target the target coordinate (holds the result).
-     */
-    transformForward(source: Coordinate, target: Coordinate): void;
-    /**
-     * Make a reverse transformation.
-     * @param source the source coordinate (holds the result).
-     * @param target the target coordinate (the coordinate to transform).
-     */
-    transformReverse(source: Coordinate, target: Coordinate): void;
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#forward
-     */
-    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#reverse
-     */
-    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-    /**
-     * Create an inverse transformation.
-     * @return the inverse transformation.
-     */
-    createInverse(): PositionVector;
-}
-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 PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).
+ *
+ * NOTE: this method is also known as "Bursa-Wolf".
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised May 2005
+ * Available at: http://www.epsg.org/
+ *
+ * Used in Europe. To convert to "Coordinate Frame Rotation" inverse the sign of the three rotations.
+ *
+ * @version 1.0 July 2005
+ */
+/** @internal */
+export declare class PositionVector extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The translation vector */
+    private _dX;
+    private _dY;
+    private _dZ;
+    /** The rotations (radians) */
+    private _rX;
+    private _rY;
+    private _rZ;
+    /** The scale factor (offset from 1.0) */
+    private _dS;
+    /**
+     * Create a new transformation.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters: ParameterValueList);
+    /**
+     * Create a new transformation.
+     * @param dX the translation vector x.
+     * @param dY the translation vector y.
+     * @param dZ the translation vector z.
+     * @param rX the first rotation (radians).
+     * @param rY the second rotation (radians).
+     * @param rZ the third rotation (radians).
+     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).
+     * @return the new transformation.
+     */
+    static create(dX: float64, dY: float64, dZ: float64, rX: float64, rY: float64, rZ: float64, dS: float64): PositionVector;
+    /**
+     * Get the 7 parameters of the method.
+     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.
+     */
+    getParameters(): Float64Array;
+    /**
+     * Make a forward transformation.
+     * @param source the source coordinate (the coordinate to transform).
+     * @param target the target coordinate (holds the result).
+     */
+    transformForward(source: Coordinate, target: Coordinate): void;
+    /**
+     * Make a reverse transformation.
+     * @param source the source coordinate (holds the result).
+     * @param target the target coordinate (the coordinate to transform).
+     */
+    transformReverse(source: Coordinate, target: Coordinate): void;
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#forward
+     */
+    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#reverse
+     */
+    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+     * Create an inverse transformation.
+     * @return the inverse transformation.
+     */
+    createInverse(): PositionVector;
+}
+export {};
 //# sourceMappingURL=PositionVector.d.ts.map

package/lib/cjs/spatial/ecrs/transformation/PositionVector.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"PositionVector.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/transformation/PositionVector.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;AAErD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,qBAAa,cAAe,SAAQ,eAAe;IAC/C,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,6BAA6B;IAC7B,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,8BAA8B;IAC9B,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,yCAAyC;IACzC,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAYjD;;;;;;;;;;OAUG;WACW,MAAM,CAAC,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,cAAc;IAY/H;;;OAGG;IACI,aAAa,IAAI,YAAY;IAYpC;;;;OAIG;IACI,gBAAgB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAUrE;;;;OAIG;IACI,gBAAgB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAIrE;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAI5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAI5F;;;OAGG;IACI,aAAa,IAAI,cAAc;CAGzC"}
+{"version":3,"file":"PositionVector.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/transformation/PositionVector.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;AAErD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,qBAAa,cAAe,SAAQ,eAAe;IAC/C,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,6BAA6B;IAC7B,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,8BAA8B;IAC9B,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,yCAAyC;IACzC,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAYjD;;;;;;;;;;OAUG;WACW,MAAM,CAAC,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,OAAO,GAAG,cAAc;IAY/H;;;OAGG;IACI,aAAa,IAAI,YAAY;IAYpC;;;;OAIG;IACI,gBAAgB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAUrE;;;;OAIG;IACI,gBAAgB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAIrE;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAI5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAI5F;;;OAGG;IACI,aAAa,IAAI,cAAc;CAGzC"}

package/lib/cjs/spatial/ecrs/transformation/PositionVector.js

@@ -1,131 +1,131 @@
-"use strict";
-/*---------------------------------------------------------------------------------------------
-* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
-* See LICENSE.md in the project root for license terms and full copyright notice.
-*--------------------------------------------------------------------------------------------*/
-/** @packageDocumentation
- * @module OrbitGT
- */
-Object.defineProperty(exports, "__esModule", { value: true });
-exports.PositionVector = void 0;
-const OperationMethod_1 = require("../OperationMethod");
-const ParameterValue_1 = require("../ParameterValue");
-const ParameterValueList_1 = require("../ParameterValueList");
-/**
- * Class PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).
- *
- * NOTE: this method is also known as "Bursa-Wolf".
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised May 2005
- * Available at: http://www.epsg.org/
- *
- * Used in Europe. To convert to "Coordinate Frame Rotation" inverse the sign of the three rotations.
- *
- * @version 1.0 July 2005
- */
-/** @internal */
-class PositionVector extends OperationMethod_1.OperationMethod {
-    /**
-     * Create a new transformation.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters) {
-        super(PositionVector.METHOD_CODE, "Position Vector 7-param. transformation", parameters);
-        /* Store the parameters */
-        this._dX = parameters.getValue(8605);
-        this._dY = parameters.getValue(8606);
-        this._dZ = parameters.getValue(8607);
-        this._rX = parameters.getValue(8608);
-        this._rY = parameters.getValue(8609);
-        this._rZ = parameters.getValue(8610);
-        this._dS = parameters.getValue(8611);
-    }
-    /**
-     * Create a new transformation.
-     * @param dX the translation vector x.
-     * @param dY the translation vector y.
-     * @param dZ the translation vector z.
-     * @param rX the first rotation (radians).
-     * @param rY the second rotation (radians).
-     * @param rZ the third rotation (radians).
-     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).
-     * @return the new transformation.
-     */
-    static create(dX, dY, dZ, rX, rY, rZ, dS) {
-        let parameters = new ParameterValueList_1.ParameterValueList();
-        parameters.add(new ParameterValue_1.ParameterValue(8605, dX, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8606, dY, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8607, dZ, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8608, rX, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8609, rY, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8610, rZ, null));
-        parameters.add(new ParameterValue_1.ParameterValue(8611, dS, null));
-        return new PositionVector(parameters);
-    }
-    /**
-     * Get the 7 parameters of the method.
-     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.
-     */
-    getParameters() {
-        let parameters = new Float64Array(7);
-        parameters[0] = this._dX;
-        parameters[1] = this._dY;
-        parameters[2] = this._dZ;
-        parameters[3] = this._rX;
-        parameters[4] = this._rY;
-        parameters[5] = this._rZ;
-        parameters[6] = this._dS;
-        return parameters;
-    }
-    /**
-     * Make a forward transformation.
-     * @param source the source coordinate (the coordinate to transform).
-     * @param target the target coordinate (holds the result).
-     */
-    transformForward(source, target) {
-        let x = source.getX();
-        let y = source.getY();
-        let z = source.getZ();
-        let m = (1.0 + this._dS);
-        target.setX(m * (1.0 * x - this._rZ * y + this._rY * z) + this._dX);
-        target.setY(m * (this._rZ * x + 1.0 * y - this._rX * z) + this._dY);
-        target.setZ(m * (-this._rY * x + this._rX * y + 1.0 * z) + this._dZ);
-    }
-    /**
-     * Make a reverse transformation.
-     * @param source the source coordinate (holds the result).
-     * @param target the target coordinate (the coordinate to transform).
-     */
-    transformReverse(source, target) {
-        this.createInverse().transformForward(target, source);
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#forward
-     */
-    forward(sourceCRS, source, targetCRS, target) {
-        this.transformForward(source, target);
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#reverse
-     */
-    reverse(sourceCRS, source, targetCRS, target) {
-        this.transformReverse(source, target);
-    }
-    /**
-     * Create an inverse transformation.
-     * @return the inverse transformation.
-     */
-    createInverse() {
-        return PositionVector.create(-this._dX, -this._dY, -this._dZ, -this._rX, -this._rY, -this._rZ, -this._dS);
-    }
-}
-exports.PositionVector = PositionVector;
-/** The code of this method */
-PositionVector.METHOD_CODE = 9606;
+"use strict";
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.PositionVector = void 0;
+const OperationMethod_1 = require("../OperationMethod");
+const ParameterValue_1 = require("../ParameterValue");
+const ParameterValueList_1 = require("../ParameterValueList");
+/**
+ * Class PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).
+ *
+ * NOTE: this method is also known as "Bursa-Wolf".
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised May 2005
+ * Available at: http://www.epsg.org/
+ *
+ * Used in Europe. To convert to "Coordinate Frame Rotation" inverse the sign of the three rotations.
+ *
+ * @version 1.0 July 2005
+ */
+/** @internal */
+class PositionVector extends OperationMethod_1.OperationMethod {
+    /**
+     * Create a new transformation.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters) {
+        super(PositionVector.METHOD_CODE, "Position Vector 7-param. transformation", parameters);
+        /* Store the parameters */
+        this._dX = parameters.getValue(8605);
+        this._dY = parameters.getValue(8606);
+        this._dZ = parameters.getValue(8607);
+        this._rX = parameters.getValue(8608);
+        this._rY = parameters.getValue(8609);
+        this._rZ = parameters.getValue(8610);
+        this._dS = parameters.getValue(8611);
+    }
+    /**
+     * Create a new transformation.
+     * @param dX the translation vector x.
+     * @param dY the translation vector y.
+     * @param dZ the translation vector z.
+     * @param rX the first rotation (radians).
+     * @param rY the second rotation (radians).
+     * @param rZ the third rotation (radians).
+     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).
+     * @return the new transformation.
+     */
+    static create(dX, dY, dZ, rX, rY, rZ, dS) {
+        let parameters = new ParameterValueList_1.ParameterValueList();
+        parameters.add(new ParameterValue_1.ParameterValue(8605, dX, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8606, dY, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8607, dZ, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8608, rX, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8609, rY, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8610, rZ, null));
+        parameters.add(new ParameterValue_1.ParameterValue(8611, dS, null));
+        return new PositionVector(parameters);
+    }
+    /**
+     * Get the 7 parameters of the method.
+     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.
+     */
+    getParameters() {
+        let parameters = new Float64Array(7);
+        parameters[0] = this._dX;
+        parameters[1] = this._dY;
+        parameters[2] = this._dZ;
+        parameters[3] = this._rX;
+        parameters[4] = this._rY;
+        parameters[5] = this._rZ;
+        parameters[6] = this._dS;
+        return parameters;
+    }
+    /**
+     * Make a forward transformation.
+     * @param source the source coordinate (the coordinate to transform).
+     * @param target the target coordinate (holds the result).
+     */
+    transformForward(source, target) {
+        let x = source.getX();
+        let y = source.getY();
+        let z = source.getZ();
+        let m = (1.0 + this._dS);
+        target.setX(m * (1.0 * x - this._rZ * y + this._rY * z) + this._dX);
+        target.setY(m * (this._rZ * x + 1.0 * y - this._rX * z) + this._dY);
+        target.setZ(m * (-this._rY * x + this._rX * y + 1.0 * z) + this._dZ);
+    }
+    /**
+     * Make a reverse transformation.
+     * @param source the source coordinate (holds the result).
+     * @param target the target coordinate (the coordinate to transform).
+     */
+    transformReverse(source, target) {
+        this.createInverse().transformForward(target, source);
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#forward
+     */
+    forward(sourceCRS, source, targetCRS, target) {
+        this.transformForward(source, target);
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#reverse
+     */
+    reverse(sourceCRS, source, targetCRS, target) {
+        this.transformReverse(source, target);
+    }
+    /**
+     * Create an inverse transformation.
+     * @return the inverse transformation.
+     */
+    createInverse() {
+        return PositionVector.create(-this._dX, -this._dY, -this._dZ, -this._rX, -this._rY, -this._rZ, -this._dS);
+    }
+}
+/** The code of this method */
+PositionVector.METHOD_CODE = 9606;
+exports.PositionVector = PositionVector;
 //# sourceMappingURL=PositionVector.js.map

package/lib/cjs/spatial/ecrs/transformation/PositionVector.js.map

@@ -1 +1 @@
-{"version":3,"file":"PositionVector.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/transformation/PositionVector.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAYH,wDAAqD;AACrD,sDAAmD;AACnD,8DAA2D;AAE3D;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,MAAa,cAAe,SAAQ,iCAAe;IAe/C;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,cAAc,CAAC,WAAW,EAAE,yCAAyC,EAAE,UAAU,CAAC,CAAC;QACzF,0BAA0B;QAC1B,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;QACrC,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;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzC,CAAC;IAED;;;;;;;;;;OAUG;IACI,MAAM,CAAC,MAAM,CAAC,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QAC1G,IAAI,UAAU,GAAuB,IAAI,uCAAkB,EAAE,CAAC;QAC9D,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,OAAO,IAAI,cAAc,CAAC,UAAU,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,aAAa;QAChB,IAAI,UAAU,GAAiB,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,OAAO,UAAU,CAAC;IACtB,CAAC;IAED;;;;OAIG;IACI,gBAAgB,CAAC,MAAkB,EAAE,MAAkB;QAC1D,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAClC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACzE,CAAC;IAED;;;;OAIG;IACI,gBAAgB,CAAC,MAAkB,EAAE,MAAkB;QAC1D,IAAI,CAAC,aAAa,EAAE,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1D,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,aAAa;QAChB,OAAO,cAAc,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC9G,CAAC;;AApHL,wCAqHC;AApHG,8BAA8B;AACP,0BAAW,GAAU,IAAI,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module OrbitGT\r\n */\r\n\r\n//package orbitgt.spatial.ecrs.transformation;\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 { ParameterValue } from \"../ParameterValue\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).\r\n *\r\n * NOTE: this method is also known as \"Bursa-Wolf\".\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 * Used in Europe. To convert to \"Coordinate Frame Rotation\" inverse the sign of the three rotations.\r\n *\r\n * @version 1.0 July 2005\r\n */\r\n/** @internal */\r\nexport class PositionVector extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9606;\r\n\r\n    /** The translation vector */\r\n    private _dX: float64;\r\n    private _dY: float64;\r\n    private _dZ: float64;\r\n    /** The rotations (radians) */\r\n    private _rX: float64;\r\n    private _rY: float64;\r\n    private _rZ: float64;\r\n    /** The scale factor (offset from 1.0) */\r\n    private _dS: float64;\r\n\r\n    /**\r\n     * Create a new transformation.\r\n     * @param parameters the values of the parameters.\r\n     */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(PositionVector.METHOD_CODE, \"Position Vector 7-param. transformation\", parameters);\r\n        /* Store the parameters */\r\n        this._dX = parameters.getValue(8605);\r\n        this._dY = parameters.getValue(8606);\r\n        this._dZ = parameters.getValue(8607);\r\n        this._rX = parameters.getValue(8608);\r\n        this._rY = parameters.getValue(8609);\r\n        this._rZ = parameters.getValue(8610);\r\n        this._dS = parameters.getValue(8611);\r\n    }\r\n\r\n    /**\r\n     * Create a new transformation.\r\n     * @param dX the translation vector x.\r\n     * @param dY the translation vector y.\r\n     * @param dZ the translation vector z.\r\n     * @param rX the first rotation (radians).\r\n     * @param rY the second rotation (radians).\r\n     * @param rZ the third rotation (radians).\r\n     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).\r\n     * @return the new transformation.\r\n     */\r\n    public static create(dX: float64, dY: float64, dZ: float64, rX: float64, rY: float64, rZ: float64, dS: float64): PositionVector {\r\n        let parameters: ParameterValueList = new ParameterValueList();\r\n        parameters.add(new ParameterValue(8605, dX, null));\r\n        parameters.add(new ParameterValue(8606, dY, null));\r\n        parameters.add(new ParameterValue(8607, dZ, null));\r\n        parameters.add(new ParameterValue(8608, rX, null));\r\n        parameters.add(new ParameterValue(8609, rY, null));\r\n        parameters.add(new ParameterValue(8610, rZ, null));\r\n        parameters.add(new ParameterValue(8611, dS, null));\r\n        return new PositionVector(parameters);\r\n    }\r\n\r\n    /**\r\n     * Get the 7 parameters of the method.\r\n     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.\r\n     */\r\n    public getParameters(): Float64Array {\r\n        let parameters: Float64Array = new Float64Array(7);\r\n        parameters[0] = this._dX;\r\n        parameters[1] = this._dY;\r\n        parameters[2] = this._dZ;\r\n        parameters[3] = this._rX;\r\n        parameters[4] = this._rY;\r\n        parameters[5] = this._rZ;\r\n        parameters[6] = this._dS;\r\n        return parameters;\r\n    }\r\n\r\n    /**\r\n     * Make a forward transformation.\r\n     * @param source the source coordinate (the coordinate to transform).\r\n     * @param target the target coordinate (holds the result).\r\n     */\r\n    public transformForward(source: Coordinate, target: Coordinate): void {\r\n        let x: float64 = source.getX();\r\n        let y: float64 = source.getY();\r\n        let z: float64 = source.getZ();\r\n        let m: float64 = (1.0 + this._dS);\r\n        target.setX(m * (1.0 * x - this._rZ * y + this._rY * z) + this._dX);\r\n        target.setY(m * (this._rZ * x + 1.0 * y - this._rX * z) + this._dY);\r\n        target.setZ(m * (-this._rY * x + this._rX * y + 1.0 * z) + this._dZ);\r\n    }\r\n\r\n    /**\r\n     * Make a reverse transformation.\r\n     * @param source the source coordinate (holds the result).\r\n     * @param target the target coordinate (the coordinate to transform).\r\n     */\r\n    public transformReverse(source: Coordinate, target: Coordinate): void {\r\n        this.createInverse().transformForward(target, source);\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        this.transformForward(source, target);\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        this.transformReverse(source, target);\r\n    }\r\n\r\n    /**\r\n     * Create an inverse transformation.\r\n     * @return the inverse transformation.\r\n     */\r\n    public createInverse(): PositionVector {\r\n        return PositionVector.create(-this._dX, -this._dY, -this._dZ, -this._rX, -this._rY, -this._rZ, -this._dS);\r\n    }\r\n}\r\n"]}
+{"version":3,"file":"PositionVector.js","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/transformation/PositionVector.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;AAC/F;;GAEG;;;AAYH,wDAAqD;AACrD,sDAAmD;AACnD,8DAA2D;AAE3D;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,MAAa,cAAe,SAAQ,iCAAe;IAe/C;;;OAGG;IACH,YAAmB,UAA8B;QAC7C,KAAK,CAAC,cAAc,CAAC,WAAW,EAAE,yCAAyC,EAAE,UAAU,CAAC,CAAC;QACzF,0BAA0B;QAC1B,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;QACrC,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;QACrC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzC,CAAC;IAED;;;;;;;;;;OAUG;IACI,MAAM,CAAC,MAAM,CAAC,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW,EAAE,EAAW;QAC1G,IAAI,UAAU,GAAuB,IAAI,uCAAkB,EAAE,CAAC;QAC9D,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,GAAG,CAAC,IAAI,+BAAc,CAAC,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QACnD,OAAO,IAAI,cAAc,CAAC,UAAU,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,aAAa;QAChB,IAAI,UAAU,GAAiB,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACnD,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC;QACzB,OAAO,UAAU,CAAC;IACtB,CAAC;IAED;;;;OAIG;IACI,gBAAgB,CAAC,MAAkB,EAAE,MAAkB;QAC1D,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,MAAM,CAAC,IAAI,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAY,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QAClC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACpE,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACzE,CAAC;IAED;;;;OAIG;IACI,gBAAgB,CAAC,MAAkB,EAAE,MAAkB;QAC1D,IAAI,CAAC,aAAa,EAAE,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1D,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,OAAO,CAAC,SAAc,EAAE,MAAkB,EAAE,SAAc,EAAE,MAAkB;QACjF,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC1C,CAAC;IAED;;;OAGG;IACI,aAAa;QAChB,OAAO,cAAc,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC9G,CAAC;;AAnHD,8BAA8B;AACP,0BAAW,GAAU,IAAI,CAAC;AAFxC,wCAAc","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.transformation;\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 { ParameterValue } from \"../ParameterValue\";\r\nimport { ParameterValueList } from \"../ParameterValueList\";\r\n\r\n/**\r\n * Class PositionVector defines a Position Vector 7-parameter transformation (see Guidance Note 2.4.3.2.1).\r\n *\r\n * NOTE: this method is also known as \"Bursa-Wolf\".\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 * Used in Europe. To convert to \"Coordinate Frame Rotation\" inverse the sign of the three rotations.\r\n *\r\n * @version 1.0 July 2005\r\n */\r\n/** @internal */\r\nexport class PositionVector extends OperationMethod {\r\n    /** The code of this method */\r\n    public static readonly METHOD_CODE: int32 = 9606;\r\n\r\n    /** The translation vector */\r\n    private _dX: float64;\r\n    private _dY: float64;\r\n    private _dZ: float64;\r\n    /** The rotations (radians) */\r\n    private _rX: float64;\r\n    private _rY: float64;\r\n    private _rZ: float64;\r\n    /** The scale factor (offset from 1.0) */\r\n    private _dS: float64;\r\n\r\n    /**\r\n     * Create a new transformation.\r\n     * @param parameters the values of the parameters.\r\n     */\r\n    public constructor(parameters: ParameterValueList) {\r\n        super(PositionVector.METHOD_CODE, \"Position Vector 7-param. transformation\", parameters);\r\n        /* Store the parameters */\r\n        this._dX = parameters.getValue(8605);\r\n        this._dY = parameters.getValue(8606);\r\n        this._dZ = parameters.getValue(8607);\r\n        this._rX = parameters.getValue(8608);\r\n        this._rY = parameters.getValue(8609);\r\n        this._rZ = parameters.getValue(8610);\r\n        this._dS = parameters.getValue(8611);\r\n    }\r\n\r\n    /**\r\n     * Create a new transformation.\r\n     * @param dX the translation vector x.\r\n     * @param dY the translation vector y.\r\n     * @param dZ the translation vector z.\r\n     * @param rX the first rotation (radians).\r\n     * @param rY the second rotation (radians).\r\n     * @param rZ the third rotation (radians).\r\n     * @param dS the scale factor (offset from 1.0) (normal units, not ppm).\r\n     * @return the new transformation.\r\n     */\r\n    public static create(dX: float64, dY: float64, dZ: float64, rX: float64, rY: float64, rZ: float64, dS: float64): PositionVector {\r\n        let parameters: ParameterValueList = new ParameterValueList();\r\n        parameters.add(new ParameterValue(8605, dX, null));\r\n        parameters.add(new ParameterValue(8606, dY, null));\r\n        parameters.add(new ParameterValue(8607, dZ, null));\r\n        parameters.add(new ParameterValue(8608, rX, null));\r\n        parameters.add(new ParameterValue(8609, rY, null));\r\n        parameters.add(new ParameterValue(8610, rZ, null));\r\n        parameters.add(new ParameterValue(8611, dS, null));\r\n        return new PositionVector(parameters);\r\n    }\r\n\r\n    /**\r\n     * Get the 7 parameters of the method.\r\n     * @return the parameters of the method {dX,dY,dZ,rX,rY,rZ,dS} in units {meter,meter,meter, radian,radian,radian, scale_offset-1}.\r\n     */\r\n    public getParameters(): Float64Array {\r\n        let parameters: Float64Array = new Float64Array(7);\r\n        parameters[0] = this._dX;\r\n        parameters[1] = this._dY;\r\n        parameters[2] = this._dZ;\r\n        parameters[3] = this._rX;\r\n        parameters[4] = this._rY;\r\n        parameters[5] = this._rZ;\r\n        parameters[6] = this._dS;\r\n        return parameters;\r\n    }\r\n\r\n    /**\r\n     * Make a forward transformation.\r\n     * @param source the source coordinate (the coordinate to transform).\r\n     * @param target the target coordinate (holds the result).\r\n     */\r\n    public transformForward(source: Coordinate, target: Coordinate): void {\r\n        let x: float64 = source.getX();\r\n        let y: float64 = source.getY();\r\n        let z: float64 = source.getZ();\r\n        let m: float64 = (1.0 + this._dS);\r\n        target.setX(m * (1.0 * x - this._rZ * y + this._rY * z) + this._dX);\r\n        target.setY(m * (this._rZ * x + 1.0 * y - this._rX * z) + this._dY);\r\n        target.setZ(m * (-this._rY * x + this._rX * y + 1.0 * z) + this._dZ);\r\n    }\r\n\r\n    /**\r\n     * Make a reverse transformation.\r\n     * @param source the source coordinate (holds the result).\r\n     * @param target the target coordinate (the coordinate to transform).\r\n     */\r\n    public transformReverse(source: Coordinate, target: Coordinate): void {\r\n        this.createInverse().transformForward(target, source);\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        this.transformForward(source, target);\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        this.transformReverse(source, target);\r\n    }\r\n\r\n    /**\r\n     * Create an inverse transformation.\r\n     * @return the inverse transformation.\r\n     */\r\n    public createInverse(): PositionVector {\r\n        return PositionVector.create(-this._dX, -this._dY, -this._dZ, -this._rX, -this._rY, -this._rZ, -this._dS);\r\n    }\r\n}\r\n"]}