gnss-js 1.2.0 → 1.3.0

package/README.md

@@ -51,6 +51,16 @@ 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

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

@@ -440,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,
@@ -450,5 +509,6 @@ export {
   computeSatPosition,
   navTimesFromEph,
   computeAllPositions,
-  computeLiveSkyPositions
+  computeLiveSkyPositions,
+  computeVisibility
 };

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,
@@ -4408,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) {
@@ -6320,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

@@ -92,12 +92,13 @@ import {
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   glonassPosition,
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-VANVV4F7.js";
+} from "./chunk-ZCNXERQ2.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -322,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,
@@ -591,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,13 @@ import {
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   glonassPosition,
   keplerPosition,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-VANVV4F7.js";
+} from "./chunk-ZCNXERQ2.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
@@ -21,6 +22,7 @@ export {
   computeDop,
   computeLiveSkyPositions,
   computeSatPosition,
+  computeVisibility,
   ecefToAzEl,
   ecefToGeodetic,
   geodeticToEcef,

package/dist/positioning.js

@@ -2,7 +2,7 @@ import {
   computeDop,
   computeSatPosition,
   ecefToAzEl
-} from "./chunk-VANVV4F7.js";
+} from "./chunk-ZCNXERQ2.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.2.0",
+  "version": "1.3.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,