ootk 4.0.9 → 5.1.0

package/src/coordinate/index.ts

@@ -1 +1,12 @@
+export * from './ClassicalElements.js';
+export * from './EquinoctialElements.js';
+export * from './FormatTle.js';
+export * from './Geodetic.js';
+export * from './ITRF.js';
+export * from './J2000.js';
+export * from './RelativeState.js';
+export * from './RIC.js';
+export * from './StateVector.js';
+export * from './TEME.js';
+export * from './Tle.js';
 export { Hill } from './Hill.js';

package/src/coordinate/tle-format-data.ts

@@ -0,0 +1,47 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+/**
+ * Represents the format data of a TLE (Two-Line Element) set. This is used to
+ * make it easier to remember the starting and ending positions of the columns
+ * containing the TLE data.
+ */
+export class TleFormatData {
+  /** The starting position of the TLE data in the source string. */
+  start: number;
+  /** The ending position of the TLE data in the source string. */
+  stop: number;
+  /** The length of the TLE data in the source string. */
+  length: number;
+
+  /**
+   * Creates a new instance of TleFormatData.
+   * @param start The starting position of the TLE data in the source string.
+   * @param end The ending position of the TLE data in the source string.
+   */
+  constructor(start: number, end: number) {
+    this.start = start - 1;
+    this.stop = end;
+    this.length = this.stop - this.start;
+  }
+}

package/src/covariance/CovarianceSample.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
  * terms of the GNU Affero General Public License as published by the Free Software
@@ -23,11 +23,13 @@ import {
   Kilometers,
   KilometersPerSecond,
   Matrix,
+  OrbitRegime,
   RelativeState,
   RIC,
+  Tle, Vec3Flat,
   Vector,
   Vector3D,
-} from 'ootk-core';
+} from '../main.js';
 import { ForceModel } from '../force/ForceModel.js';
 import { Thrust } from '../force/Thrust.js';
 import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator.js';
@@ -35,9 +37,9 @@ import { CovarianceFrame, StateCovariance } from './StateCovariance.js';
 
 // / Sigma point covariance sample.
 export class CovarianceSample {
-  private _origin: RungeKutta89Propagator;
-  private _samples: RungeKutta89Propagator[] = [];
-  private matrix_ = Matrix.zero(6, 12);
+  private readonly origin_: RungeKutta89Propagator;
+  private readonly samples_: RungeKutta89Propagator[] = [];
+  private readonly matrix_ = Matrix.zero(6, 12);
 
   /**
    * Create a new [CovarianceSample] object from an inertial state, covariance
@@ -46,19 +48,45 @@ export class CovarianceSample {
    * Two-body physics will be used if a force model is not provided.
    * @param state The origin state.
    * @param covariance The covariance.
+   * @param tle The TLE object.
    * @param originForceModel The force model for the origin state.
    * @param sampleForceModel The force model for the samples.
    */
-  constructor(state: J2000, covariance: StateCovariance, originForceModel?: ForceModel, sampleForceModel?: ForceModel) {
+  constructor(
+    state: J2000, covariance: StateCovariance, tle?: Tle, originForceModel?: ForceModel, sampleForceModel?: ForceModel,
+  ) {
     originForceModel ??= new ForceModel().setGravity();
     sampleForceModel ??= new ForceModel().setGravity();
-    this._origin = new RungeKutta89Propagator(state, originForceModel);
+
+    // Scale covariance using TLE quality and regime aging factor if TLE info is provided
+    let scale = [1, 1, 1];
+
+    if (tle) {
+      const tleAgeDays = Tle.calcElsetAge(tle.line1, new Date(), 'days');
+
+      const quality = this.evaluateTleQuality(tle);
+      const aging = this.getRegimeAgingFactor(tle, tleAgeDays);
+
+      scale = [
+        quality[0] * aging[0],
+        quality[1] * aging[1],
+        quality[2] * aging[2],
+      ];
+    }
+
+    this.origin_ = new RungeKutta89Propagator(state, originForceModel);
     const s = covariance.matrix.cholesky().elements;
     const sqrt6 = Math.sqrt(6.0);
 
     for (let i = 0; i < 6; i++) {
       for (let j = 0; j < 6; j++) {
-        s[i][j] *= sqrt6;
+      /*
+       * Apply scale[0] to R, scale[1] to I, scale[2] to C (x, y, z)
+       * Position: i = 0,1,2; Velocity: i = 3,4,5
+       */
+        const scaleIdx = i % 3;
+
+        s[i][j] *= sqrt6 * scale[scaleIdx];
       }
     }
 
@@ -92,23 +120,23 @@ export class CovarianceSample {
       if (covariance.frame === CovarianceFrame.ECI) {
         const sample = new J2000(state.epoch, state.position.add(sampleR), state.velocity.add(sampleV));
 
-        this._samples.push(new RungeKutta89Propagator(sample, sampleForceModel));
+        this.samples_.push(new RungeKutta89Propagator(sample, sampleForceModel));
       } else if (covariance.frame === CovarianceFrame.RIC) {
         const sample = new RIC(sampleR, sampleV).toJ2000(state);
 
-        this._samples.push(new RungeKutta89Propagator(sample, sampleForceModel));
+        this.samples_.push(new RungeKutta89Propagator(sample, sampleForceModel));
       }
     }
   }
 
   // / Current covariance sample epoch.
   get epoch(): Epoch {
-    return this._origin.state.epoch;
+    return this.origin_.state.epoch;
   }
 
   // / Current covariance sample origin state.
   get state(): J2000 {
-    return this._origin.state;
+    return this.origin_.state;
   }
 
   // / Rebuild covariance from sigma points.
@@ -142,16 +170,16 @@ export class CovarianceSample {
 
   // / Propagate covariance to a new epoch.
   propagate(epoch: EpochUTC): void {
-    this._origin.propagate(epoch);
-    for (const sample of this._samples) {
+    this.origin_.propagate(epoch);
+    for (const sample of this.samples_) {
       sample.propagate(epoch);
     }
   }
 
   // / Apply a maneuver to this covariance.
   maneuver(maneuver: Thrust): void {
-    this._origin.maneuver(maneuver);
-    for (const sample of this._samples) {
+    this.origin_.maneuver(maneuver);
+    for (const sample of this.samples_) {
       sample.maneuver(maneuver);
     }
   }
@@ -159,7 +187,7 @@ export class CovarianceSample {
   // / Desample covariance in J2000 frame.
   desampleJ2000(): StateCovariance {
     for (let i = 0; i < 12; i++) {
-      const state = this._samples[i].state;
+      const state = this.samples_[i].state;
 
       this.matrix_.elements[0][i] = state.position.x;
       this.matrix_.elements[1][i] = state.position.y;
@@ -175,10 +203,10 @@ export class CovarianceSample {
 
   // / Desample covariance in RIC frame.
   desampleRIC(): StateCovariance {
-    const rot = RelativeState.createMatrix(this._origin.state.position, this._origin.state.velocity);
+    const rot = RelativeState.createMatrix(this.origin_.state.position, this.origin_.state.velocity);
 
     for (let i = 0; i < 12; i++) {
-      const state = RIC.fromJ2000Matrix(this._samples[i].state, this._origin.state, rot);
+      const state = RIC.fromJ2000Matrix(this.samples_[i].state, this.origin_.state, rot);
 
       this.matrix_.elements[0][i] = state.position.x;
       this.matrix_.elements[1][i] = state.position.y;
@@ -191,4 +219,74 @@ export class CovarianceSample {
 
     return new StateCovariance(matrix, CovarianceFrame.RIC);
   }
+
+  evaluateTleQuality(tle: Tle): Vec3Flat {
+    let c = 1,
+      i = 1,
+      r = 1; // start with nominal 120 / 1000 / 100 m
+
+    /* ---- Mean–motion first derivative (rev/day²) ---- */
+    const mm1 = Math.abs(Tle.meanMoDev1(tle.line1));
+
+    if (mm1 > 1e-3) { // clear manoeuvre or very low-drag orbit
+      i *= 3.0;
+      r *= 2.0;
+      c *= 1.2;
+    } else if (mm1 > 1e-5) { // high drag but likely passive
+      i *= 1.6;
+      r *= 1.3;
+    }
+
+    /* ---- BSTAR drag term ---- */
+    const bstar = Math.abs(Tle.bstar(tle.line1));
+
+    if (bstar > 5e-4) { // <≈400 km LEO
+      i *= 1.5;
+      r *= 1.5;
+      c *= 1.1;
+    } else if (bstar > 1e-4) { // 400–600 km
+      i *= 1.2;
+      r *= 1.2;
+    }
+
+    /* ---- Eccentricity ---- */
+    const e = tle.eccentricity;
+
+    if (e > 0.05) { // HEO, GTO, cubesat transfer, etc.
+      r *= 1 + 6 * e;
+      i *= 1 + 4 * e;
+      c *= 1 + 1 * e;
+    } else if (e > 0.02) { // mildly elliptical LEO/MEO
+      r *= 1 + 3 * e;
+      i *= 1 + 2 * e;
+    }
+
+    return [r, i, c];
+  }
+
+  getRegimeAgingFactor(tle: Tle, ageDays: number): Vec3Flat {
+    const regime = tle.state.toClassicalElements().getOrbitRegime();
+    const t = Math.max(ageDays, 0) ** 1.5; // non-linear growth
+
+    switch (regime) {
+
+      case OrbitRegime.LEO:
+      // Target ~×2 (R), ×3 (I), ×2.2 (C) after 1 day
+        return [1 + 1.0 * t, 1 + 2.0 * t, 1 + 1.2 * t];
+
+      case OrbitRegime.MEO:
+      // GNSS shells – slower growth
+        return [1 + 0.4 * t, 1 + 0.9 * t, 1 + 0.6 * t];
+
+      case OrbitRegime.GEO:
+        return [1 + 0.2 * t, 1 + 0.5 * t, 1 + 0.3 * t];
+
+      case OrbitRegime.HEO:
+      // Highly elliptical transfer or Molniya
+        return [1 + 1.2 * t, 1 + 2.4 * t, 1 + 1.4 * t];
+
+      default:
+        return [1 + 0.6 * t, 1 + 1.2 * t, 1 + 0.8 * t];
+    }
+  }
 }

package/src/covariance/StateCovariance.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
  * terms of the GNU Affero General Public License as published by the Free Software
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Matrix, Vector } from 'ootk-core';
+import { Matrix, Vector } from '../main.js';
 
 /** Covariance Frame */
 export enum CovarianceFrame {

package/src/covariance/index.ts

@@ -0,0 +1,2 @@
+export * from './CovarianceSample.js';
+export * from './StateCovariance.js';

package/src/data/DataHandler.ts

@@ -0,0 +1,88 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+import { egm96Data, Egm96Entry } from './values/Egm96Data.js';
+import { hpAtmosphereData } from './values/HpAtmosphereData.js';
+import { HpAtmosphereResult } from './values/HpAtmosphereResult.js';
+import { iau1980Data, Iau1980Entry } from './values/Iau1980Data.js';
+import { leapSecondData } from './values/LeapSecondData.js';
+
+/**
+ * Astrodynamic data management singleton.
+ *
+ * This class provides access to the various data sets used by the library.
+ * It is a singleton class, and can be accessed via the static
+ * [getInstance] method.
+ */
+export class DataHandler {
+  private static readonly instance_ = new DataHandler();
+
+  private constructor() {
+    // Prevent instantiation.
+  }
+
+  /**
+   * Returns the singleton instance of the DataHandler class.
+   * @returns The singleton instance of the DataHandler class.
+   */
+  static getInstance(): DataHandler {
+    return DataHandler.instance_;
+  }
+
+  /**
+   * Retrieves the Egm96 coefficients for the given l and m values.
+   * @param l - The degree of the coefficient.
+   * @param m - The order of the coefficient.
+   * @returns The Egm96Entry object containing the coefficients.
+   */
+  getEgm96Coeffs(l: number, m: number): Egm96Entry {
+    return egm96Data.getCoeffs(l, m);
+  }
+
+  /**
+   * Retrieves the IAU 1980 coefficients for the specified row.
+   * @param row The row index of the coefficients to retrieve.
+   * @returns The IAU 1980 entry containing the coefficients.
+   */
+  getIau1980Coeffs(row: number): Iau1980Entry {
+    return iau1980Data.getCoeffs(row);
+  }
+
+  /**
+   * Retrieves the number of leap seconds for a given Julian date.
+   * @param jd The Julian date.
+   * @returns The number of leap seconds.
+   */
+  getLeapSeconds(jd: number): number {
+    return leapSecondData.getLeapSeconds(jd);
+  }
+
+  /**
+   * Retrieves the atmosphere data for a given height.
+   * @param height The height for which to retrieve the atmosphere data.
+   * @returns The atmosphere data for the given height, or null if no data is available.
+   */
+  getHpAtmosphere(height: number): HpAtmosphereResult | null {
+    return hpAtmosphereData.getAtmosphere(height);
+  }
+}

package/src/data/values/Egm96Data.ts

@@ -0,0 +1,75 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+import { factorial } from '../../utils/functions.js';
+import { egm96 } from './egm96.js';
+
+// / EGM-96 entry for l, m indexes and clm, slm values.
+export type Egm96Entry = [number, number, number, number];
+
+// / Container for EGM-96 data.
+export class Egm96Data {
+  private readonly coeffs_: Egm96Entry[];
+
+  // / Create a new [Egm96Data] object.
+  constructor(coeffs: Egm96Entry[]) {
+    this.coeffs_ = coeffs;
+  }
+
+  /**
+   * Create a new [Egm96Data] container given a list of EGM-96
+   * coefficient tuples [vals].
+   * @param vals List of EGM-96 coefficient tuples.
+   * @returns A new [Egm96Data] object.
+   */
+  static fromVals(vals: Egm96Entry[]): Egm96Data {
+    const output: Egm96Entry[] = [];
+
+    for (const v of vals) {
+      const [l, m, clm, slm] = v;
+      const k = m === 0 ? 1 : 2;
+      const a = factorial(l + m);
+      const b = factorial(l - m) * (k * (2 * l + 1));
+      const nFac = Math.sqrt(a / b);
+      const normalizedClm = clm / nFac;
+      const normalizedSlm = slm / nFac;
+
+      output.push([l, m, normalizedClm, normalizedSlm]);
+    }
+
+    return new Egm96Data(output);
+  }
+
+  // / Return de-normalized EGM-96 coefficients for a given [l] and [m] index.
+  getCoeffs(l: number, m: number): Egm96Entry {
+    return this.coeffs_[Egm96Data.index_(l, m)];
+  }
+
+  // / Return the EGM-96 index for a given [l] and [m] lookup.
+  private static index_(l: number, m: number): number {
+    return (((l - 2) * (l + 2) + l) >> 1) - 1 + m;
+  }
+}
+
+// / De-normalized EGM-96 data container.
+export const egm96Data = Egm96Data.fromVals(egm96);

package/src/data/values/HpAtmosphereData.ts

@@ -0,0 +1,76 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+import { hpAtmosphere } from './hpAtmosphere.js';
+import { HpAtmosphereResult } from './HpAtmosphereResult.js';
+
+// / Harris-Priester atmosphere entry for height, min, and max density values.
+export type HpAtmosphereEntry = [number, number, number];
+
+// / Container for Harris-Priester atmosphere data.
+export class HpAtmosphereData {
+  private readonly table_: HpAtmosphereEntry[];
+  private readonly hMin_: number;
+  private readonly hMax_: number;
+
+  constructor(table: HpAtmosphereEntry[]) {
+    this.table_ = table;
+    if (table.length === 0 || typeof table[0]?.[0] === 'undefined') {
+      throw new Error('Table must have at least one valid entry.');
+    }
+
+    this.hMin_ = table[0][0];
+    this.hMax_ = table[table.length - 1]?.[0] ?? 0;
+  }
+
+  static fromVals(vals: [number, number, number][]): HpAtmosphereData {
+    const output: HpAtmosphereEntry[] = [];
+
+    for (const v of vals) {
+      const [h, minD, maxD] = v;
+
+      output.push([h, minD, maxD]);
+    }
+
+    return new HpAtmosphereData(output);
+  }
+
+  getAtmosphere(height: number): HpAtmosphereResult | null {
+    if (height < this.hMin_ || height > this.hMax_) {
+      return null;
+    }
+    let index = 0;
+
+    while (index < this.table_.length - 2 && height > (this.table_[index + 1])[0]) {
+      index++;
+    }
+
+    return new HpAtmosphereResult(
+      height,
+      (this.table_[index]),
+      (this.table_[index + 1]),
+    );
+  }
+}
+
+export const hpAtmosphereData = HpAtmosphereData.fromVals(hpAtmosphere);

package/src/data/values/HpAtmosphereResult.ts

@@ -0,0 +1,41 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+import { HpAtmosphereEntry } from './HpAtmosphereData.js';
+
+// / Harris-Priester atmospheric density bracket.
+export class HpAtmosphereResult {
+  // / Height above Earth's surface _(km)_.
+  height: number;
+  // / Lower bound for atmospheric parameters.
+  hp0: HpAtmosphereEntry;
+  // / Upper bound for atmospheric parameters.
+  hp1: HpAtmosphereEntry;
+
+  // / Create a new [HpAtmosphereResult] object.
+  constructor(height: number, hp0: HpAtmosphereEntry, hp1: HpAtmosphereEntry) {
+    this.height = height;
+    this.hp0 = hp0;
+    this.hp1 = hp1;
+  }
+}

package/src/data/values/Iau1980Data.ts

@@ -0,0 +1,63 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * 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/>.
+ */
+
+import { iau1980 } from './iau1980.js';
+
+// / IAU a1, a2, a3, a4, a5, Ai, Bi, Ci, Di coefficients.
+export type Iau1980Entry = [number, number, number, number, number, number, number, number, number];
+
+// / Container for IAU-1980 data.
+export class Iau1980Data {
+  private readonly coeffs_: Iau1980Entry[];
+
+  // / Create a new [Iau1980Data] object.
+  constructor(coeffs: Iau1980Entry[]) {
+    this.coeffs_ = coeffs;
+  }
+
+  /**
+   * Create a new [Iau1980Data] container object from an array of IAU-1980
+   * coefficient tuples [coeffs].
+   * @param coeffs IAU-1980 coefficients.
+   * @returns A new [Iau1980Data] object.
+   */
+  static fromCoeffs(coeffs: Array<Iau1980Entry>): Iau1980Data {
+    const output: Iau1980Entry[] = [];
+
+    for (const c of coeffs) {
+      const [a1, a2, a3, a4, a5, ai, bi, ci, di] = c;
+
+      output.push([a1, a2, a3, a4, a5, ai, bi, ci, di]);
+    }
+
+    return new Iau1980Data(output);
+  }
+
+  // / Get IAU-1980 coefficients for a given row number.
+  getCoeffs(row: number): Iau1980Entry {
+    return this.coeffs_[row];
+  }
+}
+
+// / IAU-1980 data container.
+export const iau1980Data = Iau1980Data.fromCoeffs(iau1980);