ootk 6.0.4 → 6.0.6

package/dist/main.cjs

@@ -23,6 +23,8 @@ var ValidationError = class extends OotkError {
     this.value = value;
     this.name = "ValidationError";
   }
+  field;
+  value;
 };
 var ParseError = class extends OotkError {
   /**
@@ -37,6 +39,8 @@ var ParseError = class extends OotkError {
     this.line = line;
     this.name = "ParseError";
   }
+  format;
+  line;
 };
 var PropagationError = class extends OotkError {
   /**
@@ -49,6 +53,7 @@ var PropagationError = class extends OotkError {
     this.epoch = epoch;
     this.name = "PropagationError";
   }
+  epoch;
 };
 var OrbitDeterminationError = class extends OotkError {
   /**
@@ -61,6 +66,7 @@ var OrbitDeterminationError = class extends OotkError {
     this.algorithm = algorithm;
     this.name = "OrbitDeterminationError";
   }
+  algorithm;
 };
 
 // src/enums/AngularDiameterMethod.ts
@@ -79,6 +85,7 @@ var AngularDistanceMethod = /* @__PURE__ */ ((AngularDistanceMethod2) => {
 
 // src/enums/CatalogSource.ts
 var CatalogSource = /* @__PURE__ */ ((CatalogSource2) => {
+  CatalogSource2["UNKNOWN"] = "unknown";
   CatalogSource2["USSF"] = "spacetrack";
   CatalogSource2["CELESTRAK"] = "celestrak";
   CatalogSource2["CELESTRAK_SUP"] = "celestrak-sup";
@@ -277,6 +284,7 @@ var Epoch = class {
       throw new ValidationError("Epoch posix time must be a valid number", "posix", posix);
     }
   }
+  posix;
   toString() {
     return this.toDateTime().toISOString();
   }
@@ -1820,6 +1828,8 @@ 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;
   /**
@@ -2278,6 +2288,9 @@ 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.
@@ -2475,6 +2488,7 @@ var Vector = class _Vector {
     this.elements = elements;
     this.length = elements.length;
   }
+  elements;
   /**
    * The length of the vector.
    */
@@ -8595,6 +8609,12 @@ 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;
@@ -8670,6 +8690,7 @@ 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++) {
@@ -8869,6 +8890,7 @@ var VerletBlendInterpolator = class _VerletBlendInterpolator extends StateInterp
     super();
     this.ephemeris = ephemeris;
   }
+  ephemeris;
   get sizeBytes() {
     return 64 * 7 * this.ephemeris.length / 8;
   }
@@ -9028,6 +9050,7 @@ var OemParser = class _OemParser {
     let currentBlock = null;
     let currentMetadata = {};
     let metadataComments = [];
+    let userDefined = {};
     let inMeta = false;
     let inData = false;
     let inCovariance = false;
@@ -9036,6 +9059,7 @@ var OemParser = class _OemParser {
       if (line === "META_START") {
         currentMetadata = {};
         metadataComments = [];
+        userDefined = {};
         inMeta = true;
         continue;
       }
@@ -9043,6 +9067,9 @@ var OemParser = class _OemParser {
         if (metadataComments.length > 0) {
           currentMetadata.COMMENT = metadataComments;
         }
+        if (Object.keys(userDefined).length > 0) {
+          currentMetadata.USER_DEFINED = userDefined;
+        }
         currentBlock = {
           metadata: currentMetadata,
           ephemeris: []
@@ -9081,7 +9108,9 @@ var OemParser = class _OemParser {
           const eqIndex = line.indexOf("=");
           const key = line.substring(0, eqIndex).trim();
           const value = line.substring(eqIndex + 1).trim();
-          if (key === "INTERPOLATION_DEGREE") {
+          if (key.startsWith("USER_DEFINED_")) {
+            userDefined[key.substring(13)] = value;
+          } else if (key === "INTERPOLATION_DEGREE") {
             currentMetadata[key] = parseInt(value, 10);
           } else {
             currentMetadata[key] = value;
@@ -13529,6 +13558,11 @@ var Thrust = class {
       crosstrack * 1e-3
     );
   }
+  center;
+  radial;
+  intrack;
+  crosstrack;
+  durationRate;
   deltaV;
   get magnitude() {
     return this.deltaV.magnitude() * 1e3;
@@ -13567,6 +13601,9 @@ 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) {
@@ -14522,15 +14559,21 @@ ${this.line2}`;
     const line1 = _Tle.parseLine1(tleLine1);
     const line2 = _Tle.parseLine2(tleLine2);
     if (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) {
+      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");
     }
     if (line1.lineNumber1 !== 1) {
+      console.info("Line 1 line number:", line1.lineNumber1);
       throw new ParseError("First TLE line number must be 1", "TLE");
     }
     if (line2.lineNumber2 !== 2) {
+      console.info("Line 2 line number:", line2.lineNumber2);
       throw new ParseError("Second TLE line number must be 2", "TLE");
     }
     return {
@@ -14576,9 +14619,44 @@ ${this.line2}`;
     return { ...line1, ...line2 };
   }
   /**
-   * Converts a 6 digit SCC number to a 5 digit SCC alpha 5 number
-   * @param sccNum The 6 digit SCC number
-   * @returns The 5 digit SCC alpha 5 number
+   * Classifies a satellite catalog number by representation. See {@link SatNumKind}.
+   * Pure function; no side effects.
+   * @param sccNum The catalog number string.
+   * @returns The {@link SatNumKind} describing the input.
+   */
+  static classifySatNum(sccNum) {
+    if (typeof sccNum !== "string") {
+      return "invalid";
+    }
+    const trimmed = sccNum.trim();
+    if (trimmed.length === 0) {
+      return "invalid";
+    }
+    const first = trimmed[0];
+    if (trimmed.length === 5 && first in _Tle.alpha5_) {
+      return /^\d{4}$/u.test(trimmed.slice(1)) ? "alpha5" : "invalid";
+    }
+    if (!/^\d+$/u.test(trimmed)) {
+      return "invalid";
+    }
+    if (trimmed.length <= 5) {
+      return "numeric5";
+    }
+    if (trimmed.length === 6) {
+      return parseInt(trimmed, 10) <= 339999 ? "numeric6" : "extended";
+    }
+    return "extended";
+  }
+  /**
+   * Converts a 6-digit numeric SCC number to its 5-character alpha-5 form.
+   *
+   * Inputs shorter than 6 chars or already in alpha-5 form pass through unchanged.
+   *
+   * @param sccNum The SCC number to convert.
+   * @returns The 5-character alpha-5 representation.
+   * @throws {ValidationError} If `sccNum` exceeds the alpha-5 range (numeric value > 339 999,
+   *   or length > 6). For such IDs the canonical value should be kept on
+   *   `Satellite.sccNum` and the TLE column populated with the last 5 digits.
    */
   static convert6DigitToA5(sccNum) {
     if (sccNum.length < 6) {
@@ -14587,9 +14665,27 @@ ${this.line2}`;
     if (typeof sccNum[0] !== "string") {
       throw new ParseError("Invalid SCC number format", "TLE");
     }
+    if (sccNum.length > 6) {
+      throw new ValidationError(
+        "SCC number exceeds TLE alpha-5 capacity (max 339999); use Satellite.sccNum for the canonical ID",
+        "sccNum",
+        sccNum
+      );
+    }
     if (RegExp(/[A-Z]/iu, "u").test(sccNum[0])) {
       return sccNum;
     }
+    const numericValue = parseInt(sccNum, 10);
+    if (isNaN(numericValue)) {
+      throw new ValidationError("SCC number must be numeric", "sccNum", sccNum);
+    }
+    if (numericValue > 339999) {
+      throw new ValidationError(
+        "SCC number exceeds TLE alpha-5 capacity (max 339999); use Satellite.sccNum for the canonical ID",
+        "sccNum",
+        sccNum
+      );
+    }
     const rest = sccNum.slice(2, 6);
     let first = parseInt(`${sccNum[0]}${sccNum[1]}`);
     const iPlus = first >= 18 ? 1 : 0;
@@ -14598,9 +14694,17 @@ ${this.line2}`;
     return `${String.fromCharCode(first + 55)}${rest}`;
   }
   /**
-   * Converts a 5-digit SCC number to a 6-digit SCC number.
-   * @param sccNum - The 5-digit SCC number to convert.
-   * @returns The converted 6-digit SCC number.
+   * Converts a 5-character alpha-5 SCC number to its 6-digit numeric form.
+   *
+   * Inputs shorter than 5 chars pass through unchanged. Numeric inputs without
+   * an alpha-5 leading letter pass through unchanged (5-digit numeric, in-range
+   * 6-digit numeric, and extended 7+ digit numeric IDs are all identity).
+   *
+   * @param sccNum The SCC number to convert.
+   * @returns The 6-digit numeric representation, or the input itself if it
+   *   is already a numeric form.
+   * @throws {ValidationError} If `sccNum` is malformed: 6-digit numeric whose
+   *   value exceeds 339 999, or contains stray letters in a non-leading position.
    */
   static convertA5to6Digit(sccNum) {
     if (sccNum.length < 5) {
@@ -14611,10 +14715,33 @@ ${this.line2}`;
       throw new ParseError("Invalid SCC number format", "TLE");
     }
     if (values[0] in _Tle.alpha5_) {
+      if (!/^\d+$/u.test(sccNum.slice(1))) {
+        throw new ValidationError(
+          "Alpha-5 SCC number must have 4 trailing digits",
+          "sccNum",
+          sccNum
+        );
+      }
       const firstLetter = values[0];
       values[0] = _Tle.alpha5_[firstLetter];
+      return values.join("");
+    }
+    if (!/^\s*\d+\s*$/u.test(sccNum)) {
+      throw new ValidationError(
+        "SCC number must be alpha-5 or all-numeric",
+        "sccNum",
+        sccNum
+      );
     }
-    return values.join("");
+    const trimmed = sccNum.trim();
+    if (trimmed.length === 6 && parseInt(trimmed, 10) > 339999) {
+      throw new ValidationError(
+        "6-digit SCC number exceeds TLE alpha-5 capacity (max 339999)",
+        "sccNum",
+        sccNum
+      );
+    }
+    return sccNum;
   }
 };
 
@@ -14836,6 +14963,9 @@ 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) {
@@ -15218,6 +15348,8 @@ 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);
@@ -15320,6 +15452,13 @@ 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;
@@ -15788,6 +15927,9 @@ 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 = [];
@@ -16083,6 +16225,11 @@ 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);
@@ -16411,6 +16558,7 @@ var Sgp4Propagator = class extends Propagator {
     this.tle_ = tle_;
     this.cacheState_ = tle_.state.toJ2000();
   }
+  tle_;
   cacheState_;
   checkpoints_ = [];
   /**
@@ -16497,11 +16645,15 @@ var Satellite = class _Satellite extends SpaceObject {
   period;
   rightAscension;
   satrec;
-  /** The satellite catalog number as listed in the TLE. */
+  /** The canonical satellite catalog number. May be a 5-digit numeric, alpha-5,
+   * 6-digit numeric, or an extended (7+ digit) ID such as CelesTrak supplemental
+   * 9-digit IDs. */
   sccNum;
-  /** The 5 digit alpha-numeric satellite catalog number. */
+  /** The 5-character alpha-5 representation, or `null` when {@link sccNum} is
+   * an extended ID that exceeds the alpha-5 capacity (max numeric value 339 999). */
   sccNum5;
-  /** The 6 digit numeric satellite catalog number. */
+  /** The 6-digit numeric representation, or `null` when {@link sccNum} is
+   * an extended ID that exceeds the alpha-5 capacity (max numeric value 339 999). */
   sccNum6;
   tle1;
   tle2;
@@ -16683,13 +16835,42 @@ var Satellite = class _Satellite extends SpaceObject {
     return new _Satellite({ omm, name: name ?? omm.OBJECT_NAME });
   }
   // ==================== TLE/OMM Parsing ====================
+  /**
+   * Converts an OMM international designator (e.g. `"2026-114A"`) to the
+   * TLE intl-des column format (cols 10-17, e.g. `"26114A"`). Inputs that
+   * don't match the OMM `YYYY-NNNL...` pattern pass through unchanged.
+   */
+  static ommObjectIdToTleIntlDes_(objectId) {
+    if (!objectId) {
+      return "";
+    }
+    const match = objectId.match(/^\d{4}-(?<rest>\d{3}[A-Z]{1,3})$/u);
+    if (match?.groups?.rest) {
+      return objectId.slice(2, 4) + match.groups.rest;
+    }
+    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`.
+   */
+  assignAlpha5Forms_() {
+    const kind = Tle.classifySatNum(this.sccNum);
+    if (kind === "numeric5" || kind === "alpha5" || kind === "numeric6") {
+      this.sccNum5 = Tle.convert6DigitToA5(this.sccNum);
+      this.sccNum6 = Tle.convertA5to6Digit(this.sccNum5);
+    } else {
+      this.sccNum5 = null;
+      this.sccNum6 = null;
+    }
+  }
   parseTleAndUpdateOrbit_(tle1, tle2, sccNum) {
     const tleData = Tle.parse(tle1, tle2);
     this.tle1 = tle1;
     this.tle2 = tle2;
     this.sccNum = sccNum ?? tleData.satNum.toString();
-    this.sccNum5 = Tle.convert6DigitToA5(this.sccNum);
-    this.sccNum6 = Tle.convertA5to6Digit(this.sccNum5);
+    this.assignAlpha5Forms_();
     this.intlDes = tleData.intlDes;
     this.epochYear = tleData.epochYear;
     this.epochDay = tleData.epochDay;
@@ -16712,8 +16893,7 @@ var Satellite = class _Satellite extends SpaceObject {
   parseOmmAndUpdateOrbit_(omm) {
     const noradStr = String(omm.NORAD_CAT_ID);
     this.sccNum = noradStr.padStart(5, "0");
-    this.sccNum5 = Tle.convert6DigitToA5(noradStr);
-    this.sccNum6 = Tle.convertA5to6Digit(this.sccNum5);
+    this.assignAlpha5Forms_();
     this.intlDes = omm.OBJECT_ID;
     const YYYY = omm.EPOCH.slice(0, 4);
     const MM = omm.EPOCH.slice(5, 7);
@@ -16773,8 +16953,10 @@ var Satellite = class _Satellite extends SpaceObject {
       ecen: this.eccentricity,
       epochyr: this.epochYear,
       epochday: this.epochDay,
-      intl: omm.OBJECT_ID,
-      scc: this.sccNum,
+      intl: _Satellite.ommObjectIdToTleIntlDes_(omm.OBJECT_ID),
+      // Extended (7+ digit) IDs don't fit TLE cols 3-7. Truncate to the last 5
+      // digits for the TLE string; the canonical ID stays on this.sccNum.
+      scc: Tle.classifySatNum(this.sccNum) === "extended" ? this.sccNum.slice(-5) : this.sccNum,
       bstar: this.bstar,
       meanMotionDot: this.meanMoDev1,
       meanMotionDdot: this.meanMoDev2,
@@ -19119,6 +19301,13 @@ 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.
@@ -19261,6 +19450,13 @@ 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.
@@ -19818,6 +20014,8 @@ var RkCheckpoint = class {
     this.cacheState = cacheState;
     this.stepSize = stepSize;
   }
+  cacheState;
+  stepSize;
 };
 
 // src/propagator/RkResult.ts
@@ -19828,6 +20026,9 @@ var RkResult = class {
     this.error = error;
     this.newStep = newStep;
   }
+  state;
+  error;
+  newStep;
 };
 
 // src/covariance/StateCovariance.ts
@@ -19848,6 +20049,8 @@ 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;
@@ -19878,6 +20081,8 @@ 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) {
@@ -19911,6 +20116,9 @@ var BatchLeastSquaresResult = class {
     this.covariance = covariance2;
     this.rms = rms;
   }
+  state;
+  covariance;
+  rms;
 };
 
 // src/orbit-determination/BatchLeastSquaresOD.ts
@@ -19941,6 +20149,12 @@ 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. */
@@ -20041,6 +20255,7 @@ var GibbsIOD = class _GibbsIOD {
   constructor(mu = Earth.mu) {
     this.mu = mu;
   }
+  mu;
   /** Abort solve if position plane exceeds this value. */
   static coplanarThreshold_ = 5 * DEG2RAD;
   /**
@@ -20099,6 +20314,7 @@ 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);
@@ -20359,6 +20575,7 @@ 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.
@@ -21046,6 +21263,11 @@ var LevenbergMarquardtResult = class {
     this.iterations = iterations;
     this.converged = converged;
   }
+  state;
+  covariance;
+  rms;
+  iterations;
+  converged;
 };
 
 // src/orbit-determination/LevenbergMarquardtOD.ts
@@ -21076,6 +21298,12 @@ 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. */
@@ -21234,6 +21462,8 @@ var SimplexEntry = class _SimplexEntry {
     this.x_ = new Vector(points);
     this.score = this.f_(points);
   }
+  f_;
+  points;
   /**
    * The cost function score at this vertex.
    *
@@ -21766,6 +21996,8 @@ var ConjunctionAssessment = class {
     this.primaryProp = this.createPropagator(primary, false);
     this.secondaryProp = this.createPropagator(secondary, false);
   }
+  primary;
+  secondary;
   primaryProp;
   secondaryProp;
   primaryCovSample;
@@ -23211,6 +23443,9 @@ 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() {
@@ -23700,6 +23935,9 @@ var ObservationOptical = class _ObservationOptical extends Observation {
     this.observation = observation;
     this.noise_ = noise_;
   }
+  site_;
+  observation;
+  noise_;
   // / Inertial site location.
   get site() {
     return this.site_;
@@ -24182,12 +24420,12 @@ var HorizonsParser = class {
     for (let i = 0; i < lines.length; i++) {
       const line = lines[i].trim();
       if (line.startsWith("Target body name:")) {
-        targetName = this.extractValue_(line, "Target body name:");
+        targetName = this.extractTargetName_(line, "Target body name:");
         metadata.targetName = targetName;
       } else if (line.startsWith("Center body name:")) {
         centerBody = this.extractValue_(line, "Center body name:");
         metadata.centerBody = centerBody;
-      } else if (line.startsWith("Center-site name:")) {
+      } else if (line.startsWith("Center-site name:") && !centerBody) {
         centerBody = this.extractValue_(line, "Center-site name:");
         metadata.centerBody = centerBody;
       } else if (line.startsWith("Output type:")) {
@@ -24242,7 +24480,7 @@ var HorizonsParser = class {
     for (let i = 0; i < lines.length; i++) {
       const line = lines[i].trim();
       if (line.startsWith("Target body name:")) {
-        targetName = this.extractValue_(line, "Target body name:");
+        targetName = this.extractTargetName_(line, "Target body name:");
       } else if (line.startsWith("Center-site name:")) {
         observerLocation = this.extractValue_(line, "Center-site name:");
       }
@@ -24273,7 +24511,23 @@ var HorizonsParser = class {
    */
   static extractValue_(line, key) {
     const value = line.substring(key.length).trim();
-    return value.split(/\s{2,}/)[0].trim();
+    return value.split(/\s{2,}/u)[0].trim();
+  }
+  /**
+   * Extracts the target body name, stripping parenthetical descriptors and IDs.
+   *
+   * Horizons format examples:
+   * - "1 Ceres                       {source: astDys}" → "1 Ceres"
+   * - "Artemis II (spacecraft) (-1024) {source: ...}"  → "Artemis II"
+   * - "Mars (499)"                                     → "Mars"
+   */
+  static extractTargetName_(line, key) {
+    const raw = this.extractValue_(line, key);
+    const parenIdx = raw.search(/\s\(/u);
+    if (parenIdx > 0) {
+      return raw.substring(0, parenIdx).trim();
+    }
+    return raw;
   }
   /**
    * Parses a vector data line from Horizons output.
@@ -24343,7 +24597,9 @@ var HorizonsParser = class {
       return null;
     }
     const posLine = allLines[nextIndex].trim();
-    const posMatch = posLine.match(/X\s*=\s*([-\d.E+]+)\s+Y\s*=\s*([-\d.E+]+)\s+Z\s*=\s*([-\d.E+]+)/i);
+    const posMatchLabeled = posLine.match(/X\s*=\s*([-\d.E+]+)\s+Y\s*=\s*([-\d.E+]+)\s+Z\s*=\s*([-\d.E+]+)/iu);
+    const posMatchBare = !posMatchLabeled ? posLine.match(/^\s*([-\d.E+]+)\s+([-\d.E+]+)\s+([-\d.E+]+)\s*$/iu) : null;
+    const posMatch = posMatchLabeled ?? posMatchBare;
     if (!posMatch) {
       return null;
     }
@@ -24355,7 +24611,9 @@ var HorizonsParser = class {
     let velocity;
     if (nextIndex < allLines.length) {
       const velLine = allLines[nextIndex].trim();
-      const velMatch = velLine.match(/VX\s*=\s*([-\d.E+]+)\s+VY\s*=\s*([-\d.E+]+)\s+VZ\s*=\s*([-\d.E+]+)/i);
+      const velMatchLabeled = velLine.match(/VX\s*=\s*([-\d.E+]+)\s+VY\s*=\s*([-\d.E+]+)\s+VZ\s*=\s*([-\d.E+]+)/iu);
+      const velMatchBare = !velMatchLabeled ? velLine.match(/^\s*([-\d.E+]+)\s+([-\d.E+]+)\s+([-\d.E+]+)\s*$/iu) : null;
+      const velMatch = velMatchLabeled ?? velMatchBare;
       if (velMatch) {
         velocity = new Vector3D(
           parseFloat(velMatch[1]),
@@ -24365,6 +24623,17 @@ var HorizonsParser = class {
         nextIndex++;
       }
     }
+    while (nextIndex < allLines.length) {
+      const extraLine = allLines[nextIndex].trim();
+      if (extraLine === "$$EOE" || extraLine.length === 0 || /^\d+\.\d+\s*=\s*A\.D\./u.test(extraLine)) {
+        break;
+      }
+      if (/^\s*[-\d.E+]+(\s+[-\d.E+]+)*\s*$/u.test(extraLine)) {
+        nextIndex++;
+      } else {
+        break;
+      }
+    }
     return {
       data: { epoch, position, velocity },
       nextIndex
@@ -25111,6 +25380,58 @@ var OdmExporter = class _OdmExporter {
     return lines.join("\n");
   }
   // ==========================================================================
+  // OEM from State Vectors — direct J2000 array export
+  // ==========================================================================
+  /**
+   * Export an array of J2000 state vectors as OEM KVN.
+   *
+   * Unlike {@link formatOem}, which propagates a TLE-based satellite, this
+   * method serializes pre-computed state vectors directly. Useful for
+   * converting parsed ephemeris data (e.g., NASA Horizons) to CCSDS OEM.
+   *
+   * @param stateVectors - Array of J2000 state vectors to export
+   * @param metadata - Object identification metadata
+   * @param options - Export configuration
+   * @returns OEM KVN format string
+   */
+  static formatOemFromStateVectors(stateVectors, metadata, options = {}) {
+    if (stateVectors.length === 0) {
+      throw new Error("Cannot export OEM with no state vectors");
+    }
+    const refFrame = options.refFrame ?? "EME2000";
+    const centerName = options.centerName ?? "EARTH";
+    const startTime = stateVectors[0].epoch.toDateTime();
+    const stopTime = stateVectors[stateVectors.length - 1].epoch.toDateTime();
+    const lines = [];
+    lines.push("CCSDS_OEM_VERS = 2.0");
+    _OdmExporter.appendComments_(lines, options.comments);
+    lines.push(`CREATION_DATE = ${_OdmExporter.formatDateTime_(/* @__PURE__ */ new Date())}`);
+    lines.push(`ORIGINATOR = ${options.originator ?? "KeepTrack"}`);
+    if (options.messageId) {
+      lines.push(`MESSAGE_ID = ${options.messageId}`);
+    }
+    lines.push("");
+    lines.push("META_START");
+    lines.push(`OBJECT_NAME = ${metadata.objectName}`);
+    lines.push(`OBJECT_ID = ${metadata.objectId}`);
+    lines.push(`CENTER_NAME = ${centerName}`);
+    lines.push(`REF_FRAME = ${refFrame}`);
+    lines.push("TIME_SYSTEM = UTC");
+    lines.push(`START_TIME = ${_OdmExporter.formatDateTime_(startTime)}`);
+    lines.push(`STOP_TIME = ${_OdmExporter.formatDateTime_(stopTime)}`);
+    lines.push(`INTERPOLATION = ${options.interpolation ?? "LAGRANGE"}`);
+    lines.push(`INTERPOLATION_DEGREE = ${options.interpolationDegree ?? 7}`);
+    lines.push("META_STOP");
+    lines.push("");
+    for (const sv of stateVectors) {
+      const epoch = _OdmExporter.formatDateTime_(sv.epoch.toDateTime());
+      lines.push(
+        `${epoch}  ${_OdmExporter.formatNumber_(sv.position.x)}  ${_OdmExporter.formatNumber_(sv.position.y)}  ${_OdmExporter.formatNumber_(sv.position.z)}  ${_OdmExporter.formatNumber_(sv.velocity.x)}  ${_OdmExporter.formatNumber_(sv.velocity.y)}  ${_OdmExporter.formatNumber_(sv.velocity.z)}`
+      );
+    }
+    return lines.join("\n");
+  }
+  // ==========================================================================
   // OMM — Orbit Mean-Elements Message
   // ==========================================================================
   /**
@@ -26910,6 +27231,11 @@ 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_;
@@ -27431,6 +27757,11 @@ 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)
@@ -27483,6 +27814,7 @@ 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);
@@ -27723,6 +28055,8 @@ 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). */