gnss-js 2.3.0 → 2.4.0

package/dist/{chunk-TYJNN77C.js → chunk-B2XWO4RP.js}

@@ -1649,6 +1649,10 @@ var SbasProcessor = class {
   tlat = 0;
   degr = null;
   mt27 = null;
+  /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+   *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+   *  scale exponent. Keyed by PRN; cleared on a mask change. */
+  mt28 = /* @__PURE__ */ new Map();
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -1679,6 +1683,8 @@ var SbasProcessor = class {
         return this.decodeMt10(msg) ? type : -1;
       case 27:
         return this.decodeMt27(msg) ? type : -1;
+      case 28:
+        return this.decodeMt28(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -1705,6 +1711,7 @@ var SbasProcessor = class {
       return true;
     this.iodp = iodp;
     this.sats = prns.map((prn) => ({ prn }));
+    this.mt28.clear();
     this.satIdx.clear();
     this.sats.forEach((s, k) => {
       if (s.prn) this.satIdx.set(s.prn, k);
@@ -1770,7 +1777,9 @@ var SbasProcessor = class {
       iiono: u(117, 9) || 1,
       cionoRamp: u(126, 10) * 5e-6,
       rssUdre: u(136, 1) === 1,
-      rssIono: u(137, 1) === 1
+      rssIono: u(137, 1) === 1,
+      ccovariance: u(138, 7) * 0.1
+      // Table B-49: 0..12.7 m, 0.1 m
     };
     return true;
   }
@@ -1822,6 +1831,51 @@ var SbasProcessor = class {
     }
     return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
   }
+  /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+   *  two satellites per message, each an upper-triangular Cholesky factor E and a
+   *  scale exponent that together define the relative clock-ephemeris covariance
+   *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+   *  the current MT1 mask, so a mismatched IODP is dropped. */
+  decodeMt28(msg) {
+    if (this.iodp !== getBitU(msg, 14, 2)) return false;
+    for (const base of [16, 121]) {
+      const slot = getBitU(msg, base, 6);
+      const prn = this.sats[slot]?.prn;
+      const scaleExp = getBitU(msg, base + 6, 3);
+      const e = [
+        getBitU(msg, base + 9, 9),
+        getBitU(msg, base + 18, 9),
+        getBitU(msg, base + 27, 9),
+        getBitU(msg, base + 36, 9),
+        getBitS(msg, base + 45, 10),
+        getBitS(msg, base + 55, 10),
+        getBitS(msg, base + 65, 10),
+        getBitS(msg, base + 75, 10),
+        getBitS(msg, base + 85, 10),
+        getBitS(msg, base + 95, 10)
+      ];
+      if (prn && e.some((v) => v !== 0)) this.mt28.set(prn, { scaleExp, e });
+    }
+    return true;
+  }
+  /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+   *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+   *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+   *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+  mt28DeltaUdre(prn, los) {
+    const m = this.mt28.get(prn);
+    if (!m) return null;
+    const sf = 2 ** (m.scaleExp - 5);
+    const [e11, e22, e33, e44, e12, e13, e14, e23, e24, e34] = m.e;
+    const [ix, iy, iz] = los;
+    const v0 = (e11 * ix + e12 * iy + e13 * iz + e14) * sf;
+    const v1 = (e22 * iy + e23 * iz + e24) * sf;
+    const v2 = (e33 * iz + e34) * sf;
+    const v3 = e44 * sf;
+    const iCi = v0 * v0 + v1 * v1 + v2 * v2 + v3 * v3;
+    const ec = (this.degr?.ccovariance ?? 0) * sf;
+    return Math.sqrt(iCi) + ec;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -2026,7 +2080,7 @@ var SbasProcessor = class {
    * no valid (unexpired) long-term + fast correction is available. Add `dPos`
    * to the satellite ECEF position and `dClkS` to its clock bias.
    */
-  satCorrection(prn, week, tow, userLatRad, userLonRad) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad, losUnitEcef) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -2048,7 +2102,8 @@ 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 dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+    const mt28 = losUnitEcef ? this.mt28DeltaUdre(prn, losUnitEcef) : null;
+    const dudre = mt28 != null ? mt28 : userLatRad != null && userLonRad != null ? this.deltaUdre(
       userLatRad * 180 / Math.PI,
       userLonRad * 180 / Math.PI
     ) : 1;
@@ -3768,7 +3823,10 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
+        const dx = sat.x - x, dy = sat.y - y, dz = sat.z - z;
+        const dr = Math.hypot(dx, dy, dz) || 1;
+        const los = [dx / dr, dy / dr, dz / dr];
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon, los);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/index.cjs

@@ -14576,6 +14576,10 @@ var SbasProcessor = class {
   tlat = 0;
   degr = null;
   mt27 = null;
+  /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+   *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+   *  scale exponent. Keyed by PRN; cleared on a mask change. */
+  mt28 = /* @__PURE__ */ new Map();
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -14606,6 +14610,8 @@ var SbasProcessor = class {
         return this.decodeMt10(msg) ? type : -1;
       case 27:
         return this.decodeMt27(msg) ? type : -1;
+      case 28:
+        return this.decodeMt28(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -14632,6 +14638,7 @@ var SbasProcessor = class {
       return true;
     this.iodp = iodp;
     this.sats = prns.map((prn) => ({ prn }));
+    this.mt28.clear();
     this.satIdx.clear();
     this.sats.forEach((s, k) => {
       if (s.prn) this.satIdx.set(s.prn, k);
@@ -14697,7 +14704,9 @@ var SbasProcessor = class {
       iiono: u(117, 9) || 1,
       cionoRamp: u(126, 10) * 5e-6,
       rssUdre: u(136, 1) === 1,
-      rssIono: u(137, 1) === 1
+      rssIono: u(137, 1) === 1,
+      ccovariance: u(138, 7) * 0.1
+      // Table B-49: 0..12.7 m, 0.1 m
     };
     return true;
   }
@@ -14749,6 +14758,51 @@ var SbasProcessor = class {
     }
     return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
   }
+  /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+   *  two satellites per message, each an upper-triangular Cholesky factor E and a
+   *  scale exponent that together define the relative clock-ephemeris covariance
+   *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+   *  the current MT1 mask, so a mismatched IODP is dropped. */
+  decodeMt28(msg) {
+    if (this.iodp !== getBitU(msg, 14, 2)) return false;
+    for (const base of [16, 121]) {
+      const slot = getBitU(msg, base, 6);
+      const prn = this.sats[slot]?.prn;
+      const scaleExp = getBitU(msg, base + 6, 3);
+      const e = [
+        getBitU(msg, base + 9, 9),
+        getBitU(msg, base + 18, 9),
+        getBitU(msg, base + 27, 9),
+        getBitU(msg, base + 36, 9),
+        getBitS(msg, base + 45, 10),
+        getBitS(msg, base + 55, 10),
+        getBitS(msg, base + 65, 10),
+        getBitS(msg, base + 75, 10),
+        getBitS(msg, base + 85, 10),
+        getBitS(msg, base + 95, 10)
+      ];
+      if (prn && e.some((v) => v !== 0)) this.mt28.set(prn, { scaleExp, e });
+    }
+    return true;
+  }
+  /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+   *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+   *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+   *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+  mt28DeltaUdre(prn, los) {
+    const m = this.mt28.get(prn);
+    if (!m) return null;
+    const sf = 2 ** (m.scaleExp - 5);
+    const [e11, e22, e33, e44, e12, e13, e14, e23, e24, e34] = m.e;
+    const [ix, iy, iz] = los;
+    const v0 = (e11 * ix + e12 * iy + e13 * iz + e14) * sf;
+    const v1 = (e22 * iy + e23 * iz + e24) * sf;
+    const v2 = (e33 * iz + e34) * sf;
+    const v3 = e44 * sf;
+    const iCi = v0 * v0 + v1 * v1 + v2 * v2 + v3 * v3;
+    const ec = (this.degr?.ccovariance ?? 0) * sf;
+    return Math.sqrt(iCi) + ec;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -14953,7 +15007,7 @@ var SbasProcessor = class {
    * no valid (unexpired) long-term + fast correction is available. Add `dPos`
    * to the satellite ECEF position and `dClkS` to its clock bias.
    */
-  satCorrection(prn, week, tow, userLatRad, userLonRad) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad, losUnitEcef) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -14975,7 +15029,8 @@ 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 dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+    const mt28 = losUnitEcef ? this.mt28DeltaUdre(prn, losUnitEcef) : null;
+    const dudre = mt28 != null ? mt28 : userLatRad != null && userLonRad != null ? this.deltaUdre(
       userLatRad * 180 / Math.PI,
       userLonRad * 180 / Math.PI
     ) : 1;
@@ -16695,7 +16750,10 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
+        const dx = sat.x - x, dy = sat.y - y, dz = sat.z - z;
+        const dr = Math.hypot(dx, dy, dz) || 1;
+        const los = [dx / dr, dy / dr, dz / dr];
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon, los);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/index.js

@@ -209,7 +209,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-TYJNN77C.js";
+} from "./chunk-B2XWO4RP.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,

package/dist/positioning.cjs

@@ -2132,6 +2132,10 @@ var SbasProcessor = class {
   tlat = 0;
   degr = null;
   mt27 = null;
+  /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+   *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+   *  scale exponent. Keyed by PRN; cleared on a mask change. */
+  mt28 = /* @__PURE__ */ new Map();
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -2162,6 +2166,8 @@ var SbasProcessor = class {
         return this.decodeMt10(msg) ? type : -1;
       case 27:
         return this.decodeMt27(msg) ? type : -1;
+      case 28:
+        return this.decodeMt28(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -2188,6 +2194,7 @@ var SbasProcessor = class {
       return true;
     this.iodp = iodp;
     this.sats = prns.map((prn) => ({ prn }));
+    this.mt28.clear();
     this.satIdx.clear();
     this.sats.forEach((s, k) => {
       if (s.prn) this.satIdx.set(s.prn, k);
@@ -2253,7 +2260,9 @@ var SbasProcessor = class {
       iiono: u(117, 9) || 1,
       cionoRamp: u(126, 10) * 5e-6,
       rssUdre: u(136, 1) === 1,
-      rssIono: u(137, 1) === 1
+      rssIono: u(137, 1) === 1,
+      ccovariance: u(138, 7) * 0.1
+      // Table B-49: 0..12.7 m, 0.1 m
     };
     return true;
   }
@@ -2305,6 +2314,51 @@ var SbasProcessor = class {
     }
     return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
   }
+  /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+   *  two satellites per message, each an upper-triangular Cholesky factor E and a
+   *  scale exponent that together define the relative clock-ephemeris covariance
+   *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+   *  the current MT1 mask, so a mismatched IODP is dropped. */
+  decodeMt28(msg) {
+    if (this.iodp !== getBitU(msg, 14, 2)) return false;
+    for (const base of [16, 121]) {
+      const slot = getBitU(msg, base, 6);
+      const prn = this.sats[slot]?.prn;
+      const scaleExp = getBitU(msg, base + 6, 3);
+      const e = [
+        getBitU(msg, base + 9, 9),
+        getBitU(msg, base + 18, 9),
+        getBitU(msg, base + 27, 9),
+        getBitU(msg, base + 36, 9),
+        getBitS(msg, base + 45, 10),
+        getBitS(msg, base + 55, 10),
+        getBitS(msg, base + 65, 10),
+        getBitS(msg, base + 75, 10),
+        getBitS(msg, base + 85, 10),
+        getBitS(msg, base + 95, 10)
+      ];
+      if (prn && e.some((v) => v !== 0)) this.mt28.set(prn, { scaleExp, e });
+    }
+    return true;
+  }
+  /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+   *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+   *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+   *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+  mt28DeltaUdre(prn, los) {
+    const m = this.mt28.get(prn);
+    if (!m) return null;
+    const sf = 2 ** (m.scaleExp - 5);
+    const [e11, e22, e33, e44, e12, e13, e14, e23, e24, e34] = m.e;
+    const [ix, iy, iz] = los;
+    const v0 = (e11 * ix + e12 * iy + e13 * iz + e14) * sf;
+    const v1 = (e22 * iy + e23 * iz + e24) * sf;
+    const v2 = (e33 * iz + e34) * sf;
+    const v3 = e44 * sf;
+    const iCi = v0 * v0 + v1 * v1 + v2 * v2 + v3 * v3;
+    const ec = (this.degr?.ccovariance ?? 0) * sf;
+    return Math.sqrt(iCi) + ec;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -2509,7 +2563,7 @@ var SbasProcessor = class {
    * no valid (unexpired) long-term + fast correction is available. Add `dPos`
    * to the satellite ECEF position and `dClkS` to its clock bias.
    */
-  satCorrection(prn, week, tow, userLatRad, userLonRad) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad, losUnitEcef) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -2531,7 +2585,8 @@ 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 dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+    const mt28 = losUnitEcef ? this.mt28DeltaUdre(prn, losUnitEcef) : null;
+    const dudre = mt28 != null ? mt28 : userLatRad != null && userLonRad != null ? this.deltaUdre(
       userLatRad * 180 / Math.PI,
       userLonRad * 180 / Math.PI
     ) : 1;
@@ -4323,7 +4378,10 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
+        const dx = sat.x - x, dy = sat.y - y, dz = sat.z - z;
+        const dr = Math.hypot(dx, dy, dz) || 1;
+        const los = [dx / dr, dy / dr, dz / dr];
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon, los);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/positioning.d.cts

@@ -668,6 +668,8 @@ interface Degradation {
     cionoRamp: number;
     rssUdre: boolean;
     rssIono: boolean;
+    /** Cᶜᵒᵛᵃʳⁱᵃⁿᶜᵉ — the MT28 δUDRE quantization term (m), 0 when no MT28. */
+    ccovariance: number;
 }
 /** A per-satellite SBAS correction, ready to apply in an SPP solve. */
 interface SbasSatCorrection {
@@ -715,6 +717,10 @@ declare class SbasProcessor {
     private tlat;
     private degr;
     private mt27;
+    /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+     *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+     *  scale exponent. Keyed by PRN; cleared on a mask change. */
+    private mt28;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -739,6 +745,17 @@ declare class SbasProcessor {
      *  set: the highest-priority containing region wins (ties → lower δUDRE),
      *  else the outside value. 1 when no MT27 has been received. */
     deltaUdre(latDeg: number, lonDeg: number): number;
+    /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+     *  two satellites per message, each an upper-triangular Cholesky factor E and a
+     *  scale exponent that together define the relative clock-ephemeris covariance
+     *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+     *  the current MT1 mask, so a mismatched IODP is dropped. */
+    private decodeMt28;
+    /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+     *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+     *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+     *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+    private mt28DeltaUdre;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -782,7 +799,11 @@ declare class SbasProcessor {
      * 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, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number, 
+    /** Unit line-of-sight vector user→satellite in ECEF, for the MT28
+     *  geometry-dependent δUDRE (final protection-level pass). When given and
+     *  MT28 data exists for the PRN, it supersedes the MT27 region δUDRE. */
+    losUnitEcef?: readonly [number, number, number]): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -1557,7 +1578,10 @@ interface SbasCorrectionSource {
     satCorrection(prn: string, week: number, tow: number, 
     /** User geodetic lat/lon (rad) for the MT27 service-region δUDRE, when
      *  known (the final protection-level pass). Omitted during iteration. */
-    userLatRad?: number, userLonRad?: number): {
+    userLatRad?: number, userLonRad?: number, 
+    /** Unit LOS user→satellite (ECEF) for the MT28 geometry-dependent δUDRE
+     *  (final protection-level pass). */
+    losUnitEcef?: readonly [number, number, number]): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;

package/dist/positioning.d.ts

@@ -668,6 +668,8 @@ interface Degradation {
     cionoRamp: number;
     rssUdre: boolean;
     rssIono: boolean;
+    /** Cᶜᵒᵛᵃʳⁱᵃⁿᶜᵉ — the MT28 δUDRE quantization term (m), 0 when no MT28. */
+    ccovariance: number;
 }
 /** A per-satellite SBAS correction, ready to apply in an SPP solve. */
 interface SbasSatCorrection {
@@ -715,6 +717,10 @@ declare class SbasProcessor {
     private tlat;
     private degr;
     private mt27;
+    /** MT28 clock-ephemeris covariance per PRN: the upper-triangular Cholesky
+     *  factor E (10 elements [E11,E22,E33,E44,E12,E13,E14,E23,E24,E34]) and its
+     *  scale exponent. Keyed by PRN; cleared on a mask change. */
+    private mt28;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -739,6 +745,17 @@ declare class SbasProcessor {
      *  set: the highest-priority containing region wins (ties → lower δUDRE),
      *  else the outside value. 1 when no MT27 has been received. */
     deltaUdre(latDeg: number, lonDeg: number): number;
+    /** MT28 clock-ephemeris covariance (ICAO Annex 10 Vol I §3.5.4.7 / Table B-51):
+     *  two satellites per message, each an upper-triangular Cholesky factor E and a
+     *  scale exponent that together define the relative clock-ephemeris covariance
+     *  used for a geometry-dependent δUDRE (§3.5.5.6.2.5). Slots are indexed into
+     *  the current MT1 mask, so a mismatched IODP is dropped. */
+    private decodeMt28;
+    /** Geometry-dependent δUDRE multiplier from MT28 (§3.5.5.6.2.5):
+     *  δUDRE = √(Iᵀ·C·I) + ε_c, with C = Rᵀ·R, R = E·2^(scaleExp−5),
+     *  I = [î_user→sat, 1], ε_c = C_covariance·2^(scaleExp−5). Since C = RᵀR,
+     *  Iᵀ·C·I = ‖R·I‖², so we form R·I directly. Null when no MT28 for `prn`. */
+    private mt28DeltaUdre;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -782,7 +799,11 @@ declare class SbasProcessor {
      * 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, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number, 
+    /** Unit line-of-sight vector user→satellite in ECEF, for the MT28
+     *  geometry-dependent δUDRE (final protection-level pass). When given and
+     *  MT28 data exists for the PRN, it supersedes the MT27 region δUDRE. */
+    losUnitEcef?: readonly [number, number, number]): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -1557,7 +1578,10 @@ interface SbasCorrectionSource {
     satCorrection(prn: string, week: number, tow: number, 
     /** User geodetic lat/lon (rad) for the MT27 service-region δUDRE, when
      *  known (the final protection-level pass). Omitted during iteration. */
-    userLatRad?: number, userLonRad?: number): {
+    userLatRad?: number, userLonRad?: number, 
+    /** Unit LOS user→satellite (ECEF) for the MT28 geometry-dependent δUDRE
+     *  (final protection-level pass). */
+    losUnitEcef?: readonly [number, number, number]): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;

package/dist/positioning.js

@@ -39,7 +39,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-TYJNN77C.js";
+} from "./chunk-B2XWO4RP.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-JYWYQ7DW.js";
 import "./chunk-HVXYFUCB.js";

package/package.json

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