gnss-js 2.3.0 → 2.5.0

package/dist/index.d.cts

@@ -13,7 +13,7 @@ export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2Ob
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.cjs';
 export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.cjs';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.cjs';
-export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.cjs';
+export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.cjs';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.cjs';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.cjs';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';

package/dist/index.d.ts

@@ -13,7 +13,7 @@ export { RTCM2_MESSAGE_NAMES, Rtcm2Correction, Rtcm2Decoder, Rtcm2Frame, Rtcm2Ob
 export { B as BitReader, E as EphemerisInfo, R as Rtcm3DebugHandler, a as Rtcm3Decoder, b as Rtcm3Frame, d as decodeEphemeris, r as readString, c as reportDecodeError, s as setRtcm3DebugHandler } from './ephemeris-BAnBLLB9.js';
 export { IgsSsrIonoLayer, IgsSsrMessage, MessageTypeStats, ObsEpoch, ObsSatObs, ObsSignal, RTCM3_MESSAGE_NAMES, SatCn0, SignalCn0, SsrCodeBias, SsrKind, SsrMessage, SsrPhaseBias, SsrSatCorrection, StationMeta, StreamStats, SystemParams, SystemParamsMessage, createStationMeta, createStreamStats, decodeIgsSsr, decodeLegacyObs, decodeMsmFull, decodeObs, decodeSsr, decodeSystemParams, isSsrMessage, obsEpochToDate, resetGloFreqCache, rtcm3Constellation, setGloFreqNumber, updateStationMeta, updateStreamStats } from './rtcm3.js';
 export { UbxCnavResult, UbxFrame, UbxIonoUtcResult, UbxMeasurement, UbxNavOptions, UbxNavResult, UbxParseResult, UbxPvtResult, UbxRawNavOptions, UbxRawNavResult, UbxRawxEpoch, parseUbxCnav, parseUbxIonoUtc, parseUbxNav, parseUbxPvt, parseUbxRawNav, parseUbxRawx, ubxFrames } from './ubx.js';
-export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, feedGeoBlock, feedGeoL5Block, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.js';
+export { DfmcCensus, DfmcMessageCb, GeoL5Result, SBAS_L5_MT_NAMES, SBAS_L5_PREAMBLE_UW, SbasMessageCb, SbfBdsNavResult, SbfCnavEphemeris, SbfCnavResult, SbfGalEphemeris, SbfGalNavResult, SbfGeoL5Result, SbfGeoNavResult, SbfGloNavResult, SbfGpsNavResult, SbfHasMessage, SbfHasResult, SbfIonoUtcResult, SbfMeasEpoch, SbfMeasurement, SbfNavCounts, SbfNavicResult, SbfNavigation, SbfParseResult, SbfPvt, SbfPvtResult, SbfReceiverSetup, decodeSbfNavigation, extractSbfNavicSubframe, feedGeoBlock, feedGeoL5Block, feedNavicBlock, isSbasL1Preamble, isSbasL5Correction, parseSbfBdsNav, parseSbfCnav, parseSbfGalNav, parseSbfGeoL5, parseSbfGeoNav, parseSbfGloNav, parseSbfGpsNav, parseSbfHas, parseSbfIonoUtc, parseSbfMeas, parseSbfNavic, parseSbfPvt, parseSbfReceiverSetup, sbasL5MessageName, sbasL5MessageType, sbfCrc16, scanSbfFrames } from './sbf.js';
 export { NovatelEpoch, NovatelMeasurement, NovatelNavResult, NovatelParseResult, NovatelPvtResult, OEM4_HLEN, Oem4Frame, crc32, oem4Frames, parseNovatelNav, parseNovatelPvt, parseNovatelRange } from './novatel.js';
 export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, TrimbleParseResult, parseTrimble, parseTrimbleNav, trimbleFrames } from './trimble.js';
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';

package/dist/index.js

@@ -77,7 +77,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-XPQ4BU6W.js";
+} from "./chunk-7GNRNKBR.js";
 import {
   parseUbxCnav,
   parseUbxIonoUtc,
@@ -92,8 +92,10 @@ import {
   SBAS_L5_PREAMBLE_UW,
   crc16,
   decodeSbfNavigation,
+  extractSbfNavicSubframe,
   feedGeoBlock,
   feedGeoL5Block,
+  feedNavicBlock,
   isSbasL1Preamble,
   isSbasL5Correction,
   parseSbfAlmanac,
@@ -108,12 +110,13 @@ import {
   parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
+  parseSbfNavic,
   parseSbfPvt,
   parseSbfReceiverSetup,
   sbasL5MessageName,
   sbasL5MessageType,
   scanSbfFrames
-} from "./chunk-WW3DKCDN.js";
+} from "./chunk-74IVYOQV.js";
 import "./chunk-3KSDL3CO.js";
 import {
   OEM4_HLEN,
@@ -209,7 +212,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-TYJNN77C.js";
+} from "./chunk-B2XWO4RP.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,
@@ -542,9 +545,11 @@ export {
   estimateNetworkFcbs,
   estimateWidelaneFcb,
   euclidean3D,
+  extractSbfNavicSubframe,
   extractWidelaneArcs,
   feedGeoBlock,
   feedGeoL5Block,
+  feedNavicBlock,
   fetchSourcetable,
   findSatBias,
   fixPppPosition,
@@ -664,6 +669,7 @@ export {
   parseSbfIonoUtc,
   parseSbfMeas,
   parseSbfNav,
+  parseSbfNavic,
   parseSbfPvt,
   parseSbfReceiverSetup,
   parseSinexBiasDcb,

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/dist/rtcm3.cjs

@@ -564,6 +564,70 @@ function decodeQzssEphemeris(payload) {
     fitInterval
   };
 }
+function decodeNavicEphemeris(payload) {
+  if (payload.length < 60) return null;
+  const r = new BitReader(payload);
+  r.skip(12);
+  const svId = r.readU(6);
+  const week = r.readU(10);
+  const af0 = r.readS(22) * 2 ** -31;
+  const af1 = r.readS(16) * 2 ** -43;
+  const af2 = r.readS(8) * 2 ** -55;
+  const ura = r.readU(4);
+  const toc = r.readU(16) * 16;
+  const tgd = r.readS(8) * 2 ** -31;
+  const deltaN = r.readS(22) * 2 ** -41 * PI;
+  const iode = r.readU(8);
+  r.skip(10);
+  const l5Health = r.readU(1);
+  const sHealth = r.readU(1);
+  const cuc = r.readS(15) * 2 ** -28;
+  const cus = r.readS(15) * 2 ** -28;
+  const cic = r.readS(15) * 2 ** -28;
+  const cis = r.readS(15) * 2 ** -28;
+  const crc = r.readS(15) * 2 ** -4;
+  const crs = r.readS(15) * 2 ** -4;
+  const idot = r.readS(14) * 2 ** -43 * PI;
+  const m0 = r.readS(32) * 2 ** -31 * PI;
+  const toe = r.readU(16) * 16;
+  const e = r.readU(32) * 2 ** -33;
+  const sqrtA = r.readU(32) * 2 ** -19;
+  const omega0 = r.readS(32) * 2 ** -31 * PI;
+  const argPerigee = r.readS(32) * 2 ** -31 * PI;
+  const omegaDot = r.readS(22) * 2 ** -41 * PI;
+  const i0 = r.readS(32) * 2 ** -31 * PI;
+  return {
+    prn: `I${String(svId).padStart(2, "0")}`,
+    constellation: "NavIC",
+    health: l5Health << 1 | sHealth,
+    lastReceived: Date.now(),
+    messageType: 1041,
+    week,
+    ura,
+    toc,
+    toe,
+    sqrtA,
+    eccentricity: e,
+    inclination: i0,
+    omega0,
+    omegaDot,
+    argPerigee,
+    meanAnomaly: m0,
+    deltaN,
+    idot,
+    crs,
+    crc,
+    cuc,
+    cus,
+    cic,
+    cis,
+    af0,
+    af1,
+    af2,
+    iode,
+    tgd
+  };
+}
 function decodeEphemeris(frame) {
   try {
     switch (frame.messageType) {
@@ -571,6 +635,8 @@ function decodeEphemeris(frame) {
         return decodeGpsEphemeris(frame.payload);
       case 1020:
         return decodeGlonassEphemeris(frame.payload);
+      case 1041:
+        return decodeNavicEphemeris(frame.payload);
       case 1042:
         return decodeBdsEphemeris(frame.payload);
       case 1043:

package/dist/rtcm3.js

@@ -21,7 +21,7 @@ import {
   setRtcm3DebugHandler,
   updateStationMeta,
   updateStreamStats
-} from "./chunk-XPQ4BU6W.js";
+} from "./chunk-7GNRNKBR.js";
 import "./chunk-HVXYFUCB.js";
 import "./chunk-LEEU5OIO.js";
 import "./chunk-FIEWO4J4.js";