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

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

@@ -1,73 +1,73 @@
-/** @packageDocumentation
- * @module OrbitGT
- */
-declare type int32 = number;
-import { Coordinate } from "../../geom/Coordinate";
-import { CRS } from "../CRS";
-import { OperationMethod } from "../OperationMethod";
-import { ParameterValueList } from "../ParameterValueList";
-/**
- * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.
- *
- * The 'source' CRS is the geographic CRS.
- * The 'target' CRS is the projected CRS.
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised August 2006
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.4.1.1 -> 1.4.1.4
- *
- * @version 1.0 January 2007
- */
-/** @internal */
-export declare class LambertConical1SP extends OperationMethod {
-    /** The code of this method */
-    static readonly METHOD_CODE: int32;
-    /** Latitude of natural origin */
-    private _latN;
-    /** Longitude of natural origin (the central meridian) */
-    private _lonN;
-    /** Scale factor at natural origin */
-    private _k0;
-    /** False easting */
-    private _fE;
-    /** False northing */
-    private _fN;
-    /**
-     * Create a new projection.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters: ParameterValueList);
-    /**
-     * Get the sign of a number.
-     */
-    private static sign;
-    /**
-     * Calculate M.
-     */
-    private static calcM;
-    /**
-     * Calculate T.
-     */
-    private static calcT;
-    /**
-     * Calculate R.
-     */
-    private static calcR;
-    /**
-     * 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;
+import { Coordinate } from "../../geom/Coordinate";
+import { CRS } from "../CRS";
+import { OperationMethod } from "../OperationMethod";
+import { ParameterValueList } from "../ParameterValueList";
+/**
+ * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.1.1 -> 1.4.1.4
+ *
+ * @version 1.0 January 2007
+ */
+/** @internal */
+export declare class LambertConical1SP extends OperationMethod {
+    /** The code of this method */
+    static readonly METHOD_CODE: int32;
+    /** Latitude of natural origin */
+    private _latN;
+    /** Longitude of natural origin (the central meridian) */
+    private _lonN;
+    /** Scale factor at natural origin */
+    private _k0;
+    /** False easting */
+    private _fE;
+    /** False northing */
+    private _fN;
+    /**
+     * Create a new projection.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters: ParameterValueList);
+    /**
+     * Get the sign of a number.
+     */
+    private static sign;
+    /**
+     * Calculate M.
+     */
+    private static calcM;
+    /**
+     * Calculate T.
+     */
+    private static calcT;
+    /**
+     * Calculate R.
+     */
+    private static calcR;
+    /**
+     * 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=LambertConical1SP.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"LambertConical1SP.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical1SP.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,aAAK,KAAK,GAAG,MAAM,CAAC;AAIpB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,iBAAkB,SAAQ,eAAe;IAClD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,iCAAiC;IACjC,OAAO,CAAC,KAAK,CAAU;IACvB,yDAAyD;IACzD,OAAO,CAAC,KAAK,CAAU;IACvB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAUjD;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,IAAI;IAInB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAIpB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAIpB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAKpB;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAwB5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAwB/F"}
+{"version":3,"file":"LambertConical1SP.d.ts","sourceRoot":"","sources":["../../../../../src/spatial/ecrs/projection/LambertConical1SP.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,KAAK,KAAK,GAAG,MAAM,CAAC;AAIpB,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,GAAG,EAAE,MAAM,QAAQ,CAAC;AAC7B,OAAO,EAAE,eAAe,EAAE,MAAM,oBAAoB,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAE3D;;;;;;;;;;;;;;;;GAgBG;AACH,gBAAgB;AAChB,qBAAa,iBAAkB,SAAQ,eAAe;IAClD,8BAA8B;IAC9B,gBAAuB,WAAW,EAAE,KAAK,CAAQ;IAEjD,iCAAiC;IACjC,OAAO,CAAC,KAAK,CAAU;IACvB,yDAAyD;IACzD,OAAO,CAAC,KAAK,CAAU;IACvB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IACrB,oBAAoB;IACpB,OAAO,CAAC,GAAG,CAAU;IACrB,qBAAqB;IACrB,OAAO,CAAC,GAAG,CAAU;IAErB;;;OAGG;gBACgB,UAAU,EAAE,kBAAkB;IAUjD;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,IAAI;IAInB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAIpB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAIpB;;OAEG;IACH,OAAO,CAAC,MAAM,CAAC,KAAK;IAKpB;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;IAwB5F;;;OAGG;IACI,OAAO,CAAC,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,EAAE,SAAS,EAAE,GAAG,EAAE,MAAM,EAAE,UAAU,GAAG,IAAI;CAwB/F"}

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

@@ -1,130 +1,130 @@
-"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.LambertConical1SP = void 0;
-const OperationMethod_1 = require("../OperationMethod");
-/**
- * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.
- *
- * The 'source' CRS is the geographic CRS.
- * The 'target' CRS is the projected CRS.
- *
- * Based on the following document:
- *
- * Coordinate Conversions and Transformations including Formulas
- * Guidance Note Number 7, part 2
- * Revised August 2006
- * Available at: http://www.epsg.org/
- *
- * Formulas: see 1.4.1.1 -> 1.4.1.4
- *
- * @version 1.0 January 2007
- */
-/** @internal */
-class LambertConical1SP extends OperationMethod_1.OperationMethod {
-    /**
-     * Create a new projection.
-     * @param parameters the values of the parameters.
-     */
-    constructor(parameters) {
-        super(LambertConical1SP.METHOD_CODE, "Lambert Conic Conformal (1SP)", 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 sign of a number.
-     */
-    static sign(v) {
-        return (v < 0.0) ? (-1.0) : (1.0);
-    }
-    /**
-     * Calculate M.
-     */
-    static calcM(e, lat) {
-        return Math.cos(lat) / Math.pow(1.0 - Math.pow(e * Math.sin(lat), 2.0), 0.5);
-    }
-    /**
-     * Calculate T.
-     */
-    static calcT(e, lat) {
-        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);
-    }
-    /**
-     * Calculate R.
-     */
-    static calcR(a, F, t, n, k0) {
-        if (Math.abs(t) < 1.0e-6)
-            return 0.0;
-        return a * F * Math.pow(t, n) * k0;
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#forward
-     */
-    forward(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        let lon = source.getX();
-        let lat = source.getY();
-        /* Get the ellipsoid parameters */
-        let e = sourceCRS.getEllipsoid().getE();
-        let a = sourceCRS.getEllipsoid().getA();
-        /* Make the calculation */
-        let m0 = LambertConical1SP.calcM(e, this._latN);
-        let t0 = LambertConical1SP.calcT(e, this._latN);
-        let n = Math.sin(this._latN);
-        let F = m0 / (n * Math.pow(t0, n));
-        let r0 = LambertConical1SP.calcR(a, F, t0, n, this._k0);
-        let t = LambertConical1SP.calcT(e, lat);
-        let r = LambertConical1SP.calcR(a, F, t, n, this._k0);
-        let theta = n * (lon - this._lonN);
-        let E = this._fE + r * Math.sin(theta);
-        let N = this._fN + r0 - r * Math.cos(theta);
-        /* Save the position */
-        target.setX(E);
-        target.setY(N);
-        target.setZ(source.getZ()); // Keep the Z value
-    }
-    /**
-     * OperationMethod interface method.
-     * @see OperationMethod#reverse
-     */
-    reverse(sourceCRS, source, targetCRS, target) {
-        /* Get the parameters */
-        let E = target.getX();
-        let N = target.getY();
-        /* Get the ellipsoid parameters */
-        let e = sourceCRS.getEllipsoid().getE();
-        let a = sourceCRS.getEllipsoid().getA();
-        /* Make the calculation */
-        let m0 = LambertConical1SP.calcM(e, this._latN);
-        let t0 = LambertConical1SP.calcT(e, this._latN);
-        let n = Math.sin(this._latN);
-        let F = m0 / (n * Math.pow(t0, n));
-        let r0 = LambertConical1SP.calcR(a, F, t0, n, this._k0);
-        let r_ = LambertConical1SP.sign(n) * Math.pow(Math.pow(E - this._fE, 2.0) + Math.pow(r0 - (N - this._fN), 2.0), 0.5);
-        let t_ = Math.pow(r_ / (a * this._k0 * F), 1.0 / n);
-        let theta_ = Math.atan((E - this._fE) / (r0 - (N - this._fN)));
-        let lat = Math.PI * 0.5 - 2.0 * Math.atan(t_);
-        for (let i = 0; i < 7; i++)
-            lat = Math.PI * 0.5 - 2.0 * Math.atan(t_ * Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5)); // recursive formula
-        let lon = (theta_) / n + this._lonN;
-        /* Save the position */
-        source.setX(lon);
-        source.setY(lat);
-        source.setZ(target.getZ()); // Keep the Z value
-    }
-}
-exports.LambertConical1SP = LambertConical1SP;
-/** The code of this method */
-LambertConical1SP.METHOD_CODE = 9801;
+"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.LambertConical1SP = void 0;
+const OperationMethod_1 = require("../OperationMethod");
+/**
+ * Class LambertConical1SP defines a Lambert Conical Conformal map projection with 1 standard parallel.
+ *
+ * The 'source' CRS is the geographic CRS.
+ * The 'target' CRS is the projected CRS.
+ *
+ * Based on the following document:
+ *
+ * Coordinate Conversions and Transformations including Formulas
+ * Guidance Note Number 7, part 2
+ * Revised August 2006
+ * Available at: http://www.epsg.org/
+ *
+ * Formulas: see 1.4.1.1 -> 1.4.1.4
+ *
+ * @version 1.0 January 2007
+ */
+/** @internal */
+class LambertConical1SP extends OperationMethod_1.OperationMethod {
+    /**
+     * Create a new projection.
+     * @param parameters the values of the parameters.
+     */
+    constructor(parameters) {
+        super(LambertConical1SP.METHOD_CODE, "Lambert Conic Conformal (1SP)", 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 sign of a number.
+     */
+    static sign(v) {
+        return (v < 0.0) ? (-1.0) : (1.0);
+    }
+    /**
+     * Calculate M.
+     */
+    static calcM(e, lat) {
+        return Math.cos(lat) / Math.pow(1.0 - Math.pow(e * Math.sin(lat), 2.0), 0.5);
+    }
+    /**
+     * Calculate T.
+     */
+    static calcT(e, lat) {
+        return Math.tan(Math.PI * 0.25 - lat * 0.5) / Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5);
+    }
+    /**
+     * Calculate R.
+     */
+    static calcR(a, F, t, n, k0) {
+        if (Math.abs(t) < 1.0e-6)
+            return 0.0;
+        return a * F * Math.pow(t, n) * k0;
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#forward
+     */
+    forward(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        let lon = source.getX();
+        let lat = source.getY();
+        /* Get the ellipsoid parameters */
+        let e = sourceCRS.getEllipsoid().getE();
+        let a = sourceCRS.getEllipsoid().getA();
+        /* Make the calculation */
+        let m0 = LambertConical1SP.calcM(e, this._latN);
+        let t0 = LambertConical1SP.calcT(e, this._latN);
+        let n = Math.sin(this._latN);
+        let F = m0 / (n * Math.pow(t0, n));
+        let r0 = LambertConical1SP.calcR(a, F, t0, n, this._k0);
+        let t = LambertConical1SP.calcT(e, lat);
+        let r = LambertConical1SP.calcR(a, F, t, n, this._k0);
+        let theta = n * (lon - this._lonN);
+        let E = this._fE + r * Math.sin(theta);
+        let N = this._fN + r0 - r * Math.cos(theta);
+        /* Save the position */
+        target.setX(E);
+        target.setY(N);
+        target.setZ(source.getZ()); // Keep the Z value
+    }
+    /**
+     * OperationMethod interface method.
+     * @see OperationMethod#reverse
+     */
+    reverse(sourceCRS, source, targetCRS, target) {
+        /* Get the parameters */
+        let E = target.getX();
+        let N = target.getY();
+        /* Get the ellipsoid parameters */
+        let e = sourceCRS.getEllipsoid().getE();
+        let a = sourceCRS.getEllipsoid().getA();
+        /* Make the calculation */
+        let m0 = LambertConical1SP.calcM(e, this._latN);
+        let t0 = LambertConical1SP.calcT(e, this._latN);
+        let n = Math.sin(this._latN);
+        let F = m0 / (n * Math.pow(t0, n));
+        let r0 = LambertConical1SP.calcR(a, F, t0, n, this._k0);
+        let r_ = LambertConical1SP.sign(n) * Math.pow(Math.pow(E - this._fE, 2.0) + Math.pow(r0 - (N - this._fN), 2.0), 0.5);
+        let t_ = Math.pow(r_ / (a * this._k0 * F), 1.0 / n);
+        let theta_ = Math.atan((E - this._fE) / (r0 - (N - this._fN)));
+        let lat = Math.PI * 0.5 - 2.0 * Math.atan(t_);
+        for (let i = 0; i < 7; i++)
+            lat = Math.PI * 0.5 - 2.0 * Math.atan(t_ * Math.pow((1.0 - e * Math.sin(lat)) / (1.0 + e * Math.sin(lat)), e * 0.5)); // recursive formula
+        let lon = (theta_) / n + this._lonN;
+        /* Save the position */
+        source.setX(lon);
+        source.setY(lat);
+        source.setZ(target.getZ()); // Keep the Z value
+    }
+}
+/** The code of this method */
+LambertConical1SP.METHOD_CODE = 9801;
+exports.LambertConical1SP = LambertConical1SP;
 //# sourceMappingURL=LambertConical1SP.js.map

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

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