@itwin/core-orbitgt 4.0.0-dev.46 → 4.0.0-dev.48

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

@@ -1,337 +1,337 @@
-"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.Ellipsoid = void 0;
-const Registry_1 = require("./Registry");
-const Unit_1 = require("./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 */
-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_1.Unit.METER) ? a : Registry_1.Registry.getUnit(this._unitCode).toStandard(a);
-        this._invF = invF;
-        this._b = (unitCode == Unit_1.Unit.METER) ? b : Registry_1.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 + "']";
-    }
-}
-exports.Ellipsoid = Ellipsoid;
+"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.Ellipsoid = void 0;
+const Registry_1 = require("./Registry");
+const Unit_1 = require("./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 */
+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_1.Unit.METER) ? a : Registry_1.Registry.getUnit(this._unitCode).toStandard(a);
+        this._invF = invF;
+        this._b = (unitCode == Unit_1.Unit.METER) ? b : Registry_1.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 + "']";
+    }
+}
+exports.Ellipsoid = Ellipsoid;
 //# sourceMappingURL=Ellipsoid.js.map