@aguacerowx/mapsgl 0.0.57 → 0.0.58

package/src/nexrad/radarFrameGpuMatch.ts

@@ -1,111 +1,111 @@
-/**
- * Frame transforms applied in {@link MapboxRadarLayer} before GPU upload.
- * Mouse readouts must use the **same** gate/ray ordering as the texture or samples drift from pixels.
- */
-import type { DecodedRadarFrame } from './nexradArchiveCache.js';
-
-/** Shortest distance between two headings on [0, 360), in degrees. */
-function angularDistanceDeg(a: number, b: number): number {
-    let d = Math.abs(a - b) % 360;
-    if (d > 180) d = 360 - d;
-    return d;
-}
-
-/**
- * For Level-III layout products (KDP/N0H): remap incoming rays into fixed canonical azimuth bins
- * (same as MapboxRadarLayer `sortFrameToCanonicalBins`).
- */
-export function canonicalBinsRadarFrame(frame: DecodedRadarFrame): DecodedRadarFrame {
-    const nRays = frame.nRays;
-    const nGates = frame.nGates;
-    if (nRays <= 0 || nGates <= 0) return frame;
-    if (frame.rayBoundariesDeg.length < nRays + 1) return frame;
-
-    const degPerBin = 360 / nRays;
-    const bytesPerRay = nGates * 2;
-
-    const centers = new Float32Array(nRays);
-    for (let r = 0; r < nRays; r++) {
-        const lower = frame.rayBoundariesDeg[r]!;
-        const upper = frame.rayBoundariesDeg[r + 1]!;
-        const center = (lower + upper) * 0.5;
-        centers[r] = ((center % 360) + 360) % 360;
-    }
-
-    const canonicalGateData = new Uint8Array(nRays * bytesPerRay);
-    for (let i = 0; i < canonicalGateData.length; i += 2) {
-        canonicalGateData[i] = 128;
-    }
-
-    const used = new Array<boolean>(nRays).fill(false);
-    for (let bin = 0; bin < nRays; bin++) {
-        const targetDeg = ((((bin + 0.5) * degPerBin) % 360) + 360) % 360;
-        let bestR = -1;
-        let bestDist = Infinity;
-        for (let r = 0; r < nRays; r++) {
-            if (used[r]) continue;
-            const dist = angularDistanceDeg(centers[r], targetDeg);
-            if (dist < bestDist) {
-                bestDist = dist;
-                bestR = r;
-            }
-        }
-        if (bestR >= 0) {
-            used[bestR] = true;
-            canonicalGateData.set(
-                frame.gateData.subarray(bestR * bytesPerRay, (bestR + 1) * bytesPerRay),
-                bin * bytesPerRay,
-            );
-        }
-    }
-
-    const canonicalBoundaries = new Float32Array(nRays + 1);
-    for (let i = 0; i <= nRays; i++) {
-        canonicalBoundaries[i] = i * degPerBin;
-    }
-
-    return { ...frame, gateData: canonicalGateData, rayBoundariesDeg: canonicalBoundaries };
-}
-
-/** Sort rays by ascending start azimuth — same as MapboxRadarLayer `sortFrameByAzimuth`. */
-export function sortRadarFrameByAzimuth(frame: DecodedRadarFrame): DecodedRadarFrame {
-    const nRays = frame.nRays;
-    const nGates = frame.nGates;
-    const bytesPerRay = nGates * 2;
-
-    const rayOrder = Array.from({ length: nRays }, (_, i) => i).sort(
-        (a, b) => frame.rayBoundariesDeg[a] - frame.rayBoundariesDeg[b],
-    );
-
-    const sortedGateData = new Uint8Array(nRays * bytesPerRay);
-    const sortedBoundaries = new Float32Array(nRays + 1);
-
-    for (let newIdx = 0; newIdx < nRays; newIdx++) {
-        const oldIdx = rayOrder[newIdx]!;
-        sortedGateData.set(
-            frame.gateData.subarray(oldIdx * bytesPerRay, (oldIdx + 1) * bytesPerRay),
-            newIdx * bytesPerRay,
-        );
-        sortedBoundaries[newIdx] = frame.rayBoundariesDeg[oldIdx]!;
-    }
-    sortedBoundaries[nRays] =
-        frame.rayBoundariesDeg[rayOrder[nRays - 1]!]! +
-        (frame.rayBoundariesDeg[rayOrder[1]!]! - frame.rayBoundariesDeg[rayOrder[0]!]!);
-
-    return { ...frame, gateData: sortedGateData, rayBoundariesDeg: sortedBoundaries };
-}
-
-export type GpuReadoutFrameOptions = {
-    geometryLayoutKey?: string | null;
-};
-
-/**
- * Match {@link MapboxRadarLayer} `uploadFrame` gate/ray ordering for CPU readout.
- */
-export function prepareRadarFrameForGpuReadout(
-    frame: DecodedRadarFrame,
-    options?: GpuReadoutFrameOptions,
-): DecodedRadarFrame {
-    const hasLayoutKey = options?.geometryLayoutKey != null && options.geometryLayoutKey !== '';
-    return hasLayoutKey ? canonicalBinsRadarFrame(frame) : sortRadarFrameByAzimuth(frame);
-}
+/**
+ * Frame transforms applied in {@link MapboxRadarLayer} before GPU upload.
+ * Mouse readouts must use the **same** gate/ray ordering as the texture or samples drift from pixels.
+ */
+import type { DecodedRadarFrame } from './nexradArchiveCache.js';
+
+/** Shortest distance between two headings on [0, 360), in degrees. */
+function angularDistanceDeg(a: number, b: number): number {
+    let d = Math.abs(a - b) % 360;
+    if (d > 180) d = 360 - d;
+    return d;
+}
+
+/**
+ * For Level-III layout products (KDP/N0H): remap incoming rays into fixed canonical azimuth bins
+ * (same as MapboxRadarLayer `sortFrameToCanonicalBins`).
+ */
+export function canonicalBinsRadarFrame(frame: DecodedRadarFrame): DecodedRadarFrame {
+    const nRays = frame.nRays;
+    const nGates = frame.nGates;
+    if (nRays <= 0 || nGates <= 0) return frame;
+    if (frame.rayBoundariesDeg.length < nRays + 1) return frame;
+
+    const degPerBin = 360 / nRays;
+    const bytesPerRay = nGates * 2;
+
+    const centers = new Float32Array(nRays);
+    for (let r = 0; r < nRays; r++) {
+        const lower = frame.rayBoundariesDeg[r]!;
+        const upper = frame.rayBoundariesDeg[r + 1]!;
+        const center = (lower + upper) * 0.5;
+        centers[r] = ((center % 360) + 360) % 360;
+    }
+
+    const canonicalGateData = new Uint8Array(nRays * bytesPerRay);
+    for (let i = 0; i < canonicalGateData.length; i += 2) {
+        canonicalGateData[i] = 128;
+    }
+
+    const used = new Array<boolean>(nRays).fill(false);
+    for (let bin = 0; bin < nRays; bin++) {
+        const targetDeg = ((((bin + 0.5) * degPerBin) % 360) + 360) % 360;
+        let bestR = -1;
+        let bestDist = Infinity;
+        for (let r = 0; r < nRays; r++) {
+            if (used[r]) continue;
+            const dist = angularDistanceDeg(centers[r], targetDeg);
+            if (dist < bestDist) {
+                bestDist = dist;
+                bestR = r;
+            }
+        }
+        if (bestR >= 0) {
+            used[bestR] = true;
+            canonicalGateData.set(
+                frame.gateData.subarray(bestR * bytesPerRay, (bestR + 1) * bytesPerRay),
+                bin * bytesPerRay,
+            );
+        }
+    }
+
+    const canonicalBoundaries = new Float32Array(nRays + 1);
+    for (let i = 0; i <= nRays; i++) {
+        canonicalBoundaries[i] = i * degPerBin;
+    }
+
+    return { ...frame, gateData: canonicalGateData, rayBoundariesDeg: canonicalBoundaries };
+}
+
+/** Sort rays by ascending start azimuth — same as MapboxRadarLayer `sortFrameByAzimuth`. */
+export function sortRadarFrameByAzimuth(frame: DecodedRadarFrame): DecodedRadarFrame {
+    const nRays = frame.nRays;
+    const nGates = frame.nGates;
+    const bytesPerRay = nGates * 2;
+
+    const rayOrder = Array.from({ length: nRays }, (_, i) => i).sort(
+        (a, b) => frame.rayBoundariesDeg[a] - frame.rayBoundariesDeg[b],
+    );
+
+    const sortedGateData = new Uint8Array(nRays * bytesPerRay);
+    const sortedBoundaries = new Float32Array(nRays + 1);
+
+    for (let newIdx = 0; newIdx < nRays; newIdx++) {
+        const oldIdx = rayOrder[newIdx]!;
+        sortedGateData.set(
+            frame.gateData.subarray(oldIdx * bytesPerRay, (oldIdx + 1) * bytesPerRay),
+            newIdx * bytesPerRay,
+        );
+        sortedBoundaries[newIdx] = frame.rayBoundariesDeg[oldIdx]!;
+    }
+    sortedBoundaries[nRays] =
+        frame.rayBoundariesDeg[rayOrder[nRays - 1]!]! +
+        (frame.rayBoundariesDeg[rayOrder[1]!]! - frame.rayBoundariesDeg[rayOrder[0]!]!);
+
+    return { ...frame, gateData: sortedGateData, rayBoundariesDeg: sortedBoundaries };
+}
+
+export type GpuReadoutFrameOptions = {
+    geometryLayoutKey?: string | null;
+};
+
+/**
+ * Match {@link MapboxRadarLayer} `uploadFrame` gate/ray ordering for CPU readout.
+ */
+export function prepareRadarFrameForGpuReadout(
+    frame: DecodedRadarFrame,
+    options?: GpuReadoutFrameOptions,
+): DecodedRadarFrame {
+    const hasLayoutKey = options?.geometryLayoutKey != null && options.geometryLayoutKey !== '';
+    return hasLayoutKey ? canonicalBinsRadarFrame(frame) : sortRadarFrameByAzimuth(frame);
+}

package/src/satelliteDefaultColormaps.js

@@ -1,37 +1,37 @@
-/**
- * Default GOES brightness-temperature colormaps for WebGL satellite, aligned with
- * aguacero-frontend `src/lib/defaultColormaps.ts` (`SatelliteLongwaveIR` and `SatelliteWV`, K/°C units).
- *
- * Format: flat array `[norm0, '#hex', norm1, '#hex', ...]` where each **norm** is 0–100
- * (cold → hot in texture space; see `SatelliteShaderManager._generateTextureFromColormapData` and
- * frontend `normToKelvin` / `kelvinToNorm` in `utilityFunctions.ts`).
- */
-
-/** Longwave IR (C13–C16-style) — matches `SatelliteLongwaveIR.units.K.colormap` in the frontend. */
-export const DEFAULT_SATELLITE_LONGWAVE_IR_COLORMAP = [
-    5, '#000000',
-    51, '#ffffff',
-    52, '#00ffff',
-    58, '#0033cc',
-    64, '#33ff00',
-    73, '#ffee00',
-    78, '#ff0000',
-    82, '#660000',
-    88, '#808080',
-    90, '#cc66ff',
-    96, '#9933cc',
-    100, '#9933cc',
-];
-
-/** Upper / mid / lower water vapor (C08–C10) — matches `SatelliteWV.units.K.colormap` in the frontend. */
-export const DEFAULT_SATELLITE_WATER_VAPOR_COLORMAP = [
-    10, '#cece55',
-    20, '#ce5555',
-    30, '#262626',
-    45, '#ffffff',
-    55, '#55cece',
-    70, '#5555ce',
-    80, '#4d3366',
-    90, '#e673e6',
-    100, '#efe6eb',
-];
+/**
+ * Default GOES brightness-temperature colormaps for WebGL satellite, aligned with
+ * aguacero-frontend `src/lib/defaultColormaps.ts` (`SatelliteLongwaveIR` and `SatelliteWV`, K/°C units).
+ *
+ * Format: flat array `[norm0, '#hex', norm1, '#hex', ...]` where each **norm** is 0–100
+ * (cold → hot in texture space; see `SatelliteShaderManager._generateTextureFromColormapData` and
+ * frontend `normToKelvin` / `kelvinToNorm` in `utilityFunctions.ts`).
+ */
+
+/** Longwave IR (C13–C16-style) — matches `SatelliteLongwaveIR.units.K.colormap` in the frontend. */
+export const DEFAULT_SATELLITE_LONGWAVE_IR_COLORMAP = [
+    5, '#000000',
+    51, '#ffffff',
+    52, '#00ffff',
+    58, '#0033cc',
+    64, '#33ff00',
+    73, '#ffee00',
+    78, '#ff0000',
+    82, '#660000',
+    88, '#808080',
+    90, '#cc66ff',
+    96, '#9933cc',
+    100, '#9933cc',
+];
+
+/** Upper / mid / lower water vapor (C08–C10) — matches `SatelliteWV.units.K.colormap` in the frontend. */
+export const DEFAULT_SATELLITE_WATER_VAPOR_COLORMAP = [
+    10, '#cece55',
+    20, '#ce5555',
+    30, '#262626',
+    45, '#ffffff',
+    55, '#55cece',
+    70, '#5555ce',
+    80, '#4d3366',
+    90, '#e673e6',
+    100, '#efe6eb',
+];