ootk 6.0.1 → 6.0.3

package/dist/main.cjs

@@ -9849,7 +9849,7 @@ var FormatTle = class _FormatTle {
     const { inc, meanmo, rasc, argPe, meana, ecen, epochyr, epochday, intl } = tleParams;
     const scc = Tle.convert6DigitToA5(tleParams.scc);
     const epochYrStr = String(epochyr).padStart(2, "0");
-    const epochdayStr = parseFloat(String(epochday)).toFixed(8).padStart(12, "0");
+    const epochdayStr = Number(epochday).toFixed(8).padStart(12, "0");
     const incStr = _FormatTle.inclination(inc);
     const meanmoStr = _FormatTle.meanMotion(meanmo);
     const rascStr = _FormatTle.rightAscension(rasc);
@@ -9857,13 +9857,23 @@ var FormatTle = class _FormatTle {
     const meanaStr = _FormatTle.meanAnomaly(meana);
     const ecenStr = _FormatTle.eccentricity(ecen);
     const intlStr = intl.padEnd(8, " ");
-    let TLE1Ending = tleParams.sat ? tleParams.sat.tle1.substring(32, 71) : " +.00000000 +00000+0 +00000-0 0  9990";
-    TLE1Ending = TLE1Ending[1] === " " ? _FormatTle.setCharAt(TLE1Ending, 1, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[12] === " " ? _FormatTle.setCharAt(TLE1Ending, 12, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[21] === " " ? _FormatTle.setCharAt(TLE1Ending, 21, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[32] === " " ? _FormatTle.setCharAt(TLE1Ending, 32, "0") : TLE1Ending;
-    const tle1 = `1 ${scc}U ${intlStr} ${epochYrStr}${epochdayStr}${TLE1Ending}`;
-    const tle2 = `2 ${scc} ${incStr} ${rascStr} ${ecenStr} ${argPeStr} ${meanaStr} ${meanmoStr} 00010`;
+    const classification = tleParams.classification ?? "U";
+    const elementSetNo = (tleParams.elementSetNo ?? 999).toString().padStart(4, " ");
+    const ephemerisType = (tleParams.ephemerisType ?? 0).toString();
+    const revAtEpoch = (tleParams.revAtEpoch ?? 0).toString().padStart(5, " ");
+    let TLE1Ending;
+    if (tleParams.sat) {
+      TLE1Ending = tleParams.sat.tle1.substring(32, 71);
+    } else {
+      const mmDot = _FormatTle.formatMeanMotionDot(tleParams.meanMotionDot ?? 0);
+      const mmDdot = _FormatTle.formatTleExponential(tleParams.meanMotionDdot ?? 0);
+      const bstarStr = _FormatTle.formatTleExponential(tleParams.bstar ?? 0);
+      TLE1Ending = ` ${mmDot} ${mmDdot} ${bstarStr} ${ephemerisType} ${elementSetNo}0`;
+    }
+    const tle1Pre = `1 ${scc}${classification} ${intlStr} ${epochYrStr}${epochdayStr}${TLE1Ending}`;
+    const tle1 = _FormatTle.setCharAt(tle1Pre, 68, _FormatTle.tleChecksum(tle1Pre).toString());
+    const tle2Pre = `2 ${scc} ${incStr} ${rascStr} ${ecenStr} ${argPeStr} ${meanaStr} ${meanmoStr}${revAtEpoch}0`;
+    const tle2 = _FormatTle.setCharAt(tle2Pre, 68, _FormatTle.tleChecksum(tle2Pre).toString());
     return { tle1, tle2 };
   }
   /**
@@ -9877,7 +9887,7 @@ var FormatTle = class _FormatTle {
       argPe = argPe.toString();
     }
     const argPeNum = parseFloat(argPe).toFixed(4);
-    const argPe0 = argPeNum.padStart(8, "0");
+    const argPe0 = argPeNum.padStart(8, " ");
     if (argPe0.length !== 8) {
       throw new ValidationError("Argument of perigee must be 8 characters", "argPe", argPe0);
     }
@@ -9891,7 +9901,10 @@ var FormatTle = class _FormatTle {
    */
   static eccentricity(ecen) {
     if (typeof ecen === "number") {
-      ecen = ecen.toFixed(7);
+      const ecenStr = ecen.toString();
+      const dotIdx = ecenStr.indexOf(".");
+      const afterDot = dotIdx >= 0 ? ecenStr.substring(dotIdx + 1) : "";
+      ecen = `0.${afterDot.substring(0, 7).padEnd(7, "0")}`;
     }
     let ecen0 = ecen.padEnd(9, "0");
     if (ecen0[1] === ".") {
@@ -9915,7 +9928,7 @@ var FormatTle = class _FormatTle {
       inc = inc.toString();
     }
     const incNum = parseFloat(inc).toFixed(4);
-    const inc0 = incNum.padStart(8, "0");
+    const inc0 = incNum.padStart(8, " ");
     if (inc0.length !== 8) {
       throw new ValidationError("Inclination must be 8 characters", "inclination", inc0);
     }
@@ -9932,7 +9945,7 @@ var FormatTle = class _FormatTle {
       meana = meana.toString();
     }
     const meanaNum = parseFloat(meana).toFixed(4);
-    const meana0 = meanaNum.padStart(8, "0");
+    const meana0 = meanaNum.padStart(8, " ");
     if (meana0.length !== 8) {
       throw new ValidationError("Mean anomaly must be 8 characters", "meanAnomaly", meana0);
     }
@@ -9973,7 +9986,7 @@ var FormatTle = class _FormatTle {
       rasc = rasc.toString();
     }
     const rascNum = parseFloat(rasc).toFixed(4);
-    const rasc0 = rascNum.padStart(8, "0");
+    const rasc0 = rascNum.padStart(8, " ");
     if (rasc0.length !== 8) {
       throw new ValidationError("Right ascension must be 8 characters", "rightAscension", rasc0);
     }
@@ -9992,6 +10005,61 @@ var FormatTle = class _FormatTle {
     }
     return `${str.substring(0, index)}${chr}${str.substring(index + 1)}`;
   }
+  /**
+   * Format mean motion dot (first derivative / 2) for TLE line 1.
+   * Format: sign + ".NNNNNNNN" = 10 chars total.
+   * @param value - Mean motion dot value (rev/day^2)
+   * @returns Formatted 10-character string
+   */
+  static formatMeanMotionDot(value) {
+    const sign2 = value >= 0 ? " " : "-";
+    return sign2 + Math.abs(value).toFixed(8).substring(1);
+  }
+  /**
+   * Format a value in TLE exponential notation for BSTAR or mean motion ddot.
+   * Format: "sNNNNN±N" (8 chars) where mantissa has implied leading decimal point.
+   * Example: 0.00017507 → " 17507-3" (i.e., .17507 × 10^-3)
+   * @param value - The value to format
+   * @returns Formatted 8-character string
+   */
+  static formatTleExponential(value) {
+    if (value === 0) {
+      return " 00000+0";
+    }
+    const sign2 = value >= 0 ? " " : "-";
+    const absVal = Math.abs(value);
+    let exponent = Math.floor(Math.log10(absVal)) + 1;
+    let mantissa = absVal / Math.pow(10, exponent);
+    if (mantissa >= 1) {
+      mantissa /= 10;
+      exponent++;
+    }
+    if (mantissa < 0.1 && mantissa > 0) {
+      mantissa *= 10;
+      exponent--;
+    }
+    const mantissaStr = Math.round(mantissa * 1e5).toString().padStart(5, "0");
+    const expSign = exponent >= 0 ? "+" : "-";
+    const expStr = Math.abs(exponent).toString();
+    return sign2 + mantissaStr + expSign + expStr;
+  }
+  /**
+   * Compute TLE line checksum (modulo 10 sum of digits, '-' counts as 1).
+   * @param line - TLE line (first 68 characters are summed)
+   * @returns Checksum digit (0-9)
+   */
+  static tleChecksum(line) {
+    let sum = 0;
+    for (let i = 0; i < 68 && i < line.length; i++) {
+      const c = line[i];
+      if (c >= "0" && c <= "9") {
+        sum += parseInt(c);
+      } else if (c === "-") {
+        sum += 1;
+      }
+    }
+    return sum % 10;
+  }
 };
 
 // src/coordinate/RelativeState.ts
@@ -13187,13 +13255,28 @@ var ClassicalElements = class _ClassicalElements {
     const h = pos.cross(vel);
     const i = Math.acos(clamp(h.z / h.magnitude(), -1, 1));
     const n = Vector3D.zAxis.cross(h);
-    let o = Math.acos(clamp(n.x / n.magnitude(), -1, 1));
-    if (n.y < 0) {
-      o = TAU - o;
-    }
-    let w = n.angle(eVec);
-    if (eVec.z < 0) {
-      w = TAU - w;
+    const nMag = n.magnitude();
+    let o;
+    let w;
+    if (nMag < 1e-12) {
+      o = 0;
+      if (e > 1e-12) {
+        w = Math.atan2(eVec.y, eVec.x);
+        if (w < 0) {
+          w = w + TAU;
+        }
+      } else {
+        w = 0;
+      }
+    } else {
+      o = Math.acos(clamp(n.x / nMag, -1, 1));
+      if (n.y < 0) {
+        o = TAU - o;
+      }
+      w = n.angle(eVec);
+      if (eVec.z < 0) {
+        w = TAU - w;
+      }
     }
     let v = eVec.angle(pos);
     if (pos.dot(vel) < 0) {
@@ -16625,8 +16708,9 @@ var Satellite = class _Satellite extends SpaceObject {
     this.satrec = Sgp4.createSatrec(tle1, tle2);
   }
   parseOmmAndUpdateOrbit_(omm) {
-    this.sccNum = omm.NORAD_CAT_ID.padStart(5, "0");
-    this.sccNum5 = Tle.convert6DigitToA5(omm.NORAD_CAT_ID);
+    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.intlDes = omm.OBJECT_ID;
     const YYYY = omm.EPOCH.slice(0, 4);
@@ -16639,6 +16723,19 @@ var Satellite = class _Satellite extends SpaceObject {
     const dayOfYear = (epochDateObj - Date.UTC(Number(YYYY), 0, 0)) / 864e5;
     const ommParsed = {
       ...omm,
+      NORAD_CAT_ID: noradStr,
+      MEAN_MOTION: String(omm.MEAN_MOTION),
+      ECCENTRICITY: String(omm.ECCENTRICITY),
+      INCLINATION: String(omm.INCLINATION),
+      RA_OF_ASC_NODE: String(omm.RA_OF_ASC_NODE),
+      ARG_OF_PERICENTER: String(omm.ARG_OF_PERICENTER),
+      MEAN_ANOMALY: String(omm.MEAN_ANOMALY),
+      EPHEMERIS_TYPE: String(omm.EPHEMERIS_TYPE),
+      ELEMENT_SET_NO: String(omm.ELEMENT_SET_NO),
+      REV_AT_EPOCH: String(omm.REV_AT_EPOCH),
+      BSTAR: String(omm.BSTAR),
+      MEAN_MOTION_DOT: String(omm.MEAN_MOTION_DOT),
+      MEAN_MOTION_DDOT: String(omm.MEAN_MOTION_DDOT),
       epoch: {
         year: Number(YYYY),
         month: Number(MM),
@@ -16651,20 +16748,41 @@ var Satellite = class _Satellite extends SpaceObject {
     };
     this.epochYear = parseInt(YYYY.slice(2, 4));
     this.epochDay = dayOfYear;
-    this.meanMoDev1 = parseFloat(omm.MEAN_MOTION_DOT);
-    this.meanMoDev2 = parseFloat(omm.MEAN_MOTION_DDOT);
-    this.bstar = parseFloat(omm.BSTAR);
-    this.inclination = parseFloat(omm.INCLINATION);
-    this.rightAscension = parseFloat(omm.RA_OF_ASC_NODE);
-    this.eccentricity = parseFloat(omm.ECCENTRICITY);
-    this.argOfPerigee = parseFloat(omm.ARG_OF_PERICENTER);
-    this.meanAnomaly = parseFloat(omm.MEAN_ANOMALY);
-    this.meanMotion = parseFloat(omm.MEAN_MOTION);
+    this.meanMoDev1 = Number(omm.MEAN_MOTION_DOT);
+    this.meanMoDev2 = Number(omm.MEAN_MOTION_DDOT);
+    this.bstar = Number(omm.BSTAR);
+    this.inclination = Number(omm.INCLINATION);
+    this.rightAscension = Number(omm.RA_OF_ASC_NODE);
+    this.eccentricity = Number(omm.ECCENTRICITY);
+    this.argOfPerigee = Number(omm.ARG_OF_PERICENTER);
+    this.meanAnomaly = Number(omm.MEAN_ANOMALY);
+    this.meanMotion = Number(omm.MEAN_MOTION);
     this.period = 1440 / this.meanMotion;
     this.semiMajorAxis = (8681663653e-3 / this.meanMotion) ** (2 / 3);
     this.semiMinorAxis = this.semiMajorAxis * Math.sqrt(1 - this.eccentricity ** 2);
     this.apogee = this.semiMajorAxis * (1 + this.eccentricity) - 6371;
     this.perigee = this.semiMajorAxis * (1 - this.eccentricity) - 6371;
+    const { tle1, tle2 } = FormatTle.createTle({
+      inc: this.inclination,
+      meanmo: this.meanMotion,
+      rasc: this.rightAscension,
+      argPe: this.argOfPerigee,
+      meana: this.meanAnomaly,
+      ecen: this.eccentricity,
+      epochyr: this.epochYear,
+      epochday: this.epochDay,
+      intl: omm.OBJECT_ID,
+      scc: this.sccNum,
+      bstar: this.bstar,
+      meanMotionDot: this.meanMoDev1,
+      meanMotionDdot: this.meanMoDev2,
+      classification: String(omm.CLASSIFICATION_TYPE) || "U",
+      revAtEpoch: Number(omm.REV_AT_EPOCH) || 0,
+      elementSetNo: Number(omm.ELEMENT_SET_NO) || 999,
+      ephemerisType: Number(omm.EPHEMERIS_TYPE) || 0
+    });
+    this.tle1 = tle1;
+    this.tle2 = tle2;
     this.satrec = Sgp4.createSatrecFromOmm(ommParsed);
   }
   // ==================== Type Checking ====================
@@ -25180,6 +25298,61 @@ var OmmParser = class _OmmParser {
     const theory = omm.metadata.MEAN_ELEMENT_THEORY.toUpperCase();
     return theory === "SGP" || theory === "SGP4" || theory === "SGP4-XP";
   }
+  /**
+   * Parse a CelesTrak-style flat JSON OMM object into the structured ParsedOmm format.
+   *
+   * CelesTrak's JSON OMM format is a flat object with all fields at the top level,
+   * omitting metadata fields like CENTER_NAME, REF_FRAME, TIME_SYSTEM, and
+   * MEAN_ELEMENT_THEORY (which are always EARTH/TEME/UTC/SGP4 for GP data).
+   *
+   * @param omm - A flat OMM JSON object from CelesTrak
+   * @returns Parsed OMM structure with inferred metadata
+   */
+  static parseJson(omm) {
+    const result = {
+      header: {
+        CCSDS_OMM_VERS: "3.0",
+        CREATION_DATE: (/* @__PURE__ */ new Date()).toISOString(),
+        ORIGINATOR: "CelesTrak"
+      },
+      metadata: {
+        OBJECT_NAME: omm.OBJECT_NAME,
+        OBJECT_ID: omm.OBJECT_ID,
+        CENTER_NAME: "EARTH",
+        REF_FRAME: "TEME",
+        TIME_SYSTEM: "UTC",
+        MEAN_ELEMENT_THEORY: "SGP4"
+      },
+      meanElements: {
+        EPOCH: omm.EPOCH,
+        MEAN_MOTION: Number(omm.MEAN_MOTION),
+        ECCENTRICITY: Number(omm.ECCENTRICITY),
+        INCLINATION: Number(omm.INCLINATION),
+        RA_OF_ASC_NODE: Number(omm.RA_OF_ASC_NODE),
+        ARG_OF_PERICENTER: Number(omm.ARG_OF_PERICENTER),
+        MEAN_ANOMALY: Number(omm.MEAN_ANOMALY)
+      },
+      tleParameters: {
+        EPHEMERIS_TYPE: Number(omm.EPHEMERIS_TYPE),
+        CLASSIFICATION_TYPE: String(omm.CLASSIFICATION_TYPE),
+        NORAD_CAT_ID: Number(omm.NORAD_CAT_ID),
+        ELEMENT_SET_NO: Number(omm.ELEMENT_SET_NO),
+        REV_AT_EPOCH: Number(omm.REV_AT_EPOCH),
+        BSTAR: Number(omm.BSTAR),
+        MEAN_MOTION_DOT: Number(omm.MEAN_MOTION_DOT),
+        MEAN_MOTION_DDOT: Number(omm.MEAN_MOTION_DDOT)
+      }
+    };
+    return result;
+  }
+  /**
+   * Parse an array of CelesTrak-style flat JSON OMM objects.
+   * @param ommArray - Array of flat OMM JSON objects from CelesTrak
+   * @returns Array of parsed OMM structures
+   */
+  static parseJsonArray(ommArray) {
+    return ommArray.map((omm) => _OmmParser.parseJson(omm));
+  }
   /**
    * Parse all key-value pairs from OMM lines (excluding COMMENT lines).
    */

package/dist/main.d.cts

@@ -1816,6 +1816,27 @@ declare abstract class FormatTle {
      * @returns The modified string with the character set at the specified index.
      */
     static setCharAt(str: string, index: number, chr: string): string;
+    /**
+     * Format mean motion dot (first derivative / 2) for TLE line 1.
+     * Format: sign + ".NNNNNNNN" = 10 chars total.
+     * @param value - Mean motion dot value (rev/day^2)
+     * @returns Formatted 10-character string
+     */
+    static formatMeanMotionDot(value: number): string;
+    /**
+     * Format a value in TLE exponential notation for BSTAR or mean motion ddot.
+     * Format: "sNNNNN±N" (8 chars) where mantissa has implied leading decimal point.
+     * Example: 0.00017507 → " 17507-3" (i.e., .17507 × 10^-3)
+     * @param value - The value to format
+     * @returns Formatted 8-character string
+     */
+    static formatTleExponential(value: number): string;
+    /**
+     * Compute TLE line checksum (modulo 10 sum of digits, '-' counts as 1).
+     * @param line - TLE line (first 68 characters are summed)
+     * @returns Checksum digit (0-9)
+     */
+    static tleChecksum(line: string): number;
 }
 
 /**
@@ -2780,37 +2801,34 @@ declare class RIC extends RelativeState {
 }
 
 /**
- * Represents the data format for orbital elements as provided by the Omm system.
- * Each property is expressed as a string to accommodate various formatting or precision
- * requirements from the data source.
+ * Represents the data format for orbital elements as provided by the OMM system.
+ * Numeric fields accept both string and number to support CelesTrak's JSON format
+ * (which sends numbers as numbers) and other sources that send them as strings.
  */
 interface OmmDataFormat {
     OBJECT_NAME: string;
     OBJECT_ID: string;
     /** Date in YYYY-MM-DDTHH:MM:SS.SSSSSS UTC format */
     EPOCH: string;
-    MEAN_MOTION: string;
-    ECCENTRICITY: string;
-    INCLINATION: string;
-    RA_OF_ASC_NODE: string;
-    ARG_OF_PERICENTER: string;
-    MEAN_ANOMALY: string;
-    EPHEMERIS_TYPE: string;
+    MEAN_MOTION: string | number;
+    ECCENTRICITY: string | number;
+    INCLINATION: string | number;
+    RA_OF_ASC_NODE: string | number;
+    ARG_OF_PERICENTER: string | number;
+    MEAN_ANOMALY: string | number;
+    EPHEMERIS_TYPE: string | number;
     CLASSIFICATION_TYPE: string;
-    NORAD_CAT_ID: string;
-    ELEMENT_SET_NO: string;
-    REV_AT_EPOCH: string;
-    BSTAR: string;
-    MEAN_MOTION_DOT: string;
-    MEAN_MOTION_DDOT: string;
+    NORAD_CAT_ID: string | number;
+    ELEMENT_SET_NO: string | number;
+    REV_AT_EPOCH: string | number;
+    BSTAR: string | number;
+    MEAN_MOTION_DOT: string | number;
+    MEAN_MOTION_DDOT: string | number;
 }
 /**
- * Represents the parsed data format for orbital elements as provided by the Omm system.
- * Each property is expressed as a string to accommodate various formatting or precision
- * requirements from the data source.
- * The `epoch` property is an object that contains the parsed date and time values.
- * The `doy` property is the day of the year.
- * The `epoch` property is an object that contains the parsed date and time values.
+ * Represents the parsed data format for orbital elements as provided by the OMM system.
+ * String fields are preserved from the original data; the `epoch` property contains
+ * the parsed date/time values used for SGP4 initialization.
  */
 interface OmmParsedDataFormat {
     OBJECT_NAME: string;
@@ -5919,6 +5937,20 @@ type TleParams = {
     intl: string;
     /** alpha 5 satellite number */
     scc: string;
+    /** B* drag term (1/Earth radii). Used when `sat` is not provided. */
+    bstar?: number;
+    /** First derivative of mean motion / 2 (rev/day^2). Used when `sat` is not provided. */
+    meanMotionDot?: number;
+    /** Second derivative of mean motion / 6 (rev/day^3). Used when `sat` is not provided. */
+    meanMotionDdot?: number;
+    /** Classification type (default 'U'). Used when `sat` is not provided. */
+    classification?: string;
+    /** Revolution number at epoch. Used when `sat` is not provided. */
+    revAtEpoch?: number;
+    /** Element set number (default 999). Used when `sat` is not provided. */
+    elementSetNo?: number;
+    /** Ephemeris type (default 0). Used when `sat` is not provided. */
+    ephemerisType?: number;
 };
 type PositionVelocity = {
     position: Vector3D<Kilometers>;
@@ -14078,6 +14110,23 @@ declare class OmmParser {
      * @returns true if the OMM uses SGP, SGP4, or SGP4-XP mean element theory
      */
     static isTleCompatible(omm: ParsedOmm): boolean;
+    /**
+     * Parse a CelesTrak-style flat JSON OMM object into the structured ParsedOmm format.
+     *
+     * CelesTrak's JSON OMM format is a flat object with all fields at the top level,
+     * omitting metadata fields like CENTER_NAME, REF_FRAME, TIME_SYSTEM, and
+     * MEAN_ELEMENT_THEORY (which are always EARTH/TEME/UTC/SGP4 for GP data).
+     *
+     * @param omm - A flat OMM JSON object from CelesTrak
+     * @returns Parsed OMM structure with inferred metadata
+     */
+    static parseJson(omm: OmmDataFormat): ParsedOmm;
+    /**
+     * Parse an array of CelesTrak-style flat JSON OMM objects.
+     * @param ommArray - Array of flat OMM JSON objects from CelesTrak
+     * @returns Array of parsed OMM structures
+     */
+    static parseJsonArray(ommArray: OmmDataFormat[]): ParsedOmm[];
     /**
      * Parse all key-value pairs from OMM lines (excluding COMMENT lines).
      */

package/dist/main.d.ts

@@ -1816,6 +1816,27 @@ declare abstract class FormatTle {
      * @returns The modified string with the character set at the specified index.
      */
     static setCharAt(str: string, index: number, chr: string): string;
+    /**
+     * Format mean motion dot (first derivative / 2) for TLE line 1.
+     * Format: sign + ".NNNNNNNN" = 10 chars total.
+     * @param value - Mean motion dot value (rev/day^2)
+     * @returns Formatted 10-character string
+     */
+    static formatMeanMotionDot(value: number): string;
+    /**
+     * Format a value in TLE exponential notation for BSTAR or mean motion ddot.
+     * Format: "sNNNNN±N" (8 chars) where mantissa has implied leading decimal point.
+     * Example: 0.00017507 → " 17507-3" (i.e., .17507 × 10^-3)
+     * @param value - The value to format
+     * @returns Formatted 8-character string
+     */
+    static formatTleExponential(value: number): string;
+    /**
+     * Compute TLE line checksum (modulo 10 sum of digits, '-' counts as 1).
+     * @param line - TLE line (first 68 characters are summed)
+     * @returns Checksum digit (0-9)
+     */
+    static tleChecksum(line: string): number;
 }
 
 /**
@@ -2780,37 +2801,34 @@ declare class RIC extends RelativeState {
 }
 
 /**
- * Represents the data format for orbital elements as provided by the Omm system.
- * Each property is expressed as a string to accommodate various formatting or precision
- * requirements from the data source.
+ * Represents the data format for orbital elements as provided by the OMM system.
+ * Numeric fields accept both string and number to support CelesTrak's JSON format
+ * (which sends numbers as numbers) and other sources that send them as strings.
  */
 interface OmmDataFormat {
     OBJECT_NAME: string;
     OBJECT_ID: string;
     /** Date in YYYY-MM-DDTHH:MM:SS.SSSSSS UTC format */
     EPOCH: string;
-    MEAN_MOTION: string;
-    ECCENTRICITY: string;
-    INCLINATION: string;
-    RA_OF_ASC_NODE: string;
-    ARG_OF_PERICENTER: string;
-    MEAN_ANOMALY: string;
-    EPHEMERIS_TYPE: string;
+    MEAN_MOTION: string | number;
+    ECCENTRICITY: string | number;
+    INCLINATION: string | number;
+    RA_OF_ASC_NODE: string | number;
+    ARG_OF_PERICENTER: string | number;
+    MEAN_ANOMALY: string | number;
+    EPHEMERIS_TYPE: string | number;
     CLASSIFICATION_TYPE: string;
-    NORAD_CAT_ID: string;
-    ELEMENT_SET_NO: string;
-    REV_AT_EPOCH: string;
-    BSTAR: string;
-    MEAN_MOTION_DOT: string;
-    MEAN_MOTION_DDOT: string;
+    NORAD_CAT_ID: string | number;
+    ELEMENT_SET_NO: string | number;
+    REV_AT_EPOCH: string | number;
+    BSTAR: string | number;
+    MEAN_MOTION_DOT: string | number;
+    MEAN_MOTION_DDOT: string | number;
 }
 /**
- * Represents the parsed data format for orbital elements as provided by the Omm system.
- * Each property is expressed as a string to accommodate various formatting or precision
- * requirements from the data source.
- * The `epoch` property is an object that contains the parsed date and time values.
- * The `doy` property is the day of the year.
- * The `epoch` property is an object that contains the parsed date and time values.
+ * Represents the parsed data format for orbital elements as provided by the OMM system.
+ * String fields are preserved from the original data; the `epoch` property contains
+ * the parsed date/time values used for SGP4 initialization.
  */
 interface OmmParsedDataFormat {
     OBJECT_NAME: string;
@@ -5919,6 +5937,20 @@ type TleParams = {
     intl: string;
     /** alpha 5 satellite number */
     scc: string;
+    /** B* drag term (1/Earth radii). Used when `sat` is not provided. */
+    bstar?: number;
+    /** First derivative of mean motion / 2 (rev/day^2). Used when `sat` is not provided. */
+    meanMotionDot?: number;
+    /** Second derivative of mean motion / 6 (rev/day^3). Used when `sat` is not provided. */
+    meanMotionDdot?: number;
+    /** Classification type (default 'U'). Used when `sat` is not provided. */
+    classification?: string;
+    /** Revolution number at epoch. Used when `sat` is not provided. */
+    revAtEpoch?: number;
+    /** Element set number (default 999). Used when `sat` is not provided. */
+    elementSetNo?: number;
+    /** Ephemeris type (default 0). Used when `sat` is not provided. */
+    ephemerisType?: number;
 };
 type PositionVelocity = {
     position: Vector3D<Kilometers>;
@@ -14078,6 +14110,23 @@ declare class OmmParser {
      * @returns true if the OMM uses SGP, SGP4, or SGP4-XP mean element theory
      */
     static isTleCompatible(omm: ParsedOmm): boolean;
+    /**
+     * Parse a CelesTrak-style flat JSON OMM object into the structured ParsedOmm format.
+     *
+     * CelesTrak's JSON OMM format is a flat object with all fields at the top level,
+     * omitting metadata fields like CENTER_NAME, REF_FRAME, TIME_SYSTEM, and
+     * MEAN_ELEMENT_THEORY (which are always EARTH/TEME/UTC/SGP4 for GP data).
+     *
+     * @param omm - A flat OMM JSON object from CelesTrak
+     * @returns Parsed OMM structure with inferred metadata
+     */
+    static parseJson(omm: OmmDataFormat): ParsedOmm;
+    /**
+     * Parse an array of CelesTrak-style flat JSON OMM objects.
+     * @param ommArray - Array of flat OMM JSON objects from CelesTrak
+     * @returns Array of parsed OMM structures
+     */
+    static parseJsonArray(ommArray: OmmDataFormat[]): ParsedOmm[];
     /**
      * Parse all key-value pairs from OMM lines (excluding COMMENT lines).
      */

package/dist/main.js

@@ -9847,7 +9847,7 @@ var FormatTle = class _FormatTle {
     const { inc, meanmo, rasc, argPe, meana, ecen, epochyr, epochday, intl } = tleParams;
     const scc = Tle.convert6DigitToA5(tleParams.scc);
     const epochYrStr = String(epochyr).padStart(2, "0");
-    const epochdayStr = parseFloat(String(epochday)).toFixed(8).padStart(12, "0");
+    const epochdayStr = Number(epochday).toFixed(8).padStart(12, "0");
     const incStr = _FormatTle.inclination(inc);
     const meanmoStr = _FormatTle.meanMotion(meanmo);
     const rascStr = _FormatTle.rightAscension(rasc);
@@ -9855,13 +9855,23 @@ var FormatTle = class _FormatTle {
     const meanaStr = _FormatTle.meanAnomaly(meana);
     const ecenStr = _FormatTle.eccentricity(ecen);
     const intlStr = intl.padEnd(8, " ");
-    let TLE1Ending = tleParams.sat ? tleParams.sat.tle1.substring(32, 71) : " +.00000000 +00000+0 +00000-0 0  9990";
-    TLE1Ending = TLE1Ending[1] === " " ? _FormatTle.setCharAt(TLE1Ending, 1, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[12] === " " ? _FormatTle.setCharAt(TLE1Ending, 12, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[21] === " " ? _FormatTle.setCharAt(TLE1Ending, 21, "+") : TLE1Ending;
-    TLE1Ending = TLE1Ending[32] === " " ? _FormatTle.setCharAt(TLE1Ending, 32, "0") : TLE1Ending;
-    const tle1 = `1 ${scc}U ${intlStr} ${epochYrStr}${epochdayStr}${TLE1Ending}`;
-    const tle2 = `2 ${scc} ${incStr} ${rascStr} ${ecenStr} ${argPeStr} ${meanaStr} ${meanmoStr} 00010`;
+    const classification = tleParams.classification ?? "U";
+    const elementSetNo = (tleParams.elementSetNo ?? 999).toString().padStart(4, " ");
+    const ephemerisType = (tleParams.ephemerisType ?? 0).toString();
+    const revAtEpoch = (tleParams.revAtEpoch ?? 0).toString().padStart(5, " ");
+    let TLE1Ending;
+    if (tleParams.sat) {
+      TLE1Ending = tleParams.sat.tle1.substring(32, 71);
+    } else {
+      const mmDot = _FormatTle.formatMeanMotionDot(tleParams.meanMotionDot ?? 0);
+      const mmDdot = _FormatTle.formatTleExponential(tleParams.meanMotionDdot ?? 0);
+      const bstarStr = _FormatTle.formatTleExponential(tleParams.bstar ?? 0);
+      TLE1Ending = ` ${mmDot} ${mmDdot} ${bstarStr} ${ephemerisType} ${elementSetNo}0`;
+    }
+    const tle1Pre = `1 ${scc}${classification} ${intlStr} ${epochYrStr}${epochdayStr}${TLE1Ending}`;
+    const tle1 = _FormatTle.setCharAt(tle1Pre, 68, _FormatTle.tleChecksum(tle1Pre).toString());
+    const tle2Pre = `2 ${scc} ${incStr} ${rascStr} ${ecenStr} ${argPeStr} ${meanaStr} ${meanmoStr}${revAtEpoch}0`;
+    const tle2 = _FormatTle.setCharAt(tle2Pre, 68, _FormatTle.tleChecksum(tle2Pre).toString());
     return { tle1, tle2 };
   }
   /**
@@ -9875,7 +9885,7 @@ var FormatTle = class _FormatTle {
       argPe = argPe.toString();
     }
     const argPeNum = parseFloat(argPe).toFixed(4);
-    const argPe0 = argPeNum.padStart(8, "0");
+    const argPe0 = argPeNum.padStart(8, " ");
     if (argPe0.length !== 8) {
       throw new ValidationError("Argument of perigee must be 8 characters", "argPe", argPe0);
     }
@@ -9889,7 +9899,10 @@ var FormatTle = class _FormatTle {
    */
   static eccentricity(ecen) {
     if (typeof ecen === "number") {
-      ecen = ecen.toFixed(7);
+      const ecenStr = ecen.toString();
+      const dotIdx = ecenStr.indexOf(".");
+      const afterDot = dotIdx >= 0 ? ecenStr.substring(dotIdx + 1) : "";
+      ecen = `0.${afterDot.substring(0, 7).padEnd(7, "0")}`;
     }
     let ecen0 = ecen.padEnd(9, "0");
     if (ecen0[1] === ".") {
@@ -9913,7 +9926,7 @@ var FormatTle = class _FormatTle {
       inc = inc.toString();
     }
     const incNum = parseFloat(inc).toFixed(4);
-    const inc0 = incNum.padStart(8, "0");
+    const inc0 = incNum.padStart(8, " ");
     if (inc0.length !== 8) {
       throw new ValidationError("Inclination must be 8 characters", "inclination", inc0);
     }
@@ -9930,7 +9943,7 @@ var FormatTle = class _FormatTle {
       meana = meana.toString();
     }
     const meanaNum = parseFloat(meana).toFixed(4);
-    const meana0 = meanaNum.padStart(8, "0");
+    const meana0 = meanaNum.padStart(8, " ");
     if (meana0.length !== 8) {
       throw new ValidationError("Mean anomaly must be 8 characters", "meanAnomaly", meana0);
     }
@@ -9971,7 +9984,7 @@ var FormatTle = class _FormatTle {
       rasc = rasc.toString();
     }
     const rascNum = parseFloat(rasc).toFixed(4);
-    const rasc0 = rascNum.padStart(8, "0");
+    const rasc0 = rascNum.padStart(8, " ");
     if (rasc0.length !== 8) {
       throw new ValidationError("Right ascension must be 8 characters", "rightAscension", rasc0);
     }
@@ -9990,6 +10003,61 @@ var FormatTle = class _FormatTle {
     }
     return `${str.substring(0, index)}${chr}${str.substring(index + 1)}`;
   }
+  /**
+   * Format mean motion dot (first derivative / 2) for TLE line 1.
+   * Format: sign + ".NNNNNNNN" = 10 chars total.
+   * @param value - Mean motion dot value (rev/day^2)
+   * @returns Formatted 10-character string
+   */
+  static formatMeanMotionDot(value) {
+    const sign2 = value >= 0 ? " " : "-";
+    return sign2 + Math.abs(value).toFixed(8).substring(1);
+  }
+  /**
+   * Format a value in TLE exponential notation for BSTAR or mean motion ddot.
+   * Format: "sNNNNN±N" (8 chars) where mantissa has implied leading decimal point.
+   * Example: 0.00017507 → " 17507-3" (i.e., .17507 × 10^-3)
+   * @param value - The value to format
+   * @returns Formatted 8-character string
+   */
+  static formatTleExponential(value) {
+    if (value === 0) {
+      return " 00000+0";
+    }
+    const sign2 = value >= 0 ? " " : "-";
+    const absVal = Math.abs(value);
+    let exponent = Math.floor(Math.log10(absVal)) + 1;
+    let mantissa = absVal / Math.pow(10, exponent);
+    if (mantissa >= 1) {
+      mantissa /= 10;
+      exponent++;
+    }
+    if (mantissa < 0.1 && mantissa > 0) {
+      mantissa *= 10;
+      exponent--;
+    }
+    const mantissaStr = Math.round(mantissa * 1e5).toString().padStart(5, "0");
+    const expSign = exponent >= 0 ? "+" : "-";
+    const expStr = Math.abs(exponent).toString();
+    return sign2 + mantissaStr + expSign + expStr;
+  }
+  /**
+   * Compute TLE line checksum (modulo 10 sum of digits, '-' counts as 1).
+   * @param line - TLE line (first 68 characters are summed)
+   * @returns Checksum digit (0-9)
+   */
+  static tleChecksum(line) {
+    let sum = 0;
+    for (let i = 0; i < 68 && i < line.length; i++) {
+      const c = line[i];
+      if (c >= "0" && c <= "9") {
+        sum += parseInt(c);
+      } else if (c === "-") {
+        sum += 1;
+      }
+    }
+    return sum % 10;
+  }
 };
 
 // src/coordinate/RelativeState.ts
@@ -13185,13 +13253,28 @@ var ClassicalElements = class _ClassicalElements {
     const h = pos.cross(vel);
     const i = Math.acos(clamp(h.z / h.magnitude(), -1, 1));
     const n = Vector3D.zAxis.cross(h);
-    let o = Math.acos(clamp(n.x / n.magnitude(), -1, 1));
-    if (n.y < 0) {
-      o = TAU - o;
-    }
-    let w = n.angle(eVec);
-    if (eVec.z < 0) {
-      w = TAU - w;
+    const nMag = n.magnitude();
+    let o;
+    let w;
+    if (nMag < 1e-12) {
+      o = 0;
+      if (e > 1e-12) {
+        w = Math.atan2(eVec.y, eVec.x);
+        if (w < 0) {
+          w = w + TAU;
+        }
+      } else {
+        w = 0;
+      }
+    } else {
+      o = Math.acos(clamp(n.x / nMag, -1, 1));
+      if (n.y < 0) {
+        o = TAU - o;
+      }
+      w = n.angle(eVec);
+      if (eVec.z < 0) {
+        w = TAU - w;
+      }
     }
     let v = eVec.angle(pos);
     if (pos.dot(vel) < 0) {
@@ -16623,8 +16706,9 @@ var Satellite = class _Satellite extends SpaceObject {
     this.satrec = Sgp4.createSatrec(tle1, tle2);
   }
   parseOmmAndUpdateOrbit_(omm) {
-    this.sccNum = omm.NORAD_CAT_ID.padStart(5, "0");
-    this.sccNum5 = Tle.convert6DigitToA5(omm.NORAD_CAT_ID);
+    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.intlDes = omm.OBJECT_ID;
     const YYYY = omm.EPOCH.slice(0, 4);
@@ -16637,6 +16721,19 @@ var Satellite = class _Satellite extends SpaceObject {
     const dayOfYear = (epochDateObj - Date.UTC(Number(YYYY), 0, 0)) / 864e5;
     const ommParsed = {
       ...omm,
+      NORAD_CAT_ID: noradStr,
+      MEAN_MOTION: String(omm.MEAN_MOTION),
+      ECCENTRICITY: String(omm.ECCENTRICITY),
+      INCLINATION: String(omm.INCLINATION),
+      RA_OF_ASC_NODE: String(omm.RA_OF_ASC_NODE),
+      ARG_OF_PERICENTER: String(omm.ARG_OF_PERICENTER),
+      MEAN_ANOMALY: String(omm.MEAN_ANOMALY),
+      EPHEMERIS_TYPE: String(omm.EPHEMERIS_TYPE),
+      ELEMENT_SET_NO: String(omm.ELEMENT_SET_NO),
+      REV_AT_EPOCH: String(omm.REV_AT_EPOCH),
+      BSTAR: String(omm.BSTAR),
+      MEAN_MOTION_DOT: String(omm.MEAN_MOTION_DOT),
+      MEAN_MOTION_DDOT: String(omm.MEAN_MOTION_DDOT),
       epoch: {
         year: Number(YYYY),
         month: Number(MM),
@@ -16649,20 +16746,41 @@ var Satellite = class _Satellite extends SpaceObject {
     };
     this.epochYear = parseInt(YYYY.slice(2, 4));
     this.epochDay = dayOfYear;
-    this.meanMoDev1 = parseFloat(omm.MEAN_MOTION_DOT);
-    this.meanMoDev2 = parseFloat(omm.MEAN_MOTION_DDOT);
-    this.bstar = parseFloat(omm.BSTAR);
-    this.inclination = parseFloat(omm.INCLINATION);
-    this.rightAscension = parseFloat(omm.RA_OF_ASC_NODE);
-    this.eccentricity = parseFloat(omm.ECCENTRICITY);
-    this.argOfPerigee = parseFloat(omm.ARG_OF_PERICENTER);
-    this.meanAnomaly = parseFloat(omm.MEAN_ANOMALY);
-    this.meanMotion = parseFloat(omm.MEAN_MOTION);
+    this.meanMoDev1 = Number(omm.MEAN_MOTION_DOT);
+    this.meanMoDev2 = Number(omm.MEAN_MOTION_DDOT);
+    this.bstar = Number(omm.BSTAR);
+    this.inclination = Number(omm.INCLINATION);
+    this.rightAscension = Number(omm.RA_OF_ASC_NODE);
+    this.eccentricity = Number(omm.ECCENTRICITY);
+    this.argOfPerigee = Number(omm.ARG_OF_PERICENTER);
+    this.meanAnomaly = Number(omm.MEAN_ANOMALY);
+    this.meanMotion = Number(omm.MEAN_MOTION);
     this.period = 1440 / this.meanMotion;
     this.semiMajorAxis = (8681663653e-3 / this.meanMotion) ** (2 / 3);
     this.semiMinorAxis = this.semiMajorAxis * Math.sqrt(1 - this.eccentricity ** 2);
     this.apogee = this.semiMajorAxis * (1 + this.eccentricity) - 6371;
     this.perigee = this.semiMajorAxis * (1 - this.eccentricity) - 6371;
+    const { tle1, tle2 } = FormatTle.createTle({
+      inc: this.inclination,
+      meanmo: this.meanMotion,
+      rasc: this.rightAscension,
+      argPe: this.argOfPerigee,
+      meana: this.meanAnomaly,
+      ecen: this.eccentricity,
+      epochyr: this.epochYear,
+      epochday: this.epochDay,
+      intl: omm.OBJECT_ID,
+      scc: this.sccNum,
+      bstar: this.bstar,
+      meanMotionDot: this.meanMoDev1,
+      meanMotionDdot: this.meanMoDev2,
+      classification: String(omm.CLASSIFICATION_TYPE) || "U",
+      revAtEpoch: Number(omm.REV_AT_EPOCH) || 0,
+      elementSetNo: Number(omm.ELEMENT_SET_NO) || 999,
+      ephemerisType: Number(omm.EPHEMERIS_TYPE) || 0
+    });
+    this.tle1 = tle1;
+    this.tle2 = tle2;
     this.satrec = Sgp4.createSatrecFromOmm(ommParsed);
   }
   // ==================== Type Checking ====================
@@ -25178,6 +25296,61 @@ var OmmParser = class _OmmParser {
     const theory = omm.metadata.MEAN_ELEMENT_THEORY.toUpperCase();
     return theory === "SGP" || theory === "SGP4" || theory === "SGP4-XP";
   }
+  /**
+   * Parse a CelesTrak-style flat JSON OMM object into the structured ParsedOmm format.
+   *
+   * CelesTrak's JSON OMM format is a flat object with all fields at the top level,
+   * omitting metadata fields like CENTER_NAME, REF_FRAME, TIME_SYSTEM, and
+   * MEAN_ELEMENT_THEORY (which are always EARTH/TEME/UTC/SGP4 for GP data).
+   *
+   * @param omm - A flat OMM JSON object from CelesTrak
+   * @returns Parsed OMM structure with inferred metadata
+   */
+  static parseJson(omm) {
+    const result = {
+      header: {
+        CCSDS_OMM_VERS: "3.0",
+        CREATION_DATE: (/* @__PURE__ */ new Date()).toISOString(),
+        ORIGINATOR: "CelesTrak"
+      },
+      metadata: {
+        OBJECT_NAME: omm.OBJECT_NAME,
+        OBJECT_ID: omm.OBJECT_ID,
+        CENTER_NAME: "EARTH",
+        REF_FRAME: "TEME",
+        TIME_SYSTEM: "UTC",
+        MEAN_ELEMENT_THEORY: "SGP4"
+      },
+      meanElements: {
+        EPOCH: omm.EPOCH,
+        MEAN_MOTION: Number(omm.MEAN_MOTION),
+        ECCENTRICITY: Number(omm.ECCENTRICITY),
+        INCLINATION: Number(omm.INCLINATION),
+        RA_OF_ASC_NODE: Number(omm.RA_OF_ASC_NODE),
+        ARG_OF_PERICENTER: Number(omm.ARG_OF_PERICENTER),
+        MEAN_ANOMALY: Number(omm.MEAN_ANOMALY)
+      },
+      tleParameters: {
+        EPHEMERIS_TYPE: Number(omm.EPHEMERIS_TYPE),
+        CLASSIFICATION_TYPE: String(omm.CLASSIFICATION_TYPE),
+        NORAD_CAT_ID: Number(omm.NORAD_CAT_ID),
+        ELEMENT_SET_NO: Number(omm.ELEMENT_SET_NO),
+        REV_AT_EPOCH: Number(omm.REV_AT_EPOCH),
+        BSTAR: Number(omm.BSTAR),
+        MEAN_MOTION_DOT: Number(omm.MEAN_MOTION_DOT),
+        MEAN_MOTION_DDOT: Number(omm.MEAN_MOTION_DDOT)
+      }
+    };
+    return result;
+  }
+  /**
+   * Parse an array of CelesTrak-style flat JSON OMM objects.
+   * @param ommArray - Array of flat OMM JSON objects from CelesTrak
+   * @returns Array of parsed OMM structures
+   */
+  static parseJsonArray(ommArray) {
+    return ommArray.map((omm) => _OmmParser.parseJson(omm));
+  }
   /**
    * Parse all key-value pairs from OMM lines (excluding COMMENT lines).
    */

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "ootk",
-  "version": "6.0.1",
+  "version": "6.0.3",
   "type": "module",
   "description": "Orbital Object Toolkit including Multiple Propagators, Initial Orbit Determination, and Maneuver Calculations.",
   "exports": {