cog-tiler-wasm 0.3.3 → 0.3.5

package/cog-tiler.js

@@ -26,6 +26,7 @@ import proj4 from "proj4";
 import * as GeoTIFF from "geotiff";
 import geokeysToProj4 from "geotiff-geokeys-to-proj4";
 import { sampleWindowBilinear } from "./sampling.js";
+import { computeStats } from "./statistics.js";
 
 const OS = 20037508.342789244; // Web Mercator half-extent (m)
 const TILE = 256; // output tile size (px)
@@ -65,44 +66,6 @@ function dtypeOf(lv) {
   return (f || "uint") + b;
 }
 
-/** Min/max/mean/std/count/valid_percent/percentiles/histogram for a band buffer. */
-function computeStats(buf, nodata) {
-  let min = Infinity, max = -Infinity, sum = 0, sumsq = 0;
-  const valid = [];
-  for (let i = 0; i < buf.length; i++) {
-    const v = buf[i];
-    if (Number.isNaN(v)) continue;
-    if (nodata != null && v === nodata) continue;
-    if (v < min) min = v;
-    if (v > max) max = v;
-    sum += v;
-    sumsq += v * v;
-    valid.push(v);
-  }
-  const count = valid.length;
-  if (count === 0) return { count: 0, valid_percent: 0 };
-  const mean = sum / count;
-  const std = Math.sqrt(Math.max(0, sumsq / count - mean * mean));
-  valid.sort((a, b) => a - b);
-  const pct = (p) => valid[Math.min(count - 1, Math.floor((p / 100) * count))];
-  const bins = 10, span = max - min || 1, hist = new Array(bins).fill(0);
-  for (const v of valid) {
-    let k = Math.floor(((v - min) / span) * bins);
-    if (k >= bins) k = bins - 1;
-    if (k < 0) k = 0;
-    hist[k]++;
-  }
-  const edges = Array.from({ length: bins + 1 }, (_, i) => min + (span * i) / bins);
-  return {
-    min, max, mean, std, count,
-    valid_percent: (count / buf.length) * 100,
-    median: pct(50),
-    percentile_2: pct(2),
-    percentile_98: pct(98),
-    histogram: [hist, edges],
-  };
-}
-
 /** Transfer curve for a rescaled value in 0..1: stretch then gamma (mirrors the
  *  Rust `transfer`; reversal is colormap-only so it's not applied to RGB). */
 function transferCurve(t, stretch, gamma) {

package/package.json

@@ -5,7 +5,7 @@
     "Qiusheng Wu <giswqs@gmail.com>"
   ],
   "description": "Serverless Cloud Optimized GeoTIFF (COG) dynamic tiling in WebAssembly. TiTiler-style XYZ tiles, no backend.",
-  "version": "0.3.3",
+  "version": "0.3.5",
   "license": "MIT",
   "repository": {
     "type": "git",
@@ -17,6 +17,7 @@
     "cog_tiler_wasm.d.ts",
     "cog-tiler.js",
     "sampling.js",
+    "statistics.js",
     "cog-tiler.d.ts",
     "README.md",
     "LICENSE"
@@ -38,7 +39,7 @@
     }
   },
   "peerDependencies": {
-    "whitebox-wasm": "^0.4.1",
+    "whitebox-wasm": "^0.4.1 || ^0.5.1",
     "proj4": "^2.15.0",
     "geotiff": "^2.1.0 || ^3.0.0",
     "geotiff-geokeys-to-proj4": "^2024.4.13 || ^2026.8.16"

package/statistics.js

@@ -0,0 +1,44 @@
+/** Number of histogram buckets exposed through the TiTiler-style statistics API. */
+export const HISTOGRAM_BINS = 128;
+
+/** Min/max/mean/std/count/valid_percent/percentiles/histogram for a band buffer. */
+export function computeStats(buf, nodata) {
+  let min = Infinity, max = -Infinity, sum = 0, sumsq = 0;
+  const valid = [];
+  for (let i = 0; i < buf.length; i++) {
+    const v = buf[i];
+    if (Number.isNaN(v)) continue;
+    if (nodata != null && v === nodata) continue;
+    if (v < min) min = v;
+    if (v > max) max = v;
+    sum += v;
+    sumsq += v * v;
+    valid.push(v);
+  }
+  const count = valid.length;
+  if (count === 0) return { count: 0, valid_percent: 0 };
+  const mean = sum / count;
+  const std = Math.sqrt(Math.max(0, sumsq / count - mean * mean));
+  valid.sort((a, b) => a - b);
+  const pct = (p) => valid[Math.min(count - 1, Math.floor((p / 100) * count))];
+  const span = max - min || 1;
+  const histogram = new Array(HISTOGRAM_BINS).fill(0);
+  for (const v of valid) {
+    let bin = Math.floor(((v - min) / span) * HISTOGRAM_BINS);
+    if (bin >= HISTOGRAM_BINS) bin = HISTOGRAM_BINS - 1;
+    if (bin < 0) bin = 0;
+    histogram[bin]++;
+  }
+  const edges = Array.from(
+    { length: HISTOGRAM_BINS + 1 },
+    (_, i) => min + (span * i) / HISTOGRAM_BINS,
+  );
+  return {
+    min, max, mean, std, count,
+    valid_percent: (count / buf.length) * 100,
+    median: pct(50),
+    percentile_2: pct(2),
+    percentile_98: pct(98),
+    histogram: [histogram, edges],
+  };
+}