gnss-js 1.62.0 → 1.64.0

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

@@ -563,6 +563,302 @@ 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;
+var LAM_GPS_L1 = C_LIGHT / 157542e4;
+var LAM_GPS_L2 = C_LIGHT / 12276e5;
+var PPR_INVALID = -(2 ** 19);
+var PR21_INVALID = -(2 ** 13);
+function legacyLockSec(i) {
+  if (i < 24) return i;
+  if (i < 48) return 2 * i - 24;
+  if (i < 72) return 4 * i - 120;
+  if (i < 96) return 8 * i - 408;
+  if (i < 120) return 16 * i - 1176;
+  if (i < 127) return 32 * i - 3096;
+  return 937;
+}
+var two2 = (n) => n.toString().padStart(2, "0");
+function decodeLegacyObs(frame) {
+  const type = frame.messageType;
+  const gps = type >= 1001 && type <= 1004;
+  const glo = type >= 1009 && type <= 1012;
+  if (!gps && !glo) return null;
+  const ext = type === 1002 || type === 1004 || type === 1010 || type === 1012;
+  const dual = type === 1003 || type === 1004 || type === 1011 || type === 1012;
+  const system = gps ? "G" : "R";
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  r.readU(12);
+  let epochMs;
+  let nsat;
+  if (gps) {
+    epochMs = r.readU(30);
+    r.readU(1);
+    nsat = r.readU(5);
+    r.readU(1);
+    r.readU(3);
+  } else {
+    const tk = r.readU(27);
+    r.readU(1);
+    nsat = r.readU(5);
+    r.readU(1);
+    r.readU(3);
+    const gloDow = new Date(Date.now() + 3 * 36e5).getUTCDay();
+    epochMs = (gloDow << 27 | tk) >>> 0;
+  }
+  const observations = [];
+  for (let s = 0; s < nsat; s++) {
+    const satId = r.readU(6);
+    const code1 = r.readU(1);
+    const freqCh = glo ? r.readU(5) : 0;
+    const pr1raw = r.readU(gps ? 24 : 25);
+    const ppr1 = r.readS(20);
+    const lock1 = r.readU(7);
+    let amb = 0;
+    let cnr1 = 0;
+    if (ext) {
+      amb = r.readU(gps ? 8 : 7);
+      cnr1 = r.readU(8);
+    }
+    let code2 = 0;
+    let pr21 = 0;
+    let ppr2 = 0;
+    let lock2 = 0;
+    let cnr2 = 0;
+    if (dual) {
+      code2 = r.readU(2);
+      pr21 = r.readS(14);
+      ppr2 = r.readS(20);
+      lock2 = r.readU(7);
+      if (ext) cnr2 = r.readU(8);
+    }
+    if (gps && (satId < 1 || satId > 32)) continue;
+    if (glo && (satId < 1 || satId > 24)) continue;
+    const prunit = gps ? PRUNIT_GPS : PRUNIT_GLO;
+    const pr1 = pr1raw * 0.02 + amb * prunit;
+    const k = freqCh - 7;
+    const lamL1 = gps ? LAM_GPS_L1 : C_LIGHT / (GLO_F1_BASE + k * GLO_F1_STEP);
+    const lamL2 = gps ? LAM_GPS_L2 : C_LIGHT / (GLO_F2_BASE + k * GLO_F2_STEP);
+    const signals = [];
+    {
+      const sig = {
+        rinexCode: gps ? code1 ? "1W" : "1C" : code1 ? "1P" : "1C",
+        pseudorange: pr1,
+        wavelength: lamL1,
+        lockTime: legacyLockSec(lock1)
+      };
+      if (ppr1 !== PPR_INVALID) sig.phase = (pr1 + ppr1 * 5e-4) / lamL1;
+      if (ext && cnr1 > 0) sig.cn0 = cnr1 * 0.25;
+      signals.push(sig);
+    }
+    if (dual) {
+      const sig = {
+        rinexCode: gps ? code2 === 0 ? "2X" : "2W" : code2 === 0 ? "2C" : "2P",
+        wavelength: lamL2,
+        lockTime: legacyLockSec(lock2)
+      };
+      if (pr21 !== PR21_INVALID) sig.pseudorange = pr1 + pr21 * 0.02;
+      if (ppr2 !== PPR_INVALID) sig.phase = (pr1 + ppr2 * 5e-4) / lamL2;
+      if (ext && cnr2 > 0) sig.cn0 = cnr2 * 0.25;
+      signals.push(sig);
+    }
+    observations.push({
+      prn: `${system}${two2(satId)}`,
+      system,
+      signals
+    });
+  }
+  return { messageType: type, epochMs, system, observations };
+}
+
 // src/rtcm3/msm.ts
 var GPS_L1 = FREQ.G["1"];
 var GPS_L2 = FREQ.G["2"];
@@ -1097,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 {
@@ -1447,10 +1881,14 @@ export {
   Rtcm3Decoder,
   readString,
   decodeEphemeris,
+  decodeIgsSsr,
+  decodeLegacyObs,
   decodeMsmFull,
   msmEpochToDate,
   setGloFreqNumber,
   resetGloFreqCache,
+  decodeSsr,
+  isSsrMessage,
   createStationMeta,
   updateStationMeta,
   RTCM3_MESSAGE_NAMES,