gnss-js 1.48.0 → 1.49.0

package/dist/positioning.d.cts

@@ -550,6 +550,71 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/** A per-satellite SBAS correction, ready to apply in an SPP solve. */
+interface SbasSatCorrection {
+    /** Long-term ECEF position correction to add to the satellite position (m). */
+    dPos: [number, number, number];
+    /** Clock correction to add to the satellite clock bias (s) — long-term + fast. */
+    dClkS: number;
+    /** Pseudorange fast correction alone (m), for reference. */
+    prcM: number;
+    /** Variance of the correction (m²). */
+    varM2: number;
+    /** Applied long-term correction IODE (matches the broadcast ephemeris IODE). */
+    iode: number;
+}
+/** A slant ionospheric delay at a pierce point from the SBAS grid. */
+interface SbasIonoDelay {
+    /** Slant ionospheric delay on L1 (m). */
+    delayM: number;
+    /** Variance (m²). */
+    varM2: number;
+}
+declare class SbasProcessor {
+    private sats;
+    private satIdx;
+    private iodp;
+    private tlat;
+    private ion;
+    /** GEO PRN(s) whose messages have been ingested. */
+    readonly geoPrns: Set<string>;
+    /**
+     * Ingest one decoded SBAS L1 message (250-bit, MSB-first) received at the
+     * given GPS week / time-of-week. Returns the message type, or -1 if it was
+     * not a recognised/applicable correction message.
+     */
+    update(msg: Uint8Array, week: number, tow: number, geoPrn?: number): number;
+    private decodeMask;
+    private decodeFast;
+    private decodeIntegrity;
+    private decodeDegradation;
+    private decodeMixed;
+    private decodeLongCorrHalf;
+    private decodeLong0;
+    private decodeLong1;
+    private decodeIgpMask;
+    private decodeIgpDelay;
+    /** PRNs currently in the correction mask (in mask order). */
+    activeSats(): string[];
+    /** Number of ionospheric grid points with a valid delay across all bands. */
+    ionoGridPoints(): number;
+    /**
+     * Satellite range/clock correction for a PRN at a GPS week/tow, or null if
+     * no valid (unexpired) long-term + fast correction is available. Add `dPos`
+     * to the satellite ECEF position and `dClkS` to its clock bias.
+     */
+    satCorrection(prn: string, week: number, tow: number): SbasSatCorrection | null;
+    /**
+     * Slant L1 ionospheric delay (m) at the pierce point for a receiver geodetic
+     * position (lat/lon rad, height m) and a satellite az/el (rad), interpolated
+     * from the SBAS grid — or null if the pierce point isn't covered. RTKLIB
+     * `sbsioncorr`.
+     */
+    ionoDelay(week: number, tow: number, latRad: number, lonRad: number, heightM: number, azRad: number, elRad: number): SbasIonoDelay | null;
+    /** Find the four IGPs bracketing a pierce point (RTKLIB `searchigp`). */
+    private searchIgp;
+}
+
 /**
  * ANTEX phase-centre application for PPP: satellite and receiver antenna
  * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
@@ -1198,4 +1263,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 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 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

@@ -550,6 +550,71 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/** A per-satellite SBAS correction, ready to apply in an SPP solve. */
+interface SbasSatCorrection {
+    /** Long-term ECEF position correction to add to the satellite position (m). */
+    dPos: [number, number, number];
+    /** Clock correction to add to the satellite clock bias (s) — long-term + fast. */
+    dClkS: number;
+    /** Pseudorange fast correction alone (m), for reference. */
+    prcM: number;
+    /** Variance of the correction (m²). */
+    varM2: number;
+    /** Applied long-term correction IODE (matches the broadcast ephemeris IODE). */
+    iode: number;
+}
+/** A slant ionospheric delay at a pierce point from the SBAS grid. */
+interface SbasIonoDelay {
+    /** Slant ionospheric delay on L1 (m). */
+    delayM: number;
+    /** Variance (m²). */
+    varM2: number;
+}
+declare class SbasProcessor {
+    private sats;
+    private satIdx;
+    private iodp;
+    private tlat;
+    private ion;
+    /** GEO PRN(s) whose messages have been ingested. */
+    readonly geoPrns: Set<string>;
+    /**
+     * Ingest one decoded SBAS L1 message (250-bit, MSB-first) received at the
+     * given GPS week / time-of-week. Returns the message type, or -1 if it was
+     * not a recognised/applicable correction message.
+     */
+    update(msg: Uint8Array, week: number, tow: number, geoPrn?: number): number;
+    private decodeMask;
+    private decodeFast;
+    private decodeIntegrity;
+    private decodeDegradation;
+    private decodeMixed;
+    private decodeLongCorrHalf;
+    private decodeLong0;
+    private decodeLong1;
+    private decodeIgpMask;
+    private decodeIgpDelay;
+    /** PRNs currently in the correction mask (in mask order). */
+    activeSats(): string[];
+    /** Number of ionospheric grid points with a valid delay across all bands. */
+    ionoGridPoints(): number;
+    /**
+     * Satellite range/clock correction for a PRN at a GPS week/tow, or null if
+     * no valid (unexpired) long-term + fast correction is available. Add `dPos`
+     * to the satellite ECEF position and `dClkS` to its clock bias.
+     */
+    satCorrection(prn: string, week: number, tow: number): SbasSatCorrection | null;
+    /**
+     * Slant L1 ionospheric delay (m) at the pierce point for a receiver geodetic
+     * position (lat/lon rad, height m) and a satellite az/el (rad), interpolated
+     * from the SBAS grid — or null if the pierce point isn't covered. RTKLIB
+     * `sbsioncorr`.
+     */
+    ionoDelay(week: number, tow: number, latRad: number, lonRad: number, heightM: number, azRad: number, elRad: number): SbasIonoDelay | null;
+    /** Find the four IGPs bracketing a pierce point (RTKLIB `searchigp`). */
+    private searchIgp;
+}
+
 /**
  * ANTEX phase-centre application for PPP: satellite and receiver antenna
  * PCO (offset) + PCV (elevation-dependent variation), ionosphere-free
@@ -1198,4 +1263,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 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 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

@@ -1,6 +1,7 @@
 import {
   IONO_L1_M_PER_TECU,
   RtkFloatEngine,
+  SbasProcessor,
   buildPppAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,
@@ -28,18 +29,22 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-KRXX4UGR.js";
+} from "./chunk-DWLVCXA6.js";
+import "./chunk-WQX354KD.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-MIIM4LDY.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-GPEP67TG.js";
+import "./chunk-YUW33VXK.js";
+import "./chunk-V2D5KUKH.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";
 export {
   IONO_L1_M_PER_TECU,
   RtkFloatEngine,
+  SbasProcessor,
   buildPppAntenna,
   createPppCorrections,
   estimateNarrowlaneFcb,

package/dist/rinex.js

@@ -21,8 +21,6 @@ import {
   writeRinexNav4,
   writeRinexObsBlob
 } from "./chunk-6SWJTN6H.js";
-import "./chunk-YUW33VXK.js";
-import "./chunk-V2D5KUKH.js";
 import {
   SYSTEM_ORDER,
   crxDecompress,
@@ -38,6 +36,8 @@ import {
   parseSp3,
   sp3Position
 } from "./chunk-GPEP67TG.js";
+import "./chunk-YUW33VXK.js";
+import "./chunk-V2D5KUKH.js";
 import "./chunk-FIEWO4J4.js";
 export {
   EMPTY_WARNINGS,

package/dist/sbf.js

@@ -15,9 +15,9 @@ import {
   parseSbfNav,
   parseSbfReceiverSetup,
   scanSbfFrames
-} from "./chunk-MEDAWNWM.js";
+} from "./chunk-XFPDTGSV.js";
 import "./chunk-FA2UTLEQ.js";
-import "./chunk-ZIGMN3E5.js";
+import "./chunk-WQX354KD.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";
 import "./chunk-YUW33VXK.js";

package/dist/ubx.cjs

@@ -1409,6 +1409,7 @@ function parseUbxRawNav(data, opts = {}) {
         continue;
       }
       if (refWeek === 0) continue;
+      opts.onSbasMessage?.(msg, svId, refWeek, refTow);
       const eph = decodeSbasGeoNav(msg, svId, refWeek, refTow);
       if (eph) {
         const key = `${eph.prn}@${eph.tocDate.getTime()}`;

package/dist/ubx.d.cts

@@ -228,6 +228,13 @@ interface UbxRawNavOptions {
      * week and are decoded regardless); the system clock is never used.
      */
     refWeek?: number;
+    /**
+     * Called for every CRC-valid SBAS L1 C/A message (all types, not only the
+     * GEO navigation MT9), with the 250-bit message MSB-first, the GEO PRN, and
+     * the reception GPS week / time-of-week. Feed these to an
+     * {@link ../positioning/sbas.SbasProcessor} to build SBAS corrections.
+     */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
 }
 interface UbxRawNavResult {
     /**

package/dist/ubx.d.ts

@@ -228,6 +228,13 @@ interface UbxRawNavOptions {
      * week and are decoded regardless); the system clock is never used.
      */
     refWeek?: number;
+    /**
+     * Called for every CRC-valid SBAS L1 C/A message (all types, not only the
+     * GEO navigation MT9), with the 250-bit message MSB-first, the GEO PRN, and
+     * the reception GPS week / time-of-week. Feed these to an
+     * {@link ../positioning/sbas.SbasProcessor} to build SBAS corrections.
+     */
+    onSbasMessage?: (msg: Uint8Array, prn: number, week: number, tow: number) => void;
 }
 interface UbxRawNavResult {
     /**

package/dist/ubx.js

@@ -5,9 +5,9 @@ import {
   parseUbxRawNav,
   parseUbxRawx,
   ubxFrames
-} from "./chunk-EQNC24OB.js";
+} from "./chunk-DEVUKH6A.js";
 import "./chunk-FA2UTLEQ.js";
-import "./chunk-ZIGMN3E5.js";
+import "./chunk-WQX354KD.js";
 import "./chunk-YUW33VXK.js";
 import "./chunk-V2D5KUKH.js";
 import "./chunk-HVXYFUCB.js";

package/package.json

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