gnss-js 1.57.0 → 1.59.0

package/dist/{chunk-BU6YC6VK.js → chunk-7NZ4SJTJ.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,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];
@@ -1323,6 +1334,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();
@@ -1349,6 +1361,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:
@@ -1417,6 +1431,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);
@@ -1463,7 +1501,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;
   }
@@ -1487,7 +1526,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;
   }
@@ -1608,9 +1648,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
@@ -1659,14 +1706,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) {
@@ -3081,10 +3132,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;
@@ -3380,6 +3431,8 @@ export {
   solveDgnss,
   RtkFloatEngine,
   postProcessRtk,
+  sbasLongTermDeg,
+  sbasIonoDeg,
   SbasProcessor,
   niellMapping,
   gmst,

package/dist/index.cjs

@@ -303,8 +303,10 @@ __export(src_exports, {
   rhumbLine: () => rhumbLine,
   rtcm3Constellation: () => rtcm3Constellation,
   satClockCorrection: () => satClockCorrection,
+  sbasIonoDeg: () => sbasIonoDeg,
   sbasL5MessageName: () => sbasL5MessageName,
   sbasL5MessageType: () => sbasL5MessageType,
+  sbasLongTermDeg: () => sbasLongTermDeg,
   sbasProtectionLevels: () => sbasProtectionLevels,
   sbfCrc16: () => crc16,
   scanSbfFrames: () => scanSbfFrames,
@@ -13421,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 = [
@@ -13478,6 +13480,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];
@@ -13531,6 +13544,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();
@@ -13557,6 +13571,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:
@@ -13625,6 +13641,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);
@@ -13671,7 +13711,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;
   }
@@ -13695,7 +13736,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;
   }
@@ -13816,9 +13858,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
@@ -13867,14 +13916,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) {
@@ -15289,10 +15342,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;
@@ -18122,8 +18175,10 @@ function computePsdDb(centerMHz, halfSpanChips, numPoints, psdFn, f0 = 1023e3, f
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
+  sbasIonoDeg,
   sbasL5MessageName,
   sbasL5MessageType,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   sbfCrc16,
   scanSbfFrames,

package/dist/index.d.cts

@@ -18,7 +18,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.cjs';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.cjs';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.cjs';
-export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, 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 } from './positioning.cjs';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, 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 } from './positioning.cjs';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.cjs';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.cjs';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.cjs';

package/dist/index.d.ts

@@ -18,7 +18,7 @@ export { TrimbleEpoch, TrimbleFrame, TrimbleMeasurement, TrimbleNavResult, Trimb
 export { BinexEpoch, BinexMeasurement, BinexParseResult, BinexRecord, binexCrc16, binexCrc32, binexCsum8, binexRecords, decodeBinexEph, getBnxi, parseBinex } from './binex.js';
 export { AllPositionsData, AlmanacPosition, DopValues, EpochSkyData, SatAzEl, SatPoint, SatPosition, VisibilityPass, VisibilityResult, almanacEpochMs, almanacSatPosition, computeAllPositions, computeDop, computeLiveSkyPositions, computeSatPosition, computeVisibility, computeVisibilityFromPositions, dopplerHz, ecefToAzEl, ephInfoToEphemeris, glonassAlmanacEpochMs, glonassAlmanacPosition, glonassPosition, keplerPosition, maskRadForAzimuth, navTimesFromEph, rangeRate, selectEphemeris } from './orbit.js';
 export { Helmert14, REFERENCE_FRAMES, ReferenceFrame, applyHelmert, dateToEpoch, transformFrame } from './frames.js';
-export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, 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 } from './positioning.js';
+export { AntennaOffset, ArFixedSat, ArOptions, ArResult, ArSat, CorrectionResult, Degradation, DgnssOptions, DgnssSolution, EphemerisSource, FixedPositionResult, IONO_L1_M_PER_TECU, KlobucharCoeffs, LambdaResult, NetworkFcbResult, NetworkFcbSummary, NlArc, NlFcbResult, PppAntennaModel, PppArc, PppCorrectionConfig, PppCorrections, PppEpoch, PppEpochContext, PppEpochResult, PppFixState, PppOptions, PppSatObs, PppSolution, RawObservation, RtkEpochMeasurements, RtkFloatEngine, RtkFloatOptions, RtkFloatSolution, RtkMeasurement, RtkObsEpoch, RtkObsInput, RtkPostProcessOptions, RtkPostProcessResult, RtkTrackPoint, SatGeom, SbasCorrectionSource, SbasIonoDelay, SbasProcessor, SbasSatCorrection, SppOptions, SppSolution, TropoMapping, WlArc, WlFcbOptions, WlFcbResult, 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 } from './positioning.js';
 export { NtripCaster, NtripConnectionInfo, NtripNetwork, NtripStream, NtripStreamConnection, NtripVersion, Sourcetable, SourcetableEntry, connectToMountpoint, fetchSourcetable, parseSourcetable } from './ntrip.js';
 export { CompleteSatSignal, CompleteSignalStats, CompletenessAccumulator, CompletenessResult, CycleSlipAccumulator, CycleSlipEvent, CycleSlipResult, CycleSlipSignalStats, IonoAccumulator, IonoDcbResult, IonoPoint, IonoRatePoint, IonoRateSeries, IonoResult, IonoSeries, MultipathAccumulator, MultipathPoint, MultipathResult, MultipathSeries, MultipathSignalStat, QualityResult, SatDcbMap, analyzeQuality, applyIonoDcb, computeIonoRate, detrendIonoArcs, parseSinexBiasDcb } from './analysis.js';
 export { AntennaEntry, AntexFile, FrequencyData, frequencyLabel, parseAntex } from './antex.js';

package/dist/index.js

@@ -168,6 +168,8 @@ import {
   postProcessRtk,
   resolvePppAmbiguities,
   satClockCorrection,
+  sbasIonoDeg,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   solidEarthTide,
   solveDgnss,
@@ -176,7 +178,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-BU6YC6VK.js";
+} from "./chunk-7NZ4SJTJ.js";
 import "./chunk-WQX354KD.js";
 import {
   almanacEpochMs,
@@ -651,8 +653,10 @@ export {
   rhumbLine,
   rtcm3Constellation,
   satClockCorrection,
+  sbasIonoDeg,
   sbasL5MessageName,
   sbasL5MessageType,
+  sbasLongTermDeg,
   sbasProtectionLevels,
   crc16 as sbfCrc16,
   scanSbfFrames,

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,
@@ -1687,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 = [
@@ -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-7NZ4SJTJ.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/package.json

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