gnss-js 1.64.0 → 1.66.0

package/dist/{chunk-TWAC2BZF.js → chunk-2CMAFEHI.js}

@@ -6,7 +6,7 @@ import {
   computeSatPosition,
   ecefToAzEl,
   selectEphemeris
-} from "./chunk-MHGNKSSH.js";
+} from "./chunk-JYWYQ7DW.js";
 import {
   ecefToGeodetic,
   getAer,

package/dist/{chunk-MRC6LIPM.js → chunk-GH7MYXON.js}

@@ -5,7 +5,7 @@ import {
   bdsSubframeParityOk,
   galInavPageCrcOk,
   testGloString
-} from "./chunk-FA2UTLEQ.js";
+} from "./chunk-3KSDL3CO.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk

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

@@ -1531,6 +1531,30 @@ function isSsrMessage(messageType) {
   return messageType in SSR_KIND;
 }
 
+// src/rtcm3/system-params.ts
+function decodeSystemParams(frame) {
+  if (frame.messageType !== 1013) return null;
+  const r = new BitReader(frame.payload);
+  r.readU(12);
+  const referenceStationId = r.readU(12);
+  const mjd = r.readU(16);
+  const secondsOfDay = r.readU(17);
+  const nm = r.readU(5);
+  const leapSeconds = r.readU(8);
+  const messages = [];
+  for (let i = 0; i < nm; i++) {
+    messages.push({
+      messageId: r.readU(12),
+      // DF055
+      sync: r.readU(1) === 1,
+      // DF056
+      intervalS: r.readU(16) * 0.1
+      // DF057 (0.1 s)
+    });
+  }
+  return { referenceStationId, mjd, secondsOfDay, leapSeconds, messages };
+}
+
 // src/rtcm3/station.ts
 function createStationMeta() {
   return {
@@ -1544,7 +1568,8 @@ function createStationMeta() {
     receiverType: null,
     receiverFirmware: null,
     receiverSerial: null,
-    description: null
+    description: null,
+    physicalRefStationId: null
   };
 }
 function decodeStationARP(payload, msgType, meta) {
@@ -1567,6 +1592,18 @@ function decodeStationARP(payload, msgType, meta) {
     meta.antennaHeight = r.readU(16) * 1e-4;
   }
 }
+function decodeStation1032(payload, meta) {
+  if (payload.length < 19) return;
+  const r = new BitReader(payload);
+  r.skip(12);
+  meta.stationId = r.readU(12);
+  meta.physicalRefStationId = r.readU(12);
+  meta.itrf = r.readU(6);
+  const x = r.readS(38) * 1e-4;
+  const y = r.readS(38) * 1e-4;
+  const z = r.readS(38) * 1e-4;
+  if (x !== 0 || y !== 0 || z !== 0) meta.position = [x, y, z];
+}
 function decodeAntennaDescriptor(payload, msgType, meta) {
   if (payload.length < 5) return;
   const r = new BitReader(payload);
@@ -1647,6 +1684,9 @@ function updateStationMeta(meta, frame) {
       case 1006:
         decodeStationARP(frame.payload, frame.messageType, meta);
         return true;
+      case 1032:
+        decodeStation1032(frame.payload, meta);
+        return true;
       case 1007:
       case 1008:
         decodeAntennaDescriptor(frame.payload, frame.messageType, meta);
@@ -1889,6 +1929,7 @@ export {
   resetGloFreqCache,
   decodeSsr,
   isSsrMessage,
+  decodeSystemParams,
   createStationMeta,
   updateStationMeta,
   RTCM3_MESSAGE_NAMES,

package/dist/chunk-S624C5AC.js

@@ -0,0 +1,269 @@
+// src/rtcm2/index.ts
+var RTCM2_PREAMBLE = 102;
+var HAMMING = [
+  3139384448,
+  1569692224,
+  2932329728,
+  1466164864,
+  1806824256,
+  2340063680
+];
+function getbitu(buff, pos, len) {
+  let v = 0;
+  for (let i = pos; i < pos + len; i++)
+    v = v * 2 + (buff[i >> 3] >> 7 - (i & 7) & 1);
+  return v >>> 0;
+}
+function getbits(buff, pos, len) {
+  const v = getbitu(buff, pos, len);
+  return v & 1 << len - 1 ? v - 2 ** len : v;
+}
+function decodeWord(w) {
+  let word = w >>> 0;
+  if (word & 1073741824) word = (word ^ 1073741760) >>> 0;
+  let parity = 0;
+  for (let i = 0; i < 6; i++) {
+    parity <<= 1;
+    let x = (word & HAMMING[i]) >>> 6;
+    while (x) {
+      parity ^= x & 1;
+      x >>>= 1;
+    }
+  }
+  if (parity >>> 0 !== (word & 63)) return null;
+  return [word >>> 22 & 255, word >>> 14 & 255, word >>> 6 & 255];
+}
+var Rtcm2Decoder = class {
+  word = 0;
+  nbyte = 0;
+  nbit = 0;
+  len = 0;
+  buff = new Uint8Array(256);
+  /** Reset the framer (e.g. after a stream reconnect). */
+  reset() {
+    this.word = 0;
+    this.nbyte = 0;
+    this.nbit = 0;
+    this.len = 0;
+  }
+  decode(bytes) {
+    const frames = [];
+    for (const raw of bytes) {
+      if ((raw & 192) !== 64) continue;
+      let data = raw;
+      for (let i = 0; i < 6; i++, data >>= 1) {
+        this.word = (this.word << 1 | data & 1) >>> 0;
+        if (this.nbyte === 0) {
+          let preamb = this.word >>> 22 & 255;
+          if (this.word & 1073741824) preamb ^= 255;
+          if (preamb !== RTCM2_PREAMBLE) continue;
+          const d2 = decodeWord(this.word);
+          if (!d2) continue;
+          this.buff[0] = d2[0];
+          this.buff[1] = d2[1];
+          this.buff[2] = d2[2];
+          this.nbyte = 3;
+          this.nbit = 0;
+          continue;
+        }
+        if (++this.nbit < 30) continue;
+        this.nbit = 0;
+        const d = decodeWord(this.word);
+        if (!d) {
+          this.nbyte = 0;
+          this.word &= 3;
+          continue;
+        }
+        this.buff[this.nbyte] = d[0];
+        this.buff[this.nbyte + 1] = d[1];
+        this.buff[this.nbyte + 2] = d[2];
+        this.nbyte += 3;
+        if (this.nbyte === 6) this.len = (this.buff[5] >> 3) * 3 + 6;
+        if (this.nbyte < this.len) continue;
+        frames.push(this.readHeader(this.buff.slice(0, this.len)));
+        this.nbyte = 0;
+        this.word &= 3;
+      }
+    }
+    return frames;
+  }
+  readHeader(payload) {
+    return {
+      messageType: getbitu(payload, 8, 6),
+      stationId: getbitu(payload, 14, 10),
+      zcountSec: getbitu(payload, 24, 13) * 0.6,
+      sequenceNo: getbitu(payload, 37, 3),
+      lengthWords: payload[5] >> 3,
+      health: getbitu(payload, 45, 3),
+      payload
+    };
+  }
+};
+function rtcm2Station(frame) {
+  const p = frame.payload;
+  const nbit = frame.lengthWords * 24;
+  if (frame.messageType === 3) {
+    if (nbit < 96) return null;
+    return {
+      position: [
+        getbits(p, 48, 32) * 0.01,
+        getbits(p, 80, 32) * 0.01,
+        getbits(p, 112, 32) * 0.01
+      ]
+    };
+  }
+  if (frame.messageType === 22) {
+    if (nbit < 24) return null;
+    const l1 = [
+      getbits(p, 48, 8) / 25600,
+      getbits(p, 56, 8) / 25600,
+      getbits(p, 64, 8) / 25600
+    ];
+    let height;
+    if (nbit >= 48) {
+      const noh = getbits(p, 48 + 24 + 5, 1);
+      height = noh ? 0 : getbitu(p, 48 + 24 + 6, 18) / 25600;
+    }
+    return { l1Offset: l1, height };
+  }
+  return null;
+}
+function rtcm2Observation(frame) {
+  const t = frame.messageType;
+  if (t !== 18 && t !== 19) return null;
+  const p = frame.payload;
+  const nbit = frame.lengthWords * 24 + 48;
+  let i = 48;
+  if (i + 24 > nbit) return null;
+  const fRaw = getbitu(p, i, 2);
+  i += 4;
+  const microSec = getbitu(p, i, 20);
+  i += 20;
+  if (fRaw & 1) return null;
+  const freq = fRaw >> 1;
+  let sync = false;
+  const records = [];
+  while (i + 48 <= nbit) {
+    sync = getbitu(p, i, 1) === 1;
+    i += 1;
+    const code = getbitu(p, i, 1);
+    i += 1;
+    const sys = getbitu(p, i, 1);
+    i += 1;
+    let prn = getbitu(p, i, 5);
+    i += 5;
+    let loss;
+    if (t === 18) {
+      i += 3;
+      loss = getbitu(p, i, 5);
+      i += 5;
+    } else {
+      i += 8;
+    }
+    const val = t === 18 ? getbits(p, i, 32) : getbitu(p, i, 32);
+    i += 32;
+    if (prn === 0) prn = 32;
+    const system = sys ? "R" : "G";
+    const rec = {
+      prn: `${system}${String(prn).padStart(2, "0")}`,
+      system,
+      freq,
+      code: code ? "P" : "C/A",
+      lossOfLock: loss
+    };
+    if (t === 18)
+      rec.phaseCycles = -val / 256;
+    else rec.pseudorangeM = val * 0.02;
+    records.push(rec);
+  }
+  return { freq, sync, microSec, records };
+}
+function rtcm2Dgps(frame) {
+  const t = frame.messageType;
+  if (t !== 1 && t !== 9 && t !== 31) return null;
+  const system = t === 31 ? "R" : "G";
+  const p = frame.payload;
+  const nbit = frame.lengthWords * 24 + 48;
+  let i = 48;
+  const corrections = [];
+  while (i + 40 <= nbit) {
+    const fact = getbitu(p, i, 1);
+    i += 1;
+    const udre = getbitu(p, i, 2);
+    i += 2;
+    let prn = getbitu(p, i, 5);
+    i += 5;
+    const prc = getbits(p, i, 16);
+    i += 16;
+    const rrc = getbits(p, i, 8);
+    i += 8;
+    const iod = getbits(p, i, 8);
+    i += 8;
+    if (prn === 0) prn = 32;
+    corrections.push({
+      prn: `${system}${String(prn).padStart(2, "0")}`,
+      prc: prc * (fact ? 0.32 : 0.02),
+      rrc: rrc * (fact ? 0.032 : 2e-3),
+      iod,
+      udre
+    });
+  }
+  return { system, corrections };
+}
+function rtcm2Time(frame) {
+  if (frame.messageType !== 14) return null;
+  const p = frame.payload;
+  if (frame.lengthWords * 24 < 24) return null;
+  return {
+    zcountSec: getbitu(p, 24, 13) * 0.6,
+    week: getbitu(p, 48, 10),
+    hour: getbitu(p, 58, 8),
+    leapSeconds: getbitu(p, 66, 6)
+  };
+}
+function rtcm2Text(frame) {
+  if (frame.messageType !== 16) return null;
+  const p = frame.payload;
+  const nbit = frame.lengthWords * 24 + 48;
+  let s = "";
+  for (let i = 48; i + 8 <= nbit; i += 8) {
+    const c = getbitu(p, i, 8);
+    if (c === 0) break;
+    s += String.fromCharCode(c);
+  }
+  return s;
+}
+var RTCM2_MESSAGE_NAMES = {
+  1: "GPS differential corrections",
+  3: "Reference station parameters (ARP)",
+  9: "GPS partial correction set",
+  14: "GPS time of week",
+  16: "GPS special message (text)",
+  17: "GPS ephemerides",
+  18: "RTK uncorrected carrier phase",
+  19: "RTK uncorrected pseudorange",
+  22: "Extended reference station parameters",
+  23: "Antenna type definition",
+  24: "Antenna reference point (ARP)",
+  31: "GLONASS differential corrections",
+  32: "GLONASS reference station parameters",
+  34: "GLONASS partial correction set",
+  36: "GLONASS special message"
+};
+function looksLikeRtcm2(bytes) {
+  const n = Math.min(bytes.length, 4096);
+  if (n < 12) return false;
+  for (let i = 0; i < n; i++) if ((bytes[i] & 192) !== 64) return false;
+  return new Rtcm2Decoder().decode(bytes.subarray(0, n)).length >= 2;
+}
+
+export {
+  Rtcm2Decoder,
+  rtcm2Station,
+  rtcm2Observation,
+  rtcm2Dgps,
+  rtcm2Time,
+  rtcm2Text,
+  RTCM2_MESSAGE_NAMES,
+  looksLikeRtcm2
+};

package/dist/{chunk-UHLUAU6H.js → chunk-WW3DKCDN.js}

@@ -7,11 +7,14 @@ import {
   galFnavPageCrcOk,
   galInavPageCrcOk,
   testGloString
-} from "./chunk-FA2UTLEQ.js";
+} from "./chunk-3KSDL3CO.js";
 import {
   decodeSbasGeoNav,
   sbasCrcOk
 } from "./chunk-WQX354KD.js";
+import {
+  getGpsLeap
+} from "./chunk-HVXYFUCB.js";
 import {
   geodeticToEcef
 } from "./chunk-37QNKGTC.js";
@@ -26,9 +29,6 @@ import {
   getBitU,
   setBitU
 } from "./chunk-V2D5KUKH.js";
-import {
-  getGpsLeap
-} from "./chunk-HVXYFUCB.js";
 
 // src/sbf/frame.ts
 var CRC_TABLE = (() => {

package/dist/constants.js

@@ -1,10 +1,4 @@
 import "./chunk-5S5MONFK.js";
-import {
-  WGS84_ECCENTRICITY_SQUARED,
-  WGS84_FLATTENING,
-  WGS84_SEMI_MAJOR_AXIS,
-  WGS84_SEMI_MINOR_AXIS
-} from "./chunk-6FAL6P4G.js";
 import {
   DAYS_MJD2000_MJD,
   DAYS_MJD_JULIAN,
@@ -30,6 +24,12 @@ import {
   START_TAI_TIME,
   START_UNIX_TIME
 } from "./chunk-LEEU5OIO.js";
+import {
+  WGS84_ECCENTRICITY_SQUARED,
+  WGS84_FLATTENING,
+  WGS84_SEMI_MAJOR_AXIS,
+  WGS84_SEMI_MINOR_AXIS
+} from "./chunk-6FAL6P4G.js";
 import {
   ARC_GAP_FACTOR,
   BAND_LABELS,