@photivra/engine 2.21.0 → 2.22.0

package/dist/capture/jpeg-xl-tiles.js

@@ -0,0 +1,151 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * #598: JPEG XL (DNG Compression 52546, DNG 1.7.0.0) strips and tiles. The engine has no JPEG XL decoder and no runtime
+ * dependency: the caller passes one (`importDngRawImage`'s `jpegXlDecoder`), which returns each tile as a PNG. The owner chose the
+ * prebuilt jxl-oxide WebAssembly package (npm `jxl-oxide-wasm`, MIT OR Apache-2.0); `createJxlOxideDngDecoder` adapts it without
+ * importing it, so the decoder is loaded only by an application that opens such a file (THIRD_PARTY.md, "JPEG XL").
+ *
+ * The PNG is read here (W3C PNG Specification, 3rd edition: IHDR, IDAT, IEND; zlib per RFC 1950/1951 through the platform's
+ * `DecompressionStream("deflate")`; the five filter types of §9). It must be grey or RGB with no alpha, the tile's size, and 16-bit
+ * for a tile of more than 8 bits. jxl-oxide writes a sample of bit depth b at 16 bits as round(v·65535/(2^b − 1)), so the stored
+ * value is recovered exactly as round(p·(2^b − 1)/65535) (each step of 65535/(2^b − 1) ≥ 1 maps to a distinct value; measured
+ * bit-exact for 8-, 12- and 16-bit lossless tiles in test/dng-import-phone.test.ts).
+ */
+import { InvalidConfigurationError } from "../core/configuration-error.js";
+/**
+ * A {@link DngJpegXlDecoder} over jxl-oxide-wasm: `forceSrgb` off (the samples stay in the codestream's encoding), one render per
+ * tile, and the image object freed after it. Usage: `import init, * as jxl from "jxl-oxide-wasm"; await init();
+ * createJxlOxideDngDecoder(jxl)`.
+ */
+export function createJxlOxideDngDecoder(module) {
+    if (module === null || typeof module !== "object" || typeof module.JxlImage !== "function")
+        throw new InvalidConfigurationError("createJxlOxideDngDecoder: module must be the initialized jxl-oxide-wasm module.");
+    return (codestream) => {
+        const image = new module.JxlImage();
+        try {
+            image.forceSrgb = false;
+            image.feedBytes(codestream);
+            if (!image.tryInit())
+                throw new InvalidConfigurationError("DNG import: the JPEG XL tile is truncated.");
+            // encodeToPng takes the render result by value (it is freed by the call; a later free() fails).
+            return image.render().encodeToPng();
+        }
+        finally {
+            image.free();
+        }
+    };
+}
+function fail(message) { throw new InvalidConfigurationError("DNG import: " + message); }
+/** @internal The decoded tile's samples at the DNG's bit depth, interleaved (pixel, channel). */
+export async function readJpegXlTilePng(png, width, height, channels, bits) {
+    if (!(png instanceof Uint8Array) || png.length < 33 || png[0] !== 0x89 || png[1] !== 0x50 || png[2] !== 0x4e || png[3] !== 0x47)
+        fail("the JPEG XL decoder did not return a PNG.");
+    const view = new DataView(png.buffer, png.byteOffset, png.byteLength);
+    let p = 8, depth = 0, colorType = -1, w = 0, h = 0, idatBytes = 0;
+    const idat = [];
+    while (p + 12 <= png.length) {
+        const length = view.getUint32(p), type = String.fromCharCode(png[p + 4], png[p + 5], png[p + 6], png[p + 7]);
+        if (p + 12 + length > png.length)
+            fail("the decoded PNG is truncated.");
+        if (type === "IHDR") {
+            w = view.getUint32(p + 8);
+            h = view.getUint32(p + 12);
+            depth = png[p + 16];
+            colorType = png[p + 17];
+            if (png[p + 20] !== 0)
+                fail("the decoded PNG is interlaced.");
+        }
+        else if (type === "IDAT") {
+            idat.push(png.subarray(p + 8, p + 8 + length));
+            idatBytes += length;
+        }
+        else if (type === "IEND")
+            break;
+        p += 12 + length;
+    }
+    const pngChannels = colorType === 0 ? 1 : colorType === 2 ? 3 : 0;
+    if (!pngChannels)
+        fail("the decoded JPEG XL tile has an alpha or palette channel; DNG tiles are grey or RGB.");
+    if (w !== width || h !== height || pngChannels !== channels)
+        fail("the decoded JPEG XL tile does not match the strip or tile size and SamplesPerPixel.");
+    if (depth !== 16 && !(depth === 8 && bits === 8))
+        fail("the JPEG XL decoder returned " + depth + "-bit samples for a " + bits + "-bit tile.");
+    const joined = new Uint8Array(idatBytes);
+    let at = 0;
+    for (const part of idat) {
+        joined.set(part, at);
+        at += part.length;
+    }
+    const bpp = channels * depth / 8, stride = 1 + width * bpp, expected = stride * height;
+    const raw = await inflate(joined, expected);
+    const out = new Uint16Array(width * height * channels), scale = (2 ** bits - 1) / 65535;
+    let previous = new Uint8Array(stride - 1), row = new Uint8Array(stride - 1);
+    for (let y = 0; y < height; y++) {
+        const filter = raw[y * stride];
+        row.set(raw.subarray(y * stride + 1, (y + 1) * stride));
+        for (let i = 0; i < row.length; i++) {
+            const a = i >= bpp ? row[i - bpp] : 0, b = previous[i], c = i >= bpp ? previous[i - bpp] : 0;
+            let predictor = 0;
+            switch (filter) {
+                case 0: break;
+                case 1:
+                    predictor = a;
+                    break;
+                case 2:
+                    predictor = b;
+                    break;
+                case 3:
+                    predictor = (a + b) >> 1;
+                    break;
+                case 4: {
+                    const pa = Math.abs(b - c), pb = Math.abs(a - c), pc = Math.abs(a + b - 2 * c);
+                    predictor = pa <= pb && pa <= pc ? a : pb <= pc ? b : c;
+                    break;
+                }
+                default: fail("the decoded PNG has an invalid filter type.");
+            }
+            row[i] = (row[i] + predictor) & 255;
+        }
+        const o = y * width * channels;
+        if (depth === 8)
+            for (let i = 0; i < width * channels; i++)
+                out[o + i] = row[i];
+        else
+            for (let i = 0; i < width * channels; i++) {
+                const v = (row[2 * i] << 8) | row[2 * i + 1];
+                out[o + i] = bits === 16 ? v : Math.round(v * scale);
+            }
+        [previous, row] = [row, previous];
+    }
+    return out;
+}
+/** zlib inflate through the platform's DecompressionStream; the output length must be exactly `expected`. */
+async function inflate(data, expected) {
+    if (typeof DecompressionStream !== "function")
+        fail("this platform has no DecompressionStream.");
+    const out = new Uint8Array(expected);
+    let at = 0;
+    try {
+        const reader = new Blob([data]).stream().pipeThrough(new DecompressionStream("deflate")).getReader();
+        for (;;) {
+            const { done, value } = await reader.read();
+            if (done)
+                break;
+            if (at + value.length > expected) {
+                await reader.cancel();
+                fail("the decoded PNG holds more data than its size.");
+            }
+            out.set(value, at);
+            at += value.length;
+        }
+    }
+    catch (error) {
+        if (error instanceof InvalidConfigurationError)
+            throw error;
+        fail("the decoded PNG's image data is corrupt.");
+    }
+    if (at !== expected)
+        fail("the decoded PNG holds less data than its size.");
+    return out;
+}
+//# sourceMappingURL=jpeg-xl-tiles.js.map

package/dist/capture/jpeg-xl-tiles.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"jpeg-xl-tiles.js","sourceRoot":"","sources":["../../src/capture/jpeg-xl-tiles.ts"],"names":[],"mappings":"AAAA,sCAAsC;AAEtC;;;;;;;;;;;GAWG;AAEH,OAAO,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAU3E;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,MAAsB;IAC7D,IAAI,MAAM,KAAK,IAAI,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,OAAO,MAAM,CAAC,QAAQ,KAAK,UAAU;QAAE,MAAM,IAAI,yBAAyB,CAAC,iFAAiF,CAAC,CAAC;IACnN,OAAO,CAAC,UAAU,EAAE,EAAE;QACpB,MAAM,KAAK,GAAG,IAAI,MAAM,CAAC,QAAQ,EAAE,CAAC;QACpC,IAAI,CAAC;YACH,KAAK,CAAC,SAAS,GAAG,KAAK,CAAC;YACxB,KAAK,CAAC,SAAS,CAAC,UAAU,CAAC,CAAC;YAC5B,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE;gBAAE,MAAM,IAAI,yBAAyB,CAAC,4CAA4C,CAAC,CAAC;YACxG,gGAAgG;YAChG,OAAO,KAAK,CAAC,MAAM,EAAE,CAAC,WAAW,EAAE,CAAC;QACtC,CAAC;gBAAS,CAAC;YAAC,KAAK,CAAC,IAAI,EAAE,CAAC;QAAC,CAAC;IAC7B,CAAC,CAAC;AACJ,CAAC;AAED,SAAS,IAAI,CAAC,OAAe,IAAW,MAAM,IAAI,yBAAyB,CAAC,cAAc,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;AAExG,iGAAiG;AACjG,MAAM,CAAC,KAAK,UAAU,iBAAiB,CAAC,GAAe,EAAE,KAAa,EAAE,MAAc,EAAE,QAAgB,EAAE,IAAY;IACpH,IAAI,CAAC,CAAC,GAAG,YAAY,UAAU,CAAC,IAAI,GAAG,CAAC,MAAM,GAAG,EAAE,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,KAAK,IAAI;QAAE,IAAI,CAAC,2CAA2C,CAAC,CAAC;IACnL,MAAM,IAAI,GAAG,IAAI,QAAQ,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,UAAU,EAAE,GAAG,CAAC,UAAU,CAAC,CAAC;IACtE,IAAI,CAAC,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,SAAS,GAAG,CAAC,CAAC;IAClE,MAAM,IAAI,GAAiB,EAAE,CAAC;IAC9B,OAAO,CAAC,GAAG,EAAE,IAAI,GAAG,CAAC,MAAM,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,IAAI,GAAG,MAAM,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,CAAC,CAAC;QACjH,IAAI,CAAC,GAAG,EAAE,GAAG,MAAM,GAAG,GAAG,CAAC,MAAM;YAAE,IAAI,CAAC,+BAA+B,CAAC,CAAC;QACxE,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;YACpB,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;YAAC,KAAK,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAE,CAAC;YAAC,SAAS,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAE,CAAC;YACtG,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC;gBAAE,IAAI,CAAC,gCAAgC,CAAC,CAAC;QAChE,CAAC;aAAM,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;YAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;YAAC,SAAS,IAAI,MAAM,CAAC;QAAC,CAAC;aAC/F,IAAI,IAAI,KAAK,MAAM;YAAE,MAAM;QAChC,CAAC,IAAI,EAAE,GAAG,MAAM,CAAC;IACnB,CAAC;IACD,MAAM,WAAW,GAAG,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,IAAI,CAAC,WAAW;QAAE,IAAI,CAAC,sFAAsF,CAAC,CAAC;IAC/G,IAAI,CAAC,KAAK,KAAK,IAAI,CAAC,KAAK,MAAM,IAAI,WAAW,KAAK,QAAQ;QAAE,IAAI,CAAC,qFAAqF,CAAC,CAAC;IACzJ,IAAI,KAAK,KAAK,EAAE,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC,CAAC;QAAE,IAAI,CAAC,+BAA+B,GAAG,KAAK,GAAG,qBAAqB,GAAG,IAAI,GAAG,YAAY,CAAC,CAAC;IAC9I,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,SAAS,CAAC,CAAC;IACzC,IAAI,EAAE,GAAG,CAAC,CAAC;IAAC,KAAK,MAAM,IAAI,IAAI,IAAI,EAAE,CAAC;QAAC,MAAM,CAAC,GAAG,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC;QAAC,EAAE,IAAI,IAAI,CAAC,MAAM,CAAC;IAAC,CAAC;IACjF,MAAM,GAAG,GAAG,QAAQ,GAAG,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,GAAG,KAAK,GAAG,GAAG,EAAE,QAAQ,GAAG,MAAM,GAAG,MAAM,CAAC;IACvF,MAAM,GAAG,GAAG,MAAM,OAAO,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC;IAC5C,MAAM,GAAG,GAAG,IAAI,WAAW,CAAC,KAAK,GAAG,MAAM,GAAG,QAAQ,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,IAAI,IAAI,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;IACxF,IAAI,QAAQ,GAAG,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IAC5E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAE,CAAC;QAChC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,MAAM,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAE,EAAE,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,GAAG,CAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAChG,IAAI,SAAS,GAAG,CAAC,CAAC;YAClB,QAAQ,MAAM,EAAE,CAAC;gBACf,KAAK,CAAC,CAAC,CAAC,MAAM;gBACd,KAAK,CAAC;oBAAE,SAAS,GAAG,CAAC,CAAC;oBAAC,MAAM;gBAC7B,KAAK,CAAC;oBAAE,SAAS,GAAG,CAAC,CAAC;oBAAC,MAAM;gBAC7B,KAAK,CAAC;oBAAE,SAAS,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;oBAAC,MAAM;gBACxC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAAC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;oBAAC,SAAS,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAAC,MAAM;gBAAC,CAAC;gBAC3J,OAAO,CAAC,CAAC,IAAI,CAAC,6CAA6C,CAAC,CAAC;YAC/D,CAAC;YACD,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAE,GAAG,SAAS,CAAC,GAAG,GAAG,CAAC;QACvC,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,QAAQ,CAAC;QAC/B,IAAI,KAAK,KAAK,CAAC;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,QAAQ,EAAE,CAAC,EAAE;gBAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,CAAC;;YAC5E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/C,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAE,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAE,CAAC;gBAC/C,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;YACvD,CAAC;QACD,CAAC,QAAQ,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC;IACpC,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED,6GAA6G;AAC7G,KAAK,UAAU,OAAO,CAAC,IAAgB,EAAE,QAAgB;IACvD,IAAI,OAAO,mBAAmB,KAAK,UAAU;QAAE,IAAI,CAAC,2CAA2C,CAAC,CAAC;IACjG,MAAM,GAAG,GAAG,IAAI,UAAU,CAAC,QAAQ,CAAC,CAAC;IACrC,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,CAAC;QACH,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,CAAC,IAA+B,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,WAAW,CAAC,IAAI,mBAAmB,CAAC,SAAS,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QAChI,SAAS,CAAC;YACR,MAAM,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,MAAM,MAAM,CAAC,IAAI,EAAE,CAAC;YAC5C,IAAI,IAAI;gBAAE,MAAM;YAChB,IAAI,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,QAAQ,EAAE,CAAC;gBAAC,MAAM,MAAM,CAAC,MAAM,EAAE,CAAC;gBAAC,IAAI,CAAC,gDAAgD,CAAC,CAAC;YAAC,CAAC;YACpH,GAAG,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;YAAC,EAAE,IAAI,KAAK,CAAC,MAAM,CAAC;QACzC,CAAC;IACH,CAAC;IAAC,OAAO,KAAK,EAAE,CAAC;QACf,IAAI,KAAK,YAAY,yBAAyB;YAAE,MAAM,KAAK,CAAC;QAC5D,IAAI,CAAC,0CAA0C,CAAC,CAAC;IACnD,CAAC;IACD,IAAI,EAAE,KAAK,QAAQ;QAAE,IAAI,CAAC,gDAAgD,CAAC,CAAC;IAC5E,OAAO,GAAG,CAAC;AACb,CAAC","sourcesContent":["// SPDX-License-Identifier: Apache-2.0\n\n/**\n * #598: JPEG XL (DNG Compression 52546, DNG 1.7.0.0) strips and tiles. The engine has no JPEG XL decoder and no runtime\n * dependency: the caller passes one (`importDngRawImage`'s `jpegXlDecoder`), which returns each tile as a PNG. The owner chose the\n * prebuilt jxl-oxide WebAssembly package (npm `jxl-oxide-wasm`, MIT OR Apache-2.0); `createJxlOxideDngDecoder` adapts it without\n * importing it, so the decoder is loaded only by an application that opens such a file (THIRD_PARTY.md, \"JPEG XL\").\n *\n * The PNG is read here (W3C PNG Specification, 3rd edition: IHDR, IDAT, IEND; zlib per RFC 1950/1951 through the platform's\n * `DecompressionStream(\"deflate\")`; the five filter types of §9). It must be grey or RGB with no alpha, the tile's size, and 16-bit\n * for a tile of more than 8 bits. jxl-oxide writes a sample of bit depth b at 16 bits as round(v·65535/(2^b − 1)), so the stored\n * value is recovered exactly as round(p·(2^b − 1)/65535) (each step of 65535/(2^b − 1) ≥ 1 maps to a distinct value; measured\n * bit-exact for 8-, 12- and 16-bit lossless tiles in test/dng-import-phone.test.ts).\n */\n\nimport { InvalidConfigurationError } from \"../core/configuration-error.js\";\n\n/** Decodes one JPEG XL codestream (a DNG strip or tile) to PNG bytes in the codestream's own colour encoding (no colour conversion). */\nexport type DngJpegXlDecoder = (codestream: Uint8Array) => Uint8Array | Promise<Uint8Array>;\n\n/** The part of the jxl-oxide-wasm module the adapter uses (its `JxlImage` class, after the module's `init`). */\nexport interface JxlOxideModule {\n  JxlImage: new () => { feedBytes(bytes: Uint8Array): unknown; tryInit(): boolean; render(): { encodeToPng(): Uint8Array }; forceSrgb: boolean; free(): void };\n}\n\n/**\n * A {@link DngJpegXlDecoder} over jxl-oxide-wasm: `forceSrgb` off (the samples stay in the codestream's encoding), one render per\n * tile, and the image object freed after it. Usage: `import init, * as jxl from \"jxl-oxide-wasm\"; await init();\n * createJxlOxideDngDecoder(jxl)`.\n */\nexport function createJxlOxideDngDecoder(module: JxlOxideModule): DngJpegXlDecoder {\n  if (module === null || typeof module !== \"object\" || typeof module.JxlImage !== \"function\") throw new InvalidConfigurationError(\"createJxlOxideDngDecoder: module must be the initialized jxl-oxide-wasm module.\");\n  return (codestream) => {\n    const image = new module.JxlImage();\n    try {\n      image.forceSrgb = false;\n      image.feedBytes(codestream);\n      if (!image.tryInit()) throw new InvalidConfigurationError(\"DNG import: the JPEG XL tile is truncated.\");\n      // encodeToPng takes the render result by value (it is freed by the call; a later free() fails).\n      return image.render().encodeToPng();\n    } finally { image.free(); }\n  };\n}\n\nfunction fail(message: string): never { throw new InvalidConfigurationError(\"DNG import: \" + message); }\n\n/** @internal The decoded tile's samples at the DNG's bit depth, interleaved (pixel, channel). */\nexport async function readJpegXlTilePng(png: Uint8Array, width: number, height: number, channels: number, bits: number): Promise<Uint16Array> {\n  if (!(png instanceof Uint8Array) || png.length < 33 || png[0] !== 0x89 || png[1] !== 0x50 || png[2] !== 0x4e || png[3] !== 0x47) fail(\"the JPEG XL decoder did not return a PNG.\");\n  const view = new DataView(png.buffer, png.byteOffset, png.byteLength);\n  let p = 8, depth = 0, colorType = -1, w = 0, h = 0, idatBytes = 0;\n  const idat: Uint8Array[] = [];\n  while (p + 12 <= png.length) {\n    const length = view.getUint32(p), type = String.fromCharCode(png[p + 4]!, png[p + 5]!, png[p + 6]!, png[p + 7]!);\n    if (p + 12 + length > png.length) fail(\"the decoded PNG is truncated.\");\n    if (type === \"IHDR\") {\n      w = view.getUint32(p + 8); h = view.getUint32(p + 12); depth = png[p + 16]!; colorType = png[p + 17]!;\n      if (png[p + 20] !== 0) fail(\"the decoded PNG is interlaced.\");\n    } else if (type === \"IDAT\") { idat.push(png.subarray(p + 8, p + 8 + length)); idatBytes += length; }\n    else if (type === \"IEND\") break;\n    p += 12 + length;\n  }\n  const pngChannels = colorType === 0 ? 1 : colorType === 2 ? 3 : 0;\n  if (!pngChannels) fail(\"the decoded JPEG XL tile has an alpha or palette channel; DNG tiles are grey or RGB.\");\n  if (w !== width || h !== height || pngChannels !== channels) fail(\"the decoded JPEG XL tile does not match the strip or tile size and SamplesPerPixel.\");\n  if (depth !== 16 && !(depth === 8 && bits === 8)) fail(\"the JPEG XL decoder returned \" + depth + \"-bit samples for a \" + bits + \"-bit tile.\");\n  const joined = new Uint8Array(idatBytes);\n  let at = 0; for (const part of idat) { joined.set(part, at); at += part.length; }\n  const bpp = channels * depth / 8, stride = 1 + width * bpp, expected = stride * height;\n  const raw = await inflate(joined, expected);\n  const out = new Uint16Array(width * height * channels), scale = (2 ** bits - 1) / 65535;\n  let previous = new Uint8Array(stride - 1), row = new Uint8Array(stride - 1);\n  for (let y = 0; y < height; y++) {\n    const filter = raw[y * stride]!;\n    row.set(raw.subarray(y * stride + 1, (y + 1) * stride));\n    for (let i = 0; i < row.length; i++) {\n      const a = i >= bpp ? row[i - bpp]! : 0, b = previous[i]!, c = i >= bpp ? previous[i - bpp]! : 0;\n      let predictor = 0;\n      switch (filter) {\n        case 0: break;\n        case 1: predictor = a; break;\n        case 2: predictor = b; break;\n        case 3: predictor = (a + b) >> 1; break;\n        case 4: { const pa = Math.abs(b - c), pb = Math.abs(a - c), pc = Math.abs(a + b - 2 * c); predictor = pa <= pb && pa <= pc ? a : pb <= pc ? b : c; break; }\n        default: fail(\"the decoded PNG has an invalid filter type.\");\n      }\n      row[i] = (row[i]! + predictor) & 255;\n    }\n    const o = y * width * channels;\n    if (depth === 8) for (let i = 0; i < width * channels; i++) out[o + i] = row[i]!;\n    else for (let i = 0; i < width * channels; i++) {\n      const v = (row[2 * i]! << 8) | row[2 * i + 1]!;\n      out[o + i] = bits === 16 ? v : Math.round(v * scale);\n    }\n    [previous, row] = [row, previous];\n  }\n  return out;\n}\n\n/** zlib inflate through the platform's DecompressionStream; the output length must be exactly `expected`. */\nasync function inflate(data: Uint8Array, expected: number): Promise<Uint8Array> {\n  if (typeof DecompressionStream !== \"function\") fail(\"this platform has no DecompressionStream.\");\n  const out = new Uint8Array(expected);\n  let at = 0;\n  try {\n    const reader = new Blob([data as Uint8Array<ArrayBuffer>]).stream().pipeThrough(new DecompressionStream(\"deflate\")).getReader();\n    for (;;) {\n      const { done, value } = await reader.read();\n      if (done) break;\n      if (at + value.length > expected) { await reader.cancel(); fail(\"the decoded PNG holds more data than its size.\"); }\n      out.set(value, at); at += value.length;\n    }\n  } catch (error) {\n    if (error instanceof InvalidConfigurationError) throw error;\n    fail(\"the decoded PNG's image data is corrupt.\");\n  }\n  if (at !== expected) fail(\"the decoded PNG holds less data than its size.\");\n  return out;\n}\n"]}

package/dist/capture/lossy-jpeg.d.ts

@@ -0,0 +1,17 @@
+/** Decoded samples, interleaved by component within each line (line, x, component), 0..255. */
+export interface LossyJpegImage {
+    width: number;
+    height: number;
+    components: number;
+    /** "ycbcr-to-rgb" when three YCbCr components were converted to RGB; "none" otherwise. */
+    colorTransform: "ycbcr-to-rgb" | "none";
+    samples: Uint8Array;
+}
+export interface LossyJpegLimits {
+    /** Largest width x height x components accepted, checked before any sample buffer is allocated. */
+    maximumSamples: number;
+    /** An abbreviated table-specification stream (SOI, DQT/DHT, EOI) read before `bytes`, as DNG's JPEGTables tag (347) supplies. */
+    tables?: Uint8Array;
+}
+/** Decodes one DCT JPEG stream (see the module notes). Every length, count and table index is checked against the buffer. */
+export declare function parseLossyJpeg(bytes: Uint8Array, limits: LossyJpegLimits): LossyJpegImage;

package/dist/capture/lossy-jpeg.js

@@ -0,0 +1,369 @@
+// SPDX-License-Identifier: Apache-2.0
+/**
+ * Size-limited decoder for 8-bit sequential DCT JPEG (ITU-T T.81 | ISO/IEC 10918-1, processes 1 and 2: SOF0 baseline and SOF1
+ * extended sequential with Huffman coding), the variant DNG uses for Compression = 34892 ("lossy JPEG", DNG 1.4.0.0) strips and
+ * tiles. Written independently from the published standard, with no third-party decoder source. Browser-safe: plain typed arrays.
+ *
+ * Decoding follows T.81 Annex F (Huffman decoding F.2.2, DC prediction, run-length AC decoding, restart intervals F.2.2.5) and the
+ * inverse DCT of A.3.3, evaluated directly in double precision with the correctly rounded cosine table `JPEG_DCT_COSINES`, then
+ * rounded half away from zero, level-shifted by 128 and clamped to 0..255 (A.3.1). DCT decoders may differ by one code from each
+ * other (T.81 leaves the IDCT accuracy to ITU-T T.83 compliance limits); this one is deterministic.
+ *
+ * Subsampled components are interpolated linearly to full size with centred sample siting (T.871, JFIF: chroma samples sit at the
+ * centre of the luma samples they cover), so a 2:1 factor uses the weights 3/4 and 1/4. Three-component data is converted from
+ * YCbCr to RGB with the T.871 equations unless the stream says it is RGB: an APP14 "Adobe" segment with transform 0, or, without
+ * JFIF or APP14 segments, component identifiers "R", "G", "B" (the conventions ITU-T T.872 describes). Four components are
+ * returned as stored.
+ *
+ * Supported: precision 8, 1..4 components with sampling factors 1..4, interleaved and non-interleaved sequential scans, DRI restart
+ * intervals, 8- and 16-bit quantization tables, tables defined in a separate abbreviated stream (DNG JPEGTables). Rejected:
+ * progressive, hierarchical, lossless and arithmetic-coded frames, 12-bit precision, truncated data.
+ */
+import { InvalidConfigurationError } from "../core/configuration-error.js";
+import { JPEG_DCT_COSINES } from "./export-jpeg.js";
+function fail(message) { throw new InvalidConfigurationError("Lossy JPEG: " + message); }
+const ZIGZAG = /* @__PURE__ */ (() => {
+    const z = new Int32Array(64);
+    for (let d = 0, i = 0; d <= 14; d++) {
+        const points = [];
+        for (let y = 0; y < 8; y++) {
+            const x = d - y;
+            if (x >= 0 && x < 8)
+                points.push(y * 8 + x);
+        }
+        for (const p of d % 2 === 0 ? points.reverse() : points)
+            z[i++] = p;
+    }
+    return z;
+})();
+/** T.81 Annex C code generation and F.2.2.3 decoder tables. */
+function buildHuffman(counts, values) {
+    const maxCode = new Int32Array(18).fill(-1), minCode = new Int32Array(17), valPtr = new Int32Array(17);
+    let code = 0, k = 0;
+    for (let length = 1; length <= 16; length++) {
+        const n = counts[length - 1];
+        if (n) {
+            valPtr[length] = k;
+            minCode[length] = code;
+            code += n;
+            k += n;
+            maxCode[length] = code - 1;
+        }
+        if (code > (1 << length))
+            fail("over-subscribed Huffman table.");
+        code <<= 1;
+    }
+    maxCode[17] = 0x7fffffff;
+    return { maxCode, minCode, valPtr, values };
+}
+/** Decodes one DCT JPEG stream (see the module notes). Every length, count and table index is checked against the buffer. */
+export function parseLossyJpeg(bytes, limits) {
+    if (!(bytes instanceof Uint8Array))
+        fail("input must be a Uint8Array.");
+    if (limits === null || typeof limits !== "object" || !Number.isSafeInteger(limits.maximumSamples) || limits.maximumSamples < 1)
+        fail("invalid sample limit.");
+    if (limits.tables !== undefined && !(limits.tables instanceof Uint8Array))
+        fail("tables must be a Uint8Array.");
+    const quant = [], dcTables = [], acTables = [];
+    const state = { restart: 0, jfif: false, adobe: -1 };
+    let scanned = false;
+    let frame;
+    const segments = (data, tablesOnly) => {
+        const n = data.length;
+        if (n < 4 || data[0] !== 0xff || data[1] !== 0xd8)
+            fail("missing SOI marker.");
+        let p = 2;
+        for (;;) {
+            while (p < n && data[p] === 0xff && data[p + 1] === 0xff)
+                p++;
+            if (p + 2 > n || data[p] !== 0xff)
+                fail("malformed marker segment.");
+            const marker = data[p + 1];
+            if (marker === 0xd9) {
+                if (tablesOnly || scanned)
+                    return;
+                fail("no scan before EOI.");
+            }
+            if (p + 4 > n)
+                fail("malformed marker segment.");
+            const size = data[p + 2] * 256 + data[p + 3], start = p + 4, end = p + 2 + size;
+            if (size < 2 || end > n)
+                fail("segment exceeds the stream.");
+            if (marker === 0xdb) {
+                for (let q = start; q < end;) {
+                    const pq = data[q] >> 4, tq = data[q] & 15, len = pq ? 128 : 64;
+                    if (pq > 1 || tq > 3 || q + 1 + len > end)
+                        fail("invalid quantization table.");
+                    const table = new Float64Array(64);
+                    for (let i = 0; i < 64; i++) {
+                        const v = pq ? data[q + 1 + i * 2] * 256 + data[q + 2 + i * 2] : data[q + 1 + i];
+                        if (v === 0)
+                            fail("zero quantization value.");
+                        table[ZIGZAG[i]] = v;
+                    }
+                    quant[tq] = table;
+                    q += 1 + len;
+                }
+            }
+            else if (marker === 0xc4) {
+                for (let q = start; q < end;) {
+                    const tc = data[q] >> 4, th = data[q] & 15;
+                    if (tc > 1 || th > 3 || q + 17 > end)
+                        fail("invalid Huffman table.");
+                    const counts = data.subarray(q + 1, q + 17);
+                    let total = 0;
+                    for (const c of counts)
+                        total += c;
+                    if (total > 256 || q + 17 + total > end)
+                        fail("Huffman table exceeds its segment.");
+                    (tc ? acTables : dcTables)[th] = buildHuffman(counts, data.slice(q + 17, q + 17 + total));
+                    q += 17 + total;
+                }
+            }
+            else if (marker === 0xdd) {
+                if (size !== 4)
+                    fail("invalid restart interval.");
+                state.restart = data[start] * 256 + data[start + 1];
+            }
+            else if (marker === 0xe0) {
+                if (size >= 7 && data[start] === 0x4a && data[start + 1] === 0x46 && data[start + 2] === 0x49 && data[start + 3] === 0x46 && data[start + 4] === 0)
+                    state.jfif = true;
+            }
+            else if (marker === 0xee) {
+                if (size >= 14 && String.fromCharCode(...data.subarray(start, start + 5)) === "Adobe")
+                    state.adobe = data[start + 11];
+            }
+            else if (marker === 0xc0 || marker === 0xc1) {
+                if (tablesOnly || frame)
+                    fail("unexpected frame header.");
+                if (size < 8)
+                    fail("unsupported frame header.");
+                const precision = data[start], height = data[start + 1] * 256 + data[start + 2], width = data[start + 3] * 256 + data[start + 4], nf = data[start + 5];
+                if (precision !== 8)
+                    fail("only 8-bit precision is supported.");
+                if (height < 1 || width < 1 || nf < 1 || nf > 4 || size !== 8 + 3 * nf)
+                    fail("unsupported frame header.");
+                if (width * height * nf > limits.maximumSamples)
+                    fail("frame exceeds the sample limit.");
+                const raw = [];
+                for (let i = 0; i < nf; i++) {
+                    const o = start + 6 + i * 3, h = data[o + 1] >> 4, v = data[o + 1] & 15, tq = data[o + 2];
+                    if (h < 1 || h > 4 || v < 1 || v > 4 || tq > 3)
+                        fail("invalid component sampling or table.");
+                    if (raw.some((c) => c.id === data[o]))
+                        fail("duplicate component identifier.");
+                    raw.push({ id: data[o], h, v, tq });
+                }
+                const hMax = Math.max(...raw.map((c) => c.h)), vMax = Math.max(...raw.map((c) => c.v));
+                const mcusWide = Math.ceil(width / (8 * hMax)), mcusHigh = Math.ceil(height / (8 * vMax));
+                const components = raw.map((c) => {
+                    const blocksWide = mcusWide * c.h, blocksHigh = mcusHigh * c.v;
+                    return { ...c, blocksWide, blocksHigh, plane: new Uint8Array(blocksWide * blocksHigh * 64), previousDc: 0 };
+                });
+                frame = { width, height, components, hMax, vMax, mcusWide, mcusHigh };
+            }
+            else if (marker === 0xda) {
+                if (tablesOnly || !frame)
+                    fail("scan before frame header.");
+                const ns = data[start];
+                if (ns < 1 || ns > 4 || size !== 6 + 2 * ns)
+                    fail("invalid scan header.");
+                const scan = [];
+                for (let i = 0; i < ns; i++) {
+                    const c = frame.components.find((x) => x.id === data[start + 1 + i * 2]);
+                    if (!c || scan.includes(c))
+                        fail("scan references an unknown component.");
+                    const td = data[start + 2 + i * 2] >> 4, ta = data[start + 2 + i * 2] & 15;
+                    c.dc = dcTables[td];
+                    c.ac = acTables[ta];
+                    if (!c.dc || !c.ac)
+                        fail("scan references an undefined Huffman table.");
+                    if (!quant[c.tq])
+                        fail("component references an undefined quantization table.");
+                    scan.push(c);
+                }
+                const ss = data[start + 1 + 2 * ns], se = data[start + 2 + 2 * ns], a = data[start + 3 + 2 * ns];
+                if (ss !== 0 || se !== 63 || a !== 0)
+                    fail("only sequential scans are supported.");
+                if (scan.reduce((s, c) => s + (ns === 1 ? 1 : c.h * c.v), 0) > 10)
+                    fail("too many blocks per MCU.");
+                // Further scans may follow (non-interleaved sequential); EOI ends the image.
+                p = decodeScan(data, end, frame, scan, quant, state.restart);
+                scanned = true;
+                continue;
+            }
+            else if ((marker >= 0xc2 && marker <= 0xcf && marker !== 0xc4 && marker !== 0xc8 && marker !== 0xcc) || marker === 0xde) {
+                fail("only sequential Huffman DCT frames (SOF0, SOF1) are supported.");
+            }
+            else if (marker === 0xcc) {
+                fail("arithmetic coding is not supported.");
+            }
+            else if (marker < 0xc0 || marker === 0xd8 || (marker >= 0xd0 && marker <= 0xd7)) {
+                fail("unexpected marker.");
+            }
+            p = end;
+        }
+    };
+    if (limits.tables)
+        segments(limits.tables, true);
+    segments(bytes, false);
+    if (!frame)
+        fail("no image.");
+    if (frame.components.some((c) => !c.dc))
+        fail("a component has no scan.");
+    return assemble(frame, state);
+}
+function decodeScan(data, start, frame, scan, quant, restartInterval) {
+    let p = start, buffer = 0, bits = 0, padding = 0, marker = -1;
+    const fill = () => {
+        let byte = 0;
+        if (marker < 0 && p < data.length && data[p] !== 0xff)
+            byte = data[p++];
+        else if (marker < 0 && p + 1 < data.length && data[p] === 0xff && data[p + 1] === 0x00) {
+            byte = 0xff;
+            p += 2;
+        }
+        else {
+            if (marker < 0 && p + 1 < data.length)
+                marker = data[p + 1];
+            if (++padding > 8)
+                fail("entropy-coded data is truncated.");
+        }
+        buffer = ((buffer << 8) | byte) >>> 0;
+        bits += 8;
+    };
+    const readBits = (count) => {
+        if (count === 0)
+            return 0;
+        while (bits < count)
+            fill();
+        bits -= count;
+        return (buffer >>> bits) & ((1 << count) - 1);
+    };
+    const decode = (t) => {
+        let code = readBits(1), len = 1;
+        while (code > t.maxCode[len]) {
+            if (++len > 16)
+                fail("invalid Huffman code.");
+            code = (code << 1) | readBits(1);
+        }
+        const v = t.values[t.valPtr[len] + code - t.minCode[len]];
+        if (v === undefined)
+            fail("invalid Huffman code.");
+        return v;
+    };
+    const extend = (v, s) => v < (1 << (s - 1)) ? v - (1 << s) + 1 : v;
+    const coefficients = new Float64Array(64), rows = new Float64Array(64), cos = JPEG_DCT_COSINES;
+    const block = (c, bx, by) => {
+        coefficients.fill(0);
+        const t = decode(c.dc);
+        if (t > 11)
+            fail("invalid DC difference category.");
+        c.previousDc += t ? extend(readBits(t), t) : 0;
+        const q = quant[c.tq];
+        coefficients[0] = c.previousDc * q[0];
+        for (let k = 1; k < 64;) {
+            const rs = decode(c.ac), r = rs >> 4, s = rs & 15;
+            if (s === 0) {
+                if (r === 15) {
+                    k += 16;
+                    continue;
+                }
+                break;
+            }
+            k += r;
+            if (k > 63 || s > 10)
+                fail("invalid AC coefficient.");
+            const z = ZIGZAG[k];
+            coefficients[z] = extend(readBits(s), s) * q[z];
+            k++;
+        }
+        // A.3.3: s(y,x) = 1/4 Σ_u Σ_v C(u)C(v) S(v,u) cos((2x+1)uπ/16) cos((2y+1)vπ/16), rows (u) first, then columns (v).
+        for (let v = 0; v < 8; v++)
+            for (let x = 0; x < 8; x++) {
+                let sum = 0;
+                for (let u = 0; u < 8; u++)
+                    sum += (u === 0 ? Math.SQRT1_2 : 1) * coefficients[v * 8 + u] * cos[u * 8 + x];
+                rows[v * 8 + x] = sum;
+            }
+        const stride = c.blocksWide * 8, base = by * 8 * stride + bx * 8;
+        for (let y = 0; y < 8; y++)
+            for (let x = 0; x < 8; x++) {
+                let sum = 0;
+                for (let v = 0; v < 8; v++)
+                    sum += (v === 0 ? Math.SQRT1_2 : 1) * rows[v * 8 + x] * cos[v * 8 + y];
+                const value = sum / 4, rounded = value < 0 ? -Math.round(-value) : Math.round(value);
+                c.plane[base + y * stride + x] = Math.min(255, Math.max(0, rounded + 128));
+            }
+    };
+    for (const c of scan)
+        c.previousDc = 0;
+    // A single-component scan is non-interleaved: one block per MCU over the component's own (unpadded-to-MCU) block grid.
+    const single = scan.length === 1 ? scan[0] : undefined;
+    const unitsWide = single ? Math.ceil(Math.ceil(frame.width * single.h / frame.hMax) / 8) : frame.mcusWide;
+    const unitsHigh = single ? Math.ceil(Math.ceil(frame.height * single.v / frame.vMax) / 8) : frame.mcusHigh;
+    const total = unitsWide * unitsHigh;
+    for (let m = 0, restarts = 0; m < total; m++) {
+        if (restartInterval && m > 0 && m % restartInterval === 0) {
+            bits = 0;
+            buffer = 0;
+            padding = 0;
+            marker = -1;
+            if (p + 1 >= data.length || data[p] !== 0xff || data[p + 1] !== (0xd0 | (restarts & 7)))
+                fail("missing restart marker.");
+            p += 2;
+            restarts++;
+            for (const c of scan)
+                c.previousDc = 0;
+        }
+        const mx = m % unitsWide, my = Math.floor(m / unitsWide);
+        if (single)
+            block(single, mx, my);
+        else
+            for (const c of scan)
+                for (let v = 0; v < c.v; v++)
+                    for (let h = 0; h < c.h; h++)
+                        block(c, mx * c.h + h, my * c.v + v);
+    }
+    // Skip to the next marker: remaining fill bits are discarded (T.81 F.1.2.3 pads with 1-bits).
+    while (p < data.length && !(data[p] === 0xff && data[p + 1] !== 0x00 && !(data[p + 1] >= 0xd0 && data[p + 1] <= 0xd7)))
+        p++;
+    if (p + 1 >= data.length)
+        fail("missing EOI marker.");
+    return p;
+}
+function assemble(frame, state) {
+    const { width, height, components, hMax, vMax } = frame, nf = components.length, out = new Uint8Array(width * height * nf);
+    components.forEach((c, ci) => {
+        const stride = c.blocksWide * 8, cw = Math.ceil(width * c.h / hMax), ch = Math.ceil(height * c.v / vMax);
+        if (c.h === hMax && c.v === vMax) {
+            for (let y = 0; y < height; y++)
+                for (let x = 0; x < width; x++)
+                    out[(y * width + x) * nf + ci] = c.plane[y * stride + x];
+            return;
+        }
+        // Centred linear interpolation: output sample x covers component position (x + 0.5)·h/hMax − 0.5, clamped to the samples.
+        const sx = c.h / hMax, sy = c.v / vMax;
+        for (let y = 0; y < height; y++) {
+            const fy = Math.min(Math.max((y + 0.5) * sy - 0.5, 0), ch - 1), y0 = Math.floor(fy), y1 = Math.min(y0 + 1, ch - 1), wy = fy - y0;
+            for (let x = 0; x < width; x++) {
+                const fx = Math.min(Math.max((x + 0.5) * sx - 0.5, 0), cw - 1), x0 = Math.floor(fx), x1 = Math.min(x0 + 1, cw - 1), wx = fx - x0;
+                const a = c.plane[y0 * stride + x0], b = c.plane[y0 * stride + x1], d = c.plane[y1 * stride + x0], e = c.plane[y1 * stride + x1];
+                out[(y * width + x) * nf + ci] = Math.round((a * (1 - wx) + b * wx) * (1 - wy) + (d * (1 - wx) + e * wx) * wy);
+            }
+        }
+    });
+    const ids = components.map((c) => c.id), rgbIds = ids[0] === 0x52 && ids[1] === 0x47 && ids[2] === 0x42;
+    const ycc = nf === 3 && (state.adobe >= 0 ? state.adobe !== 0 : state.jfif || !rgbIds);
+    if (ycc)
+        for (let i = 0; i < out.length; i += 3) {
+            // T.871 §7: R = Y + 1.402 (Cr − 128), G = Y − 0.344136 (Cb − 128) − 0.714136 (Cr − 128), B = Y + 1.772 (Cb − 128).
+            const y = out[i], cb = out[i + 1] - 128, cr = out[i + 2] - 128;
+            out[i] = clamp8(y + 1.402 * cr);
+            out[i + 1] = clamp8(y - 0.344136 * cb - 0.714136 * cr);
+            out[i + 2] = clamp8(y + 1.772 * cb);
+        }
+    return { width, height, components: nf, colorTransform: ycc ? "ycbcr-to-rgb" : "none", samples: out };
+}
+function clamp8(v) { return Math.min(255, Math.max(0, Math.round(v))); }
+//# sourceMappingURL=lossy-jpeg.js.map