ootk 6.0.6 → 6.0.7

package/dist/main.cjs

@@ -23,8 +23,6 @@ var ValidationError = class extends OotkError {
     this.value = value;
     this.name = "ValidationError";
   }
-  field;
-  value;
 };
 var ParseError = class extends OotkError {
   /**
@@ -39,8 +37,6 @@ var ParseError = class extends OotkError {
     this.line = line;
     this.name = "ParseError";
   }
-  format;
-  line;
 };
 var PropagationError = class extends OotkError {
   /**
@@ -53,7 +49,6 @@ var PropagationError = class extends OotkError {
     this.epoch = epoch;
     this.name = "PropagationError";
   }
-  epoch;
 };
 var OrbitDeterminationError = class extends OotkError {
   /**
@@ -66,7 +61,6 @@ var OrbitDeterminationError = class extends OotkError {
     this.algorithm = algorithm;
     this.name = "OrbitDeterminationError";
   }
-  algorithm;
 };
 
 // src/enums/AngularDiameterMethod.ts
@@ -284,7 +278,6 @@ var Epoch = class {
       throw new ValidationError("Epoch posix time must be a valid number", "posix", posix);
     }
   }
-  posix;
   toString() {
     return this.toDateTime().toISOString();
   }
@@ -308,7 +301,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,
@@ -1828,8 +1821,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;
   /**
@@ -2255,6 +2246,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) => {
@@ -2281,6 +2298,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) {
@@ -2288,9 +2316,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.
@@ -2488,7 +2513,6 @@ var Vector = class _Vector {
     this.elements = elements;
     this.length = elements.length;
   }
-  elements;
   /**
    * The length of the vector.
    */
@@ -3269,7 +3293,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;
@@ -8609,12 +8633,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;
@@ -8690,7 +8708,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++) {
@@ -8890,7 +8907,6 @@ var VerletBlendInterpolator = class _VerletBlendInterpolator extends StateInterp
     super();
     this.ephemeris = ephemeris;
   }
-  ephemeris;
   get sizeBytes() {
     return 64 * 7 * this.ephemeris.length / 8;
   }
@@ -9878,7 +9894,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);
@@ -10391,7 +10410,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));
@@ -13558,11 +13577,6 @@ var Thrust = class {
       crosstrack * 1e-3
     );
   }
-  center;
-  radial;
-  intrack;
-  crosstrack;
-  durationRate;
   deltaV;
   get magnitude() {
     return this.deltaV.magnitude() * 1e3;
@@ -13601,9 +13615,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) {
@@ -13917,7 +13928,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() {
@@ -14006,7 +14017,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) {
@@ -14154,7 +14165,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 === "-") {
@@ -14162,12 +14173,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);
   }
   /**
@@ -14222,22 +14237,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;
   }
   /**
@@ -14468,6 +14486,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);
   }
@@ -14558,12 +14579,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");
@@ -14604,10 +14626,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) {
@@ -14963,9 +14986,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) {
@@ -15348,8 +15368,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);
@@ -15452,13 +15470,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;
@@ -15927,9 +15938,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 = [];
@@ -16225,11 +16233,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);
@@ -16558,7 +16561,6 @@ var Sgp4Propagator = class extends Propagator {
     this.tle_ = tle_;
     this.cacheState_ = tle_.state.toJ2000();
   }
-  tle_;
   cacheState_;
   checkpoints_ = [];
   /**
@@ -16851,15 +16853,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;
@@ -19301,13 +19318,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.
@@ -19450,13 +19460,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.
@@ -20014,8 +20017,6 @@ var RkCheckpoint = class {
     this.cacheState = cacheState;
     this.stepSize = stepSize;
   }
-  cacheState;
-  stepSize;
 };
 
 // src/propagator/RkResult.ts
@@ -20026,9 +20027,6 @@ var RkResult = class {
     this.error = error;
     this.newStep = newStep;
   }
-  state;
-  error;
-  newStep;
 };
 
 // src/covariance/StateCovariance.ts
@@ -20049,8 +20047,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;
@@ -20081,8 +20077,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) {
@@ -20116,9 +20110,6 @@ var BatchLeastSquaresResult = class {
     this.covariance = covariance2;
     this.rms = rms;
   }
-  state;
-  covariance;
-  rms;
 };
 
 // src/orbit-determination/BatchLeastSquaresOD.ts
@@ -20149,12 +20140,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. */
@@ -20255,7 +20240,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;
   /**
@@ -20314,7 +20298,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);
@@ -20575,7 +20558,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.
@@ -21263,11 +21245,6 @@ var LevenbergMarquardtResult = class {
     this.iterations = iterations;
     this.converged = converged;
   }
-  state;
-  covariance;
-  rms;
-  iterations;
-  converged;
 };
 
 // src/orbit-determination/LevenbergMarquardtOD.ts
@@ -21298,12 +21275,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. */
@@ -21462,8 +21433,6 @@ var SimplexEntry = class _SimplexEntry {
     this.x_ = new Vector(points);
     this.score = this.f_(points);
   }
-  f_;
-  points;
   /**
    * The cost function score at this vertex.
    *
@@ -21796,6 +21765,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_;
@@ -21988,6 +22033,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) {
@@ -21996,8 +22054,6 @@ var ConjunctionAssessment = class {
     this.primaryProp = this.createPropagator(primary, false);
     this.secondaryProp = this.createPropagator(secondary, false);
   }
-  primary;
-  secondary;
   primaryProp;
   secondaryProp;
   primaryCovSample;
@@ -23443,9 +23499,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() {
@@ -23935,9 +23988,6 @@ var ObservationOptical = class _ObservationOptical extends Observation {
     this.observation = observation;
     this.noise_ = noise_;
   }
-  site_;
-  observation;
-  noise_;
   // / Inertial site location.
   get site() {
     return this.site_;
@@ -27231,11 +27281,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_;
@@ -27757,11 +27802,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)
@@ -27814,7 +27854,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);
@@ -28055,8 +28094,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). */
@@ -29811,8 +29848,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 = {
@@ -29827,7 +29876,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`;
   }
 };
 
@@ -30187,6 +30236,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
@@ -30416,6 +30492,8 @@ exports.Quaternion = Quaternion;
 exports.RAD2DEG = RAD2DEG;
 exports.RADIUS_OF_EARTH = RADIUS_OF_EARTH;
 exports.RAE = RAE;
+exports.RCS_VMAG_ESTIMATE_MAX = RCS_VMAG_ESTIMATE_MAX;
+exports.RCS_VMAG_ESTIMATE_MIN = RCS_VMAG_ESTIMATE_MIN;
 exports.RIC = RIC;
 exports.RadarSensor = RadarSensor;
 exports.RadecGeocentric = RadecGeocentric;
@@ -30493,6 +30571,7 @@ exports.calcIncFromAz = calcIncFromAz;
 exports.calcInertAz = calcInertAz;
 exports.calculateFspl = calculateFspl;
 exports.calculatePropagationDelay = calculatePropagationDelay;
+exports.cappedScreeningCovarianceFromTle = cappedScreeningCovarianceFromTle;
 exports.clamp = clamp;
 exports.concat = concat;
 exports.copySign = copySign;
@@ -30517,6 +30596,7 @@ exports.eci2lla = eci2lla;
 exports.eci2rae = eci2rae;
 exports.enu2rf = enu2rf;
 exports.estimateRcs = estimateRcs;
+exports.estimateVmagFromRcs = estimateVmagFromRcs;
 exports.evalPoly = evalPoly;
 exports.factorial = factorial;
 exports.gamma = gamma;
@@ -30557,6 +30637,8 @@ exports.rae2raeOffBoresight = rae2raeOffBoresight;
 exports.rae2ruv = rae2ruv;
 exports.rae2sez = rae2sez;
 exports.relativeVelocity = relativeVelocity;
+exports.rv2tle = rv2tle;
+exports.scaleAndClampRicSigmas = scaleAndClampRicSigmas;
 exports.sec2day = sec2day;
 exports.sec2deg = sec2deg;
 exports.sec2min = sec2min;