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

package/lib/cjs/spatial/ecrs/projection/TransverseMercator.d.ts

@@ -1,109 +1,109 @@
-/** @packageDocumentation
- * @module OrbitGT
- */
-declare type int32 = number;
-declare type float64 = number;
-import { Coordinate } from "../../geom/Coordinate";
-import { CRS } from "../CRS";
-import { Ellipsoid } from "../Ellipsoid";
-import { Operation } from "../Operation";
-import { OperationMethod } from "../OperationMethod";
-import { ParameterValueList } from "../ParameterValueList";
-/**
- * Class TransverseMercator defines a Transverse Mercator projection.
- *
- * 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 August 2006
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.4.5.1: General Case
- *
- * @version 1.0 January 2007
- */
-/** @internal */
-export declare class TransverseMercator extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** The value of PI */
-    private static readonly PI;
-    /** Latitude of natural origin (radians) */
-    private _latN;
-    /** Longitude of natural origin (the central meridian) (radians) */
-    private _lonN;
-    /** The scale factor at the natural origin */
-    private _k0;
-    /** False easting */
-    private _fE;
-    /** False northing */
-    private _fN;
-    /** Forward calculation parameters */
-    private _a;
-    private _e;
-    private _e2;
-    private _e4;
-    private _e6;
-    private _e_2;
-    private _m0;
-    /** Reverse calculation parameters */
-    private _e1;
-    private _R1;
-    private _R2;
-    private _R3;
-    private _R4;
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters: ParameterValueList);
-    /**
-       * Get the power of a number (must be able to handle negative 'n' values).
-       */
-    private static pow;
-    /**
-       * Calculate M.
-       */
-    private calcM;
-    /**
-       * Initialize the projection.
-       * @param ellipsoid the ellipsoid to use.
-       * @return this projection (for convenience).
-       */
-    initializeProjection(ellipsoid: Ellipsoid): TransverseMercator;
-    /**
-       * Do the projection.
-       * @param lon the longitude (radians).
-       * @param lat the latitude (radians).
-       * @param projected the target projected coordinate (X and Y will be set).
-       */
-    toProjection(lon: float64, lat: float64, projected: Coordinate): void;
-    /**
-       * Do the inverse projection.
-       * @param x the easting.
-       * @param y the northing.
-       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).
-       */
-    toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void;
-    /**
-       * OperationMethod method.
-       * @see OperationMethod#initialize
-       */
-    initialize(operation: Operation): void;
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#forward
-       */
-    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#reverse
-       */
-    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
-}
-export {};
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+type int32 = number;
+type float64 = number;
+import { Coordinate } from "../../geom/Coordinate";
+import { CRS } from "../CRS";
+import { Ellipsoid } from "../Ellipsoid";
+import { Operation } from "../Operation";
+import { OperationMethod } from "../OperationMethod";
+import { ParameterValueList } from "../ParameterValueList";
+/**
+ * Class TransverseMercator defines a Transverse Mercator projection.
+ *
+ * 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 August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.5.1: General Case
+ *
+ * @version 1.0 January 2007
+ */
+/** @internal */
+export declare class TransverseMercator extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** The value of PI */
+    private static readonly PI;
+    /** Latitude of natural origin (radians) */
+    private _latN;
+    /** Longitude of natural origin (the central meridian) (radians) */
+    private _lonN;
+    /** The scale factor at the natural origin */
+    private _k0;
+    /** False easting */
+    private _fE;
+    /** False northing */
+    private _fN;
+    /** Forward calculation parameters */
+    private _a;
+    private _e;
+    private _e2;
+    private _e4;
+    private _e6;
+    private _e_2;
+    private _m0;
+    /** Reverse calculation parameters */
+    private _e1;
+    private _R1;
+    private _R2;
+    private _R3;
+    private _R4;
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters: ParameterValueList);
+    /**
+       * Get the power of a number (must be able to handle negative 'n' values).
+       */
+    private static pow;
+    /**
+       * Calculate M.
+       */
+    private calcM;
+    /**
+       * Initialize the projection.
+       * @param ellipsoid the ellipsoid to use.
+       * @return this projection (for convenience).
+       */
+    initializeProjection(ellipsoid: Ellipsoid): TransverseMercator;
+    /**
+       * Do the projection.
+       * @param lon the longitude (radians).
+       * @param lat the latitude (radians).
+       * @param projected the target projected coordinate (X and Y will be set).
+       */
+    toProjection(lon: float64, lat: float64, projected: Coordinate): void;
+    /**
+       * Do the inverse projection.
+       * @param x the easting.
+       * @param y the northing.
+       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).
+       */
+    toGeoGraphic(x: float64, y: float64, geographic: Coordinate): void;
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation: Operation): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS: CRS, source: Coordinate, targetCRS: CRS, target: Coordinate): void;
+}
+export {};
 //# sourceMappingURL=TransverseMercator.d.ts.map

package/lib/cjs/spatial/ecrs/projection/TransverseMercator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"TransverseMercator.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/TransverseMercator.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,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,kBAAmB,SAAQ,eAAe;IACnD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,sBAAsB;IACtB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAoB;IAE9C,2CAA2C;IAC3C,OAAO,CAAC,KAAK,CAAU;IACvB,mEAAmE;IACnE,OAAO,CAAC,KAAK,CAAU;IACvB,6CAA6C;IAC7C,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB,qCAAqC;IACrC,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,GAAG,CAAU;IACrB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAUjD;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAMlB;;SAEK;IACL,OAAO,CAAC,KAAK;IAKb;;;;SAIK;IACE,oBAAoB,CAAC,SAAS,EAAE,SAAS,GAAG,kBAAkB;IAuBrE;;;;;SAKK;IACE,YAAY,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,SAAS,EAAE,UAAU,GAAG,IAAI;IAoB5E;;;;;SAKK;IACE,YAAY,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IA8BzE;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAItD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAS5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAQ/F"}
+{"version":3,"file":"TransverseMercator.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/TransverseMercator.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,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AACzC,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,kBAAmB,SAAQ,eAAe;IACnD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,sBAAsB;IACtB,OAAO,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAoB;IAE9C,2CAA2C;IAC3C,OAAO,CAAC,KAAK,CAAU;IACvB,mEAAmE;IACnE,OAAO,CAAC,KAAK,CAAU;IACvB,6CAA6C;IAC7C,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB,qCAAqC;IACrC,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,EAAE,CAAU;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,GAAG,CAAU;IACrB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;SAGK;gBACc,UAAU,EAAE,kBAAkB;IAUjD;;SAEK;IACL,OAAO,CAAC,MAAM,CAAC,GAAG;IAMlB;;SAEK;IACL,OAAO,CAAC,KAAK;IAKb;;;;SAIK;IACE,oBAAoB,CAAC,SAAS,EAAE,SAAS,GAAG,kBAAkB;IAuBrE;;;;;SAKK;IACE,YAAY,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,SAAS,EAAE,UAAU,GAAG,IAAI;IAoB5E;;;;;SAKK;IACE,YAAY,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IA8BzE;;;SAGK;IACW,UAAU,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAItD;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAS5F;;;SAGK;IACE,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAQ/F"}

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

@@ -1,184 +1,184 @@
-"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.TransverseMercator = void 0;
-const OperationMethod_1 = require("../OperationMethod");
-/**
- * Class TransverseMercator defines a Transverse Mercator projection.
- *
- * 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 August 2006
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.4.5.1: General Case
- *
- * @version 1.0 January 2007
- */
-/** @internal */
-class TransverseMercator extends OperationMethod_1.OperationMethod {
-    /**
-       * Create a new projection.
-       * @param parameters the values of the parameters.
-       */
-    constructor(parameters) {
-        super(TransverseMercator.METHOD_CODE, "Transverse Mercator", parameters);
-        /* Store the parameters */
-        this._latN = parameters.getValue(8801);
-        this._lonN = parameters.getValue(8802);
-        this._k0 = parameters.getValue(8805);
-        this._fE = parameters.getValue(8806);
-        this._fN = parameters.getValue(8807);
-    }
-    /**
-       * Get the power of a number (must be able to handle negative 'n' values).
-       */
-    static pow(n, e) {
-        let p = 1.0;
-        for (let i = 0; i < e; i++)
-            p *= n;
-        return p;
-    }
-    /**
-       * Calculate M.
-       */
-    calcM(lat) {
-        const M = 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);
-        return this._a * M;
-    }
-    /**
-       * Initialize the projection.
-       * @param ellipsoid the ellipsoid to use.
-       * @return this projection (for convenience).
-       */
-    initializeProjection(ellipsoid) {
-        /* Prepare the forward parameters */
-        this._a = ellipsoid.getA();
-        this._e = ellipsoid.getE();
-        this._e2 = this._e * this._e;
-        this._e4 = this._e2 * this._e2;
-        this._e6 = this._e2 * this._e4;
-        this._e_2 = this._e2 / (1.0 - this._e2);
-        this._m0 = this.calcM(this._latN);
-        /* Prepare the reverse parameters */
-        const temp = Math.sqrt(1 - this._e2);
-        this._e1 = (1.0 - temp) / (1.0 + temp);
-        const e1_2 = this._e1 * this._e1;
-        const e1_3 = this._e1 * e1_2;
-        const e1_4 = e1_2 * e1_2;
-        this._R1 = (3 * this._e1 / 2 - 27 * e1_3 / 32);
-        this._R2 = (21 * e1_2 / 16 - 55 * e1_4 / 32);
-        this._R3 = (151 * e1_3 / 96);
-        this._R4 = (1097 * e1_4 / 512);
-        /* Return the projection */
-        return this;
-    }
-    /**
-       * Do the projection.
-       * @param lon the longitude (radians).
-       * @param lat the latitude (radians).
-       * @param projected the target projected coordinate (X and Y will be set).
-       */
-    toProjection(lon, lat, projected) {
-        /* Calculate parameters */
-        const tan = Math.tan(lat);
-        const T = tan * tan;
-        const cos = Math.cos(lat);
-        const C = this._e2 * cos * cos / (1.0 - this._e2);
-        const A = (lon - this._lonN) * cos;
-        const eSin = this._e * Math.sin(lat);
-        const v = this._a / Math.sqrt(1.0 - eSin * eSin);
-        const m = this.calcM(lat);
-        /* Make the calculation */
-        const temp1 = 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;
-        const E = this._fE + this._k0 * v * (temp1);
-        const temp2 = 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);
-        const N = this._fN + this._k0 * (temp2);
-        /* Save the position */
-        projected.setX(E);
-        projected.setY(N);
-    }
-    /**
-       * Do the inverse projection.
-       * @param x the easting.
-       * @param y the northing.
-       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).
-       */
-    toGeoGraphic(x, y, geographic) {
-        /* Get the parameters */
-        const E = x;
-        const N = y;
-        /* Calculate parameters */
-        const m1 = this._m0 + (N - this._fN) / this._k0;
-        const u1 = m1 / (this._a * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256));
-        const lat1 = 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);
-        const eSin1 = this._e * Math.sin(lat1);
-        const temp = Math.sqrt(1.0 - eSin1 * eSin1);
-        const v1 = this._a / temp;
-        const rho1 = this._a * (1.0 - this._e2) / TransverseMercator.pow(temp, 3);
-        const tan = Math.tan(lat1);
-        const T1 = tan * tan;
-        const T1_2 = T1 * T1;
-        const cos = Math.cos(lat1);
-        const C1 = this._e_2 * cos * cos;
-        const C1_2 = C1 * C1;
-        const D = (E - this._fE) / (v1 * this._k0);
-        const D_2 = D * D;
-        /* Make the calculation */
-        const temp1 = 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;
-        const lat = lat1 - (v1 * tan / rho1) * (temp1);
-        const temp2 = 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;
-        const lon = this._lonN + (temp2) / cos;
-        /* Save the position */
-        geographic.setX(lon);
-        geographic.setY(lat);
-    }
-    /**
-       * OperationMethod method.
-       * @see OperationMethod#initialize
-       */
-    initialize(operation) {
-        this.initializeProjection(operation.getSourceCRS().getEllipsoid()); // should this use getTargetCRS? (while testing OSTN15 RoI)
-    }
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#forward
-       */
-    forward(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        const lon = source.getX();
-        const lat = source.getY();
-        /* Do the projection */
-        this.toProjection(lon, lat, target);
-        target.setZ(source.getZ()); // Copy the Z value
-    }
-    /**
-       * OperationMethod interface method.
-       * @see OperationMethod#reverse
-       */
-    reverse(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        const E = target.getX();
-        const N = target.getY();
-        /* Do the inverse projection */
-        this.toGeoGraphic(E, N, source);
-        source.setZ(target.getZ()); // Copy the Z value
-    }
-}
-exports.TransverseMercator = TransverseMercator;
-/** The code of this method */
-TransverseMercator.METHOD_CODE = 9807;
-/** The value of PI */
-TransverseMercator.PI = Math.PI;
+"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.TransverseMercator = void 0;
+const OperationMethod_1 = require("../OperationMethod");
+/**
+ * Class TransverseMercator defines a Transverse Mercator projection.
+ *
+ * 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 August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.5.1: General Case
+ *
+ * @version 1.0 January 2007
+ */
+/** @internal */
+class TransverseMercator extends OperationMethod_1.OperationMethod {
+    /**
+       * Create a new projection.
+       * @param parameters the values of the parameters.
+       */
+    constructor(parameters) {
+        super(TransverseMercator.METHOD_CODE, "Transverse Mercator", parameters);
+        /* Store the parameters */
+        this._latN = parameters.getValue(8801);
+        this._lonN = parameters.getValue(8802);
+        this._k0 = parameters.getValue(8805);
+        this._fE = parameters.getValue(8806);
+        this._fN = parameters.getValue(8807);
+    }
+    /**
+       * Get the power of a number (must be able to handle negative 'n' values).
+       */
+    static pow(n, e) {
+        let p = 1.0;
+        for (let i = 0; i < e; i++)
+            p *= n;
+        return p;
+    }
+    /**
+       * Calculate M.
+       */
+    calcM(lat) {
+        const M = 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);
+        return this._a * M;
+    }
+    /**
+       * Initialize the projection.
+       * @param ellipsoid the ellipsoid to use.
+       * @return this projection (for convenience).
+       */
+    initializeProjection(ellipsoid) {
+        /* Prepare the forward parameters */
+        this._a = ellipsoid.getA();
+        this._e = ellipsoid.getE();
+        this._e2 = this._e * this._e;
+        this._e4 = this._e2 * this._e2;
+        this._e6 = this._e2 * this._e4;
+        this._e_2 = this._e2 / (1.0 - this._e2);
+        this._m0 = this.calcM(this._latN);
+        /* Prepare the reverse parameters */
+        const temp = Math.sqrt(1 - this._e2);
+        this._e1 = (1.0 - temp) / (1.0 + temp);
+        const e1_2 = this._e1 * this._e1;
+        const e1_3 = this._e1 * e1_2;
+        const e1_4 = e1_2 * e1_2;
+        this._R1 = (3 * this._e1 / 2 - 27 * e1_3 / 32);
+        this._R2 = (21 * e1_2 / 16 - 55 * e1_4 / 32);
+        this._R3 = (151 * e1_3 / 96);
+        this._R4 = (1097 * e1_4 / 512);
+        /* Return the projection */
+        return this;
+    }
+    /**
+       * Do the projection.
+       * @param lon the longitude (radians).
+       * @param lat the latitude (radians).
+       * @param projected the target projected coordinate (X and Y will be set).
+       */
+    toProjection(lon, lat, projected) {
+        /* Calculate parameters */
+        const tan = Math.tan(lat);
+        const T = tan * tan;
+        const cos = Math.cos(lat);
+        const C = this._e2 * cos * cos / (1.0 - this._e2);
+        const A = (lon - this._lonN) * cos;
+        const eSin = this._e * Math.sin(lat);
+        const v = this._a / Math.sqrt(1.0 - eSin * eSin);
+        const m = this.calcM(lat);
+        /* Make the calculation */
+        const temp1 = 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;
+        const E = this._fE + this._k0 * v * (temp1);
+        const temp2 = 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);
+        const N = this._fN + this._k0 * (temp2);
+        /* Save the position */
+        projected.setX(E);
+        projected.setY(N);
+    }
+    /**
+       * Do the inverse projection.
+       * @param x the easting.
+       * @param y the northing.
+       * @param geographic the target geographic coordinate (X/Lon and Y/Lat will be set) (radians).
+       */
+    toGeoGraphic(x, y, geographic) {
+        /* Get the parameters */
+        const E = x;
+        const N = y;
+        /* Calculate parameters */
+        const m1 = this._m0 + (N - this._fN) / this._k0;
+        const u1 = m1 / (this._a * (1.0 - this._e2 / 4 - 3 * this._e4 / 64 - 5 * this._e6 / 256));
+        const lat1 = 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);
+        const eSin1 = this._e * Math.sin(lat1);
+        const temp = Math.sqrt(1.0 - eSin1 * eSin1);
+        const v1 = this._a / temp;
+        const rho1 = this._a * (1.0 - this._e2) / TransverseMercator.pow(temp, 3);
+        const tan = Math.tan(lat1);
+        const T1 = tan * tan;
+        const T1_2 = T1 * T1;
+        const cos = Math.cos(lat1);
+        const C1 = this._e_2 * cos * cos;
+        const C1_2 = C1 * C1;
+        const D = (E - this._fE) / (v1 * this._k0);
+        const D_2 = D * D;
+        /* Make the calculation */
+        const temp1 = 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;
+        const lat = lat1 - (v1 * tan / rho1) * (temp1);
+        const temp2 = 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;
+        const lon = this._lonN + (temp2) / cos;
+        /* Save the position */
+        geographic.setX(lon);
+        geographic.setY(lat);
+    }
+    /**
+       * OperationMethod method.
+       * @see OperationMethod#initialize
+       */
+    initialize(operation) {
+        this.initializeProjection(operation.getSourceCRS().getEllipsoid()); // should this use getTargetCRS? (while testing OSTN15 RoI)
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#forward
+       */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const lon = source.getX();
+        const lat = source.getY();
+        /* Do the projection */
+        this.toProjection(lon, lat, target);
+        target.setZ(source.getZ()); // Copy the Z value
+    }
+    /**
+       * OperationMethod interface method.
+       * @see OperationMethod#reverse
+       */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        const E = target.getX();
+        const N = target.getY();
+        /* Do the inverse projection */
+        this.toGeoGraphic(E, N, source);
+        source.setZ(target.getZ()); // Copy the Z value
+    }
+}
+/** The code of this method */
+TransverseMercator.METHOD_CODE = 9807;
+/** The value of PI */
+TransverseMercator.PI = Math.PI;
+exports.TransverseMercator = TransverseMercator;
 //# sourceMappingURL=TransverseMercator.js.map