gnss-js 1.63.0 → 1.64.0

package/dist/{chunk-QN7Q2ALS.js → chunk-5SQLUGWE.js}

@@ -563,6 +563,192 @@ function decodeEphemeris(frame) {
   }
 }
 
+// src/rtcm3/igs-ssr.ts
+var UPDATE_INTERVAL_S = [
+  1,
+  2,
+  5,
+  10,
+  15,
+  30,
+  60,
+  120,
+  240,
+  300,
+  600,
+  900,
+  1800,
+  3600,
+  7200,
+  10800
+];
+var SYSTEM_BY_BASE = {
+  20: "G",
+  40: "R",
+  60: "E",
+  80: "J",
+  100: "C",
+  120: "S"
+};
+var KIND_BY_OFFSET = {
+  1: "orbit",
+  2: "clock",
+  3: "combined",
+  4: "highRateClock",
+  5: "codeBias",
+  6: "phaseBias",
+  7: "ura"
+};
+var two = (n) => n.toString().padStart(2, "0");
+function prnFor(system, id) {
+  if (system === "J") return `J${two(id)}`;
+  if (system === "S") return `S${two(19 + id)}`;
+  return `${system}${two(id)}`;
+}
+function ssrUraMm(v) {
+  if (v === 0) return null;
+  const cls = v >> 3 & 7;
+  const val = v & 7;
+  return 3 ** cls * (1 + val / 4) - 1;
+}
+function decodeIgsSsr(frame) {
+  if (frame.messageType !== 4076) return null;
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const version = r.readU(3);
+  const im = r.readU(8);
+  if (im === 201) return decodeIono(r, version, im);
+  const base = Math.floor(im / 20) * 20;
+  const system = SYSTEM_BY_BASE[base];
+  const kind = KIND_BY_OFFSET[im - base];
+  if (!system || !kind) return null;
+  const epochS = r.readU(20);
+  const updateIntervalS = UPDATE_INTERVAL_S[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const hasDatum = kind === "orbit" || kind === "combined";
+  const referenceDatum = hasDatum ? r.readU(1) : void 0;
+  let dispersiveConsistency;
+  let mwConsistency;
+  if (kind === "phaseBias") {
+    dispersiveConsistency = r.readU(1) === 1;
+    mwConsistency = r.readU(1) === 1;
+  }
+  const nsat = r.readU(6);
+  const readOrbit = (s) => {
+    s.iode = r.readU(8);
+    s.deltaRadial = r.readS(22) * 1e-4;
+    s.deltaAlongTrack = r.readS(20) * 4e-4;
+    s.deltaCrossTrack = r.readS(20) * 4e-4;
+    s.dotRadial = r.readS(21) * 1e-6;
+    s.dotAlongTrack = r.readS(19) * 4e-6;
+    s.dotCrossTrack = r.readS(19) * 4e-6;
+  };
+  const readClock = (s) => {
+    s.c0 = r.readS(22) * 1e-4;
+    s.c1 = r.readS(21) * 1e-6;
+    s.c2 = r.readS(27) * 2e-8;
+  };
+  const satellites = [];
+  for (let i = 0; i < nsat; i++) {
+    const id = r.readU(6);
+    const sat = { prn: prnFor(system, id) };
+    if (kind === "orbit" || kind === "combined") readOrbit(sat);
+    if (kind === "clock" || kind === "combined") readClock(sat);
+    else if (kind === "ura")
+      sat.uraMm = ssrUraMm(r.readU(6));
+    else if (kind === "highRateClock")
+      sat.highRateClock = r.readS(22) * 1e-4;
+    else if (kind === "codeBias") {
+      const n = r.readU(5);
+      const biases = [];
+      for (let b = 0; b < n; b++) {
+        const signal = r.readU(5);
+        biases.push({
+          signal,
+          signalName: `#${signal}`,
+          biasM: r.readS(14) * 0.01
+          // IDF025 0.01 m
+        });
+      }
+      sat.codeBiases = biases;
+    } else if (kind === "phaseBias") {
+      const n = r.readU(5);
+      r.readU(9);
+      r.readS(8);
+      const biases = [];
+      for (let b = 0; b < n; b++) {
+        const signal = r.readU(5);
+        const integer = r.readU(1) === 1;
+        const wideLaneGroup = r.readU(2);
+        const discontinuity = r.readU(4);
+        const biasM = r.readS(20) * 1e-4;
+        biases.push({ signal, integer, wideLaneGroup, discontinuity, biasM });
+      }
+      sat.phaseBiases = biases;
+    }
+    satellites.push(sat);
+  }
+  return {
+    messageType: 4076,
+    version,
+    igsMessageNumber: im,
+    system,
+    kind,
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    referenceDatum,
+    iodSsr,
+    providerId,
+    solutionId,
+    dispersiveConsistency,
+    mwConsistency,
+    satellites
+  };
+}
+function decodeIono(r, version, im) {
+  const epochS = r.readU(20);
+  const updateIntervalS = UPDATE_INTERVAL_S[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const qualityTecu = r.readU(9) * 0.05;
+  const nLayers = r.readU(2) + 1;
+  const layers = [];
+  for (let l = 0; l < nLayers; l++) {
+    const heightKm = r.readU(8) * 10;
+    const degree = r.readU(4) + 1;
+    const order = r.readU(4) + 1;
+    const N = degree;
+    const M = order;
+    const nCoef = (N + 1) * (N + 1) - (N - M) * (N - M + 1);
+    const nCos = nCoef;
+    let nSin = 0;
+    for (let m = 1; m <= M; m++) nSin += N - m + 1;
+    for (let i = 0; i < nCos + nSin; i++) r.readS(16);
+    layers.push({ heightKm, degree, order, coefficientCount: nCos + nSin });
+  }
+  return {
+    messageType: 4076,
+    version,
+    igsMessageNumber: im,
+    system: "I",
+    kind: "ionosphere",
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    iodSsr,
+    providerId,
+    solutionId,
+    satellites: [],
+    iono: { qualityTecu, layers }
+  };
+}
+
 // src/rtcm3/legacy-obs.ts
 var PRUNIT_GPS = 299792.458;
 var PRUNIT_GLO = 599584.916;
@@ -579,7 +765,7 @@ function legacyLockSec(i) {
   if (i < 127) return 32 * i - 3096;
   return 937;
 }
-var two = (n) => n.toString().padStart(2, "0");
+var two2 = (n) => n.toString().padStart(2, "0");
 function decodeLegacyObs(frame) {
   const type = frame.messageType;
   const gps = type >= 1001 && type <= 1004;
@@ -665,7 +851,7 @@ function decodeLegacyObs(frame) {
       signals.push(sig);
     }
     observations.push({
-      prn: `${system}${two(satId)}`,
+      prn: `${system}${two2(satId)}`,
       system,
       signals
     });
@@ -1207,6 +1393,144 @@ function resetGloFreqCache() {
   gloFreqNum.fill(-128);
 }
 
+// src/rtcm3/ssr.ts
+var UPDATE_INTERVAL_S2 = [
+  1,
+  2,
+  5,
+  10,
+  15,
+  30,
+  60,
+  120,
+  240,
+  300,
+  600,
+  900,
+  1800,
+  3600,
+  7200,
+  10800
+];
+var GPS_SIGNAL = {
+  0: "L1 C/A",
+  1: "L1 P",
+  2: "L1 Z-tracking",
+  5: "L2 C/A",
+  6: "L2 semi-codeless",
+  7: "L2C(M)",
+  8: "L2C(L)",
+  9: "L2C(M+L)",
+  10: "L2 P",
+  11: "L2 Z-tracking",
+  14: "L5 I",
+  15: "L5 Q"
+};
+var GLO_SIGNAL = {
+  0: "G1 C/A",
+  1: "G1 P",
+  2: "G2 C/A",
+  3: "G2 P"
+};
+var SSR_KIND = {
+  1057: { system: "G", kind: "orbit" },
+  1058: { system: "G", kind: "clock" },
+  1059: { system: "G", kind: "codeBias" },
+  1060: { system: "G", kind: "combined" },
+  1061: { system: "G", kind: "ura" },
+  1062: { system: "G", kind: "highRateClock" },
+  1063: { system: "R", kind: "orbit" },
+  1064: { system: "R", kind: "clock" },
+  1065: { system: "R", kind: "codeBias" },
+  1066: { system: "R", kind: "combined" },
+  1067: { system: "R", kind: "ura" },
+  1068: { system: "R", kind: "highRateClock" }
+};
+function ssrUraMm2(v) {
+  if (v === 0) return null;
+  const cls = v >> 3 & 7;
+  const val = v & 7;
+  return 3 ** cls * (1 + val / 4) - 1;
+}
+var two3 = (n) => n.toString().padStart(2, "0");
+function decodeSsr(frame) {
+  const spec = SSR_KIND[frame.messageType];
+  if (!spec) return null;
+  const { system, kind } = spec;
+  const gps = system === "G";
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const epochS = r.readU(gps ? 20 : 17);
+  const updateIntervalS = UPDATE_INTERVAL_S2[r.readU(4)];
+  const multipleMessage = r.readU(1) === 1;
+  const hasDatum = kind === "orbit" || kind === "combined";
+  const referenceDatum = hasDatum ? r.readU(1) : void 0;
+  const iodSsr = r.readU(4);
+  const providerId = r.readU(16);
+  const solutionId = r.readU(4);
+  const nsat = r.readU(6);
+  const sigName = gps ? GPS_SIGNAL : GLO_SIGNAL;
+  const satellites = [];
+  const readOrbit = (sat) => {
+    sat.deltaRadial = r.readS(22) * 1e-4;
+    sat.deltaAlongTrack = r.readS(20) * 4e-4;
+    sat.deltaCrossTrack = r.readS(20) * 4e-4;
+    sat.dotRadial = r.readS(21) * 1e-6;
+    sat.dotAlongTrack = r.readS(19) * 4e-6;
+    sat.dotCrossTrack = r.readS(19) * 4e-6;
+  };
+  const readClock = (sat) => {
+    sat.c0 = r.readS(22) * 1e-4;
+    sat.c1 = r.readS(21) * 1e-6;
+    sat.c2 = r.readS(27) * 2e-8;
+  };
+  for (let i = 0; i < nsat; i++) {
+    const satId = gps ? r.readU(6) : r.readU(5);
+    const sat = { prn: `${system}${two3(satId)}` };
+    if (kind === "orbit" || kind === "combined") {
+      sat.iode = r.readU(8);
+      readOrbit(sat);
+    }
+    if (kind === "combined") readClock(sat);
+    else if (kind === "clock") readClock(sat);
+    else if (kind === "ura")
+      sat.uraMm = ssrUraMm2(r.readU(6));
+    else if (kind === "highRateClock")
+      sat.highRateClock = r.readS(22) * 1e-4;
+    else if (kind === "codeBias") {
+      const ncb = r.readU(5);
+      const biases = [];
+      for (let c = 0; c < ncb; c++) {
+        const signal = r.readU(5);
+        const biasM = r.readS(14) * 0.01;
+        biases.push({
+          signal,
+          signalName: sigName[signal] ?? `#${signal}`,
+          biasM
+        });
+      }
+      sat.codeBiases = biases;
+    }
+    satellites.push(sat);
+  }
+  return {
+    messageType: frame.messageType,
+    system,
+    kind,
+    epochS,
+    updateIntervalS,
+    multipleMessage,
+    referenceDatum,
+    iodSsr,
+    providerId,
+    solutionId,
+    satellites
+  };
+}
+function isSsrMessage(messageType) {
+  return messageType in SSR_KIND;
+}
+
 // src/rtcm3/station.ts
 function createStationMeta() {
   return {
@@ -1557,11 +1881,14 @@ export {
   Rtcm3Decoder,
   readString,
   decodeEphemeris,
+  decodeIgsSsr,
   decodeLegacyObs,
   decodeMsmFull,
   msmEpochToDate,
   setGloFreqNumber,
   resetGloFreqCache,
+  decodeSsr,
+  isSsrMessage,
   createStationMeta,
   updateStationMeta,
   RTCM3_MESSAGE_NAMES,