gnss-js 1.58.0 → 1.60.0

package/dist/{chunk-SWVIBUHS.js → chunk-4WNIGDNR.js}

@@ -1213,7 +1213,7 @@ var D2R = Math.PI / 180;
 var R2D = 180 / Math.PI;
 var SEC_PER_WEEK = 604800;
 var MAXSBSAGEF = 30;
-var MAXSBSAGEL = 1800;
+var MAXSBSAGEL = 360;
 var MAXBAND = 10;
 var two = (n) => String(n).padStart(2, "0");
 var VAR_FCORR = [
@@ -1270,6 +1270,44 @@ var DEGF = [
   58e-4
 ];
 var degfcorr = (ai) => ai > 0 && ai <= 15 ? DEGF[ai] : 58e-4;
+var DELTA_UDRE = [
+  1,
+  1.1,
+  1.25,
+  1.5,
+  2,
+  3,
+  4,
+  5,
+  6,
+  8,
+  10,
+  20,
+  30,
+  40,
+  50,
+  100
+];
+function pointInRegion(lat, lon, rg) {
+  const poly = rg.square ? [
+    [rg.lat1, rg.lon1],
+    [rg.lat1, rg.lon2],
+    [rg.lat2, rg.lon2],
+    [rg.lat2, rg.lon1]
+  ] : [
+    [rg.lat1, rg.lon1],
+    [rg.lat2, rg.lon2],
+    [rg.lat1, rg.lon2]
+  ];
+  let inside = false;
+  for (let i = 0, j = poly.length - 1; i < poly.length; j = i++) {
+    const [yi, xi] = poly[i];
+    const [yj, xj] = poly[j];
+    if (yi > lat !== yj > lat && lon < (xj - xi) * (lat - yi) / (yj - yi) + xi)
+      inside = !inside;
+  }
+  return inside;
+}
 function sbasLongTermDeg(d, vel, ageSec) {
   if (vel) {
     if (ageSec > 0 && ageSec < d.iltcV1) return 0;
@@ -1335,6 +1373,7 @@ var SbasProcessor = class {
   iodp = -1;
   tlat = 0;
   degr = null;
+  mt27 = null;
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -1363,6 +1402,8 @@ var SbasProcessor = class {
         return this.decodeDegradation(msg) ? type : -1;
       case 10:
         return this.decodeMt10(msg) ? type : -1;
+      case 27:
+        return this.decodeMt27(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -1455,6 +1496,54 @@ var SbasProcessor = class {
     };
     return true;
   }
+  /** MT27 SBAS service message (DO-229D §A.4.4.13, Table A-20): regions with a
+   *  δUDRE that inflates σ_UDRE for users inside/outside them. */
+  decodeMt27(msg) {
+    const iods = getBitU(msg, 14, 3);
+    const svcNum = getBitU(msg, 20, 3) + 1;
+    const nRegions = getBitU(msg, 23, 3);
+    const priority = getBitU(msg, 26, 2);
+    const udreInside = getBitU(msg, 28, 4);
+    const udreOutside = getBitU(msg, 32, 4);
+    if (!this.mt27 || this.mt27.iods !== iods)
+      this.mt27 = { iods, udreOutside, msgs: /* @__PURE__ */ new Map() };
+    this.mt27.udreOutside = udreOutside;
+    const regions = [];
+    for (let r = 0; r < nRegions && r < 5; r++) {
+      const p = 36 + r * 35;
+      regions.push({
+        lat1: getBitS(msg, p, 8),
+        lon1: getBitS(msg, p + 8, 9),
+        lat2: getBitS(msg, p + 17, 8),
+        lon2: getBitS(msg, p + 25, 9),
+        square: getBitU(msg, p + 34, 1) === 1
+      });
+    }
+    this.mt27.msgs.set(svcNum, { priority, udreInside, regions });
+    return true;
+  }
+  /** δUDRE multiplier (DO-229D Table A-21) for a user location from the MT27
+   *  set: the highest-priority containing region wins (ties → lower δUDRE),
+   *  else the outside value. 1 when no MT27 has been received. */
+  deltaUdre(latDeg, lonDeg) {
+    const s = this.mt27;
+    if (!s) return 1;
+    let best = null;
+    let bestPri = -1;
+    for (const m of s.msgs.values()) {
+      for (const rg of m.regions) {
+        if (!pointInRegion(latDeg, lonDeg, rg)) continue;
+        const val = DELTA_UDRE[m.udreInside] ?? 1;
+        if (m.priority > bestPri) {
+          bestPri = m.priority;
+          best = val;
+        } else if (m.priority === bestPri && best !== null) {
+          best = Math.min(best, val);
+        }
+      }
+    }
+    return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -1626,7 +1715,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) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -1648,7 +1737,11 @@ 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 sigUdre = Math.sqrt(varfcorr(s.fcorr.udre));
+    const dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+      userLatRad * 180 / Math.PI,
+      userLonRad * 180 / Math.PI
+    ) : 1;
+    const sigUdre = Math.sqrt(varfcorr(s.fcorr.udre)) * dudre;
     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;
@@ -3387,7 +3480,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/index.cjs

@@ -13423,7 +13423,7 @@ var D2R = Math.PI / 180;
 var R2D = 180 / Math.PI;
 var SEC_PER_WEEK12 = 604800;
 var MAXSBSAGEF = 30;
-var MAXSBSAGEL = 1800;
+var MAXSBSAGEL = 360;
 var MAXBAND = 10;
 var two8 = (n) => String(n).padStart(2, "0");
 var VAR_FCORR = [
@@ -13480,6 +13480,44 @@ var DEGF = [
   58e-4
 ];
 var degfcorr = (ai) => ai > 0 && ai <= 15 ? DEGF[ai] : 58e-4;
+var DELTA_UDRE = [
+  1,
+  1.1,
+  1.25,
+  1.5,
+  2,
+  3,
+  4,
+  5,
+  6,
+  8,
+  10,
+  20,
+  30,
+  40,
+  50,
+  100
+];
+function pointInRegion(lat, lon, rg) {
+  const poly = rg.square ? [
+    [rg.lat1, rg.lon1],
+    [rg.lat1, rg.lon2],
+    [rg.lat2, rg.lon2],
+    [rg.lat2, rg.lon1]
+  ] : [
+    [rg.lat1, rg.lon1],
+    [rg.lat2, rg.lon2],
+    [rg.lat1, rg.lon2]
+  ];
+  let inside = false;
+  for (let i = 0, j = poly.length - 1; i < poly.length; j = i++) {
+    const [yi, xi] = poly[i];
+    const [yj, xj] = poly[j];
+    if (yi > lat !== yj > lat && lon < (xj - xi) * (lat - yi) / (yj - yi) + xi)
+      inside = !inside;
+  }
+  return inside;
+}
 function sbasLongTermDeg(d, vel, ageSec) {
   if (vel) {
     if (ageSec > 0 && ageSec < d.iltcV1) return 0;
@@ -13545,6 +13583,7 @@ var SbasProcessor = class {
   iodp = -1;
   tlat = 0;
   degr = null;
+  mt27 = null;
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -13573,6 +13612,8 @@ var SbasProcessor = class {
         return this.decodeDegradation(msg) ? type : -1;
       case 10:
         return this.decodeMt10(msg) ? type : -1;
+      case 27:
+        return this.decodeMt27(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -13665,6 +13706,54 @@ var SbasProcessor = class {
     };
     return true;
   }
+  /** MT27 SBAS service message (DO-229D §A.4.4.13, Table A-20): regions with a
+   *  δUDRE that inflates σ_UDRE for users inside/outside them. */
+  decodeMt27(msg) {
+    const iods = getBitU(msg, 14, 3);
+    const svcNum = getBitU(msg, 20, 3) + 1;
+    const nRegions = getBitU(msg, 23, 3);
+    const priority = getBitU(msg, 26, 2);
+    const udreInside = getBitU(msg, 28, 4);
+    const udreOutside = getBitU(msg, 32, 4);
+    if (!this.mt27 || this.mt27.iods !== iods)
+      this.mt27 = { iods, udreOutside, msgs: /* @__PURE__ */ new Map() };
+    this.mt27.udreOutside = udreOutside;
+    const regions = [];
+    for (let r = 0; r < nRegions && r < 5; r++) {
+      const p = 36 + r * 35;
+      regions.push({
+        lat1: getBitS(msg, p, 8),
+        lon1: getBitS(msg, p + 8, 9),
+        lat2: getBitS(msg, p + 17, 8),
+        lon2: getBitS(msg, p + 25, 9),
+        square: getBitU(msg, p + 34, 1) === 1
+      });
+    }
+    this.mt27.msgs.set(svcNum, { priority, udreInside, regions });
+    return true;
+  }
+  /** δUDRE multiplier (DO-229D Table A-21) for a user location from the MT27
+   *  set: the highest-priority containing region wins (ties → lower δUDRE),
+   *  else the outside value. 1 when no MT27 has been received. */
+  deltaUdre(latDeg, lonDeg) {
+    const s = this.mt27;
+    if (!s) return 1;
+    let best = null;
+    let bestPri = -1;
+    for (const m of s.msgs.values()) {
+      for (const rg of m.regions) {
+        if (!pointInRegion(latDeg, lonDeg, rg)) continue;
+        const val = DELTA_UDRE[m.udreInside] ?? 1;
+        if (m.priority > bestPri) {
+          bestPri = m.priority;
+          best = val;
+        } else if (m.priority === bestPri && best !== null) {
+          best = Math.min(best, val);
+        }
+      }
+    }
+    return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -13836,7 +13925,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) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -13858,7 +13947,11 @@ 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 sigUdre = Math.sqrt(varfcorr(s.fcorr.udre));
+    const dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+      userLatRad * 180 / Math.PI,
+      userLonRad * 180 / Math.PI
+    ) : 1;
+    const sigUdre = Math.sqrt(varfcorr(s.fcorr.udre)) * dudre;
     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;
@@ -15597,7 +15690,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/index.js

@@ -178,7 +178,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-SWVIBUHS.js";
+} from "./chunk-4WNIGDNR.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,

package/dist/positioning.cjs

@@ -1689,7 +1689,7 @@ var D2R = Math.PI / 180;
 var R2D = 180 / Math.PI;
 var SEC_PER_WEEK2 = 604800;
 var MAXSBSAGEF = 30;
-var MAXSBSAGEL = 1800;
+var MAXSBSAGEL = 360;
 var MAXBAND = 10;
 var two = (n) => String(n).padStart(2, "0");
 var VAR_FCORR = [
@@ -1746,6 +1746,44 @@ var DEGF = [
   58e-4
 ];
 var degfcorr = (ai) => ai > 0 && ai <= 15 ? DEGF[ai] : 58e-4;
+var DELTA_UDRE = [
+  1,
+  1.1,
+  1.25,
+  1.5,
+  2,
+  3,
+  4,
+  5,
+  6,
+  8,
+  10,
+  20,
+  30,
+  40,
+  50,
+  100
+];
+function pointInRegion(lat, lon, rg) {
+  const poly = rg.square ? [
+    [rg.lat1, rg.lon1],
+    [rg.lat1, rg.lon2],
+    [rg.lat2, rg.lon2],
+    [rg.lat2, rg.lon1]
+  ] : [
+    [rg.lat1, rg.lon1],
+    [rg.lat2, rg.lon2],
+    [rg.lat1, rg.lon2]
+  ];
+  let inside = false;
+  for (let i = 0, j = poly.length - 1; i < poly.length; j = i++) {
+    const [yi, xi] = poly[i];
+    const [yj, xj] = poly[j];
+    if (yi > lat !== yj > lat && lon < (xj - xi) * (lat - yi) / (yj - yi) + xi)
+      inside = !inside;
+  }
+  return inside;
+}
 function sbasLongTermDeg(d, vel, ageSec) {
   if (vel) {
     if (ageSec > 0 && ageSec < d.iltcV1) return 0;
@@ -1811,6 +1849,7 @@ var SbasProcessor = class {
   iodp = -1;
   tlat = 0;
   degr = null;
+  mt27 = null;
   ion = [];
   /** GEO PRN(s) whose messages have been ingested. */
   geoPrns = /* @__PURE__ */ new Set();
@@ -1839,6 +1878,8 @@ var SbasProcessor = class {
         return this.decodeDegradation(msg) ? type : -1;
       case 10:
         return this.decodeMt10(msg) ? type : -1;
+      case 27:
+        return this.decodeMt27(msg) ? type : -1;
       case 18:
         return this.decodeIgpMask(msg) ? type : -1;
       case 24:
@@ -1931,6 +1972,54 @@ var SbasProcessor = class {
     };
     return true;
   }
+  /** MT27 SBAS service message (DO-229D §A.4.4.13, Table A-20): regions with a
+   *  δUDRE that inflates σ_UDRE for users inside/outside them. */
+  decodeMt27(msg) {
+    const iods = getBitU(msg, 14, 3);
+    const svcNum = getBitU(msg, 20, 3) + 1;
+    const nRegions = getBitU(msg, 23, 3);
+    const priority = getBitU(msg, 26, 2);
+    const udreInside = getBitU(msg, 28, 4);
+    const udreOutside = getBitU(msg, 32, 4);
+    if (!this.mt27 || this.mt27.iods !== iods)
+      this.mt27 = { iods, udreOutside, msgs: /* @__PURE__ */ new Map() };
+    this.mt27.udreOutside = udreOutside;
+    const regions = [];
+    for (let r = 0; r < nRegions && r < 5; r++) {
+      const p = 36 + r * 35;
+      regions.push({
+        lat1: getBitS(msg, p, 8),
+        lon1: getBitS(msg, p + 8, 9),
+        lat2: getBitS(msg, p + 17, 8),
+        lon2: getBitS(msg, p + 25, 9),
+        square: getBitU(msg, p + 34, 1) === 1
+      });
+    }
+    this.mt27.msgs.set(svcNum, { priority, udreInside, regions });
+    return true;
+  }
+  /** δUDRE multiplier (DO-229D Table A-21) for a user location from the MT27
+   *  set: the highest-priority containing region wins (ties → lower δUDRE),
+   *  else the outside value. 1 when no MT27 has been received. */
+  deltaUdre(latDeg, lonDeg) {
+    const s = this.mt27;
+    if (!s) return 1;
+    let best = null;
+    let bestPri = -1;
+    for (const m of s.msgs.values()) {
+      for (const rg of m.regions) {
+        if (!pointInRegion(latDeg, lonDeg, rg)) continue;
+        const val = DELTA_UDRE[m.udreInside] ?? 1;
+        if (m.priority > bestPri) {
+          bestPri = m.priority;
+          best = val;
+        } else if (m.priority === bestPri && best !== null) {
+          best = Math.min(best, val);
+        }
+      }
+    }
+    return best ?? DELTA_UDRE[s.udreOutside] ?? 1;
+  }
   decodeMixed(msg, week, tow) {
     if (this.iodp !== getBitU(msg, 110, 2)) return false;
     const blk = getBitU(msg, 112, 2);
@@ -2102,7 +2191,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) {
+  satCorrection(prn, week, tow, userLatRad, userLonRad) {
     const idx = this.satIdx.get(prn);
     if (idx === void 0) return null;
     const s = this.sats[idx];
@@ -2124,7 +2213,11 @@ 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 sigUdre = Math.sqrt(varfcorr(s.fcorr.udre));
+    const dudre = userLatRad != null && userLonRad != null ? this.deltaUdre(
+      userLatRad * 180 / Math.PI,
+      userLonRad * 180 / Math.PI
+    ) : 1;
+    const sigUdre = Math.sqrt(varfcorr(s.fcorr.udre)) * dudre;
     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;
@@ -3935,7 +4028,7 @@ function solveSpp(pseudoranges, ephemerides, timeMs, opts = {}) {
       used.push(prn);
       azels.push(azel);
       if (sbas) {
-        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow);
+        const sc = sbas.satCorrection(prn, gpsWeek, gpsTow, rxLat, rxLon);
         const si = sbas.ionoDelay(
           gpsWeek,
           gpsTow,

package/dist/positioning.d.cts

@@ -605,6 +605,7 @@ declare class SbasProcessor {
     private iodp;
     private tlat;
     private degr;
+    private mt27;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -622,6 +623,13 @@ declare class SbasProcessor {
      *  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;
+    /** MT27 SBAS service message (DO-229D §A.4.4.13, Table A-20): regions with a
+     *  δUDRE that inflates σ_UDRE for users inside/outside them. */
+    private decodeMt27;
+    /** δUDRE multiplier (DO-229D Table A-21) for a user location from the MT27
+     *  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;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -657,7 +665,7 @@ 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): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -1292,7 +1300,10 @@ interface SppOptions {
 /** The two SBAS corrections `solveSpp` consumes — {@link ./sbas.SbasProcessor}
  *  satisfies this structurally. */
 interface SbasCorrectionSource {
-    satCorrection(prn: string, week: number, tow: number): {
+    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): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;

package/dist/positioning.d.ts

@@ -605,6 +605,7 @@ declare class SbasProcessor {
     private iodp;
     private tlat;
     private degr;
+    private mt27;
     private ion;
     /** GEO PRN(s) whose messages have been ingested. */
     readonly geoPrns: Set<string>;
@@ -622,6 +623,13 @@ declare class SbasProcessor {
      *  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;
+    /** MT27 SBAS service message (DO-229D §A.4.4.13, Table A-20): regions with a
+     *  δUDRE that inflates σ_UDRE for users inside/outside them. */
+    private decodeMt27;
+    /** δUDRE multiplier (DO-229D Table A-21) for a user location from the MT27
+     *  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;
     private decodeMixed;
     private decodeLongCorrHalf;
     private decodeLong0;
@@ -657,7 +665,7 @@ 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): SbasSatCorrection | null;
+    satCorrection(prn: string, week: number, tow: number, userLatRad?: number, userLonRad?: number): SbasSatCorrection | null;
     /** The last-decoded MT10 degradation factors, or null if none seen. */
     get degradation(): Readonly<Degradation> | null;
     /**
@@ -1292,7 +1300,10 @@ interface SppOptions {
 /** The two SBAS corrections `solveSpp` consumes — {@link ./sbas.SbasProcessor}
  *  satisfies this structurally. */
 interface SbasCorrectionSource {
-    satCorrection(prn: string, week: number, tow: number): {
+    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): {
         dPos: [number, number, number];
         dClkS: number;
         varM2: number;

package/dist/positioning.js

@@ -32,7 +32,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-SWVIBUHS.js";
+} from "./chunk-4WNIGDNR.js";
 import "./chunk-WQX354KD.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-MIIM4LDY.js";

package/package.json

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