gnss-js 1.1.0 → 1.3.0

package/README.md

@@ -51,6 +51,44 @@ const { distance } = vincenty(
 // distance in meters, angles in radians
 ```
 
+### Visibility / pass prediction
+
+```ts
+import { computeVisibility } from 'gnss-js/orbit';
+
+const vis = computeVisibility(ephemerides, rxEcef, startMs, endMs, 300, 10);
+// vis.passes: rise/peak/set per satellite; vis.pdop / vis.visibleCount
+// per epoch; vis.elevation[prn]: elevation timeline for plotting.
+```
+
+### Single-point positioning
+
+```ts
+import { solveSpp, ionoFree } from 'gnss-js/positioning';
+
+const sol = solveSpp(pseudoranges, ephemerides, epochMs);
+// sol.position (ECEF m), per-system clock biases, DOP, residuals,
+// rejected outliers. ~10 m single-frequency; use ionoFree(p1, p2, f1, f2)
+// pseudoranges for metre-level.
+```
+
+### Reference frame transformations
+
+```ts
+import { transformFrame, dateToEpoch } from 'gnss-js/frames';
+
+// ETRS89 (ETRF2000) station coordinates → ITRF2020 at epoch 2010.5
+const itrf = transformFrame(
+  [3924687.7039, 301132.7618, 5001910.7712],
+  'ETRF2000',
+  'ITRF2020',
+  2010.5
+);
+// Time-dependent 14-parameter Helmert (IERS/EPSG parameters);
+// supports ITRF88–ITRF2020, ETRF2000/2014/2020, NAD83(2011).
+// 'WGS84' is treated as an ITRF2020 alias (standard GNSS practice).
+```
+
 ### RINEX parsing
 
 ```ts

package/dist/{chunk-7NOFXKET.js → chunk-HVXYFUCB.js}

@@ -1,6 +1,4 @@
 import {
-  DAYS_MJD2000_MJD,
-  DAYS_MJD_JULIAN,
   MILLISECONDS_GPS_TAI,
   MILLISECONDS_IN_DAY,
   MILLISECONDS_IN_HOUR,
@@ -8,98 +6,13 @@ import {
   MILLISECONDS_IN_SECOND,
   MILLISECONDS_IN_WEEK,
   MILLISECONDS_TT_TAI,
-  RINEX_CODES,
   START_BDS_TIME,
   START_GAL_TIME,
   START_GLO_LEAP,
   START_GPS_TIME,
-  START_JULIAN_TAI,
-  START_NTP_TIME,
-  START_TAI_TIME
+  START_NTP_TIME
 } from "./chunk-LEEU5OIO.js";
 
-// src/types/enums.ts
-var Scale = /* @__PURE__ */ ((Scale2) => {
-  Scale2[Scale2["TAI"] = 0] = "TAI";
-  Scale2[Scale2["TT"] = 1] = "TT";
-  Scale2[Scale2["GPS"] = 2] = "GPS";
-  Scale2[Scale2["UTC"] = 3] = "UTC";
-  return Scale2;
-})(Scale || {});
-
-// src/time/functions.ts
-function getTimeOfDay(date) {
-  return date.getTime() - Date.UTC(
-    date.getUTCFullYear(),
-    date.getUTCMonth(),
-    date.getUTCDate(),
-    0,
-    0,
-    0
-  );
-}
-function getDateFromTimeOfDay(timeOfDay, dateRaw) {
-  return new Date(
-    Date.UTC(
-      dateRaw.getUTCFullYear(),
-      dateRaw.getUTCMonth(),
-      dateRaw.getUTCDate(),
-      0,
-      0,
-      0
-    ) + timeOfDay
-  );
-}
-function getDayOfYear(date) {
-  return Math.floor(
-    (date.getTime() - Date.UTC(date.getUTCFullYear())) / MILLISECONDS_IN_DAY
-  ) + 1;
-}
-function getDateFromDayOfYear(dayOfYear, date) {
-  return new Date(
-    Date.UTC(
-      date.getUTCFullYear(),
-      0,
-      1,
-      date.getUTCHours(),
-      date.getUTCMinutes(),
-      date.getUTCSeconds(),
-      date.getUTCMilliseconds()
-    ) + (dayOfYear - 1) * MILLISECONDS_IN_DAY
-  );
-}
-function getDayOfWeek(date) {
-  return date.getUTCDay();
-}
-function getDateFromDayOfWeek(dayOfWeek, dateRaw) {
-  if (dayOfWeek < 0 || dayOfWeek > 6)
-    throw new Error("Day of week must be a value between 0 and 7");
-  return new Date(
-    dateRaw.getTime() + (dayOfWeek - dateRaw.getUTCDay()) * MILLISECONDS_IN_DAY
-  );
-}
-function getDateFromHourCode(hourCode, date) {
-  const hour = RINEX_CODES.indexOf(hourCode);
-  if (hour === -1)
-    throw new Error("Hour code must be a lowercase letter between a and x");
-  const newDate = new Date(date.getTime());
-  newDate.setUTCHours(hour);
-  return newDate;
-}
-function getHourCode(date) {
-  return RINEX_CODES[date.getUTCHours()];
-}
-function getWeekOfYear(date) {
-  const target = new Date(
-    Date.UTC(date.getFullYear(), date.getMonth(), date.getDate())
-  );
-  const dayNr = (target.getUTCDay() + 6) % 7;
-  target.setUTCDate(target.getUTCDate() - dayNr + 3);
-  const firstThursday = new Date(Date.UTC(target.getUTCFullYear(), 0, 4));
-  const dayDiff = (target.getTime() - firstThursday.getTime()) / (86400 * 1e3);
-  return 1 + Math.floor(dayDiff / 7) + 1;
-}
-
 // src/time/time-difference.ts
 function getTimeDifference(startDate, finalDate) {
   return Math.floor(finalDate.getTime() - startDate.getTime());
@@ -139,6 +52,76 @@ function getTimeDifferenceFromObject(timeDifferenceObject) {
   return getTimeDifferenceFromSeconds(seconds) + getTimeDifferenceFromMinutes(minutes) + getTimeDifferenceFromHours(hours) + getTimeDifferenceFromDays(days);
 }
 
+// src/time/gnss.ts
+function getGpsTime(date) {
+  return date.getTime() - START_GPS_TIME.getTime();
+}
+function getDateFromGpsTime(gpsTime) {
+  return new Date(gpsTime + START_GPS_TIME.getTime());
+}
+function getGalTime(date) {
+  return date.getTime() - START_GAL_TIME.getTime();
+}
+function getDateFromGalTime(galTime) {
+  return new Date(galTime + START_GAL_TIME.getTime());
+}
+function getBdsTime(date) {
+  return date.getTime() - START_BDS_TIME.getTime();
+}
+function getDateFromBdsTime(bdsTime) {
+  return new Date(bdsTime + START_BDS_TIME.getTime());
+}
+function getUnixTime(date) {
+  const utc_date = getUtcDate(date);
+  return utc_date.getTime();
+}
+function getDateFromUnixTime(unixTime) {
+  return getDateFromUtc(new Date(unixTime));
+}
+function getWeekNumber(date) {
+  return Math.floor(getGpsTime(date) / MILLISECONDS_IN_WEEK);
+}
+function getTimeOfWeek(date) {
+  return getGpsTime(date) % MILLISECONDS_IN_WEEK;
+}
+function getDateFromGpsData(weekNumber, timeOfWeek) {
+  return new Date(
+    weekNumber * MILLISECONDS_IN_WEEK + timeOfWeek + START_GPS_TIME.getTime()
+  );
+}
+function getTaiDate(date) {
+  return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
+}
+function getDateFromTai(tai_date) {
+  return new Date(tai_date.getTime() - MILLISECONDS_GPS_TAI);
+}
+function getTtDate(date) {
+  return new Date(date.getTime() + MILLISECONDS_GPS_TAI + MILLISECONDS_TT_TAI);
+}
+function getDateFromTt(tt_date) {
+  return new Date(
+    tt_date.getTime() - MILLISECONDS_GPS_TAI - MILLISECONDS_TT_TAI
+  );
+}
+function getGloN4(date) {
+  return Math.floor((date.getUTCFullYear() - START_GLO_LEAP.getFullYear()) / 4);
+}
+function getGloNA(date) {
+  const n4 = getGloN4(date);
+  const init4YearPeriod = new Date(START_GLO_LEAP);
+  init4YearPeriod.setUTCFullYear(init4YearPeriod.getUTCFullYear() + n4 * 4);
+  return Math.floor(getTimeDifference(init4YearPeriod, date) / MILLISECONDS_IN_DAY) + 1;
+}
+function getDateFromGloN(n4, na, tod) {
+  const date = new Date(START_GLO_LEAP);
+  date.setFullYear(date.getUTCFullYear() + n4 * 4);
+  date.setTime(date.getTime() + (na - 1) * MILLISECONDS_IN_DAY + tod);
+  return date;
+}
+function getNtpTime(date) {
+  return getTaiDate(date).getTime() - START_NTP_TIME.getTime();
+}
+
 // src/time/utc.ts
 var LEAP_SECONDS_TABLE = [
   [3692217600, 37],
@@ -228,167 +211,7 @@ function getDateFromUtc(utc_date) {
   return gps_date;
 }
 
-// src/time/gnss.ts
-function getGpsTime(date) {
-  return date.getTime() - START_GPS_TIME.getTime();
-}
-function getDateFromGpsTime(gpsTime) {
-  return new Date(gpsTime + START_GPS_TIME.getTime());
-}
-function getGalTime(date) {
-  return date.getTime() - START_GAL_TIME.getTime();
-}
-function getDateFromGalTime(galTime) {
-  return new Date(galTime + START_GAL_TIME.getTime());
-}
-function getBdsTime(date) {
-  return date.getTime() - START_BDS_TIME.getTime();
-}
-function getDateFromBdsTime(bdsTime) {
-  return new Date(bdsTime + START_BDS_TIME.getTime());
-}
-function getUnixTime(date) {
-  const utc_date = getUtcDate(date);
-  return utc_date.getTime();
-}
-function getDateFromUnixTime(unixTime) {
-  return getDateFromUtc(new Date(unixTime));
-}
-function getWeekNumber(date) {
-  return Math.floor(getGpsTime(date) / MILLISECONDS_IN_WEEK);
-}
-function getTimeOfWeek(date) {
-  return getGpsTime(date) % MILLISECONDS_IN_WEEK;
-}
-function getDateFromGpsData(weekNumber, timeOfWeek) {
-  return new Date(
-    weekNumber * MILLISECONDS_IN_WEEK + timeOfWeek + START_GPS_TIME.getTime()
-  );
-}
-function getTaiDate(date) {
-  return new Date(date.getTime() + MILLISECONDS_GPS_TAI);
-}
-function getDateFromTai(tai_date) {
-  return new Date(tai_date.getTime() - MILLISECONDS_GPS_TAI);
-}
-function getTtDate(date) {
-  return new Date(date.getTime() + MILLISECONDS_GPS_TAI + MILLISECONDS_TT_TAI);
-}
-function getDateFromTt(tt_date) {
-  return new Date(
-    tt_date.getTime() - MILLISECONDS_GPS_TAI - MILLISECONDS_TT_TAI
-  );
-}
-function getGloN4(date) {
-  return Math.floor((date.getUTCFullYear() - START_GLO_LEAP.getFullYear()) / 4);
-}
-function getGloNA(date) {
-  const n4 = getGloN4(date);
-  const init4YearPeriod = new Date(START_GLO_LEAP);
-  init4YearPeriod.setUTCFullYear(init4YearPeriod.getUTCFullYear() + n4 * 4);
-  return Math.floor(getTimeDifference(init4YearPeriod, date) / MILLISECONDS_IN_DAY) + 1;
-}
-function getDateFromGloN(n4, na, tod) {
-  const date = new Date(START_GLO_LEAP);
-  date.setFullYear(date.getUTCFullYear() + n4 * 4);
-  date.setTime(date.getTime() + (na - 1) * MILLISECONDS_IN_DAY + tod);
-  return date;
-}
-function getNtpTime(date) {
-  return getTaiDate(date).getTime() - START_NTP_TIME.getTime();
-}
-
-// src/time/rinex.ts
-function getRINEX(date) {
-  const seconds = date.getUTCSeconds();
-  const milliseconds = date.getUTCMilliseconds();
-  const fractionalSeconds = (milliseconds / MILLISECONDS_IN_SECOND).toFixed(7).slice(1);
-  const formatted = `> ${date.getUTCFullYear()} ${String(
-    date.getUTCMonth() + 1
-  ).padStart(2, "0")} ${String(date.getUTCDate()).padStart(2, "0")} ${String(
-    date.getUTCHours()
-  ).padStart(2, "0")} ${String(date.getUTCMinutes()).padStart(
-    2,
-    "0"
-  )} ${String(seconds).padStart(2, "0")}${fractionalSeconds}`;
-  return formatted;
-}
-function getDateFromRINEX(rinex) {
-  const year = Number.parseInt(rinex.substring(2, 6));
-  const month = Number.parseInt(rinex.substring(7, 9)) - 1;
-  const day = Number.parseInt(rinex.substring(10, 12));
-  const hour = Number.parseInt(rinex.substring(13, 15));
-  const minute = Number.parseInt(rinex.substring(16, 18));
-  const second = Number.parseInt(rinex.substring(19, 21));
-  const millisecond = Number.parseInt(rinex.substring(22, 29)) / 1e4;
-  const date = new Date(
-    Date.UTC(year, month, day, hour, minute, second, millisecond)
-  );
-  return date;
-}
-
-// src/time/julian.ts
-function getJulianDate(date, scale = 1 /* TT */) {
-  if (scale === 0 /* TAI */) {
-    return (date.getTime() - START_TAI_TIME.getTime()) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
-  }
-  if (scale === 2 /* GPS */) {
-    return (date.getTime() - START_TAI_TIME.getTime() - MILLISECONDS_GPS_TAI) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
-  }
-  if (scale === 3 /* UTC */) {
-    return (date.getTime() - START_TAI_TIME.getTime() - getLeap(date) * MILLISECONDS_IN_SECOND) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
-  }
-  return (date.getTime() - START_TAI_TIME.getTime() + MILLISECONDS_TT_TAI) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
-}
-function getDateFromJulianDate(julianDate, scale = 1 /* TT */) {
-  if (scale === 0 /* TAI */) {
-    return new Date(
-      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime()
-    );
-  }
-  if (scale === 2 /* GPS */) {
-    return new Date(
-      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime() + MILLISECONDS_GPS_TAI
-    );
-  }
-  if (scale === 3 /* UTC */) {
-    const utc_date = new Date(
-      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime()
-    );
-    const tai_date = new Date(
-      utc_date.getTime() + getLeap(utc_date) * MILLISECONDS_IN_SECOND
-    );
-    const leap_seconds = getLeap(getDateFromTai(tai_date));
-    return new Date(utc_date.getTime() + leap_seconds * MILLISECONDS_IN_SECOND);
-  }
-  return new Date(
-    (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime() - MILLISECONDS_TT_TAI
-  );
-}
-function getMJD(date, scale = 1 /* TT */) {
-  return getJulianDate(date, scale) - DAYS_MJD_JULIAN;
-}
-function getDateFromMJD(mjd, scale = 1 /* TT */) {
-  return getDateFromJulianDate(mjd + DAYS_MJD_JULIAN, scale);
-}
-function getMJD2000(date, scale = 1 /* TT */) {
-  return getMJD(date, scale) - DAYS_MJD2000_MJD;
-}
-function getDateFromMJD2000(mjd2000, scale = 1 /* TT */) {
-  return getDateFromMJD(mjd2000 + DAYS_MJD2000_MJD, scale);
-}
-
 export {
-  Scale,
-  getTimeOfDay,
-  getDateFromTimeOfDay,
-  getDayOfYear,
-  getDateFromDayOfYear,
-  getDayOfWeek,
-  getDateFromDayOfWeek,
-  getDateFromHourCode,
-  getHourCode,
-  getWeekOfYear,
   getTimeDifference,
   getSecondsFromTimeDifference,
   getMinutesFromTimeDifference,
@@ -421,13 +244,5 @@ export {
   getGloN4,
   getGloNA,
   getDateFromGloN,
-  getNtpTime,
-  getRINEX,
-  getDateFromRINEX,
-  getJulianDate,
-  getDateFromJulianDate,
-  getMJD,
-  getDateFromMJD,
-  getMJD2000,
-  getDateFromMJD2000
+  getNtpTime
 };

package/dist/chunk-L5OUUNLT.js

@@ -0,0 +1,198 @@
+import {
+  getDateFromTai,
+  getLeap
+} from "./chunk-HVXYFUCB.js";
+import {
+  DAYS_MJD2000_MJD,
+  DAYS_MJD_JULIAN,
+  MILLISECONDS_GPS_TAI,
+  MILLISECONDS_IN_DAY,
+  MILLISECONDS_IN_SECOND,
+  MILLISECONDS_TT_TAI,
+  RINEX_CODES,
+  START_JULIAN_TAI,
+  START_TAI_TIME
+} from "./chunk-LEEU5OIO.js";
+
+// src/types/enums.ts
+var Scale = /* @__PURE__ */ ((Scale2) => {
+  Scale2[Scale2["TAI"] = 0] = "TAI";
+  Scale2[Scale2["TT"] = 1] = "TT";
+  Scale2[Scale2["GPS"] = 2] = "GPS";
+  Scale2[Scale2["UTC"] = 3] = "UTC";
+  return Scale2;
+})(Scale || {});
+
+// src/time/functions.ts
+function getTimeOfDay(date) {
+  return date.getTime() - Date.UTC(
+    date.getUTCFullYear(),
+    date.getUTCMonth(),
+    date.getUTCDate(),
+    0,
+    0,
+    0
+  );
+}
+function getDateFromTimeOfDay(timeOfDay, dateRaw) {
+  return new Date(
+    Date.UTC(
+      dateRaw.getUTCFullYear(),
+      dateRaw.getUTCMonth(),
+      dateRaw.getUTCDate(),
+      0,
+      0,
+      0
+    ) + timeOfDay
+  );
+}
+function getDayOfYear(date) {
+  return Math.floor(
+    (date.getTime() - Date.UTC(date.getUTCFullYear())) / MILLISECONDS_IN_DAY
+  ) + 1;
+}
+function getDateFromDayOfYear(dayOfYear, date) {
+  return new Date(
+    Date.UTC(
+      date.getUTCFullYear(),
+      0,
+      1,
+      date.getUTCHours(),
+      date.getUTCMinutes(),
+      date.getUTCSeconds(),
+      date.getUTCMilliseconds()
+    ) + (dayOfYear - 1) * MILLISECONDS_IN_DAY
+  );
+}
+function getDayOfWeek(date) {
+  return date.getUTCDay();
+}
+function getDateFromDayOfWeek(dayOfWeek, dateRaw) {
+  if (dayOfWeek < 0 || dayOfWeek > 6)
+    throw new Error("Day of week must be a value between 0 and 7");
+  return new Date(
+    dateRaw.getTime() + (dayOfWeek - dateRaw.getUTCDay()) * MILLISECONDS_IN_DAY
+  );
+}
+function getDateFromHourCode(hourCode, date) {
+  const hour = RINEX_CODES.indexOf(hourCode);
+  if (hour === -1)
+    throw new Error("Hour code must be a lowercase letter between a and x");
+  const newDate = new Date(date.getTime());
+  newDate.setUTCHours(hour);
+  return newDate;
+}
+function getHourCode(date) {
+  return RINEX_CODES[date.getUTCHours()];
+}
+function getWeekOfYear(date) {
+  const target = new Date(
+    Date.UTC(date.getFullYear(), date.getMonth(), date.getDate())
+  );
+  const dayNr = (target.getUTCDay() + 6) % 7;
+  target.setUTCDate(target.getUTCDate() - dayNr + 3);
+  const firstThursday = new Date(Date.UTC(target.getUTCFullYear(), 0, 4));
+  const dayDiff = (target.getTime() - firstThursday.getTime()) / (86400 * 1e3);
+  return 1 + Math.floor(dayDiff / 7) + 1;
+}
+
+// src/time/rinex.ts
+function getRINEX(date) {
+  const seconds = date.getUTCSeconds();
+  const milliseconds = date.getUTCMilliseconds();
+  const fractionalSeconds = (milliseconds / MILLISECONDS_IN_SECOND).toFixed(7).slice(1);
+  const formatted = `> ${date.getUTCFullYear()} ${String(
+    date.getUTCMonth() + 1
+  ).padStart(2, "0")} ${String(date.getUTCDate()).padStart(2, "0")} ${String(
+    date.getUTCHours()
+  ).padStart(2, "0")} ${String(date.getUTCMinutes()).padStart(
+    2,
+    "0"
+  )} ${String(seconds).padStart(2, "0")}${fractionalSeconds}`;
+  return formatted;
+}
+function getDateFromRINEX(rinex) {
+  const year = Number.parseInt(rinex.substring(2, 6));
+  const month = Number.parseInt(rinex.substring(7, 9)) - 1;
+  const day = Number.parseInt(rinex.substring(10, 12));
+  const hour = Number.parseInt(rinex.substring(13, 15));
+  const minute = Number.parseInt(rinex.substring(16, 18));
+  const second = Number.parseInt(rinex.substring(19, 21));
+  const millisecond = Number.parseInt(rinex.substring(22, 29)) / 1e4;
+  const date = new Date(
+    Date.UTC(year, month, day, hour, minute, second, millisecond)
+  );
+  return date;
+}
+
+// src/time/julian.ts
+function getJulianDate(date, scale = 1 /* TT */) {
+  if (scale === 0 /* TAI */) {
+    return (date.getTime() - START_TAI_TIME.getTime()) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
+  }
+  if (scale === 2 /* GPS */) {
+    return (date.getTime() - START_TAI_TIME.getTime() - MILLISECONDS_GPS_TAI) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
+  }
+  if (scale === 3 /* UTC */) {
+    return (date.getTime() - START_TAI_TIME.getTime() - getLeap(date) * MILLISECONDS_IN_SECOND) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
+  }
+  return (date.getTime() - START_TAI_TIME.getTime() + MILLISECONDS_TT_TAI) / MILLISECONDS_IN_DAY + START_JULIAN_TAI;
+}
+function getDateFromJulianDate(julianDate, scale = 1 /* TT */) {
+  if (scale === 0 /* TAI */) {
+    return new Date(
+      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime()
+    );
+  }
+  if (scale === 2 /* GPS */) {
+    return new Date(
+      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime() + MILLISECONDS_GPS_TAI
+    );
+  }
+  if (scale === 3 /* UTC */) {
+    const utc_date = new Date(
+      (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime()
+    );
+    const tai_date = new Date(
+      utc_date.getTime() + getLeap(utc_date) * MILLISECONDS_IN_SECOND
+    );
+    const leap_seconds = getLeap(getDateFromTai(tai_date));
+    return new Date(utc_date.getTime() + leap_seconds * MILLISECONDS_IN_SECOND);
+  }
+  return new Date(
+    (julianDate - START_JULIAN_TAI) * MILLISECONDS_IN_DAY + START_TAI_TIME.getTime() - MILLISECONDS_TT_TAI
+  );
+}
+function getMJD(date, scale = 1 /* TT */) {
+  return getJulianDate(date, scale) - DAYS_MJD_JULIAN;
+}
+function getDateFromMJD(mjd, scale = 1 /* TT */) {
+  return getDateFromJulianDate(mjd + DAYS_MJD_JULIAN, scale);
+}
+function getMJD2000(date, scale = 1 /* TT */) {
+  return getMJD(date, scale) - DAYS_MJD2000_MJD;
+}
+function getDateFromMJD2000(mjd2000, scale = 1 /* TT */) {
+  return getDateFromMJD(mjd2000 + DAYS_MJD2000_MJD, scale);
+}
+
+export {
+  Scale,
+  getTimeOfDay,
+  getDateFromTimeOfDay,
+  getDayOfYear,
+  getDateFromDayOfYear,
+  getDayOfWeek,
+  getDateFromDayOfWeek,
+  getDateFromHourCode,
+  getHourCode,
+  getWeekOfYear,
+  getRINEX,
+  getDateFromRINEX,
+  getJulianDate,
+  getDateFromJulianDate,
+  getMJD,
+  getDateFromMJD,
+  getMJD2000,
+  getDateFromMJD2000
+};

package/dist/{chunk-PGHDJFQK.js → chunk-ZCNXERQ2.js}

@@ -1,3 +1,6 @@
+import {
+  getGpsLeap
+} from "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic
 } from "./chunk-37QNKGTC.js";
@@ -230,7 +233,8 @@ 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 = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH) / 1e3;
@@ -403,7 +407,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);
@@ -436,6 +440,65 @@ function computeLiveSkyPositions(ephemerides, rxPos, cn0Map) {
   }
   return result;
 }
+function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, elevationMaskDeg = 10) {
+  const stepMs = stepSec * 1e3;
+  const times = [];
+  for (let t = startMs; t <= endMs; t += stepMs) times.push(t);
+  const maskRad = elevationMaskDeg * Math.PI / 180;
+  const all = computeAllPositions(ephemerides, times, rxPos);
+  const elevation = {};
+  const visibleCount = new Array(times.length).fill(0);
+  const pdop = new Array(times.length).fill(null);
+  const visiblePerEpoch = times.map(() => []);
+  for (const prn of all.prns) {
+    const series = all.positions[prn];
+    elevation[prn] = series.map((p) => p ? p.el : null);
+    for (let i = 0; i < series.length; i++) {
+      const p = series[i];
+      if (p && p.el >= maskRad) {
+        visibleCount[i]++;
+        visiblePerEpoch[i].push({ az: p.az, el: p.el });
+      }
+    }
+  }
+  for (let i = 0; i < times.length; i++) {
+    const dop = computeDop(visiblePerEpoch[i]);
+    pdop[i] = dop ? dop.pdop : null;
+  }
+  const passes = [];
+  for (const prn of all.prns) {
+    const els = elevation[prn];
+    let start = -1;
+    let peakEl = -Infinity;
+    let peakTime = 0;
+    for (let i = 0; i <= els.length; i++) {
+      const el = i < els.length ? els[i] : null;
+      const above = el !== null && el !== void 0 && el >= maskRad;
+      if (above) {
+        if (start === -1) {
+          start = i;
+          peakEl = -Infinity;
+        }
+        if (el > peakEl) {
+          peakEl = el;
+          peakTime = times[i];
+        }
+      } else if (start !== -1) {
+        passes.push({
+          prn,
+          system: prn[0],
+          rise: times[start],
+          set: times[i - 1],
+          peakTime,
+          peakEl
+        });
+        start = -1;
+      }
+    }
+  }
+  passes.sort((a, b) => a.rise - b.rise);
+  return { times, elevation, visibleCount, pdop, passes };
+}
 
 export {
   keplerPosition,
@@ -446,5 +509,6 @@ export {
   computeSatPosition,
   navTimesFromEph,
   computeAllPositions,
-  computeLiveSkyPositions
+  computeLiveSkyPositions,
+  computeVisibility
 };