ootk 3.0.0 → 4.0.0-1

package/dist/observation/Observation.d.ts

@@ -0,0 +1,89 @@
+/**
+ * @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, Matrix, Vector, Vector3D } from 'ootk-core';
+import { Propagator } from '../propagator/Propagator';
+import { RandomGaussianSource } from './../operations/RandomGaussianSource';
+import { PropagatorPairs } from './PropagatorPairs';
+/**
+ * Observation data.
+ */
+export declare abstract class Observation {
+    /** Observation epoch. */
+    abstract get epoch(): EpochUTC;
+    /** Inertial observer location. */
+    abstract get site(): J2000;
+    /** Observation noise matrix. */
+    abstract get noise(): Matrix;
+    /**
+     * Return range-normalized cross line-of-sight residual for the observation
+     * when compared against a nominal state propagator.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    clos(propagator: Propagator): number;
+    /**
+     * Return relative state residual for the observation when compared against
+     * a nominal state propagator.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    ricDiff(propagator: Propagator): Vector3D;
+    /**
+     * Convert this observation to vector form.
+     * @throws Not implemented.
+     */
+    toVector(): Vector;
+    /**
+     * Compute the state derivative matrix for this observation.
+     * @param propPairs Propagator pairs to compare against.
+     * @throws Not implemented.
+     */
+    jacobian(propPairs: PropagatorPairs): Matrix;
+    /**
+     * Compute the state residual matrix for this observation.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    residual(propagator: Propagator): Matrix;
+    /**
+     * Convert this observation's noise matrix into a covariance matrix.
+     * @returns A matrix representing the noise covariance.
+     */
+    noiseCovariance(): Matrix;
+    /**
+     * Generates a noise sample from the noise covariance matrix.
+     * @param sigma - The scaling factor for the noise covariance matrix.
+     * @returns A matrix representing the noise sample.
+     */
+    noiseSample_(sigma: number): Matrix;
+    /**
+     * Randomly sample this observation in vector form within the
+     * observation noise.
+     * @param random Random number generator.
+     * @param sigma Sigma value to scale the noise by.
+     * @returns Sampled observation.
+     */
+    sampleVector(random: RandomGaussianSource, sigma: number): Vector;
+    /**
+     * Randomly sample this observation within the observation noise, scaled to
+     * the provided sigma value.
+     * @param random Random number generator.
+     * @param sigma Sigma value to scale the noise by.
+     * @throws Not implemented.
+     */
+    sample(random: RandomGaussianSource, sigma?: number): Observation;
+}

package/dist/observation/Observation.js

@@ -0,0 +1,102 @@
+/**
+ * @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/>.
+ */
+/* eslint-disable @typescript-eslint/no-unused-vars */
+/**
+ * Observation data.
+ */
+export class Observation {
+    /**
+     * Return range-normalized cross line-of-sight residual for the observation
+     * when compared against a nominal state propagator.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    clos(propagator) {
+        throw new Error('Not implemented');
+    }
+    /**
+     * Return relative state residual for the observation when compared against
+     * a nominal state propagator.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    ricDiff(propagator) {
+        throw new Error('Not implemented');
+    }
+    /**
+     * Convert this observation to vector form.
+     * @throws Not implemented.
+     */
+    toVector() {
+        throw new Error('Not implemented');
+    }
+    /**
+     * Compute the state derivative matrix for this observation.
+     * @param propPairs Propagator pairs to compare against.
+     * @throws Not implemented.
+     */
+    jacobian(propPairs) {
+        throw new Error('Not implemented');
+    }
+    /**
+     * Compute the state residual matrix for this observation.
+     * @param propagator Propagator to compare against.
+     * @throws Not implemented.
+     */
+    residual(propagator) {
+        throw new Error('Not implemented');
+    }
+    /**
+     * Convert this observation's noise matrix into a covariance matrix.
+     * @returns A matrix representing the noise covariance.
+     */
+    noiseCovariance() {
+        return this.noise.reciprocal();
+    }
+    /**
+     * Generates a noise sample from the noise covariance matrix.
+     * @param sigma - The scaling factor for the noise covariance matrix.
+     * @returns A matrix representing the noise sample.
+     */
+    noiseSample_(sigma) {
+        return this.noiseCovariance().scale(sigma).cholesky();
+    }
+    /**
+     * Randomly sample this observation in vector form within the
+     * observation noise.
+     * @param random Random number generator.
+     * @param sigma Sigma value to scale the noise by.
+     * @returns Sampled observation.
+     */
+    sampleVector(random, sigma) {
+        const chol = this.noiseSample_(sigma);
+        const gauss = random.gaussVector(this.noise.columns);
+        const meas = this.toVector();
+        return meas.add(chol.multiplyVector(gauss));
+    }
+    /**
+     * Randomly sample this observation within the observation noise, scaled to
+     * the provided sigma value.
+     * @param random Random number generator.
+     * @param sigma Sigma value to scale the noise by.
+     * @throws Not implemented.
+     */
+    sample(random, sigma = 1.0) {
+        throw new Error('Not implemented');
+    }
+}
+//# sourceMappingURL=Observation.js.map

package/dist/observation/Observation.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Observation.js","sourceRoot":"","sources":["../../src/observation/Observation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAOH,sDAAsD;AACtD;;GAEG;AACH,MAAM,OAAgB,WAAW;IAU/B;;;;;OAKG;IACH,IAAI,CAAC,UAAsB;QACzB,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;IAED;;;;;OAKG;IACH,OAAO,CAAC,UAAsB;QAC5B,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;IAED;;;OAGG;IACH,QAAQ;QACN,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;IAED;;;;OAIG;IACH,QAAQ,CAAC,SAA0B;QACjC,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;IAED;;;;OAIG;IACH,QAAQ,CAAC,UAAsB;QAC7B,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;IAED;;;OAGG;IACH,eAAe;QACb,OAAO,IAAI,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC;IACjC,CAAC;IAED;;;;OAIG;IACH,YAAY,CAAC,KAAa;QACxB,OAAO,IAAI,CAAC,eAAe,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,QAAQ,EAAE,CAAC;IACxD,CAAC;IAED;;;;;;OAMG;IACH,YAAY,CAAC,MAA4B,EAAE,KAAa;QACtD,MAAM,IAAI,GAAG,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC;QACtC,MAAM,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QACrD,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,EAAE,CAAC;QAE7B,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,cAAc,CAAC,KAAK,CAAC,CAAC,CAAC;IAC9C,CAAC;IAED;;;;;;OAMG;IACH,MAAM,CAAC,MAA4B,EAAE,KAAK,GAAG,GAAG;QAC9C,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACrC,CAAC;CACF"}

package/dist/observation/ObservationOptical.d.ts

@@ -0,0 +1,48 @@
+/**
+ * @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, Matrix, RadecTopocentric, Radians, Vector, Vector3D } from 'ootk-core';
+import { RandomGaussianSource } from '../operations/RandomGaussianSource';
+import { Propagator } from '../propagator/Propagator';
+import { Observation } from './Observation';
+import { PropagatorPairs } from './PropagatorPairs';
+export declare class ObservationOptical extends Observation {
+    private _site;
+    observation: RadecTopocentric;
+    private _noise;
+    constructor(_site: J2000, observation: RadecTopocentric, _noise?: Matrix);
+    get site(): J2000;
+    get noise(): Matrix;
+    /**
+     * Default noise matrix _(right-ascension, declination)_.
+     * Based on the Maui Optical Site noise model.
+     */
+    static defaultNoise: Matrix;
+    get epoch(): EpochUTC;
+    toVector(): Vector;
+    clos(propagator: Propagator): number;
+    ricDiff(propagator: Propagator): Vector3D;
+    sample(random: RandomGaussianSource, sigma?: number): Observation;
+    jacobian(propPairs: PropagatorPairs): Matrix;
+    residual(propagator: Propagator): Matrix;
+    /**
+     * Create a noise matrix from right ascension and declination standard deviantions.
+     * @param raSigma Right ascension standard deviation
+     * @param decSigma Declination standard deviation
+     * @returns Noise matrix
+     */
+    static noiseFromSigmas(raSigma: Radians, decSigma: Radians): Matrix;
+}

package/dist/observation/ObservationOptical.js

@@ -0,0 +1,105 @@
+/**
+ * @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 { array2d, DEG2RAD, J2000, Matrix, RadecTopocentric, RIC, Vector, Vector3D, } from 'ootk-core';
+import { Observation } from './Observation';
+import { normalizeAngle, observationDerivative, observationNoiseFromSigmas } from './ObservationUtils';
+// / Optical observation data.
+export class ObservationOptical extends Observation {
+    _site;
+    observation;
+    _noise;
+    // / Create a new [ObservationOptical] object.
+    constructor(_site, observation, _noise = ObservationOptical.defaultNoise) {
+        super();
+        this._site = _site;
+        this.observation = observation;
+        this._noise = _noise;
+    }
+    // / Inertial site location.
+    get site() {
+        return this._site;
+    }
+    // / Noise matrix.
+    get noise() {
+        return this._noise;
+    }
+    /**
+     * Default noise matrix _(right-ascension, declination)_.
+     * Based on the Maui Optical Site noise model.
+     */
+    static defaultNoise = ObservationOptical.noiseFromSigmas(0.0037 * DEG2RAD, 0.003 * DEG2RAD);
+    get epoch() {
+        return this.observation.epoch;
+    }
+    toVector() {
+        return Vector.fromList([this.observation.rightAscension, this.observation.declination]);
+    }
+    clos(propagator) {
+        const position = propagator.propagate(this.epoch).position;
+        const offset = position.subtract(this.site.position);
+        const actual = this.observation.lineOfSight().normalize();
+        const expected = offset.normalize();
+        const slantRange = offset.magnitude();
+        const theta = actual.angle(expected);
+        if (isNaN(theta)) {
+            return 0.0;
+        }
+        return 2.0 * slantRange * Math.sin(theta * 0.5);
+    }
+    ricDiff(propagator) {
+        const r0 = this.site;
+        const r1 = propagator.propagate(this.epoch);
+        const r2 = this.observation.position(this.site, r1.position.distance(r0.position));
+        return RIC.fromJ2000(new J2000(this.epoch, r2, Vector3D.origin), r1).position;
+    }
+    sample(random, sigma = 1.0) {
+        const resultEl = this.sampleVector(random, sigma).elements;
+        return new ObservationOptical(this.site, new RadecTopocentric(this.epoch, resultEl[0], resultEl[1]), this.noise);
+    }
+    jacobian(propPairs) {
+        const result = array2d(2, 6, 0.0);
+        for (let i = 0; i < 6; i++) {
+            const step = propPairs.step(i);
+            const [high, low] = propPairs.get(i);
+            const sl = low.propagate(this.epoch);
+            const sh = high.propagate(this.epoch);
+            const ol = RadecTopocentric.fromStateVector(sl, this.site);
+            const oh = RadecTopocentric.fromStateVector(sh, this.site);
+            result[0][i] = observationDerivative(oh.rightAscension, ol.rightAscension, step, true);
+            result[1][i] = observationDerivative(oh.declination, ol.declination, step, true);
+        }
+        return new Matrix(result);
+    }
+    residual(propagator) {
+        const result = array2d(2, 1, 0.0);
+        const state = propagator.propagate(this.epoch);
+        const radec = RadecTopocentric.fromStateVector(state, this.site);
+        result[0][0] = normalizeAngle(this.observation.rightAscension, radec.rightAscension);
+        result[1][0] = normalizeAngle(this.observation.declination, radec.declination);
+        return new Matrix(result);
+    }
+    /**
+     * Create a noise matrix from right ascension and declination standard deviantions.
+     * @param raSigma Right ascension standard deviation
+     * @param decSigma Declination standard deviation
+     * @returns Noise matrix
+     */
+    static noiseFromSigmas(raSigma, decSigma) {
+        return observationNoiseFromSigmas([raSigma, decSigma]);
+    }
+}
+//# sourceMappingURL=ObservationOptical.js.map

package/dist/observation/ObservationOptical.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObservationOptical.js","sourceRoot":"","sources":["../../src/observation/ObservationOptical.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EACL,OAAO,EACP,OAAO,EAEP,KAAK,EAEL,MAAM,EACN,gBAAgB,EAEhB,GAAG,EACH,MAAM,EACN,QAAQ,GACT,MAAM,WAAW,CAAC;AAGnB,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,cAAc,EAAE,qBAAqB,EAAE,0BAA0B,EAAE,MAAM,oBAAoB,CAAC;AAGvG,8BAA8B;AAC9B,MAAM,OAAO,kBAAmB,SAAQ,WAAW;IAGvC;IACD;IACC;IAJV,8CAA8C;IAC9C,YACU,KAAY,EACb,WAA6B,EAC5B,SAAiB,kBAAkB,CAAC,YAAY;QAExD,KAAK,EAAE,CAAC;QAJA,UAAK,GAAL,KAAK,CAAO;QACb,gBAAW,GAAX,WAAW,CAAkB;QAC5B,WAAM,GAAN,MAAM,CAA0C;IAG1D,CAAC;IAED,4BAA4B;IAC5B,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;IAED,kBAAkB;IAClB,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,YAAY,GAAW,kBAAkB,CAAC,eAAe,CAC9D,MAAM,GAAG,OAAkB,EAC3B,KAAK,GAAG,OAAkB,CAC3B,CAAC;IAEF,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC;IAChC,CAAC;IAED,QAAQ;QACN,OAAO,MAAM,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,WAAW,CAAC,cAAc,EAAE,IAAI,CAAC,WAAW,CAAC,WAAW,CAAC,CAAC,CAAC;IAC1F,CAAC;IAED,IAAI,CAAC,UAAsB;QACzB,MAAM,QAAQ,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC;QAC3D,MAAM,MAAM,GAAG,QAAQ,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACrD,MAAM,MAAM,GAAG,IAAI,CAAC,WAAW,CAAC,WAAW,EAAE,CAAC,SAAS,EAAE,CAAC;QAC1D,MAAM,QAAQ,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC;QACpC,MAAM,UAAU,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC;QACtC,MAAM,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QAErC,IAAI,KAAK,CAAC,KAAK,CAAC,EAAE;YAChB,OAAO,GAAG,CAAC;SACZ;QAED,OAAO,GAAG,GAAG,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;IAClD,CAAC;IAED,OAAO,CAAC,UAAsB;QAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC;QACrB,MAAM,EAAE,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC,QAAQ,CAAC,CAAC,CAAC;QAEnF,OAAO,GAAG,CAAC,SAAS,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,EAAE,EAAE,QAAQ,CAAC,MAAuC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC;IACjH,CAAC;IAED,MAAM,CAAC,MAA4B,EAAE,KAAK,GAAG,GAAG;QAC9C,MAAM,QAAQ,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,QAAqB,CAAC;QAExE,OAAO,IAAI,kBAAkB,CAC3B,IAAI,CAAC,IAAI,EACT,IAAI,gBAAgB,CAAC,IAAI,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,EAC1D,IAAI,CAAC,KAAK,CACX,CAAC;IACJ,CAAC;IAED,QAAQ,CAAC,SAA0B;QACjC,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC/B,MAAM,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,gBAAgB,CAAC,eAAe,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;YAC3D,MAAM,EAAE,GAAG,gBAAgB,CAAC,eAAe,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;YAE3D,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,CAAC,cAAc,EAAE,EAAE,CAAC,cAAc,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;YACvF,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,CAAC,WAAW,EAAE,EAAE,CAAC,WAAW,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;SAClF;QAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;IAED,QAAQ,CAAC,UAAsB;QAC7B,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAClC,MAAM,KAAK,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,KAAK,GAAG,gBAAgB,CAAC,eAAe,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QAEjE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,WAAW,CAAC,cAAc,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;QACrF,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,WAAW,CAAC,WAAW,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC;QAE/E,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;IAED;;;;;OAKG;IACH,MAAM,CAAC,eAAe,CAAC,OAAgB,EAAE,QAAiB;QACxD,OAAO,0BAA0B,CAAC,CAAC,OAAO,EAAE,QAAQ,CAAC,CAAC,CAAC;IACzD,CAAC"}

package/dist/observation/ObservationRadar.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, Kilometers, Matrix, Radians, RAE, Vector, Vector3D } from 'ootk-core';
+import { RandomGaussianSource } from '../operations/RandomGaussianSource';
+import { Propagator } from '../propagator/Propagator';
+import { Observation } from './Observation';
+import { PropagatorPairs } from './PropagatorPairs';
+export declare class ObservationRadar extends Observation {
+    private _site;
+    observation: RAE;
+    private _noise;
+    constructor(_site: J2000, observation: RAE, _noise?: Matrix);
+    private static defaultNoise;
+    get epoch(): EpochUTC;
+    get site(): J2000;
+    get noise(): Matrix;
+    toVector(): Vector;
+    clos(propagator: Propagator): number;
+    ricDiff(propagator: Propagator): Vector3D;
+    sample(random: RandomGaussianSource, sigma?: number): Observation;
+    jacobian(propPairs: PropagatorPairs): Matrix;
+    residual(propagator: Propagator): Matrix;
+    /**
+     * Create a noise matrix from the range, azimuth, and elevation standard
+     * deviations _(kilometers/radians)_.
+     * @param rngSigma - The range standard deviation _(kilometers)_.
+     * @param azSigma - The azimuth standard deviation _(radians)_.
+     * @param elSigma - The elevation standard deviation _(radians)_.
+     * @returns The noise matrix.
+     */
+    static noiseFromSigmas(rngSigma: Kilometers, azSigma: Radians, elSigma: Radians): Matrix;
+}

package/dist/observation/ObservationRadar.js

@@ -0,0 +1,96 @@
+/**
+ * @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 { array2d, DEG2RAD, J2000, Matrix, RAE, RIC, Vector, Vector3D, } from 'ootk-core';
+import { Observation } from './Observation';
+import { normalizeAngle, observationDerivative, observationNoiseFromSigmas } from './ObservationUtils';
+// / Radar observation data.
+export class ObservationRadar extends Observation {
+    _site;
+    observation;
+    _noise;
+    // / Create a new [ObservationRadar] object.
+    constructor(_site, observation, _noise = ObservationRadar.defaultNoise) {
+        super();
+        this._site = _site;
+        this.observation = observation;
+        this._noise = _noise;
+    }
+    // / Default noise matrix.
+    static defaultNoise = ObservationRadar.noiseFromSigmas(0.32, 0.015 * DEG2RAD, 0.015 * DEG2RAD);
+    get epoch() {
+        return this.observation.epoch;
+    }
+    get site() {
+        return this._site;
+    }
+    get noise() {
+        return this._noise;
+    }
+    toVector() {
+        return Vector.fromList([this.observation.rng, this.observation.azRad, this.observation.elRad]);
+    }
+    clos(propagator) {
+        const ri = propagator.propagate(this.epoch).position.subtract(this.site.position);
+        return Math.abs(this.observation.rng - ri.magnitude());
+    }
+    ricDiff(propagator) {
+        const r0 = propagator.propagate(this.epoch);
+        const r1 = RAE.fromStateVector(r0, this.site);
+        const r2 = this.observation.position(this.site, r1.azRad, r1.elRad);
+        return RIC.fromJ2000(new J2000(this.epoch, r2, Vector3D.origin), r0).position;
+    }
+    sample(random, sigma = 1.0) {
+        const result = this.sampleVector(random, sigma).elements;
+        return new ObservationRadar(this.site, new RAE(this.observation.epoch, result[0], result[1], result[2]), this.noise);
+    }
+    jacobian(propPairs) {
+        const result = array2d(3, 6, 0.0);
+        for (let i = 0; i < 6; i++) {
+            const step = propPairs.step(i);
+            const [high, low] = propPairs.get(i);
+            const sl = low.propagate(this.epoch);
+            const sh = high.propagate(this.epoch);
+            const ol = RAE.fromStateVector(sl, this.site);
+            const oh = RAE.fromStateVector(sh, this.site);
+            result[0][i] = observationDerivative(oh.rng, ol.rng, step);
+            result[1][i] = observationDerivative(oh.azRad, ol.azRad, step, true);
+            result[2][i] = observationDerivative(oh.elRad, ol.elRad, step, true);
+        }
+        return new Matrix(result);
+    }
+    residual(propagator) {
+        const result = array2d(3, 1, 0.0);
+        const state = propagator.propagate(this.epoch);
+        const razel = RAE.fromStateVector(state, this.site);
+        result[0][0] = this.observation.rng - razel.rng;
+        result[1][0] = normalizeAngle(this.observation.azRad, razel.azRad);
+        result[2][0] = normalizeAngle(this.observation.elRad, razel.elRad);
+        return new Matrix(result);
+    }
+    /**
+     * Create a noise matrix from the range, azimuth, and elevation standard
+     * deviations _(kilometers/radians)_.
+     * @param rngSigma - The range standard deviation _(kilometers)_.
+     * @param azSigma - The azimuth standard deviation _(radians)_.
+     * @param elSigma - The elevation standard deviation _(radians)_.
+     * @returns The noise matrix.
+     */
+    static noiseFromSigmas(rngSigma, azSigma, elSigma) {
+        return observationNoiseFromSigmas([rngSigma, azSigma, elSigma]);
+    }
+}
+//# sourceMappingURL=ObservationRadar.js.map

package/dist/observation/ObservationRadar.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObservationRadar.js","sourceRoot":"","sources":["../../src/observation/ObservationRadar.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EACL,OAAO,EACP,OAAO,EAEP,KAAK,EAGL,MAAM,EAEN,GAAG,EACH,GAAG,EACH,MAAM,EACN,QAAQ,GACT,MAAM,WAAW,CAAC;AAGnB,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,cAAc,EAAE,qBAAqB,EAAE,0BAA0B,EAAE,MAAM,oBAAoB,CAAC;AAGvG,4BAA4B;AAC5B,MAAM,OAAO,gBAAiB,SAAQ,WAAW;IAE3B;IAAqB;IAA0B;IADnE,4CAA4C;IAC5C,YAAoB,KAAY,EAAS,WAAgB,EAAU,SAAiB,gBAAgB,CAAC,YAAY;QAC/G,KAAK,EAAE,CAAC;QADU,UAAK,GAAL,KAAK,CAAO;QAAS,gBAAW,GAAX,WAAW,CAAK;QAAU,WAAM,GAAN,MAAM,CAAwC;IAEjH,CAAC;IAED,0BAA0B;IAClB,MAAM,CAAC,YAAY,GACzB,gBAAgB,CAAC,eAAe,CAAC,IAAkB,EAAE,KAAK,GAAG,OAAkB,EAAE,KAAK,GAAG,OAAkB,CAAC,CAAC;IAE/G,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC;IAChC,CAAC;IAED,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,QAAQ;QACN,OAAO,MAAM,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,WAAW,CAAC,GAAG,EAAE,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;IACjG,CAAC;IAED,IAAI,CAAC,UAAsB;QACzB,MAAM,EAAE,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAElF,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC,CAAC;IACzD,CAAC;IAED,OAAO,CAAC,UAAsB;QAC5B,MAAM,EAAE,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC5C,MAAM,EAAE,GAAG,GAAG,CAAC,eAAe,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,WAAW,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC;QAEpE,OAAO,GAAG,CAAC,SAAS,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,EAAE,EAAE,QAAQ,CAAC,MAAuC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC;IACjH,CAAC;IAED,MAAM,CAAC,MAA4B,EAAE,KAAK,GAAG,GAAG;QAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,QAAQ,CAAC;QAEzD,OAAO,IAAI,gBAAgB,CACzB,IAAI,CAAC,IAAI,EACT,IAAI,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,CAAe,EAAE,MAAM,CAAC,CAAC,CAAY,EAAE,MAAM,CAAC,CAAC,CAAY,CAAC,EACpG,IAAI,CAAC,KAAK,CACX,CAAC;IACJ,CAAC;IAED,QAAQ,CAAC,SAA0B;QACjC,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC/B,MAAM,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,GAAG,CAAC,eAAe,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;YAC9C,MAAM,EAAE,GAAG,GAAG,CAAC,eAAe,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;YAE9C,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;YAC3D,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;YACrE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,qBAAqB,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;SACtE;QAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;IAED,QAAQ,CAAC,UAAsB;QAC7B,MAAM,MAAM,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAClC,MAAM,KAAK,GAAG,UAAU,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,KAAK,GAAG,GAAG,CAAC,eAAe,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QAEpD,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC;QAChD,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC;QACnE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,CAAC,CAAC;QAEnE,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,eAAe,CAAC,QAAoB,EAAE,OAAgB,EAAE,OAAgB;QAC7E,OAAO,0BAA0B,CAAC,CAAC,QAAQ,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;IAClE,CAAC"}

package/dist/observation/ObservationUtils.d.ts

@@ -0,0 +1,20 @@
+/**
+ * @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 { Matrix } from 'ootk-core';
+export declare const normalizeAngle: (a: number, b: number) => number;
+export declare const observationDerivative: (xh: number, xl: number, step: number, isAngle?: boolean) => number;
+export declare const observationNoiseFromSigmas: (sigmas: number[]) => Matrix;

package/dist/observation/ObservationUtils.js

@@ -0,0 +1,32 @@
+/**
+ * @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 { Matrix } from 'ootk-core';
+export const normalizeAngle = (a, b) => {
+    const x = a - b;
+    return Math.atan2(Math.sin(x), Math.cos(x));
+};
+export const observationDerivative = (xh, xl, step, isAngle = false) => (isAngle ? normalizeAngle(xh, xl) : xh - xl) / step;
+export const observationNoiseFromSigmas = (sigmas) => {
+    const n = sigmas.length;
+    const result = Array.from({ length: n }, () => Array(n).fill(0.0));
+    for (let i = 0; i < n; i++) {
+        const s = sigmas[i];
+        result[i][i] = 1 / (s * s);
+    }
+    return new Matrix(result);
+};
+//# sourceMappingURL=ObservationUtils.js.map

package/dist/observation/ObservationUtils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ObservationUtils.js","sourceRoot":"","sources":["../../src/observation/ObservationUtils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAEnC,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEhB,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9C,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,qBAAqB,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,IAAY,EAAE,OAAO,GAAG,KAAK,EAAU,EAAE,CACrG,CAAC,OAAO,CAAC,CAAC,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC;AAEtD,MAAM,CAAC,MAAM,0BAA0B,GAAG,CAAC,MAAgB,EAAU,EAAE;IACrE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAEnE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC1B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAEpB,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;KAC5B;IAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;AAC5B,CAAC,CAAC"}

package/dist/observation/PropagatorPairs.d.ts

@@ -0,0 +1,32 @@
+/**
+ * @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/Propagator';
+export declare class PropagatorPairs {
+    private _posStep;
+    private _velStep;
+    constructor(_posStep: number, _velStep: number);
+    private _high;
+    private _low;
+    set(index: number, high: Propagator, low: Propagator): void;
+    get(index: number): [Propagator, Propagator];
+    /**
+     * Get the step size at the provided index.
+     * @param index The index.
+     * @returns The step size.
+     */
+    step(index: number): number;
+}

package/dist/observation/PropagatorPairs.js

@@ -0,0 +1,43 @@
+/**
+ * @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/>.
+ */
+export class PropagatorPairs {
+    _posStep;
+    _velStep;
+    constructor(_posStep, _velStep) {
+        this._posStep = _posStep;
+        this._velStep = _velStep;
+        // Do nothing.
+    }
+    _high = Array(6).fill(null);
+    _low = Array(6).fill(null);
+    set(index, high, low) {
+        this._high[index] = high;
+        this._low[index] = low;
+    }
+    get(index) {
+        return [this._high[index], this._low[index]];
+    }
+    /**
+     * Get the step size at the provided index.
+     * @param index The index.
+     * @returns The step size.
+     */
+    step(index) {
+        return index < 3 ? this._posStep : this._velStep;
+    }
+}
+//# sourceMappingURL=PropagatorPairs.js.map

package/dist/observation/PropagatorPairs.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"PropagatorPairs.js","sourceRoot":"","sources":["../../src/observation/PropagatorPairs.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAIH,MAAM,OAAO,eAAe;IACN;IAA0B;IAA9C,YAAoB,QAAgB,EAAU,QAAgB;QAA1C,aAAQ,GAAR,QAAQ,CAAQ;QAAU,aAAQ,GAAR,QAAQ,CAAQ;QAC5D,cAAc;IAChB,CAAC;IAEO,KAAK,GAA0B,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACnD,IAAI,GAA0B,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAE1D,GAAG,CAAC,KAAa,EAAE,IAAgB,EAAE,GAAe;QAClD,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC;QACzB,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC;IACzB,CAAC;IAED,GAAG,CAAC,KAAa;QACf,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAE,EAAE,IAAI,CAAC,IAAI,CAAC,KAAK,CAAE,CAAC,CAAC;IACjD,CAAC;IAED;;;;OAIG;IACH,IAAI,CAAC,KAAa;QAChB,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC;IACnD,CAAC;CACF"}