ootk 3.0.0 → 4.0.0-0

package/lib/propagator/Sgp4Propagator.js

@@ -0,0 +1,55 @@
+/**
+ * @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 { Propagator } from './Propagator';
+// / SGP4 propagator.
+export class Sgp4Propagator extends Propagator {
+    // / Create a new [Sgp4Propagator] object from a [TLE].
+    constructor(tle) {
+        super();
+        this.tle = tle;
+        this._checkpoints = [];
+        this._cacheState = tle.state.toJ2000();
+    }
+    get state() {
+        return this._cacheState;
+    }
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        throw new Error('Maneuvers cannot be modelled with SGP4.');
+    }
+    maneuver(maneuver, interval = 60.0) {
+        throw new Error('Maneuvers cannot be modelled with SGP4.');
+    }
+    propagate(epoch) {
+        this._cacheState = this.tle.propagate(epoch).toJ2000();
+        return this._cacheState;
+    }
+    reset() {
+        this._cacheState = this.tle.state.toJ2000();
+    }
+    checkpoint() {
+        this._checkpoints.push(this._cacheState);
+        return this._checkpoints.length - 1;
+    }
+    clearCheckpoints() {
+        this._checkpoints = [];
+    }
+    restore(index) {
+        this._cacheState = this._checkpoints[index];
+    }
+}
+//# sourceMappingURL=Sgp4Propagator.js.map

package/lib/propagator/Sgp4Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Sgp4Propagator.js","sourceRoot":"","sources":["../../src/propagator/Sgp4Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAOH,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C,qBAAqB;AACrB,MAAM,OAAO,cAAe,SAAQ,UAAU;IAC5C,uDAAuD;IACvD,YAAoB,GAAQ;QAC1B,KAAK,EAAE,CAAC;QADU,QAAG,GAAH,GAAG,CAAK;QAMpB,iBAAY,GAAY,EAAE,CAAC;QAJjC,IAAI,CAAC,WAAW,GAAG,GAAG,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IACzC,CAAC;IAKD,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,iBAAiB,CAAC,KAAe,EAAE,MAAgB,EAAE,SAAmB,EAAE,QAAQ,GAAG,IAAI;QACvF,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,IAAI;QACxC,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,SAAS,CAAC,KAAe;QACvB,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,OAAO,EAAE,CAAC;QAEvD,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,KAAK;QACH,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IAC9C,CAAC;IAED,UAAU;QACR,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;QAEzC,OAAO,IAAI,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;IACtC,CAAC;IAED,gBAAgB;QACd,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC;IACzB,CAAC;IAED,OAAO,CAAC,KAAa;QACnB,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC;IAC9C,CAAC;CACF"}

package/lib/smoothing/ExponentialSmoothing.d.ts

@@ -0,0 +1,24 @@
+/**
+ * @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 } from 'ootk-core';
+export declare class ExponentialSmoothing {
+    private constructor();
+    static smooth(xs: number[], alpha: number): number[];
+    static smoothDouble(xs: number[], alpha: number, beta: number): number[];
+    static smoothTime(epochs: EpochUTC[], xs: number[], timeConstant: number): number[];
+}

package/lib/smoothing/ExponentialSmoothing.js

@@ -0,0 +1,82 @@
+/**
+ * @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/>.
+ */
+// / Exponential data smoothing methods.
+export class ExponentialSmoothing {
+    constructor() {
+        // disable constructor
+    }
+    /*
+     * Perform exponential smoothing on data set [xs] using the provided
+     * data smoothing factor [alpha] _(0.0 <= alpha <= 1.0)_.
+     */
+    static smooth(xs, alpha) {
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+            }
+            else {
+                ss.push(alpha * xs[i] + (1 - alpha) * ss[i - 1]);
+            }
+        }
+        return ss;
+    }
+    /*
+     * Perform exponential smoothing on data set [xs] using the provided data
+     * smoothing factor [alpha] _(0.0 <= alpha <= 1.0)_ and trend smoothing
+     * factor [beta] _(0.0 <= beta <= 1.0)_.
+     */
+    static smoothDouble(xs, alpha, beta) {
+        const bs = [];
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+                bs.push(xs[1] - xs[0]);
+            }
+            else {
+                ss.push(alpha * xs[i] + (1 - alpha) * (ss[i - 1] + bs[i - 1]));
+                bs.push(beta * (ss[i] - ss[i - 1]) + (1 - beta) * bs[i - 1]);
+            }
+        }
+        return ss;
+    }
+    /*
+     * Perform exponential smoothing on time series data set [xs] correlated
+     * with the [epochs] array using the provided [timeConstant]
+     * value _(seconds)_.
+     *
+     * Note: the [timeConstant] is the amount of time for the smoothed response
+     * of a unit step function to reach ~63.2% of the original signal.
+     */
+    static smoothTime(epochs, xs, timeConstant) {
+        const ts = epochs.map((e) => e.posix);
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+            }
+            else {
+                const a = 1 - Math.exp(-(ts[i] - ts[i - 1]) / timeConstant);
+                ss.push(a * xs[i] + (1 - a) * ss[i - 1]);
+            }
+        }
+        return ss;
+    }
+}
+//# sourceMappingURL=ExponentialSmoothing.js.map

package/lib/smoothing/ExponentialSmoothing.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ExponentialSmoothing.js","sourceRoot":"","sources":["../../src/smoothing/ExponentialSmoothing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAIH,wCAAwC;AACxC,MAAM,OAAO,oBAAoB;IAC/B;QACE,sBAAsB;IACxB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,MAAM,CAAC,EAAY,EAAE,KAAa;QACvC,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aAChB;iBAAM;gBACL,EAAE,CAAC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAClD;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,YAAY,CAAC,EAAY,EAAE,KAAa,EAAE,IAAY;QAC3D,MAAM,EAAE,GAAa,EAAE,CAAC;QACxB,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACf,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aACxB;iBAAM;gBACL,EAAE,CAAC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/D,EAAE,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAC9D;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,UAAU,CAAC,MAAkB,EAAE,EAAY,EAAE,YAAoB;QACtE,MAAM,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;QACtC,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aAChB;iBAAM;gBACL,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC;gBAE5D,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAC1C;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;CACF"}