gnss-js 1.8.0 → 1.10.0

package/dist/{chunk-W4YMQKWH.js → chunk-77VUSDNI.js}

@@ -610,6 +610,91 @@ function writeRinexNav(nav) {
   return header + records.join("\n") + "\n";
 }
 
+// src/rinex/ionex.ts
+function gridRange(l1, l2, dl) {
+  const out = [];
+  const n = Math.round((l2 - l1) / dl) + 1;
+  for (let i = 0; i < n; i++) out.push(l1 + i * dl);
+  return out;
+}
+function parseIonex(text) {
+  const lines = text.split("\n");
+  let exponent = -1;
+  let lats = [];
+  let lons = [];
+  const epochs = [];
+  const maps = [];
+  let i = 0;
+  for (; i < lines.length; i++) {
+    const line = lines[i];
+    const label = line.slice(60).trim();
+    if (label === "EXPONENT") {
+      exponent = parseInt(line.slice(0, 60).trim(), 10);
+    } else if (label === "LAT1 / LAT2 / DLAT") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lats = gridRange(a, b, d);
+    } else if (label === "LON1 / LON2 / DLON") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lons = gridRange(a, b, d);
+    } else if (label === "END OF HEADER") {
+      i++;
+      break;
+    }
+  }
+  if (lats.length === 0 || lons.length === 0) {
+    throw new Error("IONEX: missing grid definition");
+  }
+  const scale = Math.pow(10, exponent);
+  while (i < lines.length) {
+    const line = lines[i];
+    const label = line.slice(60).trim();
+    if (label === "START OF TEC MAP") {
+      const map = new Float32Array(lats.length * lons.length).fill(NaN);
+      let epochMs = 0;
+      i++;
+      while (i < lines.length) {
+        const l = lines[i];
+        const lab = l.slice(60).trim();
+        if (lab === "EPOCH OF CURRENT MAP") {
+          const f = l.slice(0, 60).trim().split(/\s+/).map(Number);
+          epochMs = Date.UTC(f[0], f[1] - 1, f[2], f[3], f[4], f[5]);
+          i++;
+        } else if (lab === "LAT/LON1/LON2/DLON/H") {
+          const lat = parseFloat(l.slice(2, 8));
+          const latIdx = lats.findIndex((v) => Math.abs(v - lat) < 1e-6);
+          i++;
+          let lonIdx = 0;
+          while (lonIdx < lons.length && i < lines.length) {
+            const row = lines[i];
+            for (let c = 0; c + 5 <= 80 && lonIdx < lons.length; c += 5) {
+              const vStr = row.slice(c, c + 5).trim();
+              if (vStr === "") continue;
+              const v = parseInt(vStr, 10);
+              if (latIdx >= 0) {
+                map[latIdx * lons.length + lonIdx] = v === 9999 ? NaN : v * scale;
+              }
+              lonIdx++;
+            }
+            i++;
+          }
+        } else if (lab === "END OF TEC MAP") {
+          i++;
+          break;
+        } else {
+          i++;
+        }
+      }
+      epochs.push(epochMs);
+      maps.push(map);
+    } else if (label === "START OF RMS MAP" || label === "END OF FILE") {
+      break;
+    } else {
+      i++;
+    }
+  }
+  return { epochs, lats, lons, maps };
+}
+
 export {
   padL,
   padR,
@@ -618,5 +703,6 @@ export {
   EMPTY_WARNINGS,
   WarningAccumulator,
   parseNavFile,
-  writeRinexNav
+  writeRinexNav,
+  parseIonex
 };

package/dist/{chunk-MN2DBDJL.js → chunk-AK6LNWZX.js}

@@ -437,15 +437,35 @@ function computeLiveSkyPositions(ephemerides, rxPos, cn0Map) {
   }
   return result;
 }
+function maskRadForAzimuth(maskDeg) {
+  if (typeof maskDeg === "number") {
+    const m = maskDeg * Math.PI / 180;
+    return () => m;
+  }
+  const n = maskDeg.length;
+  if (n === 0) return () => 0;
+  return (azRad) => {
+    const azDeg = (azRad * 180 / Math.PI % 360 + 360) % 360;
+    const pos = azDeg / 360 * n;
+    const i = Math.floor(pos) % n;
+    const f = pos - Math.floor(pos);
+    const a = maskDeg[i];
+    const b = maskDeg[(i + 1) % n];
+    return ((1 - f) * a + f * b) * Math.PI / 180;
+  };
+}
 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 maskRadFor = maskRadForAzimuth(elevationMaskDeg);
   const all = computeAllPositions(ephemerides, times, rxPos);
   const elevation = {};
   const azimuth = {};
+  const subLat = {};
+  const subLon = {};
   const visibleCount = new Array(times.length).fill(0);
+  const visibleBySystem = {};
   const pdop = new Array(times.length).fill(null);
   const gdop = new Array(times.length).fill(null);
   const hdop = new Array(times.length).fill(null);
@@ -455,10 +475,14 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     const series = all.positions[prn];
     elevation[prn] = series.map((p) => p ? p.el : null);
     azimuth[prn] = series.map((p) => p ? p.az : null);
+    subLat[prn] = series.map((p) => p ? p.lat : null);
+    subLon[prn] = series.map((p) => p ? p.lon : null);
     for (let i = 0; i < series.length; i++) {
       const p = series[i];
-      if (p && p.el >= maskRad) {
+      if (p && p.el >= maskRadFor(p.az)) {
         visibleCount[i]++;
+        const sys = prn[0];
+        (visibleBySystem[sys] ??= new Array(times.length).fill(0))[i]++;
         visiblePerEpoch[i].push({ az: p.az, el: p.el });
       }
     }
@@ -478,7 +502,8 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     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;
+      const az = all.positions[prn][i]?.az ?? 0;
+      const above = el !== null && el !== void 0 && el >= maskRadFor(az);
       if (above) {
         if (start === -1) {
           start = i;
@@ -506,7 +531,10 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     times,
     elevation,
     azimuth,
+    subLat,
+    subLon,
     visibleCount,
+    visibleBySystem,
     pdop,
     gdop,
     hdop,
@@ -526,5 +554,6 @@ export {
   computeAllPositions,
   ephInfoToEphemeris,
   computeLiveSkyPositions,
+  maskRadForAzimuth,
   computeVisibility
 };

package/dist/{chunk-GXK4MMHS.js → chunk-ZTFKRZMM.js}

@@ -2,7 +2,7 @@ import {
   computeDop,
   computeSatPosition,
   ecefToAzEl
-} from "./chunk-MN2DBDJL.js";
+} from "./chunk-AK6LNWZX.js";
 import {
   C_LIGHT,
   OMEGA_E

package/dist/index.cjs

@@ -200,6 +200,7 @@ __export(src_exports, {
   horizonDistance: () => horizonDistance,
   ionoFree: () => ionoFree,
   keplerPosition: () => keplerPosition,
+  maskRadForAzimuth: () => maskRadForAzimuth,
   msmEpochToDate: () => msmEpochToDate,
   navTimesFromEph: () => navTimesFromEph,
   nmeaCoordToDecimal: () => nmeaCoordToDecimal,
@@ -207,6 +208,7 @@ __export(src_exports, {
   padR: () => padR,
   parseAntex: () => parseAntex,
   parseCrxDataLine: () => parseCrxDataLine,
+  parseIonex: () => parseIonex,
   parseNavFile: () => parseNavFile,
   parseNmeaFile: () => parseNmeaFile,
   parseRinexStream: () => parseRinexStream,
@@ -2575,6 +2577,91 @@ function writeRinexNav(nav) {
   return header + records.join("\n") + "\n";
 }
 
+// src/rinex/ionex.ts
+function gridRange(l1, l2, dl) {
+  const out = [];
+  const n = Math.round((l2 - l1) / dl) + 1;
+  for (let i = 0; i < n; i++) out.push(l1 + i * dl);
+  return out;
+}
+function parseIonex(text) {
+  const lines = text.split("\n");
+  let exponent = -1;
+  let lats = [];
+  let lons = [];
+  const epochs = [];
+  const maps = [];
+  let i = 0;
+  for (; i < lines.length; i++) {
+    const line = lines[i];
+    const label2 = line.slice(60).trim();
+    if (label2 === "EXPONENT") {
+      exponent = parseInt(line.slice(0, 60).trim(), 10);
+    } else if (label2 === "LAT1 / LAT2 / DLAT") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lats = gridRange(a, b, d);
+    } else if (label2 === "LON1 / LON2 / DLON") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lons = gridRange(a, b, d);
+    } else if (label2 === "END OF HEADER") {
+      i++;
+      break;
+    }
+  }
+  if (lats.length === 0 || lons.length === 0) {
+    throw new Error("IONEX: missing grid definition");
+  }
+  const scale = Math.pow(10, exponent);
+  while (i < lines.length) {
+    const line = lines[i];
+    const label2 = line.slice(60).trim();
+    if (label2 === "START OF TEC MAP") {
+      const map = new Float32Array(lats.length * lons.length).fill(NaN);
+      let epochMs = 0;
+      i++;
+      while (i < lines.length) {
+        const l = lines[i];
+        const lab = l.slice(60).trim();
+        if (lab === "EPOCH OF CURRENT MAP") {
+          const f = l.slice(0, 60).trim().split(/\s+/).map(Number);
+          epochMs = Date.UTC(f[0], f[1] - 1, f[2], f[3], f[4], f[5]);
+          i++;
+        } else if (lab === "LAT/LON1/LON2/DLON/H") {
+          const lat = parseFloat(l.slice(2, 8));
+          const latIdx = lats.findIndex((v) => Math.abs(v - lat) < 1e-6);
+          i++;
+          let lonIdx = 0;
+          while (lonIdx < lons.length && i < lines.length) {
+            const row = lines[i];
+            for (let c = 0; c + 5 <= 80 && lonIdx < lons.length; c += 5) {
+              const vStr = row.slice(c, c + 5).trim();
+              if (vStr === "") continue;
+              const v = parseInt(vStr, 10);
+              if (latIdx >= 0) {
+                map[latIdx * lons.length + lonIdx] = v === 9999 ? NaN : v * scale;
+              }
+              lonIdx++;
+            }
+            i++;
+          }
+        } else if (lab === "END OF TEC MAP") {
+          i++;
+          break;
+        } else {
+          i++;
+        }
+      }
+      epochs.push(epochMs);
+      maps.push(map);
+    } else if (label2 === "START OF RMS MAP" || label2 === "END OF FILE") {
+      break;
+    } else {
+      i++;
+    }
+  }
+  return { epochs, lats, lons, maps };
+}
+
 // src/rtcm3/decoder.ts
 var debugHandler = null;
 function setRtcm3DebugHandler(fn) {
@@ -4427,15 +4514,35 @@ function computeLiveSkyPositions(ephemerides, rxPos, cn0Map) {
   }
   return result;
 }
+function maskRadForAzimuth(maskDeg) {
+  if (typeof maskDeg === "number") {
+    const m = maskDeg * Math.PI / 180;
+    return () => m;
+  }
+  const n = maskDeg.length;
+  if (n === 0) return () => 0;
+  return (azRad) => {
+    const azDeg = (azRad * 180 / Math.PI % 360 + 360) % 360;
+    const pos = azDeg / 360 * n;
+    const i = Math.floor(pos) % n;
+    const f = pos - Math.floor(pos);
+    const a = maskDeg[i];
+    const b = maskDeg[(i + 1) % n];
+    return ((1 - f) * a + f * b) * Math.PI / 180;
+  };
+}
 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 maskRadFor = maskRadForAzimuth(elevationMaskDeg);
   const all = computeAllPositions(ephemerides, times, rxPos);
   const elevation = {};
   const azimuth = {};
+  const subLat = {};
+  const subLon = {};
   const visibleCount = new Array(times.length).fill(0);
+  const visibleBySystem = {};
   const pdop = new Array(times.length).fill(null);
   const gdop = new Array(times.length).fill(null);
   const hdop = new Array(times.length).fill(null);
@@ -4445,10 +4552,14 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     const series = all.positions[prn];
     elevation[prn] = series.map((p) => p ? p.el : null);
     azimuth[prn] = series.map((p) => p ? p.az : null);
+    subLat[prn] = series.map((p) => p ? p.lat : null);
+    subLon[prn] = series.map((p) => p ? p.lon : null);
     for (let i = 0; i < series.length; i++) {
       const p = series[i];
-      if (p && p.el >= maskRad) {
+      if (p && p.el >= maskRadFor(p.az)) {
         visibleCount[i]++;
+        const sys = prn[0];
+        (visibleBySystem[sys] ??= new Array(times.length).fill(0))[i]++;
         visiblePerEpoch[i].push({ az: p.az, el: p.el });
       }
     }
@@ -4468,7 +4579,8 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     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;
+      const az = all.positions[prn][i]?.az ?? 0;
+      const above = el !== null && el !== void 0 && el >= maskRadFor(az);
       if (above) {
         if (start === -1) {
           start = i;
@@ -4496,7 +4608,10 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     times,
     elevation,
     azimuth,
+    subLat,
+    subLon,
     visibleCount,
+    visibleBySystem,
     pdop,
     gdop,
     hdop,
@@ -7353,6 +7468,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   horizonDistance,
   ionoFree,
   keplerPosition,
+  maskRadForAzimuth,
   msmEpochToDate,
   navTimesFromEph,
   nmeaCoordToDecimal,
@@ -7360,6 +7476,7 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   padR,
   parseAntex,
   parseCrxDataLine,
+  parseIonex,
   parseNavFile,
   parseNmeaFile,
   parseRinexStream,

package/dist/index.d.cts

@@ -4,12 +4,12 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.cjs';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.cjs';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.cjs';
-export { CrxField, DiffState, EMPTY_WARNINGS, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav } from './rinex.cjs';
+export { CrxField, DiffState, EMPTY_WARNINGS, IonexGrid, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, writeRinexNav } from './rinex.cjs';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.cjs';
 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, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.cjs';
+export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } from './positioning.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';

package/dist/index.d.ts

@@ -4,12 +4,12 @@ export { A as ARC_GAP_FACTOR, B as BAND_LABELS, a as BDS_SATELLITES, b as BdsSat
 export { Scale, getBdsTime, getDateFromBdsTime, getDateFromDayOfWeek, getDateFromDayOfYear, getDateFromGalTime, getDateFromGloN, getDateFromGpsData, getDateFromGpsTime, getDateFromHourCode, getDateFromJulianDate, getDateFromMJD, getDateFromMJD2000, getDateFromRINEX, getDateFromTai, getDateFromTimeOfDay, getDateFromTt, getDateFromUnixTime, getDateFromUtc, getDayOfWeek, getDayOfYear, getGalTime, getGloN4, getGloNA, getGpsLeap, getGpsTime, getHourCode, getHoursFromTimeDifference, getJulianDate, getLeap, getMJD, getMJD2000, getMinutesFromTimeDifference, getNtpTime, getRINEX, getSecondsFromTimeDifference, getTaiDate, getTimeDifference, getTimeDifferenceFromDays, getTimeDifferenceFromHours, getTimeDifferenceFromMinutes, getTimeDifferenceFromObject, getTimeDifferenceFromSeconds, getTimeOfDay, getTimeOfWeek, getTotalDaysFromTimeDifference, getTtDate, getUnixTime, getUtcDate, getWeekNumber, getWeekOfYear } from './time.js';
 export { deg2dms, deg2rad, euclidean3D, geodeticToGeohash, geodeticToMaidenhead, geodeticToUtm, greatCircleMidpoint, horizonDistance, rad2deg, rhumbLine, vincenty } from './coordinates.js';
 export { c as clampUnit, e as ecefToGeodetic, g as geodeticToEcef, a as getAer, b as getEnuDifference } from './ecef-CF0uAysr.js';
-export { CrxField, DiffState, EMPTY_WARNINGS, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav } from './rinex.js';
+export { CrxField, DiffState, EMPTY_WARNINGS, IonexGrid, RinexWarning, RinexWarnings, WarningAccumulator, WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, writeRinexNav } from './rinex.js';
 export { E as EpochSummary, P as ParseOptions, R as RinexHeader, a as RinexResult, b as RinexStats, S as SYSTEM_ORDER, c as SatObsCallback, p as parseRinexStream, s as systemCmp, d as systemName } from './parser-JPjjFgeP.js';
 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, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, navTimesFromEph, selectEphemeris } from './orbit.js';
+export { AllPositionsData, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
 export { SppOptions, SppSolution, ionoFree, satClockCorrection, solveSpp } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';

package/dist/index.js

@@ -102,7 +102,7 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-GXK4MMHS.js";
+} from "./chunk-ZTFKRZMM.js";
 import {
   computeAllPositions,
   computeDop,
@@ -113,9 +113,10 @@ import {
   ephInfoToEphemeris,
   glonassPosition,
   keplerPosition,
+  maskRadForAzimuth,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-MN2DBDJL.js";
+} from "./chunk-AK6LNWZX.js";
 import {
   getBdsTime,
   getDateFromBdsTime,
@@ -172,9 +173,10 @@ import {
   hdrLine,
   padL,
   padR,
+  parseIonex,
   parseNavFile,
   writeRinexNav
-} from "./chunk-W4YMQKWH.js";
+} from "./chunk-77VUSDNI.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -431,6 +433,7 @@ export {
   horizonDistance,
   ionoFree,
   keplerPosition,
+  maskRadForAzimuth,
   msmEpochToDate,
   navTimesFromEph,
   nmeaCoordToDecimal,
@@ -438,6 +441,7 @@ export {
   padR,
   parseAntex,
   parseCrxDataLine,
+  parseIonex,
   parseNavFile,
   parseNmeaFile,
   parseRinexStream,

package/dist/orbit.cjs

@@ -31,6 +31,7 @@ __export(orbit_exports, {
   geodeticToEcef: () => geodeticToEcef,
   glonassPosition: () => glonassPosition,
   keplerPosition: () => keplerPosition,
+  maskRadForAzimuth: () => maskRadForAzimuth,
   navTimesFromEph: () => navTimesFromEph,
   selectEphemeris: () => selectEphemeris
 });
@@ -590,15 +591,35 @@ function computeLiveSkyPositions(ephemerides, rxPos, cn0Map) {
   }
   return result;
 }
+function maskRadForAzimuth(maskDeg) {
+  if (typeof maskDeg === "number") {
+    const m = maskDeg * Math.PI / 180;
+    return () => m;
+  }
+  const n = maskDeg.length;
+  if (n === 0) return () => 0;
+  return (azRad) => {
+    const azDeg = (azRad * 180 / Math.PI % 360 + 360) % 360;
+    const pos = azDeg / 360 * n;
+    const i = Math.floor(pos) % n;
+    const f = pos - Math.floor(pos);
+    const a = maskDeg[i];
+    const b = maskDeg[(i + 1) % n];
+    return ((1 - f) * a + f * b) * Math.PI / 180;
+  };
+}
 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 maskRadFor = maskRadForAzimuth(elevationMaskDeg);
   const all = computeAllPositions(ephemerides, times, rxPos);
   const elevation = {};
   const azimuth = {};
+  const subLat = {};
+  const subLon = {};
   const visibleCount = new Array(times.length).fill(0);
+  const visibleBySystem = {};
   const pdop = new Array(times.length).fill(null);
   const gdop = new Array(times.length).fill(null);
   const hdop = new Array(times.length).fill(null);
@@ -608,10 +629,14 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     const series = all.positions[prn];
     elevation[prn] = series.map((p) => p ? p.el : null);
     azimuth[prn] = series.map((p) => p ? p.az : null);
+    subLat[prn] = series.map((p) => p ? p.lat : null);
+    subLon[prn] = series.map((p) => p ? p.lon : null);
     for (let i = 0; i < series.length; i++) {
       const p = series[i];
-      if (p && p.el >= maskRad) {
+      if (p && p.el >= maskRadFor(p.az)) {
         visibleCount[i]++;
+        const sys = prn[0];
+        (visibleBySystem[sys] ??= new Array(times.length).fill(0))[i]++;
         visiblePerEpoch[i].push({ az: p.az, el: p.el });
       }
     }
@@ -631,7 +656,8 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     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;
+      const az = all.positions[prn][i]?.az ?? 0;
+      const above = el !== null && el !== void 0 && el >= maskRadFor(az);
       if (above) {
         if (start === -1) {
           start = i;
@@ -659,7 +685,10 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
     times,
     elevation,
     azimuth,
+    subLat,
+    subLon,
     visibleCount,
+    visibleBySystem,
     pdop,
     gdop,
     hdop,
@@ -680,6 +709,7 @@ function computeVisibility(ephemerides, rxPos, startMs, endMs, stepSec = 300, el
   geodeticToEcef,
   glonassPosition,
   keplerPosition,
+  maskRadForAzimuth,
   navTimesFromEph,
   selectEphemeris
 });

package/dist/orbit.d.cts

@@ -135,8 +135,14 @@ interface VisibilityResult {
     elevation: Record<string, (number | null)[]>;
     /** Azimuth (radians) per PRN per epoch; null when no valid ephemeris. */
     azimuth: Record<string, (number | null)[]>;
+    /** Sub-satellite latitude (radians) per PRN per epoch. */
+    subLat: Record<string, (number | null)[]>;
+    /** Sub-satellite longitude (radians) per PRN per epoch. */
+    subLon: Record<string, (number | null)[]>;
     /** Number of satellites at or above the mask, per epoch. */
     visibleCount: number[];
+    /** Per-system visible counts per epoch (system letter → counts). */
+    visibleBySystem: Record<string, number[]>;
     /** PDOP per epoch (null when < 4 satellites above the mask). */
     pdop: (number | null)[];
     /** GDOP per epoch (null when < 4 satellites above the mask). */
@@ -148,6 +154,12 @@ interface VisibilityResult {
     /** Discrete above-mask passes, sorted by rise time. */
     passes: VisibilityPass[];
 }
+/**
+ * Build a mask lookup (radians) from a scalar or per-azimuth array of
+ * mask values in degrees. Array entries span 0–360° uniformly and are
+ * linearly interpolated.
+ */
+declare function maskRadForAzimuth(maskDeg: number | number[]): (azRad: number) => number;
 /**
  * Predict satellite visibility and DOP over a time window for a fixed
  * receiver location.
@@ -157,8 +169,11 @@ interface VisibilityResult {
  * @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).
+ * @param elevationMaskDeg Elevation mask in degrees (default 10) — a
+ *   scalar, or an array of per-azimuth mask values in degrees
+ *   (uniformly spanning 0–360°, e.g. 36 sectors of 10°) describing a
+ *   local horizon/obstruction profile.
  */
-declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
+declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number | number[]): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, 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, ephInfoToEphemeris, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, selectEphemeris };

package/dist/orbit.d.ts

@@ -135,8 +135,14 @@ interface VisibilityResult {
     elevation: Record<string, (number | null)[]>;
     /** Azimuth (radians) per PRN per epoch; null when no valid ephemeris. */
     azimuth: Record<string, (number | null)[]>;
+    /** Sub-satellite latitude (radians) per PRN per epoch. */
+    subLat: Record<string, (number | null)[]>;
+    /** Sub-satellite longitude (radians) per PRN per epoch. */
+    subLon: Record<string, (number | null)[]>;
     /** Number of satellites at or above the mask, per epoch. */
     visibleCount: number[];
+    /** Per-system visible counts per epoch (system letter → counts). */
+    visibleBySystem: Record<string, number[]>;
     /** PDOP per epoch (null when < 4 satellites above the mask). */
     pdop: (number | null)[];
     /** GDOP per epoch (null when < 4 satellites above the mask). */
@@ -148,6 +154,12 @@ interface VisibilityResult {
     /** Discrete above-mask passes, sorted by rise time. */
     passes: VisibilityPass[];
 }
+/**
+ * Build a mask lookup (radians) from a scalar or per-azimuth array of
+ * mask values in degrees. Array entries span 0–360° uniformly and are
+ * linearly interpolated.
+ */
+declare function maskRadForAzimuth(maskDeg: number | number[]): (azRad: number) => number;
 /**
  * Predict satellite visibility and DOP over a time window for a fixed
  * receiver location.
@@ -157,8 +169,11 @@ interface VisibilityResult {
  * @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).
+ * @param elevationMaskDeg Elevation mask in degrees (default 10) — a
+ *   scalar, or an array of per-azimuth mask values in degrees
+ *   (uniformly spanning 0–360°, e.g. 36 sectors of 10°) describing a
+ *   local horizon/obstruction profile.
  */
-declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number): VisibilityResult;
+declare function computeVisibility(ephemerides: Ephemeris[], rxPos: [number, number, number], startMs: number, endMs: number, stepSec?: number, elevationMaskDeg?: number | number[]): VisibilityResult;
 
-export { type AllPositionsData, type DopValues, type EpochSkyData, type SatAzEl, type SatPoint, type SatPosition, type VisibilityPass, type VisibilityResult, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, ecefToAzEl, ephInfoToEphemeris, 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, ephInfoToEphemeris, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, selectEphemeris };

package/dist/orbit.js

@@ -8,9 +8,10 @@ import {
   ephInfoToEphemeris,
   glonassPosition,
   keplerPosition,
+  maskRadForAzimuth,
   navTimesFromEph,
   selectEphemeris
-} from "./chunk-MN2DBDJL.js";
+} from "./chunk-AK6LNWZX.js";
 import "./chunk-HVXYFUCB.js";
 import {
   ecefToGeodetic,
@@ -31,6 +32,7 @@ export {
   geodeticToEcef,
   glonassPosition,
   keplerPosition,
+  maskRadForAzimuth,
   navTimesFromEph,
   selectEphemeris
 };

package/dist/positioning.js

@@ -2,8 +2,8 @@ import {
   ionoFree,
   satClockCorrection,
   solveSpp
-} from "./chunk-GXK4MMHS.js";
-import "./chunk-MN2DBDJL.js";
+} from "./chunk-ZTFKRZMM.js";
+import "./chunk-AK6LNWZX.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";

package/dist/rinex.cjs

@@ -30,6 +30,7 @@ __export(rinex_exports, {
   padL: () => padL,
   padR: () => padR,
   parseCrxDataLine: () => parseCrxDataLine,
+  parseIonex: () => parseIonex,
   parseNavFile: () => parseNavFile,
   parseRinexStream: () => parseRinexStream,
   systemCmp: () => systemCmp,
@@ -1449,6 +1450,91 @@ function writeRinexNav(nav) {
   );
   return header + records.join("\n") + "\n";
 }
+
+// src/rinex/ionex.ts
+function gridRange(l1, l2, dl) {
+  const out = [];
+  const n = Math.round((l2 - l1) / dl) + 1;
+  for (let i = 0; i < n; i++) out.push(l1 + i * dl);
+  return out;
+}
+function parseIonex(text) {
+  const lines = text.split("\n");
+  let exponent = -1;
+  let lats = [];
+  let lons = [];
+  const epochs = [];
+  const maps = [];
+  let i = 0;
+  for (; i < lines.length; i++) {
+    const line = lines[i];
+    const label = line.slice(60).trim();
+    if (label === "EXPONENT") {
+      exponent = parseInt(line.slice(0, 60).trim(), 10);
+    } else if (label === "LAT1 / LAT2 / DLAT") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lats = gridRange(a, b, d);
+    } else if (label === "LON1 / LON2 / DLON") {
+      const [a, b, d] = line.trim().split(/\s+/).map(Number);
+      lons = gridRange(a, b, d);
+    } else if (label === "END OF HEADER") {
+      i++;
+      break;
+    }
+  }
+  if (lats.length === 0 || lons.length === 0) {
+    throw new Error("IONEX: missing grid definition");
+  }
+  const scale = Math.pow(10, exponent);
+  while (i < lines.length) {
+    const line = lines[i];
+    const label = line.slice(60).trim();
+    if (label === "START OF TEC MAP") {
+      const map = new Float32Array(lats.length * lons.length).fill(NaN);
+      let epochMs = 0;
+      i++;
+      while (i < lines.length) {
+        const l = lines[i];
+        const lab = l.slice(60).trim();
+        if (lab === "EPOCH OF CURRENT MAP") {
+          const f = l.slice(0, 60).trim().split(/\s+/).map(Number);
+          epochMs = Date.UTC(f[0], f[1] - 1, f[2], f[3], f[4], f[5]);
+          i++;
+        } else if (lab === "LAT/LON1/LON2/DLON/H") {
+          const lat = parseFloat(l.slice(2, 8));
+          const latIdx = lats.findIndex((v) => Math.abs(v - lat) < 1e-6);
+          i++;
+          let lonIdx = 0;
+          while (lonIdx < lons.length && i < lines.length) {
+            const row = lines[i];
+            for (let c = 0; c + 5 <= 80 && lonIdx < lons.length; c += 5) {
+              const vStr = row.slice(c, c + 5).trim();
+              if (vStr === "") continue;
+              const v = parseInt(vStr, 10);
+              if (latIdx >= 0) {
+                map[latIdx * lons.length + lonIdx] = v === 9999 ? NaN : v * scale;
+              }
+              lonIdx++;
+            }
+            i++;
+          }
+        } else if (lab === "END OF TEC MAP") {
+          i++;
+          break;
+        } else {
+          i++;
+        }
+      }
+      epochs.push(epochMs);
+      maps.push(map);
+    } else if (label === "START OF RMS MAP" || label === "END OF FILE") {
+      break;
+    } else {
+      i++;
+    }
+  }
+  return { epochs, lats, lons, maps };
+}
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
   EMPTY_WARNINGS,
@@ -1461,6 +1547,7 @@ function writeRinexNav(nav) {
   padL,
   padR,
   parseCrxDataLine,
+  parseIonex,
   parseNavFile,
   parseRinexStream,
   systemCmp,

package/dist/rinex.d.cts

@@ -113,4 +113,26 @@ declare class WarningAccumulator {
  */
 declare function writeRinexNav(nav: NavResult): string;
 
-export { type CrxField, type DiffState, EMPTY_WARNINGS, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav };
+/**
+ * IONEX 1.0 parser — global ionosphere maps (GIMs) as published by
+ * the IGS analysis centres (e.g. ESA0OPSRAP_*_GIM.INX).
+ *
+ * Returns the TEC maps in TECU on a regular lat/lon grid, one map per
+ * epoch. RMS maps and the DCB auxiliary block are skipped.
+ */
+interface IonexGrid {
+    /** Map epochs (Unix ms, UTC). */
+    epochs: number[];
+    /** Grid latitudes (degrees), in file order (typically 87.5 → −87.5). */
+    lats: number[];
+    /** Grid longitudes (degrees), in file order (typically −180 → 180). */
+    lons: number[];
+    /**
+     * TEC maps in TECU: maps[epochIdx][latIdx * lons.length + lonIdx].
+     * NaN where the file marks no value (9999).
+     */
+    maps: Float32Array[];
+}
+declare function parseIonex(text: string): IonexGrid;
+
+export { type CrxField, type DiffState, EMPTY_WARNINGS, type IonexGrid, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, writeRinexNav };

package/dist/rinex.d.ts

@@ -113,4 +113,26 @@ declare class WarningAccumulator {
  */
 declare function writeRinexNav(nav: NavResult): string;
 
-export { type CrxField, type DiffState, EMPTY_WARNINGS, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, writeRinexNav };
+/**
+ * IONEX 1.0 parser — global ionosphere maps (GIMs) as published by
+ * the IGS analysis centres (e.g. ESA0OPSRAP_*_GIM.INX).
+ *
+ * Returns the TEC maps in TECU on a regular lat/lon grid, one map per
+ * epoch. RMS maps and the DCB auxiliary block are skipped.
+ */
+interface IonexGrid {
+    /** Map epochs (Unix ms, UTC). */
+    epochs: number[];
+    /** Grid latitudes (degrees), in file order (typically 87.5 → −87.5). */
+    lats: number[];
+    /** Grid longitudes (degrees), in file order (typically −180 → 180). */
+    lons: number[];
+    /**
+     * TEC maps in TECU: maps[epochIdx][latIdx * lons.length + lonIdx].
+     * NaN where the file marks no value (9999).
+     */
+    maps: Float32Array[];
+}
+declare function parseIonex(text: string): IonexGrid;
+
+export { type CrxField, type DiffState, EMPTY_WARNINGS, type IonexGrid, NavResult, RinexHeader, type RinexWarning, type RinexWarnings, WarningAccumulator, type WarningSeverity, crxDecompress, crxRepair, fmtF, hdrLine, padL, padR, parseCrxDataLine, parseIonex, writeRinexNav };

package/dist/rinex.js

@@ -5,9 +5,10 @@ import {
   hdrLine,
   padL,
   padR,
+  parseIonex,
   parseNavFile,
   writeRinexNav
-} from "./chunk-W4YMQKWH.js";
+} from "./chunk-77VUSDNI.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -29,6 +30,7 @@ export {
   padL,
   padR,
   parseCrxDataLine,
+  parseIonex,
   parseNavFile,
   parseRinexStream,
   systemCmp,

package/package.json

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