ootk 6.0.6 → 7.0.0

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!).
    *
-   * 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 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).
+   *
+   * 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;
@@ -19288,6 +19305,206 @@ var observationNoiseFromSigmas = (sigmas) => {
   return new Matrix(result);
 };
 
+// src/observation/Observation.ts
+var Observation = class {
+  /**
+   * Convert this observation's noise matrix into a covariance matrix.
+   * @returns A matrix representing the noise covariance.
+   */
+  noiseCovariance() {
+    return this.noise.reciprocal();
+  }
+  /**
+   * Generates a noise sample from the noise covariance matrix.
+   * @param sigma - The scaling factor for the noise covariance matrix.
+   * @returns A matrix representing the noise sample.
+   */
+  noiseSample_(sigma) {
+    return this.noiseCovariance().scale(sigma).cholesky();
+  }
+  /**
+   * Randomly sample this observation in vector form within the
+   * observation noise.
+   * @param random Random number generator.
+   * @param sigma Sigma value to scale the noise by.
+   * @returns Sampled observation.
+   */
+  sampleVector(random, sigma) {
+    const chol = this.noiseSample_(sigma);
+    const gauss = random.gaussVector(this.noise.columns);
+    const meas = this.toVector();
+    return meas.add(chol.multiplyVector(gauss));
+  }
+};
+
+// src/observation/ObservationOptical.ts
+var ObservationOptical = class _ObservationOptical extends Observation {
+  // / Create a new [ObservationOptical] object.
+  constructor(site_, observation, noise_ = _ObservationOptical.defaultNoise) {
+    super();
+    this.site_ = site_;
+    this.observation = observation;
+    this.noise_ = noise_;
+  }
+  // / Inertial site location.
+  get site() {
+    return this.site_;
+  }
+  // / Noise matrix.
+  get noise() {
+    return this.noise_;
+  }
+  /**
+   * Default noise matrix _(right-ascension, declination)_.
+   * Based on the Maui Optical Site noise model.
+   */
+  static defaultNoise = _ObservationOptical.noiseFromSigmas(
+    37e-4 * DEG2RAD,
+    3e-3 * DEG2RAD
+  );
+  get epoch() {
+    return this.observation.epoch;
+  }
+  toVector() {
+    return Vector.fromList([this.observation.rightAscension, this.observation.declination]);
+  }
+  clos(propagator) {
+    const position = propagator.propagate(this.epoch).position;
+    const offset = position.subtract(this.site.position);
+    const actual = this.observation.lineOfSight().normalize();
+    const expected = offset.normalize();
+    const slantRange = offset.magnitude();
+    const theta = actual.angle(expected);
+    if (isNaN(theta)) {
+      return 0;
+    }
+    return 2 * slantRange * Math.sin(theta * 0.5);
+  }
+  ricDiff(propagator) {
+    const r0 = this.site;
+    const r1 = propagator.propagate(this.epoch);
+    const r2 = this.observation.position(this.site, r1.position.distance(r0.position));
+    return RIC.fromJ2000(new J20002(this.epoch, r2, Vector3D.origin), r1).position;
+  }
+  sample(random, sigma = 1) {
+    const resultEl = this.sampleVector(random, sigma).elements;
+    return new _ObservationOptical(
+      this.site,
+      new RadecTopocentric(this.epoch, resultEl[0], resultEl[1]),
+      this.noise
+    );
+  }
+  jacobian(propPairs) {
+    const result = array2d(2, 6, 0);
+    for (let i = 0; i < 6; i++) {
+      const step = propPairs.step(i);
+      const [high, low] = propPairs.get(i);
+      const sl = low.propagate(this.epoch);
+      const sh = high.propagate(this.epoch);
+      const ol = RadecTopocentric.fromStateVector(sl, this.site);
+      const oh = RadecTopocentric.fromStateVector(sh, this.site);
+      result[0][i] = observationDerivative(oh.rightAscension, ol.rightAscension, step, true);
+      result[1][i] = observationDerivative(oh.declination, ol.declination, step, true);
+    }
+    return new Matrix(result);
+  }
+  residual(propagator) {
+    const result = array2d(2, 1, 0);
+    const state = propagator.propagate(this.epoch);
+    const radec = RadecTopocentric.fromStateVector(state, this.site);
+    result[0][0] = normalizeAngle(this.observation.rightAscension, radec.rightAscension);
+    result[1][0] = normalizeAngle(this.observation.declination, radec.declination);
+    return new Matrix(result);
+  }
+  /**
+   * Create a noise matrix from right ascension and declination standard deviantions.
+   * @param raSigma Right ascension standard deviation
+   * @param decSigma Declination standard deviation
+   * @returns Noise matrix
+   */
+  static noiseFromSigmas(raSigma, decSigma) {
+    return observationNoiseFromSigmas([raSigma, decSigma]);
+  }
+};
+
+// src/observation/ObservationRadar.ts
+var ObservationRadar = class _ObservationRadar extends Observation {
+  // / Create a new [ObservationRadar] object.
+  constructor(site_, observation, noise_ = _ObservationRadar.defaultNoise) {
+    super();
+    this.site_ = site_;
+    this.observation = observation;
+    this.noise_ = noise_;
+  }
+  // / Default noise matrix.
+  static defaultNoise = _ObservationRadar.noiseFromSigmas(0.32, 0.015 * DEG2RAD, 0.015 * DEG2RAD);
+  get epoch() {
+    return this.observation.epoch;
+  }
+  get site() {
+    return this.site_;
+  }
+  get noise() {
+    return this.noise_;
+  }
+  toVector() {
+    return Vector.fromList([this.observation.rng, this.observation.azRad, this.observation.elRad]);
+  }
+  clos(propagator) {
+    const ri = propagator.propagate(this.epoch).position.subtract(this.site.position);
+    return Math.abs(this.observation.rng - ri.magnitude());
+  }
+  ricDiff(propagator) {
+    const r0 = propagator.propagate(this.epoch);
+    const r1 = RAE.fromStateVector(r0, this.site);
+    const r2 = this.observation.position(this.site, r1.azRad, r1.elRad);
+    return RIC.fromJ2000(new J20002(this.epoch, r2, Vector3D.origin), r0).position;
+  }
+  sample(random, sigma = 1) {
+    const result = this.sampleVector(random, sigma).elements;
+    return new _ObservationRadar(
+      this.site,
+      new RAE(this.observation.epoch, result[0], result[1], result[2]),
+      this.noise
+    );
+  }
+  jacobian(propPairs) {
+    const result = array2d(3, 6, 0);
+    for (let i = 0; i < 6; i++) {
+      const step = propPairs.step(i);
+      const [high, low] = propPairs.get(i);
+      const sl = low.propagate(this.epoch);
+      const sh = high.propagate(this.epoch);
+      const ol = RAE.fromStateVector(sl, this.site);
+      const oh = RAE.fromStateVector(sh, this.site);
+      result[0][i] = observationDerivative(oh.rng, ol.rng, step);
+      result[1][i] = observationDerivative(oh.azRad, ol.azRad, step, true);
+      result[2][i] = observationDerivative(oh.elRad, ol.elRad, step, true);
+    }
+    return new Matrix(result);
+  }
+  residual(propagator) {
+    const result = array2d(3, 1, 0);
+    const state = propagator.propagate(this.epoch);
+    const razel = RAE.fromStateVector(state, this.site);
+    result[0][0] = this.observation.rng - razel.rng;
+    result[1][0] = normalizeAngle(this.observation.azRad, razel.azRad);
+    result[2][0] = normalizeAngle(this.observation.elRad, razel.elRad);
+    return new Matrix(result);
+  }
+  /**
+   * Create a noise matrix from the range, azimuth, and elevation standard
+   * deviations _(kilometers/radians)_.
+   * @param rngSigma - The range standard deviation _(kilometers)_.
+   * @param azSigma - The azimuth standard deviation _(radians)_.
+   * @param elSigma - The elevation standard deviation _(radians)_.
+   * @returns The noise matrix.
+   */
+  static noiseFromSigmas(rngSigma, azSigma, elSigma) {
+    return observationNoiseFromSigmas([rngSigma, azSigma, elSigma]);
+  }
+};
+
 // src/observation/RadecGeocentric.ts
 var RadecGeocentric = class _RadecGeocentric {
   constructor(epoch, rightAscension, declination, range, rightAscensionRate, declinationRate, rangeRate) {
@@ -19299,13 +19516,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 +19658,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 +20215,6 @@ var RkCheckpoint = class {
     this.cacheState = cacheState;
     this.stepSize = stepSize;
   }
-  cacheState;
-  stepSize;
 };
 
 // src/propagator/RkResult.ts
@@ -20024,9 +20225,6 @@ var RkResult = class {
     this.error = error;
     this.newStep = newStep;
   }
-  state;
-  error;
-  newStep;
 };
 
 // src/covariance/StateCovariance.ts
@@ -20047,8 +20245,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 +20275,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 +20308,6 @@ var BatchLeastSquaresResult = class {
     this.covariance = covariance2;
     this.rms = rms;
   }
-  state;
-  covariance;
-  rms;
 };
 
 // src/orbit-determination/BatchLeastSquaresOD.ts
@@ -20147,12 +20338,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 +20438,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 +20496,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 +20756,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 +21443,6 @@ var LevenbergMarquardtResult = class {
     this.iterations = iterations;
     this.converged = converged;
   }
-  state;
-  covariance;
-  rms;
-  iterations;
-  converged;
 };
 
 // src/orbit-determination/LevenbergMarquardtOD.ts
@@ -21296,12 +21473,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 +21631,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 +21963,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 +22231,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 +22252,6 @@ var ConjunctionAssessment = class {
     this.primaryProp = this.createPropagator(primary, false);
     this.secondaryProp = this.createPropagator(secondary, false);
   }
-  primary;
-  secondary;
   primaryProp;
   secondaryProp;
   primaryCovSample;
@@ -23400,119 +23656,6 @@ var Sensor = class _Sensor {
   }
 };
 
-// src/observation/Observation.ts
-var Observation = class {
-  /**
-   * Convert this observation's noise matrix into a covariance matrix.
-   * @returns A matrix representing the noise covariance.
-   */
-  noiseCovariance() {
-    return this.noise.reciprocal();
-  }
-  /**
-   * Generates a noise sample from the noise covariance matrix.
-   * @param sigma - The scaling factor for the noise covariance matrix.
-   * @returns A matrix representing the noise sample.
-   */
-  noiseSample_(sigma) {
-    return this.noiseCovariance().scale(sigma).cholesky();
-  }
-  /**
-   * Randomly sample this observation in vector form within the
-   * observation noise.
-   * @param random Random number generator.
-   * @param sigma Sigma value to scale the noise by.
-   * @returns Sampled observation.
-   */
-  sampleVector(random, sigma) {
-    const chol = this.noiseSample_(sigma);
-    const gauss = random.gaussVector(this.noise.columns);
-    const meas = this.toVector();
-    return meas.add(chol.multiplyVector(gauss));
-  }
-};
-
-// src/observation/ObservationRadar.ts
-var ObservationRadar = class _ObservationRadar extends Observation {
-  // / Create a new [ObservationRadar] object.
-  constructor(site_, observation, noise_ = _ObservationRadar.defaultNoise) {
-    super();
-    this.site_ = site_;
-    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() {
-    return this.observation.epoch;
-  }
-  get site() {
-    return this.site_;
-  }
-  get noise() {
-    return this.noise_;
-  }
-  toVector() {
-    return Vector.fromList([this.observation.rng, this.observation.azRad, this.observation.elRad]);
-  }
-  clos(propagator) {
-    const ri = propagator.propagate(this.epoch).position.subtract(this.site.position);
-    return Math.abs(this.observation.rng - ri.magnitude());
-  }
-  ricDiff(propagator) {
-    const r0 = propagator.propagate(this.epoch);
-    const r1 = RAE.fromStateVector(r0, this.site);
-    const r2 = this.observation.position(this.site, r1.azRad, r1.elRad);
-    return RIC.fromJ2000(new J20002(this.epoch, r2, Vector3D.origin), r0).position;
-  }
-  sample(random, sigma = 1) {
-    const result = this.sampleVector(random, sigma).elements;
-    return new _ObservationRadar(
-      this.site,
-      new RAE(this.observation.epoch, result[0], result[1], result[2]),
-      this.noise
-    );
-  }
-  jacobian(propPairs) {
-    const result = array2d(3, 6, 0);
-    for (let i = 0; i < 6; i++) {
-      const step = propPairs.step(i);
-      const [high, low] = propPairs.get(i);
-      const sl = low.propagate(this.epoch);
-      const sh = high.propagate(this.epoch);
-      const ol = RAE.fromStateVector(sl, this.site);
-      const oh = RAE.fromStateVector(sh, this.site);
-      result[0][i] = observationDerivative(oh.rng, ol.rng, step);
-      result[1][i] = observationDerivative(oh.azRad, ol.azRad, step, true);
-      result[2][i] = observationDerivative(oh.elRad, ol.elRad, step, true);
-    }
-    return new Matrix(result);
-  }
-  residual(propagator) {
-    const result = array2d(3, 1, 0);
-    const state = propagator.propagate(this.epoch);
-    const razel = RAE.fromStateVector(state, this.site);
-    result[0][0] = this.observation.rng - razel.rng;
-    result[1][0] = normalizeAngle(this.observation.azRad, razel.azRad);
-    result[2][0] = normalizeAngle(this.observation.elRad, razel.elRad);
-    return new Matrix(result);
-  }
-  /**
-   * Create a noise matrix from the range, azimuth, and elevation standard
-   * deviations _(kilometers/radians)_.
-   * @param rngSigma - The range standard deviation _(kilometers)_.
-   * @param azSigma - The azimuth standard deviation _(radians)_.
-   * @param elSigma - The elevation standard deviation _(radians)_.
-   * @returns The noise matrix.
-   */
-  static noiseFromSigmas(rngSigma, azSigma, elSigma) {
-    return observationNoiseFromSigmas([rngSigma, azSigma, elSigma]);
-  }
-};
-
 // src/sensor/RadarSensor.ts
 var RadarSensor = class extends Sensor {
   /** Radar beamwidth in degrees */
@@ -23924,99 +24067,6 @@ var MechanicalRadar = class _MechanicalRadar extends RadarSensor {
   }
 };
 
-// src/observation/ObservationOptical.ts
-var ObservationOptical = class _ObservationOptical extends Observation {
-  // / Create a new [ObservationOptical] object.
-  constructor(site_, observation, noise_ = _ObservationOptical.defaultNoise) {
-    super();
-    this.site_ = site_;
-    this.observation = observation;
-    this.noise_ = noise_;
-  }
-  site_;
-  observation;
-  noise_;
-  // / Inertial site location.
-  get site() {
-    return this.site_;
-  }
-  // / Noise matrix.
-  get noise() {
-    return this.noise_;
-  }
-  /**
-   * Default noise matrix _(right-ascension, declination)_.
-   * Based on the Maui Optical Site noise model.
-   */
-  static defaultNoise = _ObservationOptical.noiseFromSigmas(
-    37e-4 * DEG2RAD,
-    3e-3 * DEG2RAD
-  );
-  get epoch() {
-    return this.observation.epoch;
-  }
-  toVector() {
-    return Vector.fromList([this.observation.rightAscension, this.observation.declination]);
-  }
-  clos(propagator) {
-    const position = propagator.propagate(this.epoch).position;
-    const offset = position.subtract(this.site.position);
-    const actual = this.observation.lineOfSight().normalize();
-    const expected = offset.normalize();
-    const slantRange = offset.magnitude();
-    const theta = actual.angle(expected);
-    if (isNaN(theta)) {
-      return 0;
-    }
-    return 2 * slantRange * Math.sin(theta * 0.5);
-  }
-  ricDiff(propagator) {
-    const r0 = this.site;
-    const r1 = propagator.propagate(this.epoch);
-    const r2 = this.observation.position(this.site, r1.position.distance(r0.position));
-    return RIC.fromJ2000(new J20002(this.epoch, r2, Vector3D.origin), r1).position;
-  }
-  sample(random, sigma = 1) {
-    const resultEl = this.sampleVector(random, sigma).elements;
-    return new _ObservationOptical(
-      this.site,
-      new RadecTopocentric(this.epoch, resultEl[0], resultEl[1]),
-      this.noise
-    );
-  }
-  jacobian(propPairs) {
-    const result = array2d(2, 6, 0);
-    for (let i = 0; i < 6; i++) {
-      const step = propPairs.step(i);
-      const [high, low] = propPairs.get(i);
-      const sl = low.propagate(this.epoch);
-      const sh = high.propagate(this.epoch);
-      const ol = RadecTopocentric.fromStateVector(sl, this.site);
-      const oh = RadecTopocentric.fromStateVector(sh, this.site);
-      result[0][i] = observationDerivative(oh.rightAscension, ol.rightAscension, step, true);
-      result[1][i] = observationDerivative(oh.declination, ol.declination, step, true);
-    }
-    return new Matrix(result);
-  }
-  residual(propagator) {
-    const result = array2d(2, 1, 0);
-    const state = propagator.propagate(this.epoch);
-    const radec = RadecTopocentric.fromStateVector(state, this.site);
-    result[0][0] = normalizeAngle(this.observation.rightAscension, radec.rightAscension);
-    result[1][0] = normalizeAngle(this.observation.declination, radec.declination);
-    return new Matrix(result);
-  }
-  /**
-   * Create a noise matrix from right ascension and declination standard deviantions.
-   * @param raSigma Right ascension standard deviation
-   * @param decSigma Declination standard deviation
-   * @returns Noise matrix
-   */
-  static noiseFromSigmas(raSigma, decSigma) {
-    return observationNoiseFromSigmas([raSigma, decSigma]);
-  }
-};
-
 // src/sensor/OpticalSensor.ts
 var OpticalSensor = class _OpticalSensor extends Sensor {
   /** Telescope aperture in meters */
@@ -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). */
@@ -29184,6 +29221,24 @@ var LaunchTrajectoryGenerator = class _LaunchTrajectoryGenerator {
     }
     return azimuth;
   }
+  /**
+   * Compute orbital inclination from a launch azimuth and launch site latitude.
+   *
+   * Inverse of {@link computeLaunchAzimuth}: from the spherical-triangle relation
+   * `cos(i) = sin(azimuth) * cos(latitude)`. The result is the magnitude of the
+   * inclination in degrees (always in [0, 180]); the launch direction (N/S) is
+   * implied by the azimuth quadrant, not by this value.
+   *
+   * @param azimuthDeg Launch azimuth in degrees (from north, clockwise).
+   * @param latDeg Launch site latitude in degrees.
+   * @returns Orbital inclination in degrees, in [0, 180].
+   */
+  static computeInclinationFromAzimuth(azimuthDeg, latDeg) {
+    const az = azimuthDeg * DEG2RAD;
+    const lat = latDeg * DEG2RAD;
+    const cosInc = Math.max(-1, Math.min(1, Math.sin(az) * Math.cos(lat)));
+    return Math.acos(cosInc) * RAD2DEG;
+  }
   /**
    * Estimate ascent duration based on target insertion altitude.
    * Loosely calibrated: ~480s for 200km, ~540s for 400km, ~600s for 800km.
@@ -29809,8 +29864,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 +29892,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 +30252,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 +30380,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, Observation, ObservationOptical, ObservationRadar, 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 };