emf-converter 4.8.1 → 4.8.3

package/README.md

@@ -217,7 +217,7 @@ Everything is measured against output painted by Windows itself; `src/gdi-parity
 - **Colour adjustment and image effects**: all eleven EMF+ image effects are supported. Complete native tables cover every legal ColorCurve intensity. Deferred draws retain their effects. JPEG, GIF and TIFF have bundled decoders for environments without Canvas. HSL, tint, red-eye, blur, sharpen and GDI colour adjustment retain measured differences; see [the detailed limitations](docs/limitations.md) and the native fixture bounds.
 - **Pen transforms**: EMF+ pen transforms (uniform, nonuniform and skewed) match GDI+ exactly on the 17 pen-transform fixtures, including dashes, which GDI+ lays out along the path in world space at the untransformed pen width. Native fixtures verify dashes under both uniform and non-uniform scales. GDI+ rejects singular pen transforms; the converter ignores such malformed transform data. Square, round, diamond and arrow anchor caps use their native shapes; the new centered-cap fixtures match every pixel with and without antialiasing.
 - **Text**: ANSI, multi-string and small-text records now render; the Windows text fixtures retain a few glyph-edge differences and a `PolyTextOut` C1 control-glyph difference (under 0.1% of pixels; bounds in `src/emf-text-records.fixture.test.ts`). ANSI decoding uses the host's `TextDecoder` for common Windows code pages; Johab and OEM CP437 use bundled Windows mappings. Vertical `ETO_PDY` advances work with and without the `fonts` engine. Without `fonts`, SVG text measurements use deterministic estimates, so precise justification requires supplying the matching fonts.
-- **Wide pens and paths**: flat-capped GDI pens 7 px and wider can differ by a few pixels at round joins, and dashed wide Bezier curves follow `WidenPath` (which Windows' direct drawing does not quite match); at most 0.2% of pixels on the fixtures. `EMR_WIDENPATH` does not reproduce the extra inner join triangles GDI's own `WidenPath` emits (visible only when the widened outline is itself stroked). EMF+ 1-pixel antialiased lines can differ by one antialiasing sample at their ends, some closed widened outlines by one sample along an edge, and Inset or compound pens on closed figures are approximate.
+- **Wide pens and paths**: flat-capped GDI pens 7 px and wider can differ by a few pixels at round joins, and dashed wide curves follow `WidenPath` (exact for round and flat caps at one logical unit per pixel; square caps differ slightly at path ends); at most 0.2% of pixels on the fixtures. `EMR_WIDENPATH` does not reproduce the extra inner join triangles GDI's own `WidenPath` emits (visible only when the widened outline is itself stroked). EMF+ 1-pixel antialiased lines can differ by one antialiasing sample at their ends, some closed widened outlines by one sample along an edge, and Inset or compound pens on closed figures are approximate.
 - **GM_COMPATIBLE recordings**: EMF files do not record the graphics mode, and Windows plays RoundRect, Arc, Chord, Pie and null-pen Ellipse records back differently from how a GM_COMPATIBLE application drew them on screen; the converter follows Windows' playback.
 - **Small EMF+ residuals**: rotated `HighQualityBicubic` `DrawImage` edge pixels (0.14%), one-level differences at exact half-level `Blend` knots, and a few pixels of a metafile nested in `DrawImage` under a scale.
 - **WMF**: `PS_INSIDEFRAME` boxes can come out a pixel short at non-integer scales, right-to-left (`LAYOUT_RTL`) layouts differ by single pixels on mirrored diagonals, and metric map modes assume a 96 dpi reference device (Windows derives them from the physical display, so its own output varies per machine).

package/dist/browser.js

@@ -17039,13 +17039,17 @@ function expandedBlur(rgba, width, height, radius, x0, y0, x1, y1, cropped) {
   const blurred = blurReduction(radius) > 1 ? effectBlur(buf, bw, bh, radius, cropped ? { x: r, y: r, w: x1 - x0, h: y1 - y0 } : void 0, true) : gdipBlur(buf, bw, bh, radius, cropped ? r : 0, cropped, Math.ceil((y1 - y0) / (x1 - x0)));
   const w = x1 - x0;
   const out = new Uint8ClampedArray(w * (y1 - y0) * 4);
+  const halo = new Uint8ClampedArray((y1 - y0) * 4);
   for (let y = 0; y < y1 - y0; y++) {
     const from = ((y + r) * bw + r) * 4;
     for (let i = 0; i < w * 4; i++) {
       out[y * w * 4 + i] = Math.round(blurred[from + i]);
     }
+    for (let i = 0; i < 4; i++) {
+      halo[y * 4 + i] = Math.round(blurred[from + w * 4 + i]);
+    }
   }
-  return out;
+  return { out, halo };
 }
 function applyImageEffectToRect(rgba, width, height, effect, src) {
   if (width <= 0 || height <= 0 || rgba.length !== width * height * 4) {
@@ -17078,7 +17082,12 @@ function applyImageEffectToRect(rgba, width, height, effect, src) {
   };
   if (src && effect.kind === "blur" && effect.expandEdge) {
     const cropped = Math.ceil(src.x + src.w) < width && Math.ceil(src.y + src.h) < height;
-    return result(expandedBlur(rgba, width, height, effect.radius, x0, y0, x1, y1, cropped));
+    const { out: out2, halo } = expandedBlur(rgba, width, height, effect.radius, x0, y0, x1, y1, cropped);
+    const region2 = result(out2);
+    if (shownWidth === w && shownHeight === h) {
+      region2.haloRight = halo;
+    }
+    return region2;
   }
   let region = rgba;
   if (w !== width || h !== height) {
@@ -20450,6 +20459,22 @@ function cubicIntegral(t, a) {
   }
   return Math.sign(t) * v;
 }
+var FAR_FADE = [1, 1, 0.985, 0.95, 0.91, 0.85, 0.8, 0.745, 0.69, 0.63, 0.56, 0.5, 0.43, 0.35, 0.28, 0.22, 0.16, 0.11, 0.065, 0.03, 0];
+function fadeAt(d, bilinear) {
+  if (d >= 1) {
+    return 0;
+  }
+  if (bilinear) {
+    return 1 - Math.max(0, d);
+  }
+  const x = Math.max(0, d) * 20;
+  const i = Math.floor(x);
+  return FAR_FADE[i] + (FAR_FADE[i + 1] - FAR_FADE[i]) * (x - i);
+}
+function nearSourceExtent(deviceLen, srcLen) {
+  const d = Math.ceil(deviceLen - 1e-6) - srcLen;
+  return d >= -1 && d <= 1;
+}
 var BICUBIC_A = -0.5;
 var HQ_BICUBIC_A = -1;
 function axisFilter(kernel, scale, copy) {
@@ -20554,13 +20579,39 @@ function mirroredOrigin(o) {
   return f === 0 ? o : Math.floor(o) + 1 - f;
 }
 function resampleImage(rgba, width, height, spec, surface) {
+  if (spec.rightHalo && spec.rightHalo.length === height * 4) {
+    const grown = new Uint8ClampedArray((width + 1) * height * 4);
+    for (let y = 0; y < height; y++) {
+      grown.set(rgba.subarray(y * width * 4, (y + 1) * width * 4), y * (width + 1) * 4);
+      grown.set(spec.rightHalo.subarray(y * 4, y * 4 + 4), (y * (width + 1) + width) * 4);
+    }
+    return resampleImage(grown, width + 1, height, { ...spec, rightHalo: void 0 }, surface);
+  }
   if (spec.kernel === "nearest" && !spec.wrap) {
     return resampleNearest(rgba, width, height, spec, surface);
   }
   let m = snapToDeviceGrid(spec);
-  const kernel = spec.kernel;
+  let kernel = spec.kernel;
   const hq = kernel === "hq-bilinear" || kernel === "hq-bicubic";
   const axisAligned = m[1] === 0 && m[2] === 0;
+  let farFade = false;
+  let shiftX = 0;
+  let shiftY = 0;
+  if (hq && !axisAligned) {
+    const extentU = Math.hypot(m[0], m[1]) * spec.srcW;
+    const extentV = Math.hypot(m[2], m[3]) * spec.srcH;
+    if (nearSourceExtent(extentU, spec.srcW) || nearSourceExtent(extentV, spec.srcH)) {
+      kernel = kernel === "hq-bicubic" ? "bicubic" : "bilinear";
+    } else {
+      farFade = true;
+      if (spec.halfPixelOffset) {
+        const lu = Math.hypot(m[0], m[1]);
+        const lv = Math.hypot(m[2], m[3]);
+        shiftX = 0.5 * m[0] / lu + 0.5 * m[2] / lv;
+        shiftY = 0.5 * m[1] / lu + 0.5 * m[3] / lv;
+      }
+    }
+  }
   if (hq && axisAligned) {
     const ox = m[0] * spec.srcX + m[4];
     m = [m[0], 0, 0, m[3], m[4] + mirroredOrigin(ox) - ox, m[5]];
@@ -20615,6 +20666,8 @@ function resampleImage(rgba, width, height, spec, surface) {
   const vMax = spec.srcY + spec.srcH;
   const wu = [];
   const wv = [];
+  const gradU = Math.hypot(inv2[0], inv2[2]);
+  const gradV = Math.hypot(inv2[1], inv2[3]);
   const c = spec.clampArgb ?? 0;
   const edge = {
     wrap: spec.wrap,
@@ -20626,8 +20679,8 @@ function resampleImage(rgba, width, height, spec, surface) {
     const py = by0 + j + o;
     for (let i = 0; i < w; i++) {
       const px = bx0 + i + o;
-      const u = inv2[0] * px + inv2[2] * py + inv2[4];
-      const v = inv2[1] * px + inv2[3] * py + inv2[5];
+      const u = inv2[0] * (px - shiftX) + inv2[2] * (py - shiftY) + inv2[4];
+      const v = inv2[1] * (px - shiftX) + inv2[3] * (py - shiftY) + inv2[5];
       const cx = px + COVERAGE_NUDGE_X;
       const cy = py + COVERAGE_NUDGE_Y;
       const eu = inv2[0] * cx + inv2[2] * cy + inv2[4];
@@ -20635,31 +20688,35 @@ function resampleImage(rgba, width, height, spec, surface) {
       if (eu < spec.srcX || ev < spec.srcY || eu >= uMax || ev >= vMax) {
         continue;
       }
-      const cu = u;
-      const cv = v;
-      const iu0 = Math.ceil(cu - fu.radius);
-      const iu1 = Math.floor(cu + fu.radius);
-      const iv0 = Math.ceil(cv - fv.radius);
-      const iv1 = Math.floor(cv + fv.radius);
-      wu.length = 0;
-      wv.length = 0;
-      for (let t = iu0; t <= iu1; t++) {
-        wu.push(fu.weight(t, cu));
-      }
-      for (let t = iv0; t <= iv1; t++) {
-        wv.push(fv.weight(t, cv));
-      }
-      const [r, g, b, a] = blendSeparable(
-        rgba,
-        width,
-        spec.wrap ? { x0: 0, y0: 0, x1: width, y1: height } : { x0: sx0, y0: sy0, x1: sx1, y1: sy1 },
-        iu0,
-        wu,
-        iv0,
-        wv,
-        kernel === "bicubic",
-        edge
-      );
+      const box = { x0: 0, y0: 0, x1: width, y1: height };
+      const sampleAt = (cu, cv) => {
+        const iu0 = Math.ceil(cu - fu.radius);
+        const iu1 = Math.floor(cu + fu.radius);
+        const iv0 = Math.ceil(cv - fv.radius);
+        const iv1 = Math.floor(cv + fv.radius);
+        wu.length = 0;
+        wv.length = 0;
+        for (let t = iu0; t <= iu1; t++) {
+          wu.push(fu.weight(t, cu));
+        }
+        for (let t = iv0; t <= iv1; t++) {
+          wv.push(fv.weight(t, cv));
+        }
+        return blendSeparable(rgba, width, box, iu0, wu, iv0, wv, kernel === "bicubic", edge);
+      };
+      let [r, g, b, a] = sampleAt(u, v);
+      if (farFade && !spec.wrap && a > 0) {
+        const du = (uMax - u) / gradU;
+        const dv = (vMax - v) / gradV;
+        const k = 1 - (du < 1 ? fadeAt(du, kernel === "hq-bilinear") * (sampleAt(uMax, v)[3] / a) : 0) - (dv < 1 ? fadeAt(dv, kernel === "hq-bilinear") * (sampleAt(u, vMax)[3] / a) : 0);
+        if (k < 1) {
+          const f = Math.max(0, k);
+          r *= f;
+          g *= f;
+          b *= f;
+          a *= f;
+        }
+      }
       if (a <= 0) {
         continue;
       }
@@ -23486,7 +23543,8 @@ function effectedDraw(img, bitmap, effect, draw, source) {
         srcY,
         srcW,
         srcH,
-        toDevice: mulMatrix(draw.resample.toDevice, shift)
+        toDevice: mulMatrix(draw.resample.toDevice, shift),
+        rightHalo: region.haloRight
       }
     }
   };
@@ -29228,9 +29286,9 @@ function perpendicularVectors(width, dx0, dy0) {
   const vx = rounded[0], vy = rounded[1];
   return { v: flip ? [-vx, -vy] : [vx, vy], ray: flip ? [-x, -y] : [x, y] };
 }
-function squareExtension(width, dx, dy) {
+function squareExtension(width, dx, dy, scale = 1) {
   const len = Math.hypot(dx, dy);
-  const r = width / 2;
+  const r = width / 2 * scale;
   return [Math.floor(dx / len * r + 0.5), Math.floor(dy / len * r + 0.5)];
 }
 function turnSign(ax, ay, bx, by) {
@@ -29269,7 +29327,7 @@ var Outliner = class {
     const [L, R] = drawVertices(this.pen, d[0], d[1]);
     const perpendicular = curve && drawDir ? drawDir : [dx, dy];
     const vectors = perpendicularVectors(this.opts.width, perpendicular[0], perpendicular[1]);
-    return { dx, dy, L, R, v: vectors.v, vRaw: vectors.ray, e: squareExtension(this.opts.width, perpendicular[0], perpendicular[1]), curveEnd: !!drawDir, tangentLength: drawDir ? drawDir[0] * drawDir[0] + drawDir[1] * drawDir[1] : 0, curve };
+    return { dx, dy, L, R, v: vectors.v, vRaw: vectors.ray, e: squareExtension(this.opts.width, perpendicular[0], perpendicular[1]), pe: [perpendicular[0], perpendicular[1]], curveEnd: !!drawDir, tangentLength: drawDir ? drawDir[0] * drawDir[0] + drawDir[1] * drawDir[1] : 0, curve };
   }
   push(p, v) {
     this.pts.push([p[0] + v[0], p[1] + v[1]]);
@@ -29364,7 +29422,7 @@ var Outliner = class {
     return side === "R" ? s.v : [-s.v[0], -s.v[1]];
   }
   /** A cap at `p` from the `from` side of `s` round to the other side (`start`: around the back). */
-  cap(p, s, start) {
+  cap(p, s, start, scale = 1) {
     const from = start ? "L" : "R";
     const to = start ? "R" : "L";
     const { cap } = this.opts;
@@ -29382,7 +29440,8 @@ var Outliner = class {
     const sv = this.capSide(s, from);
     const ev = this.capSide(s, to);
     if (cap === "square") {
-      const e = start ? [-s.e[0], -s.e[1]] : s.e;
+      const e0 = scale === 1 ? s.e : squareExtension(this.opts.width, s.pe[0], s.pe[1], scale);
+      const e = start ? [-e0[0], -e0[1]] : e0;
       this.push(p, [sv[0] + e[0], sv[1] + e[1]]);
       this.push(p, [ev[0] + e[0], ev[1] + e[1]]);
     } else {
@@ -29478,7 +29537,7 @@ var Outliner = class {
    * lets a single point stand for a zero-length dash along `dirs[0]`;
    * `draws` (curve end tangents) picks draw vertices only.
    */
-  open(P, dirs, draws, curves) {
+  open(P, dirs, draws, curves, startScale = 1, endScale = 1) {
     if (P.length < 2 && !dirs?.[0]) {
       this.dot(P[0]);
       return;
@@ -29489,12 +29548,12 @@ var Outliner = class {
     }
     if (segs2.length === 0) {
       const s = this.seg(P[0], P[0], dirs?.[0], draws?.[0], curves?.[0]);
-      this.cap(P[0], s, true);
-      this.cap(P[0], s, false);
+      this.cap(P[0], s, true, startScale);
+      this.cap(P[0], s, false, endScale);
       this.flush();
       return;
     }
-    this.cap(P[0], segs2[0], true);
+    this.cap(P[0], segs2[0], true, startScale);
     for (let i = 0; i + 1 < segs2.length; i++) {
       const s = segs2[i];
       const t = segs2[i + 1];
@@ -29505,7 +29564,7 @@ var Outliner = class {
         this.join(P[i + 1], s, t, "R", tr < 0);
       }
     }
-    this.cap(P[P.length - 1], segs2[segs2.length - 1], false);
+    this.cap(P[P.length - 1], segs2[segs2.length - 1], false, endScale);
     for (let i = segs2.length - 1; i > 0; i--) {
       const s = segs2[i];
       const t = segs2[i - 1];
@@ -29561,7 +29620,7 @@ var Outliner = class {
     this.flush();
   }
 };
-function dashPieces(P, tangents, pattern, shorten) {
+function dashPieces(P, tangents, curveFlags, pattern, shorten, wholePixelVectors) {
   const out = [];
   const dashes = [];
   for (let i = 0; i < pattern.length; i += 2) {
@@ -29576,12 +29635,22 @@ function dashPieces(P, tangents, pattern, shorten) {
   let idx = 0;
   let left = dashes[0];
   let on = true;
-  let cur = { pts: [P[0]], dirs: [], draws: [] };
+  let cur = { pts: [P[0]], dirs: [], draws: [], curves: [], startScale: 1, endScale: 1 };
   for (let i = 0; i + 1 < P.length; i++) {
     const [x0, y0] = P[i];
     const [x1, y1] = P[i + 1];
     const dir = [x1 - x0, y1 - y0];
-    const len = Math.hypot(dir[0], dir[1]);
+    const real = Math.hypot(dir[0], dir[1]);
+    const len = wholePixelVectors ? Math.hypot(Math.floor(dir[0] / 16), Math.floor(dir[1] / 16)) * 16 : Math.hypot(dir[0], dir[1]);
+    if (len === 0) {
+      if (on && cur) {
+        cur.pts.push(P[i + 1]);
+        cur.dirs.push(dir);
+        cur.draws.push(tangents[i]);
+        cur.curves.push(curveFlags[i]);
+      }
+      continue;
+    }
     let t = 0;
     while (len - t > left || left === 0 && on) {
       t += left;
@@ -29590,10 +29659,12 @@ function dashPieces(P, tangents, pattern, shorten) {
         cur.pts.push(q);
         cur.dirs.push(dir);
         cur.draws.push(tangents[i]);
+        cur.curves.push(curveFlags[i]);
+        cur.endScale = cur.pts[0] === P[0] && q[0] === P[0][0] && q[1] === P[0][1] ? 1 : real / len;
         out.push(cur);
         cur = null;
       } else {
-        cur = { pts: [q], dirs: [], draws: [] };
+        cur = { pts: [q], dirs: [], draws: [], curves: [], startScale: real / len, endScale: 1 };
       }
       on = !on;
       idx = (idx + 1) % dashes.length;
@@ -29604,9 +29675,12 @@ function dashPieces(P, tangents, pattern, shorten) {
       cur.pts.push(P[i + 1]);
       cur.dirs.push(dir);
       cur.draws.push(tangents[i]);
+      cur.curves.push(curveFlags[i]);
+      cur.endScale = real / len;
     }
   }
   if (on && cur && cur.dirs.length > 0) {
+    cur.endScale = 1;
     out.push(cur);
   }
   return out;
@@ -29641,10 +29715,13 @@ function widenPath(path, opts) {
     if (dashed) {
       const run = closed ? [...P, P[0]] : P;
       const shorten = opts.cap === "flat" || opts.shortenDashes === false ? 0 : opts.width;
-      for (const piece of dashPieces(run, closed ? [...dirs, void 0] : dirs, opts.dashes, shorten)) {
+      const runTangents = closed ? [...dirs, void 0] : dirs;
+      const runCurves = closed ? [...curves, !!fig.roundWiden] : curves;
+      for (const piece of dashPieces(run, runTangents, runCurves, opts.dashes, shorten, !!opts.wholePixelDashVectors)) {
         const pts = [piece.pts[0]];
         const pdirs = [];
         const pdraws = [];
+        const pcurves = [];
         for (let i = 1; i < piece.pts.length; i++) {
           const q = piece.pts[i];
           const last = pts[pts.length - 1];
@@ -29652,9 +29729,10 @@ function widenPath(path, opts) {
             pts.push(q);
             pdirs.push(piece.dirs[i - 1]);
             pdraws.push(piece.draws[i - 1]);
+            pcurves.push(piece.curves[i - 1]);
           }
         }
-        outliner.open(pts, pdirs.length > 0 ? pdirs : [piece.dirs[0]], pdirs.length > 0 ? pdraws : [piece.draws[0]], pts.map(() => !!fig.roundWiden));
+        outliner.open(pts, pdirs.length > 0 ? pdirs : [piece.dirs[0]], pdirs.length > 0 ? pdraws : [piece.draws[0]], pdirs.length > 0 ? pcurves : [piece.curves[0]], piece.startScale, piece.endScale);
       }
     } else if (closed && P.length >= 3) {
       outliner.closed(P, dirs, curves);
@@ -29862,7 +29940,10 @@ function penWidenOptions(rCtx, opts = {}) {
     join: opts.roundPen ? "round" : opts.rectangle && !state.penExtended ? "miter" : !state.penExtended || joinBits === 0 ? "round" : joinBits === 4096 ? "bevel" : "miter",
     miterLimit: state.miterLimit ?? 10,
     dashes: dashes ? dashes.map((v) => v * 16) : null,
-    shortenDashes: (flags & 15) !== 7
+    shortenDashes: (flags & 15) !== 7,
+    // One logical unit per device pixel, unrotated: GDI measures dash
+    // segments in whole pixels.
+    wholePixelDashVectors: Math.abs(matrix[0] - 1) < 1e-6 && Math.abs(matrix[3] - 1) < 1e-6 && Math.abs(matrix[1]) < 1e-6 && Math.abs(matrix[2]) < 1e-6
   };
 }
 function hatchBackgroundMask(hatch) {
@@ -31402,6 +31483,15 @@ function illuminantCube(index) {
   }
   return cube;
 }
+var ILLUMINANT_MATRICES = {
+  1: [[1.64737, -0.59174, -0.05494], [0.02963, 0.96126, 881e-5], [0.2911, 0.40189, 0.30769]],
+  2: [[1.22319, -0.17765, -0.04515], [275e-5, 0.99054, 645e-5], [0.05514, 0.12795, 0.81816]],
+  3: [[1.02833, -0.01775, -0.01086], [861e-5, 0.98431, 728e-5], [-0.01616, -0.0448, 1.05995]],
+  4: [[1.17271, -0.14059, -0.03174], [-474e-5, 1.00266, 183e-5], [0.05783, 0.14663, 0.79683]],
+  5: [[1.10431, -0.08327, -0.02067], [-336e-5, 1.0025, 67e-5], [0.03344, 0.09072, 0.87702]],
+  7: [[0.92515, 0.05628, 0.01822], [318e-5, 0.9967, 34e-5], [-0.02086, -0.06817, 1.08783]],
+  8: [[1.31275, -0.26259, -0.04977], [-244e-5, 0.99695, 524e-5], [0.10854, 0.23443, 0.65827]]
+};
 function hasIlluminant(ca) {
   return illuminantCube(ca.illuminant) !== null;
 }
@@ -31472,9 +31562,9 @@ var BRIGHTNESS_STEP = 0.95625;
 var COLORFULNESS_SCALE = 1 / 100;
 var TINT_RADIANS = -Math.PI / 3;
 var RGB_TO_XYZ = [
-  [0.4124, 0.3576, 0.1805],
-  [0.2126, 0.7152, 0.0722],
-  [0.0193, 0.1192, 0.9505]
+  [0.4123908, 0.3575843, 0.1804808],
+  [0.212639, 0.7151687, 0.0721923],
+  [0.0193308, 0.1191948, 0.9505322]
 ];
 var XYZ_TO_RGB = invert3(RGB_TO_XYZ);
 var WHITE = RGB_TO_XYZ.map((row2) => row2[0] + row2[1] + row2[2]);
@@ -31502,13 +31592,13 @@ function adjustChroma(r, g, b, scale, angle) {
   const z = RGB_TO_XYZ[2][0] * r + RGB_TO_XYZ[2][1] * g + RGB_TO_XYZ[2][2] * b;
   const d = x + 15 * y + 3 * z;
   if (d <= 1e-9 || y <= 1e-9) {
-    return [r, g, b];
+    return [clamp255(r), clamp255(g), clamp255(b)];
   }
   const du = 4 * x / d - WHITE_U;
   const dv = 9 * y / d - WHITE_V;
   const cos = Math.cos(angle) * scale;
   const sin = Math.sin(angle) * scale;
-  const u = WHITE_U + du * cos - dv * sin;
+  const u = Math.max(0, WHITE_U + du * cos - dv * sin);
   const v = Math.max(1e-6, WHITE_V + du * sin + dv * cos);
   const x2 = y * 9 * u / (4 * v);
   const z2 = y * (12 - 3 * u - 20 * v) / (4 * v);
@@ -31527,6 +31617,7 @@ function adjustChroma(r, g, b, scale, angle) {
 }
 function colorAdjustmentMapper(ca) {
   const illuminant = illuminantCube(ca.illuminant);
+  const matrix = ILLUMINANT_MATRICES[ca.illuminant];
   const gammas = [ca.redGamma / 1e4, ca.greenGamma / 1e4, ca.blueGamma / 1e4];
   const black = ca.referenceBlack / 1e4 * 255;
   const span = Math.max(1e-6, ca.referenceWhite / 1e4 * 255 - black);
@@ -31555,16 +31646,17 @@ function colorAdjustmentMapper(ca) {
   }
   return (rgb2) => {
     let c = [rgb2 >> 16 & 255, rgb2 >> 8 & 255, rgb2 & 255];
-    if (illuminant) {
-      const mapped = lookupIlluminant(illuminant, c[0] / 255, c[1] / 255, c[2] / 255);
-      c = [mapped[0] * 255, mapped[1] * 255, mapped[2] * 255];
-    }
     if (chroma) {
-      if (!illuminant) {
-        c = c.map(exactPaletteLevel);
+      c = c.map(exactPaletteLevel);
+      if (matrix) {
+        const [r, g, b] = c;
+        c = matrix.map((row2) => row2[0] * r + row2[1] * g + row2[2] * b);
       }
-      c = adjustChroma(clamp255(c[0]), clamp255(c[1]), clamp255(c[2]), scale, angle);
+      c = adjustChroma(c[0], c[1], c[2], scale, angle);
       c = c.map(roundHalfUp);
+    } else if (illuminant) {
+      const mapped = lookupIlluminant(illuminant, c[0] / 255, c[1] / 255, c[2] / 255);
+      c = [mapped[0] * 255, mapped[1] * 255, mapped[2] * 255];
     }
     return curve(c[0], 0) << 16 | curve(c[1], 1) << 8 | curve(c[2], 2);
   };
@@ -31918,20 +32010,6 @@ function halftoneDespeckle(px, w, h) {
     }
   }
 }
-function halftoneSharpen(px, w, h) {
-  const src = px.slice();
-  const at = (x, y, c) => src[(Math.max(0, Math.min(h - 1, y)) * w + Math.max(0, Math.min(w - 1, x))) * 3 + c];
-  for (let y = 0; y < h; y++) {
-    for (let x = 0; x < w; x++) {
-      for (let c = 0; c < 3; c++) {
-        const v = at(x, y, c);
-        const sum = at(x - 1, y, c) + at(x + 1, y, c) + at(x, y - 1, c) + at(x, y + 1, c);
-        const out = v + Math.floor((4 * v - sum) / 8);
-        px[(y * w + x) * 3 + c] = out < 0 ? 0 : out > 255 ? 255 : out;
-      }
-    }
-  }
-}
 function halftonePairDespeckle(px, w, h, rows) {
   const n = rows ? w : h;
   const at = (k, side) => (rows ? side * w + k : k * w + side) * 3;
@@ -31954,6 +32032,46 @@ function halftoneKeepLast(px, w, h, rows) {
   }
   return out;
 }
+function halftoneReduceBoth(rect, sw, sh, dw, dh, flipX, flipY) {
+  const cols = halftoneReduceTaps(sw, dw, flipX);
+  const rows = halftoneReduceTaps(sh, dh, flipY);
+  const wide = new Int32Array(dw * sh * 3);
+  for (let y = 0; y < sh; y++) {
+    for (let x = 0; x < dw; x++) {
+      let total = 0;
+      for (const [, w] of cols[x]) total += w;
+      for (let c = 0; c < 3; c++) {
+        let sum = 0;
+        for (const [xx, w] of cols[x]) sum += rect[(y * sw + xx) * 3 + c] * w;
+        wide[(y * dw + x) * 3 + c] = Math.floor((2 * sum + total) / (2 * total));
+      }
+    }
+  }
+  const avg = new Float64Array(dw * dh * 3);
+  for (let y = 0; y < dh; y++) {
+    let total = 0;
+    for (const [, w] of rows[y]) total += w;
+    for (let x = 0; x < dw; x++) {
+      for (let c = 0; c < 3; c++) {
+        let sum = 0;
+        for (const [yy, w] of rows[y]) sum += wide[(yy * dw + x) * 3 + c] * w;
+        avg[(y * dw + x) * 3 + c] = sum / total;
+      }
+    }
+  }
+  const at = (x, y, c) => avg[(Math.max(0, Math.min(dh - 1, y)) * dw + Math.max(0, Math.min(dw - 1, x))) * 3 + c];
+  const out = new Int32Array(dw * dh * 3);
+  for (let y = 0; y < dh; y++) {
+    for (let x = 0; x < dw; x++) {
+      for (let c = 0; c < 3; c++) {
+        const v = at(x, y, c);
+        const sum = at(x - 1, y, c) + at(x + 1, y, c) + at(x, y - 1, c) + at(x, y + 1, c);
+        out[(y * dw + x) * 3 + c] = Math.max(0, Math.min(255, Math.floor(v + (4 * v - sum) / 8)));
+      }
+    }
+  }
+  return out;
+}
 function stretchHalftone(src, sx, sy, sw, sh, dw, dh, adjust, adjustAfterSampling = false, directDib = false, dither, curves) {
   const W = Math.max(0, Math.round(Math.abs(dw)));
   const H = Math.max(0, Math.round(Math.abs(dh)));
@@ -32005,11 +32123,18 @@ function stretchHalftone(src, sx, sy, sw, sh, dw, dh, adjust, adjustAfterSamplin
   }
   const reducedEven = W < SW ? SW % 2 === 0 && H % SH === 0 : SH % 2 === 0 && W % SW === 0;
   const ditherAtSource = dithered && (W >= SW && H >= SH || nativeMixed && reducedEven);
+  const extendRow = ditherAtSource && nativeMixed && W < SW && H > SH && !flipY;
+  if (extendRow) {
+    const bigger = new Int32Array((SH + 1) * SW * 3);
+    bigger.set(rect);
+    bigger.copyWithin(SH * SW * 3, (SH - 1) * SW * 3, SH * SW * 3);
+    rect = bigger;
+  }
   if (ditherAtSource) {
     ditherPixels(
       rect,
       SW,
-      SH,
+      extendRow ? SH + 1 : SH,
       dither,
       false,
       false,
@@ -32049,7 +32174,10 @@ function stretchHalftone(src, sx, sy, sw, sh, dw, dh, adjust, adjustAfterSamplin
     linear(rows, SH, H, flipY);
   }
   let out = nativeMixed ? new Float64Array(sampleW * sampleH * 3) : new Int32Array(W * H * 3);
-  for (let y = 0; y < sampleH; y++) {
+  if (reducing) {
+    out = halftoneReduceBoth(rect, SW, SH, W, H, flipX, flipY);
+  }
+  for (let y = 0; !reducing && y < sampleH; y++) {
     const ys = rows[y];
     for (let x = 0; x < sampleW; x++) {
       const xs = cols[x];
@@ -32074,15 +32202,38 @@ function stretchHalftone(src, sx, sy, sw, sh, dw, dh, adjust, adjustAfterSamplin
       out[o + 2] = round(b);
     }
   }
-  if (reducing) {
-    halftoneSharpen(out, W, H);
-  }
   if (nativeMixed) {
     const reduceHorizontal = W < SW;
+    let extension;
+    if (extendRow) {
+      const reduced = new Float64Array(sampleW * 3);
+      for (let x = 0; x < sampleW; x++) {
+        let total = 0;
+        const sums = [0, 0, 0];
+        for (const [xx, wx] of cols[x]) {
+          total += wx;
+          for (let c = 0; c < 3; c++) sums[c] += rect[(SH * SW + xx) * 3 + c] * wx;
+        }
+        for (let c = 0; c < 3; c++) reduced[x * 3 + c] = sums[c] / total;
+      }
+      extension = reduced;
+    }
     for (const horizontal of [reduceHorizontal, !reduceHorizontal]) {
       const source2 = out;
       out = new Int32Array(source2.length);
-      const at = (x, y, c) => source2[(Math.max(0, Math.min(sampleH - 1, y)) * sampleW + Math.max(0, Math.min(sampleW - 1, x))) * 3 + c];
+      const ext = extension;
+      const at = (x, y, c) => y >= sampleH && ext ? ext[Math.max(0, Math.min(sampleW - 1, x)) * 3 + c] : source2[(Math.max(0, Math.min(sampleH - 1, y)) * sampleW + Math.max(0, Math.min(sampleW - 1, x))) * 3 + c];
+      if (horizontal && ext) {
+        const next = new Int32Array(sampleW * 3);
+        for (let x = 0; x < sampleW; x++) {
+          for (let c = 0; c < 3; c++) {
+            const v = ext[x * 3 + c];
+            const sum = at(x - 1, sampleH, c) + at(x + 1, sampleH, c);
+            next[x * 3 + c] = Math.max(0, Math.min(255, Math.floor(v + (2 * v - sum) / 4)));
+          }
+        }
+        extension = next;
+      }
       for (let y = 0; y < sampleH; y++) {
         for (let x = 0; x < sampleW; x++) {
           for (let c = 0; c < 3; c++) {