gnss-js 1.50.0 → 1.51.0

package/dist/{chunk-XFPDTGSV.js → chunk-DE5IUULR.js}

@@ -705,7 +705,7 @@ function parseSbfGloNav(data) {
 }
 
 // src/sbf/rawnav-sbas.ts
-function feedGeoBlock(view, b) {
+function feedGeoBlock(view, b, onSbasMessage) {
   const tow = view.getUint32(b + 8, true);
   const wnc = view.getUint16(b + 12, true);
   const svid = view.getUint8(b + 14);
@@ -721,10 +721,11 @@ function feedGeoBlock(view, b) {
   }
   if (!sbasCrcOk(msg)) return { badCrc: true };
   const prn = Number(prnStr.slice(1)) + 100;
+  onSbasMessage?.(msg, prn, wnc, tow / 1e3);
   const eph = decodeSbasGeoNav(msg, prn, wnc, tow / 1e3);
   return eph ? { eph } : {};
 }
-function parseSbfGeoNav(data) {
+function parseSbfGeoNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -733,7 +734,7 @@ function parseSbfGeoNav(data) {
   scanSbfFrames(data, view, (id, b, len) => {
     if (id !== 4020 || len < 52) return;
     messages++;
-    const r = feedGeoBlock(view, b);
+    const r = feedGeoBlock(view, b, opts.onSbasMessage);
     if (r.eph) {
       const key = `${r.eph.prn}|${r.eph.tocDate.getTime()}`;
       if (!seen.has(key)) {
@@ -753,7 +754,7 @@ function ephKey(e) {
   const source = "source" in e ? e.source : "";
   return `${e.prn}|${e.tocDate.getTime()}|${source}`;
 }
-function decodeSbfNavigation(data) {
+function decodeSbfNavigation(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -874,7 +875,7 @@ function decodeSbfNavigation(data) {
       case 4020: {
         if (len < 52) return;
         counts.sbasRaw++;
-        const r = feedGeoBlock(view, b);
+        const r = feedGeoBlock(view, b, opts.onSbasMessage);
         if (r.eph) addEph(r.eph);
         else if (r.badCrc) counts.badFrames++;
         return;

package/dist/{chunk-VYC2FRDV.js → chunk-E3QCNJL4.js}

@@ -3074,9 +3074,11 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
     convergenceM = 1e-4,
     tgd = true,
     iono,
-    gim
+    gim,
+    sbas
   } = opts;
   const gpsTow = (timeMs - GPS_EPOCH_MS_SPP) / 1e3 % 604800;
+  const gpsWeek = Math.floor((timeMs - GPS_EPOCH_MS_SPP) / 1e3 / 604800);
   const all = [...pseudoranges.keys()].filter((p) => ephemerides.has(p));
   const minSats = (list) => 3 + new Set(list.map((p) => p[0])).size;
   if (all.length < minSats(all)) return null;
@@ -3113,7 +3115,16 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         const sat = computeSatPosition(eph, tTx - dtsClock * 1e3);
         if (!Number.isFinite(sat.x)) continue;
         const isKepler = eph.system !== "R" && eph.system !== "S";
-        const dts = dtsClock - (tgd && isKepler ? eph.tgd : 0);
+        let dts = dtsClock - (tgd && isKepler ? eph.tgd : 0);
+        if (sbas) {
+          const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+          if (sc) {
+            sat.x += sc.dPos[0];
+            sat.y += sc.dPos[1];
+            sat.z += sc.dPos[2];
+            dts += sc.dClkS;
+          }
+        }
         const travel = Math.hypot(sat.x - x9, sat.y - y2, sat.z - z2) / C_LIGHT;
         const [sx, sy, sz] = sagnac2(sat, travel);
         const rho = Math.hypot(sx - x9, sy - y2, sz - z2);
@@ -3134,9 +3145,21 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         }
         const tropo = troposphere && positionSane ? tropoDelay2(elev, rxLat, rxHgt) : 0;
         let ionoM = 0;
-        if ((gim || iono) && positionSane) {
+        if ((sbas || gim || iono) && positionSane) {
           let l1M = null;
-          if (gim)
+          if (sbas) {
+            const si = sbas.ionoDelay(
+              gpsWeek,
+              gpsTow,
+              rxLat,
+              rxLon,
+              rxHgt,
+              azim,
+              elev
+            );
+            if (si) l1M = si.delayM;
+          }
+          if (l1M === null && gim)
             l1M = gimSlantIonoDelayL1(gim, rxLat, rxLon, azim, elev, timeMs);
           if (l1M === null && iono) {
             l1M = C_LIGHT * klobucharDelay(iono, rxLat, rxLon, azim, elev, gpsTow);

package/dist/{chunk-2IM5JYKP.js → chunk-ZS5RLFWI.js}

@@ -261,7 +261,7 @@ function decodeGpsQzssEphemeris(view, p, sys) {
     tgd: view.getFloat64(p + 172, true)
   };
 }
-function parseNovatelNav(data) {
+function parseNovatelNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const stats = { badCrc: 0 };
   const ephemerides = [];
@@ -334,6 +334,7 @@ function parseNovatelNav(data) {
       const prn = view.getUint32(p + 4, true);
       const msg = new Uint8Array(32);
       msg.set(data.subarray(p + 12, p + 41));
+      opts.onSbasMessage?.(msg, prn, frame.week, frame.towMs / 1e3);
       const eph = decodeSbasGeoNav(msg, prn, frame.week, frame.towMs / 1e3);
       if (!eph) continue;
       const prev = lastGlo.get(eph.prn);

package/dist/index.cjs

@@ -7727,7 +7727,7 @@ function parseSbfGloNav(data) {
 }
 
 // src/sbf/rawnav-sbas.ts
-function feedGeoBlock(view, b) {
+function feedGeoBlock(view, b, onSbasMessage) {
   const tow = view.getUint32(b + 8, true);
   const wnc = view.getUint16(b + 12, true);
   const svid = view.getUint8(b + 14);
@@ -7743,10 +7743,11 @@ function feedGeoBlock(view, b) {
   }
   if (!sbasCrcOk(msg)) return { badCrc: true };
   const prn = Number(prnStr.slice(1)) + 100;
+  onSbasMessage?.(msg, prn, wnc, tow / 1e3);
   const eph = decodeSbasGeoNav(msg, prn, wnc, tow / 1e3);
   return eph ? { eph } : {};
 }
-function parseSbfGeoNav(data) {
+function parseSbfGeoNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -7755,7 +7756,7 @@ function parseSbfGeoNav(data) {
   scanSbfFrames(data, view, (id, b, len) => {
     if (id !== 4020 || len < 52) return;
     messages++;
-    const r = feedGeoBlock(view, b);
+    const r = feedGeoBlock(view, b, opts.onSbasMessage);
     if (r.eph) {
       const key = `${r.eph.prn}|${r.eph.tocDate.getTime()}`;
       if (!seen.has(key)) {
@@ -7775,7 +7776,7 @@ function ephKey(e) {
   const source = "source" in e ? e.source : "";
   return `${e.prn}|${e.tocDate.getTime()}|${source}`;
 }
-function decodeSbfNavigation(data) {
+function decodeSbfNavigation(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -7896,7 +7897,7 @@ function decodeSbfNavigation(data) {
       case 4020: {
         if (len < 52) return;
         counts.sbasRaw++;
-        const r = feedGeoBlock(view, b);
+        const r = feedGeoBlock(view, b, opts.onSbasMessage);
         if (r.eph) addEph(r.eph);
         else if (r.badCrc) counts.badFrames++;
         return;
@@ -9368,7 +9369,7 @@ function decodeGpsQzssEphemeris(view, p, sys) {
     tgd: view.getFloat64(p + 172, true)
   };
 }
-function parseNovatelNav(data) {
+function parseNovatelNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const stats = { badCrc: 0 };
   const ephemerides = [];
@@ -9441,6 +9442,7 @@ function parseNovatelNav(data) {
       const prn = view.getUint32(p + 4, true);
       const msg = new Uint8Array(32);
       msg.set(data.subarray(p + 12, p + 41));
+      opts.onSbasMessage?.(msg, prn, frame.week, frame.towMs / 1e3);
       const eph = decodeSbasGeoNav(msg, prn, frame.week, frame.towMs / 1e3);
       if (!eph) continue;
       const prev = lastGlo.get(eph.prn);
@@ -15002,9 +15004,11 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
     convergenceM = 1e-4,
     tgd = true,
     iono,
-    gim
+    gim,
+    sbas
   } = opts;
   const gpsTow = (timeMs - GPS_EPOCH_MS_SPP) / 1e3 % 604800;
+  const gpsWeek = Math.floor((timeMs - GPS_EPOCH_MS_SPP) / 1e3 / 604800);
   const all = [...pseudoranges.keys()].filter((p) => ephemerides.has(p));
   const minSats = (list) => 3 + new Set(list.map((p) => p[0])).size;
   if (all.length < minSats(all)) return null;
@@ -15041,7 +15045,16 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         const sat = computeSatPosition(eph, tTx - dtsClock * 1e3);
         if (!Number.isFinite(sat.x)) continue;
         const isKepler3 = eph.system !== "R" && eph.system !== "S";
-        const dts = dtsClock - (tgd && isKepler3 ? eph.tgd : 0);
+        let dts = dtsClock - (tgd && isKepler3 ? eph.tgd : 0);
+        if (sbas) {
+          const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+          if (sc) {
+            sat.x += sc.dPos[0];
+            sat.y += sc.dPos[1];
+            sat.z += sc.dPos[2];
+            dts += sc.dClkS;
+          }
+        }
         const travel = Math.hypot(sat.x - x9, sat.y - y2, sat.z - z2) / C_LIGHT;
         const [sx, sy, sz] = sagnac2(sat, travel);
         const rho = Math.hypot(sx - x9, sy - y2, sz - z2);
@@ -15062,9 +15075,21 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         }
         const tropo = troposphere && positionSane ? tropoDelay2(elev, rxLat, rxHgt) : 0;
         let ionoM = 0;
-        if ((gim || iono) && positionSane) {
+        if ((sbas || gim || iono) && positionSane) {
           let l1M = null;
-          if (gim)
+          if (sbas) {
+            const si = sbas.ionoDelay(
+              gpsWeek,
+              gpsTow,
+              rxLat,
+              rxLon,
+              rxHgt,
+              azim,
+              elev
+            );
+            if (si) l1M = si.delayM;
+          }
+          if (l1M === null && gim)
             l1M = gimSlantIonoDelayL1(gim, rxLat, rxLon, azim, elev, timeMs);
           if (l1M === null && iono) {
             l1M = C_LIGHT * klobucharDelay(iono, rxLat, rxLon, azim, elev, gpsTow);

package/dist/index.d.cts

@@ -17,7 +17,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
-export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.cjs';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.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, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';

package/dist/index.d.ts

@@ -17,7 +17,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
-export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.js';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength } from './positioning.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, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -68,7 +68,7 @@ import {
   parseSbfNav,
   parseSbfReceiverSetup,
   scanSbfFrames
-} from "./chunk-XFPDTGSV.js";
+} from "./chunk-DE5IUULR.js";
 import "./chunk-FA2UTLEQ.js";
 import {
   OEM4_HLEN,
@@ -76,7 +76,7 @@ import {
   oem4Frames,
   parseNovatelNav,
   parseNovatelRange
-} from "./chunk-2IM5JYKP.js";
+} from "./chunk-ZS5RLFWI.js";
 import {
   parseTrimble,
   parseTrimbleNav,
@@ -164,7 +164,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-VYC2FRDV.js";
+} from "./chunk-E3QCNJL4.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,

package/dist/novatel.cjs

@@ -729,7 +729,7 @@ function decodeGpsQzssEphemeris(view, p, sys) {
     tgd: view.getFloat64(p + 172, true)
   };
 }
-function parseNovatelNav(data) {
+function parseNovatelNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const stats = { badCrc: 0 };
   const ephemerides = [];
@@ -802,6 +802,7 @@ function parseNovatelNav(data) {
       const prn = view.getUint32(p + 4, true);
       const msg = new Uint8Array(32);
       msg.set(data.subarray(p + 12, p + 41));
+      opts.onSbasMessage?.(msg, prn, frame.week, frame.towMs / 1e3);
       const eph = decodeSbasGeoNav(msg, prn, frame.week, frame.towMs / 1e3);
       if (!eph) continue;
       const prev = lastGlo.get(eph.prn);

package/dist/novatel.d.cts

@@ -105,7 +105,13 @@ interface NovatelNavResult {
  * Galileo by IODNav + toe + toc per data source, BeiDou by toe + toc,
  * GLONASS by reference time and health.
  */
-declare function parseNovatelNav(data: Uint8Array): NovatelNavResult;
+declare function parseNovatelNav(data: Uint8Array, opts?: {
+    /** Called for each SBAS L1 message (all types) from RAWSBASFRAME. The
+     *  frame carries 232 bits — the correction fields fit, but the inner SBAS
+     *  CRC does not; the outer OEM4 CRC-32 guarantees the bytes. Feed an
+     *  {@link ../positioning/sbas.SbasProcessor}. */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): NovatelNavResult;
 
 /**
  * NovAtel OEM4/6/7 raw-measurement decoding (RANGE log, message 43).

package/dist/novatel.d.ts

@@ -105,7 +105,13 @@ interface NovatelNavResult {
  * Galileo by IODNav + toe + toc per data source, BeiDou by toe + toc,
  * GLONASS by reference time and health.
  */
-declare function parseNovatelNav(data: Uint8Array): NovatelNavResult;
+declare function parseNovatelNav(data: Uint8Array, opts?: {
+    /** Called for each SBAS L1 message (all types) from RAWSBASFRAME. The
+     *  frame carries 232 bits — the correction fields fit, but the inner SBAS
+     *  CRC does not; the outer OEM4 CRC-32 guarantees the bytes. Feed an
+     *  {@link ../positioning/sbas.SbasProcessor}. */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): NovatelNavResult;
 
 /**
  * NovAtel OEM4/6/7 raw-measurement decoding (RANGE log, message 43).

package/dist/novatel.js

@@ -4,7 +4,7 @@ import {
   oem4Frames,
   parseNovatelNav,
   parseNovatelRange
-} from "./chunk-2IM5JYKP.js";
+} from "./chunk-ZS5RLFWI.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-YUW33VXK.js";
 import "./chunk-V2D5KUKH.js";

package/dist/positioning.cjs

@@ -3619,9 +3619,11 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
     convergenceM = 1e-4,
     tgd = true,
     iono,
-    gim
+    gim,
+    sbas
   } = opts;
   const gpsTow = (timeMs - GPS_EPOCH_MS_SPP) / 1e3 % 604800;
+  const gpsWeek = Math.floor((timeMs - GPS_EPOCH_MS_SPP) / 1e3 / 604800);
   const all = [...pseudoranges.keys()].filter((p) => ephemerides.has(p));
   const minSats = (list) => 3 + new Set(list.map((p) => p[0])).size;
   if (all.length < minSats(all)) return null;
@@ -3658,7 +3660,16 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         const sat = computeSatPosition(eph, tTx - dtsClock * 1e3);
         if (!Number.isFinite(sat.x)) continue;
         const isKepler = eph.system !== "R" && eph.system !== "S";
-        const dts = dtsClock - (tgd && isKepler ? eph.tgd : 0);
+        let dts = dtsClock - (tgd && isKepler ? eph.tgd : 0);
+        if (sbas) {
+          const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+          if (sc) {
+            sat.x += sc.dPos[0];
+            sat.y += sc.dPos[1];
+            sat.z += sc.dPos[2];
+            dts += sc.dClkS;
+          }
+        }
         const travel = Math.hypot(sat.x - x9, sat.y - y2, sat.z - z2) / C_LIGHT;
         const [sx, sy, sz] = sagnac2(sat, travel);
         const rho = Math.hypot(sx - x9, sy - y2, sz - z2);
@@ -3679,9 +3690,21 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
         }
         const tropo = troposphere && positionSane ? tropoDelay2(elev, rxLat, rxHgt) : 0;
         let ionoM = 0;
-        if ((gim || iono) && positionSane) {
+        if ((sbas || gim || iono) && positionSane) {
           let l1M = null;
-          if (gim)
+          if (sbas) {
+            const si = sbas.ionoDelay(
+              gpsWeek,
+              gpsTow,
+              rxLat,
+              rxLon,
+              rxHgt,
+              azim,
+              elev
+            );
+            if (si) l1M = si.delayM;
+          }
+          if (l1M === null && gim)
             l1M = gimSlantIonoDelayL1(gim, rxLat, rxLon, azim, elev, timeMs);
           if (l1M === null && iono) {
             l1M = C_LIGHT * klobucharDelay(iono, rxLat, rxLon, azim, elev, gpsTow);

package/dist/positioning.d.cts

@@ -1223,6 +1223,28 @@ interface SppOptions {
      * reference the broadcast clock already matches. Default true.
      */
     tgd?: boolean;
+    /**
+     * SBAS wide-area corrections (WAAS/EGNOS/…). When given, each satellite in
+     * the SBAS PRN mask has its position + clock adjusted by the fast/long-term
+     * corrections, and — where the ionospheric pierce point is covered by the
+     * SBAS grid — the SBAS slant ionosphere replaces the Klobuchar/GIM model.
+     * Pass an {@link ./sbas.SbasProcessor} (or anything with these two methods).
+     * Combine with single-frequency pseudoranges, not iono-free input.
+     */
+    sbas?: SbasCorrectionSource;
+}
+/** The two SBAS corrections `solveSpp` consumes — {@link ./sbas.SbasProcessor}
+ *  satisfies this structurally. */
+interface SbasCorrectionSource {
+    satCorrection(prn: string, week: number, tow: number): {
+        dPos: [number, number, number];
+        dClkS: number;
+        varM2: number;
+    } | null;
+    ionoDelay(week: number, tow: number, latRad: number, lonRad: number, heightM: number, azRad: number, elRad: number): {
+        delayM: number;
+        varM2: number;
+    } | null;
 }
 interface SppSolution {
     /** Receiver position, ECEF meters. */
@@ -1263,4 +1285,4 @@ declare function ionoFree(p1: number, p2: number, f1: number, f2: number): numbe
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasCorrectionSource, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -1223,6 +1223,28 @@ interface SppOptions {
      * reference the broadcast clock already matches. Default true.
      */
     tgd?: boolean;
+    /**
+     * SBAS wide-area corrections (WAAS/EGNOS/…). When given, each satellite in
+     * the SBAS PRN mask has its position + clock adjusted by the fast/long-term
+     * corrections, and — where the ionospheric pierce point is covered by the
+     * SBAS grid — the SBAS slant ionosphere replaces the Klobuchar/GIM model.
+     * Pass an {@link ./sbas.SbasProcessor} (or anything with these two methods).
+     * Combine with single-frequency pseudoranges, not iono-free input.
+     */
+    sbas?: SbasCorrectionSource;
+}
+/** The two SBAS corrections `solveSpp` consumes — {@link ./sbas.SbasProcessor}
+ *  satisfies this structurally. */
+interface SbasCorrectionSource {
+    satCorrection(prn: string, week: number, tow: number): {
+        dPos: [number, number, number];
+        dClkS: number;
+        varM2: number;
+    } | null;
+    ionoDelay(week: number, tow: number, latRad: number, lonRad: number, heightM: number, azRad: number, elRad: number): {
+        delayM: number;
+        varM2: number;
+    } | null;
 }
 interface SppSolution {
     /** Receiver position, ECEF meters. */
@@ -1263,4 +1285,4 @@ declare function ionoFree(p1: number, p2: number, f1: number, f2: number): numbe
  */
 declare function solveSpp(pseudoranges: Map<string, number>, ephemerides: Map<string, Ephemeris>, timeMs: number, opts?: SppOptions): SppSolution | null;
 
-export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type DgnssOptions, type DgnssSolution, type EphemerisSource, type FixedPositionResult, IONO_L1_M_PER_TECU, type KlobucharCoeffs, type LambdaResult, type NetworkFcbResult, type NetworkFcbSummary, type NlArc, type NlFcbResult, type PppAntennaModel, type PppArc, type PppCorrectionConfig, type PppCorrections, type PppEpoch, type PppEpochContext, type PppEpochResult, type PppFixState, type PppOptions, type PppSatObs, type PppSolution, type RawObservation, type RtkEpochMeasurements, RtkFloatEngine, type RtkFloatOptions, type RtkFloatSolution, type RtkMeasurement, type RtkObsEpoch, type RtkObsInput, type RtkPostProcessOptions, type RtkPostProcessResult, type RtkTrackPoint, type SatGeom, type SbasCorrectionSource, type SbasIonoDelay, SbasProcessor, type SbasSatCorrection, type SppOptions, type SppSolution, type TropoMapping, type WlArc, type WlFcbOptions, type WlFcbResult, type WlObs, buildPppAntenna, createPppCorrections, estimateNarrowlaneFcb, estimateNetworkFcbs, estimateWidelaneFcb, extractWidelaneArcs, fixPppPosition, gimSlantIonoDelayL1, gimVerticalTec, gmst, ionoFree, klobucharDelay, lambdaReduction, lambdaSearch, moonEcef, niellMapping, nlWavelength, postProcessRtk, resolvePppAmbiguities, satClockCorrection, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -29,7 +29,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-VYC2FRDV.js";
+} from "./chunk-E3QCNJL4.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-MIIM4LDY.js";

package/dist/sbf.cjs

@@ -1981,7 +1981,7 @@ function decodeSbasGeoNav(msg, prn, week, tow) {
 }
 
 // src/sbf/rawnav-sbas.ts
-function feedGeoBlock(view, b) {
+function feedGeoBlock(view, b, onSbasMessage) {
   const tow = view.getUint32(b + 8, true);
   const wnc = view.getUint16(b + 12, true);
   const svid = view.getUint8(b + 14);
@@ -1997,10 +1997,11 @@ function feedGeoBlock(view, b) {
   }
   if (!sbasCrcOk(msg)) return { badCrc: true };
   const prn = Number(prnStr.slice(1)) + 100;
+  onSbasMessage?.(msg, prn, wnc, tow / 1e3);
   const eph = decodeSbasGeoNav(msg, prn, wnc, tow / 1e3);
   return eph ? { eph } : {};
 }
-function parseSbfGeoNav(data) {
+function parseSbfGeoNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -2009,7 +2010,7 @@ function parseSbfGeoNav(data) {
   scanSbfFrames(data, view, (id, b, len) => {
     if (id !== 4020 || len < 52) return;
     messages++;
-    const r = feedGeoBlock(view, b);
+    const r = feedGeoBlock(view, b, opts.onSbasMessage);
     if (r.eph) {
       const key = `${r.eph.prn}|${r.eph.tocDate.getTime()}`;
       if (!seen.has(key)) {
@@ -2029,7 +2030,7 @@ function ephKey(e) {
   const source = "source" in e ? e.source : "";
   return `${e.prn}|${e.tocDate.getTime()}|${source}`;
 }
-function decodeSbfNavigation(data) {
+function decodeSbfNavigation(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
   const seen = /* @__PURE__ */ new Set();
@@ -2150,7 +2151,7 @@ function decodeSbfNavigation(data) {
       case 4020: {
         if (len < 52) return;
         counts.sbasRaw++;
-        const r = feedGeoBlock(view, b);
+        const r = feedGeoBlock(view, b, opts.onSbasMessage);
         if (r.eph) addEph(r.eph);
         else if (r.badCrc) counts.badFrames++;
         return;

package/dist/sbf.d.cts

@@ -265,7 +265,11 @@ interface SbfNavigation {
  * F/NAV, GLONASS strings, BeiDou D1/D2, and the Klobuchar/NeQuick iono +
  * GPS-UTC leap-second blocks. See the module header for merge semantics.
  */
-declare function decodeSbfNavigation(data: Uint8Array): SbfNavigation;
+declare function decodeSbfNavigation(data: Uint8Array, opts?: {
+    /** Called for every CRC-valid SBAS L1 message (all types) — feed an
+     *  {@link ../positioning/sbas.SbasProcessor} to build SBAS corrections. */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): SbfNavigation;
 
 /**
  * Septentrio SBF raw Galileo navigation-page blocks: GALRawINAV (4023)
@@ -372,7 +376,7 @@ declare function parseSbfGloNav(data: Uint8Array): SbfGloNavResult;
  * message type 9. Returns the ephemeris, `{ badCrc: true }` on a failed
  * CRC, or `{}` for any other (non–type-9) message.
  */
-declare function feedGeoBlock(view: DataView, b: number): {
+declare function feedGeoBlock(view: DataView, b: number, onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void): {
     eph?: GlonassEphemeris;
     badCrc?: boolean;
 };
@@ -391,7 +395,9 @@ interface SbfGeoNavResult {
  * the same blocks through {@link feedGeoBlock}; this standalone parser
  * mirrors the per-constellation `parseSbf*` helpers.
  */
-declare function parseSbfGeoNav(data: Uint8Array): SbfGeoNavResult;
+declare function parseSbfGeoNav(data: Uint8Array, opts?: {
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): SbfGeoNavResult;
 
 /** Constellations HAS corrects: GPS (GNSS ID 0), Galileo (GNSS ID 2). */
 type HasSystem = 'G' | 'E';

package/dist/sbf.d.ts

@@ -265,7 +265,11 @@ interface SbfNavigation {
  * F/NAV, GLONASS strings, BeiDou D1/D2, and the Klobuchar/NeQuick iono +
  * GPS-UTC leap-second blocks. See the module header for merge semantics.
  */
-declare function decodeSbfNavigation(data: Uint8Array): SbfNavigation;
+declare function decodeSbfNavigation(data: Uint8Array, opts?: {
+    /** Called for every CRC-valid SBAS L1 message (all types) — feed an
+     *  {@link ../positioning/sbas.SbasProcessor} to build SBAS corrections. */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): SbfNavigation;
 
 /**
  * Septentrio SBF raw Galileo navigation-page blocks: GALRawINAV (4023)
@@ -372,7 +376,7 @@ declare function parseSbfGloNav(data: Uint8Array): SbfGloNavResult;
  * message type 9. Returns the ephemeris, `{ badCrc: true }` on a failed
  * CRC, or `{}` for any other (non–type-9) message.
  */
-declare function feedGeoBlock(view: DataView, b: number): {
+declare function feedGeoBlock(view: DataView, b: number, onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void): {
     eph?: GlonassEphemeris;
     badCrc?: boolean;
 };
@@ -391,7 +395,9 @@ interface SbfGeoNavResult {
  * the same blocks through {@link feedGeoBlock}; this standalone parser
  * mirrors the per-constellation `parseSbf*` helpers.
  */
-declare function parseSbfGeoNav(data: Uint8Array): SbfGeoNavResult;
+declare function parseSbfGeoNav(data: Uint8Array, opts?: {
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
+}): SbfGeoNavResult;
 
 /** Constellations HAS corrects: GPS (GNSS ID 0), Galileo (GNSS ID 2). */
 type HasSystem = 'G' | 'E';

package/dist/sbf.js

@@ -15,7 +15,7 @@ import {
   parseSbfNav,
   parseSbfReceiverSetup,
   scanSbfFrames
-} from "./chunk-XFPDTGSV.js";
+} from "./chunk-DE5IUULR.js";
 import "./chunk-FA2UTLEQ.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-37QNKGTC.js";

package/package.json

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