gnss-js 1.56.0 → 1.58.0

package/dist/positioning.cjs

@@ -43,6 +43,8 @@ __export(positioning_exports, {
   postProcessRtk: () => postProcessRtk,
   resolvePppAmbiguities: () => resolvePppAmbiguities,
   satClockCorrection: () => satClockCorrection,
+  sbasIonoDeg: () => sbasIonoDeg,
+  sbasLongTermDeg: () => sbasLongTermDeg,
   sbasProtectionLevels: () => sbasProtectionLevels,
   solidEarthTide: () => solidEarthTide,
   solveDgnss: () => solveDgnss,
@@ -1744,6 +1746,17 @@ var DEGF = [
   58e-4
 ];
 var degfcorr = (ai) => ai > 0 && ai <= 15 ? DEGF[ai] : 58e-4;
+function sbasLongTermDeg(d, vel, ageSec) {
+  if (vel) {
+    if (ageSec > 0 && ageSec < d.iltcV1) return 0;
+    return d.cltcLsb + d.cltcV1 * Math.max(0, -ageSec, ageSec - d.iltcV1);
+  }
+  return d.cltcV0 * Math.floor(Math.abs(ageSec) / d.iltcV0);
+}
+function sbasIonoDeg(d, ageSec) {
+  const dt = Math.abs(ageSec);
+  return d.cionoStep * Math.floor(dt / d.iiono) + d.cionoRamp * dt;
+}
 var x1 = [-75, -65, -55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 65, 75, 85];
 var x2 = [-55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55];
 var x3 = [-75, -65, -55, -50, -45, -40, -35, -30, -25, -20, -15, -10, -5, 0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 65, 75];
@@ -1797,6 +1810,7 @@ var SbasProcessor = class {
   satIdx = /* @__PURE__ */ new Map();
   iodp = -1;
   tlat = 0;
+  degr = null;
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -1823,6 +1837,8 @@ var SbasProcessor = class {
         return this.decodeIntegrity(msg) ? type : -1;
       case 7:
         return this.decodeDegradation(msg) ? type : -1;
+      case 10:
+        return this.decodeMt10(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -1891,6 +1907,30 @@ var SbasProcessor = class {
     }
     return true;
   }
+  /** MT10 degradation factors (DO-229D Table A-9). Fields follow the 8-bit
+   *  preamble + 6-bit type (start bit 14), MSB-first, with the listed LSBs.
+   *  Iiono / Iltc_v0 of 0 must be read as 1 (Table A-9, Note 3). */
+  decodeMt10(msg) {
+    const u = (p, n) => getBitU(msg, p, n);
+    this.degr = {
+      brrc: u(14, 10) * 2e-3,
+      cltcLsb: u(24, 10) * 2e-3,
+      cltcV1: u(34, 10) * 5e-5,
+      iltcV1: u(44, 9),
+      cltcV0: u(53, 10) * 2e-3,
+      iltcV0: u(63, 9) || 1,
+      cgeoLsb: u(72, 10) * 5e-4,
+      cgeoV: u(82, 10) * 5e-5,
+      igeo: u(92, 9),
+      cer: u(101, 6) * 0.5,
+      cionoStep: u(107, 10) * 1e-3,
+      iiono: u(117, 9) || 1,
+      cionoRamp: u(126, 10) * 5e-6,
+      rssUdre: u(136, 1) === 1,
+      rssIono: u(137, 1) === 1
+    };
+    return true;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -1937,7 +1977,8 @@ var SbasProcessor = class {
       dvel: [0, 0, 0],
       daf0: getBitS(msg, p + 41, 10) * P2_31,
       daf1: 0,
-      t0s: gpsSeconds(week, tow)
+      t0s: gpsSeconds(week, tow),
+      vel: false
     };
     return true;
   }
@@ -1961,7 +2002,8 @@ var SbasProcessor = class {
       ],
       daf0: getBitS(msg, p + 47, 11) * P2_31,
       daf1: getBitS(msg, p + 82, 8) * P2_39,
-      t0s: gpsSeconds(week, tow + t)
+      t0s: gpsSeconds(week, tow + t),
+      vel: true
     };
     return true;
   }
@@ -2082,9 +2124,16 @@ var SbasProcessor = class {
     const tf = now - s.fcorr.t0s + this.tlat;
     if (Math.abs(tf) > MAXSBSAGEF || s.fcorr.udre >= 15) return null;
     const prc = s.fcorr.prc;
-    const varM2 = varfcorr(s.fcorr.udre) + degfcorr(s.fcorr.ai) * tf * tf / 2;
+    const sigUdre = Math.sqrt(varfcorr(s.fcorr.udre));
+    const efc = degfcorr(s.fcorr.ai) * tf * tf / 2;
+    const eltc = this.degr && s.lcorr && s.lcorr.t0s !== 0 ? sbasLongTermDeg(this.degr, s.lcorr.vel, now - s.lcorr.t0s) : 0;
+    const varM2 = this.degr?.rssUdre ? sigUdre * sigUdre + efc * efc + eltc * eltc : (sigUdre + efc + eltc) ** 2;
     return { dPos: drs, dClkS: dclk + prc / C_LIGHT, prcM: prc, varM2, iode };
   }
+  /** The last-decoded MT10 degradation factors, or null if none seen. */
+  get degradation() {
+    return this.degr;
+  }
   /**
    * 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
@@ -2133,14 +2182,18 @@ var SbasProcessor = class {
     } else return null;
     const now = gpsSeconds(week, tow);
     let delay = 0;
-    let varM2 = 0;
+    let varVert = 0;
     for (let i = 0; i < 4; i++) {
       const g = pts[i];
       if (!g) continue;
       delay += w[i] * g.delay;
-      varM2 += w[i] * varicorr(g.give) * 9e-8 * Math.abs(now - g.t0s);
+      const sigGive = Math.sqrt(varicorr(g.give));
+      const d = this.degr;
+      const eiono = d ? sbasIonoDeg(d, now - g.t0s) : 0;
+      const sig2 = d?.rssIono ? sigGive * sigGive + eiono * eiono : (sigGive + eiono) ** 2;
+      varVert += w[i] * sig2;
     }
-    return { delayM: delay * fp, varM2: varM2 * fp * fp };
+    return { delayM: delay * fp, varM2: varVert * fp * fp };
   }
   /** Find the four IGPs bracketing a pierce point (RTKLIB `searchigp`). */
   searchIgp(latP, lonP) {
@@ -3627,10 +3680,10 @@ function sagnac2(pos, travelTimeS) {
 function sbasAirTropoVar(elRad) {
   const elDeg = elRad * 180 / Math.PI;
   const sigMp = 0.13 + 0.53 * Math.exp(-elDeg / 10);
-  const sigNoise = 0.15 + 0.43 * Math.exp(-elDeg / 6.9);
+  const NOISE_DIVG_AADA = 0.36;
   const m = 1.001 / Math.sqrt(2001e-6 + Math.sin(elRad) * Math.sin(elRad));
   const sigTropo = 0.12 * m;
-  return sigMp * sigMp + sigNoise * sigNoise + sigTropo * sigTropo;
+  return sigMp * sigMp + NOISE_DIVG_AADA * NOISE_DIVG_AADA + sigTropo * sigTropo;
 }
 function invertMatrix(A) {
   const n = A.length;
@@ -3939,6 +3992,8 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
   postProcessRtk,
   resolvePppAmbiguities,
   satClockCorrection,
+  sbasIonoDeg,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   solidEarthTide,
   solveDgnss,

package/dist/positioning.d.cts

@@ -550,6 +550,35 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/**
+ * Long-term correction degradation ε_ltc (m) — DO-229D §A.4.5.1.3, eqs A-54
+ * (velocity code 1) and A-55 (velocity code 0). `ageSec` = t − t₀ (s).
+ */
+declare function sbasLongTermDeg(d: Degradation, vel: boolean, ageSec: number): number;
+/**
+ * Ionospheric-correction degradation ε_iono (m) — DO-229D §A.4.5.2, eq A-59.
+ * `ageSec` = t − t_iono (s) of the grid point.
+ */
+declare function sbasIonoDeg(d: Degradation, ageSec: number): number;
+/** MT10 degradation parameters (DO-229D Table A-9), needed for the residual
+ *  variances σ²_flt / σ²_ionogrid (§A.4.5, §J.2.2). */
+interface Degradation {
+    brrc: number;
+    cltcLsb: number;
+    cltcV1: number;
+    iltcV1: number;
+    cltcV0: number;
+    iltcV0: number;
+    cgeoLsb: number;
+    cgeoV: number;
+    igeo: number;
+    cer: number;
+    cionoStep: number;
+    iiono: number;
+    cionoRamp: number;
+    rssUdre: boolean;
+    rssIono: boolean;
+}
 /** 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). */
@@ -575,6 +604,7 @@ declare class SbasProcessor {
     private satIdx;
     private iodp;
     private tlat;
+    private degr;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -588,6 +618,10 @@ declare class SbasProcessor {
     private decodeFast;
     private decodeIntegrity;
     private decodeDegradation;
+    /** MT10 degradation factors (DO-229D Table A-9). Fields follow the 8-bit
+     *  preamble + 6-bit type (start bit 14), MSB-first, with the listed LSBs.
+     *  Iiono / Iltc_v0 of 0 must be read as 1 (Table A-9, Note 3). */
+    private decodeMt10;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -624,6 +658,8 @@ declare class SbasProcessor {
      * to the satellite ECEF position and `dClkS` to its clock bias.
      */
     satCorrection(prn: string, week: number, tow: number): SbasSatCorrection | null;
+    /** The last-decoded MT10 degradation factors, or null if none seen. */
+    get degradation(): Readonly<Degradation> | 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
@@ -1334,4 +1370,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 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 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, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, 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, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.d.ts

@@ -550,6 +550,35 @@ declare function lambdaSearch(a: Float64Array, Q: Float64Array, n: number, m?: n
  */
 declare function lambdaReduction(Q: Float64Array, n: number): Float64Array | null;
 
+/**
+ * Long-term correction degradation ε_ltc (m) — DO-229D §A.4.5.1.3, eqs A-54
+ * (velocity code 1) and A-55 (velocity code 0). `ageSec` = t − t₀ (s).
+ */
+declare function sbasLongTermDeg(d: Degradation, vel: boolean, ageSec: number): number;
+/**
+ * Ionospheric-correction degradation ε_iono (m) — DO-229D §A.4.5.2, eq A-59.
+ * `ageSec` = t − t_iono (s) of the grid point.
+ */
+declare function sbasIonoDeg(d: Degradation, ageSec: number): number;
+/** MT10 degradation parameters (DO-229D Table A-9), needed for the residual
+ *  variances σ²_flt / σ²_ionogrid (§A.4.5, §J.2.2). */
+interface Degradation {
+    brrc: number;
+    cltcLsb: number;
+    cltcV1: number;
+    iltcV1: number;
+    cltcV0: number;
+    iltcV0: number;
+    cgeoLsb: number;
+    cgeoV: number;
+    igeo: number;
+    cer: number;
+    cionoStep: number;
+    iiono: number;
+    cionoRamp: number;
+    rssUdre: boolean;
+    rssIono: boolean;
+}
 /** 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). */
@@ -575,6 +604,7 @@ declare class SbasProcessor {
     private satIdx;
     private iodp;
     private tlat;
+    private degr;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -588,6 +618,10 @@ declare class SbasProcessor {
     private decodeFast;
     private decodeIntegrity;
     private decodeDegradation;
+    /** MT10 degradation factors (DO-229D Table A-9). Fields follow the 8-bit
+     *  preamble + 6-bit type (start bit 14), MSB-first, with the listed LSBs.
+     *  Iiono / Iltc_v0 of 0 must be read as 1 (Table A-9, Note 3). */
+    private decodeMt10;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -624,6 +658,8 @@ declare class SbasProcessor {
      * to the satellite ECEF position and `dClkS` to its clock bias.
      */
     satCorrection(prn: string, week: number, tow: number): SbasSatCorrection | null;
+    /** The last-decoded MT10 degradation factors, or null if none seen. */
+    get degradation(): Readonly<Degradation> | 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
@@ -1334,4 +1370,4 @@ declare function sbasProtectionLevels(sats: readonly {
  */
 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 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, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };
+export { type AntennaOffset, type ArFixedSat, type ArOptions, type ArResult, type ArSat, type CorrectionResult, type Degradation, 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, sbasIonoDeg, sbasLongTermDeg, sbasProtectionLevels, solidEarthTide, solveDgnss, solvePpp, solveSpp, sunEcef, toRtkEpoch, wlWavelength };

package/dist/positioning.js

@@ -22,6 +22,8 @@ import {
   postProcessRtk,
   resolvePppAmbiguities,
   satClockCorrection,
+  sbasIonoDeg,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   solidEarthTide,
   solveDgnss,
@@ -30,7 +32,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-BU6YC6VK.js";
+} from "./chunk-SWVIBUHS.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-MIIM4LDY.js";
@@ -66,6 +68,8 @@ export {
   postProcessRtk,
   resolvePppAmbiguities,
   satClockCorrection,
+  sbasIonoDeg,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   solidEarthTide,
   solveDgnss,

package/dist/sbf.cjs

@@ -20,12 +20,18 @@ var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: tru
 // src/sbf/index.ts
 var sbf_exports = {};
 __export(sbf_exports, {
+  SBAS_L5_MT_NAMES: () => SBAS_L5_MT_NAMES,
+  SBAS_L5_PREAMBLE_UW: () => SBAS_L5_PREAMBLE_UW,
   decodeSbfNavigation: () => decodeSbfNavigation,
   feedGeoBlock: () => feedGeoBlock,
+  feedGeoL5Block: () => feedGeoL5Block,
+  isSbasL1Preamble: () => isSbasL1Preamble,
+  isSbasL5Correction: () => isSbasL5Correction,
   parseSbfAlmanac: () => parseSbfAlmanac,
   parseSbfBdsNav: () => parseSbfBdsNav,
   parseSbfCnav: () => parseSbfCnav,
   parseSbfGalNav: () => parseSbfGalNav,
+  parseSbfGeoL5: () => parseSbfGeoL5,
   parseSbfGeoNav: () => parseSbfGeoNav,
   parseSbfGloNav: () => parseSbfGloNav,
   parseSbfGpsNav: () => parseSbfGpsNav,
@@ -35,6 +41,8 @@ __export(sbf_exports, {
   parseSbfNav: () => parseSbfNav,
   parseSbfPvt: () => parseSbfPvt,
   parseSbfReceiverSetup: () => parseSbfReceiverSetup,
+  sbasL5MessageName: () => sbasL5MessageName,
+  sbasL5MessageType: () => sbasL5MessageType,
   sbfCrc16: () => crc16,
   scanSbfFrames: () => scanSbfFrames
 });
@@ -2033,13 +2041,42 @@ function decodeSbasGeoNav(msg, prn, week, tow) {
   };
 }
 
+// src/navbits/sbas-l5.ts
+var SBAS_L5_PREAMBLE_UW = 6056250;
+function isSbasL1Preamble(msg) {
+  return msg[0] === 83 || msg[0] === 154 || msg[0] === 198;
+}
+function sbasL5MessageType(msg) {
+  return getBitU(msg, 4, 6);
+}
+var SBAS_L5_MT_NAMES = {
+  0: "Do Not Use (test)",
+  31: "SBAS satellite mask",
+  32: "Satellite clock-ephemeris corrections + covariance",
+  34: "Integrity (DFREI/DFRECI)",
+  35: "Integrity (DFREI/DFRECI)",
+  36: "Integrity (DFREI/DFRECI)",
+  37: "Degradation parameters + DFREI scale table",
+  39: "SBAS satellite clock/ephemeris/covariance-1",
+  40: "SBAS satellite clock/ephemeris/covariance-2",
+  42: "SNT-to-UTC offset",
+  47: "SBAS satellite almanacs",
+  62: "Reserved",
+  63: "Null message"
+};
+function sbasL5MessageName(type) {
+  return SBAS_L5_MT_NAMES[type] ?? (type >= 1 && type <= 61 ? "Spare" : `Type ${type}`);
+}
+function isSbasL5Correction(type) {
+  return type === 31 || type === 32 || type >= 34 && type <= 37 || type === 39 || type === 40 || type === 42 || type === 47;
+}
+
 // src/sbf/rawnav-sbas.ts
-function feedGeoBlock(view, b, onSbasMessage) {
+function readGeoMessage(view, b) {
   const tow = view.getUint32(b + 8, true);
   const wnc = view.getUint16(b + 12, true);
-  const svid = view.getUint8(b + 14);
-  const prnStr = svidToPrn(svid);
-  if (!prnStr || prnStr[0] !== "S") return {};
+  const prnStr = svidToPrn(view.getUint8(b + 14));
+  if (!prnStr || prnStr[0] !== "S") return null;
   const msg = new Uint8Array(32);
   for (let k = 0; k < 8; k++) {
     const w = view.getUint32(b + 20 + 4 * k, true);
@@ -2048,12 +2085,29 @@ function feedGeoBlock(view, b, onSbasMessage) {
     msg[4 * k + 2] = w >>> 8 & 255;
     msg[4 * k + 3] = w & 255;
   }
-  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 { msg, prn: Number(prnStr.slice(1)) + 100, wnc, towSec: tow / 1e3 };
+}
+function feedGeoBlock(view, b, onSbasMessage) {
+  const g = readGeoMessage(view, b);
+  if (!g) return {};
+  if (!sbasCrcOk(g.msg)) return { badCrc: true };
+  onSbasMessage?.(g.msg, g.prn, g.wnc, g.towSec);
+  const eph = decodeSbasGeoNav(g.msg, g.prn, g.wnc, g.towSec);
   return eph ? { eph } : {};
 }
+function feedGeoL5Block(view, b, onSbasMessage, onDfmcMessage) {
+  const g = readGeoMessage(view, b);
+  if (!g) return {};
+  if (!sbasCrcOk(g.msg)) return { badCrc: true };
+  if (isSbasL1Preamble(g.msg)) {
+    onSbasMessage?.(g.msg, g.prn, g.wnc, g.towSec);
+    const eph = decodeSbasGeoNav(g.msg, g.prn, g.wnc, g.towSec);
+    return eph ? { kind: "l1", eph } : { kind: "l1" };
+  }
+  const dfmcType = sbasL5MessageType(g.msg);
+  onDfmcMessage?.(g.msg, g.prn, g.wnc, g.towSec, dfmcType);
+  return { kind: "l5", dfmcType };
+}
 function parseSbfGeoNav(data, opts = {}) {
   const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
   const ephemerides = [];
@@ -2076,6 +2130,47 @@ function parseSbfGeoNav(data, opts = {}) {
   });
   return { ephemerides, badCrc, messages };
 }
+function parseSbfGeoL5(data, opts = {}) {
+  const view = new DataView(data.buffer, data.byteOffset, data.byteLength);
+  const ephemerides = [];
+  const seen = /* @__PURE__ */ new Set();
+  const byType = {};
+  const dfmcPrns = /* @__PURE__ */ new Set();
+  let legacyMessages = 0;
+  let dfmcMessages = 0;
+  let badCrc = 0;
+  let messages = 0;
+  scanSbfFrames(data, view, (id, b, len) => {
+    if (id !== 4021 || len < 52) return;
+    messages++;
+    const r = feedGeoL5Block(view, b, opts.onSbasMessage, opts.onDfmcMessage);
+    if (r.badCrc) {
+      badCrc++;
+    } else if (r.kind === "l1") {
+      legacyMessages++;
+      if (r.eph) {
+        const key = `${r.eph.prn}|${r.eph.tocDate.getTime()}`;
+        if (!seen.has(key)) {
+          seen.add(key);
+          ephemerides.push(r.eph);
+        }
+      }
+    } else if (r.kind === "l5" && r.dfmcType !== void 0) {
+      dfmcMessages++;
+      byType[r.dfmcType] = (byType[r.dfmcType] ?? 0) + 1;
+      const prn = view.getUint8(b + 14);
+      const s = svidToPrn(prn);
+      if (s) dfmcPrns.add(s);
+    }
+  });
+  return {
+    ephemerides,
+    legacyMessages,
+    dfmc: { byType, prns: [...dfmcPrns].sort(), messages: dfmcMessages },
+    badCrc,
+    messages
+  };
+}
 
 // src/sbf/navigation.ts
 var F4_DNU4 = -2e10;
@@ -2098,8 +2193,11 @@ function decodeSbfNavigation(data, opts = {}) {
     gloRaw: 0,
     bdsRaw: 0,
     sbasRaw: 0,
+    sbasL5Raw: 0,
     badFrames: 0
   };
+  const dfmcByType = {};
+  const dfmcPrns = /* @__PURE__ */ new Set();
   const cnavAsm = newCnavAssemblers();
   const galAsm = newGalAssemblers();
   const gloAsm = new GloStringAssembler();
@@ -2209,6 +2307,25 @@ function decodeSbfNavigation(data, opts = {}) {
         else if (r.badCrc) counts.badFrames++;
         return;
       }
+      case 4021: {
+        if (len < 52) return;
+        counts.sbasL5Raw++;
+        const r = feedGeoL5Block(
+          view,
+          b,
+          opts.onSbasMessage,
+          opts.onDfmcMessage
+        );
+        if (r.badCrc) counts.badFrames++;
+        else if (r.kind === "l1") {
+          if (r.eph) addEph(r.eph);
+        } else if (r.kind === "l5" && r.dfmcType !== void 0) {
+          dfmcByType[r.dfmcType] = (dfmcByType[r.dfmcType] ?? 0) + 1;
+          const s = svidToPrn(view.getUint8(b + 14));
+          if (s) dfmcPrns.add(s);
+        }
+        return;
+      }
       /* ---- iono / UTC ---- */
       case 5893:
       // GPSIon
@@ -2233,7 +2350,12 @@ function decodeSbfNavigation(data, opts = {}) {
         return;
     }
   });
-  return { ephemerides, cnav, ionoCorrections, leapSeconds, counts };
+  const dfmc = {
+    byType: dfmcByType,
+    prns: [...dfmcPrns].sort(),
+    messages: Object.values(dfmcByType).reduce((a, v) => a + v, 0)
+  };
+  return { ephemerides, cnav, ionoCorrections, leapSeconds, dfmc, counts };
 }
 
 // src/navbits/rs255.ts
@@ -3428,12 +3550,18 @@ function parseSbfMeas(data) {
 }
 // Annotate the CommonJS export names for ESM import in node:
 0 && (module.exports = {
+  SBAS_L5_MT_NAMES,
+  SBAS_L5_PREAMBLE_UW,
   decodeSbfNavigation,
   feedGeoBlock,
+  feedGeoL5Block,
+  isSbasL1Preamble,
+  isSbasL5Correction,
   parseSbfAlmanac,
   parseSbfBdsNav,
   parseSbfCnav,
   parseSbfGalNav,
+  parseSbfGeoL5,
   parseSbfGeoNav,
   parseSbfGloNav,
   parseSbfGpsNav,
@@ -3443,6 +3571,8 @@ function parseSbfMeas(data) {
   parseSbfNav,
   parseSbfPvt,
   parseSbfReceiverSetup,
+  sbasL5MessageName,
+  sbasL5MessageType,
   sbfCrc16,
   scanSbfFrames
 });