@photivra/engine 1.4.0 → 2.0.0

package/dist/output/sdr-rendering.js

@@ -1,39 +1,119 @@
 // SPDX-License-Identifier: Apache-2.0
 /**
- * Module boundary and integration notes.
- * Validates every explicit rendering/encoding policy; unknown HDR/display modes fail.
- * Validates declared image state and a bounded finite inline raster; never resolves WB/color.
+ * The SDR (standard dynamic range) reference renderer: linear sRGB → rendering exposure → tone
+ * curve → gamut handling → sRGB transfer → integer codes.
+ *
+ * It takes linear RGB that already has colour and white balance (WB) applied. Unknown HDR or
+ * display modes are rejected.
  * @see docs/SDR_RENDERING.md for equations, coordinate/unit conventions, blockers and support limits.
  */
 import { encodeSrgbComponent } from "../color/srgb-transfer.js";
+import { resolveCaptureColorModel } from "../color/capture-color.js";
 import { requireAllowlistedRecord, requirePublicOpaqueId } from "../core/record-validation.js";
 import { calculatedResult } from "../core/calculation-result.js";
 import { InvalidConfigurationError } from "../core/configuration-error.js";
 import { InvalidScientificResultError } from "../core/validation.js";
 export const SDR_RENDERING_SCHEMA_VERSION = "0.1.0";
+/**
+ * Version of the `generic-sdr-rendering-reference` model. 0.2.0 added luminance tone curves and fixed
+ * round-to-nearest, ties up (`Math.round`; the old `Math.floor(x+.5)` rounded 0.49999999999999994 up).
+ */
+export const SDR_RENDERING_MODEL_VERSION = "0.2.0";
 function object(v, keys) {
     return requireAllowlistedRecord(v, keys, "Invalid SDR fields.");
 }
 function publicId(v) {
     return requirePublicOpaqueId(v, "Invalid public SDR identity.");
 }
-/** Validates every explicit rendering/encoding policy; unknown HDR/display modes fail. */
+function bounded(value, minimumExclusive, maximum) {
+    if (typeof value !== "number" || !Number.isFinite(value) || value <= minimumExclusive || value > maximum) {
+        throw new InvalidConfigurationError("Unsupported SDR tone-curve parameters.");
+    }
+    return value;
+}
+function parseToneCurveParameters(curve, value) {
+    if (curve === "identity" || curve === "positive-reinhard-per-channel") {
+        if (value !== undefined)
+            throw new InvalidConfigurationError("Per-channel tone curves take no parameters.");
+        return undefined;
+    }
+    if (curve === "extended-reinhard-luminance") {
+        const p = requireAllowlistedRecord(value, ["whiteLuminance"], "Invalid SDR tone-curve parameters.");
+        return { whiteLuminance: bounded(p.whiteLuminance, 0, 65536) };
+    }
+    if (curve === "naka-rushton-luminance") {
+        const p = requireAllowlistedRecord(value, ["exponent", "semiSaturation", "whiteLuminance"], "Invalid SDR tone-curve parameters.");
+        const exponent = bounded(p.exponent, 0, 4);
+        if (exponent < .25)
+            throw new InvalidConfigurationError("Unsupported SDR tone-curve parameters.");
+        return { exponent, semiSaturation: bounded(p.semiSaturation, 0, 64), whiteLuminance: bounded(p.whiteLuminance, 0, 65536) };
+    }
+    throw new InvalidConfigurationError("Unsupported SDR rendering profile.");
+}
+/** Linear-sRGB/D65 relative-luminance weights: the Y row of the engine's sRGB to XYZ matrix (W3C CSS Color 4 sRGB; also IEC 61966-2-1). */
+const LUMINANCE = resolveCaptureColorModel().cameraRgbToXyz[1];
+/**
+ * Internal shared tone operator. The reference renderer, the native SDR kernel and RAW development all
+ * call this one function, so every path uses the same per-pixel operation order. Not exported from the package.
+ */
+export function prepareSdrPixelTone(profile) {
+    const curve = profile.toneCurve, p = profile.toneCurveParameters;
+    if (curve === "identity")
+        return (r, g, b, out) => { out[0] = r; out[1] = g; out[2] = b; };
+    if (curve === "positive-reinhard-per-channel")
+        return (r, g, b, out) => {
+            out[0] = r > 0 ? r / (1 + r) : r;
+            out[1] = g > 0 ? g / (1 + g) : g;
+            out[2] = b > 0 ? b / (1 + b) : b;
+        };
+    let map;
+    if (curve === "extended-reinhard-luminance" && p && !("exponent" in p)) {
+        const white2 = p.whiteLuminance * p.whiteLuminance;
+        // Written as L/(1+L)·(1+L/Lw²) so a very large finite L cannot overflow an intermediate product.
+        map = (l) => l / (1 + l) * (1 + l / white2);
+    }
+    else if (curve === "naka-rushton-luminance" && p && "exponent" in p) {
+        // g(L) = Lⁿ/(Lⁿ+σⁿ) = 1/(1+(σ/L)ⁿ): no overflow for large L.
+        const q = p;
+        const n = q.exponent, sigma = q.semiSaturation, g = (l) => 1 / (1 + (sigma / l) ** n), white = g(q.whiteLuminance);
+        map = (l) => g(l) / white;
+    }
+    else
+        throw new InvalidConfigurationError("Unsupported SDR rendering profile.");
+    const kr = LUMINANCE[0], kg = LUMINANCE[1], kb = LUMINANCE[2];
+    return (r, g, b, out) => {
+        const l = kr * r + kg * g + kb * b;
+        if (!(l > 0)) {
+            out[0] = r;
+            out[1] = g;
+            out[2] = b;
+            return;
+        }
+        const ratio = map(l) / l, mr = r * ratio, mg = g * ratio, mb = b * ratio;
+        if (!Number.isFinite(mr) || !Number.isFinite(mg) || !Number.isFinite(mb))
+            throw new InvalidScientificResultError("Tone mapping produced a non-finite sample.");
+        out[0] = mr;
+        out[1] = mg;
+        out[2] = mb;
+    };
+}
+/** Validates a rendering and encoding profile. Unknown HDR or display modes are rejected. */
 export function parseSdrRenderingProfile(value) {
-    const r = object(value, ["schemaVersion", "profileId", "profileVersion", "renderingExposureEv", "toneCurve", "gamutHandling",
+    const r = object(value, ["schemaVersion", "profileId", "profileVersion", "renderingExposureEv", "toneCurve", "toneCurveParameters", "gamutHandling",
         "outputDynamicRange", "transferFunction", "bitDepth", "rounding", "dither"]);
+    const toneCurveParameters = parseToneCurveParameters(r.toneCurve, r.toneCurveParameters);
     if (r.schemaVersion !== SDR_RENDERING_SCHEMA_VERSION || typeof r.renderingExposureEv !== "number" ||
         !Number.isFinite(r.renderingExposureEv) || Math.abs(r.renderingExposureEv) > 32 ||
-        (r.toneCurve !== "identity" && r.toneCurve !== "positive-reinhard-per-channel") ||
         (r.gamutHandling !== "clip-components" && r.gamutHandling !== "reject-out-of-range") ||
         r.outputDynamicRange !== "sdr" || r.transferFunction !== "srgb" || (r.bitDepth !== 8 && r.bitDepth !== 16) ||
         r.rounding !== "nearest-ties-up" || r.dither !== "none")
         throw new InvalidConfigurationError("Unsupported SDR rendering profile.");
     return { schemaVersion: SDR_RENDERING_SCHEMA_VERSION, profileId: publicId(r.profileId), profileVersion: publicId(r.profileVersion),
         renderingExposureEv: Object.is(r.renderingExposureEv, -0) ? 0 : r.renderingExposureEv, toneCurve: r.toneCurve,
-        gamutHandling: r.gamutHandling, outputDynamicRange: r.outputDynamicRange, transferFunction: r.transferFunction,
+        ...(toneCurveParameters === undefined ? {} : { toneCurveParameters }), gamutHandling: r.gamutHandling, outputDynamicRange: r.outputDynamicRange, transferFunction: r.transferFunction,
         bitDepth: r.bitDepth, rounding: r.rounding, dither: r.dither };
 }
-/** Validates declared image state and a bounded finite inline raster; never resolves WB/color. */
+/** Validates the declared image state and a bounded, finite, inline image. It never applies white balance or colour conversion. */
 export function parseSdrRenderingInput(value) {
     const r = object(value, ["sourceImageStateId", "inputImageState", "inputColorSpace", "whiteBalanceHandling",
         "pixelWidth", "pixelHeight", "referenceWhiteValue", "samples", "profile"]);
@@ -54,7 +134,8 @@ export function parseSdrRenderingInput(value) {
 }
 /**
  * Standalone SDR rendering: rendering exposure → tone → gamut → sRGB transfer → quantization.
- * No physical capture, metering, WB/color/correction integration or display adaptation runs.
+ * No physical capture, metering, white balance, colour conversion, lens correction or display
+ * adaptation runs.
  */
 export function calculateSdrRendering(input) {
     const value = parseSdrRenderingInput(input), p = value.profile;
@@ -62,19 +143,26 @@ export function calculateSdrRendering(input) {
     if (!Number.isFinite(scale) || scale <= 0)
         throw new InvalidScientificResultError("Rendering scale is non-finite/unrepresentable.");
     let negative = 0, above = 0, changed = 0, low = 0, high = 0;
-    const tone = value.samples.map((sample) => {
-        const exposed = sample * scale;
-        if (!Number.isFinite(exposed))
-            throw new InvalidScientificResultError("Rendering exposure produced a non-finite sample.");
-        if (exposed < 0)
-            negative++;
-        if (exposed > 1)
-            above++;
-        const mapped = p.toneCurve === "positive-reinhard-per-channel" && exposed > 0 ? exposed / (1 + exposed) : exposed;
-        if (mapped !== exposed)
-            changed++;
-        return Object.is(mapped, -0) ? 0 : mapped;
-    });
+    const toneOf = prepareSdrPixelTone(p), exposed = [0, 0, 0], mapped = [0, 0, 0], tone = [];
+    for (let i = 0; i < value.samples.length; i += 3) {
+        for (let c = 0; c < 3; c++) {
+            const e = value.samples[i + c] * scale;
+            if (!Number.isFinite(e))
+                throw new InvalidScientificResultError("Rendering exposure produced a non-finite sample.");
+            if (e < 0)
+                negative++;
+            if (e > 1)
+                above++;
+            exposed[c] = e;
+        }
+        toneOf(exposed[0], exposed[1], exposed[2], mapped);
+        for (let c = 0; c < 3; c++) {
+            const m = mapped[c];
+            if (m !== exposed[c])
+                changed++;
+            tone.push(m === 0 ? 0 : m);
+        }
+    }
     const linear = tone.map((sample) => {
         if (sample < 0 || sample > 1) {
             if (p.gamutHandling === "reject-out-of-range")
@@ -90,7 +178,8 @@ export function calculateSdrRendering(input) {
     });
     const encoded = linear.map(encodeSrgbComponent);
     const maximum = 2 ** p.bitDepth - 1;
-    const codes = encoded.map((sample) => Math.floor(sample * maximum + .5));
+    // Nearest, ties up: encoded samples are in [0,1], where Math.round is exact ties-up rounding.
+    const codes = encoded.map((sample) => Math.round(sample * maximum));
     return calculatedResult({ schemaVersion: SDR_RENDERING_SCHEMA_VERSION, sourceImageStateId: value.sourceImageStateId,
         profile: p, pixelWidth: value.pixelWidth, pixelHeight: value.pixelHeight, channelOrder: ["red", "green", "blue"],
         whiteBalanceHandling: value.whiteBalanceHandling, toneMappedLinearSamples: tone, outputLinearSamples: linear,
@@ -99,8 +188,8 @@ export function calculateSdrRendering(input) {
             whitePointXy: [.3127, .3290], referenceWhiteLuminanceCdM2: 80, codeMinimum: 0, codeMaximum: maximum },
         diagnostics: { renderingNegativeSampleCount: negative, renderingAboveReferenceSampleCount: above,
             toneChangedSampleCount: changed, gamutClippedLowSampleCount: low, gamutClippedHighSampleCount: high, captureSaturation: "not-consumed" },
-        displayAdaptation: { kind: "external-platform", applied: false, displayCapabilities: "unknown" } }, "generic-sdr-rendering-reference", SDR_RENDERING_SCHEMA_VERSION, ["Deterministic first-party rendering choice, not a calibrated camera look or scene-light model",
-        "Per-channel positive Reinhard can change chromaticity; component clipping is not perceptual gamut mapping",
+        displayAdaptation: { kind: "external-platform", applied: false, displayCapabilities: "unknown" } }, "generic-sdr-rendering-reference", SDR_RENDERING_MODEL_VERSION, ["Deterministic first-party rendering choice, not a calibrated camera look or scene-light model",
+        "Per-channel positive Reinhard can change chromaticity; luminance curves preserve chromaticity until component clipping, which is not perceptual gamut mapping",
         "Input WB/color completion is a producer declaration; capture/correction adapter and production activation remain pending",
         "sRGB reference encoding assumptions do not establish actual display luminance, saturation or viewing conditions"]);
 }

package/dist/output/sdr-rendering.js.map

@@ -1 +1 @@
-{"version":3,"file":"sdr-rendering.js","sourceRoot":"","sources":["../../src/output/sdr-rendering.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;GAKG;AAEH,OAAO,EAAE,mBAAmB,EAAE,MAAM,2BAA2B,CAAC;AAChE,OAAO,EAAE,wBAAwB,EAAE,qBAAqB,EAAE,MAAM,8BAA8B,CAAC;AAC/F,OAAO,EAAE,gBAAgB,EAA0B,MAAM,+BAA+B,CAAC;AACzF,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,4BAA4B,EAAE,MAAM,uBAAuB,CAAC;AAErE,MAAM,CAAC,MAAM,4BAA4B,GAAG,OAAgB,CAAC;AA4B7D,SAAS,MAAM,CAAC,CAAU,EAAE,IAAuB;IACjD,OAAO,wBAAwB,CAAC,CAAC,EAAE,IAAI,EAAE,qBAAqB,CAAC,CAAC;AAClE,CAAC;AACD,SAAS,QAAQ,CAAC,CAAU;IAC1B,OAAO,qBAAqB,CAAC,CAAC,EAAE,8BAA8B,CAAC,CAAC;AAClE,CAAC;AACD,0FAA0F;AAC1F,MAAM,UAAU,wBAAwB,CAAC,KAAc;IACrD,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,eAAe,EAAE,WAAW,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,WAAW,EAAE,eAAe;QAC1H,oBAAoB,EAAE,kBAAkB,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC,CAAC;IAC/E,IAAI,CAAC,CAAC,aAAa,KAAK,4BAA4B,IAAI,OAAO,CAAC,CAAC,mBAAmB,KAAK,QAAQ;QAC7F,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,mBAAmB,CAAC,GAAG,EAAE;QAC/E,CAAC,CAAC,CAAC,SAAS,KAAK,UAAU,IAAI,CAAC,CAAC,SAAS,KAAK,+BAA+B,CAAC;QAC/E,CAAC,CAAC,CAAC,aAAa,KAAK,iBAAiB,IAAI,CAAC,CAAC,aAAa,KAAK,qBAAqB,CAAC;QACpF,CAAC,CAAC,kBAAkB,KAAK,KAAK,IAAI,CAAC,CAAC,gBAAgB,KAAK,MAAM,IAAI,CAAC,CAAC,CAAC,QAAQ,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,EAAE,CAAC;QAC1G,CAAC,CAAC,QAAQ,KAAK,iBAAiB,IAAI,CAAC,CAAC,MAAM,KAAK,MAAM;QAAE,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;IACvI,OAAO,EAAE,aAAa,EAAE,4BAA4B,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,EAAE,cAAc,EAAE,QAAQ,CAAC,CAAC,CAAC,cAAc,CAAC;QAChI,mBAAmB,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,mBAAmB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,mBAAmB,EAAE,SAAS,EAAE,CAAC,CAAC,SAAS;QAC7G,aAAa,EAAE,CAAC,CAAC,aAAa,EAAE,kBAAkB,EAAE,CAAC,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,CAAC,gBAAgB;QAC9G,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC;AACnE,CAAC;AACD,kGAAkG;AAClG,MAAM,UAAU,sBAAsB,CAAC,KAAc;IACnD,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,oBAAoB,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,sBAAsB;QACzG,YAAY,EAAE,aAAa,EAAE,qBAAqB,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;IAC7E,MAAM,OAAO,GAAG,CAAC,CAAC,OAAO,CAAC;IAC1B,IAAI,CAAC,CAAC,eAAe,KAAK,8BAA8B,IAAI,CAAC,CAAC,eAAe,KAAK,iBAAiB;QAC/F,CAAC,CAAC,CAAC,oBAAoB,KAAK,0BAA0B,IAAI,CAAC,CAAC,oBAAoB,KAAK,cAAc,CAAC;QACpG,OAAO,CAAC,CAAC,UAAU,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,UAAU,IAAI,CAAC;QAC5F,OAAO,CAAC,CAAC,WAAW,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC,WAAW,IAAI,CAAC;QAC/F,OAAO,CAAC,CAAC,mBAAmB,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,mBAAmB,IAAI,CAAC;QAClH,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,CAAC,CAAC,UAAU,GAAC,CAAC,CAAC,WAAW,GAAC,CAAC,IAAI,OAAO,CAAC,MAAM,GAAG,MAAM;QACrG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC;QAC9E,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,sCAAsC,CAAC,CAAC;IACnJ,OAAO,EAAE,kBAAkB,EAAE,QAAQ,CAAC,CAAC,CAAC,kBAAkB,CAAC,EAAE,eAAe,EAAE,CAAC,CAAC,eAAe,EAAE,eAAe,EAAE,CAAC,CAAC,eAAe;QACjI,oBAAoB,EAAE,CAAC,CAAC,oBAAoB,EAAE,UAAU,EAAE,CAAC,CAAC,UAAU,EAAE,WAAW,EAAE,CAAC,CAAC,WAAW;QAClG,mBAAmB,EAAE,CAAC,CAAC,mBAAmB,EAAE,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzG,OAAO,EAAE,wBAAwB,CAAC,CAAC,CAAC,OAAO,CAAC,EAAE,CAAC;AACnD,CAAC;AAmCD;;;GAGG;AACH,MAAM,UAAU,qBAAqB,CAAC,KAAwB;IAC5D,MAAM,KAAK,GAAG,sBAAsB,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC;IAC/D,MAAM,KAAK,GAAG,CAAC,IAAE,CAAC,CAAC,mBAAmB,GAAC,KAAK,CAAC,mBAAmB,CAAC;IACjE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;QAAE,MAAM,IAAI,4BAA4B,CAAC,gDAAgD,CAAC,CAAC;IACpI,IAAI,QAAQ,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,CAAC;IAC5D,MAAM,IAAI,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QACxC,MAAM,OAAO,GAAG,MAAM,GAAC,KAAK,CAAC;QAC7B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC;YAAE,MAAM,IAAI,4BAA4B,CAAC,kDAAkD,CAAC,CAAC;QAC1H,IAAI,OAAO,GAAG,CAAC;YAAE,QAAQ,EAAE,CAAC;QAC5B,IAAI,OAAO,GAAG,CAAC;YAAE,KAAK,EAAE,CAAC;QACzB,MAAM,MAAM,GAAG,CAAC,CAAC,SAAS,KAAK,+BAA+B,IAAI,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,GAAC,CAAC,CAAC,GAAC,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC;QAC9G,IAAI,MAAM,KAAK,OAAO;YAAE,OAAO,EAAE,CAAC;QAClC,OAAO,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IAC5C,CAAC,CAAC,CAAC;IACH,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QACjC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,aAAa,KAAK,qBAAqB;gBAAE,MAAM,IAAI,yBAAyB,CAAC,mDAAmD,CAAC,CAAC;YACxI,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;gBAAC,GAAG,EAAE,CAAC;gBAAC,OAAO,CAAC,CAAC;YAAC,CAAC;YACpC,IAAI,EAAE,CAAC;YAAC,OAAO,CAAC,CAAC;QACnB,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC,CAAC;IACH,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,mBAAmB,CAAC,CAAC;IAChD,MAAM,OAAO,GAAG,CAAC,IAAE,CAAC,CAAC,QAAQ,GAAC,CAAC,CAAC;IAChC,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAC,OAAO,GAAC,EAAE,CAAC,CAAC,CAAC;IACrE,OAAO,gBAAgB,CAAC,EAAE,aAAa,EAAE,4BAA4B,EAAE,kBAAkB,EAAE,KAAK,CAAC,kBAAkB;QACjH,OAAO,EAAE,CAAC,EAAE,UAAU,EAAE,KAAK,CAAC,UAAU,EAAE,WAAW,EAAE,KAAK,CAAC,WAAW,EAAE,YAAY,EAAE,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC;QAChH,oBAAoB,EAAE,KAAK,CAAC,oBAAoB,EAAE,uBAAuB,EAAE,IAAI,EAAE,mBAAmB,EAAE,MAAM;QAC5G,oBAAoB,EAAE,OAAO,EAAE,cAAc,EAAE,KAAK;QACpD,cAAc,EAAE,EAAE,UAAU,EAAE,qBAAqB,EAAE,UAAU,EAAE,MAAM,EAAE,WAAW,EAAE,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,IAAI,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE;YAC5I,YAAY,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,2BAA2B,EAAE,EAAE,EAAE,WAAW,EAAE,CAAC,EAAE,WAAW,EAAE,OAAO,EAAE;QACvG,WAAW,EAAE,EAAE,4BAA4B,EAAE,QAAQ,EAAE,kCAAkC,EAAE,KAAK;YAC9F,sBAAsB,EAAE,OAAO,EAAE,0BAA0B,EAAE,GAAG,EAAE,2BAA2B,EAAE,IAAI,EAAE,iBAAiB,EAAE,cAAc,EAAE;QAC1I,iBAAiB,EAAE,EAAE,IAAI,EAAE,mBAAmB,EAAE,OAAO,EAAE,KAAK,EAAE,mBAAmB,EAAE,SAAS,EAAE,EAAE,EACpG,iCAAiC,EAAE,4BAA4B,EAC/D,CAAC,+FAA+F;QAC9F,2GAA2G;QAC3G,0HAA0H;QAC1H,iHAAiH,CAAC,CAAC,CAAC;AACxH,CAAC"}
+{"version":3,"file":"sdr-rendering.js","sourceRoot":"","sources":["../../src/output/sdr-rendering.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;GAOG;AAEH,OAAO,EAAE,mBAAmB,EAAE,MAAM,2BAA2B,CAAC;AAChE,OAAO,EAAE,wBAAwB,EAAE,MAAM,2BAA2B,CAAC;AACrE,OAAO,EAAE,wBAAwB,EAAE,qBAAqB,EAAE,MAAM,8BAA8B,CAAC;AAC/F,OAAO,EAAE,gBAAgB,EAA0B,MAAM,+BAA+B,CAAC;AACzF,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAC3E,OAAO,EAAE,4BAA4B,EAAE,MAAM,uBAAuB,CAAC;AAErE,MAAM,CAAC,MAAM,4BAA4B,GAAG,OAAgB,CAAC;AAC7D;;;GAGG;AACH,MAAM,CAAC,MAAM,2BAA2B,GAAG,OAAgB,CAAC;AA6C5D,SAAS,MAAM,CAAC,CAAU,EAAE,IAAuB;IACjD,OAAO,wBAAwB,CAAC,CAAC,EAAE,IAAI,EAAE,qBAAqB,CAAC,CAAC;AAClE,CAAC;AACD,SAAS,QAAQ,CAAC,CAAU;IAC1B,OAAO,qBAAqB,CAAC,CAAC,EAAE,8BAA8B,CAAC,CAAC;AAClE,CAAC;AACD,SAAS,OAAO,CAAC,KAAc,EAAE,gBAAwB,EAAE,OAAe;IACxE,IAAI,OAAO,KAAK,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,gBAAgB,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;QACzG,MAAM,IAAI,yBAAyB,CAAC,wCAAwC,CAAC,CAAC;IAChF,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AACD,SAAS,wBAAwB,CAAC,KAAc,EAAE,KAAc;IAC9D,IAAI,KAAK,KAAK,UAAU,IAAI,KAAK,KAAK,+BAA+B,EAAE,CAAC;QACtE,IAAI,KAAK,KAAK,SAAS;YAAE,MAAM,IAAI,yBAAyB,CAAC,6CAA6C,CAAC,CAAC;QAC5G,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,IAAI,KAAK,KAAK,6BAA6B,EAAE,CAAC;QAC5C,MAAM,CAAC,GAAG,wBAAwB,CAAC,KAAK,EAAE,CAAC,gBAAgB,CAAC,EAAE,oCAAoC,CAAC,CAAC;QACpG,OAAO,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC,CAAC,cAAc,EAAE,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;IACjE,CAAC;IACD,IAAI,KAAK,KAAK,wBAAwB,EAAE,CAAC;QACvC,MAAM,CAAC,GAAG,wBAAwB,CAAC,KAAK,EAAE,CAAC,UAAU,EAAE,gBAAgB,EAAE,gBAAgB,CAAC,EAAE,oCAAoC,CAAC,CAAC;QAClI,MAAM,QAAQ,GAAG,OAAO,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC3C,IAAI,QAAQ,GAAG,GAAG;YAAE,MAAM,IAAI,yBAAyB,CAAC,wCAAwC,CAAC,CAAC;QAClG,OAAO,EAAE,QAAQ,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC,CAAC,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,cAAc,EAAE,OAAO,CAAC,CAAC,CAAC,cAAc,EAAE,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;IAC7H,CAAC;IACD,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;AAC5E,CAAC;AACD,2IAA2I;AAC3I,MAAM,SAAS,GAAG,wBAAwB,EAAE,CAAC,cAAc,CAAC,CAAC,CAAuB,CAAC;AAGrF;;;GAGG;AACH,MAAM,UAAU,mBAAmB,CAAC,OAAuE;IACzG,MAAM,KAAK,GAAG,OAAO,CAAC,SAAS,EAAE,CAAC,GAAG,OAAO,CAAC,mBAAmB,CAAC;IACjE,IAAI,KAAK,KAAK,UAAU;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3F,IAAI,KAAK,KAAK,+BAA+B;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE;YACrE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3F,CAAC,CAAC;IACF,IAAI,GAA0B,CAAC;IAC/B,IAAI,KAAK,KAAK,6BAA6B,IAAI,CAAC,IAAI,CAAC,CAAC,UAAU,IAAI,CAAC,CAAC,EAAE,CAAC;QACvE,MAAM,MAAM,GAAG,CAAC,CAAC,cAAc,GAAC,CAAC,CAAC,cAAc,CAAC;QACjD,iGAAiG;QACjG,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,GAAC,MAAM,CAAC,CAAC;IACpD,CAAC;SAAM,IAAI,KAAK,KAAK,wBAAwB,IAAI,CAAC,IAAI,UAAU,IAAI,CAAC,EAAE,CAAC;QACtE,6DAA6D;QAC7D,MAAM,CAAC,GAAG,CAA0D,CAAC;QACrE,MAAM,CAAC,GAAG,CAAC,CAAC,QAAQ,EAAE,KAAK,GAAG,CAAC,CAAC,cAAc,EAAE,CAAC,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,KAAK,GAAC,CAAC,CAAC,IAAE,CAAC,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC;QAC3H,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,KAAK,CAAC;IAC1C,CAAC;;QAAM,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;IACjF,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,SAAS,CAAC,CAAC,CAAE,EAAE,EAAE,GAAG,SAAS,CAAC,CAAC,CAAE,CAAC;IACjE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE;QACtB,MAAM,CAAC,GAAG,EAAE,GAAC,CAAC,GAAC,EAAE,GAAC,CAAC,GAAC,EAAE,GAAC,CAAC,CAAC;QACzB,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAAC,OAAO;QAAC,CAAC;QAC7D,MAAM,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,GAAC,CAAC,EAAE,EAAE,GAAG,CAAC,GAAC,KAAK,EAAE,EAAE,GAAG,CAAC,GAAC,KAAK,EAAE,EAAE,GAAG,CAAC,GAAC,KAAK,CAAC;QACjE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;YAAE,MAAM,IAAI,4BAA4B,CAAC,4CAA4C,CAAC,CAAC;QAC/J,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACxC,CAAC,CAAC;AACJ,CAAC;AACD,6FAA6F;AAC7F,MAAM,UAAU,wBAAwB,CAAC,KAAc;IACrD,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,eAAe,EAAE,WAAW,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,WAAW,EAAE,qBAAqB,EAAE,eAAe;QACjJ,oBAAoB,EAAE,kBAAkB,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC,CAAC;IAC/E,MAAM,mBAAmB,GAAG,wBAAwB,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC,mBAAmB,CAAC,CAAC;IACzF,IAAI,CAAC,CAAC,aAAa,KAAK,4BAA4B,IAAI,OAAO,CAAC,CAAC,mBAAmB,KAAK,QAAQ;QAC7F,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,mBAAmB,CAAC,GAAG,EAAE;QAC/E,CAAC,CAAC,CAAC,aAAa,KAAK,iBAAiB,IAAI,CAAC,CAAC,aAAa,KAAK,qBAAqB,CAAC;QACpF,CAAC,CAAC,kBAAkB,KAAK,KAAK,IAAI,CAAC,CAAC,gBAAgB,KAAK,MAAM,IAAI,CAAC,CAAC,CAAC,QAAQ,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,EAAE,CAAC;QAC1G,CAAC,CAAC,QAAQ,KAAK,iBAAiB,IAAI,CAAC,CAAC,MAAM,KAAK,MAAM;QAAE,MAAM,IAAI,yBAAyB,CAAC,oCAAoC,CAAC,CAAC;IACvI,OAAO,EAAE,aAAa,EAAE,4BAA4B,EAAE,SAAS,EAAE,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,EAAE,cAAc,EAAE,QAAQ,CAAC,CAAC,CAAC,cAAc,CAAC;QAChI,mBAAmB,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,mBAAmB,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,mBAAmB,EAAE,SAAS,EAAE,CAAC,CAAC,SAAyB;QAC7H,GAAG,CAAC,mBAAmB,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,mBAAmB,EAAE,CAAC,EAAE,aAAa,EAAE,CAAC,CAAC,aAAa,EAAE,kBAAkB,EAAE,CAAC,CAAC,kBAAkB,EAAE,gBAAgB,EAAE,CAAC,CAAC,gBAAgB;QACrL,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC;AACnE,CAAC;AACD,mIAAmI;AACnI,MAAM,UAAU,sBAAsB,CAAC,KAAc;IACnD,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,oBAAoB,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,sBAAsB;QACzG,YAAY,EAAE,aAAa,EAAE,qBAAqB,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC,CAAC;IAC7E,MAAM,OAAO,GAAG,CAAC,CAAC,OAAO,CAAC;IAC1B,IAAI,CAAC,CAAC,eAAe,KAAK,8BAA8B,IAAI,CAAC,CAAC,eAAe,KAAK,iBAAiB;QAC/F,CAAC,CAAC,CAAC,oBAAoB,KAAK,0BAA0B,IAAI,CAAC,CAAC,oBAAoB,KAAK,cAAc,CAAC;QACpG,OAAO,CAAC,CAAC,UAAU,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,UAAU,IAAI,CAAC;QAC5F,OAAO,CAAC,CAAC,WAAW,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC,WAAW,IAAI,CAAC;QAC/F,OAAO,CAAC,CAAC,mBAAmB,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC,mBAAmB,IAAI,CAAC;QAClH,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,MAAM,KAAK,CAAC,CAAC,UAAU,GAAC,CAAC,CAAC,WAAW,GAAC,CAAC,IAAI,OAAO,CAAC,MAAM,GAAG,MAAM;QACrG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,OAAO,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC;QAC9E,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,yBAAyB,CAAC,sCAAsC,CAAC,CAAC;IACnJ,OAAO,EAAE,kBAAkB,EAAE,QAAQ,CAAC,CAAC,CAAC,kBAAkB,CAAC,EAAE,eAAe,EAAE,CAAC,CAAC,eAAe,EAAE,eAAe,EAAE,CAAC,CAAC,eAAe;QACjI,oBAAoB,EAAE,CAAC,CAAC,oBAAoB,EAAE,UAAU,EAAE,CAAC,CAAC,UAAU,EAAE,WAAW,EAAE,CAAC,CAAC,WAAW;QAClG,mBAAmB,EAAE,CAAC,CAAC,mBAAmB,EAAE,OAAO,EAAE,OAAO,CAAC,GAAG,CAAC,CAAC,CAAS,EAAE,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzG,OAAO,EAAE,wBAAwB,CAAC,CAAC,CAAC,OAAO,CAAC,EAAE,CAAC;AACnD,CAAC;AAwCD;;;;GAIG;AACH,MAAM,UAAU,qBAAqB,CAAC,KAAwB;IAC5D,MAAM,KAAK,GAAG,sBAAsB,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC;IAC/D,MAAM,KAAK,GAAG,CAAC,IAAE,CAAC,CAAC,mBAAmB,GAAC,KAAK,CAAC,mBAAmB,CAAC;IACjE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC;QAAE,MAAM,IAAI,4BAA4B,CAAC,gDAAgD,CAAC,CAAC;IACpI,IAAI,QAAQ,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,CAAC;IAC5D,MAAM,MAAM,GAAG,mBAAmB,CAAC,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,IAAI,GAAa,EAAE,CAAC;IACpG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;QACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,KAAK,CAAC,OAAO,CAAC,CAAC,GAAC,CAAC,CAAE,GAAC,KAAK,CAAC;YACpC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;gBAAE,MAAM,IAAI,4BAA4B,CAAC,kDAAkD,CAAC,CAAC;YACpH,IAAI,CAAC,GAAG,CAAC;gBAAE,QAAQ,EAAE,CAAC;YACtB,IAAI,CAAC,GAAG,CAAC;gBAAE,KAAK,EAAE,CAAC;YACnB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,CAAC;QACD,MAAM,CAAC,OAAO,CAAC,CAAC,CAAE,EAAE,OAAO,CAAC,CAAC,CAAE,EAAE,OAAO,CAAC,CAAC,CAAE,EAAE,MAAM,CAAC,CAAC;QACtD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC;YACrB,IAAI,CAAC,KAAK,OAAO,CAAC,CAAC,CAAC;gBAAE,OAAO,EAAE,CAAC;YAChC,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,CAAC;IACH,CAAC;IACD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE;QACjC,IAAI,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,aAAa,KAAK,qBAAqB;gBAAE,MAAM,IAAI,yBAAyB,CAAC,mDAAmD,CAAC,CAAC;YACxI,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;gBAAC,GAAG,EAAE,CAAC;gBAAC,OAAO,CAAC,CAAC;YAAC,CAAC;YACpC,IAAI,EAAE,CAAC;YAAC,OAAO,CAAC,CAAC;QACnB,CAAC;QACD,OAAO,MAAM,CAAC;IAChB,CAAC,CAAC,CAAC;IACH,MAAM,OAAO,GAAG,MAAM,CAAC,GAAG,CAAC,mBAAmB,CAAC,CAAC;IAChD,MAAM,OAAO,GAAG,CAAC,IAAE,CAAC,CAAC,QAAQ,GAAC,CAAC,CAAC;IAChC,8FAA8F;IAC9F,MAAM,KAAK,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAC,OAAO,CAAC,CAAC,CAAC;IAClE,OAAO,gBAAgB,CAAC,EAAE,aAAa,EAAE,4BAA4B,EAAE,kBAAkB,EAAE,KAAK,CAAC,kBAAkB;QACjH,OAAO,EAAE,CAAC,EAAE,UAAU,EAAE,KAAK,CAAC,UAAU,EAAE,WAAW,EAAE,KAAK,CAAC,WAAW,EAAE,YAAY,EAAE,CAAC,KAAK,EAAE,OAAO,EAAE,MAAM,CAAC;QAChH,oBAAoB,EAAE,KAAK,CAAC,oBAAoB,EAAE,uBAAuB,EAAE,IAAI,EAAE,mBAAmB,EAAE,MAAM;QAC5G,oBAAoB,EAAE,OAAO,EAAE,cAAc,EAAE,KAAK;QACpD,cAAc,EAAE,EAAE,UAAU,EAAE,qBAAqB,EAAE,UAAU,EAAE,MAAM,EAAE,WAAW,EAAE,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,IAAI,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE;YAC5I,YAAY,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,EAAE,2BAA2B,EAAE,EAAE,EAAE,WAAW,EAAE,CAAC,EAAE,WAAW,EAAE,OAAO,EAAE;QACvG,WAAW,EAAE,EAAE,4BAA4B,EAAE,QAAQ,EAAE,kCAAkC,EAAE,KAAK;YAC9F,sBAAsB,EAAE,OAAO,EAAE,0BAA0B,EAAE,GAAG,EAAE,2BAA2B,EAAE,IAAI,EAAE,iBAAiB,EAAE,cAAc,EAAE;QAC1I,iBAAiB,EAAE,EAAE,IAAI,EAAE,mBAAmB,EAAE,OAAO,EAAE,KAAK,EAAE,mBAAmB,EAAE,SAAS,EAAE,EAAE,EACpG,iCAAiC,EAAE,2BAA2B,EAC9D,CAAC,+FAA+F;QAC9F,+JAA+J;QAC/J,0HAA0H;QAC1H,iHAAiH,CAAC,CAAC,CAAC;AACxH,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * The SDR (standard dynamic range) reference renderer: linear sRGB → rendering exposure → tone\n * curve → gamut handling → sRGB transfer → integer codes.\n *\n * It takes linear RGB that already has colour and white balance (WB) applied. Unknown HDR or\n * display modes are rejected.\n * @see docs/SDR_RENDERING.md for equations, coordinate/unit conventions, blockers and support limits.\n */\n\nimport { encodeSrgbComponent } from \"../color/srgb-transfer.js\";\nimport { resolveCaptureColorModel } from \"../color/capture-color.js\";\nimport { requireAllowlistedRecord, requirePublicOpaqueId } from \"../core/record-validation.js\";\nimport { calculatedResult, type CalculationResult } from \"../core/calculation-result.js\";\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\nimport { InvalidScientificResultError } from \"../core/validation.js\";\n\nexport const SDR_RENDERING_SCHEMA_VERSION = \"0.1.0\" as const;\n/**\n * Version of the `generic-sdr-rendering-reference` model. 0.2.0 added luminance tone curves and fixed\n * round-to-nearest, ties up (`Math.round`; the old `Math.floor(x+.5)` rounded 0.49999999999999994 up).\n */\nexport const SDR_RENDERING_MODEL_VERSION = \"0.2.0\" as const;\n/**\n * Tone curves. A per-channel curve acts on each exposed component separately. A luminance curve\n * scales all three components of a pixel by f(L)/L, where L is the linear-sRGB/D65 relative\n * luminance, so the colour (chromaticity) is kept before gamut handling.\n * - `identity`: f(x)=x.\n * - `positive-reinhard-per-channel`: f(x)=x/(1+x) for x>0 (Reinhard et al. 2002, eq. 3), per channel.\n * - `extended-reinhard-luminance`: f(L)=L(1+L/Lw²)/(1+L) (Reinhard et al. 2002, eq. 4); L=Lw maps to 1.\n * - `naka-rushton-luminance`: f(L)=g(L)/g(Lw), g(L)=Lⁿ/(Lⁿ+σⁿ) (Naka & Rushton 1966; Reinhard & Devlin\n *   2005): a sigmoid in log luminance with contrast n and mid-response σ, i.e. a film-like S-curve.\n */\nexport type SdrToneCurve = \"identity\" | \"positive-reinhard-per-channel\" | \"extended-reinhard-luminance\" | \"naka-rushton-luminance\";\n/** Required only by luminance curves; all values are relative linear luminance after rendering exposure. */\nexport type SdrToneCurveParameters =\n  | { whiteLuminance: number }\n  | { exponent: number; semiSaturation: number; whiteLuminance: number };\n/** Rendering choices made by the engine user. They are independent of physical exposure and display calibration. */\nexport interface SdrRenderingProfile {\n  schemaVersion: typeof SDR_RENDERING_SCHEMA_VERSION;\n  profileId: string;\n  profileVersion: string;\n  renderingExposureEv: number;\n  toneCurve: SdrToneCurve;\n  /** Present exactly for `extended-reinhard-luminance` ({whiteLuminance}) and `naka-rushton-luminance`. */\n  toneCurveParameters?: SdrToneCurveParameters;\n  gamutHandling: \"clip-components\" | \"reject-out-of-range\";\n  outputDynamicRange: \"sdr\";\n  transferFunction: \"srgb\";\n  bitDepth: 8 | 16;\n  rounding: \"nearest-ties-up\";\n  dither: \"none\";\n}\n/** Only linear-sRGB/D65 data that has been explicitly colour-transformed, with all three channels at each pixel, is accepted. */\nexport interface SdrRenderingInput {\n  sourceImageStateId: string;\n  inputImageState: \"color-transformed-linear-rgb\";\n  inputColorSpace: \"linear-srgb-d65\";\n  whiteBalanceHandling: \"already-applied-upstream\" | \"not-required\";\n  pixelWidth: number;\n  pixelHeight: number;\n  referenceWhiteValue: number;\n  /** Relative linear RGB, row-major tightly interleaved red/green/blue. */\n  samples: readonly number[];\n  profile: SdrRenderingProfile;\n}\nfunction object(v: unknown, keys: readonly string[]): Record<string, unknown> {\n  return requireAllowlistedRecord(v, keys, \"Invalid SDR fields.\");\n}\nfunction publicId(v: unknown): string {\n  return requirePublicOpaqueId(v, \"Invalid public SDR identity.\");\n}\nfunction bounded(value: unknown, minimumExclusive: number, maximum: number): number {\n  if (typeof value !== \"number\" || !Number.isFinite(value) || value <= minimumExclusive || value > maximum) {\n    throw new InvalidConfigurationError(\"Unsupported SDR tone-curve parameters.\");\n  }\n  return value;\n}\nfunction parseToneCurveParameters(curve: unknown, value: unknown): SdrToneCurveParameters | undefined {\n  if (curve === \"identity\" || curve === \"positive-reinhard-per-channel\") {\n    if (value !== undefined) throw new InvalidConfigurationError(\"Per-channel tone curves take no parameters.\");\n    return undefined;\n  }\n  if (curve === \"extended-reinhard-luminance\") {\n    const p = requireAllowlistedRecord(value, [\"whiteLuminance\"], \"Invalid SDR tone-curve parameters.\");\n    return { whiteLuminance: bounded(p.whiteLuminance, 0, 65536) };\n  }\n  if (curve === \"naka-rushton-luminance\") {\n    const p = requireAllowlistedRecord(value, [\"exponent\", \"semiSaturation\", \"whiteLuminance\"], \"Invalid SDR tone-curve parameters.\");\n    const exponent = bounded(p.exponent, 0, 4);\n    if (exponent < .25) throw new InvalidConfigurationError(\"Unsupported SDR tone-curve parameters.\");\n    return { exponent, semiSaturation: bounded(p.semiSaturation, 0, 64), whiteLuminance: bounded(p.whiteLuminance, 0, 65536) };\n  }\n  throw new InvalidConfigurationError(\"Unsupported SDR rendering profile.\");\n}\n/** Linear-sRGB/D65 relative-luminance weights: the Y row of the engine's sRGB to XYZ matrix (W3C CSS Color 4 sRGB; also IEC 61966-2-1). */\nconst LUMINANCE = resolveCaptureColorModel().cameraRgbToXyz[1]! as readonly number[];\n/** Writes the three tone-mapped samples of one exposed pixel. Exposed values are finite. */\nexport type SdrPixelTone = (red: number, green: number, blue: number, out: Float64Array | number[]) => void;\n/**\n * Internal shared tone operator. The reference renderer, the native SDR kernel and RAW development all\n * call this one function, so every path uses the same per-pixel operation order. Not exported from the package.\n */\nexport function prepareSdrPixelTone(profile: Pick<SdrRenderingProfile, \"toneCurve\" | \"toneCurveParameters\">): SdrPixelTone {\n  const curve = profile.toneCurve, p = profile.toneCurveParameters;\n  if (curve === \"identity\") return (r, g, b, out) => { out[0] = r; out[1] = g; out[2] = b; };\n  if (curve === \"positive-reinhard-per-channel\") return (r, g, b, out) => {\n    out[0] = r > 0 ? r/(1+r) : r; out[1] = g > 0 ? g/(1+g) : g; out[2] = b > 0 ? b/(1+b) : b;\n  };\n  let map: (l: number) => number;\n  if (curve === \"extended-reinhard-luminance\" && p && !(\"exponent\" in p)) {\n    const white2 = p.whiteLuminance*p.whiteLuminance;\n    // Written as L/(1+L)·(1+L/Lw²) so a very large finite L cannot overflow an intermediate product.\n    map = (l: number): number => l/(1+l)*(1+l/white2);\n  } else if (curve === \"naka-rushton-luminance\" && p && \"exponent\" in p) {\n    // g(L) = Lⁿ/(Lⁿ+σⁿ) = 1/(1+(σ/L)ⁿ): no overflow for large L.\n    const q = p as Extract<SdrToneCurveParameters, { exponent: number }>;\n    const n = q.exponent, sigma = q.semiSaturation, g = (l: number): number => 1/(1+(sigma/l)**n), white = g(q.whiteLuminance);\n    map = (l: number): number => g(l)/white;\n  } else throw new InvalidConfigurationError(\"Unsupported SDR rendering profile.\");\n  const kr = LUMINANCE[0]!, kg = LUMINANCE[1]!, kb = LUMINANCE[2]!;\n  return (r, g, b, out) => {\n    const l = kr*r+kg*g+kb*b;\n    if (!(l > 0)) { out[0] = r; out[1] = g; out[2] = b; return; }\n    const ratio = map(l)/l, mr = r*ratio, mg = g*ratio, mb = b*ratio;\n    if (!Number.isFinite(mr) || !Number.isFinite(mg) || !Number.isFinite(mb)) throw new InvalidScientificResultError(\"Tone mapping produced a non-finite sample.\");\n    out[0] = mr; out[1] = mg; out[2] = mb;\n  };\n}\n/** Validates a rendering and encoding profile. Unknown HDR or display modes are rejected. */\nexport function parseSdrRenderingProfile(value: unknown): SdrRenderingProfile {\n  const r = object(value, [\"schemaVersion\", \"profileId\", \"profileVersion\", \"renderingExposureEv\", \"toneCurve\", \"toneCurveParameters\", \"gamutHandling\",\n    \"outputDynamicRange\", \"transferFunction\", \"bitDepth\", \"rounding\", \"dither\"]);\n  const toneCurveParameters = parseToneCurveParameters(r.toneCurve, r.toneCurveParameters);\n  if (r.schemaVersion !== SDR_RENDERING_SCHEMA_VERSION || typeof r.renderingExposureEv !== \"number\" ||\n      !Number.isFinite(r.renderingExposureEv) || Math.abs(r.renderingExposureEv) > 32 ||\n      (r.gamutHandling !== \"clip-components\" && r.gamutHandling !== \"reject-out-of-range\") ||\n      r.outputDynamicRange !== \"sdr\" || r.transferFunction !== \"srgb\" || (r.bitDepth !== 8 && r.bitDepth !== 16) ||\n      r.rounding !== \"nearest-ties-up\" || r.dither !== \"none\") throw new InvalidConfigurationError(\"Unsupported SDR rendering profile.\");\n  return { schemaVersion: SDR_RENDERING_SCHEMA_VERSION, profileId: publicId(r.profileId), profileVersion: publicId(r.profileVersion),\n    renderingExposureEv: Object.is(r.renderingExposureEv, -0) ? 0 : r.renderingExposureEv, toneCurve: r.toneCurve as SdrToneCurve,\n    ...(toneCurveParameters === undefined ? {} : { toneCurveParameters }), gamutHandling: r.gamutHandling, outputDynamicRange: r.outputDynamicRange, transferFunction: r.transferFunction,\n    bitDepth: r.bitDepth, rounding: r.rounding, dither: r.dither };\n}\n/** Validates the declared image state and a bounded, finite, inline image. It never applies white balance or colour conversion. */\nexport function parseSdrRenderingInput(value: unknown): SdrRenderingInput {\n  const r = object(value, [\"sourceImageStateId\", \"inputImageState\", \"inputColorSpace\", \"whiteBalanceHandling\",\n    \"pixelWidth\", \"pixelHeight\", \"referenceWhiteValue\", \"samples\", \"profile\"]);\n  const samples = r.samples;\n  if (r.inputImageState !== \"color-transformed-linear-rgb\" || r.inputColorSpace !== \"linear-srgb-d65\" ||\n      (r.whiteBalanceHandling !== \"already-applied-upstream\" && r.whiteBalanceHandling !== \"not-required\") ||\n      typeof r.pixelWidth !== \"number\" || !Number.isSafeInteger(r.pixelWidth) || r.pixelWidth <= 0 ||\n      typeof r.pixelHeight !== \"number\" || !Number.isSafeInteger(r.pixelHeight) || r.pixelHeight <= 0 ||\n      typeof r.referenceWhiteValue !== \"number\" || !Number.isFinite(r.referenceWhiteValue) || r.referenceWhiteValue <= 0 ||\n      !Array.isArray(samples) || samples.length !== r.pixelWidth*r.pixelHeight*3 || samples.length > 262144 ||\n      Array.from({ length: samples.length }, (_, i) => i in samples).includes(false) ||\n      samples.some((v) => typeof v !== \"number\" || !Number.isFinite(v))) throw new InvalidConfigurationError(\"Unsupported SDR input domain/raster.\");\n  return { sourceImageStateId: publicId(r.sourceImageStateId), inputImageState: r.inputImageState, inputColorSpace: r.inputColorSpace,\n    whiteBalanceHandling: r.whiteBalanceHandling, pixelWidth: r.pixelWidth, pixelHeight: r.pixelHeight,\n    referenceWhiteValue: r.referenceWhiteValue, samples: samples.map((v: number) => Object.is(v, -0) ? 0 : v),\n    profile: parseSdrRenderingProfile(r.profile) };\n}\nexport interface SdrRenderingResult {\n  schemaVersion: typeof SDR_RENDERING_SCHEMA_VERSION;\n  sourceImageStateId: string;\n  profile: SdrRenderingProfile;\n  pixelWidth: number;\n  pixelHeight: number;\n  channelOrder: readonly [\"red\", \"green\", \"blue\"];\n  whiteBalanceHandling: SdrRenderingInput[\"whiteBalanceHandling\"];\n  /** Relative normalized linear samples after exposure/tone, before gamut handling. */\n  toneMappedLinearSamples: readonly number[];\n  outputLinearSamples: readonly number[];\n  outputEncodedSamples: readonly number[];\n  integerSamples: readonly number[];\n  outputEncoding: {\n    imageState: \"output-referred-sdr\";\n    colorSpace: \"srgb\";\n    primariesXy: { red: readonly [number, number]; green: readonly [number, number]; blue: readonly [number, number] };\n    whitePointXy: readonly [number, number];\n    referenceWhiteLuminanceCdM2: 80;\n    codeMinimum: 0;\n    codeMaximum: number;\n  };\n  diagnostics: {\n    renderingNegativeSampleCount: number;\n    renderingAboveReferenceSampleCount: number;\n    toneChangedSampleCount: number;\n    gamutClippedLowSampleCount: number;\n    gamutClippedHighSampleCount: number;\n    /**\n     * `not-consumed`: no upstream saturation data (full well, ADC or RAW) is available, and it is never\n     * guessed from the output.\n     * `consumed-in-raw-development`: RAW development already clipped or reconstructed flagged and\n     * white-level samples before white balance (see its highlight diagnostics).\n     */\n    captureSaturation: \"not-consumed\" | \"consumed-in-raw-development\";\n  };\n  /** Describes the encoded output for the platform's colour management. No real display is assessed. */\n  displayAdaptation: { kind: \"external-platform\"; applied: false; displayCapabilities: \"unknown\" };\n}\n/**\n * Standalone SDR rendering: rendering exposure → tone → gamut → sRGB transfer → quantization.\n * No physical capture, metering, white balance, colour conversion, lens correction or display\n * adaptation runs.\n */\nexport function calculateSdrRendering(input: SdrRenderingInput): CalculationResult<SdrRenderingResult> {\n  const value = parseSdrRenderingInput(input), p = value.profile;\n  const scale = 2**p.renderingExposureEv/value.referenceWhiteValue;\n  if (!Number.isFinite(scale) || scale <= 0) throw new InvalidScientificResultError(\"Rendering scale is non-finite/unrepresentable.\");\n  let negative = 0, above = 0, changed = 0, low = 0, high = 0;\n  const toneOf = prepareSdrPixelTone(p), exposed = [0, 0, 0], mapped = [0, 0, 0], tone: number[] = [];\n  for (let i = 0; i < value.samples.length; i += 3) {\n    for (let c = 0; c < 3; c++) {\n      const e = value.samples[i+c]!*scale;\n      if (!Number.isFinite(e)) throw new InvalidScientificResultError(\"Rendering exposure produced a non-finite sample.\");\n      if (e < 0) negative++;\n      if (e > 1) above++;\n      exposed[c] = e;\n    }\n    toneOf(exposed[0]!, exposed[1]!, exposed[2]!, mapped);\n    for (let c = 0; c < 3; c++) {\n      const m = mapped[c]!;\n      if (m !== exposed[c]) changed++;\n      tone.push(m === 0 ? 0 : m);\n    }\n  }\n  const linear = tone.map((sample) => {\n    if (sample < 0 || sample > 1) {\n      if (p.gamutHandling === \"reject-out-of-range\") throw new InvalidConfigurationError(\"Tone output exceeds the declared SDR gamut range.\");\n      if (sample < 0) { low++; return 0; }\n      high++; return 1;\n    }\n    return sample;\n  });\n  const encoded = linear.map(encodeSrgbComponent);\n  const maximum = 2**p.bitDepth-1;\n  // Nearest, ties up: encoded samples are in [0,1], where Math.round is exact ties-up rounding.\n  const codes = encoded.map((sample) => Math.round(sample*maximum));\n  return calculatedResult({ schemaVersion: SDR_RENDERING_SCHEMA_VERSION, sourceImageStateId: value.sourceImageStateId,\n    profile: p, pixelWidth: value.pixelWidth, pixelHeight: value.pixelHeight, channelOrder: [\"red\", \"green\", \"blue\"],\n    whiteBalanceHandling: value.whiteBalanceHandling, toneMappedLinearSamples: tone, outputLinearSamples: linear,\n    outputEncodedSamples: encoded, integerSamples: codes,\n    outputEncoding: { imageState: \"output-referred-sdr\", colorSpace: \"srgb\", primariesXy: { red: [.64, .33], green: [.30, .60], blue: [.15, .06] },\n      whitePointXy: [.3127, .3290], referenceWhiteLuminanceCdM2: 80, codeMinimum: 0, codeMaximum: maximum },\n    diagnostics: { renderingNegativeSampleCount: negative, renderingAboveReferenceSampleCount: above,\n      toneChangedSampleCount: changed, gamutClippedLowSampleCount: low, gamutClippedHighSampleCount: high, captureSaturation: \"not-consumed\" },\n    displayAdaptation: { kind: \"external-platform\", applied: false, displayCapabilities: \"unknown\" } },\n  \"generic-sdr-rendering-reference\", SDR_RENDERING_MODEL_VERSION,\n  [\"Deterministic first-party rendering choice, not a calibrated camera look or scene-light model\",\n    \"Per-channel positive Reinhard can change chromaticity; luminance curves preserve chromaticity until component clipping, which is not perceptual gamut mapping\",\n    \"Input WB/color completion is a producer declaration; capture/correction adapter and production activation remain pending\",\n    \"sRGB reference encoding assumptions do not establish actual display luminance, saturation or viewing conditions\"]);\n}\n"]}

package/dist/output/subject-framing-crop.d.ts

@@ -1,8 +1,6 @@
 /**
- * Module boundary and integration notes.
- * Calculates the same-aspect-ratio crop needed for a subject to occupy a requested fraction of frame
- * height. This calculation assumes a centered/positionable crop with enough spatial margin around the
- * subject. It does not validate subject position against image edges.
+ * Subject framing: how much to crop (keeping the aspect ratio) so a subject fills a chosen share
+ * of the frame height.
  * @see docs/PROCESSED_OUTPUT.md for equations, coordinate/unit conventions, blockers and support
  * limits.
  */
@@ -33,12 +31,10 @@ export interface SubjectFramingCrop {
     cropped: boolean;
 }
 /**
- * Calculates the same-aspect-ratio crop needed for a subject to occupy a
- * requested fraction of frame height.
+ * Calculates the crop (keeping the aspect ratio) that makes a subject fill a
+ * requested fraction of the frame height.
  *
- * This calculation assumes a centered/positionable crop with enough spatial
- * margin around the subject. It does not validate subject position against
- * image edges.
+ * It assumes the crop can be centred or moved freely, with enough room around
+ * the subject. It does not check the subject's position against the image edges.
  */
 export declare function calculateSubjectFramingCrop(input: CalculateSubjectFramingCropInput): CalculationResult<SubjectFramingCrop>;
-//# sourceMappingURL=subject-framing-crop.d.ts.map

package/dist/output/subject-framing-crop.js

@@ -1,23 +1,23 @@
 // SPDX-License-Identifier: Apache-2.0
 /**
- * Module boundary and integration notes.
- * Calculates the same-aspect-ratio crop needed for a subject to occupy a requested fraction of frame
- * height. This calculation assumes a centered/positionable crop with enough spatial margin around the
- * subject. It does not validate subject position against image edges.
+ * Subject framing: how much to crop (keeping the aspect ratio) so a subject fills a chosen share
+ * of the frame height.
  * @see docs/PROCESSED_OUTPUT.md for equations, coordinate/unit conventions, blockers and support
  * limits.
  */
 import { calculatedResult } from "../core/calculation-result.js";
 import { InvalidScientificInputError, requirePositiveFinite, requirePositiveInteger } from "../core/validation.js";
+import { floorWholePixels } from "./crop.js";
+import { requireInputKeys } from "../core/record-validation.js";
 /**
- * Calculates the same-aspect-ratio crop needed for a subject to occupy a
- * requested fraction of frame height.
+ * Calculates the crop (keeping the aspect ratio) that makes a subject fill a
+ * requested fraction of the frame height.
  *
- * This calculation assumes a centered/positionable crop with enough spatial
- * margin around the subject. It does not validate subject position against
- * image edges.
+ * It assumes the crop can be centred or moved freely, with enough room around
+ * the subject. It does not check the subject's position against the image edges.
  */
 export function calculateSubjectFramingCrop(input) {
+    requireInputKeys(input, "calculateSubjectFramingCrop", ["pixelWidth", "pixelHeight", "subjectHeightPixels", "targetSubjectHeightFraction"]);
     requirePositiveInteger("pixelWidth", input.pixelWidth);
     requirePositiveInteger("pixelHeight", input.pixelHeight);
     requirePositiveFinite("subjectHeightPixels", input.subjectHeightPixels);
@@ -29,8 +29,18 @@ export function calculateSubjectFramingCrop(input) {
     const cropFactor = requestedCropHeight >= input.pixelHeight
         ? 1
         : input.pixelHeight / requestedCropHeight;
-    const pixelHeight = Math.max(1, Math.floor(input.pixelHeight / cropFactor));
-    const pixelWidth = Math.max(1, Math.floor(input.pixelWidth / cropFactor));
+    // Height is the limiting axis. Compute it straight from the requested crop
+    // height, not through H / (H / requested), whose rounding could lose a pixel.
+    // It is the largest whole-pixel height within the requested crop, so the
+    // subject fills at least the target fraction. It is never smaller than the
+    // subject, so rounding alone never clips the subject. Width follows from that
+    // integer height and the source aspect ratio.
+    const pixelHeight = cropFactor === 1
+        ? input.pixelHeight
+        : Math.min(input.pixelHeight, Math.max(floorWholePixels(requestedCropHeight), Math.ceil(input.subjectHeightPixels)));
+    const pixelWidth = cropFactor === 1
+        ? input.pixelWidth
+        : Math.max(1, Math.min(input.pixelWidth, Math.round((pixelHeight * input.pixelWidth) / input.pixelHeight)));
     const subjectHeightFraction = input.subjectHeightPixels / pixelHeight;
     return calculatedResult({
         cropFactor,
@@ -40,8 +50,8 @@ export function calculateSubjectFramingCrop(input) {
         subjectHeightFraction,
         subjectClipped: subjectHeightFraction > 1,
         cropped: cropFactor > 1
-    }, "subject-height-framing-crop", "1.0.0", [
-        "Original aspect ratio is preserved",
+    }, "subject-height-framing-crop", "2.0.0", [
+        "Crop height is the largest whole-pixel height within the requested crop, but never smaller than the subject; width is the nearest integer to the original aspect ratio",
         "Crop can be positioned around the subject",
         "Subject edge/margin constraints are not modeled",
         "subjectHeightFraction may exceed 1 when the subject is clipped by frame height"

package/dist/output/subject-framing-crop.js.map

@@ -1 +1 @@
-{"version":3,"file":"subject-framing-crop.js","sourceRoot":"","sources":["../../src/output/subject-framing-crop.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;GAOG;AAEH,OAAO,EACL,gBAAgB,EAEjB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EACL,2BAA2B,EAC3B,qBAAqB,EACrB,sBAAsB,EACvB,MAAM,uBAAuB,CAAC;AA6B/B;;;;;;;GAOG;AACH,MAAM,UAAU,2BAA2B,CACzC,KAAuC;IAEvC,sBAAsB,CAAC,YAAY,EAAE,KAAK,CAAC,UAAU,CAAC,CAAC;IACvD,sBAAsB,CAAC,aAAa,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC;IACzD,qBAAqB,CAAC,qBAAqB,EAAE,KAAK,CAAC,mBAAmB,CAAC,CAAC;IACxE,qBAAqB,CACnB,6BAA6B,EAC7B,KAAK,CAAC,2BAA2B,CAClC,CAAC;IAEF,IAAI,KAAK,CAAC,2BAA2B,GAAG,CAAC,EAAE,CAAC;QAC1C,MAAM,IAAI,2BAA2B,CACnC,8DAA8D,CAC/D,CAAC;IACJ,CAAC;IAED,MAAM,mBAAmB,GACvB,KAAK,CAAC,mBAAmB,GAAG,KAAK,CAAC,2BAA2B,CAAC;IAChE,MAAM,UAAU,GACd,mBAAmB,IAAI,KAAK,CAAC,WAAW;QACtC,CAAC,CAAC,CAAC;QACH,CAAC,CAAC,KAAK,CAAC,WAAW,GAAG,mBAAmB,CAAC;IAE9C,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAC1B,CAAC,EACD,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,WAAW,GAAG,UAAU,CAAC,CAC3C,CAAC;IACF,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CACzB,CAAC,EACD,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,UAAU,GAAG,UAAU,CAAC,CAC1C,CAAC;IACF,MAAM,qBAAqB,GACzB,KAAK,CAAC,mBAAmB,GAAG,WAAW,CAAC;IAE1C,OAAO,gBAAgB,CACrB;QACE,UAAU;QACV,UAAU;QACV,WAAW;QACX,UAAU,EAAE,CAAC,UAAU,GAAG,WAAW,CAAC,GAAG,SAAS;QAClD,qBAAqB;QACrB,cAAc,EAAE,qBAAqB,GAAG,CAAC;QACzC,OAAO,EAAE,UAAU,GAAG,CAAC;KACxB,EACD,6BAA6B,EAC7B,OAAO,EACP;QACE,oCAAoC;QACpC,2CAA2C;QAC3C,iDAAiD;QACjD,gFAAgF;KACjF,CACF,CAAC;AACJ,CAAC"}
+{"version":3,"file":"subject-framing-crop.js","sourceRoot":"","sources":["../../src/output/subject-framing-crop.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;GAKG;AAEH,OAAO,EACL,gBAAgB,EAEjB,MAAM,+BAA+B,CAAC;AACvC,OAAO,EACL,2BAA2B,EAC3B,qBAAqB,EACrB,sBAAsB,EACvB,MAAM,uBAAuB,CAAC;AAC/B,OAAO,EAAE,gBAAgB,EAAE,MAAM,WAAW,CAAC;AAC7C,OAAO,EAAE,gBAAgB,EAAE,MAAM,8BAA8B,CAAC;AA6BhE;;;;;;GAMG;AACH,MAAM,UAAU,2BAA2B,CACzC,KAAuC;IAEvC,gBAAgB,CAAC,KAAK,EAAE,6BAA6B,EAAE,CAAC,YAAY,EAAE,aAAa,EAAE,qBAAqB,EAAE,6BAA6B,CAAC,CAAC,CAAC;IAC5I,sBAAsB,CAAC,YAAY,EAAE,KAAK,CAAC,UAAU,CAAC,CAAC;IACvD,sBAAsB,CAAC,aAAa,EAAE,KAAK,CAAC,WAAW,CAAC,CAAC;IACzD,qBAAqB,CAAC,qBAAqB,EAAE,KAAK,CAAC,mBAAmB,CAAC,CAAC;IACxE,qBAAqB,CACnB,6BAA6B,EAC7B,KAAK,CAAC,2BAA2B,CAClC,CAAC;IAEF,IAAI,KAAK,CAAC,2BAA2B,GAAG,CAAC,EAAE,CAAC;QAC1C,MAAM,IAAI,2BAA2B,CACnC,8DAA8D,CAC/D,CAAC;IACJ,CAAC;IAED,MAAM,mBAAmB,GACvB,KAAK,CAAC,mBAAmB,GAAG,KAAK,CAAC,2BAA2B,CAAC;IAChE,MAAM,UAAU,GACd,mBAAmB,IAAI,KAAK,CAAC,WAAW;QACtC,CAAC,CAAC,CAAC;QACH,CAAC,CAAC,KAAK,CAAC,WAAW,GAAG,mBAAmB,CAAC;IAE9C,2EAA2E;IAC3E,8EAA8E;IAC9E,yEAAyE;IACzE,2EAA2E;IAC3E,8EAA8E;IAC9E,8CAA8C;IAC9C,MAAM,WAAW,GACf,UAAU,KAAK,CAAC;QACd,CAAC,CAAC,KAAK,CAAC,WAAW;QACnB,CAAC,CAAC,IAAI,CAAC,GAAG,CACN,KAAK,CAAC,WAAW,EACjB,IAAI,CAAC,GAAG,CACN,gBAAgB,CAAC,mBAAmB,CAAC,EACrC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,mBAAmB,CAAC,CACrC,CACF,CAAC;IACR,MAAM,UAAU,GACd,UAAU,KAAK,CAAC;QACd,CAAC,CAAC,KAAK,CAAC,UAAU;QAClB,CAAC,CAAC,IAAI,CAAC,GAAG,CACN,CAAC,EACD,IAAI,CAAC,GAAG,CACN,KAAK,CAAC,UAAU,EAChB,IAAI,CAAC,KAAK,CAAC,CAAC,WAAW,GAAG,KAAK,CAAC,UAAU,CAAC,GAAG,KAAK,CAAC,WAAW,CAAC,CACjE,CACF,CAAC;IACR,MAAM,qBAAqB,GACzB,KAAK,CAAC,mBAAmB,GAAG,WAAW,CAAC;IAE1C,OAAO,gBAAgB,CACrB;QACE,UAAU;QACV,UAAU;QACV,WAAW;QACX,UAAU,EAAE,CAAC,UAAU,GAAG,WAAW,CAAC,GAAG,SAAS;QAClD,qBAAqB;QACrB,cAAc,EAAE,qBAAqB,GAAG,CAAC;QACzC,OAAO,EAAE,UAAU,GAAG,CAAC;KACxB,EACD,6BAA6B,EAC7B,OAAO,EACP;QACE,wKAAwK;QACxK,2CAA2C;QAC3C,iDAAiD;QACjD,gFAAgF;KACjF,CACF,CAAC;AACJ,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * Subject framing: how much to crop (keeping the aspect ratio) so a subject fills a chosen share\n * of the frame height.\n * @see docs/PROCESSED_OUTPUT.md for equations, coordinate/unit conventions, blockers and support\n * limits.\n */\n\nimport {\n  calculatedResult,\n  type CalculationResult\n} from \"../core/calculation-result.js\";\nimport {\n  InvalidScientificInputError,\n  requirePositiveFinite,\n  requirePositiveInteger\n} from \"../core/validation.js\";\nimport { floorWholePixels } from \"./crop.js\";\nimport { requireInputKeys } from \"../core/record-validation.js\";\n\nexport interface CalculateSubjectFramingCropInput {\n  /** Source image width in pixels. */\n  pixelWidth: number;\n  /** Source image height in pixels. */\n  pixelHeight: number;\n  /** Projected subject height in source-image pixels. */\n  subjectHeightPixels: number;\n  /** Desired subject height as a fraction of output frame height, >0 and <=1. */\n  targetSubjectHeightFraction: number;\n}\n\nexport interface SubjectFramingCrop {\n  cropFactor: number;\n  pixelWidth: number;\n  pixelHeight: number;\n  megapixels: number;\n  /**\n   * Actual subject height divided by cropped frame height.\n   * Values greater than 1 mean the subject is taller than the output frame.\n   */\n  subjectHeightFraction: number;\n  /** Whether the projected subject exceeds the cropped frame height. */\n  subjectClipped: boolean;\n  /** Whether achieving the target required an additional crop. */\n  cropped: boolean;\n}\n\n/**\n * Calculates the crop (keeping the aspect ratio) that makes a subject fill a\n * requested fraction of the frame height.\n *\n * It assumes the crop can be centred or moved freely, with enough room around\n * the subject. It does not check the subject's position against the image edges.\n */\nexport function calculateSubjectFramingCrop(\n  input: CalculateSubjectFramingCropInput\n): CalculationResult<SubjectFramingCrop> {\n  requireInputKeys(input, \"calculateSubjectFramingCrop\", [\"pixelWidth\", \"pixelHeight\", \"subjectHeightPixels\", \"targetSubjectHeightFraction\"]);\n  requirePositiveInteger(\"pixelWidth\", input.pixelWidth);\n  requirePositiveInteger(\"pixelHeight\", input.pixelHeight);\n  requirePositiveFinite(\"subjectHeightPixels\", input.subjectHeightPixels);\n  requirePositiveFinite(\n    \"targetSubjectHeightFraction\",\n    input.targetSubjectHeightFraction\n  );\n\n  if (input.targetSubjectHeightFraction > 1) {\n    throw new InvalidScientificInputError(\n      \"targetSubjectHeightFraction must be less than or equal to 1.\"\n    );\n  }\n\n  const requestedCropHeight =\n    input.subjectHeightPixels / input.targetSubjectHeightFraction;\n  const cropFactor =\n    requestedCropHeight >= input.pixelHeight\n      ? 1\n      : input.pixelHeight / requestedCropHeight;\n\n  // Height is the limiting axis. Compute it straight from the requested crop\n  // height, not through H / (H / requested), whose rounding could lose a pixel.\n  // It is the largest whole-pixel height within the requested crop, so the\n  // subject fills at least the target fraction. It is never smaller than the\n  // subject, so rounding alone never clips the subject. Width follows from that\n  // integer height and the source aspect ratio.\n  const pixelHeight =\n    cropFactor === 1\n      ? input.pixelHeight\n      : Math.min(\n          input.pixelHeight,\n          Math.max(\n            floorWholePixels(requestedCropHeight),\n            Math.ceil(input.subjectHeightPixels)\n          )\n        );\n  const pixelWidth =\n    cropFactor === 1\n      ? input.pixelWidth\n      : Math.max(\n          1,\n          Math.min(\n            input.pixelWidth,\n            Math.round((pixelHeight * input.pixelWidth) / input.pixelHeight)\n          )\n        );\n  const subjectHeightFraction =\n    input.subjectHeightPixels / pixelHeight;\n\n  return calculatedResult(\n    {\n      cropFactor,\n      pixelWidth,\n      pixelHeight,\n      megapixels: (pixelWidth * pixelHeight) / 1_000_000,\n      subjectHeightFraction,\n      subjectClipped: subjectHeightFraction > 1,\n      cropped: cropFactor > 1\n    },\n    \"subject-height-framing-crop\",\n    \"2.0.0\",\n    [\n      \"Crop height is the largest whole-pixel height within the requested crop, but never smaller than the subject; width is the nearest integer to the original aspect ratio\",\n      \"Crop can be positioned around the subject\",\n      \"Subject edge/margin constraints are not modeled\",\n      \"subjectHeightFraction may exceed 1 when the subject is clipped by frame height\"\n    ]\n  );\n}\n"]}

package/dist/processing/demosaic.d.ts

@@ -0,0 +1,90 @@
+import { type ColorMatrix3x3 } from "../color/colorimetry-math.js";
+import { type SensorColorSamplingProfile } from "../sensor/color-sampling.js";
+import { type RawReconstructionBorderHandling } from "../sensor/raw-reconstruction.js";
+import { type RawFrameReconstructionPhaseProfile } from "../capture/raw-frame-reconstruction.js";
+import { type ProcessingProvenance, type ProcessingProvider, type ProcessingRgbImage, type ProcessingTask } from "./processing-task.js";
+/** A single-channel linear CFA mosaic in absolute native site coordinates (row-major, +Y down). */
+export interface CfaMosaicPlane {
+    width: number;
+    height: number;
+    captureModeId: string;
+    /** Periodic mosaic anchored at site (0, 0) of this plane. */
+    colorSamplingProfile: SensorColorSamplingProfile;
+    /** Linear values, e.g. black-subtracted normalized RAW or merged HDR radiance. */
+    samples: Float32Array;
+}
+export interface CfaPlaneReconstructionInput {
+    plane: CfaMosaicPlane;
+    phaseProfiles: readonly RawFrameReconstructionPhaseProfile[];
+    /** Default `reject`: every neighbourhood must lie inside the plane. */
+    borderHandling?: RawReconstructionBorderHandling;
+    /** Output region in plane coordinates; default the whole plane. */
+    region?: {
+        x: number;
+        y: number;
+        width: number;
+        height: number;
+    };
+    maximumWorkingBytes?: number;
+}
+export interface CfaPlaneReconstructionOutput {
+    width: number;
+    height: number;
+    region: {
+        x: number;
+        y: number;
+        width: number;
+        height: number;
+    };
+    channelIds: readonly string[];
+    /** Interleaved `channelIds.length` channels, row-major over the region. */
+    data: Float32Array;
+    borderHandling: RawReconstructionBorderHandling;
+    provenance: ProcessingProvenance;
+}
+/**
+ * Runs explicit linear reconstruction kernels, one per CFA phase, over a typed CFA plane.
+ * - With `mirror` borders, a site outside the plane is read from its mirrored (reflect-101) position.
+ *   This is accepted only if the mirrored site has the same CFA channel, which is always true for a
+ *   2×2 Bayer pattern.
+ * - With `reject`, out-of-plane reads throw (fail closed), like the reference resolver.
+ * Values are linear and signed. No clamping, white balance or colour transform is applied.
+ */
+export declare function createCfaPlaneReconstructionTask(input: CfaPlaneReconstructionInput, provider: ProcessingProvider): ProcessingTask<CfaPlaneReconstructionOutput>;
+export interface AhdDemosaicInput {
+    plane: CfaMosaicPlane;
+    /** CFA channel ids of the red, green and blue sites; output order is red, green, blue. */
+    channelRoles: {
+        red: string;
+        green: string;
+        blue: string;
+    };
+    /**
+     * Camera-linear RGB → CIE XYZ, used only for the homogeneity metric (in CIELAB). Default: the
+     * sRGB/D65 matrix, so the plane is treated as linear sRGB for the metric. Output values are never
+     * transformed.
+     */
+    cameraToXyz?: ColorMatrix3x3;
+    maximumWorkingBytes?: number;
+}
+export interface AhdDemosaicOutput extends ProcessingRgbImage {
+    /** Per pixel: 0 horizontal, 1 vertical, 2 equal homogeneity (average of both). */
+    direction: Uint8Array;
+    provenance: ProcessingProvenance;
+}
+/**
+ * Adaptive homogeneity-directed (AHD) demosaicing after Hirakawa & Parks, IEEE TIP 14(3), 2005:
+ * 1. green is interpolated separately along rows and columns with the 5-tap directional filter
+ *    [−1/4, 1/2, 1/2, 1/2, −1/4] (bilinear green plus the own-channel second difference);
+ * 2. red and blue are filled from bilinearly interpolated colour differences (R−G, B−G) against each
+ *    directional green, giving two full-colour candidates;
+ * 3. both candidates are converted to CIELAB. At each pixel, the homogeneity is the number of
+ *    4-neighbours within adaptive luminance and chrominance tolerances ε_L, ε_C. Each tolerance is the
+ *    smaller of the two directions' largest differences along their own direction;
+ * 4. homogeneity is summed over 3×3, and the more homogeneous candidate is chosen (a tie averages
+ *    both).
+ * The paper's optional iterative median step for reducing artifacts is not applied. Mirrored
+ * (reflect-101) borders keep each site's Bayer channel. The method is non-linear: the output is not a
+ * linear function of the mosaic.
+ */
+export declare function createAhdDemosaicTask(input: AhdDemosaicInput, provider: ProcessingProvider): ProcessingTask<AhdDemosaicOutput>;