@aeroswarm/weather 0.1.1 → 0.1.2

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,6 +2,7 @@ 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 "./time.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,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,213 @@ 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/index.ts
 function decodeWeatherMessage(message, reference = /* @__PURE__ */ new Date()) {
   const reports = [];
@@ -572,13 +779,27 @@ export {
   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/package.json

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