@photivra/engine 1.4.0 → 2.0.0

package/dist/processing/denoise.js

@@ -0,0 +1,214 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * Image-domain denoising, sharpening and synthetic film grain on linear or developed RGB:
+ * - Bilateral filter (Tomasi & Manduchi, ICCV 1998): Gaussian spatial and Gaussian range weights,
+ *   range distance Euclidean over the three channels of the supplied data.
+ * - Non-local means (Buades, Coll & Morel, CVPR 2005): weights exp(−max(d² − 2σ², 0)/h²), where d² is
+ *   the mean squared difference of (2f+1)² patches over all channels. This is the uniform-patch form
+ *   with the known-noise offset from Buades et al., IPOL 2011; σ = 0 gives the 2005 form with a
+ *   uniform patch kernel. The centre pixel gets the largest weight in its search window.
+ * - Unsharp masking: out = in + amount·(in − G_σ ⊗ in) where |in − G_σ ⊗ in| > threshold.
+ * - Synthetic film grain: an aesthetic, deterministic, signal-dependent texture. It is NOT a physical
+ *   film model (no silver-halide grain statistics, dye clouds or film characteristic curve).
+ * Noise in linear RAW-derived data depends on the signal level, but these filters use one global set of
+ * parameters.
+ * @see docs/PROCESSING.md for equations, sources and limits.
+ */
+import { InvalidConfigurationError } from "../core/configuration-error.js";
+import { createDeterministicUniformStream, deriveNativeSiteStreamSeedUint32, sampleDeterministicStandardNormal } from "../sensor/stochastic-seed.js";
+import { convolveRows, convolveSeparable, gaussianKernel, reflect } from "./image-ops.js";
+import { createProcessingTask, finiteRows, maximumWorkingBytes, optionalFiniteIn, optionalIntegerIn, processingProvenance, requireEnum, requireFiniteIn, requireImage, requireOptions } from "./processing-task.js";
+/** Brute-force bilateral filter (window radius ⌈3σ_s⌉, reflect-101 borders). */
+export function createBilateralFilterTask(input, provider) {
+    const r = requireOptions(input, ["image", "spatialSigmaPixels", "rangeSigma", "maximumWorkingBytes"], "createBilateralFilterTask input");
+    const image = requireImage(r.image, "image"), sigmaS = requireFiniteIn(r.spatialSigmaPixels, "spatialSigmaPixels", 0, 64, true);
+    const sigmaR = requireFiniteIn(r.rangeSigma, "rangeSigma", 0, 1e30, true), W = image.width, H = image.height, rad = Math.ceil(3 * sigmaS);
+    const spatial = new Float64Array((2 * rad + 1) ** 2);
+    for (let dy = -rad; dy <= rad; dy++)
+        for (let dx = -rad; dx <= rad; dx++)
+            spatial[(dy + rad) * (2 * rad + 1) + dx + rad] = Math.exp(-(dx * dx + dy * dy) / (2 * sigmaS * sigmaS));
+    return createProcessingTask({
+        provider, totalTileCount: 2 * H, estimatedWorkingBytes: W * H * 12, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),
+        body: function* () {
+            yield* finiteRows(image.data, W, H, 3, "image");
+            const d = image.data, out = new Float32Array(3 * W * H), k = 1 / (2 * sigmaR * sigmaR);
+            for (let y = 0; y < H; y++) {
+                for (let x = 0; x < W; x++) {
+                    const i = 3 * (y * W + x), r0 = d[i], g0 = d[i + 1], b0 = d[i + 2];
+                    let ws = 0, sr = 0, sg = 0, sb = 0;
+                    for (let dy = -rad; dy <= rad; dy++) {
+                        const row = reflect(y + dy, H) * W;
+                        for (let dx = -rad; dx <= rad; dx++) {
+                            const j = 3 * (row + reflect(x + dx, W)), R = d[j], G = d[j + 1], B = d[j + 2];
+                            const w = spatial[(dy + rad) * (2 * rad + 1) + dx + rad] * Math.exp(-((R - r0) ** 2 + (G - g0) ** 2 + (B - b0) ** 2) * k);
+                            ws += w;
+                            sr += w * R;
+                            sg += w * G;
+                            sb += w * B;
+                        }
+                    }
+                    out[i] = sr / ws;
+                    out[i + 1] = sg / ws;
+                    out[i + 2] = sb / ws;
+                }
+                yield;
+            }
+            return { width: W, height: H, data: out, provenance: processingProvenance("denoise-bilateral-tomasi-manduchi-1998", [
+                    "Gaussian spatial (sigma " + sigmaS + " px, radius " + rad + ") times Gaussian range weight on Euclidean RGB distance (sigma " + sigmaR + ").",
+                    "Range distance in the supplied data's units, not CIELAB; reflect-101 borders; one global parameter set (noise is not modelled per pixel)."
+                ]) };
+        }
+    });
+}
+/** Non-local means (Buades, Coll & Morel 2005), offset-major with box-filtered patch distances. */
+export function createNonLocalMeansTask(input, provider) {
+    const r = requireOptions(input, ["image", "filteringParameter", "noiseStandardDeviation", "patchRadius", "searchRadius", "maximumWorkingBytes"], "createNonLocalMeansTask input");
+    const image = requireImage(r.image, "image"), h = requireFiniteIn(r.filteringParameter, "filteringParameter", 0, 1e30, true);
+    const sigma = optionalFiniteIn(r.noiseStandardDeviation, 0, "noiseStandardDeviation", 0, 1e30);
+    const f = optionalIntegerIn(r.patchRadius, 1, "patchRadius", 0, 16), sr = optionalIntegerIn(r.searchRadius, 5, "searchRadius", 1, 32);
+    const W = image.width, H = image.height, n = W * H, offsets = (2 * sr + 1) ** 2 - 1;
+    const box = new Float64Array(2 * f + 1).fill(1 / (2 * f + 1));
+    return createProcessingTask({
+        provider, totalTileCount: H + offsets * (2 * H + convolveRows(H)) + H, estimatedWorkingBytes: n * (12 + 4 * 5 + 12),
+        maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),
+        body: function* () {
+            yield* finiteRows(image.data, W, H, 3, "image");
+            const d = image.data, acc = new Float32Array(3 * n), wsum = new Float32Array(n), wmax = new Float32Array(n);
+            const diff = new Float32Array(n), tmp = new Float32Array(n), dist = new Float32Array(n), inv = 1 / (h * h), offset2 = 2 * sigma * sigma;
+            for (let oy = -sr; oy <= sr; oy++)
+                for (let ox = -sr; ox <= sr; ox++) {
+                    if (ox === 0 && oy === 0)
+                        continue;
+                    for (let y = 0; y < H; y++) {
+                        const yy = reflect(y + oy, H) * W;
+                        for (let x = 0; x < W; x++) {
+                            const i = 3 * (y * W + x), j = 3 * (yy + reflect(x + ox, W));
+                            diff[y * W + x] = ((d[i] - d[j]) ** 2 + (d[i + 1] - d[j + 1]) ** 2 + (d[i + 2] - d[j + 2]) ** 2) / 3;
+                        }
+                        yield;
+                    }
+                    yield* convolveSeparable(diff, W, H, 1, box, tmp, dist);
+                    for (let y = 0; y < H; y++) {
+                        const yy = reflect(y + oy, H) * W;
+                        for (let x = 0; x < W; x++) {
+                            const p = y * W + x, j = 3 * (yy + reflect(x + ox, W)), w = Math.exp(-Math.max(dist[p] - offset2, 0) * inv);
+                            acc[3 * p] = acc[3 * p] + w * d[j];
+                            acc[3 * p + 1] = acc[3 * p + 1] + w * d[j + 1];
+                            acc[3 * p + 2] = acc[3 * p + 2] + w * d[j + 2];
+                            wsum[p] = wsum[p] + w;
+                            if (w > wmax[p])
+                                wmax[p] = w;
+                        }
+                        yield;
+                    }
+                }
+            const out = new Float32Array(3 * n);
+            for (let y = 0; y < H; y++) {
+                for (let p = y * W; p < (y + 1) * W; p++) {
+                    const wc = wmax[p], total = wsum[p] + wc;
+                    for (let c = 0; c < 3; c++)
+                        out[3 * p + c] = total > 0 ? (acc[3 * p + c] + wc * d[3 * p + c]) / total : d[3 * p + c];
+                }
+                yield;
+            }
+            return { width: W, height: H, data: out, provenance: processingProvenance("denoise-non-local-means-buades-2005", [
+                    "Weights exp(-max(d^2 - 2 sigma^2, 0)/h^2); d^2 = mean squared difference over a uniform (2f+1)^2 patch and all channels.",
+                    "Search window (2r+1)^2 with reflect-101 borders; centre weight = largest other weight.",
+                    "One global h and sigma; signal-dependent RAW noise needs a variance-stabilizing transform not applied here."
+                ]) };
+        }
+    });
+}
+/** Unsharp masking: in + amount·(in − G_σ ⊗ in) where the detail exceeds the threshold. */
+export function createUnsharpMaskTask(input, provider) {
+    const r = requireOptions(input, ["image", "radiusSigmaPixels", "amount", "threshold", "maximumWorkingBytes"], "createUnsharpMaskTask input");
+    const image = requireImage(r.image, "image"), sigma = requireFiniteIn(r.radiusSigmaPixels, "radiusSigmaPixels", 0, 256, true);
+    const amount = requireFiniteIn(r.amount, "amount", 0, 100), threshold = optionalFiniteIn(r.threshold, 0, "threshold", 0, 1e30);
+    const W = image.width, H = image.height, kernel = gaussianKernel(sigma);
+    return createProcessingTask({
+        provider, totalTileCount: H + convolveRows(H) + H, estimatedWorkingBytes: W * H * 36, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),
+        body: function* () {
+            yield* finiteRows(image.data, W, H, 3, "image");
+            const d = image.data, blur = new Float32Array(d.length), tmp = new Float32Array(d.length), out = new Float32Array(d.length);
+            yield* convolveSeparable(d, W, H, 3, kernel, tmp, blur);
+            for (let y = 0; y < H; y++) {
+                for (let i = 3 * y * W; i < 3 * (y + 1) * W; i++) {
+                    const detail = d[i] - blur[i];
+                    out[i] = Math.abs(detail) > threshold ? d[i] + amount * detail : d[i];
+                }
+                yield;
+            }
+            return { width: W, height: H, data: out, provenance: processingProvenance("sharpen-unsharp-mask", [
+                    "out = in + amount (in - G_sigma * in) per component where |in - G_sigma * in| > threshold; sigma " + sigma + " px, reflect-101 borders.",
+                    "Output is not clipped; overshoot (halo) at edges is inherent to unsharp masking."
+                ]) };
+        }
+    });
+}
+/** Random-stream id for synthetic grain. Different from every RAW charge, read-noise and fixed-pattern stream. */
+export const SYNTHETIC_FILM_GRAIN_STREAM_ID = 64;
+/**
+ * Adds deterministic synthetic grain:
+ * 1. Unit-variance Gaussian white noise per pixel. It is seeded by the absolute pixel index, so tiling
+ *    never changes it.
+ * 2. Optionally blurred by a Gaussian with σ = grain size, then renormalized to unit variance.
+ * 3. Scaled by a(v) = intensity·2√(v(1 − v)): zero at black and white, `intensity` at mid-grey.
+ * 4. Clipped to [0, 1].
+ * It is labelled synthetic and aesthetic. It is never a film model.
+ */
+export function createSyntheticFilmGrainTask(input, provider) {
+    const r = requireOptions(input, ["image", "seedUint32", "intensity", "grainSizePixels", "mode", "maximumWorkingBytes"], "createSyntheticFilmGrainTask input");
+    const image = requireImage(r.image, "image");
+    if (typeof r.seedUint32 !== "number" || !Number.isInteger(r.seedUint32) || r.seedUint32 < 0 || r.seedUint32 > 0xffffffff) {
+        throw new InvalidConfigurationError("seedUint32 must be a uint32 integer.");
+    }
+    const seed = r.seedUint32, intensity = requireFiniteIn(r.intensity, "intensity", 0, 1), size = optionalFiniteIn(r.grainSizePixels, 0.7, "grainSizePixels", 0, 64);
+    const mode = requireEnum(r.mode, ["monochrome", "per-channel"], "mode", "monochrome"), channels = mode === "monochrome" ? 1 : 3;
+    const W = image.width, H = image.height, n = W * H, kernel = gaussianKernel(size);
+    let k2 = 0;
+    for (const v of kernel)
+        k2 += v * v;
+    return createProcessingTask({
+        provider, totalTileCount: 3 * H + convolveRows(H), estimatedWorkingBytes: n * channels * 12 + n * 12, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),
+        body: function* () {
+            yield* finiteRows(image.data, W, H, 3, "image", 0);
+            const noise = new Float32Array(n * channels), tmp = new Float32Array(n * channels), field = new Float32Array(n * channels);
+            for (let y = 0; y < H; y++) {
+                for (let p = y * W; p < (y + 1) * W; p++) {
+                    const stream = createDeterministicUniformStream(deriveNativeSiteStreamSeedUint32(seed, p, SYNTHETIC_FILM_GRAIN_STREAM_ID));
+                    for (let c = 0; c < channels; c++)
+                        noise[p * channels + c] = sampleDeterministicStandardNormal(stream);
+                }
+                yield;
+            }
+            // Blurred white noise has variance (Σk²)² for a separable kernel; dividing by Σk² restores unit variance.
+            yield* convolveSeparable(noise, W, H, channels, kernel, tmp, field);
+            const d = image.data, out = new Float32Array(3 * n);
+            let clipped = 0;
+            for (let y = 0; y < H; y++) {
+                for (let p = y * W; p < (y + 1) * W; p++)
+                    for (let c = 0; c < 3; c++) {
+                        const v = Math.min(d[3 * p + c], 1), a = intensity * 2 * Math.sqrt(v * (1 - v));
+                        let o = v + a * field[p * channels + (channels === 1 ? 0 : c)] / k2;
+                        if (o < 0) {
+                            o = 0;
+                            clipped++;
+                        }
+                        else if (o > 1) {
+                            o = 1;
+                            clipped++;
+                        }
+                        out[3 * p + c] = o;
+                    }
+                yield;
+            }
+            return { width: W, height: H, data: out, syntheticAesthetic: true, physicalFilmModel: false, clippedSampleCount: clipped,
+                provenance: processingProvenance("synthetic-aesthetic-film-grain", [
+                    "Synthetic aesthetic texture only: not a physical film, silver-halide or dye-cloud model.",
+                    "Deterministic per-pixel Gaussian noise (seed " + seed + ", stream " + SYNTHETIC_FILM_GRAIN_STREAM_ID + "), Gaussian-correlated (sigma " + size + " px), unit variance.",
+                    "Amplitude intensity*2*sqrt(v(1-v)) on encoded values; output clipped to [0, 1]."
+                ]) };
+        }
+    });
+}
+//# sourceMappingURL=denoise.js.map

package/dist/processing/denoise.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"denoise.js","sourceRoot":"","sources":["../../src/processing/denoise.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,gCAAgC,EAAE,gCAAgC,EAAE,iCAAiC,EAAE,MAAM,8BAA8B,CAAC;AACrJ,OAAO,EAAE,YAAY,EAAE,iBAAiB,EAAE,cAAc,EAAE,OAAO,EAAE,MAAM,gBAAgB,CAAC;AAC1F,OAAO,EACL,oBAAoB,EAAE,UAAU,EAAE,mBAAmB,EAAE,gBAAgB,EAAE,iBAAiB,EAAE,oBAAoB,EAAE,WAAW,EAAE,eAAe,EAC9I,YAAY,EAAE,cAAc,EAE7B,MAAM,sBAAsB,CAAC;AAc9B,gFAAgF;AAChF,MAAM,UAAU,yBAAyB,CAAC,KAA2B,EAAE,QAA4B;IACjG,MAAM,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,CAAC,OAAO,EAAE,oBAAoB,EAAE,YAAY,EAAE,qBAAqB,CAAC,EAAE,iCAAiC,CAAC,CAAC;IACzI,MAAM,KAAK,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,OAAO,CAAC,EAAE,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC,kBAAkB,EAAE,oBAAoB,EAAE,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;IAChI,MAAM,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC,UAAU,EAAE,YAAY,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;IAC1I,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IACrD,KAAK,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,IAAI,GAAG,EAAE,EAAE,EAAE;QAAE,KAAK,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,IAAI,GAAG,EAAE,EAAE,EAAE;YAAE,OAAO,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;IAClL,OAAO,oBAAoB,CAAC;QAC1B,QAAQ,EAAE,cAAc,EAAE,CAAC,GAAG,CAAC,EAAE,qBAAqB,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,mBAAmB,EAAE,mBAAmB,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACnI,IAAI,EAAE,QAAQ,CAAC;YACb,KAAK,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;YAChD,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC;YACvF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;oBACtE,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;oBACnC,KAAK,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,IAAI,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;wBACpC,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;wBACnC,KAAK,IAAI,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,IAAI,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;4BACpC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;4BAClF,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;4BAC3H,EAAE,IAAI,CAAC,CAAC;4BAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;4BAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;4BAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;wBACjD,CAAC;oBACH,CAAC;oBACD,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;oBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;oBAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC;gBAC/D,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,UAAU,EAAE,oBAAoB,CAAC,wCAAwC,EAAE;oBAClH,0BAA0B,GAAG,MAAM,GAAG,cAAc,GAAG,GAAG,GAAG,iEAAiE,GAAG,MAAM,GAAG,IAAI;oBAC9I,2IAA2I;iBAAC,CAAC,EAAE,CAAC;QACpJ,CAAC;KACF,CAAC,CAAC;AACL,CAAC;AAeD,mGAAmG;AACnG,MAAM,UAAU,uBAAuB,CAAC,KAAyB,EAAE,QAA4B;IAC7F,MAAM,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,CAAC,OAAO,EAAE,oBAAoB,EAAE,wBAAwB,EAAE,aAAa,EAAE,cAAc,EAAE,qBAAqB,CAAC,EAC7I,+BAA+B,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,OAAO,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,kBAAkB,EAAE,oBAAoB,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IAC7H,MAAM,KAAK,GAAG,gBAAgB,CAAC,CAAC,CAAC,sBAAsB,EAAE,CAAC,EAAE,wBAAwB,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IAC/F,MAAM,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,WAAW,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,iBAAiB,CAAC,CAAC,CAAC,YAAY,EAAE,CAAC,EAAE,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IACtI,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACpF,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC9D,OAAO,oBAAoB,CAAC;QAC1B,QAAQ,EAAE,cAAc,EAAE,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,qBAAqB,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;QACnH,mBAAmB,EAAE,mBAAmB,CAAC,CAAC,CAAC,mBAAmB,CAAC;QAC/D,IAAI,EAAE,QAAQ,CAAC;YACb,KAAK,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;YAChD,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;YAC5G,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC;YACxI,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE;gBAAE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBACrE,IAAI,EAAE,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC;wBAAE,SAAS;oBACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;wBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;4BAC7D,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;wBAC7G,CAAC;wBACD,KAAK,CAAC;oBACR,CAAC;oBACD,KAAK,CAAC,CAAC,iBAAiB,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;oBACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;wBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,OAAO,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;4BAC7G,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;4BAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;4BAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC;4BACzI,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC;4BAAC,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,CAAE;gCAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;wBACxD,CAAC;wBACD,KAAK,CAAC;oBACR,CAAC;gBACH,CAAC;YACD,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBACzC,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAE,EAAE,KAAK,GAAG,IAAI,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC;oBAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;wBAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;gBAC1H,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,UAAU,EAAE,oBAAoB,CAAC,qCAAqC,EAAE;oBAC/G,0HAA0H;oBAC1H,wFAAwF;oBACxF,6GAA6G;iBAAC,CAAC,EAAE,CAAC;QACtH,CAAC;KACF,CAAC,CAAC;AACL,CAAC;AAWD,2FAA2F;AAC3F,MAAM,UAAU,qBAAqB,CAAC,KAAuB,EAAE,QAA4B;IACzF,MAAM,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,CAAC,OAAO,EAAE,mBAAmB,EAAE,QAAQ,EAAE,WAAW,EAAE,qBAAqB,CAAC,EAAE,6BAA6B,CAAC,CAAC;IAC7I,MAAM,KAAK,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,OAAO,CAAC,EAAE,KAAK,GAAG,eAAe,CAAC,CAAC,CAAC,iBAAiB,EAAE,mBAAmB,EAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;IAC9H,MAAM,MAAM,GAAG,eAAe,CAAC,CAAC,CAAC,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,GAAG,CAAC,EAAE,SAAS,GAAG,gBAAgB,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IAC/H,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,MAAM,GAAG,cAAc,CAAC,KAAK,CAAC,CAAC;IACxE,OAAO,oBAAoB,CAAC;QAC1B,QAAQ,EAAE,cAAc,EAAE,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,qBAAqB,EAAE,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,mBAAmB,EAAE,mBAAmB,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACrJ,IAAI,EAAE,QAAQ,CAAC;YACb,KAAK,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;YAChD,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YAC5H,KAAK,CAAC,CAAC,iBAAiB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;YACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBACjD,MAAM,MAAM,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC,CAAE,CAAC;oBAChC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,GAAG,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAC1E,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,UAAU,EAAE,oBAAoB,CAAC,sBAAsB,EAAE;oBAChG,mGAAmG,GAAG,KAAK,GAAG,2BAA2B;oBACzI,kFAAkF;iBAAC,CAAC,EAAE,CAAC;QAC3F,CAAC;KACF,CAAC,CAAC;AACL,CAAC;AAED,kHAAkH;AAClH,MAAM,CAAC,MAAM,8BAA8B,GAAG,EAAW,CAAC;AAuB1D;;;;;;;;GAQG;AACH,MAAM,UAAU,4BAA4B,CAAC,KAA8B,EAAE,QAA4B;IACvG,MAAM,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,CAAC,OAAO,EAAE,YAAY,EAAE,WAAW,EAAE,iBAAiB,EAAE,MAAM,EAAE,qBAAqB,CAAC,EAAE,oCAAoC,CAAC,CAAC;IAC9J,MAAM,KAAK,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC;IAC7C,IAAI,OAAO,CAAC,CAAC,UAAU,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,UAAU,GAAG,CAAC,IAAI,CAAC,CAAC,UAAU,GAAG,UAAU,EAAE,CAAC;QACzH,MAAM,IAAI,yBAAyB,CAAC,sCAAsC,CAAC,CAAC;IAC9E,CAAC;IACD,MAAM,IAAI,GAAG,CAAC,CAAC,UAAU,EAAE,SAAS,GAAG,eAAe,CAAC,CAAC,CAAC,SAAS,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,GAAG,gBAAgB,CAAC,CAAC,CAAC,eAAe,EAAE,GAAG,EAAE,iBAAiB,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;IAClK,MAAM,IAAI,GAAG,WAAW,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,YAAY,EAAE,aAAa,CAAU,EAAE,MAAM,EAAE,YAAY,CAAC,EAAE,QAAQ,GAAG,IAAI,KAAK,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACzI,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,MAAM,GAAG,cAAc,CAAC,IAAI,CAAC,CAAC;IAClF,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,KAAK,MAAM,CAAC,IAAI,MAAM;QAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;IACpC,OAAO,oBAAoB,CAAC;QAC1B,QAAQ,EAAE,cAAc,EAAE,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,EAAE,qBAAqB,EAAE,CAAC,GAAG,QAAQ,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,EAAE,mBAAmB,EAAE,mBAAmB,CAAC,CAAC,CAAC,mBAAmB,CAAC;QACrK,IAAI,EAAE,QAAQ,CAAC;YACb,KAAK,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;YACnD,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC,EAAE,KAAK,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC;YAC3H,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBACzC,MAAM,MAAM,GAAG,gCAAgC,CAAC,gCAAgC,CAAC,IAAI,EAAE,CAAC,EAAE,8BAA8B,CAAC,CAAC,CAAC;oBAC3H,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE;wBAAE,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,iCAAiC,CAAC,MAAM,CAAC,CAAC;gBACzG,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,0GAA0G;YAC1G,KAAK,CAAC,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,EAAE,MAAM,EAAE,GAAG,EAAE,KAAK,CAAC,CAAC;YACpE,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACpD,IAAI,OAAO,GAAG,CAAC,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBACrE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;wBACjF,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,QAAQ,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC;wBACrE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;4BAAC,CAAC,GAAG,CAAC,CAAC;4BAAC,OAAO,EAAE,CAAC;wBAAC,CAAC;6BAAM,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;4BAAC,CAAC,GAAG,CAAC,CAAC;4BAAC,OAAO,EAAE,CAAC;wBAAC,CAAC;wBACtE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;oBACrB,CAAC;gBACD,KAAK,CAAC;YACR,CAAC;YACD,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,kBAAkB,EAAE,IAAI,EAAE,iBAAiB,EAAE,KAAK,EAAE,kBAAkB,EAAE,OAAO;gBACtH,UAAU,EAAE,oBAAoB,CAAC,gCAAgC,EAAE;oBACjE,0FAA0F;oBAC1F,+CAA+C,GAAG,IAAI,GAAG,WAAW,GAAG,8BAA8B,GAAG,gCAAgC,GAAG,IAAI,GAAG,sBAAsB;oBACxK,iFAAiF;iBAAC,CAAC,EAAE,CAAC;QAC5F,CAAC;KACF,CAAC,CAAC;AACL,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Image-domain denoising, sharpening and synthetic film grain on linear or developed RGB:\n * - Bilateral filter (Tomasi & Manduchi, ICCV 1998): Gaussian spatial and Gaussian range weights,\n *   range distance Euclidean over the three channels of the supplied data.\n * - Non-local means (Buades, Coll & Morel, CVPR 2005): weights exp(−max(d² − 2σ², 0)/h²), where d² is\n *   the mean squared difference of (2f+1)² patches over all channels. This is the uniform-patch form\n *   with the known-noise offset from Buades et al., IPOL 2011; σ = 0 gives the 2005 form with a\n *   uniform patch kernel. The centre pixel gets the largest weight in its search window.\n * - Unsharp masking: out = in + amount·(in − G_σ ⊗ in) where |in − G_σ ⊗ in| > threshold.\n * - Synthetic film grain: an aesthetic, deterministic, signal-dependent texture. It is NOT a physical\n *   film model (no silver-halide grain statistics, dye clouds or film characteristic curve).\n * Noise in linear RAW-derived data depends on the signal level, but these filters use one global set of\n * parameters.\n * @see docs/PROCESSING.md for equations, sources and limits.\n */\n\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { createDeterministicUniformStream, deriveNativeSiteStreamSeedUint32, sampleDeterministicStandardNormal } from \"../sensor/stochastic-seed.js\";\nimport { convolveRows, convolveSeparable, gaussianKernel, reflect } from \"./image-ops.js\";\nimport {\n  createProcessingTask, finiteRows, maximumWorkingBytes, optionalFiniteIn, optionalIntegerIn, processingProvenance, requireEnum, requireFiniteIn,\n  requireImage, requireOptions,\n  type ProcessingBody, type ProcessingProvenance, type ProcessingProvider, type ProcessingRgbImage, type ProcessingTask\n} from \"./processing-task.js\";\n\nexport interface ProcessedRgbImage extends ProcessingRgbImage {\n  provenance: ProcessingProvenance;\n}\n\nexport interface BilateralFilterInput {\n  image: ProcessingRgbImage;\n  spatialSigmaPixels: number;\n  /** Range σ in the units of the image values. */\n  rangeSigma: number;\n  maximumWorkingBytes?: number;\n}\n\n/** Brute-force bilateral filter (window radius ⌈3σ_s⌉, reflect-101 borders). */\nexport function createBilateralFilterTask(input: BilateralFilterInput, provider: ProcessingProvider): ProcessingTask<ProcessedRgbImage> {\n  const r = requireOptions(input, [\"image\", \"spatialSigmaPixels\", \"rangeSigma\", \"maximumWorkingBytes\"], \"createBilateralFilterTask input\");\n  const image = requireImage(r.image, \"image\"), sigmaS = requireFiniteIn(r.spatialSigmaPixels, \"spatialSigmaPixels\", 0, 64, true);\n  const sigmaR = requireFiniteIn(r.rangeSigma, \"rangeSigma\", 0, 1e30, true), W = image.width, H = image.height, rad = Math.ceil(3 * sigmaS);\n  const spatial = new Float64Array((2 * rad + 1) ** 2);\n  for (let dy = -rad; dy <= rad; dy++) for (let dx = -rad; dx <= rad; dx++) spatial[(dy + rad) * (2 * rad + 1) + dx + rad] = Math.exp(-(dx * dx + dy * dy) / (2 * sigmaS * sigmaS));\n  return createProcessingTask({\n    provider, totalTileCount: 2 * H, estimatedWorkingBytes: W * H * 12, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),\n    body: function* (): ProcessingBody<ProcessedRgbImage> {\n      yield* finiteRows(image.data, W, H, 3, \"image\");\n      const d = image.data, out = new Float32Array(3 * W * H), k = 1 / (2 * sigmaR * sigmaR);\n      for (let y = 0; y < H; y++) {\n        for (let x = 0; x < W; x++) {\n          const i = 3 * (y * W + x), r0 = d[i]!, g0 = d[i + 1]!, b0 = d[i + 2]!;\n          let ws = 0, sr = 0, sg = 0, sb = 0;\n          for (let dy = -rad; dy <= rad; dy++) {\n            const row = reflect(y + dy, H) * W;\n            for (let dx = -rad; dx <= rad; dx++) {\n              const j = 3 * (row + reflect(x + dx, W)), R = d[j]!, G = d[j + 1]!, B = d[j + 2]!;\n              const w = spatial[(dy + rad) * (2 * rad + 1) + dx + rad]! * Math.exp(-((R - r0) ** 2 + (G - g0) ** 2 + (B - b0) ** 2) * k);\n              ws += w; sr += w * R; sg += w * G; sb += w * B;\n            }\n          }\n          out[i] = sr / ws; out[i + 1] = sg / ws; out[i + 2] = sb / ws;\n        }\n        yield;\n      }\n      return { width: W, height: H, data: out, provenance: processingProvenance(\"denoise-bilateral-tomasi-manduchi-1998\", [\n        \"Gaussian spatial (sigma \" + sigmaS + \" px, radius \" + rad + \") times Gaussian range weight on Euclidean RGB distance (sigma \" + sigmaR + \").\",\n        \"Range distance in the supplied data's units, not CIELAB; reflect-101 borders; one global parameter set (noise is not modelled per pixel).\"]) };\n    }\n  });\n}\n\nexport interface NonLocalMeansInput {\n  image: ProcessingRgbImage;\n  /** Filtering parameter h (image units). */\n  filteringParameter: number;\n  /** Known noise standard deviation σ (image units); default 0. */\n  noiseStandardDeviation?: number;\n  /** Patch radius f (patch (2f+1)²); default 1. */\n  patchRadius?: number;\n  /** Search-window radius r ((2r+1)² candidates); default 5. */\n  searchRadius?: number;\n  maximumWorkingBytes?: number;\n}\n\n/** Non-local means (Buades, Coll & Morel 2005), offset-major with box-filtered patch distances. */\nexport function createNonLocalMeansTask(input: NonLocalMeansInput, provider: ProcessingProvider): ProcessingTask<ProcessedRgbImage> {\n  const r = requireOptions(input, [\"image\", \"filteringParameter\", \"noiseStandardDeviation\", \"patchRadius\", \"searchRadius\", \"maximumWorkingBytes\"],\n    \"createNonLocalMeansTask input\");\n  const image = requireImage(r.image, \"image\"), h = requireFiniteIn(r.filteringParameter, \"filteringParameter\", 0, 1e30, true);\n  const sigma = optionalFiniteIn(r.noiseStandardDeviation, 0, \"noiseStandardDeviation\", 0, 1e30);\n  const f = optionalIntegerIn(r.patchRadius, 1, \"patchRadius\", 0, 16), sr = optionalIntegerIn(r.searchRadius, 5, \"searchRadius\", 1, 32);\n  const W = image.width, H = image.height, n = W * H, offsets = (2 * sr + 1) ** 2 - 1;\n  const box = new Float64Array(2 * f + 1).fill(1 / (2 * f + 1));\n  return createProcessingTask({\n    provider, totalTileCount: H + offsets * (2 * H + convolveRows(H)) + H, estimatedWorkingBytes: n * (12 + 4 * 5 + 12),\n    maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),\n    body: function* (): ProcessingBody<ProcessedRgbImage> {\n      yield* finiteRows(image.data, W, H, 3, \"image\");\n      const d = image.data, acc = new Float32Array(3 * n), wsum = new Float32Array(n), wmax = new Float32Array(n);\n      const diff = new Float32Array(n), tmp = new Float32Array(n), dist = new Float32Array(n), inv = 1 / (h * h), offset2 = 2 * sigma * sigma;\n      for (let oy = -sr; oy <= sr; oy++) for (let ox = -sr; ox <= sr; ox++) {\n        if (ox === 0 && oy === 0) continue;\n        for (let y = 0; y < H; y++) {\n          const yy = reflect(y + oy, H) * W;\n          for (let x = 0; x < W; x++) {\n            const i = 3 * (y * W + x), j = 3 * (yy + reflect(x + ox, W));\n            diff[y * W + x] = ((d[i]! - d[j]!) ** 2 + (d[i + 1]! - d[j + 1]!) ** 2 + (d[i + 2]! - d[j + 2]!) ** 2) / 3;\n          }\n          yield;\n        }\n        yield* convolveSeparable(diff, W, H, 1, box, tmp, dist);\n        for (let y = 0; y < H; y++) {\n          const yy = reflect(y + oy, H) * W;\n          for (let x = 0; x < W; x++) {\n            const p = y * W + x, j = 3 * (yy + reflect(x + ox, W)), w = Math.exp(-Math.max(dist[p]! - offset2, 0) * inv);\n            acc[3 * p] = acc[3 * p]! + w * d[j]!; acc[3 * p + 1] = acc[3 * p + 1]! + w * d[j + 1]!; acc[3 * p + 2] = acc[3 * p + 2]! + w * d[j + 2]!;\n            wsum[p] = wsum[p]! + w; if (w > wmax[p]!) wmax[p] = w;\n          }\n          yield;\n        }\n      }\n      const out = new Float32Array(3 * n);\n      for (let y = 0; y < H; y++) {\n        for (let p = y * W; p < (y + 1) * W; p++) {\n          const wc = wmax[p]!, total = wsum[p]! + wc;\n          for (let c = 0; c < 3; c++) out[3 * p + c] = total > 0 ? (acc[3 * p + c]! + wc * d[3 * p + c]!) / total : d[3 * p + c]!;\n        }\n        yield;\n      }\n      return { width: W, height: H, data: out, provenance: processingProvenance(\"denoise-non-local-means-buades-2005\", [\n        \"Weights exp(-max(d^2 - 2 sigma^2, 0)/h^2); d^2 = mean squared difference over a uniform (2f+1)^2 patch and all channels.\",\n        \"Search window (2r+1)^2 with reflect-101 borders; centre weight = largest other weight.\",\n        \"One global h and sigma; signal-dependent RAW noise needs a variance-stabilizing transform not applied here.\"]) };\n    }\n  });\n}\n\nexport interface UnsharpMaskInput {\n  image: ProcessingRgbImage;\n  radiusSigmaPixels: number;\n  amount: number;\n  /** Detail magnitude below which a component is left unchanged; default 0. */\n  threshold?: number;\n  maximumWorkingBytes?: number;\n}\n\n/** Unsharp masking: in + amount·(in − G_σ ⊗ in) where the detail exceeds the threshold. */\nexport function createUnsharpMaskTask(input: UnsharpMaskInput, provider: ProcessingProvider): ProcessingTask<ProcessedRgbImage> {\n  const r = requireOptions(input, [\"image\", \"radiusSigmaPixels\", \"amount\", \"threshold\", \"maximumWorkingBytes\"], \"createUnsharpMaskTask input\");\n  const image = requireImage(r.image, \"image\"), sigma = requireFiniteIn(r.radiusSigmaPixels, \"radiusSigmaPixels\", 0, 256, true);\n  const amount = requireFiniteIn(r.amount, \"amount\", 0, 100), threshold = optionalFiniteIn(r.threshold, 0, \"threshold\", 0, 1e30);\n  const W = image.width, H = image.height, kernel = gaussianKernel(sigma);\n  return createProcessingTask({\n    provider, totalTileCount: H + convolveRows(H) + H, estimatedWorkingBytes: W * H * 36, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),\n    body: function* (): ProcessingBody<ProcessedRgbImage> {\n      yield* finiteRows(image.data, W, H, 3, \"image\");\n      const d = image.data, blur = new Float32Array(d.length), tmp = new Float32Array(d.length), out = new Float32Array(d.length);\n      yield* convolveSeparable(d, W, H, 3, kernel, tmp, blur);\n      for (let y = 0; y < H; y++) {\n        for (let i = 3 * y * W; i < 3 * (y + 1) * W; i++) {\n          const detail = d[i]! - blur[i]!;\n          out[i] = Math.abs(detail) > threshold ? d[i]! + amount * detail : d[i]!;\n        }\n        yield;\n      }\n      return { width: W, height: H, data: out, provenance: processingProvenance(\"sharpen-unsharp-mask\", [\n        \"out = in + amount (in - G_sigma * in) per component where |in - G_sigma * in| > threshold; sigma \" + sigma + \" px, reflect-101 borders.\",\n        \"Output is not clipped; overshoot (halo) at edges is inherent to unsharp masking.\"]) };\n    }\n  });\n}\n\n/** Random-stream id for synthetic grain. Different from every RAW charge, read-noise and fixed-pattern stream. */\nexport const SYNTHETIC_FILM_GRAIN_STREAM_ID = 64 as const;\n\nexport interface SyntheticFilmGrainInput {\n  /** Display-referred (encoded) RGB with values in [0, 1]. */\n  image: ProcessingRgbImage;\n  /** uint32 seed; equal seeds reproduce the grain exactly. */\n  seedUint32: number;\n  /** Grain standard deviation at mid-grey (value 0.5), in encoded units. */\n  intensity: number;\n  /** Spatial correlation σ of the grain in pixels; 0 = per-pixel white grain. Default 0.7. */\n  grainSizePixels?: number;\n  /** `monochrome` (default) adds one grain field to all channels; `per-channel` uses three. */\n  mode?: \"monochrome\" | \"per-channel\";\n  maximumWorkingBytes?: number;\n}\n\nexport interface SyntheticFilmGrainOutput extends ProcessedRgbImage {\n  /** Always true: the grain is an aesthetic texture, not a physical film model. */\n  syntheticAesthetic: true;\n  physicalFilmModel: false;\n  clippedSampleCount: number;\n}\n\n/**\n * Adds deterministic synthetic grain:\n * 1. Unit-variance Gaussian white noise per pixel. It is seeded by the absolute pixel index, so tiling\n *    never changes it.\n * 2. Optionally blurred by a Gaussian with σ = grain size, then renormalized to unit variance.\n * 3. Scaled by a(v) = intensity·2√(v(1 − v)): zero at black and white, `intensity` at mid-grey.\n * 4. Clipped to [0, 1].\n * It is labelled synthetic and aesthetic. It is never a film model.\n */\nexport function createSyntheticFilmGrainTask(input: SyntheticFilmGrainInput, provider: ProcessingProvider): ProcessingTask<SyntheticFilmGrainOutput> {\n  const r = requireOptions(input, [\"image\", \"seedUint32\", \"intensity\", \"grainSizePixels\", \"mode\", \"maximumWorkingBytes\"], \"createSyntheticFilmGrainTask input\");\n  const image = requireImage(r.image, \"image\");\n  if (typeof r.seedUint32 !== \"number\" || !Number.isInteger(r.seedUint32) || r.seedUint32 < 0 || r.seedUint32 > 0xffffffff) {\n    throw new InvalidConfigurationError(\"seedUint32 must be a uint32 integer.\");\n  }\n  const seed = r.seedUint32, intensity = requireFiniteIn(r.intensity, \"intensity\", 0, 1), size = optionalFiniteIn(r.grainSizePixels, 0.7, \"grainSizePixels\", 0, 64);\n  const mode = requireEnum(r.mode, [\"monochrome\", \"per-channel\"] as const, \"mode\", \"monochrome\"), channels = mode === \"monochrome\" ? 1 : 3;\n  const W = image.width, H = image.height, n = W * H, kernel = gaussianKernel(size);\n  let k2 = 0;\n  for (const v of kernel) k2 += v * v;\n  return createProcessingTask({\n    provider, totalTileCount: 3 * H + convolveRows(H), estimatedWorkingBytes: n * channels * 12 + n * 12, maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),\n    body: function* (): ProcessingBody<SyntheticFilmGrainOutput> {\n      yield* finiteRows(image.data, W, H, 3, \"image\", 0);\n      const noise = new Float32Array(n * channels), tmp = new Float32Array(n * channels), field = new Float32Array(n * channels);\n      for (let y = 0; y < H; y++) {\n        for (let p = y * W; p < (y + 1) * W; p++) {\n          const stream = createDeterministicUniformStream(deriveNativeSiteStreamSeedUint32(seed, p, SYNTHETIC_FILM_GRAIN_STREAM_ID));\n          for (let c = 0; c < channels; c++) noise[p * channels + c] = sampleDeterministicStandardNormal(stream);\n        }\n        yield;\n      }\n      // Blurred white noise has variance (Σk²)² for a separable kernel; dividing by Σk² restores unit variance.\n      yield* convolveSeparable(noise, W, H, channels, kernel, tmp, field);\n      const d = image.data, out = new Float32Array(3 * n);\n      let clipped = 0;\n      for (let y = 0; y < H; y++) {\n        for (let p = y * W; p < (y + 1) * W; p++) for (let c = 0; c < 3; c++) {\n          const v = Math.min(d[3 * p + c]!, 1), a = intensity * 2 * Math.sqrt(v * (1 - v));\n          let o = v + a * field[p * channels + (channels === 1 ? 0 : c)]! / k2;\n          if (o < 0) { o = 0; clipped++; } else if (o > 1) { o = 1; clipped++; }\n          out[3 * p + c] = o;\n        }\n        yield;\n      }\n      return { width: W, height: H, data: out, syntheticAesthetic: true, physicalFilmModel: false, clippedSampleCount: clipped,\n        provenance: processingProvenance(\"synthetic-aesthetic-film-grain\", [\n          \"Synthetic aesthetic texture only: not a physical film, silver-halide or dye-cloud model.\",\n          \"Deterministic per-pixel Gaussian noise (seed \" + seed + \", stream \" + SYNTHETIC_FILM_GRAIN_STREAM_ID + \"), Gaussian-correlated (sigma \" + size + \" px), unit variance.\",\n          \"Amplitude intensity*2*sqrt(v(1-v)) on encoded values; output clipped to [0, 1].\"]) };\n    }\n  });\n}\n"]}

package/dist/processing/focus-stacking.d.ts

@@ -0,0 +1,17 @@
+import { type ProcessingProvenance, type ProcessingProvider, type ProcessingRgbImage, type ProcessingTask } from "./processing-task.js";
+export interface FocusStackInput {
+    /** Registered frames of one scene, ordered by focus distance (any consistent order). */
+    frames: readonly ProcessingRgbImage[];
+    /** σ (pixels) of the Gaussian window averaging the squared Laplacian; default 2. */
+    sharpnessSigmaPixels?: number;
+    /** Pyramid levels for seam-free blending; default down to a 1-pixel side (at most 10). */
+    pyramidLevels?: number;
+    maximumWorkingBytes?: number;
+}
+export interface FocusStackOutput extends ProcessingRgbImage {
+    /** Index of the sharpest frame per pixel (the blending selection). */
+    selectedFrameIndex: Uint8Array;
+    provenance: ProcessingProvenance;
+}
+/** Laplacian-energy focus stacking with Laplacian-pyramid blending. */
+export declare function createFocusStackTask(input: FocusStackInput, provider: ProcessingProvider): ProcessingTask<FocusStackOutput>;

package/dist/processing/focus-stacking.js

@@ -0,0 +1,80 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * Focus stacking: combines aligned (registered) frames focused at different distances into one sharp
+ * image.
+ *
+ * - Sharpness for each frame and pixel is the locally averaged Laplacian energy S_k = G_σ ⊗ (∇²Y_k)²
+ *   of the Rec.709-weighted luma Y. This is a Laplacian-based focus measure (see Nayar & Nakagawa 1994,
+ *   and the survey by Pertuz, Puig & García 2013).
+ * - Each pixel picks the sharpest frame; a tie picks the lowest index.
+ * - The selection masks blend the frames with Laplacian pyramids (Burt & Adelson 1983; pyramid image
+ *   fusion as in Burt & Kolczynski 1993). This hides the seams between selections at every scale.
+ *
+ * Not done: frame alignment, magnification (focus-breathing) correction and depth estimation. The
+ * index map only says which frame was picked; it is not calibrated depth.
+ * @see docs/PROCESSING.md for equations, sources and limits.
+ */
+import { autoLevels, collapse, collapseRows, convolveRows, convolveSeparable, gaussianKernel, gaussianPyramid, gaussianPyramidRows, laplacianPyramid, laplacianPyramidRows, luminanceWeights, pyramidBytes, reflect } from "./image-ops.js";
+import { createProcessingTask, finiteRows, maximumWorkingBytes, optionalFiniteIn, optionalIntegerIn, processingProvenance, requireImageList, requireOptions } from "./processing-task.js";
+/** Laplacian-energy focus stacking with Laplacian-pyramid blending. */
+export function createFocusStackTask(input, provider) {
+    const r = requireOptions(input, ["frames", "sharpnessSigmaPixels", "pyramidLevels", "maximumWorkingBytes"], "createFocusStackTask input");
+    const frames = requireImageList(r.frames, "frames", 1), W = frames[0].width, H = frames[0].height, n = W * H, K = frames.length;
+    const sigma = optionalFiniteIn(r.sharpnessSigmaPixels, 2, "sharpnessSigmaPixels", 0, 1000);
+    const levels = optionalIntegerIn(r.pyramidLevels, autoLevels(W, H, 10), "pyramidLevels", 1, autoLevels(W, H, 32));
+    const kernel = gaussianKernel(sigma);
+    const total = K * H + K * (H + convolveRows(H)) + K * (laplacianPyramidRows(W, H, levels) + gaussianPyramidRows(W, H, levels) + H) + collapseRows(W, H, levels);
+    return createProcessingTask({
+        provider, totalTileCount: total, estimatedWorkingBytes: n * 4 * 4 + n + 2 * pyramidBytes(W, H, 3, levels) + 2 * pyramidBytes(W, H, 1, levels),
+        maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),
+        body: function* () {
+            for (const [k, f] of frames.entries())
+                yield* finiteRows(f.data, W, H, 3, "frames[" + k + "]");
+            const [lr, lg, lb] = luminanceWeights("linear-srgb-d65");
+            const energy = new Float32Array(n), tmp = new Float32Array(n), smooth = new Float32Array(n), best = new Float32Array(n).fill(-1);
+            const index = new Uint8Array(n);
+            for (const [k, f] of frames.entries()) {
+                const d = f.data;
+                const y709 = (x, y) => { const i = 3 * (reflect(y, H) * W + reflect(x, W)); return lr * d[i] + lg * d[i + 1] + lb * d[i + 2]; };
+                for (let y = 0; y < H; y++) {
+                    for (let x = 0; x < W; x++) {
+                        const lap = y709(x - 1, y) + y709(x + 1, y) + y709(x, y - 1) + y709(x, y + 1) - 4 * y709(x, y);
+                        energy[y * W + x] = lap * lap;
+                    }
+                    yield;
+                }
+                yield* convolveSeparable(energy, W, H, 1, kernel, tmp, smooth);
+                for (let i = 0; i < n; i++)
+                    if (smooth[i] > best[i]) {
+                        best[i] = smooth[i];
+                        index[i] = k;
+                    }
+            }
+            let acc = null, sizes = [];
+            const mask = new Float32Array(n);
+            for (const [k, f] of frames.entries()) {
+                for (let y = 0; y < H; y++) {
+                    for (let i = y * W; i < (y + 1) * W; i++)
+                        mask[i] = index[i] === k ? 1 : 0;
+                    yield;
+                }
+                const lap = yield* laplacianPyramid(f.data, W, H, 3, levels), gw = yield* gaussianPyramid(mask, W, H, 1, levels);
+                sizes = lap.sizes;
+                acc ??= lap.sizes.map((s) => new Float32Array(s.width * s.height * 3));
+                for (let l = 0; l < levels; l++) {
+                    const a = acc[l], L = lap.levels[l], G = gw.levels[l];
+                    for (let i = 0; i < G.length; i++)
+                        for (let c = 0; c < 3; c++)
+                            a[3 * i + c] = a[3 * i + c] + G[i] * L[3 * i + c];
+                }
+            }
+            const data = yield* collapse({ levels: acc, sizes, channels: 3 });
+            return { width: W, height: H, data, selectedFrameIndex: index, provenance: processingProvenance("focus-stack-laplacian-energy-pyramid", [
+                    "Sharpness: Gaussian-averaged squared 3x3 Laplacian of Rec.709-weighted luma; argmax per pixel, ties to the lowest index.",
+                    "One-hot masks blended with Burt-Adelson Laplacian pyramids (reflect-101 borders); values are not clipped.",
+                    "Frames must be registered; no alignment, focus-breathing correction or calibrated depth is produced."
+                ]) };
+        }
+    });
+}
+//# sourceMappingURL=focus-stacking.js.map

package/dist/processing/focus-stacking.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"focus-stacking.js","sourceRoot":"","sources":["../../src/processing/focus-stacking.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EACL,UAAU,EAAE,QAAQ,EAAE,YAAY,EAAE,YAAY,EAAE,iBAAiB,EAAE,cAAc,EAAE,eAAe,EAAE,mBAAmB,EAAE,gBAAgB,EAC3I,oBAAoB,EAAE,gBAAgB,EAAE,YAAY,EAAE,OAAO,EAC9D,MAAM,gBAAgB,CAAC;AACxB,OAAO,EACL,oBAAoB,EAAE,UAAU,EAAE,mBAAmB,EAAE,gBAAgB,EAAE,iBAAiB,EAAE,oBAAoB,EAAE,gBAAgB,EAAE,cAAc,EAEnJ,MAAM,sBAAsB,CAAC;AAkB9B,uEAAuE;AACvE,MAAM,UAAU,oBAAoB,CAAC,KAAsB,EAAE,QAA4B;IACvF,MAAM,CAAC,GAAG,cAAc,CAAC,KAAK,EAAE,CAAC,QAAQ,EAAE,sBAAsB,EAAE,eAAe,EAAE,qBAAqB,CAAC,EAAE,4BAA4B,CAAC,CAAC;IAC1I,MAAM,MAAM,GAAG,gBAAgB,CAAC,CAAC,CAAC,MAAM,EAAE,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC,KAAK,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IAClI,MAAM,KAAK,GAAG,gBAAgB,CAAC,CAAC,CAAC,oBAAoB,EAAE,CAAC,EAAE,sBAAsB,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IAC3F,MAAM,MAAM,GAAG,iBAAiB,CAAC,CAAC,CAAC,aAAa,EAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,eAAe,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IAClH,MAAM,MAAM,GAAG,cAAc,CAAC,KAAK,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,mBAAmB,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;IAChK,OAAO,oBAAoB,CAAC;QAC1B,QAAQ,EAAE,cAAc,EAAE,KAAK,EAAE,qBAAqB,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,GAAG,YAAY,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC;QAC7I,mBAAmB,EAAE,mBAAmB,CAAC,CAAC,CAAC,mBAAmB,CAAC;QAC/D,IAAI,EAAE,QAAQ,CAAC;YACb,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE;gBAAE,KAAK,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;YAC/F,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,gBAAgB,CAAC,iBAAiB,CAAC,CAAC;YACzD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YACjI,MAAM,KAAK,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;YAChC,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,EAAE,CAAC;gBACtC,MAAM,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;gBACjB,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC,CAAC,CAAC;gBAC3J,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;wBAC/F,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;oBAChC,CAAC;oBACD,KAAK,CAAC;gBACR,CAAC;gBACD,KAAK,CAAC,CAAC,iBAAiB,CAAC,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC;gBAC/D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,IAAI,MAAM,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC,CAAE,EAAE,CAAC;wBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC;wBAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBAAC,CAAC;YAChG,CAAC;YACD,IAAI,GAAG,GAA0B,IAAI,EAAE,KAAK,GAAwC,EAAE,CAAC;YACvF,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;YACjC,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,MAAM,CAAC,OAAO,EAAE,EAAE,CAAC;gBACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;wBAAE,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3E,KAAK,CAAC;gBACR,CAAC;gBACD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,eAAe,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;gBACjH,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;gBAClB,GAAG,KAAK,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,YAAY,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;gBACvE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBAChC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAE,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC,CAAC,CAAE,CAAC;oBACzD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;wBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;4BAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;gBACtH,CAAC;YACH,CAAC;YACD,MAAM,IAAI,GAAG,KAAK,CAAC,CAAC,QAAQ,CAAC,EAAE,MAAM,EAAE,GAAI,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC;YACnE,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,kBAAkB,EAAE,KAAK,EAAE,UAAU,EAAE,oBAAoB,CAAC,sCAAsC,EAAE;oBACtI,0HAA0H;oBAC1H,2GAA2G;oBAC3G,sGAAsG;iBAAC,CAAC,EAAE,CAAC;QAC/G,CAAC;KACF,CAAC,CAAC;AACL,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Focus stacking: combines aligned (registered) frames focused at different distances into one sharp\n * image.\n *\n * - Sharpness for each frame and pixel is the locally averaged Laplacian energy S_k = G_σ ⊗ (∇²Y_k)²\n *   of the Rec.709-weighted luma Y. This is a Laplacian-based focus measure (see Nayar & Nakagawa 1994,\n *   and the survey by Pertuz, Puig & García 2013).\n * - Each pixel picks the sharpest frame; a tie picks the lowest index.\n * - The selection masks blend the frames with Laplacian pyramids (Burt & Adelson 1983; pyramid image\n *   fusion as in Burt & Kolczynski 1993). This hides the seams between selections at every scale.\n *\n * Not done: frame alignment, magnification (focus-breathing) correction and depth estimation. The\n * index map only says which frame was picked; it is not calibrated depth.\n * @see docs/PROCESSING.md for equations, sources and limits.\n */\n\nimport {\n  autoLevels, collapse, collapseRows, convolveRows, convolveSeparable, gaussianKernel, gaussianPyramid, gaussianPyramidRows, laplacianPyramid,\n  laplacianPyramidRows, luminanceWeights, pyramidBytes, reflect\n} from \"./image-ops.js\";\nimport {\n  createProcessingTask, finiteRows, maximumWorkingBytes, optionalFiniteIn, optionalIntegerIn, processingProvenance, requireImageList, requireOptions,\n  type ProcessingBody, type ProcessingProvenance, type ProcessingProvider, type ProcessingRgbImage, type ProcessingTask\n} from \"./processing-task.js\";\n\nexport interface FocusStackInput {\n  /** Registered frames of one scene, ordered by focus distance (any consistent order). */\n  frames: readonly ProcessingRgbImage[];\n  /** σ (pixels) of the Gaussian window averaging the squared Laplacian; default 2. */\n  sharpnessSigmaPixels?: number;\n  /** Pyramid levels for seam-free blending; default down to a 1-pixel side (at most 10). */\n  pyramidLevels?: number;\n  maximumWorkingBytes?: number;\n}\n\nexport interface FocusStackOutput extends ProcessingRgbImage {\n  /** Index of the sharpest frame per pixel (the blending selection). */\n  selectedFrameIndex: Uint8Array;\n  provenance: ProcessingProvenance;\n}\n\n/** Laplacian-energy focus stacking with Laplacian-pyramid blending. */\nexport function createFocusStackTask(input: FocusStackInput, provider: ProcessingProvider): ProcessingTask<FocusStackOutput> {\n  const r = requireOptions(input, [\"frames\", \"sharpnessSigmaPixels\", \"pyramidLevels\", \"maximumWorkingBytes\"], \"createFocusStackTask input\");\n  const frames = requireImageList(r.frames, \"frames\", 1), W = frames[0]!.width, H = frames[0]!.height, n = W * H, K = frames.length;\n  const sigma = optionalFiniteIn(r.sharpnessSigmaPixels, 2, \"sharpnessSigmaPixels\", 0, 1000);\n  const levels = optionalIntegerIn(r.pyramidLevels, autoLevels(W, H, 10), \"pyramidLevels\", 1, autoLevels(W, H, 32));\n  const kernel = gaussianKernel(sigma);\n  const total = K * H + K * (H + convolveRows(H)) + K * (laplacianPyramidRows(W, H, levels) + gaussianPyramidRows(W, H, levels) + H) + collapseRows(W, H, levels);\n  return createProcessingTask({\n    provider, totalTileCount: total, estimatedWorkingBytes: n * 4 * 4 + n + 2 * pyramidBytes(W, H, 3, levels) + 2 * pyramidBytes(W, H, 1, levels),\n    maximumWorkingBytes: maximumWorkingBytes(r.maximumWorkingBytes),\n    body: function* (): ProcessingBody<FocusStackOutput> {\n      for (const [k, f] of frames.entries()) yield* finiteRows(f.data, W, H, 3, \"frames[\" + k + \"]\");\n      const [lr, lg, lb] = luminanceWeights(\"linear-srgb-d65\");\n      const energy = new Float32Array(n), tmp = new Float32Array(n), smooth = new Float32Array(n), best = new Float32Array(n).fill(-1);\n      const index = new Uint8Array(n);\n      for (const [k, f] of frames.entries()) {\n        const d = f.data;\n        const y709 = (x: number, y: number): number => { const i = 3 * (reflect(y, H) * W + reflect(x, W)); return lr * d[i]! + lg * d[i + 1]! + lb * d[i + 2]!; };\n        for (let y = 0; y < H; y++) {\n          for (let x = 0; x < W; x++) {\n            const lap = y709(x - 1, y) + y709(x + 1, y) + y709(x, y - 1) + y709(x, y + 1) - 4 * y709(x, y);\n            energy[y * W + x] = lap * lap;\n          }\n          yield;\n        }\n        yield* convolveSeparable(energy, W, H, 1, kernel, tmp, smooth);\n        for (let i = 0; i < n; i++) if (smooth[i]! > best[i]!) { best[i] = smooth[i]!; index[i] = k; }\n      }\n      let acc: Float32Array[] | null = null, sizes: { width: number; height: number }[] = [];\n      const mask = new Float32Array(n);\n      for (const [k, f] of frames.entries()) {\n        for (let y = 0; y < H; y++) {\n          for (let i = y * W; i < (y + 1) * W; i++) mask[i] = index[i] === k ? 1 : 0;\n          yield;\n        }\n        const lap = yield* laplacianPyramid(f.data, W, H, 3, levels), gw = yield* gaussianPyramid(mask, W, H, 1, levels);\n        sizes = lap.sizes;\n        acc ??= lap.sizes.map((s) => new Float32Array(s.width * s.height * 3));\n        for (let l = 0; l < levels; l++) {\n          const a = acc[l]!, L = lap.levels[l]!, G = gw.levels[l]!;\n          for (let i = 0; i < G.length; i++) for (let c = 0; c < 3; c++) a[3 * i + c] = a[3 * i + c]! + G[i]! * L[3 * i + c]!;\n        }\n      }\n      const data = yield* collapse({ levels: acc!, sizes, channels: 3 });\n      return { width: W, height: H, data, selectedFrameIndex: index, provenance: processingProvenance(\"focus-stack-laplacian-energy-pyramid\", [\n        \"Sharpness: Gaussian-averaged squared 3x3 Laplacian of Rec.709-weighted luma; argmax per pixel, ties to the lowest index.\",\n        \"One-hot masks blended with Burt-Adelson Laplacian pyramids (reflect-101 borders); values are not clipped.\",\n        \"Frames must be registered; no alignment, focus-breathing correction or calibrated depth is produced.\"]) };\n    }\n  });\n}\n"]}

package/dist/processing/hdr.d.ts

@@ -0,0 +1,54 @@
+import { type ProcessingLinearColorSpace, type ProcessingProvenance, type ProcessingProvider, type ProcessingRgbImage, type ProcessingTask } from "./processing-task.js";
+export interface HdrBracketFrame {
+    exposureTimeSeconds: number;
+    /** Linear RAW codes, one per sample. */
+    codes: Uint16Array | Float32Array;
+    blackLevelCode: number;
+    /** Codes at or above this level are saturated and excluded. */
+    whiteLevelCode: number;
+    /** Optional per-sample flags. Non-zero excludes the sample (for example physical saturation before the ADC). */
+    saturationFlags?: Uint8Array;
+    /** System conversion gain K, electrons per code (includes the frame's ISO gain). */
+    conversionGainElectronsPerCode: number;
+    /** RMS read noise in electrons for this frame's readout regime. */
+    readNoiseElectrons: number;
+}
+export interface HdrBracketMergeInput {
+    /** Samples per frame (e.g. width·height for a mosaic, width·height·3 for RGB). */
+    sampleCount: number;
+    frames: readonly HdrBracketFrame[];
+    /** Variance reweighting iterations after the shot-noise-limited start; default 2. */
+    iterations?: number;
+    maximumWorkingBytes?: number;
+}
+export interface HdrBracketMergeOutput {
+    /** Scene-linear photo-electron rate per sample, e⁻/s (relative scene radiance; not absolute). */
+    radianceElectronsPerSecond: Float32Array;
+    /** Model standard deviation of each estimate, e⁻/s. */
+    standardDeviationElectronsPerSecond: Float32Array;
+    /** Frames that contributed; 0 marks a sample saturated in every frame (value is a lower bound). */
+    contributingFrameCount: Uint8Array;
+    saturatedInAllFramesCount: number;
+    provenance: ProcessingProvenance;
+}
+/** Noise-optimal merge of linear RAW exposure brackets into a scene-linear photo-electron rate. */
+export declare function createHdrBracketMergeTask(input: HdrBracketMergeInput, provider: ProcessingProvider): ProcessingTask<HdrBracketMergeOutput>;
+/** The only output transfer in 2.0: unclipped float32 linear values (PQ and HLG are not shipped). */
+export type HdrTransfer = "linear-float32";
+export interface HdrImageEncodingInput {
+    image: ProcessingRgbImage;
+    inputColorSpace: ProcessingLinearColorSpace;
+    transfer: HdrTransfer;
+    maximumWorkingBytes?: number;
+}
+export interface HdrImageEncodingOutput {
+    width: number;
+    height: number;
+    colorSpace: "rec2020";
+    transfer: HdrTransfer;
+    /** Interleaved linear BT.2020 RGB (relative, unclipped). */
+    linear: Float32Array;
+    provenance: ProcessingProvenance;
+}
+/** Converts linear RGB to unclipped float32 linear BT.2020 (relative values, no transfer function or quantization). */
+export declare function createHdrImageEncodingTask(input: HdrImageEncodingInput, provider: ProcessingProvider): ProcessingTask<HdrImageEncodingOutput>;