gnss-js 1.1.0 → 1.3.0

package/dist/index.cjs

@@ -110,6 +110,7 @@ __export(src_exports, {
   computePsdDb: () => computePsdDb,
   computeSatPosition: () => computeSatPosition,
   computeStats: () => computeStats2,
+  computeVisibility: () => computeVisibility,
   connectToMountpoint: () => connectToMountpoint,
   createStationMeta: () => createStationMeta,
   createStreamStats: () => createStreamStats,
@@ -4201,7 +4202,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_EPOCH2 = START_GPS_TIME.getTime();
   const gpsSeconds = (timeMs - GPS_EPOCH2) / 1e3;
@@ -4374,7 +4376,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);
@@ -4407,6 +4409,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 };
+}
 
 // src/ntrip/index.ts
 function parseStreamEntry(fields) {
@@ -6319,6 +6380,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   computePsdDb,
   computeSatPosition,
   computeStats,
+  computeVisibility,
   connectToMountpoint,
   createStationMeta,
   createStreamStats,

package/dist/index.d.cts

@@ -9,7 +9,7 @@ export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResul
 export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, p as parseNavFile } from './nav-BAI1a9vK.cjs';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.cjs';
 export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
-export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.cjs';
+export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';

package/dist/index.d.ts

@@ -9,7 +9,7 @@ export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResul
 export { E as Ephemeris, G as GlonassEphemeris, K as KeplerEphemeris, N as NavHeader, a as NavResult, p as parseNavFile } from './nav-BAI1a9vK.js';
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-Ohl6NAfw.js';
 export { MessageTypeStats, MsmEpoch, MsmSatObs, MsmSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, StationMeta, StreamStats, createStationMeta, createStreamStats, decodeMsmFull, msmEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
-export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.js';
+export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, analyzeQuality } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -1,3 +1,8 @@
+import {
+  connectToMountpoint,
+  fetchSourcetable,
+  parseSourcetable
+} from "./chunk-BJHTBYNG.js";
 import {
   CompletenessAccumulator,
   CycleSlipAccumulator,
@@ -51,40 +56,69 @@ import {
 } 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,
+  computeVisibility,
+  ecefToAzEl,
+  glonassPosition,
+  keplerPosition,
+  navTimesFromEph,
+  selectEphemeris
+} from "./chunk-ZCNXERQ2.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,29 +127,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 {
+  WGS84_ECCENTRICITY_SQUARED,
+  WGS84_FLATTENING,
+  WGS84_SEMI_MAJOR_AXIS,
+  WGS84_SEMI_MINOR_AXIS
+} from "./chunk-6FAL6P4G.js";
 import {
   EMPTY_WARNINGS,
   WarningAccumulator,
@@ -152,53 +184,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-JDO3LEPC.js";
-import {
-  ARC_GAP_FACTOR,
-  BAND_LABELS,
-  BDS_SATELLITES,
-  C_LIGHT,
-  DEFAULT_ELEV_MASK_DEG,
-  DUAL_FREQ_PAIRS,
-  FREQ,
-  GLO_CHANNEL_FALLBACK,
-  GLO_F1_BASE,
-  GLO_F1_STEP,
-  GLO_F2_BASE,
-  GLO_F2_STEP,
-  GLO_F3,
-  SYSTEM_NAMES,
-  SYS_SHORT,
-  buildGloChannelMap,
-  buildObsIndices,
-  formatUTCTime,
-  getFreq,
-  gloFreq
-} from "./chunk-W5WKEV7U.js";
-import {
-  computeAllPositions,
-  computeDop,
-  computeLiveSkyPositions,
-  computeSatPosition,
-  ecefToAzEl,
-  glonassPosition,
-  keplerPosition,
-  navTimesFromEph,
-  selectEphemeris
-} from "./chunk-PGHDJFQK.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";
+} from "./chunk-7EEWQ5DU.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,
@@ -225,10 +211,27 @@ import {
   START_UNIX_TIME
 } from "./chunk-LEEU5OIO.js";
 import {
-  connectToMountpoint,
-  fetchSourcetable,
-  parseSourcetable
-} from "./chunk-BJHTBYNG.js";
+  ARC_GAP_FACTOR,
+  BAND_LABELS,
+  BDS_SATELLITES,
+  C_LIGHT,
+  DEFAULT_ELEV_MASK_DEG,
+  DUAL_FREQ_PAIRS,
+  FREQ,
+  GLO_CHANNEL_FALLBACK,
+  GLO_F1_BASE,
+  GLO_F1_STEP,
+  GLO_F2_BASE,
+  GLO_F2_STEP,
+  GLO_F3,
+  SYSTEM_NAMES,
+  SYS_SHORT,
+  buildGloChannelMap,
+  buildObsIndices,
+  formatUTCTime,
+  getFreq,
+  gloFreq
+} from "./chunk-W5WKEV7U.js";
 export {
   ARC_GAP_FACTOR,
   BAND_LABELS,
@@ -320,6 +323,7 @@ export {
   computePsdDb,
   computeSatPosition,
   computeStats,
+  computeVisibility,
   connectToMountpoint,
   createStationMeta,
   createStreamStats,

package/dist/orbit.cjs

@@ -24,6 +24,7 @@ __export(orbit_exports, {
   computeDop: () => computeDop,
   computeLiveSkyPositions: () => computeLiveSkyPositions,
   computeSatPosition: () => computeSatPosition,
+  computeVisibility: () => computeVisibility,
   ecefToAzEl: () => ecefToAzEl,
   ecefToGeodetic: () => ecefToGeodetic,
   geodeticToEcef: () => geodeticToEcef,
@@ -72,8 +73,93 @@ function ecefToGeodetic(x, y, z) {
 }
 
 // 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;
@@ -299,7 +385,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;
@@ -472,7 +559,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);
@@ -505,12 +592,72 @@ 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 };
+}
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
   computeAllPositions,
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   ecefToGeodetic,
   geodeticToEcef,

package/dist/orbit.d.cts

@@ -113,5 +113,42 @@ declare function computeAllPositions(ephemerides: Ephemeris[], times: number[],
  * Returns array of SatAzEl for all satellites with valid ephemeris.
  */
 declare function computeLiveSkyPositions(ephemerides: Map<string, EphemerisInfo>, rxPos: [number, number, number], cn0Map?: Map<string, number>): SatAzEl[];
+/** A single above-mask visibility interval for one satellite. */
+interface VisibilityPass {
+    prn: string;
+    system: string;
+    /** Rise time (Unix ms) — first sample at or above the mask. */
+    rise: number;
+    /** Set time (Unix ms) — last sample at or above the mask. */
+    set: number;
+    /** Time of peak elevation (Unix ms). */
+    peakTime: number;
+    /** Peak elevation (radians). */
+    peakEl: number;
+}
+interface VisibilityResult {
+    /** Sample epochs (Unix ms). */
+    times: number[];
+    /** Elevation (radians) per PRN per epoch; null when below −0.05 rad / no eph. */
+    elevation: Record<string, (number | null)[]>;
+    /** Number of satellites at or above the mask, per epoch. */
+    visibleCount: number[];
+    /** PDOP per epoch (null when < 4 satellites above the mask). */
+    pdop: (number | null)[];
+    /** Discrete above-mask passes, sorted by rise time. */
+    passes: VisibilityPass[];
+}
+/**
+ * Predict satellite visibility and DOP over a time window for a fixed
+ * receiver location.
+ *
+ * @param ephemerides Broadcast ephemerides covering the window.
+ * @param rxPos Receiver ECEF position (meters).
+ * @param startMs Window start (Unix ms).
+ * @param endMs Window end (Unix ms).
+ * @param stepSec Sample spacing in seconds (default 300).
+ * @param elevationMaskDeg Elevation mask in degrees (default 10).
+ */
+declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };
+export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };

package/dist/orbit.d.ts

@@ -113,5 +113,42 @@ declare function computeAllPositions(ephemerides: Ephemeris[], times: number[],
  * Returns array of SatAzEl for all satellites with valid ephemeris.
  */
 declare function computeLiveSkyPositions(ephemerides: Map<string, EphemerisInfo>, rxPos: [number, number, number], cn0Map?: Map<string, number>): SatAzEl[];
+/** A single above-mask visibility interval for one satellite. */
+interface VisibilityPass {
+    prn: string;
+    system: string;
+    /** Rise time (Unix ms) — first sample at or above the mask. */
+    rise: number;
+    /** Set time (Unix ms) — last sample at or above the mask. */
+    set: number;
+    /** Time of peak elevation (Unix ms). */
+    peakTime: number;
+    /** Peak elevation (radians). */
+    peakEl: number;
+}
+interface VisibilityResult {
+    /** Sample epochs (Unix ms). */
+    times: number[];
+    /** Elevation (radians) per PRN per epoch; null when below −0.05 rad / no eph. */
+    elevation: Record<string, (number | null)[]>;
+    /** Number of satellites at or above the mask, per epoch. */
+    visibleCount: number[];
+    /** PDOP per epoch (null when < 4 satellites above the mask). */
+    pdop: (number | null)[];
+    /** Discrete above-mask passes, sorted by rise time. */
+    passes: VisibilityPass[];
+}
+/**
+ * Predict satellite visibility and DOP over a time window for a fixed
+ * receiver location.
+ *
+ * @param ephemerides Broadcast ephemerides covering the window.
+ * @param rxPos Receiver ECEF position (meters).
+ * @param startMs Window start (Unix ms).
+ * @param endMs Window end (Unix ms).
+ * @param stepSec Sample spacing in seconds (default 300).
+ * @param elevationMaskDeg Elevation mask in degrees (default 10).
+ */
+declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };
+export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris };

package/dist/orbit.js

@@ -3,12 +3,14 @@ import {
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   glonassPosition,
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-PGHDJFQK.js";
+} from "./chunk-ZCNXERQ2.js";
+import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
   geodeticToEcef
@@ -20,6 +22,7 @@ export {
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   ecefToGeodetic,
   geodeticToEcef,