ootk 4.0.9 → 5.1.0

package/src/utils/constants.ts

@@ -0,0 +1,113 @@
+import { Days, Degrees, Milliseconds, Minutes, Radians, Seconds } from '../types/types.js';
+
+/**
+ * Full circle in radians (PI * 2)
+ *
+ * https://tauday.com/tau-manifesto
+ */
+export const TAU = (2.0 * Math.PI) as Radians;
+
+/**
+ * Represents half of the mathematical constant PI.
+ */
+export const halfPi = (0.5 * Math.PI) as Radians;
+
+/**
+ * Converts degrees to radians.
+ */
+export const DEG2RAD = (Math.PI / 180.0) as Radians;
+
+/**
+ * Converts radians to degrees.
+ */
+export const RAD2DEG = (180.0 / Math.PI) as Degrees;
+
+/**
+ * Conversion factor from seconds to degrees.
+ */
+export const sec2deg = (1.0 / 60.0 / 60.0) as Degrees;
+
+/**
+ * Conversion factor from seconds to days.
+ */
+export const sec2day: number = sec2deg / 24.0;
+
+/**
+ * Conversion factor from arcseconds to radians.
+ */
+export const asec2rad = (sec2deg * DEG2RAD) as Radians;
+
+/**
+ * Convert ten-thousandths of an arcsecond to radians.
+ */
+export const ttasec2rad = (asec2rad / 10000.0) as Radians;
+
+/**
+ * Convert milliarcseconds to radians.
+ */
+export const masec2rad = (asec2rad / 1000.0) as Radians;
+
+/**
+ * The angular velocity of the Earth in radians per second.
+ */
+export const angularVelocityOfEarth = 7.292115e-5;
+
+/**
+ * Astronomical unit in kilometers.
+ */
+export const astronomicalUnit = 149597870.0;
+
+// / Convert milliseconds to seconds.
+export const msec2sec = 1e-3 as Seconds;
+
+// / Speed of light.
+export const cMPerSec = 299792458;
+export const cKmPerSec = 299792458 / 1000;
+export const cKmPerMs = 299792458 / 1000 / 1000;
+
+// / Milliseconds per day.
+export const MS_PER_DAY = 86400000;
+
+// / Seconds per day.
+export const secondsPerDay = 86400.0;
+
+// / Convert seconds to minutes.
+export const sec2min = (1.0 / 60.0) as Minutes;
+
+// / Seconds per sidereal day.
+export const secondsPerSiderealDay = 86164.0905;
+
+// / Seconds per week.
+export const secondsPerWeek: number = secondsPerDay * 7.0;
+
+/**
+ * Half the number of radians in a circle.
+ */
+export const PI = Math.PI as Radians;
+
+export const x2o3 = 2.0 / 3.0;
+export const temp4 = 1.5e-12;
+
+/**
+ * The number of minutes in a day.
+ */
+export const MINUTES_PER_DAY = 1440 as Minutes;
+
+/**
+ * The number of milliseconds in a day.
+ */
+export const MILLISECONDS_TO_DAYS = 1.15741e-8;
+
+/**
+ * The number of milliseconds in a day.
+ */
+export const MILLISECONDS_PER_DAY = <Days>1000 * 60 * 60 * 24;
+
+/**
+ * The number of milliseconds in a second.
+ */
+export const MILLISECONDS_PER_SECOND = <Milliseconds>1000;
+
+export const RADIUS_OF_EARTH = 6371; // Radius of Earth in kilometers
+
+export const earthGravityParam = 398600.4415;

package/src/utils/create-covariance-from-tle.ts

@@ -0,0 +1,75 @@
+import { Tle, CovarianceSample, CovarianceFrame, StateCovariance } from '../main.js';
+
+/**
+ * Creates a 6x6 state covariance matrix from a TLE
+ * @param tleLine1 The first line of the TLE
+ * @param tleLine2 The second line of the TLE
+ * @param frame The covariance frame (CovarianceFrame.ECI or CovarianceFrame.RIC)
+ * @param sigmaScale Scaling factor for the sigmas (default: 1.0)
+ * @returns A StateCovariance object containing the 6x6 covariance matrix
+ */
+export function createCovarianceFromTle(
+  tleLine1: string,
+  tleLine2: string,
+  frame: CovarianceFrame = CovarianceFrame.RIC,
+  sigmaScale: number = 1.0,
+): StateCovariance {
+  // Parse the TLE and get the state vector
+  const tle = new Tle(tleLine1, tleLine2);
+  const temeState = tle.state;
+  const j2000State = temeState.toJ2000();
+
+  // eslint-disable-next-line no-console
+  console.log('J2000 State:', j2000State);
+
+  /*
+   * Determine appropriate sigma values based on TLE
+   * These are rough estimates and can be adjusted based on your needs
+   */
+  const positionSigma = 1.0 * sigmaScale; // km
+  const velocitySigma = 0.001 * sigmaScale; // km/s
+
+  // Create a diagonal covariance matrix with sigma values
+  const sigmas = [
+    positionSigma, positionSigma, positionSigma,
+    velocitySigma, velocitySigma, velocitySigma,
+  ];
+
+  // Generate the covariance matrix in the desired frame
+  return StateCovariance.fromSigmas(sigmas, frame);
+}
+
+/**
+ * Creates a sample-based covariance from a TLE with more realistic uncertainties
+ * @param tleLine1 The first line of the TLE
+ * @param tleLine2 The second line of the TLE
+ * @param frame The covariance frame (CovarianceFrame.ECI or CovarianceFrame.RIC)
+ * @returns A StateCovariance object
+ */
+export function createSampleCovarianceFromTle(
+  tleLine1: string,
+  tleLine2: string,
+  frame: CovarianceFrame = CovarianceFrame.RIC,
+): StateCovariance {
+  // Parse the TLE and get the state vector
+  const tle = new Tle(tleLine1, tleLine2);
+  const temeState = tle.state;
+  const j2000State = temeState.toJ2000();
+
+  /*
+   * Create initial covariance with basic sigma values
+   * Position uncertainties are higher in-track than radial/cross-track
+   * for most space catalog objects
+   */
+  const sigmas = frame === CovarianceFrame.RIC
+    ? [0.12, 1.0, 0.1, 0.00012, 0.001, 0.0001] // RIC frame: [R,I,C,Rdot,Idot,Cdot]
+    : [0.6, 0.6, 0.6, 0.0006, 0.0006, 0.0006]; // ECI frame: [x,y,z,vx,vy,vz]
+
+  const covariance = StateCovariance.fromSigmas(sigmas, frame);
+
+  // Create a covariance sample that will be used to generate a more realistic covariance
+  const sample = new CovarianceSample(j2000State, covariance, tle);
+
+  // Return the desample in the appropriate frame
+  return frame === CovarianceFrame.RIC ? sample.desampleRIC() : sample.desampleJ2000();
+}