@saber-usa/node-common 1.7.1 → 1.7.2

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@saber-usa/node-common",
-  "version": "1.7.1",
+  "version": "1.7.2",
   "description": "Common node functions for Saber",
   "main": "src/index.js",
   "type": "module",

package/src/astro.js

@@ -1539,6 +1539,53 @@ const keplerianToCartesian = (elset, mu = 398600.4418) => {
     }
 };
 
+const cartesianToElsetElements = (pv, epoch) => {
+    // sma, eccentricity, inclination, raan, argp, trueAnomaly
+    const kepl = cartesianToKeplerian(
+        multiply(posToArray(pv.position), 1000.0), // km to m
+        multiply(posToArray(pv.velocity), 1000.0), // km/s to m/s
+    );
+
+    const elset = {
+        SemiMajorAxis: kepl.a,
+        Eccentricity: kepl.e,
+        Inclination: kepl.i,
+        Raan: kepl.raan,
+        ArgOfPerigee: kepl.w,
+    };
+
+    // Mean motion in radians per second
+    const mu = 3.986004418e14; // (m^3)/(s^2) WGS-84 Earth Mu
+    const meanMotion = Math.sqrt(mu/Math.pow((elset.SemiMajorAxis*1000.0), 3));
+
+    elset.MeanMotion = meanMotion / (2*Math.PI) * 60 * 60 * 24; // rads/s to revs per day
+
+    const trueAnomaly = kepl.f;
+
+    // Compute Eccentric Anomaly from True Anomaly and Eccentricity
+    const E = 2 * Math.atan2(Math.tan(trueAnomaly * DEG2RAD)
+        * Math.sqrt(1-elset.Eccentricity), Math.sqrt(1+elset.Eccentricity));
+
+    // Compute Mean Anomaly from Eccentric Anomaly and Eccentricity
+    elset.MeanAnomaly = (E - (elset.Eccentricity)*Math.sin(E)) * RAD2DEG;
+
+    elset.Period = 2*Math.PI/meanMotion / 60; // period in minutes
+    elset.Apogee = elset.SemiMajorAxis * (1 + elset.Eccentricity); // km
+    elset.Perigee = elset.SemiMajorAxis * (1 - elset.Eccentricity); // km
+
+    elset.Longitude = wrapToRange(
+        eciToGeodetic(
+            pv.position,
+            gstime(new Date(parseDate(epoch))),
+        ).longitude*RAD2DEG,
+        0, 360);
+
+    elset.LonDriftDegreesPerDay = null;
+    elset.BStar = null;
+
+    return elset;
+};
+
 /**
  * Get LEO RPO data for a given target satellite and a set of potential threat satellites.
  * @param {String} line1, line 1 of the target satellite
@@ -3131,6 +3178,7 @@ export {REGIMES,
     planeChangePureInclinationDeltaV,
     cartesianToKeplerian,
     keplerianToCartesian,
+    cartesianToElsetElements,
     getLeoRpoData,
     getGeoRpoData,
     getGeoShadowZones,