ootk 6.0.6 → 6.0.7

package/dist/main.js

@@ -21,8 +21,6 @@ var ValidationError = class extends OotkError {
     this.value = value;
     this.name = "ValidationError";
   }
-  field;
-  value;
 };
 var ParseError = class extends OotkError {
   /**
@@ -37,8 +35,6 @@ var ParseError = class extends OotkError {
     this.line = line;
     this.name = "ParseError";
   }
-  format;
-  line;
 };
 var PropagationError = class extends OotkError {
   /**
@@ -51,7 +47,6 @@ var PropagationError = class extends OotkError {
     this.epoch = epoch;
     this.name = "PropagationError";
   }
-  epoch;
 };
 var OrbitDeterminationError = class extends OotkError {
   /**
@@ -64,7 +59,6 @@ var OrbitDeterminationError = class extends OotkError {
     this.algorithm = algorithm;
     this.name = "OrbitDeterminationError";
   }
-  algorithm;
 };
 
 // src/enums/AngularDiameterMethod.ts
@@ -282,7 +276,6 @@ var Epoch = class {
       throw new ValidationError("Epoch posix time must be a valid number", "posix", posix);
     }
   }
-  posix;
   toString() {
     return this.toDateTime().toISOString();
   }
@@ -306,7 +299,7 @@ var Epoch = class {
     const currentDateObj = this.toDateTime();
     const epochYear = currentDateObj.getUTCFullYear().toString().slice(2, 4);
     const epochDay = this.getDayOfYear_(currentDateObj);
-    const timeOfDay = (currentDateObj.getUTCHours() * 60 + currentDateObj.getUTCMinutes()) / 1440;
+    const timeOfDay = (currentDateObj.getUTCHours() * 3600 + currentDateObj.getUTCMinutes() * 60 + currentDateObj.getUTCSeconds() + currentDateObj.getUTCMilliseconds() / 1e3) / 86400;
     const epochDayStr = (epochDay + timeOfDay).toFixed(8).padStart(12, "0");
     return {
       epochYr: epochYear,
@@ -1826,8 +1819,6 @@ var EpochGPS = class _EpochGPS {
       throw new ValidationError("GPS seconds must be between 0 and 604799", "seconds", seconds);
     }
   }
-  week;
-  seconds;
   /** Cached GPS reference epoch (1980-01-06T00:00:00.000Z) */
   static reference_ = null;
   /**
@@ -2253,6 +2244,32 @@ function createSampleCovarianceFromTle(tleLine1, tleLine2, frame = "ric" /* RIC
   const sample = new CovarianceSample(j2000State, covariance2, tle);
   return frame === "ric" /* RIC */ ? sample.desampleRIC() : sample.desampleJ2000();
 }
+var DEFAULT_SCREENING_SIGMA_CAPS = [1200, 1e3, 5e3];
+function cappedScreeningCovarianceFromTle(tleLine1, tleLine2, confidenceLevel, caps = DEFAULT_SCREENING_SIGMA_CAPS) {
+  let elements;
+  try {
+    elements = createSampleCovarianceFromTle(tleLine1, tleLine2).matrix.elements;
+  } catch {
+    return buildFallbackScreeningCovariance_(caps);
+  }
+  const radialVar = elements[0][0];
+  const intrackVar = elements[1][1];
+  const crosstrackVar = elements[2][2];
+  elements[0][0] = Math.min(Math.sqrt(radialVar) * confidenceLevel, caps[0]);
+  elements[1][1] = Math.min(Math.sqrt(crosstrackVar) * confidenceLevel, caps[1]);
+  elements[2][2] = Math.min(Math.sqrt(intrackVar) * confidenceLevel, caps[2]);
+  if (!elements[0][0] || !elements[1][1] || !elements[2][2]) {
+    return buildFallbackScreeningCovariance_(caps);
+  }
+  return new StateCovariance(new Matrix(elements), "eci" /* ECI */);
+}
+function buildFallbackScreeningCovariance_(caps) {
+  const fallback = Array.from({ length: 6 }, () => new Array(6).fill(0));
+  fallback[0][0] = caps[0];
+  fallback[1][1] = caps[1];
+  fallback[2][2] = caps[2];
+  return new StateCovariance(new Matrix(fallback), "eci" /* ECI */);
+}
 
 // src/utils/jacobian.ts
 var jacobian = (f, m, x0, step = 1e-5) => {
@@ -2279,6 +2296,17 @@ function linearDistance(pos1, pos2) {
   return Math.sqrt((pos1.x - pos2.x) ** 2 + (pos1.y - pos2.y) ** 2 + (pos1.z - pos2.z) ** 2);
 }
 
+// src/utils/photometry.ts
+var RCS_VMAG_ESTIMATE_MIN = -5;
+var RCS_VMAG_ESTIMATE_MAX = 15;
+function estimateVmagFromRcs(rcs) {
+  if (typeof rcs !== "number" || !Number.isFinite(rcs) || rcs <= 0) {
+    return null;
+  }
+  const vmag = -1.3 - 2.5 * Math.log10(rcs);
+  return Math.min(Math.max(vmag, RCS_VMAG_ESTIMATE_MIN), RCS_VMAG_ESTIMATE_MAX);
+}
+
 // src/operations/Vector3D.ts
 var Vector3D = class _Vector3D {
   constructor(x, y, z) {
@@ -2286,9 +2314,6 @@ var Vector3D = class _Vector3D {
     this.y = y;
     this.z = z;
   }
-  x;
-  y;
-  z;
   /**
    * Create a new Vector3D object from the first three elements of a Vector
    * object.
@@ -2486,7 +2511,6 @@ var Vector = class _Vector {
     this.elements = elements;
     this.length = elements.length;
   }
-  elements;
   /**
    * The length of the vector.
    */
@@ -3267,7 +3291,7 @@ var Random = class {
   }
 };
 
-// ../../../node_modules/astronomy-engine/esm/astronomy.js
+// node_modules/astronomy-engine/esm/astronomy.js
 var C_AUDAY = 173.1446326846693;
 var KM_PER_AU2 = 14959787069098932e-8;
 var DEG2RAD2 = 0.017453292519943295;
@@ -8607,12 +8631,6 @@ var CubicSpline = class {
     this.p1 = p1;
     this.m1 = m1;
   }
-  t0;
-  p0;
-  m0;
-  t1;
-  p1;
-  m1;
   // / Interpolate position at the provided time [t] _(POSIX seconds)_.
   position_(t) {
     const t2 = t * t;
@@ -8688,7 +8706,6 @@ var CubicSplineInterpolator = class _CubicSplineInterpolator extends StateInterp
     super();
     this.splines_ = splines_;
   }
-  splines_;
   static fromEphemeris(ephemeris) {
     const splines = [];
     for (let i = 0; i < ephemeris.length - 1; i++) {
@@ -8888,7 +8905,6 @@ var VerletBlendInterpolator = class _VerletBlendInterpolator extends StateInterp
     super();
     this.ephemeris = ephemeris;
   }
-  ephemeris;
   get sizeBytes() {
     return 64 * 7 * this.ephemeris.length / 8;
   }
@@ -9876,7 +9892,10 @@ var FormatTle = class _FormatTle {
    */
   static createTle(tleParams) {
     const { inc, meanmo, rasc, argPe, meana, ecen, epochyr, epochday, intl } = tleParams;
-    const scc = Tle.convert6DigitToA5(tleParams.scc);
+    let scc = Tle.convert6DigitToA5(tleParams.scc);
+    if (/^\d{1,4}$/u.test(scc)) {
+      scc = scc.padStart(5, "0");
+    }
     const epochYrStr = String(epochyr).padStart(2, "0");
     const epochdayStr = Number(epochday).toFixed(8).padStart(12, "0");
     const incStr = _FormatTle.inclination(inc);
@@ -10389,7 +10408,7 @@ var Sgp4 = class _Sgp4 {
       error: 0 /* NO_ERROR */
     };
     const xpdotp = 1440 / (2 * PI);
-    satrec.satnum = tleLine1.substring(2, 7);
+    satrec.satnum = tleLine1.charAt(7) === "V" ? "0" : tleLine1.substring(2, 7);
     satrec.epochyr = parseInt(tleLine1.substring(18, 20));
     satrec.epochdays = parseFloat(tleLine1.substring(20, 32));
     satrec.ndot = parseFloat(tleLine1.substring(33, 43));
@@ -13556,11 +13575,6 @@ var Thrust = class {
       crosstrack * 1e-3
     );
   }
-  center;
-  radial;
-  intrack;
-  crosstrack;
-  durationRate;
   deltaV;
   get magnitude() {
     return this.deltaV.magnitude() * 1e3;
@@ -13599,9 +13613,6 @@ var Hill = class _Hill {
     this.semimajorAxis_ = semimajorAxis;
     this.meanMotion_ = Earth.smaToMeanMotion(this.semimajorAxis_);
   }
-  epoch;
-  position;
-  velocity;
   semimajorAxis_;
   meanMotion_;
   static fromState(origin, radialPosition, intrackPosition, nodeVelocity, nodeOffsetTime) {
@@ -13915,7 +13926,7 @@ var Tle = class _Tle {
     this.line1 = line1;
     this.line2 = line2;
     this.epoch = _Tle.parseEpoch_(line1.substring(18, 32));
-    this.satnum = parseInt(_Tle.convertA5to6Digit(line1.substring(2, 7)));
+    this.satnum = _Tle.satNum(this.line1);
     this.satrec_ = Sgp4.createSatrec(line1, line2, gravConst, opsMode);
   }
   toString() {
@@ -14004,7 +14015,7 @@ ${this.line2}`;
       epochYearFull = 1900 + epochYearShort;
     }
     const epochJday = epochDayOfYear + epochYearFull * 365;
-    const currentJday = getDayOfYear() + currentYearFull * 365;
+    const currentJday = getDayOfYear(nowInput) + currentYearFull * 365;
     const currentTime = (nowInput.getUTCHours() * 3600 + nowInput.getUTCMinutes() * 60 + nowInput.getUTCSeconds()) / 86400;
     const daysOld = currentJday + currentTime - epochJday;
     switch (outputUnits) {
@@ -14152,7 +14163,7 @@ ${this.line2}`;
     const BSTAR_PART_3 = _Tle.bstar_.start + 6;
     const BSTAR_PART_4 = _Tle.bstar_.stop - 1;
     const bstarSymbol = tleLine1.substring(_Tle.bstar_.start, BSTAR_PART_2);
-    let bstar1 = parseFloat(`0.${tleLine1.substring(BSTAR_PART_2, BSTAR_PART_3)}`);
+    const mantissaFraction = tleLine1.substring(BSTAR_PART_2, BSTAR_PART_3);
     const exponentSymbol = tleLine1.substring(BSTAR_PART_3, BSTAR_PART_4);
     let exponent = parseInt(tleLine1.substring(BSTAR_PART_4, _Tle.bstar_.stop));
     if (exponentSymbol === "-") {
@@ -14160,12 +14171,16 @@ ${this.line2}`;
     } else if (exponentSymbol !== "+" && exponentSymbol !== " " && exponentSymbol !== "0") {
       throw new ParseError(`Invalid BSTAR exponent symbol: ${exponentSymbol}`, "TLE");
     }
-    bstar1 *= 10 ** exponent;
+    let sign2 = 1;
+    let integerPart = "0";
     if (bstarSymbol === "-") {
-      bstar1 *= -1;
-    } else if (bstarSymbol === "+" || bstarSymbol === " " || bstarSymbol === "0") ; else {
+      sign2 = -1;
+    } else if (bstarSymbol === "+" || bstarSymbol === " " || bstarSymbol === "0") ; else if (/^\d$/u.test(bstarSymbol)) {
+      integerPart = bstarSymbol;
+    } else {
       throw new ParseError(`Invalid BSTAR symbol: ${bstarSymbol}`, "TLE");
     }
+    const bstar1 = sign2 * parseFloat(`${integerPart}.${mantissaFraction}`) * 10 ** exponent;
     return toPrecision(bstar1, 14);
   }
   /**
@@ -14220,22 +14235,25 @@ ${this.line2}`;
   }
   /**
    * Private value - used by United States Space Force to reference the orbit model used to generate the Tle. Will
-   * always be seen as zero externally (e.g. by "us", unless you are "them" - in which case, hello!).
+   * almost always be seen as zero externally (e.g. by "us", unless you are "them" - in which case, hello!).
+   *
+   * A value of 1 is tolerated in addition to 0: the field is informational and does not change how SGP4/SDP4
+   * propagate the element set, and a handful of archival TLEs carry a 1 (originally an SGP marker).
    *
-   * Starting in 2024, this may contain a 4 if the Tle was generated using the new SGP4-XP model. Until the source code
-   * is released, there is no way to support that format in JavaScript or TypeScript.
+   * A value of 4 indicates the SGP4-XP model. Until that source code is released there is no way to support that
+   * format in JavaScript or TypeScript, so it is rejected explicitly. Any other value is treated as malformed.
    * @example 0
    * @param tleLine1 The first line of the Tle to parse.
-   * @returns The ephemeris type.
+   * @returns The ephemeris type (0 or 1).
    */
   static ephemerisType(tleLine1) {
     const ephemerisType = parseInt(tleLine1.substring(_Tle.ephemerisType_.start, _Tle.ephemerisType_.stop));
-    if (ephemerisType !== 0 && ephemerisType !== 4) {
-      throw new ParseError(`Invalid ephemeris type: ${ephemerisType}`, "TLE");
-    }
     if (ephemerisType === 4) {
       throw new ParseError("SGP4-XP ephemeris type is not supported", "TLE");
     }
+    if (ephemerisType !== 0 && ephemerisType !== 1) {
+      throw new ParseError(`Invalid ephemeris type: ${ephemerisType}`, "TLE");
+    }
     return ephemerisType;
   }
   /**
@@ -14466,6 +14484,9 @@ ${this.line2}`;
    */
   static satNum(tleLine) {
     const satNumStr = tleLine.substring(_Tle.satNum_.start, _Tle.satNum_.stop);
+    if (satNumStr.trim().length === 0) {
+      return NaN;
+    }
     const sixDigitSatNum = _Tle.convertA5to6Digit(satNumStr);
     return parseInt(sixDigitSatNum);
   }
@@ -14556,12 +14577,13 @@ ${this.line2}`;
   static parse(tleLine1, tleLine2) {
     const line1 = _Tle.parseLine1(tleLine1);
     const line2 = _Tle.parseLine2(tleLine2);
-    if (line1.satNum !== line2.satNum) {
+    const isVimpel = line1.classification === "V";
+    if (!isVimpel && line1.satNum !== line2.satNum) {
       console.info("Line 1 satNum:", line1.satNum);
       console.info("Line 2 satNum:", line2.satNum);
       throw new ParseError("Satellite numbers do not match between TLE lines", "TLE");
     }
-    if (line1.satNumRaw !== line2.satNumRaw) {
+    if (!isVimpel && line1.satNumRaw !== line2.satNumRaw) {
       console.info("Line 1 satNumRaw:", line1.satNumRaw);
       console.info("Line 2 satNumRaw:", line2.satNumRaw);
       throw new ParseError("Raw satellite numbers do not match between TLE lines", "TLE");
@@ -14602,10 +14624,11 @@ ${this.line2}`;
   static parseAll(tleLine1, tleLine2) {
     const line1 = _Tle.parseLine1(tleLine1);
     const line2 = _Tle.parseLine2(tleLine2);
-    if (line1.satNum !== line2.satNum) {
+    const isVimpel = line1.classification === "V";
+    if (!isVimpel && line1.satNum !== line2.satNum) {
       throw new ParseError("Satellite numbers do not match between TLE lines", "TLE");
     }
-    if (line1.satNumRaw !== line2.satNumRaw) {
+    if (!isVimpel && line1.satNumRaw !== line2.satNumRaw) {
       throw new ParseError("Raw satellite numbers do not match between TLE lines", "TLE");
     }
     if (line1.lineNumber1 !== 1) {
@@ -14961,9 +14984,6 @@ var SolarRadiationPressure = class _SolarRadiationPressure extends Force {
     this.area = area;
     this.reflectCoeff = reflectCoeff;
   }
-  mass;
-  area;
-  reflectCoeff;
   // / Solar pressure _(N/m²)_;
   static _kRef = 456e-8 * KM_PER_AU ** 2;
   acceleration(state) {
@@ -15346,8 +15366,6 @@ var ThirdBodyGravity = class _ThirdBodyGravity {
     this.moon = moon;
     this.sun = sun;
   }
-  moon;
-  sun;
   static moonGravity_(state) {
     const rMoon = Moon.eci(state.epoch.toDateTime());
     const aNum = rMoon.subtract(state.position);
@@ -15450,13 +15468,6 @@ var RAE = class _RAE {
     this.azRateRad = azRateRad;
     this.elRateRad = elRateRad;
   }
-  epoch;
-  rng;
-  azRad;
-  elRad;
-  rngRate;
-  azRateRad;
-  elRateRad;
   // / Create a new [Razel] object, using degrees for the angular values.
   static fromDegrees(epoch, range, azimuth, elevation, rangeRate, azimuthRate, elevationRate) {
     const azimuthRateRad = azimuthRate ? azimuthRate * DEG2RAD : void 0;
@@ -15925,9 +15936,6 @@ var RungeKuttaAdaptive = class _RungeKuttaAdaptive extends Propagator {
     this._cacheState = this.initState_;
     this.tolerance_ = Math.max(_RungeKuttaAdaptive._minTolerance, Math.abs(tolerance_));
   }
-  initState_;
-  forceModel_;
-  tolerance_;
   // / Initial state vector.
   _cacheState;
   _checkpoints = [];
@@ -16223,11 +16231,6 @@ var RungeKutta4Propagator = class extends Propagator {
     this.checkpoints_ = checkpoints_;
     this.stepSize_ = Math.abs(stepSize_);
   }
-  initState_;
-  forceModel_;
-  stepSize_;
-  cacheState_;
-  checkpoints_;
   // / Set the integrator step size to the provided number of [seconds].
   setStepSize(seconds) {
     this.stepSize_ = Math.abs(seconds);
@@ -16556,7 +16559,6 @@ var Sgp4Propagator = class extends Propagator {
     this.tle_ = tle_;
     this.cacheState_ = tle_.state.toJ2000();
   }
-  tle_;
   cacheState_;
   checkpoints_ = [];
   /**
@@ -16849,15 +16851,30 @@ var Satellite = class _Satellite extends SpaceObject {
     return objectId;
   }
   /**
-   * Derives {@link sccNum5} and {@link sccNum6} from {@link sccNum}. When the
-   * canonical ID is an extended (7+ digit) value that exceeds the TLE alpha-5
-   * capacity (max 339 999), both forms are set to `null`.
+   * Normalizes {@link sccNum} to the display-canonical numeric form and derives
+   * {@link sccNum5} and {@link sccNum6}. The class invariant is that
+   * `Satellite.sccNum` is always numeric — never an alpha-5 string. Alpha-5
+   * inputs ("T0001") are converted to their 6-digit numeric equivalent
+   * ("270001"); the alpha-5 form is preserved on {@link sccNum5}.
+   *
+   * Extended (7+ digit) IDs that exceed the alpha-5 capacity (max 339 999)
+   * leave sccNum5/sccNum6 set to `null`.
+   *
+   * Invalid input (empty string, malformed token) is passed through unchanged
+   * so callers can store placeholder sccNums on notional / debris stubs.
    */
   assignAlpha5Forms_() {
+    try {
+      this.sccNum = Tle.convertA5to6Digit(this.sccNum);
+    } catch {
+    }
+    if (/^0+\d/u.test(this.sccNum)) {
+      this.sccNum = this.sccNum.replace(/^0+/u, "");
+    }
     const kind = Tle.classifySatNum(this.sccNum);
-    if (kind === "numeric5" || kind === "alpha5" || kind === "numeric6") {
+    if (kind === "numeric5" || kind === "numeric6") {
       this.sccNum5 = Tle.convert6DigitToA5(this.sccNum);
-      this.sccNum6 = Tle.convertA5to6Digit(this.sccNum5);
+      this.sccNum6 = this.sccNum;
     } else {
       this.sccNum5 = null;
       this.sccNum6 = null;
@@ -19299,13 +19316,6 @@ var RadecGeocentric = class _RadecGeocentric {
     this.declinationRate = declinationRate;
     this.rangeRate = rangeRate;
   }
-  epoch;
-  rightAscension;
-  declination;
-  range;
-  rightAscensionRate;
-  declinationRate;
-  rangeRate;
   /**
    * Creates a RadecGeocentric object from the given parameters in degrees.
    * @param epoch - The epoch in UTC.
@@ -19448,13 +19458,6 @@ var RadecTopocentric = class _RadecTopocentric {
     this.declinationRate = declinationRate;
     this.rangeRate = rangeRate;
   }
-  epoch;
-  rightAscension;
-  declination;
-  range;
-  rightAscensionRate;
-  declinationRate;
-  rangeRate;
   /**
    * Create a new RadecTopocentric object, using degrees for the angular values.
    * @param epoch UTC epoch.
@@ -20012,8 +20015,6 @@ var RkCheckpoint = class {
     this.cacheState = cacheState;
     this.stepSize = stepSize;
   }
-  cacheState;
-  stepSize;
 };
 
 // src/propagator/RkResult.ts
@@ -20024,9 +20025,6 @@ var RkResult = class {
     this.error = error;
     this.newStep = newStep;
   }
-  state;
-  error;
-  newStep;
 };
 
 // src/covariance/StateCovariance.ts
@@ -20047,8 +20045,6 @@ var StateCovariance = class _StateCovariance {
     this.matrix = matrix;
     this.frame = frame;
   }
-  matrix;
-  frame;
   // / Create a new [StateCovariance] object from 1-sigma values.
   static fromSigmas(sigmas, frame) {
     const n = sigmas.length;
@@ -20079,8 +20075,6 @@ var PropagatorPairs = class {
     this.posStep_ = posStep_;
     this.velStep_ = velStep_;
   }
-  posStep_;
-  velStep_;
   _high = Array(6).fill(null);
   _low = Array(6).fill(null);
   set(index, high, low) {
@@ -20114,9 +20108,6 @@ var BatchLeastSquaresResult = class {
     this.covariance = covariance2;
     this.rms = rms;
   }
-  state;
-  covariance;
-  rms;
 };
 
 // src/orbit-determination/BatchLeastSquaresOD.ts
@@ -20147,12 +20138,6 @@ var BatchLeastSquaresOD = class _BatchLeastSquaresOD {
     this.propagator_ = new RungeKutta89Propagator(this.apriori_, this.forceModel_);
     this.nominal_ = this.propagator_.propagate(this.start_);
   }
-  observations_;
-  apriori_;
-  forceModel_;
-  posStep_;
-  velStep_;
-  fastDerivatives_;
   /** Propagator pair cache, for generating observation Jacobians. */
   propPairs_;
   /**  Nominal state propagator. */
@@ -20253,7 +20238,6 @@ var GibbsIOD = class _GibbsIOD {
   constructor(mu = Earth.mu) {
     this.mu = mu;
   }
-  mu;
   /** Abort solve if position plane exceeds this value. */
   static coplanarThreshold_ = 5 * DEG2RAD;
   /**
@@ -20312,7 +20296,6 @@ var HerrickGibbsIOD = class {
   constructor(mu = Earth.mu) {
     this.mu = mu;
   }
-  mu;
   // / Attempt to create a state estimate from three inertial position vectors.
   solve(r1, t1, r2, t2, r3, t3) {
     const dt31 = t3.difference(t1);
@@ -20573,7 +20556,6 @@ var LambertIOD = class _LambertIOD {
   constructor(mu = Earth.mu) {
     this.mu = mu;
   }
-  mu;
   /**
    * Try to guess the short path argument given an [interceptor] and
    * [target] state.
@@ -21261,11 +21243,6 @@ var LevenbergMarquardtResult = class {
     this.iterations = iterations;
     this.converged = converged;
   }
-  state;
-  covariance;
-  rms;
-  iterations;
-  converged;
 };
 
 // src/orbit-determination/LevenbergMarquardtOD.ts
@@ -21296,12 +21273,6 @@ var LevenbergMarquardtOD = class _LevenbergMarquardtOD {
     const propagator = new RungeKutta89Propagator(this.apriori_, this.forceModel_);
     this.nominal_ = propagator.propagate(this.start_);
   }
-  observations_;
-  apriori_;
-  forceModel_;
-  posStep_;
-  velStep_;
-  fastDerivatives_;
   /** Propagator pair cache, for generating observation Jacobians. */
   propPairs_;
   /**  State estimate during solve. */
@@ -21460,8 +21431,6 @@ var SimplexEntry = class _SimplexEntry {
     this.x_ = new Vector(points);
     this.score = this.f_(points);
   }
-  f_;
-  points;
   /**
    * The cost function score at this vertex.
    *
@@ -21794,6 +21763,82 @@ var ModifiedGoodingIOD = class {
   }
 };
 
+// src/orbit-determination/Rv2Tle.ts
+var angleDiff_ = (a, b) => {
+  let diff = (a - b) % TAU;
+  if (diff > Math.PI) {
+    diff -= TAU;
+  } else if (diff < -Math.PI) {
+    diff += TAU;
+  }
+  return diff;
+};
+var wrapTau_ = (angle) => {
+  let wrapped = angle % TAU;
+  if (wrapped < 0) {
+    wrapped += TAU;
+  }
+  return wrapped;
+};
+var toElements_ = (epoch, position, velocity) => new J20002(
+  epoch,
+  new Vector3D(position.x, position.y, position.z),
+  new Vector3D(velocity.x, velocity.y, velocity.z)
+).toClassicalElements();
+var rv2tle = (epoch, position, velocity, options = {}) => {
+  const maxIterations = options.maxIterations ?? 15;
+  const toleranceKm = options.toleranceKm ?? 1e-4;
+  const epochUtc = EpochUTC.fromDateTime(epoch);
+  const target = toElements_(epochUtc, position, velocity);
+  let mean2 = new ClassicalElements({
+    epoch: epochUtc,
+    semimajorAxis: target.semimajorAxis,
+    eccentricity: target.eccentricity,
+    inclination: target.inclination,
+    rightAscension: target.rightAscension,
+    argPerigee: target.argPerigee,
+    trueAnomaly: target.trueAnomaly
+  });
+  let best = null;
+  for (let iteration = 1; iteration <= maxIterations; iteration++) {
+    let tle;
+    let sv;
+    try {
+      tle = Tle.fromClassicalElements(mean2);
+      const tsinceMin = (epoch.getTime() - tle.epoch.toDateTime().getTime()) / 6e4;
+      sv = Sgp4.propagate(Sgp4.createSatrec(tle.line1, tle.line2), tsinceMin);
+    } catch {
+      return best;
+    }
+    if (!sv.position || typeof sv.position === "boolean" || !sv.velocity || typeof sv.velocity === "boolean") {
+      return best;
+    }
+    const achievedPos = sv.position;
+    const achievedVel = sv.velocity;
+    const dx = achievedPos.x - position.x;
+    const dy = achievedPos.y - position.y;
+    const dz = achievedPos.z - position.z;
+    const positionErrorKm = Math.sqrt(dx * dx + dy * dy + dz * dz);
+    if (!best || positionErrorKm < best.positionErrorKm) {
+      best = { tle1: tle.line1, tle2: tle.line2, positionErrorKm, iterations: iteration };
+    }
+    if (positionErrorKm < toleranceKm) {
+      break;
+    }
+    const achieved = toElements_(epochUtc, achievedPos, achievedVel);
+    mean2 = new ClassicalElements({
+      epoch: epochUtc,
+      semimajorAxis: mean2.semimajorAxis + (target.semimajorAxis - achieved.semimajorAxis),
+      eccentricity: Math.min(Math.max(mean2.eccentricity + (target.eccentricity - achieved.eccentricity), 0), 0.999),
+      inclination: Math.min(Math.max(mean2.inclination + angleDiff_(target.inclination, achieved.inclination), 0), Math.PI),
+      rightAscension: wrapTau_(mean2.rightAscension + angleDiff_(target.rightAscension, achieved.rightAscension)),
+      argPerigee: wrapTau_(mean2.argPerigee + angleDiff_(target.argPerigee, achieved.argPerigee)),
+      trueAnomaly: wrapTau_(mean2.trueAnomaly + angleDiff_(target.trueAnomaly, achieved.trueAnomaly))
+    });
+  }
+  return best;
+};
+
 // src/covariance/CovarianceSample.ts
 var CovarianceSample = class {
   origin_;
@@ -21986,6 +22031,19 @@ var CovarianceSample = class {
   }
 };
 
+// src/covariance/displayRadii.ts
+function scaleAndClampRicSigmas(sigmas, confidence, caps) {
+  const { radial, inTrack, crossTrack } = sigmas;
+  if (![radial, inTrack, crossTrack].every((v) => Number.isFinite(v) && v >= 0)) {
+    return null;
+  }
+  return {
+    radial: Math.min(radial * confidence, caps.radial),
+    inTrack: Math.min(inTrack * confidence, caps.inTrack),
+    crossTrack: Math.min(crossTrack * confidence, caps.crossTrack)
+  };
+}
+
 // src/conjunction/ConjunctionAssessment.ts
 var ConjunctionAssessment = class {
   constructor(primary, secondary) {
@@ -21994,8 +22052,6 @@ var ConjunctionAssessment = class {
     this.primaryProp = this.createPropagator(primary, false);
     this.secondaryProp = this.createPropagator(secondary, false);
   }
-  primary;
-  secondary;
   primaryProp;
   secondaryProp;
   primaryCovSample;
@@ -23441,9 +23497,6 @@ var ObservationRadar = class _ObservationRadar extends Observation {
     this.observation = observation;
     this.noise_ = noise_;
   }
-  site_;
-  observation;
-  noise_;
   // / Default noise matrix.
   static defaultNoise = _ObservationRadar.noiseFromSigmas(0.32, 0.015 * DEG2RAD, 0.015 * DEG2RAD);
   get epoch() {
@@ -23933,9 +23986,6 @@ var ObservationOptical = class _ObservationOptical extends Observation {
     this.observation = observation;
     this.noise_ = noise_;
   }
-  site_;
-  observation;
-  noise_;
   // / Inertial site location.
   get site() {
     return this.site_;
@@ -27229,11 +27279,6 @@ var ChebyshevCoefficients = class _ChebyshevCoefficients {
     this.cyd_ = _ChebyshevCoefficients._derivative(a, b, cy_);
     this.czd_ = _ChebyshevCoefficients._derivative(a, b, cz_);
   }
-  a;
-  b;
-  cx_;
-  cy_;
-  cz_;
   cxd_;
   cyd_;
   czd_;
@@ -27755,11 +27800,6 @@ var TwoBurnOrbitTransfer = class _TwoBurnOrbitTransfer {
     this.vTransB = vTransB;
     this.tTrans = tTrans;
   }
-  vInit;
-  vFinal;
-  vTransA;
-  vTransB;
-  tTrans;
   /**
    * Calculates the parameters for a Hohmann transfer orbit between two circular orbits.
    * @param rInit The initial radius of the orbit. (km)
@@ -27812,7 +27852,6 @@ var ChebyshevCompressor = class _ChebyshevCompressor {
   constructor(interpolator_) {
     this.interpolator_ = interpolator_;
   }
-  interpolator_;
   /** Returns the cosine of π times x. */
   static cosPi_(x) {
     return Math.cos(Math.PI * x);
@@ -28053,8 +28092,6 @@ var SimpleLinearRegression = class _SimpleLinearRegression {
     this.ys = ys;
     this.update();
   }
-  xs;
-  ys;
   /** Line slope (m in y = mx + b). */
   slope_ = 0;
   /** Y-axis intercept (b in y = mx + b). */
@@ -29809,8 +29846,20 @@ var OrbitFinder = class _OrbitFinder {
     this.debug_(`Found perigee at altitude ${perigeeParams.altitude} km`);
     return perigeeParams;
   }
+  /**
+   * Returns the 5-char satnum to embed in the synthesized TLEs.
+   *
+   * Why not {@link Satellite.sccNum}? For extended (7+ digit) IDs sccNum is
+   * the full 9-digit canonical value, which would push the TLE past 69 chars
+   * and break Sgp4.createSatrec. The input tle1 always has a valid 5-char
+   * satnum at columns 2-6, and SGP4 itself doesn't care what value lives
+   * there as long as line 1 and line 2 agree — so use that.
+   */
+  tleSatNum_() {
+    return this.sat_.tle1.substring(2, 7);
+  }
   generateTle1_() {
-    return `1 ${this.sat_.sccNum}U ${this.sat_.tle1.substring(9, 17)} ${this.sat_.tle1.substring(18, 32)}${this.sat_.tle1.substring(32, 71)}`;
+    return `1 ${this.tleSatNum_()}U ${this.sat_.tle1.substring(9, 17)} ${this.sat_.tle1.substring(18, 32)}${this.sat_.tle1.substring(32, 71)}`;
   }
   generateTle2_(newParams) {
     const mergedParams = {
@@ -29825,7 +29874,7 @@ var OrbitFinder = class _OrbitFinder {
     const argPer = mergedParams.argOfPerigee.toFixed(4).padStart(8, "0");
     const meanA = mergedParams.meanAnomaly.toFixed(4).padStart(8, "0");
     const meanMo = this.sat_.tle2.substring(52, 63);
-    return `2 ${this.sat_.sccNum} ${inc} ${raan} ${ecc} ${argPer} ${meanA} ${meanMo}    10`;
+    return `2 ${this.tleSatNum_()} ${inc} ${raan} ${ecc} ${argPer} ${meanA} ${meanMo}    10`;
   }
 };
 
@@ -30185,6 +30234,33 @@ var RepeatGroundTrack = class _RepeatGroundTrack {
  * 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/>.
  */
+/**
+ * @author Theodore Kruczek.
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ *
+ * @license MIT License
+ *
+ * @Copyright (c) 2025 Theodore Kruczek
+ *
+ * Permission is hereby granted, free of charge, to any person obtaining a copy
+ * of this software and associated documentation files (the "Software"), to deal
+ * in the Software without restriction, including without limitation the rights
+ * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+ * copies of the Software, and to permit persons to whom the Software is
+ * furnished to do so, subject to the following conditions:
+ *
+ * The above copyright notice and this permission notice shall be included in all
+ * copies or substantial portions of the Software.
+ *
+ * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+ * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+ * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+ * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+ * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+ * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+ * SOFTWARE.
+ */
 /**
  * @author Theodore Kruczek
  * @description Orbital Object ToolKit (ootk) is a collection of tools for working
@@ -30286,4 +30362,4 @@ astronomy-engine/esm/astronomy.js:
    *)
 */
 
-export { AccessCalculator, AngularDiameterMethod, AngularDistanceMethod, Antenna, AtmosphericDrag, BaseObject, BatchLeastSquaresOD, BatchLeastSquaresResult, Beacon, BoxMuller, CatalogScreener, CatalogSource, CdmExporter, CdmParser, CelestialBody, CelestialBodyType, CenterBody, CenterBodyMu, ChebyshevCoefficients, ChebyshevCompressor, ChebyshevInterpolator, ClassicalElements, CommDeviceType, CommLink, CommunicationDevice, ConjunctionAssessment, ConjunctionEvent, ConstellationGenerator, ContactScheduler, ContactSelectionStrategy, CovarianceFrame, CovarianceRealism, CovarianceSample, CubicSpline, CubicSplineInterpolator, DEFAULT_INTERPOLATOR, DEFAULT_LAGRANGE_ORDER, DEG2RAD, DataHandler, DormandPrince54Propagator, DownhillSimplex, DynamicGroundObject, Earth, EarthGravity, EphemerisBody, EphemerisSatellite, Epoch, EpochGPS, EpochTAI, EpochTDB, EpochTT, EpochUTC, EpochWindow, EquinoctialElements, EulerAngles, FieldInterpolator, FieldOfView, Force, ForceModel, FormatTle, FovFrame, FovShape, Geodetic, GibbsIOD, GoldenSection, GoodingIOD, Gravity, GroundObject, GroundStation, HerrickGibbsIOD, Hill, History, HorizonsParser, ITRF, Interpolator, InterpolatorType, J20002 as J2000, Jupiter, KM_PER_AU, KeplerPropagator, LagrangeInterpolator, LambertIOD, LandObject, LaserRangingSensor, LaunchTrajectoryGenerator, LevenbergMarquardtOD, LevenbergMarquardtResult, MILLISECONDS_PER_DAY, MILLISECONDS_PER_SECOND, MILLISECONDS_TO_DAYS, MINUTES_PER_DAY, MS_PER_DAY, Marker, Mars, Matrix, MechanicalRadar, Mercury, ModifiedGoodingIOD, ModulationType, Moon, MoonBody, Neptune, OdmExporter, OemParser, OmmParser, OotkError, OpticalSensor, OrbitDeterminationError, OrbitFinder, OrbitRegime, PI, ParseError, PassType, PassiveRFSensor, PayloadStatus, PhasedArrayRadar, PlaneChangeBurn, PlanetBody, Pluto, PolynomialRegression, ProbabilityOfCollision, PropagationError, Propagator, PropagatorType, Quaternion, RAD2DEG, RADIUS_OF_EARTH, RAE, RIC, RadarSensor, RadecGeocentric, RadecTopocentric, Random, RandomGaussianSource, Receiver, RelativeState, RepeatGroundTrack, RkCheckpoint, RkResult, RungeKutta4Propagator, RungeKutta89Propagator, RungeKuttaAdaptive, SLR_WAVELENGTHS, SPEED_OF_LIGHT_KM_S, Satellite, Saturn, ScheduledContact, ScreeningFilter, SegmentedLagrangeInterpolator, Sensor, SensorType, Sgp4, Sgp4Error, Sgp4ErrorCode, Sgp4OpsMode, Sgp4Propagator, SimpleLinearRegression, SolarRadiationPressure, SolarSystem, SpaceObject, SpaceObjectType, Star, StateCovariance, StateInterpolator, StateVector2 as StateVector, Sun, SunBody, SunStatus, TAU, TEME, ThirdBodyGravity, Thrust, TimeStamped, Tle, Transmitter, Transponder, TwoBurnOrbitTransfer, Uranus, ValidationError, Vector, Vector3D, Venus, VerletBlendInterpolator, VisualizationHelpers, Waypoint, ZoomValue, acoth, acsch, angularDiameter, angularDistance, angularVelocityOfEarth, array2d, asec2rad, asech, azel2uv, bodyTypeLookup, boresightFrameFromAzElRoll, cKmPerMs, cKmPerSec, cMPerSec, calcGmst, calcIncFromAz, calcInertAz, calculateFspl, calculatePropagationDelay, clamp, concat, copySign, covariance, createCovarianceFromTle, createSampleCovarianceFromTle, createVec, csch, dbmToDbw, dbwToDbm, dbwToWatts, deg2rad, derivative, dopplerFactor, earthGravityParam, ecef2eci, ecef2enu, ecef2rae, ecefRad2rae, eci2ecef, eci2lla, eci2rae, enu2rf, estimateRcs, evalPoly, factorial, gamma, getDayOfYear, getDegLat, getDegLon, getRadLat, getRadLon, halfPi, isLeapYear, jacobian, jday, linearDistance, linearInterpolate, lla2ecef, lla2eci, lla2sez, llaRad2ecef, log10, mag2db, masec2rad, matchHalfPlane, mean, msec2sec, newtonM, newtonNu, normalizeAngle, observationDerivative, observationNoiseFromSigmas, parseCenterBody, rad2deg, radecToPosition, radecToVelocity, rae2ecef, rae2eci, rae2enu, rae2raeOffBoresight, rae2ruv, rae2sez, relativeVelocity, sec2day, sec2deg, sec2min, sech, secondsPerDay, secondsPerSiderealDay, secondsPerWeek, sez2rae, sign, spaceObjType2Str, std, temp4, toPrecision, ttasec2rad, uv2azel, wattsToDbw, wrapAngle, x2o3 };
+export { AccessCalculator, AngularDiameterMethod, AngularDistanceMethod, Antenna, AtmosphericDrag, BaseObject, BatchLeastSquaresOD, BatchLeastSquaresResult, Beacon, BoxMuller, CatalogScreener, CatalogSource, CdmExporter, CdmParser, CelestialBody, CelestialBodyType, CenterBody, CenterBodyMu, ChebyshevCoefficients, ChebyshevCompressor, ChebyshevInterpolator, ClassicalElements, CommDeviceType, CommLink, CommunicationDevice, ConjunctionAssessment, ConjunctionEvent, ConstellationGenerator, ContactScheduler, ContactSelectionStrategy, CovarianceFrame, CovarianceRealism, CovarianceSample, CubicSpline, CubicSplineInterpolator, DEFAULT_INTERPOLATOR, DEFAULT_LAGRANGE_ORDER, DEG2RAD, DataHandler, DormandPrince54Propagator, DownhillSimplex, DynamicGroundObject, Earth, EarthGravity, EphemerisBody, EphemerisSatellite, Epoch, EpochGPS, EpochTAI, EpochTDB, EpochTT, EpochUTC, EpochWindow, EquinoctialElements, EulerAngles, FieldInterpolator, FieldOfView, Force, ForceModel, FormatTle, FovFrame, FovShape, Geodetic, GibbsIOD, GoldenSection, GoodingIOD, Gravity, GroundObject, GroundStation, HerrickGibbsIOD, Hill, History, HorizonsParser, ITRF, Interpolator, InterpolatorType, J20002 as J2000, Jupiter, KM_PER_AU, KeplerPropagator, LagrangeInterpolator, LambertIOD, LandObject, LaserRangingSensor, LaunchTrajectoryGenerator, LevenbergMarquardtOD, LevenbergMarquardtResult, MILLISECONDS_PER_DAY, MILLISECONDS_PER_SECOND, MILLISECONDS_TO_DAYS, MINUTES_PER_DAY, MS_PER_DAY, Marker, Mars, Matrix, MechanicalRadar, Mercury, ModifiedGoodingIOD, ModulationType, Moon, MoonBody, Neptune, OdmExporter, OemParser, OmmParser, OotkError, OpticalSensor, OrbitDeterminationError, OrbitFinder, OrbitRegime, PI, ParseError, PassType, PassiveRFSensor, PayloadStatus, PhasedArrayRadar, PlaneChangeBurn, PlanetBody, Pluto, PolynomialRegression, ProbabilityOfCollision, PropagationError, Propagator, PropagatorType, Quaternion, RAD2DEG, RADIUS_OF_EARTH, RAE, RCS_VMAG_ESTIMATE_MAX, RCS_VMAG_ESTIMATE_MIN, RIC, RadarSensor, RadecGeocentric, RadecTopocentric, Random, RandomGaussianSource, Receiver, RelativeState, RepeatGroundTrack, RkCheckpoint, RkResult, RungeKutta4Propagator, RungeKutta89Propagator, RungeKuttaAdaptive, SLR_WAVELENGTHS, SPEED_OF_LIGHT_KM_S, Satellite, Saturn, ScheduledContact, ScreeningFilter, SegmentedLagrangeInterpolator, Sensor, SensorType, Sgp4, Sgp4Error, Sgp4ErrorCode, Sgp4OpsMode, Sgp4Propagator, SimpleLinearRegression, SolarRadiationPressure, SolarSystem, SpaceObject, SpaceObjectType, Star, StateCovariance, StateInterpolator, StateVector2 as StateVector, Sun, SunBody, SunStatus, TAU, TEME, ThirdBodyGravity, Thrust, TimeStamped, Tle, Transmitter, Transponder, TwoBurnOrbitTransfer, Uranus, ValidationError, Vector, Vector3D, Venus, VerletBlendInterpolator, VisualizationHelpers, Waypoint, ZoomValue, acoth, acsch, angularDiameter, angularDistance, angularVelocityOfEarth, array2d, asec2rad, asech, azel2uv, bodyTypeLookup, boresightFrameFromAzElRoll, cKmPerMs, cKmPerSec, cMPerSec, calcGmst, calcIncFromAz, calcInertAz, calculateFspl, calculatePropagationDelay, cappedScreeningCovarianceFromTle, clamp, concat, copySign, covariance, createCovarianceFromTle, createSampleCovarianceFromTle, createVec, csch, dbmToDbw, dbwToDbm, dbwToWatts, deg2rad, derivative, dopplerFactor, earthGravityParam, ecef2eci, ecef2enu, ecef2rae, ecefRad2rae, eci2ecef, eci2lla, eci2rae, enu2rf, estimateRcs, estimateVmagFromRcs, evalPoly, factorial, gamma, getDayOfYear, getDegLat, getDegLon, getRadLat, getRadLon, halfPi, isLeapYear, jacobian, jday, linearDistance, linearInterpolate, lla2ecef, lla2eci, lla2sez, llaRad2ecef, log10, mag2db, masec2rad, matchHalfPlane, mean, msec2sec, newtonM, newtonNu, normalizeAngle, observationDerivative, observationNoiseFromSigmas, parseCenterBody, rad2deg, radecToPosition, radecToVelocity, rae2ecef, rae2eci, rae2enu, rae2raeOffBoresight, rae2ruv, rae2sez, relativeVelocity, rv2tle, scaleAndClampRicSigmas, sec2day, sec2deg, sec2min, sech, secondsPerDay, secondsPerSiderealDay, secondsPerWeek, sez2rae, sign, spaceObjType2Str, std, temp4, toPrecision, ttasec2rad, uv2azel, wattsToDbw, wrapAngle, x2o3 };