gnss-js 1.0.1 → 1.2.0

package/dist/index.js

@@ -1,9 +1,14 @@
+import {
+  connectToMountpoint,
+  fetchSourcetable,
+  parseSourcetable
+} from "./chunk-BJHTBYNG.js";
 import {
   CompletenessAccumulator,
   CycleSlipAccumulator,
   MultipathAccumulator,
   analyzeQuality
-} from "./chunk-SO3POWWR.js";
+} from "./chunk-G3N4S3DM.js";
 import {
   frequencyLabel,
   parseAntex
@@ -48,43 +53,71 @@ import {
   phiBOCc,
   phiBOCs,
   phiBPSK
-} from "./chunk-HKN3PUGN.js";
+} from "./chunk-5SPJH4MG.js";
 import {
   Scale,
-  getBdsTime,
-  getDateFromBdsTime,
   getDateFromDayOfWeek,
   getDateFromDayOfYear,
-  getDateFromGalTime,
-  getDateFromGloN,
-  getDateFromGpsData,
-  getDateFromGpsTime,
   getDateFromHourCode,
   getDateFromJulianDate,
   getDateFromMJD,
   getDateFromMJD2000,
   getDateFromRINEX,
-  getDateFromTai,
   getDateFromTimeOfDay,
+  getDayOfWeek,
+  getDayOfYear,
+  getHourCode,
+  getJulianDate,
+  getMJD,
+  getMJD2000,
+  getRINEX,
+  getTimeOfDay,
+  getWeekOfYear
+} from "./chunk-L5OUUNLT.js";
+import {
+  deg2dms,
+  deg2rad,
+  euclidean3D,
+  geodeticToGeohash,
+  geodeticToMaidenhead,
+  geodeticToUtm,
+  greatCircleMidpoint,
+  horizonDistance,
+  rad2deg,
+  rhumbLine,
+  vincenty
+} from "./chunk-MIIM4LDY.js";
+import {
+  computeAllPositions,
+  computeDop,
+  computeLiveSkyPositions,
+  computeSatPosition,
+  ecefToAzEl,
+  glonassPosition,
+  keplerPosition,
+  navTimesFromEph,
+  selectEphemeris
+} from "./chunk-VANVV4F7.js";
+import {
+  getBdsTime,
+  getDateFromBdsTime,
+  getDateFromGalTime,
+  getDateFromGloN,
+  getDateFromGpsData,
+  getDateFromGpsTime,
+  getDateFromTai,
   getDateFromTt,
   getDateFromUnixTime,
   getDateFromUtc,
-  getDayOfWeek,
-  getDayOfYear,
   getGalTime,
   getGloN4,
   getGloNA,
   getGpsLeap,
   getGpsTime,
-  getHourCode,
   getHoursFromTimeDifference,
-  getJulianDate,
   getLeap,
-  getMJD,
-  getMJD2000,
   getMinutesFromTimeDifference,
   getNtpTime,
-  getRINEX,
   getSecondsFromTimeDifference,
   getTaiDate,
   getTimeDifference,
@@ -93,54 +126,27 @@ import {
   getTimeDifferenceFromMinutes,
   getTimeDifferenceFromObject,
   getTimeDifferenceFromSeconds,
-  getTimeOfDay,
   getTimeOfWeek,
   getTotalDaysFromTimeDifference,
   getTtDate,
   getUnixTime,
   getUtcDate,
-  getWeekNumber,
-  getWeekOfYear
-} from "./chunk-7NOFXKET.js";
+  getWeekNumber
+} from "./chunk-HVXYFUCB.js";
 import {
-  deg2dms,
-  deg2rad,
-  euclidean3D,
-  geodeticToGeohash,
-  geodeticToMaidenhead,
-  geodeticToUtm,
-  greatCircleMidpoint,
-  horizonDistance,
-  rad2deg,
-  rhumbLine,
-  vincenty
-} from "./chunk-MIIM4LDY.js";
+  clampUnit,
+  ecefToGeodetic,
+  geodeticToEcef,
+  getAer,
+  getEnuDifference
+} from "./chunk-37QNKGTC.js";
 import "./chunk-5S5MONFK.js";
 import {
-  DAYS_MJD2000_MJD,
-  DAYS_MJD_JULIAN,
-  MILLISECONDS_GPS_TAI,
-  MILLISECONDS_IN_DAY,
-  MILLISECONDS_IN_HOUR,
-  MILLISECONDS_IN_MINUTE,
-  MILLISECONDS_IN_SECOND,
-  MILLISECONDS_IN_WEEK,
-  MILLISECONDS_TT_TAI,
-  RINEX_CODES,
-  SECONDS_IN_DAY,
-  SECONDS_IN_HOUR,
-  SECONDS_IN_MINUTE,
-  SECONDS_IN_WEEK,
-  START_BDS_TIME,
-  START_GAL_TIME,
-  START_GLO_LEAP,
-  START_GPS_TIME,
-  START_JULIAN_TAI,
-  START_MJD_UNIX_SECONDS,
-  START_NTP_TIME,
-  START_TAI_TIME,
-  START_UNIX_TIME
-} from "./chunk-LEEU5OIO.js";
+  WGS84_ECCENTRICITY_SQUARED,
+  WGS84_FLATTENING,
+  WGS84_SEMI_MAJOR_AXIS,
+  WGS84_SEMI_MINOR_AXIS
+} from "./chunk-6FAL6P4G.js";
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
@@ -159,7 +165,7 @@ import {
   parseRinexStream,
   systemCmp,
   systemName
-} from "./chunk-HBLU2EJ4.js";
+} from "./chunk-OZCYOM5D.js";
 import {
   BitReader,
   RTCM3_MESSAGE_NAMES,
@@ -170,11 +176,39 @@ import {
   decodeMsmFull,
   msmEpochToDate,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rtcm3Constellation,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-WP2JFDLA.js";
+} from "./chunk-7EEWQ5DU.js";
+import {
+  DAYS_MJD2000_MJD,
+  DAYS_MJD_JULIAN,
+  MILLISECONDS_GPS_TAI,
+  MILLISECONDS_IN_DAY,
+  MILLISECONDS_IN_HOUR,
+  MILLISECONDS_IN_MINUTE,
+  MILLISECONDS_IN_SECOND,
+  MILLISECONDS_IN_WEEK,
+  MILLISECONDS_TT_TAI,
+  RINEX_CODES,
+  SECONDS_IN_DAY,
+  SECONDS_IN_HOUR,
+  SECONDS_IN_MINUTE,
+  SECONDS_IN_WEEK,
+  START_BDS_TIME,
+  START_GAL_TIME,
+  START_GLO_LEAP,
+  START_GPS_TIME,
+  START_JULIAN_TAI,
+  START_MJD_UNIX_SECONDS,
+  START_NTP_TIME,
+  START_TAI_TIME,
+  START_UNIX_TIME
+} from "./chunk-LEEU5OIO.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -196,36 +230,7 @@ import {
   formatUTCTime,
   getFreq,
   gloFreq
-} from "./chunk-LWNTWBHB.js";
-import {
-  computeAllPositions,
-  computeDop,
-  computeLiveSkyPositions,
-  computeSatPosition,
-  ecefToAzEl,
-  glonassPosition,
-  keplerPosition,
-  navTimesFromEph,
-  selectEphemeris
-} from "./chunk-4YN353Q7.js";
-import {
-  clampUnit,
-  ecefToGeodetic,
-  geodeticToEcef,
-  getAer,
-  getEnuDifference
-} from "./chunk-37QNKGTC.js";
-import {
-  WGS84_ECCENTRICITY_SQUARED,
-  WGS84_FLATTENING,
-  WGS84_SEMI_MAJOR_AXIS,
-  WGS84_SEMI_MINOR_AXIS
-} from "./chunk-6FAL6P4G.js";
-import {
-  connectToMountpoint,
-  fetchSourcetable,
-  parseSourcetable
-} from "./chunk-BJHTBYNG.js";
+} from "./chunk-W5WKEV7U.js";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -413,10 +418,13 @@ export {
   phiBPSK,
   rad2deg,
   readString,
+  reportDecodeError,
   resetGloFreqCache,
   rhumbLine,
   rtcm3Constellation,
   selectEphemeris,
+  setGloFreqNumber,
+  setRtcm3DebugHandler,
   systemCmp,
   systemName,
   updateStationMeta,

package/dist/orbit.cjs

@@ -71,12 +71,103 @@ function ecefToGeodetic(x, y, z) {
   return [lat, lon, alt];
 }
 
+// src/constants/time.ts
+var MILLISECONDS_IN_SECOND = 1e3;
+var MILLISECONDS_GPS_TAI = 19e3;
+var START_GPS_TIME = /* @__PURE__ */ new Date("1980-01-06T00:00:00Z");
+var START_BDS_TIME = /* @__PURE__ */ new Date("2006-01-01T00:00:14Z");
+var START_NTP_TIME = /* @__PURE__ */ new Date("1900-01-01T00:00:00Z");
+
+// src/time/gnss.ts
+function getTaiDate(date) {
+  return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
+}
+function getNtpTime(date) {
+  return getTaiDate(date).getTime() - START_NTP_TIME.getTime();
+}
+
+// src/time/utc.ts
+var LEAP_SECONDS_TABLE = [
+  [3692217600, 37],
+  // 1 Jan 2017
+  [3644697600, 36],
+  // 1 Jul 2015
+  [3550089600, 35],
+  // 1 Jul 2012
+  [3439756800, 34],
+  // 1 Jan 2009
+  [3345062400, 33],
+  // 1 Jan 2006
+  [3124137600, 32],
+  // 1 Jan 1999
+  [3076704e3, 31],
+  // 1 Jul 1997
+  [3029443200, 30],
+  // 1 Jan 1996
+  [2982009600, 29],
+  // 1 Jul 1994
+  [2950473600, 28],
+  // 1 Jul 1993
+  [2918937600, 27],
+  // 1 Jul 1992
+  [2871676800, 26],
+  // 1 Jan 1991
+  [2840140800, 25],
+  // 1 Jan 1990
+  [2776982400, 24],
+  // 1 Jan 1988
+  [2698012800, 23],
+  // 1 Jul 1985
+  [2634854400, 22],
+  // 1 Jul 1983
+  [2603318400, 21],
+  // 1 Jul 1982
+  [2571782400, 20],
+  // 1 Jul 1981
+  [2524521600, 19],
+  // 1 Jan 1980
+  [2492985600, 18],
+  // 1 Jan 1979
+  [2461449600, 17],
+  // 1 Jan 1978
+  [2429913600, 16],
+  // 1 Jan 1977
+  [2398291200, 15],
+  // 1 Jan 1976
+  [2366755200, 14],
+  // 1 Jan 1975
+  [2335219200, 13],
+  // 1 Jan 1974
+  [2303683200, 12],
+  // 1 Jan 1973
+  [2287785600, 11],
+  // 1 Jul 1972
+  [2272060800, 10]
+  // 1 Jan 1972
+];
+function getLeap(date) {
+  const ntp_time = getNtpTime(date);
+  for (const [timestamp, leapSeconds] of LEAP_SECONDS_TABLE) {
+    if (ntp_time / MILLISECONDS_IN_SECOND - leapSeconds >= timestamp) {
+      return leapSeconds;
+    }
+  }
+  return 8;
+}
+function getGpsLeap(date) {
+  const leap_seconds = getLeap(date);
+  if (leap_seconds < 0) return 0;
+  return leap_seconds - 19;
+}
+
 // src/orbit/index.ts
 var GM_GPS = 3986005e8;
 var GM_GAL = 3986004418e5;
 var GM_BDS = 3986004418e5;
 var GM_GLO = 39860044e7;
 var OMEGA_E = 72921151467e-15;
+var BDS_GEO_MAX_INCLINATION_RAD = 0.1;
+var BDS_GEO_MIN_SEMIMAJOR_AXIS_M = 4e7;
 var AE_GLO = 6378136;
 var J2_GLO = 108263e-8;
 var TWO_PI = 2 * Math.PI;
@@ -114,7 +205,7 @@ function keplerPosition(eph, t) {
   const ik = eph.i0 + dik + eph.idot * tk;
   const xp = rk * Math.cos(uk);
   const yp = rk * Math.sin(uk);
-  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < 0.1 && a > 4e7;
+  const isBdsGeo = eph.system === "C" && Math.abs(eph.i0) < BDS_GEO_MAX_INCLINATION_RAD && a > BDS_GEO_MIN_SEMIMAJOR_AXIS_M;
   if (isBdsGeo) {
     const omegak2 = eph.omega0 + eph.omegaDot * tk - OMEGA_E * eph.toe;
     const cosO2 = Math.cos(omegak2);
@@ -293,11 +384,13 @@ function selectEphemeris(ephemerides, prn, timeMs) {
 }
 function computeSatPosition(eph, timeMs) {
   if (eph.system === "R" || eph.system === "S") {
-    return glonassPosition(eph, timeMs / 1e3);
+    const leapMs = getGpsLeap(new Date(timeMs)) * 1e3;
+    return glonassPosition(eph, (timeMs - leapMs) / 1e3);
   }
-  const GPS_EPOCH = Date.UTC(1980, 0, 6);
+  const GPS_EPOCH = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH) / 1e3;
-  const tow = gpsSeconds % (7 * 86400);
+  let tow = gpsSeconds % (7 * 86400);
+  if (eph.system === "C") tow = (tow - 14 + 604800) % 604800;
   return keplerPosition(eph, tow);
 }
 function navTimesFromEph(ephs) {
@@ -379,8 +472,8 @@ function selectBest(ephs, timeMs) {
   }
   return best;
 }
-var GPS_EPOCH_MS = Date.UTC(1980, 0, 6);
-var BDS_EPOCH_MS = Date.UTC(2006, 0, 1);
+var GPS_EPOCH_MS = START_GPS_TIME.getTime();
+var BDS_EPOCH_MS = START_BDS_TIME.getTime();
 var GAL_EPOCH_MS = GPS_EPOCH_MS;
 function ephInfoToEphemeris(info) {
   const sys = info.prn.charAt(0);
@@ -421,6 +514,9 @@ function ephInfoToEphemeris(info) {
   }
   if (info.sqrtA === void 0 || info.eccentricity === void 0 || info.inclination === void 0 || info.omega0 === void 0 || info.argPerigee === void 0 || info.meanAnomaly === void 0 || info.toe === void 0 || info.week === void 0)
     return null;
+  if (info.eccentricity < 0 || info.eccentricity >= 1 || info.sqrtA < 3e3 || info.sqrtA > 9e3 || info.toe < 0 || info.toe > 604800) {
+    return null;
+  }
   let epochMs;
   const tocSec = info.toc ?? info.toe;
   if (sys === "C") {
@@ -462,7 +558,7 @@ function ephInfoToEphemeris(info) {
   };
 }
 function computeLiveSkyPositions(ephemerides, rxPos, cn0Map) {
-  const now = Date.now();
+  const now = Date.now() + getGpsLeap(/* @__PURE__ */ new Date()) * 1e3;
   const result = [];
   const [rxX, rxY, rxZ] = rxPos;
   const [rxLat, rxLon] = ecefToGeodetic(rxX, rxY, rxZ);

package/dist/orbit.d.cts

@@ -1,5 +1,5 @@
 import { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris } from './nav-BAI1a9vK.cjs';
-import { E as EphemerisInfo } from './ephemeris-C10stHhM.cjs';
+import { E as EphemerisInfo } from './ephemeris-Ohl6NAfw.cjs';
 export { e as ecefToGeodetic, g as geodeticToEcef } from './ecef-CF0uAysr.cjs';
 
 /**

package/dist/orbit.d.ts

@@ -1,5 +1,5 @@
 import { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris } from './nav-BAI1a9vK.js';
-import { E as EphemerisInfo } from './ephemeris-C10stHhM.js';
+import { E as EphemerisInfo } from './ephemeris-Ohl6NAfw.js';
 export { e as ecefToGeodetic, g as geodeticToEcef } from './ecef-CF0uAysr.js';
 
 /**

package/dist/orbit.js

@@ -8,12 +8,14 @@ import {
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-4YN353Q7.js";
+} from "./chunk-VANVV4F7.js";
+import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
   geodeticToEcef
 } from "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
+import "./chunk-LEEU5OIO.js";
 export {
   computeAllPositions,
   computeDop,