gnss-js 0.1.2 → 1.0.1

package/dist/nav-BAI1a9vK.d.cts

@@ -0,0 +1,108 @@
+/**
+ * RINEX 3.x/4.x navigation file parser.
+ * Parses GPS, Galileo, BeiDou (Keplerian), and GLONASS (state vector) ephemerides.
+ */
+interface NavHeader {
+    version: number;
+    type: string;
+    leapSeconds: number | null;
+    ionoCorrections: Record<string, number[]>;
+}
+/** GPS / Galileo / BeiDou / QZSS / NavIC Keplerian ephemeris. */
+interface KeplerEphemeris {
+    system: 'G' | 'E' | 'C' | 'J' | 'I';
+    /** Satellite PRN, e.g. "G14". */
+    prn: string;
+    /** Epoch of clock in seconds of GPS week. */
+    toc: number;
+    /** Epoch of clock as a Date. */
+    tocDate: Date;
+    /** SV clock bias in seconds. */
+    af0: number;
+    /** SV clock drift in s/s. */
+    af1: number;
+    /** SV clock drift rate in s/s². */
+    af2: number;
+    /** Issue of data (ephemeris). */
+    iode: number;
+    /** Amplitude of sine correction to orbit radius in meters. */
+    crs: number;
+    /** Mean motion difference in rad/s. */
+    deltaN: number;
+    /** Mean anomaly at reference time in radians. */
+    m0: number;
+    /** Amplitude of cosine correction to argument of latitude in radians. */
+    cuc: number;
+    /** Eccentricity (dimensionless). */
+    e: number;
+    /** Amplitude of sine correction to argument of latitude in radians. */
+    cus: number;
+    /** Square root of semi-major axis in m^(1/2). */
+    sqrtA: number;
+    /** Time of ephemeris in seconds of GPS week. */
+    toe: number;
+    /** Amplitude of cosine correction to inclination in radians. */
+    cic: number;
+    /** Right ascension of ascending node at reference time in radians. */
+    omega0: number;
+    /** Amplitude of sine correction to inclination in radians. */
+    cis: number;
+    /** Inclination at reference time in radians. */
+    i0: number;
+    /** Amplitude of cosine correction to orbit radius in meters. */
+    crc: number;
+    /** Argument of perigee in radians. */
+    omega: number;
+    /** Rate of right ascension in rad/s. */
+    omegaDot: number;
+    /** Rate of inclination in rad/s. */
+    idot: number;
+    /** GPS week number. */
+    week: number;
+    /** SV health flag. */
+    svHealth: number;
+    /** Total group delay in seconds. */
+    tgd: number;
+}
+/** GLONASS / SBAS state-vector ephemeris. Positions in km, velocities in km/s, accelerations in km/s². */
+interface GlonassEphemeris {
+    system: 'R' | 'S';
+    prn: string;
+    tocDate: Date;
+    /** Clock bias (seconds). */
+    tauN: number;
+    /** Relative frequency bias (dimensionless). */
+    gammaN: number;
+    /** Message frame time (seconds of UTC day). */
+    messageFrameTime: number;
+    /** X position in km (PZ-90). */
+    x: number;
+    /** X velocity in km/s. */
+    xDot: number;
+    /** X acceleration in km/s² (luni-solar). */
+    xAcc: number;
+    /** Y position in km (PZ-90). */
+    y: number;
+    /** Y velocity in km/s. */
+    yDot: number;
+    /** Y acceleration in km/s² (luni-solar). */
+    yAcc: number;
+    /** Z position in km (PZ-90). */
+    z: number;
+    /** Z velocity in km/s. */
+    zDot: number;
+    /** Z acceleration in km/s² (luni-solar). */
+    zAcc: number;
+    /** Health flag (0 = healthy). */
+    health: number;
+    /** GLONASS frequency channel number (-7 to +6). */
+    freqNum: number;
+}
+type Ephemeris = KeplerEphemeris | GlonassEphemeris;
+interface NavResult {
+    header: NavHeader;
+    ephemerides: Ephemeris[];
+}
+declare function parseNavFile(text: string): NavResult;
+
+export { type Ephemeris as E, type GlonassEphemeris as G, type KeplerEphemeris as K, type NavHeader as N, type NavResult as a, parseNavFile as p };

package/dist/nav-BAI1a9vK.d.ts

@@ -0,0 +1,108 @@
+/**
+ * RINEX 3.x/4.x navigation file parser.
+ * Parses GPS, Galileo, BeiDou (Keplerian), and GLONASS (state vector) ephemerides.
+ */
+interface NavHeader {
+    version: number;
+    type: string;
+    leapSeconds: number | null;
+    ionoCorrections: Record<string, number[]>;
+}
+/** GPS / Galileo / BeiDou / QZSS / NavIC Keplerian ephemeris. */
+interface KeplerEphemeris {
+    system: 'G' | 'E' | 'C' | 'J' | 'I';
+    /** Satellite PRN, e.g. "G14". */
+    prn: string;
+    /** Epoch of clock in seconds of GPS week. */
+    toc: number;
+    /** Epoch of clock as a Date. */
+    tocDate: Date;
+    /** SV clock bias in seconds. */
+    af0: number;
+    /** SV clock drift in s/s. */
+    af1: number;
+    /** SV clock drift rate in s/s². */
+    af2: number;
+    /** Issue of data (ephemeris). */
+    iode: number;
+    /** Amplitude of sine correction to orbit radius in meters. */
+    crs: number;
+    /** Mean motion difference in rad/s. */
+    deltaN: number;
+    /** Mean anomaly at reference time in radians. */
+    m0: number;
+    /** Amplitude of cosine correction to argument of latitude in radians. */
+    cuc: number;
+    /** Eccentricity (dimensionless). */
+    e: number;
+    /** Amplitude of sine correction to argument of latitude in radians. */
+    cus: number;
+    /** Square root of semi-major axis in m^(1/2). */
+    sqrtA: number;
+    /** Time of ephemeris in seconds of GPS week. */
+    toe: number;
+    /** Amplitude of cosine correction to inclination in radians. */
+    cic: number;
+    /** Right ascension of ascending node at reference time in radians. */
+    omega0: number;
+    /** Amplitude of sine correction to inclination in radians. */
+    cis: number;
+    /** Inclination at reference time in radians. */
+    i0: number;
+    /** Amplitude of cosine correction to orbit radius in meters. */
+    crc: number;
+    /** Argument of perigee in radians. */
+    omega: number;
+    /** Rate of right ascension in rad/s. */
+    omegaDot: number;
+    /** Rate of inclination in rad/s. */
+    idot: number;
+    /** GPS week number. */
+    week: number;
+    /** SV health flag. */
+    svHealth: number;
+    /** Total group delay in seconds. */
+    tgd: number;
+}
+/** GLONASS / SBAS state-vector ephemeris. Positions in km, velocities in km/s, accelerations in km/s². */
+interface GlonassEphemeris {
+    system: 'R' | 'S';
+    prn: string;
+    tocDate: Date;
+    /** Clock bias (seconds). */
+    tauN: number;
+    /** Relative frequency bias (dimensionless). */
+    gammaN: number;
+    /** Message frame time (seconds of UTC day). */
+    messageFrameTime: number;
+    /** X position in km (PZ-90). */
+    x: number;
+    /** X velocity in km/s. */
+    xDot: number;
+    /** X acceleration in km/s² (luni-solar). */
+    xAcc: number;
+    /** Y position in km (PZ-90). */
+    y: number;
+    /** Y velocity in km/s. */
+    yDot: number;
+    /** Y acceleration in km/s² (luni-solar). */
+    yAcc: number;
+    /** Z position in km (PZ-90). */
+    z: number;
+    /** Z velocity in km/s. */
+    zDot: number;
+    /** Z acceleration in km/s² (luni-solar). */
+    zAcc: number;
+    /** Health flag (0 = healthy). */
+    health: number;
+    /** GLONASS frequency channel number (-7 to +6). */
+    freqNum: number;
+}
+type Ephemeris = KeplerEphemeris | GlonassEphemeris;
+interface NavResult {
+    header: NavHeader;
+    ephemerides: Ephemeris[];
+}
+declare function parseNavFile(text: string): NavResult;
+
+export { type Ephemeris as E, type GlonassEphemeris as G, type KeplerEphemeris as K, type NavHeader as N, type NavResult as a, parseNavFile as p };

package/dist/nmea.cjs

@@ -0,0 +1,247 @@
+"use strict";
+var __defProp = Object.defineProperty;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
+var __getOwnPropNames = Object.getOwnPropertyNames;
+var __hasOwnProp = Object.prototype.hasOwnProperty;
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
+
+// src/nmea/index.ts
+var nmea_exports = {};
+__export(nmea_exports, {
+  computeStats: () => computeStats,
+  nmeaCoordToDecimal: () => nmeaCoordToDecimal,
+  parseNmeaFile: () => parseNmeaFile,
+  verifyChecksum: () => verifyChecksum
+});
+module.exports = __toCommonJS(nmea_exports);
+function nmeaCoordToDecimal(raw, direction) {
+  if (!raw || !direction) return NaN;
+  const isLon = direction === "E" || direction === "W";
+  const degLen = isLon ? 3 : 2;
+  const degrees = parseInt(raw.substring(0, degLen), 10);
+  const minutes = parseFloat(raw.substring(degLen));
+  const decimal = degrees + minutes / 60;
+  return direction === "S" || direction === "W" ? -decimal : decimal;
+}
+function verifyChecksum(sentence) {
+  const starIdx = sentence.indexOf("*");
+  if (starIdx === -1) return false;
+  const body = sentence.substring(sentence.indexOf("$") + 1, starIdx);
+  let xor = 0;
+  for (let i = 0; i < body.length; i++) xor ^= body.charCodeAt(i);
+  const expected = parseInt(sentence.substring(starIdx + 1, starIdx + 3), 16);
+  return xor === expected;
+}
+function buildDate(timeStr, dateStr) {
+  if (!timeStr || timeStr.length < 6) return null;
+  const h = parseInt(timeStr.substring(0, 2), 10);
+  const m = parseInt(timeStr.substring(2, 4), 10);
+  const s = parseFloat(timeStr.substring(4));
+  if (isNaN(h) || isNaN(m) || isNaN(s)) return null;
+  let year = 2e3, month = 0, day = 1;
+  if (dateStr && dateStr.length >= 6) {
+    day = parseInt(dateStr.substring(0, 2), 10);
+    month = parseInt(dateStr.substring(2, 4), 10) - 1;
+    const yy = parseInt(dateStr.substring(4, 6), 10);
+    year = yy < 80 ? 2e3 + yy : 1900 + yy;
+  }
+  const ms = Math.round(s % 1 * 1e3);
+  return new Date(Date.UTC(year, month, day, h, m, Math.floor(s), ms));
+}
+function parseGGA(fields) {
+  if (fields.length < 10) return null;
+  const quality = parseInt(fields[6], 10);
+  if (isNaN(quality) || quality === 0) return null;
+  const lat = nmeaCoordToDecimal(fields[2], fields[3]);
+  const lon = nmeaCoordToDecimal(fields[4], fields[5]);
+  if (isNaN(lat) || isNaN(lon)) return null;
+  return {
+    time: buildDate(fields[1]),
+    timeKey: fields[1],
+    lat,
+    lon,
+    alt: fields[9] ? parseFloat(fields[9]) : null,
+    satellites: fields[7] ? parseInt(fields[7], 10) : null,
+    fixQuality: quality
+  };
+}
+function parseRMC(fields) {
+  if (fields.length < 10) return null;
+  if (fields[2] !== "A") return null;
+  const lat = nmeaCoordToDecimal(fields[3], fields[4]);
+  const lon = nmeaCoordToDecimal(fields[5], fields[6]);
+  if (isNaN(lat) || isNaN(lon)) return null;
+  return {
+    time: buildDate(fields[1], fields[9]),
+    timeKey: fields[1],
+    lat,
+    lon,
+    speed: fields[7] ? parseFloat(fields[7]) : null,
+    course: fields[8] ? parseFloat(fields[8]) : null,
+    hasDate: true
+  };
+}
+function computeStats(fixes) {
+  const valid = fixes.filter((f) => f.fixQuality > 0);
+  const withTime = valid.filter((f) => f.time !== null);
+  const withSat = valid.filter((f) => f.satellites !== null).map((f) => f.satellites);
+  let startTime = null;
+  let endTime = null;
+  let duration = null;
+  if (withTime.length >= 2) {
+    startTime = withTime[0].time;
+    endTime = withTime[withTime.length - 1].time;
+    duration = (endTime.getTime() - startTime.getTime()) / 1e3;
+  } else if (withTime.length === 1) {
+    startTime = endTime = withTime[0].time;
+  }
+  let totalDistance = null;
+  if (valid.length >= 2) {
+    let dist = 0;
+    for (let i = 1; i < valid.length; i++) {
+      const prev = valid[i - 1], curr = valid[i];
+      dist += haversine(prev.lat, prev.lon, curr.lat, curr.lon);
+    }
+    totalDistance = dist;
+  }
+  const speeds = valid.filter((f) => f.speed !== null).map((f) => f.speed * 1.852);
+  const maxSpeed = speeds.length > 0 ? Math.max(...speeds) : null;
+  let cep = null;
+  let drms2 = null;
+  let hRms = null;
+  let vRms = null;
+  if (valid.length >= 2) {
+    const toRad = (d) => d * Math.PI / 180;
+    const R = 6378137;
+    const e2 = 0.00669437999014;
+    const ecef = valid.map((f) => {
+      const lat = toRad(f.lat), lon = toRad(f.lon), h = f.alt ?? 0;
+      const N = R / Math.sqrt(1 - e2 * Math.sin(lat) ** 2);
+      return [
+        (N + h) * Math.cos(lat) * Math.cos(lon),
+        (N + h) * Math.cos(lat) * Math.sin(lon),
+        ((1 - e2) * N + h) * Math.sin(lat)
+      ];
+    });
+    const n = ecef.length;
+    const mx = ecef.reduce((s, p) => s + p[0], 0) / n;
+    const my = ecef.reduce((s, p) => s + p[1], 0) / n;
+    const mz = ecef.reduce((s, p) => s + p[2], 0) / n;
+    const refLon = Math.atan2(my, mx);
+    const refP = Math.sqrt(mx * mx + my * my);
+    const refLat = Math.atan2(mz, refP * (1 - e2));
+    const sinLat = Math.sin(refLat), cosLat = Math.cos(refLat);
+    const sinLon = Math.sin(refLon), cosLon = Math.cos(refLon);
+    const enu = ecef.map((p) => {
+      const dx = p[0] - mx, dy = p[1] - my, dz = p[2] - mz;
+      return {
+        e: -sinLon * dx + cosLon * dy,
+        n: -cosLon * sinLat * dx - sinLon * sinLat * dy + cosLat * dz,
+        u: cosLon * cosLat * dx + sinLon * cosLat * dy + sinLat * dz
+      };
+    });
+    const hDist = enu.map((p) => Math.sqrt(p.e * p.e + p.n * p.n));
+    hDist.sort((a, b) => a - b);
+    const medIdx = Math.floor(hDist.length * 0.5);
+    cep = hDist[medIdx];
+    const hSqSum = enu.reduce((s, p) => s + p.e * p.e + p.n * p.n, 0);
+    hRms = Math.sqrt(hSqSum / n);
+    drms2 = 2 * hRms;
+    const vSqSum = enu.reduce((s, p) => s + p.u * p.u, 0);
+    vRms = Math.sqrt(vSqSum / n);
+  }
+  return {
+    totalFixes: fixes.length,
+    validFixes: valid.length,
+    duration,
+    startTime,
+    endTime,
+    totalDistance,
+    maxSpeed,
+    cep,
+    drms2,
+    hRms,
+    vRms,
+    avgSatellites: withSat.length > 0 ? withSat.reduce((a, b) => a + b, 0) / withSat.length : null
+  };
+}
+function haversine(lat1, lon1, lat2, lon2) {
+  const toRad = (d) => d * Math.PI / 180;
+  const dLat = toRad(lat2 - lat1);
+  const dLon = toRad(lon2 - lon1);
+  const a = Math.sin(dLat / 2) ** 2 + Math.cos(toRad(lat1)) * Math.cos(toRad(lat2)) * Math.sin(dLon / 2) ** 2;
+  return 6371e3 * 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
+}
+function parseNmeaFile(content) {
+  const lines = content.split(/\r?\n/);
+  const fixMap = /* @__PURE__ */ new Map();
+  const fixOrder = [];
+  let lastDate;
+  for (const raw of lines) {
+    const line = raw.trim();
+    if (!line.startsWith("$")) continue;
+    if (line.includes("*") && !verifyChecksum(line)) continue;
+    const fields = line.split("*")[0].split(",");
+    const sentenceId = fields[0].substring(3);
+    let partial = null;
+    if (sentenceId === "GGA") partial = parseGGA(fields);
+    else if (sentenceId === "RMC") partial = parseRMC(fields);
+    if (!partial) continue;
+    const key = partial.timeKey || `line-${fixOrder.length}`;
+    if (partial.hasDate && partial.time) {
+      const d = partial.time;
+      lastDate = `${String(d.getUTCDate()).padStart(2, "0")}${String(d.getUTCMonth() + 1).padStart(2, "0")}${String(d.getUTCFullYear() % 100).padStart(2, "0")}`;
+    }
+    const existing = fixMap.get(key);
+    if (existing) {
+      if (partial.alt !== void 0 && partial.alt !== null)
+        existing.alt = partial.alt;
+      if (partial.satellites !== void 0 && partial.satellites !== null)
+        existing.satellites = partial.satellites;
+      if (partial.fixQuality !== void 0)
+        existing.fixQuality = partial.fixQuality;
+      if (partial.speed !== void 0 && partial.speed !== null)
+        existing.speed = partial.speed;
+      if (partial.course !== void 0 && partial.course !== null)
+        existing.course = partial.course;
+      if (partial.hasDate && partial.time) existing.time = partial.time;
+    } else {
+      let time = partial.time;
+      if (!partial.hasDate && time && lastDate) {
+        time = buildDate(partial.timeKey, lastDate);
+      }
+      fixMap.set(key, {
+        time,
+        lat: partial.lat ?? 0,
+        lon: partial.lon ?? 0,
+        alt: partial.alt ?? null,
+        satellites: partial.satellites ?? null,
+        fixQuality: partial.fixQuality ?? 0,
+        speed: partial.speed ?? null,
+        course: partial.course ?? null
+      });
+      fixOrder.push(key);
+    }
+  }
+  const fixes = fixOrder.map((k) => fixMap.get(k));
+  return { fixes, stats: computeStats(fixes) };
+}
+// Annotate the CommonJS export names for ESM import in node:
+0 && (module.exports = {
+  computeStats,
+  nmeaCoordToDecimal,
+  parseNmeaFile,
+  verifyChecksum
+});

package/dist/nmea.d.cts

@@ -0,0 +1,37 @@
+interface NmeaFix {
+    time: Date | null;
+    lat: number;
+    lon: number;
+    alt: number | null;
+    satellites: number | null;
+    fixQuality: number;
+    speed: number | null;
+    course: number | null;
+}
+interface NmeaStats {
+    totalFixes: number;
+    validFixes: number;
+    duration: number | null;
+    startTime: Date | null;
+    endTime: Date | null;
+    totalDistance: number | null;
+    maxSpeed: number | null;
+    cep: number | null;
+    drms2: number | null;
+    hRms: number | null;
+    vRms: number | null;
+    avgSatellites: number | null;
+}
+interface NmeaTrack {
+    fixes: NmeaFix[];
+    stats: NmeaStats;
+}
+/** Convert NMEA latitude/longitude (ddmm.mmmm / dddmm.mmmm) to decimal degrees. */
+declare function nmeaCoordToDecimal(raw: string, direction: string): number;
+/** Verify XOR checksum between '$' and '*'. */
+declare function verifyChecksum(sentence: string): boolean;
+declare function computeStats(fixes: NmeaFix[]): NmeaStats;
+/** Parse an NMEA log file and return a track with fixes and statistics. */
+declare function parseNmeaFile(content: string): NmeaTrack;
+
+export { type NmeaFix, type NmeaStats, type NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum };

package/dist/nmea.d.ts

@@ -0,0 +1,37 @@
+interface NmeaFix {
+    time: Date | null;
+    lat: number;
+    lon: number;
+    alt: number | null;
+    satellites: number | null;
+    fixQuality: number;
+    speed: number | null;
+    course: number | null;
+}
+interface NmeaStats {
+    totalFixes: number;
+    validFixes: number;
+    duration: number | null;
+    startTime: Date | null;
+    endTime: Date | null;
+    totalDistance: number | null;
+    maxSpeed: number | null;
+    cep: number | null;
+    drms2: number | null;
+    hRms: number | null;
+    vRms: number | null;
+    avgSatellites: number | null;
+}
+interface NmeaTrack {
+    fixes: NmeaFix[];
+    stats: NmeaStats;
+}
+/** Convert NMEA latitude/longitude (ddmm.mmmm / dddmm.mmmm) to decimal degrees. */
+declare function nmeaCoordToDecimal(raw: string, direction: string): number;
+/** Verify XOR checksum between '$' and '*'. */
+declare function verifyChecksum(sentence: string): boolean;
+declare function computeStats(fixes: NmeaFix[]): NmeaStats;
+/** Parse an NMEA log file and return a track with fixes and statistics. */
+declare function parseNmeaFile(content: string): NmeaTrack;
+
+export { type NmeaFix, type NmeaStats, type NmeaTrack, computeStats, nmeaCoordToDecimal, parseNmeaFile, verifyChecksum };

package/dist/nmea.js

@@ -0,0 +1,12 @@
+import {
+  computeStats,
+  nmeaCoordToDecimal,
+  parseNmeaFile,
+  verifyChecksum
+} from "./chunk-YDMYO3YK.js";
+export {
+  computeStats,
+  nmeaCoordToDecimal,
+  parseNmeaFile,
+  verifyChecksum
+};