ootk 5.1.1 → 6.0.0-rc.2

package/dist/conjunction/ConjunctionAssessment.js

@@ -0,0 +1,281 @@
+/**
+ * @author Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * 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 '../time/EpochUTC.js';
+import { RIC } from '../coordinate/RIC.js';
+import { StateCovariance, CovarianceFrame } from '../covariance/StateCovariance.js';
+import { CovarianceSample } from '../covariance/CovarianceSample.js';
+import { Sgp4Propagator } from '../propagator/Sgp4Propagator.js';
+import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator.js';
+import { ForceModel } from '../force/ForceModel.js';
+import { GoldenSection } from '../optimize/GoldenSection.js';
+import { Matrix } from '../operations/Matrix.js';
+import { ConjunctionEvent } from './ConjunctionEvent.js';
+import { ProbabilityOfCollision } from './ProbabilityOfCollision.js';
+/**
+ * High accuracy conjunction assessment workflow.
+ *
+ * This class provides a comprehensive workflow for assessing conjunctions between
+ * space objects using:
+ * - High accuracy propagators (SGP4 or numerical integrators)
+ * - Covariance matrices based on historical TLE accuracy
+ * - Time of Closest Approach (TCA) finding using optimization
+ * - Probability of collision calculation using Chan's 2D method
+ *
+ * @example
+ * ```typescript
+ * const primaryTle = new Tle(line1, line2);
+ * const secondaryTle = new Tle(line1, line2);
+ *
+ * const assessment = new ConjunctionAssessment(
+ *   { tle: primaryTle, radius: 0.01 as Kilometers },
+ *   { tle: secondaryTle, radius: 0.01 as Kilometers },
+ * );
+ *
+ * const event = assessment.assess({
+ *   startTime: EpochUTC.fromDateTime(new Date('2025-01-01T00:00:00Z')),
+ *   endTime: EpochUTC.fromDateTime(new Date('2025-01-02T00:00:00Z')),
+ *   useHighFidelity: true,
+ *   propagateCovariance: true,
+ * });
+ *
+ * console.log(event.toString());
+ * ```
+ */
+export class ConjunctionAssessment {
+    primary;
+    secondary;
+    primaryProp;
+    secondaryProp;
+    primaryCovSample;
+    secondaryCovSample;
+    constructor(primary, secondary) {
+        this.primary = primary;
+        this.secondary = secondary;
+        // Initialize propagators (will be replaced in assess() if needed)
+        this.primaryProp = this.createPropagator(primary, false);
+        this.secondaryProp = this.createPropagator(secondary, false);
+    }
+    /**
+     * Performs conjunction assessment over the specified time window.
+     *
+     * @param options Assessment configuration options
+     * @returns ConjunctionEvent with TCA, miss distance, and probability of collision
+     */
+    assess(options) {
+        const { startTime, endTime, useHighFidelity = false, forceModel, propagateCovariance = false, tcaTolerance = 0.001, searchStepSize = 60.0, } = options;
+        // Create propagators based on options
+        this.primaryProp = this.primary.propagator ?? this.createPropagator(this.primary, useHighFidelity, forceModel);
+        this.secondaryProp =
+            this.secondary.propagator ?? this.createPropagator(this.secondary, useHighFidelity, forceModel);
+        // Initialize covariance samples if requested
+        if (propagateCovariance) {
+            this.initializeCovarianceSamples(useHighFidelity, forceModel);
+        }
+        // Find Time of Closest Approach (TCA)
+        const tca = this.findTCA(startTime, endTime, searchStepSize, tcaTolerance);
+        // Propagate states to TCA
+        const primaryState = this.primaryProp.propagate(tca);
+        const secondaryState = this.secondaryProp.propagate(tca);
+        // Compute relative state in RIC frame
+        const relativeState = RIC.fromJ2000(secondaryState, primaryState);
+        // Extract RIC components
+        const radialDistance = Math.abs(relativeState.position.x);
+        const intrackDistance = Math.abs(relativeState.position.y);
+        const crosstrackDistance = Math.abs(relativeState.position.z);
+        const missDistance = relativeState.range;
+        const relativeVelocity = relativeState.velocity.magnitude();
+        // Compute combined covariance and Pc if available
+        let combinedCovariance;
+        let probabilityOfCollision;
+        if (this.primaryCovSample && this.secondaryCovSample && this.primary.radius && this.secondary.radius) {
+            // Propagate covariance samples to TCA
+            this.primaryCovSample.propagate(tca);
+            this.secondaryCovSample.propagate(tca);
+            // Get covariances in RIC frame
+            const primaryCov = this.primaryCovSample.covarianceRIC();
+            const secondaryCov = this.secondaryCovSample.covarianceRIC();
+            // Combine covariances
+            combinedCovariance = ProbabilityOfCollision.combineCovarianceMatrices(primaryCov, secondaryCov);
+            // Calculate probability of collision
+            const combinedRadius = (this.primary.radius + this.secondary.radius);
+            probabilityOfCollision = ProbabilityOfCollision.calculate(relativeState.position, relativeState.velocity, combinedCovariance, combinedRadius);
+        }
+        else if (this.primary.covariance && this.secondary.covariance && this.primary.radius && this.secondary.radius) {
+            // Use provided covariances without propagation
+            // Transform to RIC if needed
+            let primaryCovRIC = this.primary.covariance;
+            let secondaryCovRIC = this.secondary.covariance;
+            if (this.primary.covariance.frame === CovarianceFrame.ECI) {
+                primaryCovRIC = this.transformCovarianceToRIC(this.primary.covariance, primaryState);
+            }
+            if (this.secondary.covariance.frame === CovarianceFrame.ECI) {
+                secondaryCovRIC = this.transformCovarianceToRIC(this.secondary.covariance, secondaryState);
+            }
+            combinedCovariance = ProbabilityOfCollision.combineCovarianceMatrices(primaryCovRIC, secondaryCovRIC);
+            const combinedRadius = (this.primary.radius + this.secondary.radius);
+            probabilityOfCollision = ProbabilityOfCollision.calculate(relativeState.position, relativeState.velocity, combinedCovariance, combinedRadius);
+        }
+        return new ConjunctionEvent(tca, primaryState, secondaryState, relativeState, missDistance, radialDistance, intrackDistance, crosstrackDistance, relativeVelocity, combinedCovariance, probabilityOfCollision, this.primary.radius, this.secondary.radius);
+    }
+    /**
+     * Finds the Time of Closest Approach (TCA) using golden section search.
+     *
+     * @param startTime Search window start
+     * @param endTime Search window end
+     * @param stepSize Initial search step size in seconds
+     * @param tolerance TCA search tolerance in seconds
+     * @returns TCA epoch
+     */
+    findTCA(startTime, endTime, stepSize, tolerance) {
+        // Coarse search to find approximate TCA region
+        let minRange = Infinity;
+        let minEpoch = startTime;
+        let current = startTime;
+        while (current.posix <= endTime.posix) {
+            const primary = this.primaryProp.propagate(current);
+            const secondary = this.secondaryProp.propagate(current);
+            const ric = RIC.fromJ2000(secondary, primary);
+            const range = ric.range;
+            if (range < minRange) {
+                minRange = range;
+                minEpoch = current;
+            }
+            current = current.roll(stepSize);
+        }
+        // Fine search using golden section optimization
+        const searchWindow = 2 * stepSize;
+        const lowerBound = Math.max(startTime.posix, minEpoch.posix - searchWindow);
+        const upperBound = Math.min(endTime.posix, minEpoch.posix + searchWindow);
+        const tcaPosix = GoldenSection.search((posix) => {
+            const epoch = new EpochUTC(posix);
+            const primary = this.primaryProp.propagate(epoch);
+            const secondary = this.secondaryProp.propagate(epoch);
+            const ric = RIC.fromJ2000(secondary, primary);
+            return ric.range;
+        }, lowerBound, upperBound, { tolerance });
+        return new EpochUTC(tcaPosix);
+    }
+    /**
+     * Creates a propagator for a space object.
+     *
+     * @param obj Space object
+     * @param useHighFidelity Whether to use high-fidelity propagation
+     * @param forceModel Optional force model for numerical propagation
+     * @returns Propagator instance
+     */
+    createPropagator(obj, useHighFidelity, forceModel) {
+        if (obj.tle) {
+            if (useHighFidelity) {
+                // Convert TLE to state and use RK89
+                const state = obj.tle.toJ2000();
+                const fm = forceModel ?? new ForceModel().setGravity();
+                return new RungeKutta89Propagator(state, fm);
+            }
+            return new Sgp4Propagator(obj.tle);
+        }
+        else if (obj.state) {
+            if (useHighFidelity) {
+                const fm = forceModel ?? new ForceModel().setGravity();
+                return new RungeKutta89Propagator(obj.state, fm);
+            }
+            // For low-fidelity, still use RK89 but with simpler force model
+            const fm = new ForceModel().setGravity();
+            return new RungeKutta89Propagator(obj.state, fm);
+        }
+        throw new Error('Space object must have either a TLE or state vector');
+    }
+    /**
+     * Initializes covariance samples for both objects.
+     *
+     * @param useHighFidelity Whether to use high-fidelity propagation
+     * @param forceModel Optional force model
+     */
+    initializeCovarianceSamples(useHighFidelity, forceModel) {
+        const fm = forceModel ?? new ForceModel().setGravity();
+        // Primary covariance
+        if (this.primary.tle) {
+            const state = this.primary.tle.toJ2000();
+            const covariance = this.primary.covariance ??
+                StateCovariance.fromSigmas([1.0, 1.0, 1.0, 0.001, 0.001, 0.001], CovarianceFrame.RIC);
+            this.primaryCovSample = new CovarianceSample(state, covariance, this.primary.tle, fm, fm);
+        }
+        else if (this.primary.state && this.primary.covariance) {
+            this.primaryCovSample = new CovarianceSample(this.primary.state, this.primary.covariance, undefined, fm, fm);
+        }
+        // Secondary covariance
+        if (this.secondary.tle) {
+            const state = this.secondary.tle.toJ2000();
+            const covariance = this.secondary.covariance ??
+                StateCovariance.fromSigmas([1.0, 1.0, 1.0, 0.001, 0.001, 0.001], CovarianceFrame.RIC);
+            this.secondaryCovSample = new CovarianceSample(state, covariance, this.secondary.tle, fm, fm);
+        }
+        else if (this.secondary.state && this.secondary.covariance) {
+            this.secondaryCovSample = new CovarianceSample(this.secondary.state, this.secondary.covariance, undefined, fm, fm);
+        }
+    }
+    /**
+     * Transforms ECI covariance to RIC frame.
+     *
+     * @param covariance ECI covariance
+     * @param state State vector for RIC frame definition
+     * @returns RIC covariance
+     */
+    transformCovarianceToRIC(covariance, state) {
+        // Create RIC transformation matrix (3x3)
+        const ricMatrix = this.createRICTransformMatrix(state.position, state.velocity);
+        // Build 6x6 transformation matrix (block diagonal with 3x3 rotation)
+        const transform = this.build6x6Transform(ricMatrix);
+        // Transform covariance: C_ric = T * C_eci * T^T
+        const covRIC = transform.multiply(covariance.matrix).multiply(transform.transpose());
+        return new StateCovariance(covRIC, CovarianceFrame.RIC);
+    }
+    /**
+     * Creates the 3x3 RIC transformation matrix.
+     *
+     * @param position Position vector
+     * @param velocity Velocity vector
+     * @returns 3x3 RIC transformation matrix
+     */
+    createRICTransformMatrix(position, velocity) {
+        const ru = position.normalize();
+        const cu = position.cross(velocity).normalize();
+        const iu = cu.cross(ru).normalize();
+        return new Matrix([
+            [ru.x, ru.y, ru.z],
+            [iu.x, iu.y, iu.z],
+            [cu.x, cu.y, cu.z],
+        ]);
+    }
+    /**
+     * Builds a 6x6 transformation matrix from a 3x3 rotation matrix.
+     *
+     * @param rot 3x3 rotation matrix
+     * @returns 6x6 transformation matrix
+     */
+    build6x6Transform(rot) {
+        const elements = [
+            [rot.elements[0][0], rot.elements[0][1], rot.elements[0][2], 0, 0, 0],
+            [rot.elements[1][0], rot.elements[1][1], rot.elements[1][2], 0, 0, 0],
+            [rot.elements[2][0], rot.elements[2][1], rot.elements[2][2], 0, 0, 0],
+            [0, 0, 0, rot.elements[0][0], rot.elements[0][1], rot.elements[0][2]],
+            [0, 0, 0, rot.elements[1][0], rot.elements[1][1], rot.elements[1][2]],
+            [0, 0, 0, rot.elements[2][0], rot.elements[2][1], rot.elements[2][2]],
+        ];
+        return new Matrix(elements);
+    }
+}
+//# sourceMappingURL=ConjunctionAssessment.js.map

package/dist/conjunction/ConjunctionAssessment.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ConjunctionAssessment.js","sourceRoot":"","sources":["../../src/conjunction/ConjunctionAssessment.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAE/C,OAAO,EAAE,GAAG,EAAE,MAAM,sBAAsB,CAAC;AAE3C,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,kCAAkC,CAAC;AACpF,OAAO,EAAE,gBAAgB,EAAE,MAAM,mCAAmC,CAAC;AAErE,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,sBAAsB,EAAE,MAAM,yCAAyC,CAAC;AACjF,OAAO,EAAE,UAAU,EAAE,MAAM,wBAAwB,CAAC;AACpD,OAAO,EAAE,aAAa,EAAE,MAAM,8BAA8B,CAAC;AAC7D,OAAO,EAAE,MAAM,EAAE,MAAM,yBAAyB,CAAC;AAEjD,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,OAAO,EAAE,sBAAsB,EAAE,MAAM,6BAA6B,CAAC;AAoDrE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6BG;AACH,MAAM,OAAO,qBAAqB;IAOtB;IACA;IAPF,WAAW,CAAa;IACxB,aAAa,CAAa;IAC1B,gBAAgB,CAAoB;IACpC,kBAAkB,CAAoB;IAE9C,YACU,OAAoB,EACpB,SAAsB;QADtB,YAAO,GAAP,OAAO,CAAa;QACpB,cAAS,GAAT,SAAS,CAAa;QAE9B,kEAAkE;QAClE,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,gBAAgB,CAAC,OAAO,EAAE,KAAK,CAAC,CAAC;QACzD,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,gBAAgB,CAAC,SAAS,EAAE,KAAK,CAAC,CAAC;IAC/D,CAAC;IAED;;;;;OAKG;IACH,MAAM,CAAC,OAAqC;QAC1C,MAAM,EACJ,SAAS,EACT,OAAO,EACP,eAAe,GAAG,KAAK,EACvB,UAAU,EACV,mBAAmB,GAAG,KAAK,EAC3B,YAAY,GAAG,KAAgB,EAC/B,cAAc,GAAG,IAAe,GACjC,GAAG,OAAO,CAAC;QAEZ,sCAAsC;QACtC,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,IAAI,IAAI,CAAC,gBAAgB,CAAC,IAAI,CAAC,OAAO,EAAE,eAAe,EAAE,UAAU,CAAC,CAAC;QAC/G,IAAI,CAAC,aAAa;YAChB,IAAI,CAAC,SAAS,CAAC,UAAU,IAAI,IAAI,CAAC,gBAAgB,CAAC,IAAI,CAAC,SAAS,EAAE,eAAe,EAAE,UAAU,CAAC,CAAC;QAElG,6CAA6C;QAC7C,IAAI,mBAAmB,EAAE,CAAC;YACxB,IAAI,CAAC,2BAA2B,CAAC,eAAe,EAAE,UAAU,CAAC,CAAC;QAChE,CAAC;QAED,sCAAsC;QACtC,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,SAAS,EAAE,OAAO,EAAE,cAAc,EAAE,YAAY,CAAC,CAAC;QAE3E,0BAA0B;QAC1B,MAAM,YAAY,GAAG,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QACrD,MAAM,cAAc,GAAG,IAAI,CAAC,aAAa,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;QAEzD,sCAAsC;QACtC,MAAM,aAAa,GAAG,GAAG,CAAC,SAAS,CAAC,cAAc,EAAE,YAAY,CAAC,CAAC;QAElE,yBAAyB;QACzB,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC,CAAe,CAAC;QACxE,MAAM,eAAe,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC,CAAe,CAAC;QACzE,MAAM,kBAAkB,GAAG,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC,CAAe,CAAC;QAC5E,MAAM,YAAY,GAAG,aAAa,CAAC,KAAK,CAAC;QACzC,MAAM,gBAAgB,GAAG,aAAa,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAE5D,kDAAkD;QAClD,IAAI,kBAA+C,CAAC;QACpD,IAAI,sBAA0C,CAAC;QAE/C,IAAI,IAAI,CAAC,gBAAgB,IAAI,IAAI,CAAC,kBAAkB,IAAI,IAAI,CAAC,OAAO,CAAC,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,CAAC;YACrG,sCAAsC;YACtC,IAAI,CAAC,gBAAgB,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;YACrC,IAAI,CAAC,kBAAkB,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;YAEvC,+BAA+B;YAC/B,MAAM,UAAU,GAAG,IAAI,CAAC,gBAAgB,CAAC,aAAa,EAAE,CAAC;YACzD,MAAM,YAAY,GAAG,IAAI,CAAC,kBAAkB,CAAC,aAAa,EAAE,CAAC;YAE7D,sBAAsB;YACtB,kBAAkB,GAAG,sBAAsB,CAAC,yBAAyB,CAAC,UAAU,EAAE,YAAY,CAAC,CAAC;YAEhG,qCAAqC;YACrC,MAAM,cAAc,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAe,CAAC;YAEnF,sBAAsB,GAAG,sBAAsB,CAAC,SAAS,CACvD,aAAa,CAAC,QAAQ,EACtB,aAAa,CAAC,QAAQ,EACtB,kBAAkB,EAClB,cAAc,CACf,CAAC;QACJ,CAAC;aAAM,IAAI,IAAI,CAAC,OAAO,CAAC,UAAU,IAAI,IAAI,CAAC,SAAS,CAAC,UAAU,IAAI,IAAI,CAAC,OAAO,CAAC,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,CAAC;YAChH,+CAA+C;YAC/C,6BAA6B;YAC7B,IAAI,aAAa,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC;YAC5C,IAAI,eAAe,GAAG,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC;YAEhD,IAAI,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,KAAK,KAAK,eAAe,CAAC,GAAG,EAAE,CAAC;gBAC1D,aAAa,GAAG,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,OAAO,CAAC,UAAU,EAAE,YAAY,CAAC,CAAC;YACvF,CAAC;YACD,IAAI,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,KAAK,KAAK,eAAe,CAAC,GAAG,EAAE,CAAC;gBAC5D,eAAe,GAAG,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,SAAS,CAAC,UAAU,EAAE,cAAc,CAAC,CAAC;YAC7F,CAAC;YAED,kBAAkB,GAAG,sBAAsB,CAAC,yBAAyB,CAAC,aAAa,EAAE,eAAe,CAAC,CAAC;YAEtG,MAAM,cAAc,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAe,CAAC;YAEnF,sBAAsB,GAAG,sBAAsB,CAAC,SAAS,CACvD,aAAa,CAAC,QAAQ,EACtB,aAAa,CAAC,QAAQ,EACtB,kBAAkB,EAClB,cAAc,CACf,CAAC;QACJ,CAAC;QAED,OAAO,IAAI,gBAAgB,CACzB,GAAG,EACH,YAAY,EACZ,cAAc,EACd,aAAa,EACb,YAAY,EACZ,cAAc,EACd,eAAe,EACf,kBAAkB,EAClB,gBAAgB,EAChB,kBAAkB,EAClB,sBAAsB,EACtB,IAAI,CAAC,OAAO,CAAC,MAAM,EACnB,IAAI,CAAC,SAAS,CAAC,MAAM,CACtB,CAAC;IACJ,CAAC;IAED;;;;;;;;OAQG;IACK,OAAO,CAAC,SAAmB,EAAE,OAAiB,EAAE,QAAiB,EAAE,SAAkB;QAC3F,+CAA+C;QAC/C,IAAI,QAAQ,GAAG,QAAQ,CAAC;QACxB,IAAI,QAAQ,GAAG,SAAS,CAAC;QACzB,IAAI,OAAO,GAAG,SAAS,CAAC;QAExB,OAAO,OAAO,CAAC,KAAK,IAAI,OAAO,CAAC,KAAK,EAAE,CAAC;YACtC,MAAM,OAAO,GAAG,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;YACpD,MAAM,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;YACxD,MAAM,GAAG,GAAG,GAAG,CAAC,SAAS,CAAC,SAAS,EAAE,OAAO,CAAC,CAAC;YAC9C,MAAM,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;YAExB,IAAI,KAAK,GAAG,QAAQ,EAAE,CAAC;gBACrB,QAAQ,GAAG,KAAK,CAAC;gBACjB,QAAQ,GAAG,OAAO,CAAC;YACrB,CAAC;YAED,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACnC,CAAC;QAED,gDAAgD;QAChD,MAAM,YAAY,GAAG,CAAC,GAAG,QAAQ,CAAC;QAClC,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,KAAK,EAAE,QAAQ,CAAC,KAAK,GAAG,YAAY,CAAC,CAAC;QAC5E,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,EAAE,QAAQ,CAAC,KAAK,GAAG,YAAY,CAAC,CAAC;QAE1E,MAAM,QAAQ,GAAG,aAAa,CAAC,MAAM,CACnC,CAAC,KAAK,EAAE,EAAE;YACR,MAAM,KAAK,GAAG,IAAI,QAAQ,CAAC,KAAgB,CAAC,CAAC;YAC7C,MAAM,OAAO,GAAG,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;YAClD,MAAM,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;YACtD,MAAM,GAAG,GAAG,GAAG,CAAC,SAAS,CAAC,SAAS,EAAE,OAAO,CAAC,CAAC;YAE9C,OAAO,GAAG,CAAC,KAAK,CAAC;QACnB,CAAC,EACD,UAAU,EACV,UAAU,EACV,EAAE,SAAS,EAAE,CACd,CAAC;QAEF,OAAO,IAAI,QAAQ,CAAC,QAAmB,CAAC,CAAC;IAC3C,CAAC;IAED;;;;;;;OAOG;IACK,gBAAgB,CAAC,GAAgB,EAAE,eAAwB,EAAE,UAAuB;QAC1F,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC;YACZ,IAAI,eAAe,EAAE,CAAC;gBACpB,oCAAoC;gBACpC,MAAM,KAAK,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC;gBAChC,MAAM,EAAE,GAAG,UAAU,IAAI,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;gBAEvD,OAAO,IAAI,sBAAsB,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;YAC/C,CAAC;YAED,OAAO,IAAI,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACrC,CAAC;aAAM,IAAI,GAAG,CAAC,KAAK,EAAE,CAAC;YACrB,IAAI,eAAe,EAAE,CAAC;gBACpB,MAAM,EAAE,GAAG,UAAU,IAAI,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;gBAEvD,OAAO,IAAI,sBAAsB,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;YACnD,CAAC;YAED,gEAAgE;YAChE,MAAM,EAAE,GAAG,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;YAEzC,OAAO,IAAI,sBAAsB,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACnD,CAAC;QAED,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,CAAC;IACzE,CAAC;IAED;;;;;OAKG;IACK,2BAA2B,CAAC,eAAwB,EAAE,UAAuB;QACnF,MAAM,EAAE,GAAG,UAAU,IAAI,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;QAEvD,qBAAqB;QACrB,IAAI,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,CAAC;YACrB,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC;YACzC,MAAM,UAAU,GACd,IAAI,CAAC,OAAO,CAAC,UAAU;gBACvB,eAAe,CAAC,UAAU,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;YAExF,IAAI,CAAC,gBAAgB,GAAG,IAAI,gBAAgB,CAAC,KAAK,EAAE,UAAU,EAAE,IAAI,CAAC,OAAO,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC5F,CAAC;aAAM,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,IAAI,IAAI,CAAC,OAAO,CAAC,UAAU,EAAE,CAAC;YACzD,IAAI,CAAC,gBAAgB,GAAG,IAAI,gBAAgB,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,IAAI,CAAC,OAAO,CAAC,UAAU,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC/G,CAAC;QAED,uBAAuB;QACvB,IAAI,IAAI,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC;YACvB,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,OAAO,EAAE,CAAC;YAC3C,MAAM,UAAU,GACd,IAAI,CAAC,SAAS,CAAC,UAAU;gBACzB,eAAe,CAAC,UAAU,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;YAExF,IAAI,CAAC,kBAAkB,GAAG,IAAI,gBAAgB,CAAC,KAAK,EAAE,UAAU,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAChG,CAAC;aAAM,IAAI,IAAI,CAAC,SAAS,CAAC,KAAK,IAAI,IAAI,CAAC,SAAS,CAAC,UAAU,EAAE,CAAC;YAC7D,IAAI,CAAC,kBAAkB,GAAG,IAAI,gBAAgB,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,IAAI,CAAC,SAAS,CAAC,UAAU,EAAE,SAAS,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QACrH,CAAC;IACH,CAAC;IAED;;;;;;OAMG;IACK,wBAAwB,CAAC,UAA2B,EAAE,KAAY;QACxE,yCAAyC;QACzC,MAAM,SAAS,GAAG,IAAI,CAAC,wBAAwB,CAAC,KAAK,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC,CAAC;QAEhF,qEAAqE;QACrE,MAAM,SAAS,GAAG,IAAI,CAAC,iBAAiB,CAAC,SAAS,CAAC,CAAC;QAEpD,gDAAgD;QAChD,MAAM,MAAM,GAAG,SAAS,CAAC,QAAQ,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,QAAQ,CAAC,SAAS,CAAC,SAAS,EAAE,CAAC,CAAC;QAErF,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;IAC1D,CAAC;IAED;;;;;;OAMG;IACK,wBAAwB,CAAC,QAAa,EAAE,QAAa;QAC3D,MAAM,EAAE,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;QAChC,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,SAAS,EAAE,CAAC;QAChD,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,SAAS,EAAE,CAAC;QAEpC,OAAO,IAAI,MAAM,CAAC;YAChB,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;YAClB,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;YAClB,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;SACnB,CAAC,CAAC;IACL,CAAC;IAED;;;;;OAKG;IACK,iBAAiB,CAAC,GAAW;QACnC,MAAM,QAAQ,GAAG;YACf,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACrE,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACrE,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;YACrE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SACtE,CAAC;QAEF,OAAO,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC;IAC9B,CAAC;CACF"}

package/dist/conjunction/ConjunctionEvent.d.ts

@@ -0,0 +1,104 @@
+/**
+ * @author Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * 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 type { EpochUTC } from '../time/EpochUTC.js';
+import type { J2000 } from '../coordinate/J2000.js';
+import type { RIC } from '../coordinate/RIC.js';
+import type { StateCovariance } from '../covariance/StateCovariance.js';
+import type { Kilometers, KilometersPerSecond } from '../main.js';
+/**
+ * Represents the result of a conjunction assessment between two space objects.
+ * Contains all relevant information about the close approach event.
+ */
+export declare class ConjunctionEvent {
+    /** Time of Closest Approach (TCA) */
+    tca: EpochUTC;
+    /** Primary object state at TCA in J2000 frame */
+    primaryState: J2000;
+    /** Secondary object state at TCA in J2000 frame */
+    secondaryState: J2000;
+    /** Relative state in RIC frame (relative to primary) */
+    relativeState: RIC;
+    /** Total miss distance at TCA (km) */
+    missDistance: Kilometers;
+    /** Radial component of miss distance (km) */
+    radialDistance: Kilometers;
+    /** Intrack component of miss distance (km) */
+    intrackDistance: Kilometers;
+    /** Crosstrack component of miss distance (km) */
+    crosstrackDistance: Kilometers;
+    /** Relative velocity magnitude at TCA (km/s) */
+    relativeVelocity: KilometersPerSecond;
+    /** Combined position covariance matrix in RIC frame (optional) */
+    combinedCovariance?: StateCovariance | undefined;
+    /** Probability of collision (optional, 0-1) */
+    probabilityOfCollision?: number | undefined;
+    /** Hard body radius for primary object (km, optional) */
+    primaryRadius?: Kilometers | undefined;
+    /** Hard body radius for secondary object (km, optional) */
+    secondaryRadius?: Kilometers | undefined;
+    constructor(
+    /** Time of Closest Approach (TCA) */
+    tca: EpochUTC, 
+    /** Primary object state at TCA in J2000 frame */
+    primaryState: J2000, 
+    /** Secondary object state at TCA in J2000 frame */
+    secondaryState: J2000, 
+    /** Relative state in RIC frame (relative to primary) */
+    relativeState: RIC, 
+    /** Total miss distance at TCA (km) */
+    missDistance: Kilometers, 
+    /** Radial component of miss distance (km) */
+    radialDistance: Kilometers, 
+    /** Intrack component of miss distance (km) */
+    intrackDistance: Kilometers, 
+    /** Crosstrack component of miss distance (km) */
+    crosstrackDistance: Kilometers, 
+    /** Relative velocity magnitude at TCA (km/s) */
+    relativeVelocity: KilometersPerSecond, 
+    /** Combined position covariance matrix in RIC frame (optional) */
+    combinedCovariance?: StateCovariance | undefined, 
+    /** Probability of collision (optional, 0-1) */
+    probabilityOfCollision?: number | undefined, 
+    /** Hard body radius for primary object (km, optional) */
+    primaryRadius?: Kilometers | undefined, 
+    /** Hard body radius for secondary object (km, optional) */
+    secondaryRadius?: Kilometers | undefined);
+    /**
+     * Returns a formatted string representation of the conjunction event.
+     * @returns A multi-line string with conjunction details.
+     */
+    toString(): string;
+    /**
+     * Checks if this is a high-risk conjunction based on miss distance and Pc.
+     * @param distanceThreshold Miss distance threshold in km (default: 1.0 km)
+     * @param pcThreshold Probability of collision threshold (default: 1e-4)
+     * @returns True if the conjunction exceeds risk thresholds.
+     */
+    isHighRisk(distanceThreshold?: Kilometers, pcThreshold?: number): boolean;
+    /**
+     * Gets the Mahalanobis distance if covariance is available.
+     * This is the miss distance normalized by the combined covariance.
+     * @returns The Mahalanobis distance (unitless), or undefined if no covariance.
+     */
+    getMahalanobisDistance(): number | undefined;
+    /**
+     * Extracts the 3x3 position covariance from a 6x6 state covariance matrix.
+     * @param stateCov 6x6 state covariance matrix
+     * @returns 3x3 position covariance matrix
+     */
+    private extractPositionCovariance;
+}

package/dist/conjunction/ConjunctionEvent.js

@@ -0,0 +1,156 @@
+/**
+ * @author Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * 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 '../operations/Matrix.js';
+/**
+ * Represents the result of a conjunction assessment between two space objects.
+ * Contains all relevant information about the close approach event.
+ */
+export class ConjunctionEvent {
+    tca;
+    primaryState;
+    secondaryState;
+    relativeState;
+    missDistance;
+    radialDistance;
+    intrackDistance;
+    crosstrackDistance;
+    relativeVelocity;
+    combinedCovariance;
+    probabilityOfCollision;
+    primaryRadius;
+    secondaryRadius;
+    constructor(
+    /** Time of Closest Approach (TCA) */
+    tca, 
+    /** Primary object state at TCA in J2000 frame */
+    primaryState, 
+    /** Secondary object state at TCA in J2000 frame */
+    secondaryState, 
+    /** Relative state in RIC frame (relative to primary) */
+    relativeState, 
+    /** Total miss distance at TCA (km) */
+    missDistance, 
+    /** Radial component of miss distance (km) */
+    radialDistance, 
+    /** Intrack component of miss distance (km) */
+    intrackDistance, 
+    /** Crosstrack component of miss distance (km) */
+    crosstrackDistance, 
+    /** Relative velocity magnitude at TCA (km/s) */
+    relativeVelocity, 
+    /** Combined position covariance matrix in RIC frame (optional) */
+    combinedCovariance, 
+    /** Probability of collision (optional, 0-1) */
+    probabilityOfCollision, 
+    /** Hard body radius for primary object (km, optional) */
+    primaryRadius, 
+    /** Hard body radius for secondary object (km, optional) */
+    secondaryRadius) {
+        this.tca = tca;
+        this.primaryState = primaryState;
+        this.secondaryState = secondaryState;
+        this.relativeState = relativeState;
+        this.missDistance = missDistance;
+        this.radialDistance = radialDistance;
+        this.intrackDistance = intrackDistance;
+        this.crosstrackDistance = crosstrackDistance;
+        this.relativeVelocity = relativeVelocity;
+        this.combinedCovariance = combinedCovariance;
+        this.probabilityOfCollision = probabilityOfCollision;
+        this.primaryRadius = primaryRadius;
+        this.secondaryRadius = secondaryRadius;
+    }
+    /**
+     * Returns a formatted string representation of the conjunction event.
+     * @returns A multi-line string with conjunction details.
+     */
+    toString() {
+        const lines = [
+            '[Conjunction Event]',
+            `  TCA: ${this.tca.toString()}`,
+            `  Miss Distance: ${this.missDistance.toFixed(6)} km`,
+            `    Radial:     ${this.radialDistance.toFixed(6)} km`,
+            `    Intrack:    ${this.intrackDistance.toFixed(6)} km`,
+            `    Crosstrack: ${this.crosstrackDistance.toFixed(6)} km`,
+            `  Relative Velocity: ${this.relativeVelocity.toFixed(6)} km/s`,
+        ];
+        if (this.probabilityOfCollision !== undefined) {
+            lines.push(`  Probability of Collision: ${this.probabilityOfCollision.toExponential(6)}`);
+        }
+        if (this.primaryRadius !== undefined && this.secondaryRadius !== undefined) {
+            const combinedRadius = this.primaryRadius + this.secondaryRadius;
+            lines.push(`  Combined Hard Body Radius: ${combinedRadius.toFixed(3)} km`);
+        }
+        return lines.join('\n');
+    }
+    /**
+     * Checks if this is a high-risk conjunction based on miss distance and Pc.
+     * @param distanceThreshold Miss distance threshold in km (default: 1.0 km)
+     * @param pcThreshold Probability of collision threshold (default: 1e-4)
+     * @returns True if the conjunction exceeds risk thresholds.
+     */
+    isHighRisk(distanceThreshold = 1.0, pcThreshold = 1e-4) {
+        const distanceRisk = this.missDistance < distanceThreshold;
+        const pcRisk = this.probabilityOfCollision !== undefined && this.probabilityOfCollision > pcThreshold;
+        return distanceRisk || pcRisk;
+    }
+    /**
+     * Gets the Mahalanobis distance if covariance is available.
+     * This is the miss distance normalized by the combined covariance.
+     * @returns The Mahalanobis distance (unitless), or undefined if no covariance.
+     */
+    getMahalanobisDistance() {
+        if (!this.combinedCovariance) {
+            return undefined;
+        }
+        // Extract position-only covariance (first 3x3 block)
+        const posCovariance = this.extractPositionCovariance(this.combinedCovariance.matrix);
+        const relativePosition = [
+            this.relativeState.position.x,
+            this.relativeState.position.y,
+            this.relativeState.position.z,
+        ];
+        try {
+            // Compute Mahalanobis distance: sqrt(r^T * C^-1 * r)
+            const covInv = posCovariance.inverse();
+            const temp = covInv.multiplyVector(relativePosition);
+            let mahalanobis = 0;
+            for (let i = 0; i < 3; i++) {
+                mahalanobis += relativePosition[i] * temp[i];
+            }
+            return Math.sqrt(mahalanobis);
+        }
+        catch {
+            // Covariance may be singular
+            return undefined;
+        }
+    }
+    /**
+     * Extracts the 3x3 position covariance from a 6x6 state covariance matrix.
+     * @param stateCov 6x6 state covariance matrix
+     * @returns 3x3 position covariance matrix
+     */
+    extractPositionCovariance(stateCov) {
+        const elements = [
+            [stateCov.elements[0][0], stateCov.elements[0][1], stateCov.elements[0][2]],
+            [stateCov.elements[1][0], stateCov.elements[1][1], stateCov.elements[1][2]],
+            [stateCov.elements[2][0], stateCov.elements[2][1], stateCov.elements[2][2]],
+        ];
+        return new Matrix(elements);
+    }
+}
+//# sourceMappingURL=ConjunctionEvent.js.map

package/dist/conjunction/ConjunctionEvent.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ConjunctionEvent.js","sourceRoot":"","sources":["../../src/conjunction/ConjunctionEvent.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAOH,OAAO,EAAE,MAAM,EAAE,MAAM,yBAAyB,CAAC;AAEjD;;;GAGG;AACH,MAAM,OAAO,gBAAgB;IAGlB;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IAEA;IA1BT;IACE,qCAAqC;IAC9B,GAAa;IACpB,iDAAiD;IAC1C,YAAmB;IAC1B,mDAAmD;IAC5C,cAAqB;IAC5B,wDAAwD;IACjD,aAAkB;IACzB,sCAAsC;IAC/B,YAAwB;IAC/B,6CAA6C;IACtC,cAA0B;IACjC,8CAA8C;IACvC,eAA2B;IAClC,iDAAiD;IAC1C,kBAA8B;IACrC,gDAAgD;IACzC,gBAAqC;IAC5C,kEAAkE;IAC3D,kBAAoC;IAC3C,+CAA+C;IACxC,sBAA+B;IACtC,yDAAyD;IAClD,aAA0B;IACjC,2DAA2D;IACpD,eAA4B;QAxB5B,QAAG,GAAH,GAAG,CAAU;QAEb,iBAAY,GAAZ,YAAY,CAAO;QAEnB,mBAAc,GAAd,cAAc,CAAO;QAErB,kBAAa,GAAb,aAAa,CAAK;QAElB,iBAAY,GAAZ,YAAY,CAAY;QAExB,mBAAc,GAAd,cAAc,CAAY;QAE1B,oBAAe,GAAf,eAAe,CAAY;QAE3B,uBAAkB,GAAlB,kBAAkB,CAAY;QAE9B,qBAAgB,GAAhB,gBAAgB,CAAqB;QAErC,uBAAkB,GAAlB,kBAAkB,CAAkB;QAEpC,2BAAsB,GAAtB,sBAAsB,CAAS;QAE/B,kBAAa,GAAb,aAAa,CAAa;QAE1B,oBAAe,GAAf,eAAe,CAAa;IAClC,CAAC;IAEJ;;;OAGG;IACH,QAAQ;QACN,MAAM,KAAK,GAAG;YACZ,qBAAqB;YACrB,UAAU,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,EAAE;YAC/B,oBAAoB,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK;YACrD,mBAAmB,IAAI,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK;YACtD,mBAAmB,IAAI,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK;YACvD,mBAAmB,IAAI,CAAC,kBAAkB,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK;YAC1D,wBAAwB,IAAI,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO;SAChE,CAAC;QAEF,IAAI,IAAI,CAAC,sBAAsB,KAAK,SAAS,EAAE,CAAC;YAC9C,KAAK,CAAC,IAAI,CAAC,+BAA+B,IAAI,CAAC,sBAAsB,CAAC,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QAC5F,CAAC;QAED,IAAI,IAAI,CAAC,aAAa,KAAK,SAAS,IAAI,IAAI,CAAC,eAAe,KAAK,SAAS,EAAE,CAAC;YAC3E,MAAM,cAAc,GAAG,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,eAAe,CAAC;YAEjE,KAAK,CAAC,IAAI,CAAC,gCAAgC,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;QAC7E,CAAC;QAED,OAAO,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAC1B,CAAC;IAED;;;;;OAKG;IACH,UAAU,CAAC,oBAAgC,GAAiB,EAAE,cAAsB,IAAI;QACtF,MAAM,YAAY,GAAG,IAAI,CAAC,YAAY,GAAG,iBAAiB,CAAC;QAC3D,MAAM,MAAM,GAAG,IAAI,CAAC,sBAAsB,KAAK,SAAS,IAAI,IAAI,CAAC,sBAAsB,GAAG,WAAW,CAAC;QAEtG,OAAO,YAAY,IAAI,MAAM,CAAC;IAChC,CAAC;IAED;;;;OAIG;IACH,sBAAsB;QACpB,IAAI,CAAC,IAAI,CAAC,kBAAkB,EAAE,CAAC;YAC7B,OAAO,SAAS,CAAC;QACnB,CAAC;QAED,qDAAqD;QACrD,MAAM,aAAa,GAAG,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,kBAAkB,CAAC,MAAM,CAAC,CAAC;QACrF,MAAM,gBAAgB,GAAG;YACvB,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;YAC7B,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;YAC7B,IAAI,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;SAC9B,CAAC;QAEF,IAAI,CAAC;YACH,qDAAqD;YACrD,MAAM,MAAM,GAAG,aAAa,CAAC,OAAO,EAAE,CAAC;YACvC,MAAM,IAAI,GAAG,MAAM,CAAC,cAAc,CAAC,gBAAgB,CAAC,CAAC;YACrD,IAAI,WAAW,GAAG,CAAC,CAAC;YAEpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,WAAW,IAAI,gBAAgB,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YAC/C,CAAC;YAED,OAAO,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;QAChC,CAAC;QAAC,MAAM,CAAC;YACP,6BAA6B;YAC7B,OAAO,SAAS,CAAC;QACnB,CAAC;IACH,CAAC;IAED;;;;OAIG;IACK,yBAAyB,CAAC,QAAgB;QAChD,MAAM,QAAQ,GAAG;YACf,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3E,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3E,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC5E,CAAC;QAEF,OAAO,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC;IAC9B,CAAC;CACF"}

package/dist/conjunction/ProbabilityOfCollision.d.ts

@@ -0,0 +1,102 @@
+/**
+ * @author Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * 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 { Vector3D } from '../operations/Vector3D.js';
+import type { StateCovariance } from '../covariance/StateCovariance.js';
+import type { Kilometers } from '../main.js';
+/**
+ * Probability of Collision calculator using Chan's 2D method.
+ *
+ * This implementation projects the combined covariance matrix onto the
+ * encounter plane (B-plane) perpendicular to the relative velocity vector,
+ * then computes the probability that the relative position lies within
+ * the combined hard body radius.
+ *
+ * Reference: Chan, F. K. (2008). "Spacecraft Collision Probability"
+ */
+export declare class ProbabilityOfCollision {
+    /**
+     * Computes probability of collision using Chan's 2D method.
+     *
+     * @param relativePosition Relative position vector in RIC frame (km)
+     * @param relativeVelocity Relative velocity vector in RIC frame (km/s)
+     * @param combinedCovariance Combined position covariance in RIC frame (6x6)
+     * @param combinedRadius Combined hard body radius (km)
+     * @returns Probability of collision (0 to 1)
+     */
+    static calculate(relativePosition: Vector3D<Kilometers>, relativeVelocity: Vector3D, combinedCovariance: StateCovariance, combinedRadius: Kilometers): number;
+    /**
+     * Calculates 2D probability of collision in the encounter plane.
+     *
+     * Uses the analytical solution for 2D Gaussian probability within a circle.
+     *
+     * @param x X-coordinate in encounter plane (km)
+     * @param y Y-coordinate in encounter plane (km)
+     * @param cov2D 2x2 covariance matrix in encounter plane
+     * @param radius Combined hard body radius (km)
+     * @returns Probability of collision (0 to 1)
+     */
+    private static calculatePc2D;
+    /**
+     * Chan's analytical approximation for 2D Pc.
+     *
+     * Reference: Chan, F. K. (2008). "Spacecraft Collision Probability"
+     *
+     * @param mahalanobisDistance Mahalanobis distance (sqrt of d2)
+     * @param radius Combined hard body radius
+     * @param sigma1 Larger semi-axis of covariance ellipse
+     * @param sigma2 Smaller semi-axis of covariance ellipse
+     * @returns Probability of collision (0 to 1)
+     */
+    private static chanPc2D;
+    /**
+     * Approximates the complementary error function erfc(x).
+     *
+     * Uses Abramowitz and Stegun approximation (formula 7.1.26).
+     *
+     * @param x Input value
+     * @returns erfc(x)
+     */
+    private static approximateErfc;
+    /**
+     * Fallback 3D probability calculation for cases with very low relative velocity.
+     *
+     * Uses simple spherical approximation.
+     *
+     * @param relativePosition Relative position vector
+     * @param posCovariance 3x3 position covariance matrix
+     * @param combinedRadius Combined hard body radius
+     * @returns Probability of collision (0 to 1)
+     */
+    private static calculate3D;
+    /**
+     * Combines two covariance matrices (primary and secondary).
+     *
+     * The combined covariance is simply the sum of the two covariances,
+     * assuming they are independent.
+     *
+     * @param primaryCov Primary object covariance
+     * @param secondaryCov Secondary object covariance
+     * @returns Combined covariance
+     */
+    static combineCovarianceMatrices(primaryCov: StateCovariance, secondaryCov: StateCovariance): StateCovariance;
+    /**
+     * Extracts the 3x3 position covariance from a 6x6 state covariance matrix.
+     * @param stateCov 6x6 state covariance matrix
+     * @returns 3x3 position covariance matrix
+     */
+    private static extractPositionCovariance;
+}