@itwin/core-orbitgt 3.0.0-dev.100

package/lib/esm/spatial/ecrs/Ellipsoid.d.ts

@@ -0,0 +1,189 @@
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+declare type int32 = number;
+declare type float64 = number;
+import { Coordinate } from "../geom/Coordinate";
+/**
+ * Class Ellipsoid defines the parameters of an earth ellipsoid.
+ *
+ * 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/
+ *
+ * Geocentric coordinates are defined as follows:
+ * The point (0,0,0) denotes the center of the ellipsoid. The z-axis is defined as being parallel to the earth rotational axis, pointing towards north.
+ * The x-axis intersects the ellipsoid at the 0 deg latitude, 0 deg longitude point.
+ *
+ * @version 1.0 July 2005
+ */
+/** @internal */
+export declare class Ellipsoid {
+    /** The code */
+    private _code;
+    /** The name */
+    private _name;
+    /** The code of the unit of measure */
+    private _unitCode;
+    /** Semi-major axis (meter) */
+    private _a;
+    /** Semi-minor axis (meter) (derived) */
+    private _b;
+    /** Flattening (derived) */
+    private _f;
+    /** Inverse flattening */
+    private _invF;
+    /** Eccentricity (derived) */
+    private _e;
+    /** Eccentricity squared (derived) */
+    private _e2;
+    /**
+     * Create a new ellipsoid.
+     * @param code the code.
+     * @param name the name.
+     * @param unitCode the code of the unit of measure.
+     * @param a the semi-major axis.
+     * @param invF the inverse flattening (value like 300, not like 1/300).
+     * @param b the semi-minor axis.
+     */
+    constructor(code: int32, name: string, unitCode: int32, a: float64, invF: float64, b: float64);
+    /**
+     * Get the code.
+     * @return the code.
+     */
+    getCode(): int32;
+    /**
+     * Get the name.
+     * @return the name.
+     */
+    getName(): string;
+    /**
+     * Get the code of the unit of measure.
+     * @return the code of the unit of measure.
+     */
+    getUnitCode(): int32;
+    /**
+     * Get the semi-major axis (in meter).
+     * @return the semi-major axis.
+     */
+    getA(): float64;
+    /**
+     * Get the semi-minor axis (in meter).
+     * @return the semi-minor axis.
+     */
+    getB(): float64;
+    /**
+     * Get the flattening.
+     * @return the flattening.
+     */
+    getF(): float64;
+    /**
+     * Get the inverse flattening.
+     * @return the inverse flattening.
+     */
+    getInvF(): float64;
+    /**
+     * Get the eccentricity.
+     * @return the eccentricity.
+     */
+    getE(): float64;
+    /**
+     * Get the radius of curvature of the ellipsoid in the plane of the meridian at a given latitude (how latitude curves).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getMeridianRadius(lat: float64): float64;
+    /**
+     * Get the radius of the circle defined by all points at a certain latitude (a circle of latitude, also called a parallel) (a at latitude 0, 0 at latitude 90) (how longitude curves).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getParallelRadius(lat: float64): float64;
+    /**
+     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getPrimeVerticalRadius(lat: float64): float64;
+    /**
+     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).
+     * @param lat the latitude (in radians).
+     * @param sinLat the sinus of the latitude.
+     * @return a radius (in meter).
+     */
+    private getPrimeVerticalRadius2;
+    /**
+     * Get the number of meter per radian of longitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per radian.
+     */
+    getMeterPerRadOfLon(lat: float64): float64;
+    /**
+     * Get the number of meter per degree of longitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per degree.
+     */
+    getMeterPerDegreeOfLon(lat: float64): float64;
+    /**
+     * Get the number of meter per radian of latitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per radian.
+     */
+    getMeterPerRadOfLat(lat: float64): float64;
+    /**
+     * Get the number of meter per degree of latitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per degree.
+     */
+    getMeterPerDegreeOfLat(lat: float64): float64;
+    /**
+     * Calculate the secant.
+     * @param v an angle (in radians).
+     * @return the secant.
+     */
+    private static sec;
+    /**
+     * Convert 3D geographic coordinates to 3D geocentric coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geographic the geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter).
+     * @param geocentric the new geocentric coordinates (in meter) (the result object).
+     */
+    toGeoCentric(geographic: Coordinate, geocentric: Coordinate): void;
+    /**
+     * Convert 3D geographic coordinates to 3D geocentric coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geographic the geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter).
+     * @param geocentric the new geocentric coordinates (in meter) (the result object).
+     */
+    toGeoCentricDeg(geographic: Coordinate, geocentric: Coordinate): void;
+    /**
+     * Convert 3D geocentric coordinates to 3D geographic coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geocentric the geocentric coordinates (in meter).
+     * @param geographic the new geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter) (the result object).
+     */
+    toGeoGraphic(geocentric: Coordinate, geographic: Coordinate): void;
+    /**
+     * Convert 3D geocentric coordinates to 3D geographic coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geocentric the geocentric coordinates (in meter).
+     * @param geographic the new geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter) (the result object).
+     */
+    toGeoGraphicDeg(geocentric: Coordinate, geographic: Coordinate): void;
+    /**
+     * Check if another ellipsoid is compatible with this one.
+     * @param other the other ellipsoid.
+     * @return true if compatible.
+     */
+    isCompatible(other: Ellipsoid): boolean;
+    /**
+     * The standard toString method.
+     * @see Object#toString
+     */
+    toString(): string;
+}
+export {};
+//# sourceMappingURL=Ellipsoid.d.ts.map

package/lib/esm/spatial/ecrs/Ellipsoid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"Ellipsoid.d.ts","sourceRoot":"","sources":["../../../../src/spatial/ecrs/Ellipsoid.ts"],"names":[],"mappings":"AAIA;;GAEG;AAMH,aAAK,KAAK,GAAG,MAAM,CAAC;AAEpB,aAAK,OAAO,GAAG,MAAM,CAAC;AAEtB,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAIhD;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,qBAAa,SAAS;IAClB,eAAe;IACf,OAAO,CAAC,KAAK,CAAQ;IACrB,eAAe;IACf,OAAO,CAAC,KAAK,CAAS;IACtB,sCAAsC;IACtC,OAAO,CAAC,SAAS,CAAQ;IACzB,8BAA8B;IAC9B,OAAO,CAAC,EAAE,CAAU;IACpB,wCAAwC;IACxC,OAAO,CAAC,EAAE,CAAU;IACpB,2BAA2B;IAC3B,OAAO,CAAC,EAAE,CAAU;IACpB,yBAAyB;IACzB,OAAO,CAAC,KAAK,CAAU;IACvB,6BAA6B;IAC7B,OAAO,CAAC,EAAE,CAAU;IACpB,qCAAqC;IACrC,OAAO,CAAC,GAAG,CAAU;IAErB;;;;;;;;OAQG;gBACgB,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,EAAE,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO;IAgBpG;;;OAGG;IACI,OAAO,IAAI,KAAK;IAIvB;;;OAGG;IACI,OAAO,IAAI,MAAM;IAIxB;;;OAGG;IACI,WAAW,IAAI,KAAK;IAI3B;;;OAGG;IACI,IAAI,IAAI,OAAO;IAItB;;;OAGG;IACI,IAAI,IAAI,OAAO;IAItB;;;OAGG;IACI,IAAI,IAAI,OAAO;IAItB;;;OAGG;IACI,OAAO,IAAI,OAAO;IAIzB;;;OAGG;IACI,IAAI,IAAI,OAAO;IAItB;;;;OAIG;IACI,iBAAiB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAM/C;;;;OAIG;IACI,iBAAiB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAM/C;;;;OAIG;IACI,sBAAsB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAKpD;;;;;OAKG;IACH,OAAO,CAAC,uBAAuB;IAK/B;;;;OAIG;IACI,mBAAmB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAIjD;;;;OAIG;IACI,sBAAsB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAIpD;;;;OAIG;IACI,mBAAmB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAIjD;;;;OAIG;IACI,sBAAsB,CAAC,GAAG,EAAE,OAAO,GAAG,OAAO;IAIpD;;;;OAIG;IACH,OAAO,CAAC,MAAM,CAAC,GAAG;IAIlB;;;;;OAKG;IACI,YAAY,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IAkBzE;;;;;OAKG;IACI,eAAe,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IAkB5E;;;;;OAKG;IACI,YAAY,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IAyCzE;;;;;OAKG;IACI,eAAe,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,EAAE,UAAU,GAAG,IAAI;IAuB5E;;;;OAIG;IACI,YAAY,CAAC,KAAK,EAAE,SAAS,GAAG,OAAO;IAQ9C;;;OAGG;IACI,QAAQ,IAAI,MAAM;CAG5B"}

package/lib/esm/spatial/ecrs/Ellipsoid.js

@@ -0,0 +1,333 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+import { Registry } from "./Registry";
+import { Unit } from "./Unit";
+/**
+ * Class Ellipsoid defines the parameters of an earth ellipsoid.
+ *
+ * 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/
+ *
+ * Geocentric coordinates are defined as follows:
+ * The point (0,0,0) denotes the center of the ellipsoid. The z-axis is defined as being parallel to the earth rotational axis, pointing towards north.
+ * The x-axis intersects the ellipsoid at the 0 deg latitude, 0 deg longitude point.
+ *
+ * @version 1.0 July 2005
+ */
+/** @internal */
+export class Ellipsoid {
+    /**
+     * Create a new ellipsoid.
+     * @param code the code.
+     * @param name the name.
+     * @param unitCode the code of the unit of measure.
+     * @param a the semi-major axis.
+     * @param invF the inverse flattening (value like 300, not like 1/300).
+     * @param b the semi-minor axis.
+     */
+    constructor(code, name, unitCode, a, invF, b) {
+        /* Store parameters */
+        this._code = code;
+        this._name = name;
+        this._unitCode = unitCode;
+        this._a = (unitCode == Unit.METER) ? a : Registry.getUnit(this._unitCode).toStandard(a);
+        this._invF = invF;
+        this._b = (unitCode == Unit.METER) ? b : Registry.getUnit(this._unitCode).toStandard(b);
+        /* Derive parameters */
+        if (this._invF == 0.0)
+            this._invF = a / (a - b);
+        this._f = (1.0 / this._invF);
+        if (this._b == 0.0)
+            this._b = this._a * (1.0 - this._f);
+        this._e = Math.sqrt(2.0 * this._f - this._f * this._f);
+        this._e2 = this._e * this._e; // or = (a*a-b*b)/(a*a);
+    }
+    /**
+     * Get the code.
+     * @return the code.
+     */
+    getCode() {
+        return this._code;
+    }
+    /**
+     * Get the name.
+     * @return the name.
+     */
+    getName() {
+        return this._name;
+    }
+    /**
+     * Get the code of the unit of measure.
+     * @return the code of the unit of measure.
+     */
+    getUnitCode() {
+        return this._unitCode;
+    }
+    /**
+     * Get the semi-major axis (in meter).
+     * @return the semi-major axis.
+     */
+    getA() {
+        return this._a;
+    }
+    /**
+     * Get the semi-minor axis (in meter).
+     * @return the semi-minor axis.
+     */
+    getB() {
+        return this._b;
+    }
+    /**
+     * Get the flattening.
+     * @return the flattening.
+     */
+    getF() {
+        return this._f;
+    }
+    /**
+     * Get the inverse flattening.
+     * @return the inverse flattening.
+     */
+    getInvF() {
+        return this._invF;
+    }
+    /**
+     * Get the eccentricity.
+     * @return the eccentricity.
+     */
+    getE() {
+        return this._e;
+    }
+    /**
+     * Get the radius of curvature of the ellipsoid in the plane of the meridian at a given latitude (how latitude curves).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getMeridianRadius(lat) {
+        // Formula: see http://en.wikipedia.org/wiki/Latitude
+        let t = this._e * Math.sin(lat);
+        return this._a * (1.0 - this._e2) / Math.pow(1.0 - t * t, 3.5);
+    }
+    /**
+     * Get the radius of the circle defined by all points at a certain latitude (a circle of latitude, also called a parallel) (a at latitude 0, 0 at latitude 90) (how longitude curves).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getParallelRadius(lat) {
+        // Formula: see http://en.wikipedia.org/wiki/Longitude
+        let t = this._e * Math.sin(lat);
+        return this._a * Math.cos(lat) / Math.sqrt(1.0 - t * t);
+    }
+    /**
+     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).
+     * @param lat the latitude (in radians).
+     * @return a radius (in meter).
+     */
+    getPrimeVerticalRadius(lat) {
+        let t = this._e * Math.sin(lat);
+        return this._a / Math.sqrt(1.0 - t * t);
+    }
+    /**
+     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).
+     * @param lat the latitude (in radians).
+     * @param sinLat the sinus of the latitude.
+     * @return a radius (in meter).
+     */
+    getPrimeVerticalRadius2(lat, sinLat) {
+        let t = this._e * sinLat;
+        return this._a / Math.sqrt(1.0 - t * t);
+    }
+    /**
+     * Get the number of meter per radian of longitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per radian.
+     */
+    getMeterPerRadOfLon(lat) {
+        return this.getParallelRadius(lat);
+    }
+    /**
+     * Get the number of meter per degree of longitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per degree.
+     */
+    getMeterPerDegreeOfLon(lat) {
+        return this.getMeterPerRadOfLon(lat) * Math.PI / 180.0;
+    }
+    /**
+     * Get the number of meter per radian of latitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per radian.
+     */
+    getMeterPerRadOfLat(lat) {
+        return this.getMeridianRadius(lat);
+    }
+    /**
+     * Get the number of meter per degree of latitude.
+     * @param lat the latitude (in radians).
+     * @return the number of meter per degree.
+     */
+    getMeterPerDegreeOfLat(lat) {
+        return this.getMeterPerRadOfLat(lat) * Math.PI / 180.0;
+    }
+    /**
+     * Calculate the secant.
+     * @param v an angle (in radians).
+     * @return the secant.
+     */
+    static sec(v) {
+        return 1.0 / Math.cos(v);
+    }
+    /**
+     * Convert 3D geographic coordinates to 3D geocentric coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geographic the geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter).
+     * @param geocentric the new geocentric coordinates (in meter) (the result object).
+     */
+    toGeoCentric(geographic, geocentric) {
+        /* Get the parameters */
+        let lon = geographic.getX();
+        let lat = geographic.getY();
+        let h = geographic.getZ();
+        /* Calculate */
+        let sinLat = Math.sin(lat);
+        let v = this.getPrimeVerticalRadius2(lat, sinLat);
+        let s = (v + h) * Math.cos(lat);
+        let x = s * Math.cos(lon);
+        let y = s * Math.sin(lon);
+        let z = ((1.0 - this._e2) * v + h) * sinLat;
+        /* Store the new coordinates */
+        geocentric.setX(x);
+        geocentric.setY(y);
+        geocentric.setZ(z);
+    }
+    /**
+     * Convert 3D geographic coordinates to 3D geocentric coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geographic the geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter).
+     * @param geocentric the new geocentric coordinates (in meter) (the result object).
+     */
+    toGeoCentricDeg(geographic, geocentric) {
+        /* Get the parameters */
+        let lon = geographic.getX() / 180.0 * Math.PI;
+        let lat = geographic.getY() / 180.0 * Math.PI;
+        let h = geographic.getZ();
+        /* Calculate */
+        let sinLat = Math.sin(lat);
+        let v = this.getPrimeVerticalRadius2(lat, sinLat);
+        let s = (v + h) * Math.cos(lat);
+        let x = s * Math.cos(lon);
+        let y = s * Math.sin(lon);
+        let z = ((1.0 - this._e2) * v + h) * sinLat;
+        /* Store the new coordinates */
+        geocentric.setX(x);
+        geocentric.setY(y);
+        geocentric.setZ(z);
+    }
+    /**
+     * Convert 3D geocentric coordinates to 3D geographic coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geocentric the geocentric coordinates (in meter).
+     * @param geographic the new geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter) (the result object).
+     */
+    toGeoGraphic(geocentric, geographic) {
+        /* Get the parameters */
+        let x = geocentric.getX();
+        let y = geocentric.getY();
+        let z = geocentric.getZ();
+        /* Calculate */
+        let r = Math.sqrt(x * x + y * y);
+        let ir = (1.0 / r);
+        let lat = Math.atan(z * ir);
+        for (let i = 0; i < 7; i++) {
+            let sinLat = Math.sin(lat);
+            let vi = this.getPrimeVerticalRadius2(lat, sinLat);
+            lat = Math.atan((z + this._e2 * vi * sinLat) * ir);
+        }
+        let lon = Math.atan2(y, x);
+        let v = this.getPrimeVerticalRadius(lat);
+        let h = x * Ellipsoid.sec(lon) * Ellipsoid.sec(lat) - v;
+        /* Store the new coordinates */
+        geographic.setX(lon);
+        geographic.setY(lat);
+        geographic.setZ(h);
+        //
+        // loop lat error:
+        //
+        // i=0: 0.0031239393816600014 // km level
+        // i=1: 6.753891540589585E-6 // m level
+        // i=2: 1.4568768968992174E-8 // mm level
+        // i=3: 3.1426083957342144E-11 // um level
+        // i=4: 6.783462680459706E-14 // nm level
+        // i=5: 1.1102230246251565E-16 // atomic level
+        // i=6: 0.0
+        //
+        // EPSG example for the WGS84 datum/ellipsoid:
+        // X = 3771793.97;
+        // Y =  140253.34;
+        // Z = 5124304.35;
+        // lat = 53.8093944; (degrees)
+        // lon =  2.12955; (degrees)
+        // h = 73.001873;
+    }
+    /**
+     * Convert 3D geocentric coordinates to 3D geographic coordinates.
+     * See the EPSG Guidance Note, 2.2.1
+     * @param geocentric the geocentric coordinates (in meter).
+     * @param geographic the new geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter) (the result object).
+     */
+    toGeoGraphicDeg(geocentric, geographic) {
+        /* Get the parameters */
+        let x = geocentric.getX();
+        let y = geocentric.getY();
+        let z = geocentric.getZ();
+        /* Calculate */
+        let r = Math.sqrt(x * x + y * y);
+        let ir = (1.0 / r);
+        let lat = Math.atan(z * ir);
+        for (let i = 0; i < 7; i++) {
+            let sinLat = Math.sin(lat);
+            let vi = this.getPrimeVerticalRadius2(lat, sinLat);
+            lat = Math.atan((z + this._e2 * vi * sinLat) * ir);
+        }
+        let lon = Math.atan2(y, x);
+        let v = this.getPrimeVerticalRadius(lat);
+        let h = x * Ellipsoid.sec(lon) * Ellipsoid.sec(lat) - v;
+        /* Store the new coordinates */
+        geographic.setX(lon / Math.PI * 180.0);
+        geographic.setY(lat / Math.PI * 180.0);
+        geographic.setZ(h);
+    }
+    /**
+     * Check if another ellipsoid is compatible with this one.
+     * @param other the other ellipsoid.
+     * @return true if compatible.
+     */
+    isCompatible(other) {
+        if (other._code == this._code)
+            return true;
+        if (other._unitCode != this._unitCode)
+            return false;
+        if (Math.abs(other._a - this._a) > 0.001)
+            return false;
+        if (Math.abs(other._b - this._b) > 0.001)
+            return false;
+        return true;
+    }
+    /**
+     * The standard toString method.
+     * @see Object#toString
+     */
+    toString() {
+        return "[Ellipsoid:code=" + this._code + ",name='" + this._name + "']";
+    }
+}
+//# sourceMappingURL=Ellipsoid.js.map

package/lib/esm/spatial/ecrs/Ellipsoid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Ellipsoid.js","sourceRoot":"","sources":["../../../../src/spatial/ecrs/Ellipsoid.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAWH,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC;AAE9B;;;;;;;;;;;;;;;GAeG;AACH,gBAAgB;AAChB,MAAM,OAAO,SAAS;IAoBlB;;;;;;;;OAQG;IACH,YAAmB,IAAW,EAAE,IAAY,EAAE,QAAe,EAAE,CAAU,EAAE,IAAa,EAAE,CAAU;QAChG,sBAAsB;QACtB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;QAClB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;QAClB,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC;QAC1B,IAAI,CAAC,EAAE,GAAG,CAAC,QAAQ,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QACxF,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;QAClB,IAAI,CAAC,EAAE,GAAG,CAAC,QAAQ,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QACxF,uBAAuB;QACvB,IAAI,IAAI,CAAC,KAAK,IAAI,GAAG;YAAE,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAChD,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QAC7B,IAAI,IAAI,CAAC,EAAE,IAAI,GAAG;YAAE,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACxD,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACvD,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,wBAAwB;IAC1D,CAAC;IAED;;;OAGG;IACI,OAAO;QACV,OAAO,IAAI,CAAC,KAAK,CAAC;IACtB,CAAC;IAED;;;OAGG;IACI,OAAO;QACV,OAAO,IAAI,CAAC,KAAK,CAAC;IACtB,CAAC;IAED;;;OAGG;IACI,WAAW;QACd,OAAO,IAAI,CAAC,SAAS,CAAC;IAC1B,CAAC;IAED;;;OAGG;IACI,IAAI;QACP,OAAO,IAAI,CAAC,EAAE,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,IAAI;QACP,OAAO,IAAI,CAAC,EAAE,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,IAAI;QACP,OAAO,IAAI,CAAC,EAAE,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,OAAO;QACV,OAAO,IAAI,CAAC,KAAK,CAAC;IACtB,CAAC;IAED;;;OAGG;IACI,IAAI;QACP,OAAO,IAAI,CAAC,EAAE,CAAC;IACnB,CAAC;IAED;;;;OAIG;IACI,iBAAiB,CAAC,GAAY;QACjC,qDAAqD;QACrD,IAAI,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,OAAO,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;IACnE,CAAC;IAED;;;;OAIG;IACI,iBAAiB,CAAC,GAAY;QACjC,sDAAsD;QACtD,IAAI,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED;;;;OAIG;IACI,sBAAsB,CAAC,GAAY;QACtC,IAAI,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED;;;;;OAKG;IACK,uBAAuB,CAAC,GAAY,EAAE,MAAe;QACzD,IAAI,CAAC,GAAY,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC;QAClC,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED;;;;OAIG;IACI,mBAAmB,CAAC,GAAY;QACnC,OAAO,IAAI,CAAC,iBAAiB,CAAC,GAAG,CAAC,CAAC;IACvC,CAAC;IAED;;;;OAIG;IACI,sBAAsB,CAAC,GAAY;QACtC,OAAO,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC;IAC3D,CAAC;IAED;;;;OAIG;IACI,mBAAmB,CAAC,GAAY;QACnC,OAAO,IAAI,CAAC,iBAAiB,CAAC,GAAG,CAAC,CAAC;IACvC,CAAC;IAED;;;;OAIG;IACI,sBAAsB,CAAC,GAAY;QACtC,OAAO,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC;IAC3D,CAAC;IAED;;;;OAIG;IACK,MAAM,CAAC,GAAG,CAAC,CAAU;QACzB,OAAO,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC7B,CAAC;IAED;;;;;OAKG;IACI,YAAY,CAAC,UAAsB,EAAE,UAAsB;QAC9D,wBAAwB;QACxB,IAAI,GAAG,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACrC,IAAI,GAAG,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACrC,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,eAAe;QACf,IAAI,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACpC,IAAI,CAAC,GAAY,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;QAC3D,IAAI,CAAC,GAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,IAAI,CAAC,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,IAAI,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QACrD,+BAA+B;QAC/B,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;;;;OAKG;IACI,eAAe,CAAC,UAAsB,EAAE,UAAsB;QACjE,wBAAwB;QACxB,IAAI,GAAG,GAAY,UAAU,CAAC,IAAI,EAAE,GAAG,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QACvD,IAAI,GAAG,GAAY,UAAU,CAAC,IAAI,EAAE,GAAG,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC;QACvD,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,eAAe;QACf,IAAI,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACpC,IAAI,CAAC,GAAY,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;QAC3D,IAAI,CAAC,GAAY,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACzC,IAAI,CAAC,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,IAAI,CAAC,GAAY,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,IAAI,CAAC,GAAY,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QACrD,+BAA+B;QAC/B,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;;;;OAKG;IACI,YAAY,CAAC,UAAsB,EAAE,UAAsB;QAC9D,wBAAwB;QACxB,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,eAAe;QACf,IAAI,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1C,IAAI,EAAE,GAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC5B,IAAI,GAAG,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,IAAI,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACpC,IAAI,EAAE,GAAY,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;YAC5D,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC;SACtD;QACD,IAAI,GAAG,GAAY,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,CAAC,GAAY,IAAI,CAAC,sBAAsB,CAAC,GAAG,CAAC,CAAC;QAClD,IAAI,CAAC,GAAY,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjE,+BAA+B;QAC/B,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,EAAE;QACF,kBAAkB;QAClB,EAAE;QACF,yCAAyC;QACzC,uCAAuC;QACvC,yCAAyC;QACzC,0CAA0C;QAC1C,yCAAyC;QACzC,8CAA8C;QAC9C,WAAW;QACX,EAAE;QACF,8CAA8C;QAC9C,kBAAkB;QAClB,kBAAkB;QAClB,kBAAkB;QAClB,8BAA8B;QAC9B,4BAA4B;QAC5B,iBAAiB;IACrB,CAAC;IAED;;;;;OAKG;IACI,eAAe,CAAC,UAAsB,EAAE,UAAsB;QACjE,wBAAwB;QACxB,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,IAAI,CAAC,GAAY,UAAU,CAAC,IAAI,EAAE,CAAC;QACnC,eAAe;QACf,IAAI,CAAC,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1C,IAAI,EAAE,GAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QAC5B,IAAI,GAAG,GAAY,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAW,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAChC,IAAI,MAAM,GAAY,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACpC,IAAI,EAAE,GAAY,IAAI,CAAC,uBAAuB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;YAC5D,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC;SACtD;QACD,IAAI,GAAG,GAAY,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpC,IAAI,CAAC,GAAY,IAAI,CAAC,sBAAsB,CAAC,GAAG,CAAC,CAAC;QAClD,IAAI,CAAC,GAAY,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACjE,+BAA+B;QAC/B,UAAU,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC;QACvC,UAAU,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC;QACvC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED;;;;OAIG;IACI,YAAY,CAAC,KAAgB;QAChC,IAAI,KAAK,CAAC,KAAK,IAAI,IAAI,CAAC,KAAK;YAAE,OAAO,IAAI,CAAC;QAC3C,IAAI,KAAK,CAAC,SAAS,IAAI,IAAI,CAAC,SAAS;YAAE,OAAO,KAAK,CAAC;QACpD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,KAAK;YAAE,OAAO,KAAK,CAAC;QACvD,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,KAAK;YAAE,OAAO,KAAK,CAAC;QACvD,OAAO,IAAI,CAAC;IAChB,CAAC;IAED;;;OAGG;IACI,QAAQ;QACX,OAAO,kBAAkB,GAAG,IAAI,CAAC,KAAK,GAAG,SAAS,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;IAC3E,CAAC;CACJ","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;\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 { Registry } from \"./Registry\";\r\nimport { Unit } from \"./Unit\";\r\n\r\n/**\r\n * Class Ellipsoid defines the parameters of an earth ellipsoid.\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 * Geocentric coordinates are defined as follows:\r\n * The point (0,0,0) denotes the center of the ellipsoid. The z-axis is defined as being parallel to the earth rotational axis, pointing towards north.\r\n * The x-axis intersects the ellipsoid at the 0 deg latitude, 0 deg longitude point.\r\n *\r\n * @version 1.0 July 2005\r\n */\r\n/** @internal */\r\nexport class Ellipsoid {\r\n    /** The code */\r\n    private _code: int32;\r\n    /** The name */\r\n    private _name: string;\r\n    /** The code of the unit of measure */\r\n    private _unitCode: int32;\r\n    /** Semi-major axis (meter) */\r\n    private _a: float64;\r\n    /** Semi-minor axis (meter) (derived) */\r\n    private _b: float64;\r\n    /** Flattening (derived) */\r\n    private _f: float64;\r\n    /** Inverse flattening */\r\n    private _invF: float64;\r\n    /** Eccentricity (derived) */\r\n    private _e: float64;\r\n    /** Eccentricity squared (derived) */\r\n    private _e2: float64;\r\n\r\n    /**\r\n     * Create a new ellipsoid.\r\n     * @param code the code.\r\n     * @param name the name.\r\n     * @param unitCode the code of the unit of measure.\r\n     * @param a the semi-major axis.\r\n     * @param invF the inverse flattening (value like 300, not like 1/300).\r\n     * @param b the semi-minor axis.\r\n     */\r\n    public constructor(code: int32, name: string, unitCode: int32, a: float64, invF: float64, b: float64) {\r\n        /* Store parameters */\r\n        this._code = code;\r\n        this._name = name;\r\n        this._unitCode = unitCode;\r\n        this._a = (unitCode == Unit.METER) ? a : Registry.getUnit(this._unitCode).toStandard(a);\r\n        this._invF = invF;\r\n        this._b = (unitCode == Unit.METER) ? b : Registry.getUnit(this._unitCode).toStandard(b);\r\n        /* Derive parameters */\r\n        if (this._invF == 0.0) this._invF = a / (a - b);\r\n        this._f = (1.0 / this._invF);\r\n        if (this._b == 0.0) this._b = this._a * (1.0 - this._f);\r\n        this._e = Math.sqrt(2.0 * this._f - this._f * this._f);\r\n        this._e2 = this._e * this._e; // or = (a*a-b*b)/(a*a);\r\n    }\r\n\r\n    /**\r\n     * Get the code.\r\n     * @return the code.\r\n     */\r\n    public getCode(): int32 {\r\n        return this._code;\r\n    }\r\n\r\n    /**\r\n     * Get the name.\r\n     * @return the name.\r\n     */\r\n    public getName(): string {\r\n        return this._name;\r\n    }\r\n\r\n    /**\r\n     * Get the code of the unit of measure.\r\n     * @return the code of the unit of measure.\r\n     */\r\n    public getUnitCode(): int32 {\r\n        return this._unitCode;\r\n    }\r\n\r\n    /**\r\n     * Get the semi-major axis (in meter).\r\n     * @return the semi-major axis.\r\n     */\r\n    public getA(): float64 {\r\n        return this._a;\r\n    }\r\n\r\n    /**\r\n     * Get the semi-minor axis (in meter).\r\n     * @return the semi-minor axis.\r\n     */\r\n    public getB(): float64 {\r\n        return this._b;\r\n    }\r\n\r\n    /**\r\n     * Get the flattening.\r\n     * @return the flattening.\r\n     */\r\n    public getF(): float64 {\r\n        return this._f;\r\n    }\r\n\r\n    /**\r\n     * Get the inverse flattening.\r\n     * @return the inverse flattening.\r\n     */\r\n    public getInvF(): float64 {\r\n        return this._invF;\r\n    }\r\n\r\n    /**\r\n     * Get the eccentricity.\r\n     * @return the eccentricity.\r\n     */\r\n    public getE(): float64 {\r\n        return this._e;\r\n    }\r\n\r\n    /**\r\n     * Get the radius of curvature of the ellipsoid in the plane of the meridian at a given latitude (how latitude curves).\r\n     * @param lat the latitude (in radians).\r\n     * @return a radius (in meter).\r\n     */\r\n    public getMeridianRadius(lat: float64): float64 {\r\n        // Formula: see http://en.wikipedia.org/wiki/Latitude\r\n        let t: float64 = this._e * Math.sin(lat);\r\n        return this._a * (1.0 - this._e2) / Math.pow(1.0 - t * t, 3.5);\r\n    }\r\n\r\n    /**\r\n     * Get the radius of the circle defined by all points at a certain latitude (a circle of latitude, also called a parallel) (a at latitude 0, 0 at latitude 90) (how longitude curves).\r\n     * @param lat the latitude (in radians).\r\n     * @return a radius (in meter).\r\n     */\r\n    public getParallelRadius(lat: float64): float64 {\r\n        // Formula: see http://en.wikipedia.org/wiki/Longitude\r\n        let t: float64 = this._e * Math.sin(lat);\r\n        return this._a * Math.cos(lat) / Math.sqrt(1.0 - t * t);\r\n    }\r\n\r\n    /**\r\n     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).\r\n     * @param lat the latitude (in radians).\r\n     * @return a radius (in meter).\r\n     */\r\n    public getPrimeVerticalRadius(lat: float64): float64 {\r\n        let t: float64 = this._e * Math.sin(lat);\r\n        return this._a / Math.sqrt(1.0 - t * t);\r\n    }\r\n\r\n    /**\r\n     * Get the radius of curvature of the ellipsoid perpendicular to the meridian at a given latitude (a at latitude 0, a*a/b at latitude 90).\r\n     * @param lat the latitude (in radians).\r\n     * @param sinLat the sinus of the latitude.\r\n     * @return a radius (in meter).\r\n     */\r\n    private getPrimeVerticalRadius2(lat: float64, sinLat: float64): float64 {\r\n        let t: float64 = this._e * sinLat;\r\n        return this._a / Math.sqrt(1.0 - t * t);\r\n    }\r\n\r\n    /**\r\n     * Get the number of meter per radian of longitude.\r\n     * @param lat the latitude (in radians).\r\n     * @return the number of meter per radian.\r\n     */\r\n    public getMeterPerRadOfLon(lat: float64): float64 {\r\n        return this.getParallelRadius(lat);\r\n    }\r\n\r\n    /**\r\n     * Get the number of meter per degree of longitude.\r\n     * @param lat the latitude (in radians).\r\n     * @return the number of meter per degree.\r\n     */\r\n    public getMeterPerDegreeOfLon(lat: float64): float64 {\r\n        return this.getMeterPerRadOfLon(lat) * Math.PI / 180.0;\r\n    }\r\n\r\n    /**\r\n     * Get the number of meter per radian of latitude.\r\n     * @param lat the latitude (in radians).\r\n     * @return the number of meter per radian.\r\n     */\r\n    public getMeterPerRadOfLat(lat: float64): float64 {\r\n        return this.getMeridianRadius(lat);\r\n    }\r\n\r\n    /**\r\n     * Get the number of meter per degree of latitude.\r\n     * @param lat the latitude (in radians).\r\n     * @return the number of meter per degree.\r\n     */\r\n    public getMeterPerDegreeOfLat(lat: float64): float64 {\r\n        return this.getMeterPerRadOfLat(lat) * Math.PI / 180.0;\r\n    }\r\n\r\n    /**\r\n     * Calculate the secant.\r\n     * @param v an angle (in radians).\r\n     * @return the secant.\r\n     */\r\n    private static sec(v: float64): float64 {\r\n        return 1.0 / Math.cos(v);\r\n    }\r\n\r\n    /**\r\n     * Convert 3D geographic coordinates to 3D geocentric coordinates.\r\n     * See the EPSG Guidance Note, 2.2.1\r\n     * @param geographic the geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter).\r\n     * @param geocentric the new geocentric coordinates (in meter) (the result object).\r\n     */\r\n    public toGeoCentric(geographic: Coordinate, geocentric: Coordinate): void {\r\n        /* Get the parameters */\r\n        let lon: float64 = geographic.getX();\r\n        let lat: float64 = geographic.getY();\r\n        let h: float64 = geographic.getZ();\r\n        /* Calculate */\r\n        let sinLat: float64 = Math.sin(lat);\r\n        let v: float64 = this.getPrimeVerticalRadius2(lat, sinLat);\r\n        let s: float64 = (v + h) * Math.cos(lat);\r\n        let x: float64 = s * Math.cos(lon);\r\n        let y: float64 = s * Math.sin(lon);\r\n        let z: float64 = ((1.0 - this._e2) * v + h) * sinLat;\r\n        /* Store the new coordinates */\r\n        geocentric.setX(x);\r\n        geocentric.setY(y);\r\n        geocentric.setZ(z);\r\n    }\r\n\r\n    /**\r\n     * Convert 3D geographic coordinates to 3D geocentric coordinates.\r\n     * See the EPSG Guidance Note, 2.2.1\r\n     * @param geographic the geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter).\r\n     * @param geocentric the new geocentric coordinates (in meter) (the result object).\r\n     */\r\n    public toGeoCentricDeg(geographic: Coordinate, geocentric: Coordinate): void {\r\n        /* Get the parameters */\r\n        let lon: float64 = geographic.getX() / 180.0 * Math.PI;\r\n        let lat: float64 = geographic.getY() / 180.0 * Math.PI;\r\n        let h: float64 = geographic.getZ();\r\n        /* Calculate */\r\n        let sinLat: float64 = Math.sin(lat);\r\n        let v: float64 = this.getPrimeVerticalRadius2(lat, sinLat);\r\n        let s: float64 = (v + h) * Math.cos(lat);\r\n        let x: float64 = s * Math.cos(lon);\r\n        let y: float64 = s * Math.sin(lon);\r\n        let z: float64 = ((1.0 - this._e2) * v + h) * sinLat;\r\n        /* Store the new coordinates */\r\n        geocentric.setX(x);\r\n        geocentric.setY(y);\r\n        geocentric.setZ(z);\r\n    }\r\n\r\n    /**\r\n     * Convert 3D geocentric coordinates to 3D geographic coordinates.\r\n     * See the EPSG Guidance Note, 2.2.1\r\n     * @param geocentric the geocentric coordinates (in meter).\r\n     * @param geographic the new geographic coordinates (lon(x) and lat(y) in radians, height(z) in meter) (the result object).\r\n     */\r\n    public toGeoGraphic(geocentric: Coordinate, geographic: Coordinate): void {\r\n        /* Get the parameters */\r\n        let x: float64 = geocentric.getX();\r\n        let y: float64 = geocentric.getY();\r\n        let z: float64 = geocentric.getZ();\r\n        /* Calculate */\r\n        let r: float64 = Math.sqrt(x * x + y * y);\r\n        let ir: float64 = (1.0 / r);\r\n        let lat: float64 = Math.atan(z * ir);\r\n        for (let i: number = 0; i < 7; i++) {\r\n            let sinLat: float64 = Math.sin(lat);\r\n            let vi: float64 = this.getPrimeVerticalRadius2(lat, sinLat);\r\n            lat = Math.atan((z + this._e2 * vi * sinLat) * ir);\r\n        }\r\n        let lon: float64 = Math.atan2(y, x);\r\n        let v: float64 = this.getPrimeVerticalRadius(lat);\r\n        let h: float64 = x * Ellipsoid.sec(lon) * Ellipsoid.sec(lat) - v;\r\n        /* Store the new coordinates */\r\n        geographic.setX(lon);\r\n        geographic.setY(lat);\r\n        geographic.setZ(h);\r\n        //\r\n        // loop lat error:\r\n        //\r\n        // i=0: 0.0031239393816600014 // km level\r\n        // i=1: 6.753891540589585E-6 // m level\r\n        // i=2: 1.4568768968992174E-8 // mm level\r\n        // i=3: 3.1426083957342144E-11 // um level\r\n        // i=4: 6.783462680459706E-14 // nm level\r\n        // i=5: 1.1102230246251565E-16 // atomic level\r\n        // i=6: 0.0\r\n        //\r\n        // EPSG example for the WGS84 datum/ellipsoid:\r\n        // X = 3771793.97;\r\n        // Y =  140253.34;\r\n        // Z = 5124304.35;\r\n        // lat = 53.8093944; (degrees)\r\n        // lon =  2.12955; (degrees)\r\n        // h = 73.001873;\r\n    }\r\n\r\n    /**\r\n     * Convert 3D geocentric coordinates to 3D geographic coordinates.\r\n     * See the EPSG Guidance Note, 2.2.1\r\n     * @param geocentric the geocentric coordinates (in meter).\r\n     * @param geographic the new geographic coordinates (lon(x) and lat(y) in degrees, height(z) in meter) (the result object).\r\n     */\r\n    public toGeoGraphicDeg(geocentric: Coordinate, geographic: Coordinate): void {\r\n        /* Get the parameters */\r\n        let x: float64 = geocentric.getX();\r\n        let y: float64 = geocentric.getY();\r\n        let z: float64 = geocentric.getZ();\r\n        /* Calculate */\r\n        let r: float64 = Math.sqrt(x * x + y * y);\r\n        let ir: float64 = (1.0 / r);\r\n        let lat: float64 = Math.atan(z * ir);\r\n        for (let i: number = 0; i < 7; i++) {\r\n            let sinLat: float64 = Math.sin(lat);\r\n            let vi: float64 = this.getPrimeVerticalRadius2(lat, sinLat);\r\n            lat = Math.atan((z + this._e2 * vi * sinLat) * ir);\r\n        }\r\n        let lon: float64 = Math.atan2(y, x);\r\n        let v: float64 = this.getPrimeVerticalRadius(lat);\r\n        let h: float64 = x * Ellipsoid.sec(lon) * Ellipsoid.sec(lat) - v;\r\n        /* Store the new coordinates */\r\n        geographic.setX(lon / Math.PI * 180.0);\r\n        geographic.setY(lat / Math.PI * 180.0);\r\n        geographic.setZ(h);\r\n    }\r\n\r\n    /**\r\n     * Check if another ellipsoid is compatible with this one.\r\n     * @param other the other ellipsoid.\r\n     * @return true if compatible.\r\n     */\r\n    public isCompatible(other: Ellipsoid): boolean {\r\n        if (other._code == this._code) return true;\r\n        if (other._unitCode != this._unitCode) return false;\r\n        if (Math.abs(other._a - this._a) > 0.001) return false;\r\n        if (Math.abs(other._b - this._b) > 0.001) return false;\r\n        return true;\r\n    }\r\n\r\n    /**\r\n     * The standard toString method.\r\n     * @see Object#toString\r\n     */\r\n    public toString(): string {\r\n        return \"[Ellipsoid:code=\" + this._code + \",name='\" + this._name + \"']\";\r\n    }\r\n}\r\n"]}

package/lib/esm/spatial/ecrs/OnlineEngine.d.ts

@@ -0,0 +1,49 @@
+/** @packageDocumentation
+ * @module OrbitGT
+ */
+import { CRSEngine } from "../crs/CRSEngine";
+import { Bounds } from "../geom/Bounds";
+import { Coordinate } from "../geom/Coordinate";
+/**
+ * Class OnlineEngine implements a CRS engine that reads CRS declarations from an online server.
+ */
+/** @internal */
+export declare class OnlineEngine extends CRSEngine {
+    /** The optional online registry */
+    private _onlineRegistry;
+    /**
+       * Create a new online engine.
+       * @return the new engine.
+       */
+    static create(): Promise<OnlineEngine>;
+    /**
+       * Create a new engine.
+       */
+    private constructor();
+    /**
+       * Set the online registry to the engine.
+       * @param onlineRegistry the registry to add.
+       */
+    private setOnlineRegistry;
+    /**
+       * CRSEngine method.
+       */
+    prepareForArea(crs: string, area: Bounds): Promise<Bounds>;
+    /**
+       * CRSEngine method.
+       */
+    transformPoint(point: Coordinate, sourceCRS: string, targetCRS: string): Coordinate;
+    /**
+       * CRSEngine method.
+       */
+    isGeocentricCRS(crs: string): boolean;
+    /**
+       * CRSEngine method.
+       */
+    isGeographicCRS(crs: string): boolean;
+    /**
+       * CRSEngine method.
+       */
+    isProjectedCRS(crs: string): boolean;
+}
+//# sourceMappingURL=OnlineEngine.d.ts.map

package/lib/esm/spatial/ecrs/OnlineEngine.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"OnlineEngine.d.ts","sourceRoot":"","sources":["../../../../src/spatial/ecrs/OnlineEngine.ts"],"names":[],"mappings":"AAIA;;GAEG;AAWH,OAAO,EAAE,SAAS,EAAE,MAAM,kBAAkB,CAAC;AAC7C,OAAO,EAAE,MAAM,EAAE,MAAM,gBAAgB,CAAC;AACxC,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAMhD;;GAEG;AACH,gBAAgB;AAChB,qBAAa,YAAa,SAAQ,SAAS;IACvC,mCAAmC;IACnC,OAAO,CAAC,eAAe,CAAiB;IAExC;;;SAGK;WACe,MAAM,IAAI,OAAO,CAAC,YAAY,CAAC;IAMnD;;SAEK;IACL,OAAO;IAKP;;;SAGK;YACS,iBAAiB;IAQ/B;;SAEK;IACiB,cAAc,CAAC,GAAG,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAWhF;;SAEK;IACW,cAAc,CAAC,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,UAAU;IAUnG;;SAEK;IACW,eAAe,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO;IAMrD;;SAEK;IACW,eAAe,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO;IAMrD;;SAEK;IACW,cAAc,CAAC,GAAG,EAAE,MAAM,GAAG,OAAO;CAKvD"}