ootk 4.0.9 → 5.1.0

package/README.md

@@ -1,10 +1,8 @@
 # ootk
 
-[![build](https://img.shields.io/github/workflow/status/thkruz/ootk/CI?style=flat-square)](https://github.com/thkruz/ootk/actions?query=workflow%3ACI)
 ![Size](https://img.shields.io/github/languages/code-size/thkruz/ootk?style=flat-square)
 [![Release](https://img.shields.io/github/v/release/thkruz/ootk?style=flat-square)](https://www.npmjs.com/package/ootk)
 [![Issues](https://img.shields.io/github/issues/thkruz/ootk?style=flat-square)](https://github.com/thkruz/ootk/issues)
-[![Coverage](https://img.shields.io/codecov/c/github/thkruz/ootk?style=flat-square)](https://codecov.io/gh/thkruz/ootk)
 [![License](https://img.shields.io/github/license/thkruz/ootk?style=flat-square)](LICENSE.MD)
 
 > An Orbital Object Toolkit in Your Web Browser
@@ -15,19 +13,6 @@ developed to simplify the math and let you focus on using the results.
 Most of the functionality was originally written for [KeepTrack](https://github.com/thkruz/keeptrack.space) and then
 later refactored into this library for others to use.
 
-## :blue_book: Table of Contents
-
-- [Installation](#Installation)
-- [Usage](#Usage)
-  - [Loading the Library](#Loading-the-Library)
-  - [Propagating a TLE](#Propagating-a-TLE)
-  - [Creating a Satellite](#Creating-a-Satellite)
-  - [Creating a Sensor](#Creating-a-Sensor)
-- [Contributing](#Contributing)
-- [Building](#Building)
-- [Contributors](#Contributors)
-- [License](#License)
-
 ## :wrench: Installation
 
 Install the library with [NPM](https://www.npmjs.com/):

package/dist/body/Celestial.d.ts

@@ -0,0 +1,96 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { AzEl, Degrees, Radians } from '../main.js';
+/**
+ * Celestial is a static class that provides methods for calculating the position of celestial objects such as the Sun,
+ * Moon, and planets in the sky. To create an instance of a Celestial object, use the Star class.
+ */
+export declare class Celestial {
+    private constructor();
+    /**
+     * Calculates the azimuth and elevation of a celestial object at a given date, latitude,
+     * longitude, right ascension, and declination.
+     * @param date - The date for which to calculate the azimuth and elevation.
+     * @param lat - The latitude of the observer.
+     * @param lon - The longitude of the observer.
+     * @param ra - The right ascension of the celestial object.
+     * @param dec - The declination of the celestial object.
+     * @returns An object containing the azimuth and elevation in degrees.
+     */
+    static azEl(date: Date, lat: Degrees, lon: Degrees, ra: Radians, dec: Radians): AzEl<Degrees>;
+    /**
+     * Atmospheric refraction in astronomy, refers to the bending of light as it passes through the Earth's
+     * atmosphere. This effect is most noticeable for celestial objects like stars and planets when they are
+     * close to the horizon. Here's a breakdown of how it works:
+     *
+     * Actual Position: Due to this bending of light, the apparent position of a celestial object is slightly
+     * different from its true position in the sky. When a star or planet is near the horizon, the effect is more
+     * pronounced because the light path passes through more of the Earth's atmosphere, which increases the amount of
+     * bending.
+     *
+     * A familiar example of atmospheric refraction is observed during sunrise and sunset. The Sun appears to
+     * be above the horizon when it is actually just below it. This is because the light from the Sun is bent
+     * upwards as it passes through the atmosphere.
+     * @param h - elevation
+     * @returns refraction
+     */
+    static atmosphericRefraction(h: Radians): Radians;
+    /**
+     * Calculate the declination. Similar to latitude on Earth, declination is another celestial coordinate.
+     * It measures how far north or south an object is from the celestial equator
+     * @param l - ecliptic longitude
+     * @param b - ecliptic latitude
+     * @returns declination
+     */
+    static declination(l: number, b: number): Radians;
+    /**
+     * Calculate the right ascension. This is a celestial coordinate used to determine the position of objects
+     * in the sky. It's analogous to longitude on Earth. Right Ascension indicates how far east an object is
+     * from the vernal equinox along the celestial equator.
+     * @param l - ecliptic longitude
+     * @param b - ecliptic latitude
+     * @returns right ascension
+     */
+    static rightAscension(l: number, b: number): Radians;
+    /**
+     * Calculate the elevation. Elevation, or altitude, is the angle between an object in the sky and the
+     * observer's local horizon. It's commonly expressed in degrees, where 0 degrees is right at the horizon
+     * and 90 degrees is directly overhead (the zenith), but we are using radians to support trigonometric
+     * functions like Math.sin() and Math.cos().
+     * @param H - siderealTime
+     * @param phi - latitude
+     * @param dec - The declination of the sun
+     * @returns elevation
+     */
+    static elevation(H: number, phi: Radians, dec: Radians): Radians;
+    /**
+     * Calculate the azimuth. This is a compass direction measurement. Azimuth measures the angle along
+     * the horizon from a specific reference direction (usually true north) to the point where a vertical
+     * line from the object intersects the horizon.
+     * @param H - siderealTime
+     * @param phi - latitude
+     * @param dec - The declination of the sun
+     * @returns azimuth in rad
+     */
+    static azimuth(H: number, phi: Radians, dec: Radians): Radians;
+}

package/dist/body/Celestial.js

@@ -0,0 +1,124 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { RAD2DEG, Sun } from '../main.js';
+/**
+ * Celestial is a static class that provides methods for calculating the position of celestial objects such as the Sun,
+ * Moon, and planets in the sky. To create an instance of a Celestial object, use the Star class.
+ */
+export class Celestial {
+    constructor() {
+        // disable constructor
+    }
+    /**
+     * Calculates the azimuth and elevation of a celestial object at a given date, latitude,
+     * longitude, right ascension, and declination.
+     * @param date - The date for which to calculate the azimuth and elevation.
+     * @param lat - The latitude of the observer.
+     * @param lon - The longitude of the observer.
+     * @param ra - The right ascension of the celestial object.
+     * @param dec - The declination of the celestial object.
+     * @returns An object containing the azimuth and elevation in degrees.
+     */
+    static azEl(date, lat, lon, ra, dec) {
+        const c = {
+            ra,
+            dec,
+            dist: 0,
+        };
+        const azEl = Sun.azEl(date, lat, lon, c);
+        const el = (azEl.el + Celestial.atmosphericRefraction(azEl.el)); // elevation correction for refraction
+        return {
+            az: (azEl.az * RAD2DEG),
+            el: (el * RAD2DEG),
+        };
+    }
+    /**
+     * Atmospheric refraction in astronomy, refers to the bending of light as it passes through the Earth's
+     * atmosphere. This effect is most noticeable for celestial objects like stars and planets when they are
+     * close to the horizon. Here's a breakdown of how it works:
+     *
+     * Actual Position: Due to this bending of light, the apparent position of a celestial object is slightly
+     * different from its true position in the sky. When a star or planet is near the horizon, the effect is more
+     * pronounced because the light path passes through more of the Earth's atmosphere, which increases the amount of
+     * bending.
+     *
+     * A familiar example of atmospheric refraction is observed during sunrise and sunset. The Sun appears to
+     * be above the horizon when it is actually just below it. This is because the light from the Sun is bent
+     * upwards as it passes through the atmosphere.
+     * @param h - elevation
+     * @returns refraction
+     */
+    static atmosphericRefraction(h) {
+        if (h < 0) {
+            h = 0;
+        }
+        return (0.0002967 / Math.tan(h + 0.00312536 / (h + 0.08901179)));
+    }
+    /**
+     * Calculate the declination. Similar to latitude on Earth, declination is another celestial coordinate.
+     * It measures how far north or south an object is from the celestial equator
+     * @param l - ecliptic longitude
+     * @param b - ecliptic latitude
+     * @returns declination
+     */
+    static declination(l, b) {
+        return Math.asin(Math.sin(b) * Math.cos(Sun.e) + Math.cos(b) * Math.sin(Sun.e) * Math.sin(l));
+    }
+    /**
+     * Calculate the right ascension. This is a celestial coordinate used to determine the position of objects
+     * in the sky. It's analogous to longitude on Earth. Right Ascension indicates how far east an object is
+     * from the vernal equinox along the celestial equator.
+     * @param l - ecliptic longitude
+     * @param b - ecliptic latitude
+     * @returns right ascension
+     */
+    static rightAscension(l, b) {
+        return Math.atan2(Math.sin(l) * Math.cos(Sun.e) - Math.tan(b) * Math.sin(Sun.e), Math.cos(l));
+    }
+    /**
+     * Calculate the elevation. Elevation, or altitude, is the angle between an object in the sky and the
+     * observer's local horizon. It's commonly expressed in degrees, where 0 degrees is right at the horizon
+     * and 90 degrees is directly overhead (the zenith), but we are using radians to support trigonometric
+     * functions like Math.sin() and Math.cos().
+     * @param H - siderealTime
+     * @param phi - latitude
+     * @param dec - The declination of the sun
+     * @returns elevation
+     */
+    static elevation(H, phi, dec) {
+        return Math.asin(Math.sin(phi) * Math.sin(dec) + Math.cos(phi) * Math.cos(dec) * Math.cos(H));
+    }
+    /**
+     * Calculate the azimuth. This is a compass direction measurement. Azimuth measures the angle along
+     * the horizon from a specific reference direction (usually true north) to the point where a vertical
+     * line from the object intersects the horizon.
+     * @param H - siderealTime
+     * @param phi - latitude
+     * @param dec - The declination of the sun
+     * @returns azimuth in rad
+     */
+    static azimuth(H, phi, dec) {
+        return (Math.PI + Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(dec) * Math.cos(phi)));
+    }
+}
+//# sourceMappingURL=Celestial.js.map

package/dist/body/Celestial.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Celestial.js","sourceRoot":"","sources":["../../src/body/Celestial.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAA6B,OAAO,EAAkB,GAAG,EAAE,MAAM,YAAY,CAAC;AAErF;;;GAGG;AACH,MAAM,OAAO,SAAS;IACpB;QACE,sBAAsB;IACxB,CAAC;IAED;;;;;;;;;OASG;IACH,MAAM,CAAC,IAAI,CAAC,IAAU,EAAE,GAAY,EAAE,GAAY,EAAE,EAAW,EAAE,GAAY;QAC3E,MAAM,CAAC,GAAU;YACf,EAAE;YACF,GAAG;YACH,IAAI,EAAE,CAAe;SACtB,CAAC;QACF,MAAM,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;QAEzC,MAAM,EAAE,GAAY,CAAC,IAAI,CAAC,EAAE,GAAG,SAAS,CAAC,qBAAqB,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,sCAAsC;QAEhH,OAAO;YACL,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAY;YAClC,EAAE,EAAE,CAAC,EAAE,GAAG,OAAO,CAAY;SAC9B,CAAC;IACJ,CAAC;IAED;;;;;;;;;;;;;;;OAeG;IACH,MAAM,CAAC,qBAAqB,CAAC,CAAU;QACrC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,CAAC,GAAY,CAAC,CAAC;QACjB,CAAC;QAED,OAAgB,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAC5E,CAAC;IAED;;;;;;OAMG;IACH,MAAM,CAAC,WAAW,CAAC,CAAS,EAAE,CAAS;QACrC,OAAgB,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACzG,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,CAAS;QACxC,OAAgB,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACzG,CAAC;IAED;;;;;;;;;OASG;IACH,MAAM,CAAC,SAAS,CAAC,CAAS,EAAE,GAAY,EAAE,GAAY;QACpD,OAAgB,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACzG,CAAC;IAED;;;;;;;;OAQG;IACH,MAAM,CAAC,OAAO,CAAC,CAAS,EAAE,GAAY,EAAE,GAAY;QAClD,OAAgB,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACnH,CAAC;CACF"}

package/dist/body/Earth.d.ts

@@ -0,0 +1,68 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, Kilometers, RadiansPerSecond, Vector3D } from '../main.js';
+import { NutationAngles } from './NutationAngles.js';
+import { PrecessionAngles } from './PrecessionAngles.js';
+export declare class Earth {
+    private constructor();
+    static readonly mu: number;
+    static readonly radiusEquator: Kilometers;
+    static readonly flattening: number;
+    static readonly radiusPolar: Kilometers;
+    static readonly radiusMean: Kilometers;
+    static readonly eccentricitySquared: number;
+    static readonly j2: number;
+    static readonly j3: number;
+    static readonly j4: number;
+    static readonly j5: number;
+    static readonly j6: number;
+    static readonly rotation: Vector3D<RadiansPerSecond>;
+    static smaToMeanMotion(semimajorAxis: Kilometers): RadiansPerSecond;
+    /**
+     * Converts revolutions per day to semi-major axis.
+     * @param rpd - The number of revolutions per day.
+     * @returns The semi-major axis value.
+     */
+    static revsPerDayToSma(rpd: number): number;
+    static precession(epoch: EpochUTC): PrecessionAngles;
+    static nutation(epoch: EpochUTC): NutationAngles;
+    static smaToDrift(semimajorAxis: number): number;
+    static smaToDriftDegrees(semimajorAxis: number): number;
+    static driftToSemimajorAxis(driftRate: number): number;
+    static driftDegreesToSma(driftRate: number): number;
+    /**
+     * Calculates the diameter of the Earth based on the satellite position.
+     * @param satPos The position of the satellite.
+     * @returns The diameter of the Earth.
+     */
+    static diameter(satPos: Vector3D): number;
+    private static readonly zetaPoly_;
+    private static readonly thetaPoly_;
+    private static readonly zedPoly_;
+    private static readonly moonAnomPoly_;
+    private static readonly sunAnomPoly_;
+    private static readonly moonLatPoly_;
+    private static readonly sunElongPoly_;
+    private static readonly moonRaanPoly_;
+    private static readonly meanEpsilonPoly_;
+}

package/dist/body/Earth.js

@@ -0,0 +1,198 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { angularDiameter, AngularDiameterMethod, asec2rad, DataHandler, DEG2RAD, earthGravityParam, evalPoly, RAD2DEG, secondsPerDay, secondsPerSiderealDay, TAU, ttasec2rad, Vector3D, } from '../main.js';
+// / Earth metrics and operations.
+export class Earth {
+    constructor() {
+        // disable constructor
+    }
+    // / Earth gravitational parameter _(km²/s³)_.
+    static mu = earthGravityParam;
+    // / Earth equatorial radius.
+    static radiusEquator = 6378.1363;
+    // / Earth coefficient of flattening _(unitless)_.
+    static flattening = 1.0 / 298.257223563;
+    // / Earth polar radius.
+    static radiusPolar = (Earth.radiusEquator * (1.0 - Earth.flattening));
+    // / Earth mean radius.
+    static radiusMean = ((2.0 * Earth.radiusEquator + Earth.radiusPolar) / 3.0);
+    // / Earth eccentricity squared _(unitless)_.
+    static eccentricitySquared = Earth.flattening * (2.0 - Earth.flattening);
+    // / Earth J2 effect coefficient _(unitless)_.
+    static j2 = 1.08262668355315e-3;
+    // / Earth J3 effect coefficient _(unitless)_.
+    static j3 = -2.53265648533224e-6;
+    // / Earth J4 effect coefficient _(unitless)_.
+    static j4 = -1.619621591367e-6;
+    // / Earth J5 effect coefficient _(unitless)_.
+    static j5 = -2.27296082868698e-7;
+    // / Earth J6 effect coefficient _(unitless)_.
+    static j6 = 5.40681239107085e-7;
+    // / Earth rotation vector _(rad/s)_.
+    static rotation = new Vector3D(0, 0, 7.292115146706979e-5);
+    // / Calculate mean motion _(rad/s)_ from a given [semimajorAxis] _(km)_.
+    static smaToMeanMotion(semimajorAxis) {
+        return Math.sqrt(Earth.mu / (semimajorAxis * semimajorAxis * semimajorAxis));
+    }
+    /**
+     * Converts revolutions per day to semi-major axis.
+     * @param rpd - The number of revolutions per day.
+     * @returns The semi-major axis value.
+     */
+    static revsPerDayToSma(rpd) {
+        return Earth.mu ** (1 / 3) / ((TAU * rpd) / secondsPerDay) ** (2 / 3);
+    }
+    // / Calculate Earth [PrecessionAngles] at a given UTC [epoch].
+    static precession(epoch) {
+        const t = epoch.toTT().toJulianCenturies();
+        const zeta = evalPoly(t, Earth.zetaPoly_);
+        const theta = evalPoly(t, Earth.thetaPoly_);
+        const zed = evalPoly(t, Earth.zedPoly_);
+        return { zeta: zeta, theta: theta, zed: zed };
+    }
+    // / Calculate Earth [NutationAngles] for a given UTC [epoch].
+    static nutation(epoch) {
+        const t = epoch.toTT().toJulianCenturies();
+        const moonAnom = evalPoly(t, Earth.moonAnomPoly_);
+        const sunAnom = evalPoly(t, Earth.sunAnomPoly_);
+        const moonLat = evalPoly(t, Earth.moonLatPoly_);
+        const sunElong = evalPoly(t, Earth.sunElongPoly_);
+        const moonRaan = evalPoly(t, Earth.moonRaanPoly_);
+        let deltaPsi = 0.0;
+        let deltaEpsilon = 0.0;
+        const dh = DataHandler.getInstance();
+        for (let i = 0; i < 4; i++) {
+            const [a1, a2, a3, a4, a5, ai, bi, ci, di] = dh.getIau1980Coeffs(i);
+            const arg = a1 * moonAnom + a2 * sunAnom + a3 * moonLat + a4 * sunElong + a5 * moonRaan;
+            const sinC = ai + bi * t;
+            const cosC = ci + di * t;
+            deltaPsi += sinC * Math.sin(arg);
+            deltaEpsilon += cosC * Math.cos(arg);
+        }
+        deltaPsi *= ttasec2rad;
+        deltaEpsilon *= ttasec2rad;
+        const meanEpsilon = evalPoly(t, Earth.meanEpsilonPoly_);
+        const epsilon = meanEpsilon + deltaEpsilon;
+        const eqEq = deltaPsi * Math.cos(meanEpsilon) +
+            0.00264 * asec2rad * Math.sin(moonRaan) +
+            0.000063 * asec2rad * Math.sin(2.0 * moonRaan);
+        const gast = epoch.gmstAngle() + eqEq;
+        return {
+            dPsi: deltaPsi,
+            dEps: deltaEpsilon,
+            mEps: meanEpsilon,
+            eps: epsilon,
+            eqEq: eqEq,
+            gast: gast,
+        };
+    }
+    // / Convert a [semimajorAxis] _(km)_ to an eastward drift rate _(rad/day)_.
+    static smaToDrift(semimajorAxis) {
+        const t = (TAU * Math.sqrt(semimajorAxis ** 3 / Earth.mu)) / secondsPerSiderealDay;
+        return (1.0 - t) * TAU;
+    }
+    // / Convert a [semimajorAxis] _(km)_ to an eastward drift rate _(°/day)_.
+    static smaToDriftDegrees(semimajorAxis) {
+        return Earth.smaToDrift(semimajorAxis) * RAD2DEG;
+    }
+    // / Convert an eastward [driftRate] _(rad/day)_ to a semimajor-axis _(km)_.
+    static driftToSemimajorAxis(driftRate) {
+        const t = (-driftRate / TAU + 1) * secondsPerSiderealDay;
+        return ((Earth.mu * t * t) / (4 * Math.PI * Math.PI)) ** (1 / 3);
+    }
+    // / Convert an eastward [driftRate] _(°/day)_ to a semimajor-axis _(km)_.
+    static driftDegreesToSma(driftRate) {
+        return Earth.driftToSemimajorAxis(DEG2RAD * driftRate);
+    }
+    /**
+     * Calculates the diameter of the Earth based on the satellite position.
+     * @param satPos The position of the satellite.
+     * @returns The diameter of the Earth.
+     */
+    static diameter(satPos) {
+        return angularDiameter(Earth.radiusEquator * 2, satPos.magnitude(), AngularDiameterMethod.Sphere);
+    }
+    // / Earth precession `zeta` polynomial coefficients.
+    static zetaPoly_ = Float64Array.from([
+        0.017998 * asec2rad,
+        0.30188 * asec2rad,
+        2306.2181 * asec2rad,
+        0.0,
+    ]);
+    // / Earth precession `theta` polynomial coefficients.
+    static thetaPoly_ = Float64Array.from([
+        -0.041833 * asec2rad,
+        -0.42665 * asec2rad,
+        2004.3109 * asec2rad,
+        0.0,
+    ]);
+    // / Earth precession `zed` polynomial coefficients.
+    static zedPoly_ = Float64Array.from([
+        0.018203 * asec2rad,
+        1.09468 * asec2rad,
+        2306.2181 * asec2rad,
+        0,
+    ]);
+    static moonAnomPoly_ = Float64Array.from([
+        1.4343e-5 * DEG2RAD,
+        0.0088553 * DEG2RAD,
+        (1325.0 * 360.0 + 198.8675605) * DEG2RAD,
+        134.96340251 * DEG2RAD,
+    ]);
+    // / Earth nutation Sun anomaly polynomial coefficients.
+    static sunAnomPoly_ = Float64Array.from([
+        3.8e-8 * DEG2RAD,
+        -0.0001537 * DEG2RAD,
+        (99.0 * 360.0 + 359.0502911) * DEG2RAD,
+        357.52910918 * DEG2RAD,
+    ]);
+    // / Earth nutation Moon latitude polynomial coefficients.
+    static moonLatPoly_ = Float64Array.from([
+        -2.88e-7 * DEG2RAD,
+        -0.003542 * DEG2RAD,
+        (1342.0 * 360.0 + 82.0174577) * DEG2RAD,
+        93.27209062 * DEG2RAD,
+    ]);
+    // / Earth nutation Sun elongation polynomial coefficients.
+    static sunElongPoly_ = Float64Array.from([
+        1.831e-6 * DEG2RAD,
+        -0.0017696 * DEG2RAD,
+        (1236.0 * 360.0 + 307.1114469) * DEG2RAD,
+        297.85019547 * DEG2RAD,
+    ]);
+    // / Earth nutation Moon right-ascension polynomial coefficients.
+    static moonRaanPoly_ = Float64Array.from([
+        2.139e-6 * DEG2RAD,
+        0.0020756 * DEG2RAD,
+        -(5.0 * 360.0 + 134.1361851) * DEG2RAD,
+        125.04455501 * DEG2RAD,
+    ]);
+    // / Earth nutation mean epsilon polynomial coefficients.
+    static meanEpsilonPoly_ = Float64Array.from([
+        0.001813 * asec2rad,
+        -0.00059 * asec2rad,
+        -46.815 * asec2rad,
+        84381.448 * asec2rad,
+    ]);
+}
+//# sourceMappingURL=Earth.js.map

package/dist/body/Earth.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Earth.js","sourceRoot":"","sources":["../../src/body/Earth.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EACL,eAAe,EACf,qBAAqB,EACrB,QAAQ,EACR,WAAW,EACX,OAAO,EACP,iBAAiB,EAEjB,QAAQ,EAER,OAAO,EAGP,aAAa,EACb,qBAAqB,EACrB,GAAG,EACH,UAAU,EACV,QAAQ,GACT,MAAM,YAAY,CAAC;AAIpB,kCAAkC;AAClC,MAAM,OAAO,KAAK;IAChB;QACE,sBAAsB;IACxB,CAAC;IAED,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,iBAAiB,CAAC;IAE/C,6BAA6B;IAC7B,MAAM,CAAU,aAAa,GAAG,SAAuB,CAAC;IAExD,kDAAkD;IAClD,MAAM,CAAU,UAAU,GAAW,GAAG,GAAG,aAAa,CAAC;IAEzD,wBAAwB;IACxB,MAAM,CAAU,WAAW,GAAG,CAAC,KAAK,CAAC,aAAa,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,UAAU,CAAC,CAAe,CAAC;IAE7F,uBAAuB;IACvB,MAAM,CAAU,UAAU,GAAG,CAAC,CAAC,GAAG,GAAG,KAAK,CAAC,aAAa,GAAG,KAAK,CAAC,WAAW,CAAC,GAAG,GAAG,CAAe,CAAC;IAEnG,6CAA6C;IAC7C,MAAM,CAAU,mBAAmB,GAAW,KAAK,CAAC,UAAU,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,UAAU,CAAC,CAAC;IAE1F,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,mBAAmB,CAAC;IAEjD,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,CAAC,mBAAmB,CAAC;IAElD,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,CAAC,iBAAiB,CAAC;IAEhD,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,CAAC,mBAAmB,CAAC;IAElD,8CAA8C;IAC9C,MAAM,CAAU,EAAE,GAAW,mBAAmB,CAAC;IAEjD,qCAAqC;IACrC,MAAM,CAAU,QAAQ,GAAG,IAAI,QAAQ,CACrC,CAAqB,EACrB,CAAqB,EACrB,oBAAwC,CACzC,CAAC;IAEF,yEAAyE;IACzE,MAAM,CAAC,eAAe,CAAC,aAAyB;QAC9C,OAAO,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,aAAa,GAAG,aAAa,GAAG,aAAa,CAAC,CAAqB,CAAC;IACnG,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,eAAe,CAAC,GAAW;QAChC,OAAO,KAAK,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxE,CAAC;IAED,+DAA+D;IAC/D,MAAM,CAAC,UAAU,CAAC,KAAe;QAC/B,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,iBAAiB,EAAE,CAAC;QAC3C,MAAM,IAAI,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,SAAS,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,UAAU,CAAC,CAAC;QAC5C,MAAM,GAAG,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,QAAQ,CAAC,CAAC;QAExC,OAAO,EAAE,IAAI,EAAE,IAAe,EAAE,KAAK,EAAE,KAAgB,EAAE,GAAG,EAAE,GAAc,EAAE,CAAC;IACjF,CAAC;IAED,8DAA8D;IAC9D,MAAM,CAAC,QAAQ,CAAC,KAAe;QAC7B,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,CAAC,iBAAiB,EAAE,CAAC;QAC3C,MAAM,QAAQ,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,aAAa,CAAC,CAAC;QAClD,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;QAChD,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,YAAY,CAAC,CAAC;QAChD,MAAM,QAAQ,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,aAAa,CAAC,CAAC;QAClD,MAAM,QAAQ,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,aAAa,CAAC,CAAC;QAClD,IAAI,QAAQ,GAAG,GAAG,CAAC;QACnB,IAAI,YAAY,GAAG,GAAG,CAAC;QACvB,MAAM,EAAE,GAAG,WAAW,CAAC,WAAW,EAAE,CAAC;QAErC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;YACpE,MAAM,GAAG,GAAG,EAAE,GAAG,QAAQ,GAAG,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,OAAO,GAAG,EAAE,GAAG,QAAQ,GAAG,EAAE,GAAG,QAAQ,CAAC;YACxF,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YACzB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAEzB,QAAQ,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACjC,YAAY,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvC,CAAC;QACD,QAAQ,IAAI,UAAU,CAAC;QACvB,YAAY,IAAI,UAAU,CAAC;QAC3B,MAAM,WAAW,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,gBAAgB,CAAC,CAAC;QACxD,MAAM,OAAO,GAAG,WAAW,GAAG,YAAY,CAAC;QAC3C,MAAM,IAAI,GACR,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC;YAChC,OAAO,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC;YACvC,QAAQ,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC;QACjD,MAAM,IAAI,GAAG,KAAK,CAAC,SAAS,EAAE,GAAG,IAAI,CAAC;QAEtC,OAAO;YACL,IAAI,EAAE,QAAmB;YACzB,IAAI,EAAE,YAAuB;YAC7B,IAAI,EAAE,WAAsB;YAC5B,GAAG,EAAE,OAAkB;YACvB,IAAI,EAAE,IAAe;YACrB,IAAI,EAAE,IAAe;SACtB,CAAC;IACJ,CAAC;IAED,4EAA4E;IAC5E,MAAM,CAAC,UAAU,CAAC,aAAqB;QACrC,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,aAAa,IAAI,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,qBAAqB,CAAC;QAEnF,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;IACzB,CAAC;IAED,0EAA0E;IAC1E,MAAM,CAAC,iBAAiB,CAAC,aAAqB;QAC5C,OAAO,KAAK,CAAC,UAAU,CAAC,aAAa,CAAC,GAAG,OAAO,CAAC;IACnD,CAAC;IAED,4EAA4E;IAC5E,MAAM,CAAC,oBAAoB,CAAC,SAAiB;QAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,SAAS,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,qBAAqB,CAAC;QAEzD,OAAO,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACnE,CAAC;IAED,0EAA0E;IAC1E,MAAM,CAAC,iBAAiB,CAAC,SAAiB;QACxC,OAAO,KAAK,CAAC,oBAAoB,CAAC,OAAO,GAAG,SAAS,CAAC,CAAC;IACzD,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,QAAQ,CAAC,MAAgB;QAC9B,OAAO,eAAe,CAAC,KAAK,CAAC,aAAa,GAAG,CAAC,EAAE,MAAM,CAAC,SAAS,EAAE,EAAE,qBAAqB,CAAC,MAAM,CAAC,CAAC;IACpG,CAAC;IAED,qDAAqD;IAC7C,MAAM,CAAU,SAAS,GAAiB,YAAY,CAAC,IAAI,CAAC;QAClE,QAAQ,GAAG,QAAQ;QACnB,OAAO,GAAG,QAAQ;QAClB,SAAS,GAAG,QAAQ;QACpB,GAAG;KACJ,CAAC,CAAC;IAEH,sDAAsD;IAC9C,MAAM,CAAU,UAAU,GAAiB,YAAY,CAAC,IAAI,CAAC;QACnE,CAAC,QAAQ,GAAG,QAAQ;QACpB,CAAC,OAAO,GAAG,QAAQ;QACnB,SAAS,GAAG,QAAQ;QACpB,GAAG;KACJ,CAAC,CAAC;IAEH,oDAAoD;IAC5C,MAAM,CAAU,QAAQ,GAAiB,YAAY,CAAC,IAAI,CAAC;QACjE,QAAQ,GAAG,QAAQ;QACnB,OAAO,GAAG,QAAQ;QAClB,SAAS,GAAG,QAAQ;QACpB,CAAC;KACF,CAAC,CAAC;IAEK,MAAM,CAAU,aAAa,GAAiB,YAAY,CAAC,IAAI,CAAC;QACtE,SAAS,GAAG,OAAO;QACnB,SAAS,GAAG,OAAO;QACnB,CAAC,MAAM,GAAG,KAAK,GAAG,WAAW,CAAC,GAAG,OAAO;QACxC,YAAY,GAAG,OAAO;KACvB,CAAC,CAAC;IAEH,wDAAwD;IAChD,MAAM,CAAU,YAAY,GAAiB,YAAY,CAAC,IAAI,CAAC;QACrE,MAAM,GAAG,OAAO;QAChB,CAAC,SAAS,GAAG,OAAO;QACpB,CAAC,IAAI,GAAG,KAAK,GAAG,WAAW,CAAC,GAAG,OAAO;QACtC,YAAY,GAAG,OAAO;KACvB,CAAC,CAAC;IAEH,0DAA0D;IAClD,MAAM,CAAU,YAAY,GAAiB,YAAY,CAAC,IAAI,CAAC;QACrE,CAAC,OAAO,GAAG,OAAO;QAClB,CAAC,QAAQ,GAAG,OAAO;QACnB,CAAC,MAAM,GAAG,KAAK,GAAG,UAAU,CAAC,GAAG,OAAO;QACvC,WAAW,GAAG,OAAO;KACtB,CAAC,CAAC;IAEH,2DAA2D;IACnD,MAAM,CAAU,aAAa,GAAiB,YAAY,CAAC,IAAI,CAAC;QACtE,QAAQ,GAAG,OAAO;QAClB,CAAC,SAAS,GAAG,OAAO;QACpB,CAAC,MAAM,GAAG,KAAK,GAAG,WAAW,CAAC,GAAG,OAAO;QACxC,YAAY,GAAG,OAAO;KACvB,CAAC,CAAC;IAEH,iEAAiE;IACzD,MAAM,CAAU,aAAa,GAAiB,YAAY,CAAC,IAAI,CAAC;QACtE,QAAQ,GAAG,OAAO;QAClB,SAAS,GAAG,OAAO;QACnB,CAAC,CAAC,GAAG,GAAG,KAAK,GAAG,WAAW,CAAC,GAAG,OAAO;QACtC,YAAY,GAAG,OAAO;KACvB,CAAC,CAAC;IAEH,yDAAyD;IACjD,MAAM,CAAU,gBAAgB,GAAiB,YAAY,CAAC,IAAI,CAAC;QACzE,QAAQ,GAAG,QAAQ;QACnB,CAAC,OAAO,GAAG,QAAQ;QACnB,CAAC,MAAM,GAAG,QAAQ;QAClB,SAAS,GAAG,QAAQ;KACrB,CAAC,CAAC"}