@aeroswarm/weather 0.1.1 → 0.1.3

package/README.md

@@ -25,3 +25,16 @@ const snapshot = animationSnapshot({ observation: metar, forecast: taf, at: new
 
 METAR/TAF describe one point (the aerodrome), not a weather field, so the animation is a
 representation around the airport, not the real position of clouds or rain.
+
+## Gridded fields
+
+`grid.ts` turns wind and scalar data into grids, whatever the source:
+
+- `windGridFromStations(stations, { bounds, nx, ny })` interpolates the wind of several stations (IDW);
+  `stationWindsFromMetars(reports, locate)` takes them from decoded METARs.
+- `parseGrib2Json(json)` reads `grib2json` output (convert GRIB binary first with `grib2json` or
+  `wgrib2`); `parseGridJson(json)` reads a plain `{ bounds, nx, ny, u, v }` or `{ ..., values }`.
+- `stationWindsFromModelPoints(response)` reads multi-location model responses (Open-Meteo style).
+- `sampleWind`, `sampleScalar` (bilinear) and `blendWindGrids` (between two forecast hours).
+
+Row 0 is north, wind is in m/s (`u` east, `v` north), cells without data are NaN.

package/dist/grid.d.ts

@@ -0,0 +1,116 @@
+import type { Metar } from "./types.js";
+/**
+ * Medan cuaca berbentuk grid (angin, awan, hujan, suhu) dari sumber mana pun: interpolasi banyak
+ * stasiun METAR, JSON grid dari layanan luar, atau hasil konversi GRIB. Semua murni, tanpa jaringan.
+ *
+ * Konvensi: baris 0 = utara, kolom 0 = barat. Angin dalam m/s: `u` ke timur, `v` ke utara.
+ * Sel tanpa data berisi NaN.
+ */
+export interface GridBounds {
+    west: number;
+    south: number;
+    east: number;
+    north: number;
+}
+export interface ScalarGrid {
+    bounds: GridBounds;
+    nx: number;
+    ny: number;
+    values: Float32Array;
+}
+export interface WindGrid {
+    bounds: GridBounds;
+    nx: number;
+    ny: number;
+    u: Float32Array;
+    v: Float32Array;
+}
+export interface GeoPoint {
+    lat: number;
+    lon: number;
+}
+export interface StationWind extends GeoPoint {
+    speedKt: number;
+    /** Arah ASAL angin (derajat); null untuk VRB (dianggap tenang). */
+    directionDeg: number | null;
+}
+export interface StationValue extends GeoPoint {
+    value: number;
+}
+export interface InterpolationOptions {
+    bounds: GridBounds;
+    nx: number;
+    ny: number;
+    /** Pangkat bobot 1/d^p. Bawaan 2. */
+    power?: number;
+    /** Sel yang lebih jauh dari stasiun mana pun dari jarak ini (km) dibiarkan tanpa data. Bawaan 400. */
+    maxDistanceKm?: number;
+}
+export declare const speedOf: (u: number, v: number) => number;
+/** Arah ASAL angin (derajat, 0-360) dari komponen u/v. */
+export declare function directionFromUV(u: number, v: number): number;
+/** Komponen u/v (m/s) dari kecepatan dan arah ASAL angin. */
+export declare function uvFromSpeedDirection(speedMs: number, fromDeg: number): {
+    u: number;
+    v: number;
+};
+export declare function pointOf(bounds: GridBounds, nx: number, ny: number, col: number, row: number): GeoPoint;
+/**
+ * Interpolasi bobot-jarak-terbalik (IDW) nilai di banyak titik ke grid. Tepat di sebuah stasiun
+ * hasilnya persis nilai stasiun itu. Sel yang jauh dari semua stasiun diberi NaN.
+ */
+export declare function interpolateScalar(points: readonly StationValue[], options: InterpolationOptions): ScalarGrid;
+/** Grid angin dari angin di banyak stasiun (mis. METAR beberapa bandara). */
+export declare function windGridFromStations(stations: readonly StationWind[], options: InterpolationOptions): WindGrid;
+/** Nilai grid di (lon, lat) dengan interpolasi bilinear; null di luar grid atau tanpa data. */
+export declare function sampleScalar(grid: ScalarGrid, lon: number, lat: number): number | null;
+/** Angin (m/s) di (lon, lat); null di luar grid atau tanpa data. */
+export declare function sampleWind(grid: WindGrid, lon: number, lat: number): {
+    u: number;
+    v: number;
+} | null;
+/** Campur dua grid angin yang bentuknya sama (mis. dua jam prakiraan), `t` 0-1. */
+export declare function blendWindGrids(a: WindGrid, b: WindGrid, t: number): WindGrid;
+/**
+ * Baca JSON hasil `grib2json` (format yang juga dipakai leaflet-velocity): dua rekaman, komponen angin
+ * u (parameter 2) dan v (parameter 3). GRIB biner dikonversi dulu ke JSON ini (grib2json/wgrib2).
+ */
+export declare function parseGrib2Json(input: unknown): WindGrid;
+/**
+ * Angin stasiun dari beberapa METAR/SPECI. `locate` memberi koordinat dari kode ICAO (mis. dari
+ * tabel bandara milik aplikasi); stasiun tanpa koordinat dilewati.
+ */
+export declare function stationWindsFromMetars(reports: readonly Metar[], locate: (icao: string) => GeoPoint | undefined): StationWind[];
+/**
+ * Format JSON sederhana bawaan: `{ bounds: {west,south,east,north}, nx, ny, u, v }` untuk angin atau
+ * `{ bounds, nx, ny, values }` untuk skalar (awan, hujan, suhu). Larik boleh berisi null untuk "tanpa data".
+ */
+export declare function parseGridJson(input: unknown): WindGrid | ScalarGrid;
+declare const OPEN_METEO_SPEED_TO_MS: Record<string, number>;
+/**
+ * Angin dari respons banyak lokasi layanan model bergaya Open-Meteo: larik `{latitude, longitude,
+ * current|hourly: {wind_speed_10m, wind_direction_10m}}`. `hourIndex` memilih jam dari `hourly`.
+ * `speedUnit` = satuan kecepatan di respons (bawaan km/h, seperti Open-Meteo).
+ */
+export declare function stationWindsFromModelPoints(response: unknown, options?: {
+    hourIndex?: number;
+    speedUnit?: keyof typeof OPEN_METEO_SPEED_TO_MS;
+}): StationWind[];
+export type DecodedGrid = {
+    kind: "wind";
+    grid: WindGrid;
+} | {
+    kind: "scalar";
+    grid: ScalarGrid;
+};
+/**
+ * Kenali dan baca JSON medan cuaca dari sumber mana pun: `grib2json`, JSON grid bawaan (`u`/`v` atau
+ * `values`), atau respons banyak lokasi layanan model (diinterpolasi ke grid `fallback`, wajib
+ * diberikan untuk bentuk terakhir ini). Melempar galat yang menjelaskan bila bentuknya tidak dikenal.
+ */
+export declare function decodeGridJson(input: unknown, fallback?: InterpolationOptions & {
+    hourIndex?: number;
+    speedUnit?: keyof typeof OPEN_METEO_SPEED_TO_MS;
+}): DecodedGrid;
+export {};
+//# sourceMappingURL=grid.d.ts.map

package/dist/grid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"grid.d.ts","sourceRoot":"","sources":["../src/grid.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,YAAY,CAAC;AAExC;;;;;;GAMG;AAEH,MAAM,WAAW,UAAU;IACzB,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,MAAM,CAAC;IACd,IAAI,EAAE,MAAM,CAAC;IACb,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,UAAU;IACzB,MAAM,EAAE,UAAU,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,MAAM,EAAE,YAAY,CAAC;CACtB;AAED,MAAM,WAAW,QAAQ;IACvB,MAAM,EAAE,UAAU,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,CAAC,EAAE,YAAY,CAAC;IAChB,CAAC,EAAE,YAAY,CAAC;CACjB;AAED,MAAM,WAAW,QAAQ;IACvB,GAAG,EAAE,MAAM,CAAC;IACZ,GAAG,EAAE,MAAM,CAAC;CACb;AAED,MAAM,WAAW,WAAY,SAAQ,QAAQ;IAC3C,OAAO,EAAE,MAAM,CAAC;IAChB,mEAAmE;IACnE,YAAY,EAAE,MAAM,GAAG,IAAI,CAAC;CAC7B;AAED,MAAM,WAAW,YAAa,SAAQ,QAAQ;IAC5C,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,oBAAoB;IACnC,MAAM,EAAE,UAAU,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,qCAAqC;IACrC,KAAK,CAAC,EAAE,MAAM,CAAC;IACf,sGAAsG;IACtG,aAAa,CAAC,EAAE,MAAM,CAAC;CACxB;AAKD,eAAO,MAAM,OAAO,GAAI,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,MAA0B,CAAC;AAE1E,0DAA0D;AAC1D,wBAAgB,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE5D;AAED,6DAA6D;AAC7D,wBAAgB,oBAAoB,CAClC,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,GACd;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,CAG1B;AAED,wBAAgB,OAAO,CAAC,MAAM,EAAE,UAAU,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,QAAQ,CAKtG;AAQD;;;GAGG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,SAAS,YAAY,EAAE,EAC/B,OAAO,EAAE,oBAAoB,GAC5B,UAAU,CA2BZ;AAED,6EAA6E;AAC7E,wBAAgB,oBAAoB,CAClC,QAAQ,EAAE,SAAS,WAAW,EAAE,EAChC,OAAO,EAAE,oBAAoB,GAC5B,QAAQ,CAUV;AA2BD,+FAA+F;AAC/F,wBAAgB,YAAY,CAAC,IAAI,EAAE,UAAU,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM,GAAG,IAAI,CAEtF;AAED,oEAAoE;AACpE,wBAAgB,UAAU,CAAC,IAAI,EAAE,QAAQ,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAA;CAAE,GAAG,IAAI,CAIpG;AAED,mFAAmF;AACnF,wBAAgB,cAAc,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,GAAG,QAAQ,CAK5E;AAgCD;;;GAGG;AACH,wBAAgB,cAAc,CAAC,KAAK,EAAE,OAAO,GAAG,QAAQ,CAuCvD;AAED;;;GAGG;AACH,wBAAgB,sBAAsB,CACpC,OAAO,EAAE,SAAS,KAAK,EAAE,EACzB,MAAM,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,GAAG,SAAS,GAC7C,WAAW,EAAE,CAcf;AAED;;;GAGG;AACH,wBAAgB,aAAa,CAAC,KAAK,EAAE,OAAO,GAAG,QAAQ,GAAG,UAAU,CAiBnE;AAED,QAAA,MAAM,sBAAsB,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAKlD,CAAC;AAEF;;;;GAIG;AACH,wBAAgB,2BAA2B,CACzC,QAAQ,EAAE,OAAO,EACjB,OAAO,GAAE;IAAE,SAAS,CAAC,EAAE,MAAM,CAAC;IAAC,SAAS,CAAC,EAAE,MAAM,OAAO,sBAAsB,CAAA;CAAO,GACpF,WAAW,EAAE,CAuBf;AAED,MAAM,MAAM,WAAW,GACnB;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,QAAQ,CAAA;CAAE,GAChC;IAAE,IAAI,EAAE,QAAQ,CAAC;IAAC,IAAI,EAAE,UAAU,CAAA;CAAE,CAAC;AAEzC;;;;GAIG;AACH,wBAAgB,cAAc,CAC5B,KAAK,EAAE,OAAO,EACd,QAAQ,CAAC,EAAE,oBAAoB,GAAG;IAAE,SAAS,CAAC,EAAE,MAAM,CAAC;IAAC,SAAS,CAAC,EAAE,MAAM,OAAO,sBAAsB,CAAA;CAAE,GACxG,WAAW,CAiBb"}

package/dist/grid.test.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=grid.test.d.ts.map

package/dist/grid.test.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"grid.test.d.ts","sourceRoot":"","sources":["../src/grid.test.ts"],"names":[],"mappings":""}

package/dist/index.d.ts

@@ -2,8 +2,10 @@ import type { Report } from "./types.js";
 export * from "./aftn.js";
 export * from "./animation.js";
 export * from "./conditions.js";
+export * from "./grid.js";
 export * from "./metar.js";
 export * from "./taf.js";
+export * from "./store.js";
 export * from "./time.js";
 export * from "./types.js";
 /**

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AAEzC,cAAc,WAAW,CAAC;AAC1B,cAAc,gBAAgB,CAAC;AAC/B,cAAc,iBAAiB,CAAC;AAChC,cAAc,YAAY,CAAC;AAC3B,cAAc,UAAU,CAAC;AACzB,cAAc,WAAW,CAAC;AAC1B,cAAc,YAAY,CAAC;AAE3B;;;GAGG;AACH,wBAAgB,oBAAoB,CAClC,OAAO,EAAE,MAAM,EACf,SAAS,GAAE,IAAiB,GAC3B,MAAM,EAAE,CAUV"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAGA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AAEzC,cAAc,WAAW,CAAC;AAC1B,cAAc,gBAAgB,CAAC;AAC/B,cAAc,iBAAiB,CAAC;AAChC,cAAc,WAAW,CAAC;AAC1B,cAAc,YAAY,CAAC;AAC3B,cAAc,UAAU,CAAC;AACzB,cAAc,YAAY,CAAC;AAC3B,cAAc,WAAW,CAAC;AAC1B,cAAc,YAAY,CAAC;AAE3B;;;GAGG;AACH,wBAAgB,oBAAoB,CAClC,OAAO,EAAE,MAAM,EACf,SAAS,GAAE,IAAiB,GAC3B,MAAM,EAAE,CAUV"}

package/dist/index.js

@@ -559,6 +559,427 @@ function animationSnapshot(input) {
   return null;
 }
 
+// src/grid.ts
+var KT_TO_MS = 0.514444;
+var KM_PER_DEG_LAT = 111.32;
+var speedOf = (u, v) => Math.hypot(u, v);
+function directionFromUV(u, v) {
+  return (Math.atan2(-u, -v) * 180 / Math.PI + 360) % 360;
+}
+function uvFromSpeedDirection(speedMs, fromDeg) {
+  const rad = fromDeg * Math.PI / 180;
+  return { u: -speedMs * Math.sin(rad), v: -speedMs * Math.cos(rad) };
+}
+function pointOf(bounds, nx, ny, col, row) {
+  return {
+    lon: bounds.west + (bounds.east - bounds.west) * col / Math.max(1, nx - 1),
+    lat: bounds.north - (bounds.north - bounds.south) * row / Math.max(1, ny - 1)
+  };
+}
+function distanceKm(a, b) {
+  const dLat = (a.lat - b.lat) * KM_PER_DEG_LAT;
+  const dLon = (a.lon - b.lon) * KM_PER_DEG_LAT * Math.cos((a.lat + b.lat) / 2 * (Math.PI / 180));
+  return Math.hypot(dLat, dLon);
+}
+function interpolateScalar(points, options) {
+  const { bounds, nx, ny } = options;
+  const power = options.power ?? 2;
+  const maxKm = options.maxDistanceKm ?? 400;
+  const values = new Float32Array(nx * ny).fill(Number.NaN);
+  for (let row = 0; row < ny; row++) {
+    for (let col = 0; col < nx; col++) {
+      const cell = pointOf(bounds, nx, ny, col, row);
+      let weightSum = 0;
+      let sum = 0;
+      let exact = null;
+      for (const p of points) {
+        const d = distanceKm(cell, p);
+        if (d > maxKm) continue;
+        if (d < 0.5) {
+          exact = p.value;
+          break;
+        }
+        const w = 1 / d ** power;
+        weightSum += w;
+        sum += w * p.value;
+      }
+      values[row * nx + col] = exact ?? (weightSum > 0 ? sum / weightSum : Number.NaN);
+    }
+  }
+  return { bounds, nx, ny, values };
+}
+function windGridFromStations(stations, options) {
+  const components = stations.map((s) => {
+    const speedMs = s.speedKt * KT_TO_MS;
+    const { u: u2, v: v2 } = s.directionDeg === null ? { u: 0, v: 0 } : uvFromSpeedDirection(speedMs, s.directionDeg);
+    return { lat: s.lat, lon: s.lon, u: u2, v: v2 };
+  });
+  const u = interpolateScalar(components.map((c) => ({ lat: c.lat, lon: c.lon, value: c.u })), options);
+  const v = interpolateScalar(components.map((c) => ({ lat: c.lat, lon: c.lon, value: c.v })), options);
+  return { bounds: options.bounds, nx: options.nx, ny: options.ny, u: u.values, v: v.values };
+}
+function bilinear(values, bounds, nx, ny, lon, lat) {
+  if (lon < bounds.west || lon > bounds.east || lat < bounds.south || lat > bounds.north) return null;
+  const fx = (lon - bounds.west) / (bounds.east - bounds.west || 1) * (nx - 1);
+  const fy = (bounds.north - lat) / (bounds.north - bounds.south || 1) * (ny - 1);
+  const x0 = Math.min(nx - 1, Math.floor(fx));
+  const y0 = Math.min(ny - 1, Math.floor(fy));
+  const x1 = Math.min(nx - 1, x0 + 1);
+  const y1 = Math.min(ny - 1, y0 + 1);
+  const tx = fx - x0;
+  const ty = fy - y0;
+  const a = values[y0 * nx + x0];
+  const b = values[y0 * nx + x1];
+  const c = values[y1 * nx + x0];
+  const d = values[y1 * nx + x1];
+  const result = (a * (1 - tx) + b * tx) * (1 - ty) + (c * (1 - tx) + d * tx) * ty;
+  return Number.isNaN(result) ? null : result;
+}
+function sampleScalar(grid, lon, lat) {
+  return bilinear(grid.values, grid.bounds, grid.nx, grid.ny, lon, lat);
+}
+function sampleWind(grid, lon, lat) {
+  const u = bilinear(grid.u, grid.bounds, grid.nx, grid.ny, lon, lat);
+  const v = bilinear(grid.v, grid.bounds, grid.nx, grid.ny, lon, lat);
+  return u === null || v === null ? null : { u, v };
+}
+function blendWindGrids(a, b, t) {
+  if (a.nx !== b.nx || a.ny !== b.ny) throw new Error("Grid angin harus berukuran sama");
+  const mix = (x, y) => Float32Array.from(x, (value, i) => value + (y[i] - value) * t);
+  return { bounds: a.bounds, nx: a.nx, ny: a.ny, u: mix(a.u, b.u), v: mix(a.v, b.v) };
+}
+var isRecord = (value) => typeof value === "object" && value !== null;
+function isGrib2Json(input) {
+  return Array.isArray(input) && input.length >= 2 && input.every((r) => isRecord(r) && isRecord(r.header) && Array.isArray(r.data));
+}
+function parseGrib2Json(input) {
+  if (!isGrib2Json(input)) throw new Error("Bukan JSON grib2json: perlu larik rekaman {header, data}");
+  const uRecord = input.find((r) => r.header.parameterNumber === 2 && (r.header.parameterCategory ?? 2) === 2);
+  const vRecord = input.find((r) => r.header.parameterNumber === 3 && (r.header.parameterCategory ?? 2) === 2);
+  if (!uRecord || !vRecord) throw new Error("Komponen angin u (2) dan v (3) tidak ditemukan");
+  const { nx, ny, lo1, la1 } = uRecord.header;
+  if (uRecord.data.length !== nx * ny || vRecord.data.length !== nx * ny) {
+    throw new Error("Jumlah data tidak sama dengan nx * ny");
+  }
+  const dx = uRecord.header.dx ?? ((uRecord.header.lo2 ?? lo1) - lo1) / Math.max(1, nx - 1);
+  const dy = Math.abs(uRecord.header.dy ?? ((uRecord.header.la2 ?? la1) - la1) / Math.max(1, ny - 1));
+  const lo2 = uRecord.header.lo2 ?? lo1 + dx * (nx - 1);
+  const la2 = uRecord.header.la2 ?? la1 - dy * (ny - 1);
+  const norm = (lon) => lon > 180 ? lon - 360 : lon;
+  const north = Math.max(la1, la2);
+  const south = Math.min(la1, la2);
+  const flip = la1 < la2;
+  const orient = (data) => {
+    const out = new Float32Array(nx * ny);
+    for (let row = 0; row < ny; row++) {
+      const source = flip ? ny - 1 - row : row;
+      for (let col = 0; col < nx; col++) {
+        const value = data[source * nx + col];
+        out[row * nx + col] = typeof value === "number" ? value : Number.NaN;
+      }
+    }
+    return out;
+  };
+  return {
+    bounds: { west: norm(lo1), south, east: norm(lo2), north },
+    nx,
+    ny,
+    u: orient(uRecord.data),
+    v: orient(vRecord.data)
+  };
+}
+function stationWindsFromMetars(reports, locate) {
+  const out = [];
+  for (const report of reports) {
+    const place = locate(report.station);
+    if (!place) continue;
+    const { wind } = report.conditions;
+    out.push({
+      lat: place.lat,
+      lon: place.lon,
+      speedKt: wind.speedKt,
+      directionDeg: wind.variable ? null : wind.directionDeg
+    });
+  }
+  return out;
+}
+function parseGridJson(input) {
+  if (!isRecord(input) || !isRecord(input.bounds)) throw new Error("JSON grid tidak valid: bounds hilang");
+  const { west, south, east, north } = input.bounds;
+  const nx = Number(input.nx);
+  const ny = Number(input.ny);
+  if (![west, south, east, north].every((n) => typeof n === "number") || !(nx > 0) || !(ny > 0)) {
+    throw new Error("JSON grid tidak valid: bounds/nx/ny");
+  }
+  const bounds = { west, south, east, north };
+  const toArray = (value) => {
+    if (!Array.isArray(value) || value.length !== nx * ny) throw new Error("Ukuran larik tidak sama dengan nx * ny");
+    return Float32Array.from(value, (n) => typeof n === "number" ? n : Number.NaN);
+  };
+  if (input.u !== void 0 && input.v !== void 0) {
+    return { bounds, nx, ny, u: toArray(input.u), v: toArray(input.v) };
+  }
+  return { bounds, nx, ny, values: toArray(input.values) };
+}
+var OPEN_METEO_SPEED_TO_MS = {
+  "km/h": 1 / 3.6,
+  "m/s": 1,
+  kn: KT_TO_MS,
+  mph: 0.44704
+};
+function stationWindsFromModelPoints(response, options = {}) {
+  const list = Array.isArray(response) ? response : [response];
+  const toMs = OPEN_METEO_SPEED_TO_MS[options.speedUnit ?? "km/h"] ?? 1 / 3.6;
+  const out = [];
+  for (const item of list) {
+    if (!isRecord(item)) continue;
+    const lat = Number(item.latitude);
+    const lon = Number(item.longitude);
+    const source = isRecord(item.current) ? item.current : isRecord(item.hourly) ? item.hourly : null;
+    if (!source || Number.isNaN(lat) || Number.isNaN(lon)) continue;
+    const pick = (value) => Array.isArray(value) ? Number(value[options.hourIndex ?? 0]) : Number(value);
+    const speed = pick(source.wind_speed_10m);
+    const direction = pick(source.wind_direction_10m);
+    if (Number.isNaN(speed)) continue;
+    out.push({
+      lat,
+      lon,
+      speedKt: speed * toMs / KT_TO_MS,
+      directionDeg: Number.isNaN(direction) ? null : direction
+    });
+  }
+  return out;
+}
+function decodeGridJson(input, fallback) {
+  if (isGrib2Json(input)) return { kind: "wind", grid: parseGrib2Json(input) };
+  if (isRecord(input) && isRecord(input.bounds)) {
+    const grid = parseGridJson(input);
+    return "u" in grid ? { kind: "wind", grid } : { kind: "scalar", grid };
+  }
+  const items = Array.isArray(input) ? input : [input];
+  if (items.some((item) => isRecord(item) && item.latitude !== void 0 && item.longitude !== void 0)) {
+    if (!fallback) throw new Error("Respons banyak lokasi perlu bounds, nx, dan ny untuk dijadikan grid");
+    const stations = stationWindsFromModelPoints(input, fallback);
+    if (stations.length === 0) throw new Error("Tidak ada titik dengan wind_speed_10m pada respons");
+    return { kind: "wind", grid: windGridFromStations(stations, fallback) };
+  }
+  throw new Error("Format medan tidak dikenal: harus grib2json, JSON grid {bounds,nx,ny,u,v|values}, atau respons banyak lokasi");
+}
+
+// src/store.ts
+var MAX_OBSERVATIONS_PER_STATION = 48;
+var MAX_FORECASTS_PER_STATION = 8;
+var rainRate = (intensity) => intensity * intensity * 15;
+var hhmm = (d) => `${String(d.getUTCHours()).padStart(2, "0")}:${String(d.getUTCMinutes()).padStart(2, "0")}`;
+var WeatherStore = class {
+  constructor(locate) {
+    this.locate = locate;
+  }
+  locate;
+  observations = /* @__PURE__ */ new Map();
+  forecasts = /* @__PURE__ */ new Map();
+  listeners = /* @__PURE__ */ new Set();
+  /** Naik tiap data berubah; berguna untuk memicu render ulang. */
+  version = 0;
+  /** Terima pesan mentah (laporan polos, buletin, atau pesan AFTN; boleh banyak laporan). */
+  ingest(message, options = {}) {
+    const receivedAt = options.receivedAt ?? /* @__PURE__ */ new Date();
+    const manual = options.manual === true;
+    const result = { metars: 0, tafs: 0, skipped: 0, stations: [] };
+    const touched = /* @__PURE__ */ new Set();
+    for (const raw of extractReports(message)) {
+      if (raw.kind === "TAF") {
+        const taf = parseTaf(raw.text, receivedAt);
+        if (!taf) {
+          result.skipped++;
+          continue;
+        }
+        this.add(this.forecasts, taf.station, { report: taf, manual, receivedAt }, MAX_FORECASTS_PER_STATION);
+        touched.add(taf.station);
+        result.tafs++;
+      } else {
+        const metar = parseMetar(raw.text, receivedAt);
+        if (!metar) {
+          result.skipped++;
+          continue;
+        }
+        this.add(this.observations, metar.station, { report: metar, manual, receivedAt }, MAX_OBSERVATIONS_PER_STATION);
+        touched.add(metar.station);
+        result.metars++;
+      }
+    }
+    result.stations = [...touched];
+    if (result.metars + result.tafs > 0) this.changed();
+    return result;
+  }
+  /** Hapus semua data manual (seluruh stasiun, atau satu stasiun). */
+  clearManual(station) {
+    let removed = false;
+    for (const store of [this.observations, this.forecasts]) {
+      for (const [icao, list] of store) {
+        if (station && icao !== station) continue;
+        const kept = list.filter((s) => !s.manual);
+        if (kept.length !== list.length) removed = true;
+        store.set(icao, kept);
+      }
+    }
+    if (removed) this.changed();
+  }
+  /** Hapus semua data (otomatis dan manual). */
+  clear() {
+    this.observations.clear();
+    this.forecasts.clear();
+    this.changed();
+  }
+  hasManual() {
+    for (const store of [this.observations, this.forecasts]) {
+      for (const list of store.values()) if (list.some((s) => s.manual)) return true;
+    }
+    return false;
+  }
+  stations() {
+    return [.../* @__PURE__ */ new Set([...this.observations.keys(), ...this.forecasts.keys()])].filter(
+      (icao) => (this.observations.get(icao)?.length ?? 0) + (this.forecasts.get(icao)?.length ?? 0) > 0
+    );
+  }
+  onChange(listener) {
+    this.listeners.add(listener);
+    return () => this.listeners.delete(listener);
+  }
+  /** Cuaca satu stasiun pada waktu `at`; null bila tidak ada data. */
+  snapshotAt(station, at, options = {}) {
+    return animationSnapshot({
+      observation: this.pickObservation(station, at, options.maxObservationAgeMs),
+      forecast: this.pickForecast(station, at),
+      at,
+      maxObservationAgeMs: options.maxObservationAgeMs,
+      includeProbable: options.includeProbable
+    });
+  }
+  /** Rentang waktu dan pembagian periode (BASE/BECMG/TEMPO/FM) dari TAF, untuk timeline. */
+  timeline() {
+    const segments = [];
+    const observations = [];
+    let from = null;
+    let to = null;
+    const widen = (a, b) => {
+      from = from === null ? a.getTime() : Math.min(from, a.getTime());
+      to = to === null ? b.getTime() : Math.max(to, b.getTime());
+    };
+    for (const [station, list] of this.forecasts) {
+      const taf = this.latest(list, (t) => t.issuedAt);
+      if (!taf || taf.cancelled) continue;
+      widen(taf.validFrom, taf.validTo);
+      for (const period of taf.periods) {
+        if (period.change === "BASE") continue;
+        const label = period.change === "FM" ? `FM ${hhmm(period.from)}` : `${period.probability ? `PROB${period.probability} ` : ""}${period.change} ${hhmm(period.from)}-${hhmm(period.to)}`;
+        segments.push({ station, change: period.change, from: period.from, to: period.to, label });
+      }
+    }
+    for (const list of this.observations.values()) {
+      for (const source of list) {
+        observations.push(source.report.observedAt);
+        widen(source.report.observedAt, source.report.observedAt);
+      }
+    }
+    if (from === null || to === null) return null;
+    observations.sort((a, b) => a.getTime() - b.getTime());
+    return { from: new Date(from), to: new Date(to), segments, observations };
+  }
+  /**
+   * Medan angin, awan, hujan, dan suhu pada waktu `at`, dari semua stasiun yang punya data dan
+   * koordinat, diinterpolasi ke grid. Mengembalikan null bila tidak ada stasiun yang bisa dipakai.
+   */
+  fieldsAt(at, grid, options = {}) {
+    const winds = [];
+    const clouds = [];
+    const rains = [];
+    const temperatures = [];
+    const used = [];
+    for (const station of this.stations()) {
+      const place = this.locate(station);
+      if (!place) continue;
+      const snapshot = this.snapshotAt(station, at, options);
+      if (!snapshot) continue;
+      used.push(station);
+      const { state } = snapshot;
+      winds.push({
+        lat: place.lat,
+        lon: place.lon,
+        speedKt: state.wind.speedKt,
+        directionDeg: state.wind.fromDeg
+      });
+      clouds.push({ lat: place.lat, lon: place.lon, value: state.clouds.total * 100 });
+      rains.push({
+        lat: place.lat,
+        lon: place.lon,
+        value: state.precipitation ? rainRate(state.precipitation.intensity) : 0
+      });
+      const observed = this.pickObservation(station, at);
+      if (observed?.temperatureC !== null && observed?.temperatureC !== void 0) {
+        temperatures.push({ lat: place.lat, lon: place.lon, value: observed.temperatureC });
+      }
+    }
+    if (used.length === 0) return null;
+    return {
+      at,
+      stations: used,
+      wind: windGridFromStations(winds, grid),
+      cloud: interpolateScalar(clouds, grid),
+      rain: interpolateScalar(rains, grid),
+      temperature: temperatures.length > 0 ? interpolateScalar(temperatures, grid) : null
+    };
+  }
+  // ---- internal -------------------------------------------------------------
+  changed() {
+    this.version++;
+    for (const listener of this.listeners) listener();
+  }
+  add(store, station, source, limit) {
+    const list = store.get(station) ?? [];
+    const key = (s) => {
+      const report = s.report;
+      return (report.observedAt ?? report.issuedAt)?.getTime() ?? 0;
+    };
+    const next = list.filter((s) => !(s.manual === source.manual && key(s) === key(source)));
+    next.push(source);
+    next.sort((a, b) => key(a) - key(b));
+    store.set(station, next.slice(-limit));
+  }
+  latest(list, timeOf) {
+    const pool = list.some((s) => s.manual) ? list.filter((s) => s.manual) : list;
+    let best = null;
+    for (const source of pool) {
+      if (!best || timeOf(source.report).getTime() >= timeOf(best.report).getTime()) best = source;
+    }
+    return best?.report ?? null;
+  }
+  pickObservation(station, at, maxAgeMs = 3 * 36e5) {
+    const list = this.observations.get(station) ?? [];
+    const usable = (s) => {
+      const age = at.getTime() - s.report.observedAt.getTime();
+      return age >= 0 && age <= maxAgeMs;
+    };
+    const manual = list.filter((s) => s.manual && usable(s));
+    const pool = manual.length > 0 ? manual : list.filter(usable);
+    let best = null;
+    for (const source of pool) {
+      if (!best || source.report.observedAt.getTime() >= best.report.observedAt.getTime()) best = source;
+    }
+    if (best) return best.report;
+    const before = list.filter((s) => s.report.observedAt.getTime() <= at.getTime());
+    return before.length > 0 ? before[before.length - 1].report : null;
+  }
+  pickForecast(station, at) {
+    const list = this.forecasts.get(station) ?? [];
+    const covers = (s) => !s.report.cancelled && at.getTime() >= s.report.validFrom.getTime() && at.getTime() < s.report.validTo.getTime();
+    const manual = list.filter((s) => s.manual && covers(s));
+    const pool = manual.length > 0 ? manual : list.filter(covers);
+    if (pool.length === 0) return null;
+    return this.latest(pool, (t) => t.issuedAt);
+  }
+};
+
 // src/index.ts
 function decodeWeatherMessage(message, reference = /* @__PURE__ */ new Date()) {
   const reports = [];
@@ -569,16 +990,31 @@ function decodeWeatherMessage(message, reference = /* @__PURE__ */ new Date()) {
   return reports;
 }
 export {
+  WeatherStore,
   animationFromConditions,
   animationSnapshot,
   blendConditions,
+  blendWindGrids,
+  decodeGridJson,
   decodeWeatherMessage,
+  directionFromUV,
   emptyConditions,
   extractReports,
   forecastAt,
+  interpolateScalar,
   mergeConditions,
   parseConditionTokens,
+  parseGrib2Json,
+  parseGridJson,
   parseMetar,
   parseTaf,
-  resolveDayTime
+  pointOf,
+  resolveDayTime,
+  sampleScalar,
+  sampleWind,
+  speedOf,
+  stationWindsFromMetars,
+  stationWindsFromModelPoints,
+  uvFromSpeedDirection,
+  windGridFromStations
 };

package/dist/store.d.ts

@@ -0,0 +1,91 @@
+import { type AnimationSnapshot } from "./animation.js";
+import { type GeoPoint, type GridBounds, type ScalarGrid, type WindGrid } from "./grid.js";
+import type { TafChange } from "./types.js";
+export type Locate = (icao: string) => GeoPoint | undefined;
+export interface IngestOptions {
+    /** Waktu pesan diterima (menentukan bulan/tahun laporan). Bawaan: sekarang. */
+    receivedAt?: Date;
+    /** Data manual selalu menimpa data otomatis untuk stasiun dan waktu yang sama. */
+    manual?: boolean;
+}
+export interface IngestResult {
+    metars: number;
+    tafs: number;
+    /** Laporan yang dikenali kata kuncinya tetapi tidak terbaca. */
+    skipped: number;
+    stations: string[];
+}
+export interface TimelineSegment {
+    station: string;
+    change: TafChange;
+    from: Date;
+    to: Date;
+    label: string;
+}
+export interface Timeline {
+    from: Date;
+    to: Date;
+    segments: TimelineSegment[];
+    /** Waktu pengamatan (METAR/SPECI) yang ada. */
+    observations: Date[];
+}
+export interface WeatherFields {
+    at: Date;
+    /** Stasiun yang punya data pada waktu itu. */
+    stations: string[];
+    wind: WindGrid;
+    /** Tutupan awan 0-100. */
+    cloud: ScalarGrid;
+    /** Curah hujan (mm/jam), perkiraan dari intensitas pada laporan. */
+    rain: ScalarGrid;
+    /** Suhu (°C); null bila tidak ada pengamatan bersuhu. */
+    temperature: ScalarGrid | null;
+}
+export interface FieldGridOptions {
+    bounds: GridBounds;
+    nx: number;
+    ny: number;
+    maxDistanceKm?: number;
+}
+/**
+ * Menampung laporan cuaca (METAR, SPECI, TAF) dari mana pun, lalu menjawab "bagaimana cuacanya di
+ * waktu T": per stasiun maupun sebagai medan grid. Data manual selalu menimpa data otomatis.
+ */
+export declare class WeatherStore {
+    private readonly locate;
+    private readonly observations;
+    private readonly forecasts;
+    private readonly listeners;
+    /** Naik tiap data berubah; berguna untuk memicu render ulang. */
+    version: number;
+    constructor(locate: Locate);
+    /** Terima pesan mentah (laporan polos, buletin, atau pesan AFTN; boleh banyak laporan). */
+    ingest(message: string, options?: IngestOptions): IngestResult;
+    /** Hapus semua data manual (seluruh stasiun, atau satu stasiun). */
+    clearManual(station?: string): void;
+    /** Hapus semua data (otomatis dan manual). */
+    clear(): void;
+    hasManual(): boolean;
+    stations(): string[];
+    onChange(listener: () => void): () => void;
+    /** Cuaca satu stasiun pada waktu `at`; null bila tidak ada data. */
+    snapshotAt(station: string, at: Date, options?: {
+        maxObservationAgeMs?: number;
+        includeProbable?: boolean;
+    }): AnimationSnapshot | null;
+    /** Rentang waktu dan pembagian periode (BASE/BECMG/TEMPO/FM) dari TAF, untuk timeline. */
+    timeline(): Timeline | null;
+    /**
+     * Medan angin, awan, hujan, dan suhu pada waktu `at`, dari semua stasiun yang punya data dan
+     * koordinat, diinterpolasi ke grid. Mengembalikan null bila tidak ada stasiun yang bisa dipakai.
+     */
+    fieldsAt(at: Date, grid: FieldGridOptions, options?: {
+        includeProbable?: boolean;
+    }): WeatherFields | null;
+    private changed;
+    private add;
+    private latest;
+    private pickObservation;
+    private pickForecast;
+}
+//# sourceMappingURL=store.d.ts.map

package/dist/store.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"store.d.ts","sourceRoot":"","sources":["../src/store.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,KAAK,iBAAiB,EAAqB,MAAM,gBAAgB,CAAC;AAC3E,OAAO,EACL,KAAK,QAAQ,EACb,KAAK,UAAU,EAEf,KAAK,UAAU,EAGf,KAAK,QAAQ,EAEd,MAAM,WAAW,CAAC;AAGnB,OAAO,KAAK,EAAc,SAAS,EAAE,MAAM,YAAY,CAAC;AAExD,MAAM,MAAM,MAAM,GAAG,CAAC,IAAI,EAAE,MAAM,KAAK,QAAQ,GAAG,SAAS,CAAC;AAE5D,MAAM,WAAW,aAAa;IAC5B,+EAA+E;IAC/E,UAAU,CAAC,EAAE,IAAI,CAAC;IAClB,kFAAkF;IAClF,MAAM,CAAC,EAAE,OAAO,CAAC;CAClB;AAED,MAAM,WAAW,YAAY;IAC3B,MAAM,EAAE,MAAM,CAAC;IACf,IAAI,EAAE,MAAM,CAAC;IACb,gEAAgE;IAChE,OAAO,EAAE,MAAM,CAAC;IAChB,QAAQ,EAAE,MAAM,EAAE,CAAC;CACpB;AAED,MAAM,WAAW,eAAe;IAC9B,OAAO,EAAE,MAAM,CAAC;IAChB,MAAM,EAAE,SAAS,CAAC;IAClB,IAAI,EAAE,IAAI,CAAC;IACX,EAAE,EAAE,IAAI,CAAC;IACT,KAAK,EAAE,MAAM,CAAC;CACf;AAED,MAAM,WAAW,QAAQ;IACvB,IAAI,EAAE,IAAI,CAAC;IACX,EAAE,EAAE,IAAI,CAAC;IACT,QAAQ,EAAE,eAAe,EAAE,CAAC;IAC5B,+CAA+C;IAC/C,YAAY,EAAE,IAAI,EAAE,CAAC;CACtB;AAED,MAAM,WAAW,aAAa;IAC5B,EAAE,EAAE,IAAI,CAAC;IACT,8CAA8C;IAC9C,QAAQ,EAAE,MAAM,EAAE,CAAC;IACnB,IAAI,EAAE,QAAQ,CAAC;IACf,0BAA0B;IAC1B,KAAK,EAAE,UAAU,CAAC;IAClB,oEAAoE;IACpE,IAAI,EAAE,UAAU,CAAC;IACjB,yDAAyD;IACzD,WAAW,EAAE,UAAU,GAAG,IAAI,CAAC;CAChC;AAED,MAAM,WAAW,gBAAgB;IAC/B,MAAM,EAAE,UAAU,CAAC;IACnB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,aAAa,CAAC,EAAE,MAAM,CAAC;CACxB;AAiBD;;;GAGG;AACH,qBAAa,YAAY;IAOX,OAAO,CAAC,QAAQ,CAAC,MAAM;IANnC,OAAO,CAAC,QAAQ,CAAC,YAAY,CAAsC;IACnE,OAAO,CAAC,QAAQ,CAAC,SAAS,CAAoC;IAC9D,OAAO,CAAC,QAAQ,CAAC,SAAS,CAAyB;IACnD,iEAAiE;IACjE,OAAO,SAAK;gBAEiB,MAAM,EAAE,MAAM;IAE3C,2FAA2F;IAC3F,MAAM,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,GAAE,aAAkB,GAAG,YAAY;IAiClE,oEAAoE;IACpE,WAAW,CAAC,OAAO,CAAC,EAAE,MAAM,GAAG,IAAI;IAanC,8CAA8C;IAC9C,KAAK,IAAI,IAAI;IAMb,SAAS,IAAI,OAAO;IAOpB,QAAQ,IAAI,MAAM,EAAE;IAMpB,QAAQ,CAAC,QAAQ,EAAE,MAAM,IAAI,GAAG,MAAM,IAAI;IAK1C,oEAAoE;IACpE,UAAU,CACR,OAAO,EAAE,MAAM,EACf,EAAE,EAAE,IAAI,EACR,OAAO,GAAE;QAAE,mBAAmB,CAAC,EAAE,MAAM,CAAC;QAAC,eAAe,CAAC,EAAE,OAAO,CAAA;KAAO,GACxE,iBAAiB,GAAG,IAAI;IAU3B,0FAA0F;IAC1F,QAAQ,IAAI,QAAQ,GAAG,IAAI;IAkC3B;;;OAGG;IACH,QAAQ,CAAC,EAAE,EAAE,IAAI,EAAE,IAAI,EAAE,gBAAgB,EAAE,OAAO,GAAE;QAAE,eAAe,CAAC,EAAE,OAAO,CAAA;KAAO,GAAG,aAAa,GAAG,IAAI;IA6C7G,OAAO,CAAC,OAAO;IAKf,OAAO,CAAC,GAAG;IAkBX,OAAO,CAAC,MAAM;IAUd,OAAO,CAAC,eAAe;IAmBvB,OAAO,CAAC,YAAY;CASrB"}

package/dist/store.test.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=store.test.d.ts.map

package/dist/store.test.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"store.test.d.ts","sourceRoot":"","sources":["../src/store.test.ts"],"names":[],"mappings":""}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@aeroswarm/weather",
-  "version": "0.1.1",
+  "version": "0.1.3",
   "description": "Decode raw AFTN/METAR/SPECI/TAF messages and turn them into animation parameters.",
   "type": "module",
   "main": "./dist/index.js",