ootk 2.0.1 → 4.0.0-0

package/lib/orbit_determination/GoodingIOD.js

@@ -0,0 +1,244 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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 { Earth, J2000, TAU, Vector3D } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator';
+import { GibbsIOD } from './GibbsIOD';
+import { LambertIOD } from './LambertIOD';
+/**
+ * Gooding angles-only initial orbit determination.
+ *
+ * Used for orbit determination from three optical observations.
+ */
+export class GoodingIOD {
+    constructor(o1, o2, o3, mu = Earth.mu) {
+        this._vObserverPosition1 = Vector3D.origin;
+        this._vObserverPosition2 = Vector3D.origin;
+        this._vObserverPosition3 = Vector3D.origin;
+        this._r = 0.0;
+        this._v = 0.0;
+        this._t = 0.0;
+        this._r1 = 0.0;
+        this._r2 = 0.0;
+        this._r3 = 0.0;
+        this._rho1 = 0.0;
+        this._rho2 = 0.0;
+        this._rho3 = 0.0;
+        this._d1 = 0.0;
+        this._d3 = 0.0;
+        this._facFiniteDiff = 0.0;
+        this._forceModel = new ForceModel().setGravity(1.0);
+        this._mu = mu;
+        this._o1 = o1;
+        this._o2 = o2;
+        this._o3 = o3;
+    }
+    _getPositionOnLoS2({ e1, r01, e3, r03, t13, t12, nRev, posigrade, }) {
+        const p1 = this._vObserverPosition1.add(e1.scale(r01));
+        this._r1 = p1.magnitude();
+        const p3 = this._vObserverPosition3.add(e3.scale(r03));
+        this._r3 = p3.magnitude();
+        const p13 = p1.cross(p3);
+        let th = Math.atan2(p13.magnitude(), p1.dot(p3));
+        if (!posigrade) {
+            th = TAU - th;
+        }
+        const v1 = new Float64Array(2);
+        const exitflag = LambertIOD.solve(this._r1, this._r3, th, t13, nRev, v1);
+        if (exitflag) {
+            const pn = p1.cross(p3);
+            const pt = pn.cross(p1);
+            let rt = pt.magnitude();
+            if (!posigrade) {
+                rt = -rt;
+            }
+            const vel1 = p1.scale(v1[0] / this._r1).add(pt.scale(v1[1] / rt));
+            const p2 = new RungeKutta89Propagator(new J2000(this._o1.epoch, p1, vel1), this._forceModel).propagate(this._o1.epoch.roll(t12)).position;
+            return p2;
+        }
+        return null;
+    }
+    _modifyIterate(lineOfSight1, lineOfSight3) {
+        const r13 = this._vObserverPosition3.subtract(this._vObserverPosition1);
+        this._d1 = r13.dot(lineOfSight1);
+        this._d3 = r13.dot(lineOfSight3);
+        const d2 = lineOfSight1.dot(lineOfSight3);
+        const d4 = 1.0 - d2 * d2;
+        this._rho1 = Math.max((this._d1 - this._d3 * d2) / d4, 0.0);
+        this._rho3 = Math.max((this._d1 * d2 - this._d3) / d4, 0.0);
+    }
+    _computeDerivatives({ x, y, lineOfSight1, lineOfSight3, pin, ein, t13, t12, nrev, direction, fd, gd, }) {
+        const p = pin.normalize();
+        const en = ein.normalize();
+        const dx = this._facFiniteDiff * x;
+        const dy = this._facFiniteDiff * y;
+        const cm1 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x - dx,
+            e3: lineOfSight3,
+            r03: y,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this._vObserverPosition2);
+        const fm1 = p.dot(cm1);
+        const gm1 = en.dot(cm1);
+        const cp1 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x + dx,
+            e3: lineOfSight3,
+            r03: y,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this._vObserverPosition2);
+        const fp1 = p.dot(cp1);
+        const gp1 = en.dot(cp1);
+        const fx = (fp1 - fm1) / (2.0 * dx);
+        const gx = (gp1 - gm1) / (2.0 * dx);
+        const cm3 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x,
+            e3: lineOfSight3,
+            r03: y - dy,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this._vObserverPosition2);
+        const fm3 = p.dot(cm3);
+        const gm3 = en.dot(cm3);
+        const cp3 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x,
+            e3: lineOfSight3,
+            r03: y + dy,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this._vObserverPosition2);
+        const fp3 = p.dot(cp3);
+        const gp3 = en.dot(cp3);
+        const fy = (fp3 - fm3) / (2.0 * dy);
+        const gy = (gp3 - gm3) / (2.0 * dy);
+        fd[0] = fx;
+        fd[1] = fy;
+        gd[0] = gx;
+        gd[1] = gy;
+    }
+    solve(r1Init, r3Init, nRev = 0, direction = true) {
+        const lineOfSight1 = this._o1.observation.lineOfSight();
+        const lineOfSight2 = this._o2.observation.lineOfSight();
+        const lineOfSight3 = this._o3.observation.lineOfSight();
+        this._r = Math.max(r1Init, r3Init);
+        this._v = Math.sqrt(this._mu / this._r);
+        this._t = this._r / this._v;
+        this._vObserverPosition1 = this._o1.site.position.scale(1.0 / this._r);
+        this._vObserverPosition2 = this._o2.site.position.scale(1.0 / this._r);
+        this._vObserverPosition3 = this._o3.site.position.scale(1.0 / this._r);
+        const maxiter = 100;
+        this._solveRangeProblem({
+            rho1init: r1Init / this._r,
+            rho3init: r3Init / this._r,
+            t13: this._o3.epoch.difference(this._o1.epoch) / this._t,
+            t12: this._o2.epoch.difference(this._o1.epoch) / this._t,
+            nrev: nRev,
+            direction,
+            lineOfSight1,
+            lineOfSight2,
+            lineOfSight3,
+            maxIterations: maxiter,
+        });
+        const gibbs = new GibbsIOD(this._mu);
+        const p1 = this._vObserverPosition1.add(lineOfSight1.scale(this._rho1)).scale(this._r);
+        const p2 = this._vObserverPosition2.add(lineOfSight2.scale(this._rho2)).scale(this._r);
+        const p3 = this._vObserverPosition3.add(lineOfSight3.scale(this._rho3)).scale(this._r);
+        return gibbs.solve(p1, p2, p3, this._o2.epoch, this._o3.epoch);
+    }
+    _solveRangeProblem({ rho1init, rho3init, t13, t12, nrev, direction, lineOfSight1, lineOfSight2, lineOfSight3, maxIterations, }) {
+        const arbf = 1e-6;
+        const cvtol = 1e-14;
+        this._rho1 = rho1init;
+        this._rho3 = rho3init;
+        let iter = 0;
+        let stoppingCriterion = 10.0 * cvtol;
+        while (iter < maxIterations && Math.abs(stoppingCriterion) > cvtol) {
+            this._facFiniteDiff = arbf;
+            const p2 = this._getPositionOnLoS2({
+                e1: lineOfSight1,
+                r01: this._rho1,
+                e3: lineOfSight3,
+                r03: this._rho3,
+                t13,
+                t12,
+                nRev: nrev,
+                posigrade: direction,
+            });
+            if (p2 === null) {
+                this._modifyIterate(lineOfSight1, lineOfSight3);
+            }
+            else {
+                this._r2 = p2.magnitude();
+                const c = p2.subtract(this._vObserverPosition2);
+                this._rho2 = c.magnitude();
+                const cr = lineOfSight2.dot(c);
+                const u = lineOfSight2.cross(c);
+                const p = u.cross(lineOfSight2).normalize();
+                const ent = lineOfSight2.cross(p);
+                const enr = ent.magnitude();
+                if (enr === 0.0) {
+                    return;
+                }
+                const en = ent.normalize();
+                const fc = p.dot(c);
+                const fd = new Float64Array(2);
+                const gd = new Float64Array(2);
+                this._computeDerivatives({
+                    x: this._rho1,
+                    y: this._rho3,
+                    lineOfSight1,
+                    lineOfSight3,
+                    pin: p,
+                    ein: en,
+                    t13,
+                    t12,
+                    nrev,
+                    direction,
+                    fd,
+                    gd,
+                });
+                const fr1 = fd[0];
+                const fr3 = fd[1];
+                const gr1 = gd[0];
+                const gr3 = gd[1];
+                const detj = fr1 * gr3 - fr3 * gr1;
+                this._d3 = (-gr3 * fc) / detj;
+                this._d1 = (gr1 * fc) / detj;
+                this._rho1 += this._d3;
+                this._rho3 += this._d1;
+                const den = Math.max(cr, this._r2);
+                stoppingCriterion = fc / den;
+            }
+            ++iter;
+        }
+    }
+}
+//# sourceMappingURL=GoodingIOD.js.map

package/lib/orbit_determination/GoodingIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"GoodingIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/GoodingIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AACxD,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AACjD,OAAO,EAAE,sBAAsB,EAAE,MAAM,sCAAsC,CAAC;AAE9E,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAyC1C;;;;GAIG;AACH,MAAM,OAAO,UAAU;IAsBrB,YAAY,EAAsB,EAAE,EAAsB,EAAE,EAAsB,EAAE,KAAa,KAAK,CAAC,EAAE;QAjBjG,wBAAmB,GAAa,QAAQ,CAAC,MAAM,CAAC;QAChD,wBAAmB,GAAa,QAAQ,CAAC,MAAM,CAAC;QAChD,wBAAmB,GAAa,QAAQ,CAAC,MAAM,CAAC;QAChD,OAAE,GAAG,GAAG,CAAC;QACT,OAAE,GAAG,GAAG,CAAC;QACT,OAAE,GAAG,GAAG,CAAC;QACT,QAAG,GAAG,GAAG,CAAC;QACV,QAAG,GAAG,GAAG,CAAC;QACV,QAAG,GAAG,GAAG,CAAC;QACV,UAAK,GAAG,GAAG,CAAC;QACZ,UAAK,GAAG,GAAG,CAAC;QACZ,UAAK,GAAG,GAAG,CAAC;QACZ,QAAG,GAAG,GAAG,CAAC;QACV,QAAG,GAAG,GAAG,CAAC;QACV,mBAAc,GAAG,GAAG,CAAC;QACZ,gBAAW,GAAe,IAAI,UAAU,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;QAG1E,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;IAChB,CAAC;IAEO,kBAAkB,CAAC,EACzB,EAAE,EACF,GAAG,EACH,EAAE,EACF,GAAG,EACH,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,GACe;QACxB,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvD,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAE1B,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvD,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAE1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QAEzB,IAAI,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,SAAS,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,IAAI,CAAC,SAAS,EAAE;YACd,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;SACf;QAED,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,QAAQ,GAAG,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC;QAEzE,IAAI,QAAQ,EAAE;YACZ,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YAExB,IAAI,EAAE,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;YAExB,IAAI,CAAC,SAAS,EAAE;gBACd,EAAE,GAAG,CAAC,EAAE,CAAC;aACV;YAED,MAAM,IAAI,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YAElE,MAAM,EAAE,GAAG,IAAI,sBAAsB,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC,SAAS,CACpG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CACzB,CAAC,QAAQ,CAAC;YAEX,OAAO,EAAE,CAAC;SACX;QAED,OAAO,IAAI,CAAC;IACd,CAAC;IAEO,cAAc,CAAC,YAAsB,EAAE,YAAsB;QACnE,MAAM,GAAG,GAAG,IAAI,CAAC,mBAAmB,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAExE,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACjC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QAEzB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC;QAC5D,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC;IAC9D,CAAC;IAEO,mBAAmB,CAAC,EAC1B,CAAC,EACD,CAAC,EACD,YAAY,EACZ,YAAY,EACZ,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,EACT,EAAE,EACF,EAAE,GACuB;QACzB,MAAM,CAAC,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAE3B,MAAM,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QACnC,MAAM,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,CAAC,CAAC;QAEnC,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAE;YACX,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAE;YACX,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QAEpC,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAE;YACX,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAE;YACX,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QAEpC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACb,CAAC;IAED,KAAK,CAAC,MAAc,EAAE,MAAc,EAAE,IAAI,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI;QAC9D,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAAE,CAAC;QACxD,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAAE,CAAC;QACxD,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAAE,CAAC;QAExD,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QACnC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACxC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAE5B,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACvE,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACvE,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAEvE,MAAM,OAAO,GAAG,GAAG,CAAC;QAEpB,IAAI,CAAC,kBAAkB,CAAC;YACtB,QAAQ,EAAE,MAAM,GAAG,IAAI,CAAC,EAAE;YAC1B,QAAQ,EAAE,MAAM,GAAG,IAAI,CAAC,EAAE;YAC1B,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE;YACxD,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAE;YACxD,IAAI,EAAE,IAAI;YACV,SAAS;YACT,YAAY;YACZ,YAAY;YACZ,YAAY;YACZ,aAAa,EAAE,OAAO;SACvB,CAAC,CAAC;QACH,MAAM,KAAK,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAEvF,OAAO,KAAK,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IACjE,CAAC;IAEO,kBAAkB,CAAC,EACzB,QAAQ,EACR,QAAQ,EACR,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,EACT,YAAY,EACZ,YAAY,EACZ,YAAY,EACZ,aAAa,GACW;QACxB,MAAM,IAAI,GAAG,IAAI,CAAC;QAClB,MAAM,KAAK,GAAG,KAAK,CAAC;QAEpB,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC;QACtB,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC;QAEtB,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,iBAAiB,GAAG,IAAI,GAAG,KAAK,CAAC;QAErC,OAAO,IAAI,GAAG,aAAa,IAAI,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC,GAAG,KAAK,EAAE;YAClE,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC;YAE3B,MAAM,EAAE,GAAG,IAAI,CAAC,kBAAkB,CAAC;gBACjC,EAAE,EAAE,YAAY;gBAChB,GAAG,EAAE,IAAI,CAAC,KAAK;gBACf,EAAE,EAAE,YAAY;gBAChB,GAAG,EAAE,IAAI,CAAC,KAAK;gBACf,GAAG;gBACH,GAAG;gBACH,IAAI,EAAE,IAAI;gBACV,SAAS,EAAE,SAAS;aACrB,CAAC,CAAC;YAEH,IAAI,EAAE,KAAK,IAAI,EAAE;gBACf,IAAI,CAAC,cAAc,CAAC,YAAY,EAAE,YAAY,CAAC,CAAC;aACjD;iBAAM;gBACL,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;gBAC1B,MAAM,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;gBAEhD,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAE/B,MAAM,CAAC,GAAG,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAChC,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,SAAS,EAAE,CAAC;gBAC5C,MAAM,GAAG,GAAG,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAElC,MAAM,GAAG,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;gBAE5B,IAAI,GAAG,KAAK,GAAG,EAAE;oBACf,OAAO;iBACR;gBAED,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;gBAE3B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAEpB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;gBAC/B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;gBAE/B,IAAI,CAAC,mBAAmB,CAAC;oBACvB,CAAC,EAAE,IAAI,CAAC,KAAK;oBACb,CAAC,EAAE,IAAI,CAAC,KAAK;oBACb,YAAY;oBACZ,YAAY;oBACZ,GAAG,EAAE,CAAC;oBACN,GAAG,EAAE,EAAE;oBACP,GAAG;oBACH,GAAG;oBACH,IAAI;oBACJ,SAAS;oBACT,EAAE;oBACF,EAAE;iBACH,CAAC,CAAC;gBAEH,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;gBAEnC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;gBAC9B,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;gBAE7B,IAAI,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC;gBACvB,IAAI,CAAC,KAAK,IAAI,IAAI,CAAC,GAAG,CAAC;gBAEvB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;gBAEnC,iBAAiB,GAAG,EAAE,GAAG,GAAG,CAAC;aAC9B;YAED,EAAE,IAAI,CAAC;SACR;IACH,CAAC;CACF"}

package/{dist/objects/space-object.d.ts → lib/orbit_determination/HerrickGibbsIOD.d.ts}

@@ -1,13 +1,8 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
- *
- * @file The SpaceObject class is meant to provide properties to satellites/stars that
- * do not apply to ground based objects.
+ * @author @thkruz Theodore Kruczek
  *
  * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 Theodore Kruczek
+ * @Copyright (c) 2020-2024 Theodore Kruczek
  *
  * 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
@@ -20,19 +15,19 @@
  * 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 { BaseObject } from './base-object';
-import { SpaceObjectType } from '../types/types';
-interface ObjectInfo {
-    name?: string;
-    type?: SpaceObjectType;
-    rcs?: number;
-    vmag?: number;
-}
-export declare class SpaceObject extends BaseObject {
-    rcs: number;
-    vmag: number;
-    constructor(info: ObjectInfo);
-    getRcs(): number;
-    getVmag(): number;
+import { EpochUTC, J2000, Vector3D } from 'ootk-core';
+/**
+ * Herrik-Gibbs 3-position initial orbit determination.
+ *
+ * Possibly better than regular Gibbs IOD for closely spaced position
+ * vectors (less than 5°).
+ */
+export declare class HerrickGibbsIOD {
+    mu: number;
+    /**
+     * Create a new [HerrickGibbsIOD] object with optional
+     * gravitational parameter [mu].
+     */
+    constructor(mu?: number);
+    solve(r1: Vector3D, t1: EpochUTC, r2: Vector3D, t2: EpochUTC, r3: Vector3D, t3: EpochUTC): J2000;
 }
-export {};

package/lib/orbit_determination/HerrickGibbsIOD.js

@@ -0,0 +1,49 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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 { Earth, J2000 } from 'ootk-core';
+/**
+ * Herrik-Gibbs 3-position initial orbit determination.
+ *
+ * Possibly better than regular Gibbs IOD for closely spaced position
+ * vectors (less than 5°).
+ */
+export class HerrickGibbsIOD {
+    /**
+     * Create a new [HerrickGibbsIOD] object with optional
+     * gravitational parameter [mu].
+     */
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+        // Nothing to do here.
+    }
+    // / Attempt to create a state estimate from three inertial position vectors.
+    solve(r1, t1, r2, t2, r3, t3) {
+        const dt31 = t3.difference(t1);
+        const dt32 = t3.difference(t2);
+        const dt21 = t2.difference(t1);
+        const r1m = r1.magnitude();
+        const r2m = r2.magnitude();
+        const r3m = r3.magnitude();
+        const vA = r1.scale(-dt32 * (1.0 / (dt21 * dt31) + this.mu / (12.0 * r1m * r1m * r1m)));
+        const vB = r2.scale((dt32 - dt21) * (1.0 / (dt21 * dt32) + this.mu / (12.0 * r2m * r2m * r2m)));
+        const vC = r3.scale(dt21 * (1.0 / (dt32 * dt31) + this.mu / (12.0 * r3m * r3m * r3m)));
+        const v2 = vA.add(vB).add(vC);
+        return new J2000(t2, r2, v2);
+    }
+}
+//# sourceMappingURL=HerrickGibbsIOD.js.map

package/lib/orbit_determination/HerrickGibbsIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"HerrickGibbsIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/HerrickGibbsIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAY,KAAK,EAAY,MAAM,WAAW,CAAC;AAE7D;;;;;GAKG;AACH,MAAM,OAAO,eAAe;IAC1B;;;OAGG;IACH,YAAmB,KAAa,KAAK,CAAC,EAAE;QAArB,OAAE,GAAF,EAAE,CAAmB;QACtC,sBAAsB;IACxB,CAAC;IAED,6EAA6E;IAC7E,KAAK,CAAC,EAAY,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY;QACtF,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACxF,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAChG,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAE9B,OAAO,IAAI,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAC/B,CAAC;CACF"}

package/lib/orbit_determination/LambertIOD.d.ts

@@ -0,0 +1,46 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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, J2000, Vector3D } from 'ootk-core';
+export declare class LambertIOD {
+    mu: number;
+    /**
+     * Create a new [LambertIOD] object with optional
+     * gravitational parameter [mu].
+     */
+    constructor(mu?: number);
+    /**
+     * Try to guess the short path argument given an [interceptor] and
+     * [target] state.
+     */
+    static useShortPath(interceptor: J2000, target: J2000): boolean;
+    private static timeOfFlight_;
+    /**
+     * Attempt to solve output velocity [v1] _(km/s)_ given radii [r1] and
+     * [r2] _(canonical)_, sweep angle [dth] _(rad)_, time of flight [tau]
+     * _(canonical)_, and number of revolutions _(mRev)_.
+     */
+    static solve(r1: number, r2: number, dth: number, tau: number, mRev: number, v1: Float64Array): boolean;
+    /**
+     * Estimate a state vector for inertial position [p1] _(km)_ given the
+     * two epoch and positions.
+     */
+    estimate(p1: Vector3D, p2: Vector3D, t1: EpochUTC, t2: EpochUTC, { posigrade, nRev }?: {
+        posigrade?: boolean;
+        nRev?: number;
+    }): J2000 | null;
+}

package/lib/orbit_determination/LambertIOD.js

@@ -0,0 +1,199 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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 { Earth, J2000 } from 'ootk-core';
+// / Lambert two-position and time initial orbit determination.
+export class LambertIOD {
+    /**
+     * Create a new [LambertIOD] object with optional
+     * gravitational parameter [mu].
+     */
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+        // Nothing to do here.
+    }
+    /**
+     * Try to guess the short path argument given an [interceptor] and
+     * [target] state.
+     */
+    static useShortPath(interceptor, target) {
+        const transN = interceptor.position.cross(target.position);
+        const h = interceptor.position.cross(interceptor.velocity);
+        return h.dot(transN) >= 0;
+    }
+    static timeOfFlight_(x, longway, mrev, minSma, speri, chord) {
+        const a = minSma / (1.0 - x * x);
+        let tof;
+        if (Math.abs(x) < 1) {
+            const beta = longway * 2.0 * Math.asin(Math.sqrt((speri - chord) / (2.0 * a)));
+            const alpha = 2.0 * Math.acos(x);
+            tof = a * Math.sqrt(a) * (alpha - Math.sin(alpha) - (beta - Math.sin(beta)) + 2.0 * Math.PI * mrev);
+        }
+        else {
+            const alpha = 2.0 * Math.acosh(x);
+            const beta = longway * 2.0 * Math.asinh(Math.sqrt((speri - chord) / (-2.0 * a)));
+            tof = -a * Math.sqrt(-a) * (Math.sinh(alpha) - alpha - (Math.sinh(beta) - beta));
+        }
+        return tof;
+    }
+    /**
+     * Attempt to solve output velocity [v1] _(km/s)_ given radii [r1] and
+     * [r2] _(canonical)_, sweep angle [dth] _(rad)_, time of flight [tau]
+     * _(canonical)_, and number of revolutions _(mRev)_.
+     */
+    static solve(r1, r2, dth, tau, mRev, v1) {
+        const leftBranch = dth < Math.PI;
+        let longway = 1;
+        if (dth > Math.PI) {
+            longway = -1;
+        }
+        const m = Math.abs(mRev);
+        const rtof = Math.abs(tau);
+        const theta = dth;
+        const chord = Math.sqrt(r1 * r1 + r2 * r2 - 2.0 * r1 * r2 * Math.cos(theta));
+        const speri = 0.5 * (r1 + r2 + chord);
+        const minSma = 0.5 * speri;
+        const lambda = longway * Math.sqrt(1.0 - chord / speri);
+        const logt = Math.log(rtof);
+        let in1;
+        let in2;
+        let x1;
+        let x2;
+        if (m === 0) {
+            in1 = -0.6523333;
+            in2 = 0.6523333;
+            x1 = Math.log(1.0 + in1);
+            x2 = Math.log(1.0 + in2);
+        }
+        else {
+            if (!leftBranch) {
+                in1 = -0.523334;
+                in2 = -0.223334;
+            }
+            else {
+                in1 = 0.723334;
+                in2 = 0.523334;
+            }
+            x1 = Math.tan((in1 * Math.PI) / 2);
+            x2 = Math.tan((in2 * Math.PI) / 2);
+        }
+        const tof1 = LambertIOD.timeOfFlight_(in1, longway, m, minSma, speri, chord);
+        const tof2 = LambertIOD.timeOfFlight_(in2, longway, m, minSma, speri, chord);
+        let y1;
+        let y2;
+        if (m === 0) {
+            y1 = Math.log(tof1) - logt;
+            y2 = Math.log(tof2) - logt;
+        }
+        else {
+            y1 = tof1 - rtof;
+            y2 = tof2 - rtof;
+        }
+        let err = 1e20;
+        let iterations = 0;
+        const tol = 1e-13;
+        const maxiter = 50;
+        let xnew = 0.0;
+        while (err > tol && iterations < maxiter) {
+            xnew = (x1 * y2 - y1 * x2) / (y2 - y1);
+            let xt;
+            if (m === 0) {
+                xt = Math.exp(xnew) - 1.0;
+            }
+            else {
+                xt = (Math.atan(xnew) * 2.0) / Math.PI;
+            }
+            const tof = LambertIOD.timeOfFlight_(xt, longway, m, minSma, speri, chord);
+            let ynew;
+            if (m === 0) {
+                ynew = Math.log(tof) - logt;
+            }
+            else {
+                ynew = tof - rtof;
+            }
+            x1 = x2;
+            x2 = xnew;
+            y1 = y2;
+            y2 = ynew;
+            err = Math.abs(x1 - xnew);
+            ++iterations;
+        }
+        if (err > tol) {
+            return false;
+        }
+        let x;
+        if (m === 0) {
+            x = Math.exp(xnew) - 1.0;
+        }
+        else {
+            x = (Math.atan(xnew) * 2.0) / Math.PI;
+        }
+        const sma = minSma / (1.0 - x * x);
+        let eta;
+        if (x < 1) {
+            const alfa = 2.0 * Math.acos(x);
+            const beta = longway * 2.0 * Math.asin(Math.sqrt((speri - chord) / (2.0 * sma)));
+            const psi = (alfa - beta) / 2.0;
+            const sinPsi = Math.sin(psi);
+            const etaSq = (2.0 * sma * sinPsi * sinPsi) / speri;
+            eta = Math.sqrt(etaSq);
+        }
+        else {
+            const gamma = 2.0 * Math.acosh(x);
+            const delta = longway * 2.0 * Math.asinh(Math.sqrt((chord - speri) / (2.0 * sma)));
+            const psi = (gamma - delta) / 2.0;
+            const sinhPsi = Math.sinh(psi);
+            const etaSq = (-2.0 * sma * sinhPsi * sinhPsi) / speri;
+            eta = Math.sqrt(etaSq);
+        }
+        const vr1 = (1.0 / eta) * Math.sqrt(1.0 / minSma) * ((2.0 * lambda * minSma) / r1 - (lambda + x * eta));
+        const vt1 = (1.0 / eta) * Math.sqrt(1.0 / minSma) * Math.sqrt(r2 / r1) * Math.sin(dth / 2.0);
+        v1[0] = vr1;
+        v1[1] = vt1;
+        return true;
+    }
+    /**
+     * Estimate a state vector for inertial position [p1] _(km)_ given the
+     * two epoch and positions.
+     */
+    estimate(p1, p2, t1, t2, { posigrade = true, nRev = 0 } = {}) {
+        const r1 = p1.magnitude();
+        const r2 = p2.magnitude();
+        const tof = t2.difference(t1);
+        const r = Math.max(r1, r2);
+        const v = Math.sqrt(this.mu / r);
+        const t = r / v;
+        let dth = p1.angle(p2);
+        if (!posigrade) {
+            dth = 2 * Math.PI - dth;
+        }
+        const vDep = new Float64Array(2);
+        const exitFlag = LambertIOD.solve(r1 / r, r2 / r, dth, tof / t, nRev, vDep);
+        if (exitFlag) {
+            const pn = p1.cross(p2);
+            const pt = pn.cross(p1);
+            let rt = pt.magnitude();
+            if (!posigrade) {
+                rt = -rt;
+            }
+            const vel1 = p1.scale((v * vDep[0]) / r1).add(pt.scale((v * vDep[1]) / rt));
+            return new J2000(t1, p1, vel1);
+        }
+        return null;
+    }
+}
+//# sourceMappingURL=LambertIOD.js.map

package/lib/orbit_determination/LambertIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"LambertIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/LambertIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAY,KAAK,EAAY,MAAM,WAAW,CAAC;AAE7D,+DAA+D;AAC/D,MAAM,OAAO,UAAU;IACrB;;;OAGG;IACH,YAAmB,KAAa,KAAK,CAAC,EAAE;QAArB,OAAE,GAAF,EAAE,CAAmB;QACtC,sBAAsB;IACxB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,YAAY,CAAC,WAAkB,EAAE,MAAa;QACnD,MAAM,MAAM,GAAG,WAAW,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAG,WAAW,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,CAAC,QAAQ,CAAC,CAAC;QAE3D,OAAO,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAC5B,CAAC;IAEO,MAAM,CAAC,aAAa,CAC1B,CAAS,EACT,OAAe,EACf,IAAY,EACZ,MAAc,EACd,KAAa,EACb,KAAa;QAEb,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACjC,IAAI,GAAW,CAAC;QAEhB,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE;YACnB,MAAM,IAAI,GAAG,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/E,MAAM,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAEjC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;SACrG;aAAM;YACL,MAAM,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAClC,MAAM,IAAI,GAAG,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YAEjF,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;SAClF;QAED,OAAO,GAAG,CAAC;IACb,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,KAAK,CAAC,EAAU,EAAE,EAAU,EAAE,GAAW,EAAE,GAAW,EAAE,IAAY,EAAE,EAAgB;QAC3F,MAAM,UAAU,GAAG,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC;QACjC,IAAI,OAAO,GAAG,CAAC,CAAC;QAEhB,IAAI,GAAG,GAAG,IAAI,CAAC,EAAE,EAAE;YACjB,OAAO,GAAG,CAAC,CAAC,CAAC;SACd;QAED,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACzB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC3B,MAAM,KAAK,GAAG,GAAG,CAAC;QAElB,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7E,MAAM,KAAK,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;QAEtC,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QAC3B,MAAM,MAAM,GAAG,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC;QACxD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAE5B,IAAI,GAAW,CAAC;QAChB,IAAI,GAAW,CAAC;QAChB,IAAI,EAAU,CAAC;QACf,IAAI,EAAU,CAAC;QAEf,IAAI,CAAC,KAAK,CAAC,EAAE;YACX,GAAG,GAAG,CAAC,SAAS,CAAC;YACjB,GAAG,GAAG,SAAS,CAAC;YAChB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YACzB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;SAC1B;aAAM;YACL,IAAI,CAAC,UAAU,EAAE;gBACf,GAAG,GAAG,CAAC,QAAQ,CAAC;gBAChB,GAAG,GAAG,CAAC,QAAQ,CAAC;aACjB;iBAAM;gBACL,GAAG,GAAG,QAAQ,CAAC;gBACf,GAAG,GAAG,QAAQ,CAAC;aAChB;YACD,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACnC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;SACpC;QACD,MAAM,IAAI,GAAG,UAAU,CAAC,aAAa,CAAC,GAAG,EAAE,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAC7E,MAAM,IAAI,GAAG,UAAU,CAAC,aAAa,CAAC,GAAG,EAAE,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;QAE7E,IAAI,EAAU,CAAC;QACf,IAAI,EAAU,CAAC;QAEf,IAAI,CAAC,KAAK,CAAC,EAAE;YACX,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;YAC3B,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;SAC5B;aAAM;YACL,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;YACjB,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;SAClB;QACD,IAAI,GAAG,GAAG,IAAI,CAAC;QACf,IAAI,UAAU,GAAG,CAAC,CAAC;QACnB,MAAM,GAAG,GAAG,KAAK,CAAC;QAClB,MAAM,OAAO,GAAG,EAAE,CAAC;QACnB,IAAI,IAAI,GAAG,GAAG,CAAC;QAEf,OAAO,GAAG,GAAG,GAAG,IAAI,UAAU,GAAG,OAAO,EAAE;YACxC,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YACvC,IAAI,EAAU,CAAC;YAEf,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;aAC3B;iBAAM;gBACL,EAAE,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;aACxC;YAED,MAAM,GAAG,GAAG,UAAU,CAAC,aAAa,CAAC,EAAE,EAAE,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;YAE3E,IAAI,IAAY,CAAC;YAEjB,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;aAC7B;iBAAM;gBACL,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;aACnB;YAED,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,IAAI,CAAC;YACV,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,IAAI,CAAC;YAEV,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;YAC1B,EAAE,UAAU,CAAC;SACd;QAED,IAAI,GAAG,GAAG,GAAG,EAAE;YACb,OAAO,KAAK,CAAC;SACd;QAED,IAAI,CAAS,CAAC;QAEd,IAAI,CAAC,KAAK,CAAC,EAAE;YACX,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC;SAC1B;aAAM;YACL,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;SACvC;QAED,MAAM,GAAG,GAAG,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAEnC,IAAI,GAAW,CAAC;QAEhB,IAAI,CAAC,GAAG,CAAC,EAAE;YACT,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAChC,MAAM,IAAI,GAAG,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YACjF,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,GAAG,CAAC;YAChC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAC7B,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,KAAK,CAAC;YAEpD,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;SACxB;aAAM;YACL,MAAM,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAClC,MAAM,KAAK,GAAG,OAAO,GAAG,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YACnF,MAAM,GAAG,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,GAAG,CAAC;YAClC,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YAC/B,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,OAAO,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC;YAEvD,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;SACxB;QAED,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACxG,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAE7F,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACZ,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAEZ,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;OAGG;IACH,QAAQ,CACN,EAAY,EACZ,EAAY,EACZ,EAAY,EACZ,EAAY,EACZ,EAAE,SAAS,GAAG,IAAI,EAAE,IAAI,GAAG,CAAC,KAA6C,EAAE;QAE3E,MAAM,EAAE,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC1B,MAAM,EAAE,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC1B,MAAM,GAAG,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAE9B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAEhB,IAAI,GAAG,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QAEvB,IAAI,CAAC,SAAS,EAAE;YACd,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;SACzB;QAED,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,QAAQ,GAAG,UAAU,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAE5E,IAAI,QAAQ,EAAE;YACZ,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACxB,IAAI,EAAE,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;YAExB,IAAI,CAAC,SAAS,EAAE;gBACd,EAAE,GAAG,CAAC,EAAE,CAAC;aACV;YACD,MAAM,IAAI,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YAE5E,OAAO,IAAI,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;SAChC;QAED,OAAO,IAAI,CAAC;IACd,CAAC;CACF"}

package/lib/orbit_determination/ModifiedGoodingIOD.d.ts

@@ -0,0 +1,41 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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 { J2000 } from 'ootk-core';
+import { ObservationOptical } from './../observation/ObservationOptical';
+type SolveOptions = {
+    nRev?: number;
+    direction?: boolean;
+    posSearch?: number;
+    velSearch?: number;
+    tolerance?: number;
+    printIter?: boolean;
+};
+/**
+ * Gooding angles-only initial orbit determination.
+ *
+ * Used for orbit determination from multiple optical observations.
+ */
+export declare class ModifiedGoodingIOD {
+    private _observations;
+    private _mu;
+    constructor(observations: ObservationOptical[], mu?: number);
+    private _createInitial;
+    private _createErrorFunction;
+    solve(r0: number, rN: number, { nRev, direction, posSearch, velSearch, tolerance, printIter, }: SolveOptions): J2000;
+}
+export {};