@12ui/design 0.2.71 → 0.2.72

package/dist/scene/image-sheet.js

@@ -0,0 +1,118 @@
+import path from 'node:path';
+import sharp from 'sharp';
+import { SceneArgumentError } from './numbers.js';
+import { encodePng, readRgb, resizeRgb } from './raster.js';
+/**
+ * The two pictures the self-review loop looks at. A port of `sheet.py`.
+ *
+ * Sheet: SOURCE first, then each capture, in 640 px wide tiles with a large
+ * label on each (the file name; --note is added to the SOURCE tile).
+ * Blend (2x2 of 900 px tiles): SOURCE | RENDER above, 50% BLEND | DIFFERENCE x3
+ * below. Labels are Pango text rendered by sharp's bundled libvips: no font file
+ * and no network.
+ */
+const escapeMarkup = (text) => text
+    .replaceAll('&', '&')
+    .replaceAll('<', '&lt;')
+    .replaceAll('>', '&gt;');
+/** A white label on a black box, as a PNG to composite at (x - 6, y - 3) of the text origin. */
+const label = async (text, fontPx) => {
+    const { data, info } = await sharp({
+        text: {
+            text: `<span foreground="white">${escapeMarkup(text)}</span>`,
+            font: `sans ${fontPx}`,
+            dpi: 72,
+            rgba: true,
+        },
+    }).png().toBuffer({ resolveWithObject: true });
+    const width = info.width + 12;
+    const height = info.height + 6;
+    const input = await sharp({
+        create: { width, height, channels: 4, background: { r: 0, g: 0, b: 0, alpha: 1 } },
+    }).composite([{ input: data, left: 6, top: 3 }]).png().toBuffer();
+    return { input, width, height };
+};
+const composeGrid = async (args) => {
+    const rows = Math.ceil(args.tiles.length / args.columns);
+    const width = args.columns * args.tileWidth;
+    const height = rows * args.tileHeight;
+    const layers = [];
+    for (const [index, tile] of args.tiles.entries()) {
+        const left = (index % args.columns) * args.tileWidth;
+        const top = Math.floor(index / args.columns) * args.tileHeight;
+        const resized = await resizeRgb(tile.image, args.tileWidth, args.tileHeight);
+        layers.push({ input: await encodePng(resized), left, top });
+        const tag = await label(tile.name, args.fontPx);
+        layers.push({
+            input: tag.input,
+            left: Math.min(left + args.labelOffset - 6, Math.max(width - tag.width, 0)),
+            top: top + args.labelOffset - 3,
+        });
+    }
+    await sharp({ create: { width, height, channels: 3, background: args.background } })
+        .composite(layers)
+        .png()
+        .toFile(args.outFile);
+    return { width, height };
+};
+/** Python's round() for the tile height: ties go to the even neighbour. */
+const roundHalfEven = (value) => {
+    const floor = Math.floor(value);
+    const fraction = value - floor;
+    if (Math.abs(fraction - 0.5) > 1e-9)
+        return Math.round(value);
+    return floor % 2 === 0 ? floor : floor + 1;
+};
+const captureName = (file) => path.basename(file).replace(/\.[^.]*$/u, '');
+export const writeContactSheet = async (args) => {
+    const columns = args.columns ?? 3;
+    if (!Number.isInteger(columns) || columns < 1)
+        throw new SceneArgumentError('--cols must be a positive integer');
+    const source = await readRgb(args.sourceFile);
+    const tileWidth = 640;
+    const tileHeight = roundHalfEven((source.height * tileWidth) / source.width);
+    const tiles = [
+        { name: `SOURCE${args.note ? `  ${args.note}` : ''}`, image: source },
+        ...await Promise.all(args.captureFiles.map(async (file) => ({ name: captureName(file), image: await readRgb(file) }))),
+    ];
+    return composeGrid({
+        tiles,
+        columns,
+        tileWidth,
+        tileHeight,
+        background: { r: 24, g: 24, b: 24 },
+        labelOffset: 10,
+        fontPx: 22,
+        outFile: args.outFile,
+    });
+};
+export const writeBlend = async (args) => {
+    const source = await readRgb(args.sourceFile);
+    const render = await resizeRgb(await readRgb(args.renderFile), source.width, source.height);
+    const blend = Buffer.alloc(source.data.length);
+    const difference = Buffer.alloc(source.data.length);
+    for (let index = 0; index < source.data.length; index += 1) {
+        const a = source.data[index], b = render.data[index];
+        blend[index] = Math.round((a + b) / 2);
+        difference[index] = Math.min(255, Math.abs(a - b) * 3);
+    }
+    const tileWidth = 900;
+    const tileHeight = roundHalfEven((source.height * tileWidth) / source.width);
+    const same = { width: source.width, height: source.height };
+    return composeGrid({
+        tiles: [
+            { name: 'SOURCE', image: source },
+            { name: 'RENDER', image: render },
+            { name: '50% BLEND', image: { ...same, data: blend } },
+            { name: 'DIFFERENCE x3', image: { ...same, data: difference } },
+        ],
+        columns: 2,
+        tileWidth,
+        tileHeight,
+        background: { r: 0, g: 0, b: 0 },
+        labelOffset: 12,
+        fontPx: 26,
+        outFile: args.outFile,
+    });
+};
+//# sourceMappingURL=image-sheet.js.map

package/dist/scene/image-sheet.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"image-sheet.js","sourceRoot":"","sources":["../../src/scene/image-sheet.ts"],"names":[],"mappings":"AAAA,OAAO,IAAI,MAAM,WAAW,CAAC;AAE7B,OAAO,KAA8B,MAAM,OAAO,CAAC;AAEnD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClD,OAAO,EAAE,SAAS,EAAE,OAAO,EAAE,SAAS,EAAiB,MAAM,aAAa,CAAC;AAE3E;;;;;;;;GAQG;AAEH,MAAM,YAAY,GAAG,CAAC,IAAY,EAAU,EAAE,CAAC,IAAI;KAChD,UAAU,CAAC,GAAG,EAAE,OAAO,CAAC;KACxB,UAAU,CAAC,GAAG,EAAE,MAAM,CAAC;KACvB,UAAU,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;AAE3B,gGAAgG;AAChG,MAAM,KAAK,GAAG,KAAK,EAAE,IAAY,EAAE,MAAc,EAA6D,EAAE;IAC9G,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,MAAM,KAAK,CAAC;QACjC,IAAI,EAAE;YACJ,IAAI,EAAE,4BAA4B,YAAY,CAAC,IAAI,CAAC,SAAS;YAC7D,IAAI,EAAE,QAAQ,MAAM,EAAE;YACtB,GAAG,EAAE,EAAE;YACP,IAAI,EAAE,IAAI;SACX;KACF,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,EAAE,iBAAiB,EAAE,IAAI,EAAE,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC;IAC9B,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;IAC/B,MAAM,KAAK,GAAG,MAAM,KAAK,CAAC;QACxB,MAAM,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,UAAU,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,EAAE;KACnF,CAAC,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,QAAQ,EAAE,CAAC;IAClE,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC;AAClC,CAAC,CAAC;AAIF,MAAM,WAAW,GAAG,KAAK,EAAE,IAS1B,EAA8C,EAAE;IAC/C,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;IACzD,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,GAAG,IAAI,CAAC,SAAS,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,GAAG,IAAI,CAAC,UAAU,CAAC;IACtC,MAAM,MAAM,GAAqB,EAAE,CAAC;IACpC,KAAK,MAAM,CAAC,KAAK,EAAE,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,OAAO,EAAE,EAAE,CAAC;QACjD,MAAM,IAAI,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC;QACrD,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;QAC/D,MAAM,OAAO,GAAG,MAAM,SAAS,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,UAAU,CAAC,CAAC;QAC7E,MAAM,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC,OAAO,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC;QAC5D,MAAM,GAAG,GAAG,MAAM,KAAK,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAChD,MAAM,CAAC,IAAI,CAAC;YACV,KAAK,EAAE,GAAG,CAAC,KAAK;YAChB,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;YAC3E,GAAG,EAAE,GAAG,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC;SAChC,CAAC,CAAC;IACL,CAAC;IACD,MAAM,KAAK,CAAC,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,UAAU,EAAE,IAAI,CAAC,UAAU,EAAE,EAAE,CAAC;SACjF,SAAS,CAAC,MAAM,CAAC;SACjB,GAAG,EAAE;SACL,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IACxB,OAAO,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC;AAC3B,CAAC,CAAC;AAEF,2EAA2E;AAC3E,MAAM,aAAa,GAAG,CAAC,KAAa,EAAU,EAAE;IAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAChC,MAAM,QAAQ,GAAG,KAAK,GAAG,KAAK,CAAC;IAC/B,IAAI,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,GAAG,IAAI;QAAE,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;IAC9D,OAAO,KAAK,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;AAC7C,CAAC,CAAC;AAEF,MAAM,WAAW,GAAG,CAAC,IAAY,EAAU,EAAE,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,OAAO,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC;AAE3F,MAAM,CAAC,MAAM,iBAAiB,GAAG,KAAK,EAAE,IAMvC,EAA8C,EAAE;IAC/C,MAAM,OAAO,GAAG,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC;IAClC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,IAAI,OAAO,GAAG,CAAC;QAAE,MAAM,IAAI,kBAAkB,CAAC,mCAAmC,CAAC,CAAC;IACjH,MAAM,MAAM,GAAG,MAAM,OAAO,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;IAC9C,MAAM,SAAS,GAAG,GAAG,CAAC;IACtB,MAAM,UAAU,GAAG,aAAa,CAAC,CAAC,MAAM,CAAC,MAAM,GAAG,SAAS,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAW;QACpB,EAAE,IAAI,EAAE,SAAS,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE;QACrE,GAAG,MAAM,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,WAAW,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,MAAM,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;KACvH,CAAC;IACF,OAAO,WAAW,CAAC;QACjB,KAAK;QACL,OAAO;QACP,SAAS;QACT,UAAU;QACV,UAAU,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE;QACnC,WAAW,EAAE,EAAE;QACf,MAAM,EAAE,EAAE;QACV,OAAO,EAAE,IAAI,CAAC,OAAO;KACtB,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,MAAM,CAAC,MAAM,UAAU,GAAG,KAAK,EAAE,IAIhC,EAA8C,EAAE;IAC/C,MAAM,MAAM,GAAG,MAAM,OAAO,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;IAC9C,MAAM,MAAM,GAAG,MAAM,SAAS,CAAC,MAAM,OAAO,CAAC,IAAI,CAAC,UAAU,CAAC,EAAE,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5F,MAAM,KAAK,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC/C,MAAM,UAAU,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACpD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3D,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAE,EAAE,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,KAAK,CAAE,CAAC;QACvD,KAAK,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvC,UAAU,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACzD,CAAC;IACD,MAAM,SAAS,GAAG,GAAG,CAAC;IACtB,MAAM,UAAU,GAAG,aAAa,CAAC,CAAC,MAAM,CAAC,MAAM,GAAG,SAAS,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAC7E,MAAM,IAAI,GAAG,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,CAAC;IAC5D,OAAO,WAAW,CAAC;QACjB,KAAK,EAAE;YACL,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE;YACjC,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE;YACjC,EAAE,IAAI,EAAE,WAAW,EAAE,KAAK,EAAE,EAAE,GAAG,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE;YACtD,EAAE,IAAI,EAAE,eAAe,EAAE,KAAK,EAAE,EAAE,GAAG,IAAI,EAAE,IAAI,EAAE,UAAU,EAAE,EAAE;SAChE;QACD,OAAO,EAAE,CAAC;QACV,SAAS;QACT,UAAU;QACV,UAAU,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE;QAChC,WAAW,EAAE,EAAE;QACf,MAAM,EAAE,EAAE;QACV,OAAO,EAAE,IAAI,CAAC,OAAO;KACtB,CAAC,CAAC;AACL,CAAC,CAAC"}

package/dist/scene/least-squares.d.ts

@@ -0,0 +1,31 @@
+/**
+ * Box-constrained Levenberg-Marquardt with a numeric Jacobian, and a multistart
+ * around it. A line-for-line port of the proven Python solver
+ * (`solve_camera.py` lm/multistart): the same damping schedule, the same
+ * forward-difference step, the same clamping, so a spec solves to the same
+ * minimum in both.
+ */
+export type Residuals = (parameters: Float64Array) => Float64Array;
+export declare class SingularMatrixError extends Error {
+    constructor();
+}
+/**
+ * Solves `matrix * x = rhs` by Gaussian elimination with partial pivoting.
+ * An exactly zero pivot is a singular system, as LAPACK reports it.
+ */
+export declare const solveLinearSystem: (matrix: readonly Float64Array[], rhs: Float64Array) => Float64Array;
+export type LeastSquaresResult = {
+    parameters: Float64Array;
+    cost: number;
+};
+export declare const levenbergMarquardt: (residuals: Residuals, start: ArrayLike<number>, low: ArrayLike<number>, high: ArrayLike<number>, iterations?: number) => LeastSquaresResult;
+/**
+ * Every start gets a short coarse run; the best `keep` are refined fully and
+ * the lowest cost wins. Ties keep the earlier start, as Python's stable sort did.
+ */
+export declare const multistart: (residuals: Residuals, starts: readonly ArrayLike<number>[], low: ArrayLike<number>, high: ArrayLike<number>, { keep, coarse, refine }?: {
+    keep?: number;
+    coarse?: number;
+    refine?: number;
+}) => LeastSquaresResult;
+//# sourceMappingURL=least-squares.d.ts.map

package/dist/scene/least-squares.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"least-squares.d.ts","sourceRoot":"","sources":["../../src/scene/least-squares.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAEH,MAAM,MAAM,SAAS,GAAG,CAAC,UAAU,EAAE,YAAY,KAAK,YAAY,CAAC;AAEnE,qBAAa,mBAAoB,SAAQ,KAAK;;CAK7C;AAgBD;;;GAGG;AACH,eAAO,MAAM,iBAAiB,GAAI,QAAQ,SAAS,YAAY,EAAE,EAAE,KAAK,YAAY,KAAG,YAoCtF,CAAC;AAEF,MAAM,MAAM,kBAAkB,GAAG;IAAE,UAAU,EAAE,YAAY,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CAAC;AAE5E,eAAO,MAAM,kBAAkB,GAC7B,WAAW,SAAS,EACpB,OAAO,SAAS,CAAC,MAAM,CAAC,EACxB,KAAK,SAAS,CAAC,MAAM,CAAC,EACtB,MAAM,SAAS,CAAC,MAAM,CAAC,EACvB,mBAAe,KACd,kBAmEF,CAAC;AAEF;;;GAGG;AACH,eAAO,MAAM,UAAU,GACrB,WAAW,SAAS,EACpB,QAAQ,SAAS,SAAS,CAAC,MAAM,CAAC,EAAE,EACpC,KAAK,SAAS,CAAC,MAAM,CAAC,EACtB,MAAM,SAAS,CAAC,MAAM,CAAC,EACvB,2BAAyC;IAAE,IAAI,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,MAAM,CAAA;CAAO,KAChG,kBAeF,CAAC"}

package/dist/scene/least-squares.js

@@ -0,0 +1,173 @@
+/**
+ * Box-constrained Levenberg-Marquardt with a numeric Jacobian, and a multistart
+ * around it. A line-for-line port of the proven Python solver
+ * (`solve_camera.py` lm/multistart): the same damping schedule, the same
+ * forward-difference step, the same clamping, so a spec solves to the same
+ * minimum in both.
+ */
+export class SingularMatrixError extends Error {
+    constructor() {
+        super('singular matrix');
+        this.name = 'SingularMatrixError';
+    }
+}
+const clamp = (values, low, high) => {
+    const out = new Float64Array(values.length);
+    for (let index = 0; index < values.length; index += 1) {
+        out[index] = Math.min(Math.max(values[index], low[index]), high[index]);
+    }
+    return out;
+};
+const sumOfSquares = (values) => {
+    let total = 0;
+    for (const value of values)
+        total += value * value;
+    return total;
+};
+/**
+ * Solves `matrix * x = rhs` by Gaussian elimination with partial pivoting.
+ * An exactly zero pivot is a singular system, as LAPACK reports it.
+ */
+export const solveLinearSystem = (matrix, rhs) => {
+    const size = rhs.length;
+    const a = matrix.map((row) => Float64Array.from(row));
+    const b = Float64Array.from(rhs);
+    for (let column = 0; column < size; column += 1) {
+        let pivotRow = column;
+        let pivotSize = Math.abs(a[column][column]);
+        for (let row = column + 1; row < size; row += 1) {
+            const candidate = Math.abs(a[row][column]);
+            if (candidate > pivotSize) {
+                pivotSize = candidate;
+                pivotRow = row;
+            }
+        }
+        if (pivotSize === 0)
+            throw new SingularMatrixError();
+        if (pivotRow !== column) {
+            [a[column], a[pivotRow]] = [a[pivotRow], a[column]];
+            [b[column], b[pivotRow]] = [b[pivotRow], b[column]];
+        }
+        const pivot = a[column];
+        for (let row = column + 1; row < size; row += 1) {
+            const target = a[row];
+            const factor = target[column] / pivot[column];
+            if (factor === 0)
+                continue;
+            for (let k = column; k < size; k += 1)
+                target[k] -= factor * pivot[k];
+            b[row] -= factor * b[column];
+        }
+    }
+    const x = new Float64Array(size);
+    for (let row = size - 1; row >= 0; row -= 1) {
+        let total = b[row];
+        const current = a[row];
+        for (let k = row + 1; k < size; k += 1)
+            total -= current[k] * x[k];
+        x[row] = total / current[row];
+    }
+    return x;
+};
+export const levenbergMarquardt = (residuals, start, low, high, iterations = 60) => {
+    let p = clamp(start, low, high);
+    let lambda = 1e-2;
+    let r = residuals(p);
+    let cost = sumOfSquares(r);
+    const n = p.length;
+    for (let iteration = 0; iteration < iterations; iteration += 1) {
+        const m = r.length;
+        // Forward-difference Jacobian, column by column.
+        const jacobian = [];
+        for (let i = 0; i < n; i += 1) {
+            const step = 1e-6 * Math.max(1, Math.abs(p[i]));
+            const q = Float64Array.from(p);
+            q[i] += step;
+            const rq = residuals(q);
+            const column = new Float64Array(m);
+            for (let k = 0; k < m; k += 1)
+                column[k] = (rq[k] - r[k]) / step;
+            jacobian.push(column);
+        }
+        const normal = Array.from({ length: n }, () => new Float64Array(n));
+        const gradient = new Float64Array(n);
+        for (let i = 0; i < n; i += 1) {
+            const ci = jacobian[i];
+            let g = 0;
+            for (let k = 0; k < m; k += 1)
+                g += ci[k] * r[k];
+            gradient[i] = g;
+            for (let j = i; j < n; j += 1) {
+                const cj = jacobian[j];
+                let total = 0;
+                for (let k = 0; k < m; k += 1)
+                    total += ci[k] * cj[k];
+                normal[i][j] = total;
+                normal[j][i] = total;
+            }
+        }
+        let improved = false;
+        for (let attempt = 0; attempt < 12; attempt += 1) {
+            const damped = normal.map((row, i) => {
+                const copy = Float64Array.from(row);
+                copy[i] += lambda * (row[i] + 1e-9);
+                return copy;
+            });
+            let step;
+            try {
+                step = solveLinearSystem(damped, gradient.map((value) => -value));
+            }
+            catch (error) {
+                if (error instanceof SingularMatrixError) {
+                    lambda *= 10;
+                    continue;
+                }
+                throw error;
+            }
+            const candidate = clamp(p.map((value, index) => value + step[index]), low, high);
+            const rq = residuals(candidate);
+            const cq = sumOfSquares(rq);
+            if (cq < cost) {
+                p = candidate;
+                r = rq;
+                cost = cq;
+                lambda = Math.max(lambda / 3, 1e-9);
+                improved = true;
+                break;
+            }
+            lambda *= 4;
+        }
+        if (!improved || cost < 1e-10)
+            break;
+    }
+    return { parameters: p, cost };
+};
+/**
+ * Every start gets a short coarse run; the best `keep` are refined fully and
+ * the lowest cost wins. Ties keep the earlier start, as Python's stable sort did.
+ */
+export const multistart = (residuals, starts, low, high, { keep = 6, coarse = 25, refine = 200 } = {}) => {
+    const scored = starts.map((start, index) => ({
+        index,
+        ...levenbergMarquardt(residuals, start, low, high, coarse),
+    }));
+    scored.sort((left, right) => (compareCost(left.cost, right.cost) || left.index - right.index));
+    let best = null;
+    for (const candidate of scored.slice(0, keep)) {
+        const refined = levenbergMarquardt(residuals, candidate.parameters, low, high, refine);
+        if (best === null || refined.cost < best.cost)
+            best = refined;
+    }
+    if (best === null)
+        throw new Error('multistart needs at least one start');
+    return best;
+};
+/** NaN sorts last, like a failed start should. */
+const compareCost = (left, right) => {
+    if (Number.isNaN(left))
+        return Number.isNaN(right) ? 0 : 1;
+    if (Number.isNaN(right))
+        return -1;
+    return left - right;
+};
+//# sourceMappingURL=least-squares.js.map

package/dist/scene/least-squares.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"least-squares.js","sourceRoot":"","sources":["../../src/scene/least-squares.ts"],"names":[],"mappings":"AAAA;;;;;;GAMG;AAIH,MAAM,OAAO,mBAAoB,SAAQ,KAAK;IAC5C;QACE,KAAK,CAAC,iBAAiB,CAAC,CAAC;QACzB,IAAI,CAAC,IAAI,GAAG,qBAAqB,CAAC;IACpC,CAAC;CACF;AAED,MAAM,KAAK,GAAG,CAAC,MAAyB,EAAE,GAAsB,EAAE,IAAuB,EAAgB,EAAE;IACzG,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5C,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,CAAC,MAAM,EAAE,KAAK,IAAI,CAAC,EAAE,CAAC;QACtD,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,KAAK,CAAE,EAAE,GAAG,CAAC,KAAK,CAAE,CAAC,EAAE,IAAI,CAAC,KAAK,CAAE,CAAC,CAAC;IAC7E,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC,CAAC;AAEF,MAAM,YAAY,GAAG,CAAC,MAAoB,EAAU,EAAE;IACpD,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,KAAK,IAAI,MAAM;QAAE,KAAK,IAAI,KAAK,GAAG,KAAK,CAAC;IACnD,OAAO,KAAK,CAAC;AACf,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,MAA+B,EAAE,GAAiB,EAAgB,EAAE;IACpG,MAAM,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC;IACxB,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACtD,MAAM,CAAC,GAAG,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACjC,KAAK,IAAI,MAAM,GAAG,CAAC,EAAE,MAAM,GAAG,IAAI,EAAE,MAAM,IAAI,CAAC,EAAE,CAAC;QAChD,IAAI,QAAQ,GAAG,MAAM,CAAC;QACtB,IAAI,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAE,CAAC,MAAM,CAAE,CAAC,CAAC;QAC9C,KAAK,IAAI,GAAG,GAAG,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;YAChD,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAE,CAAC,MAAM,CAAE,CAAC,CAAC;YAC7C,IAAI,SAAS,GAAG,SAAS,EAAE,CAAC;gBAC1B,SAAS,GAAG,SAAS,CAAC;gBACtB,QAAQ,GAAG,GAAG,CAAC;YACjB,CAAC;QACH,CAAC;QACD,IAAI,SAAS,KAAK,CAAC;YAAE,MAAM,IAAI,mBAAmB,EAAE,CAAC;QACrD,IAAI,QAAQ,KAAK,MAAM,EAAE,CAAC;YACxB,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAE,EAAE,CAAC,CAAC,MAAM,CAAE,CAAC,CAAC;YACtD,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAE,EAAE,CAAC,CAAC,MAAM,CAAE,CAAC,CAAC;QACxD,CAAC;QACD,MAAM,KAAK,GAAG,CAAC,CAAC,MAAM,CAAE,CAAC;QACzB,KAAK,IAAI,GAAG,GAAG,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;YAChD,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,CAAE,CAAC;YACvB,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAE,GAAG,KAAK,CAAC,MAAM,CAAE,CAAC;YAChD,IAAI,MAAM,KAAK,CAAC;gBAAE,SAAS;YAC3B,KAAK,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,IAAI,CAAC;gBAAE,MAAM,CAAC,CAAC,CAAE,IAAI,MAAM,GAAG,KAAK,CAAC,CAAC,CAAE,CAAC;YACxE,CAAC,CAAC,GAAG,CAAE,IAAI,MAAM,GAAG,CAAC,CAAC,MAAM,CAAE,CAAC;QACjC,CAAC;IACH,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC;IACjC,KAAK,IAAI,GAAG,GAAG,IAAI,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAC5C,IAAI,KAAK,GAAG,CAAC,CAAC,GAAG,CAAE,CAAC;QACpB,MAAM,OAAO,GAAG,CAAC,CAAC,GAAG,CAAE,CAAC;QACxB,KAAK,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,IAAI,CAAC;YAAE,KAAK,IAAI,OAAO,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;QACrE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,OAAO,CAAC,GAAG,CAAE,CAAC;IACjC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAIF,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAChC,SAAoB,EACpB,KAAwB,EACxB,GAAsB,EACtB,IAAuB,EACvB,UAAU,GAAG,EAAE,EACK,EAAE;IACtB,IAAI,CAAC,GAAG,KAAK,CAAC,KAAK,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;IAChC,IAAI,MAAM,GAAG,IAAI,CAAC;IAClB,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;IACrB,IAAI,IAAI,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,KAAK,IAAI,SAAS,GAAG,CAAC,EAAE,SAAS,GAAG,UAAU,EAAE,SAAS,IAAI,CAAC,EAAE,CAAC;QAC/D,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,iDAAiD;QACjD,MAAM,QAAQ,GAAmB,EAAE,CAAC;QACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;YAC9B,MAAM,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,CAAC;YACjD,MAAM,CAAC,GAAG,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAC/B,CAAC,CAAC,CAAC,CAAE,IAAI,IAAI,CAAC;YACd,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YACxB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;YACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC;gBAAE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC,GAAG,IAAI,CAAC;YACnE,QAAQ,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACxB,CAAC;QACD,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,MAAM,QAAQ,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC;YACxB,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC;gBAAE,CAAC,IAAI,EAAE,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;YACnD,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;gBAC9B,MAAM,EAAE,GAAG,QAAQ,CAAC,CAAC,CAAE,CAAC;gBACxB,IAAI,KAAK,GAAG,CAAC,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC;oBAAE,KAAK,IAAI,EAAE,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC,CAAC,CAAE,CAAC;gBACxD,MAAM,CAAC,CAAC,CAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;gBACtB,MAAM,CAAC,CAAC,CAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;YACxB,CAAC;QACH,CAAC;QACD,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,KAAK,IAAI,OAAO,GAAG,CAAC,EAAE,OAAO,GAAG,EAAE,EAAE,OAAO,IAAI,CAAC,EAAE,CAAC;YACjD,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE;gBACnC,MAAM,IAAI,GAAG,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;gBACpC,IAAI,CAAC,CAAC,CAAE,IAAI,MAAM,GAAG,CAAC,GAAG,CAAC,CAAC,CAAE,GAAG,IAAI,CAAC,CAAC;gBACtC,OAAO,IAAI,CAAC;YACd,CAAC,CAAC,CAAC;YACH,IAAI,IAAkB,CAAC;YACvB,IAAI,CAAC;gBACH,IAAI,GAAG,iBAAiB,CAAC,MAAM,EAAE,QAAQ,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;YACpE,CAAC;YAAC,OAAO,KAAK,EAAE,CAAC;gBACf,IAAI,KAAK,YAAY,mBAAmB,EAAE,CAAC;oBACzC,MAAM,IAAI,EAAE,CAAC;oBACb,SAAS;gBACX,CAAC;gBACD,MAAM,KAAK,CAAC;YACd,CAAC;YACD,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAE,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;YAClF,MAAM,EAAE,GAAG,SAAS,CAAC,SAAS,CAAC,CAAC;YAChC,MAAM,EAAE,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC;YAC5B,IAAI,EAAE,GAAG,IAAI,EAAE,CAAC;gBACd,CAAC,GAAG,SAAS,CAAC;gBACd,CAAC,GAAG,EAAE,CAAC;gBACP,IAAI,GAAG,EAAE,CAAC;gBACV,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC;gBACpC,QAAQ,GAAG,IAAI,CAAC;gBAChB,MAAM;YACR,CAAC;YACD,MAAM,IAAI,CAAC,CAAC;QACd,CAAC;QACD,IAAI,CAAC,QAAQ,IAAI,IAAI,GAAG,KAAK;YAAE,MAAM;IACvC,CAAC;IACD,OAAO,EAAE,UAAU,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC;AACjC,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,UAAU,GAAG,CACxB,SAAoB,EACpB,MAAoC,EACpC,GAAsB,EACtB,IAAuB,EACvB,EAAE,IAAI,GAAG,CAAC,EAAE,MAAM,GAAG,EAAE,EAAE,MAAM,GAAG,GAAG,KAA0D,EAAE,EAC7E,EAAE;IACtB,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;QAC3C,KAAK;QACL,GAAG,kBAAkB,CAAC,SAAS,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,CAAC;KAC3D,CAAC,CAAC,CAAC;IACJ,MAAM,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC,CAC3B,WAAW,CAAC,IAAI,CAAC,IAAI,EAAE,KAAK,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAC/D,CAAC,CAAC;IACH,IAAI,IAAI,GAA8B,IAAI,CAAC;IAC3C,KAAK,MAAM,SAAS,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC;QAC9C,MAAM,OAAO,GAAG,kBAAkB,CAAC,SAAS,EAAE,SAAS,CAAC,UAAU,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,CAAC,CAAC;QACvF,IAAI,IAAI,KAAK,IAAI,IAAI,OAAO,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI;YAAE,IAAI,GAAG,OAAO,CAAC;IAChE,CAAC;IACD,IAAI,IAAI,KAAK,IAAI;QAAE,MAAM,IAAI,KAAK,CAAC,qCAAqC,CAAC,CAAC;IAC1E,OAAO,IAAI,CAAC;AACd,CAAC,CAAC;AAEF,kDAAkD;AAClD,MAAM,WAAW,GAAG,CAAC,IAAY,EAAE,KAAa,EAAU,EAAE;IAC1D,IAAI,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,IAAI,MAAM,CAAC,KAAK,CAAC,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC,CAAC;IACnC,OAAO,IAAI,GAAG,KAAK,CAAC;AACtB,CAAC,CAAC"}

package/dist/scene/numbers.d.ts

@@ -0,0 +1,17 @@
+/** A malformed scene argument: a refusal, raised before any work starts. */
+export declare class SceneArgumentError extends Error {
+    constructor(message: string);
+}
+/** Round to `digits` decimals for JSON output (the Python helpers' `round`). */
+export declare const roundTo: (value: number, digits: number) => number;
+/** Fixed-point text, like Python's `f"{value:.Nf}"`. */
+export declare const fixed: (value: number, digits: number) => string;
+/** Parses "a,b" or "a,b,c" into numbers, naming the argument when it is malformed. */
+export declare const parseNumberTuple: (text: string, count: number, what: string) => number[];
+/** Parses a box "x0,y0,x1,y1" into integers (fractions truncate, as Python's int() did). */
+export declare const parseBox: (text: string, what: string) => [number, number, number, number];
+/** Splits a `;`-separated list, dropping empty entries. */
+export declare const splitList: (text: string) => string[];
+/** "WxH" into two positive integers. */
+export declare const parseSize: (text: string, what: string) => [number, number];
+//# sourceMappingURL=numbers.d.ts.map

package/dist/scene/numbers.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"numbers.d.ts","sourceRoot":"","sources":["../../src/scene/numbers.ts"],"names":[],"mappings":"AAAA,4EAA4E;AAC5E,qBAAa,kBAAmB,SAAQ,KAAK;gBAC/B,OAAO,EAAE,MAAM;CAI5B;AAED,gFAAgF;AAChF,eAAO,MAAM,OAAO,GAAI,OAAO,MAAM,EAAE,QAAQ,MAAM,KAAG,MAGvD,CAAC;AAEF,wDAAwD;AACxD,eAAO,MAAM,KAAK,GAAI,OAAO,MAAM,EAAE,QAAQ,MAAM,KAAG,MAA+B,CAAC;AAEtF,sFAAsF;AACtF,eAAO,MAAM,gBAAgB,GAAI,MAAM,MAAM,EAAE,OAAO,MAAM,EAAE,MAAM,MAAM,KAAG,MAAM,EAOlF,CAAC;AAEF,4FAA4F;AAC5F,eAAO,MAAM,QAAQ,GAAI,MAAM,MAAM,EAAE,MAAM,MAAM,KAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAGpF,CAAC;AAEF,2DAA2D;AAC3D,eAAO,MAAM,SAAS,GAAI,MAAM,MAAM,KAAG,MAAM,EAAgE,CAAC;AAEhH,wCAAwC;AACxC,eAAO,MAAM,SAAS,GAAI,MAAM,MAAM,EAAE,MAAM,MAAM,KAAG,CAAC,MAAM,EAAE,MAAM,CAMrE,CAAC"}

package/dist/scene/numbers.js

@@ -0,0 +1,39 @@
+/** A malformed scene argument: a refusal, raised before any work starts. */
+export class SceneArgumentError extends Error {
+    constructor(message) {
+        super(message);
+        this.name = 'SceneArgumentError';
+    }
+}
+/** Round to `digits` decimals for JSON output (the Python helpers' `round`). */
+export const roundTo = (value, digits) => {
+    const rounded = Number(value.toFixed(digits));
+    return Object.is(rounded, -0) ? 0 : rounded;
+};
+/** Fixed-point text, like Python's `f"{value:.Nf}"`. */
+export const fixed = (value, digits) => value.toFixed(digits);
+/** Parses "a,b" or "a,b,c" into numbers, naming the argument when it is malformed. */
+export const parseNumberTuple = (text, count, what) => {
+    const parts = text.split(',').map((part) => part.trim());
+    const values = parts.map(Number);
+    if (parts.length !== count || parts.some((part) => part === '') || values.some((value) => !Number.isFinite(value))) {
+        throw new SceneArgumentError(`${what} must be ${count} comma-separated numbers, got "${text}"`);
+    }
+    return values;
+};
+/** Parses a box "x0,y0,x1,y1" into integers (fractions truncate, as Python's int() did). */
+export const parseBox = (text, what) => {
+    const [x0, y0, x1, y1] = parseNumberTuple(text, 4, what).map(Math.trunc);
+    return [x0, y0, x1, y1];
+};
+/** Splits a `;`-separated list, dropping empty entries. */
+export const splitList = (text) => text.split(';').map((part) => part.trim()).filter(Boolean);
+/** "WxH" into two positive integers. */
+export const parseSize = (text, what) => {
+    const match = /^(\d+)x(\d+)$/u.exec(text.trim());
+    if (!match || Number(match[1]) < 1 || Number(match[2]) < 1) {
+        throw new SceneArgumentError(`${what} must be <width>x<height> in pixels, got "${text}"`);
+    }
+    return [Number(match[1]), Number(match[2])];
+};
+//# sourceMappingURL=numbers.js.map

package/dist/scene/numbers.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"numbers.js","sourceRoot":"","sources":["../../src/scene/numbers.ts"],"names":[],"mappings":"AAAA,4EAA4E;AAC5E,MAAM,OAAO,kBAAmB,SAAQ,KAAK;IAC3C,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,oBAAoB,CAAC;IACnC,CAAC;CACF;AAED,gFAAgF;AAChF,MAAM,CAAC,MAAM,OAAO,GAAG,CAAC,KAAa,EAAE,MAAc,EAAU,EAAE;IAC/D,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC;IAC9C,OAAO,MAAM,CAAC,EAAE,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC;AAC9C,CAAC,CAAC;AAEF,wDAAwD;AACxD,MAAM,CAAC,MAAM,KAAK,GAAG,CAAC,KAAa,EAAE,MAAc,EAAU,EAAE,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;AAEtF,sFAAsF;AACtF,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAAC,IAAY,EAAE,KAAa,EAAE,IAAY,EAAY,EAAE;IACtF,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;IACzD,MAAM,MAAM,GAAG,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;IACjC,IAAI,KAAK,CAAC,MAAM,KAAK,KAAK,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;QACnH,MAAM,IAAI,kBAAkB,CAAC,GAAG,IAAI,YAAY,KAAK,kCAAkC,IAAI,GAAG,CAAC,CAAC;IAClG,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC,CAAC;AAEF,4FAA4F;AAC5F,MAAM,CAAC,MAAM,QAAQ,GAAG,CAAC,IAAY,EAAE,IAAY,EAAoC,EAAE;IACvF,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,gBAAgB,CAAC,IAAI,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAqC,CAAC;IAC7G,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAC1B,CAAC,CAAC;AAEF,2DAA2D;AAC3D,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,IAAY,EAAY,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC;AAEhH,wCAAwC;AACxC,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,IAAY,EAAE,IAAY,EAAoB,EAAE;IACxE,MAAM,KAAK,GAAG,gBAAgB,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;IACjD,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;QAC3D,MAAM,IAAI,kBAAkB,CAAC,GAAG,IAAI,6CAA6C,IAAI,GAAG,CAAC,CAAC;IAC5F,CAAC;IACD,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC9C,CAAC,CAAC"}

package/dist/scene/placement.d.ts

@@ -0,0 +1,43 @@
+import { type Camera, type ImageSize, type Vec3 } from './camera-model.js';
+/**
+ * Placing objects through a solved camera: pixel rays onto a plane, world
+ * points back to pixels, and the pose of a rectangle from its four corners.
+ * Ports of `solve_camera.py ray|project|rect`, with the sized rectangle fixed:
+ * the Python helper crashed whenever both sides were known.
+ */
+export type RayLine = {
+    u: number;
+    v: number;
+    point: Vec3;
+    behind: boolean;
+};
+export declare const castRays: (camera: Camera, size: ImageSize, pixels: ReadonlyArray<[number, number]>, planeY: number) => RayLine[];
+export declare const formatRayLine: (line: RayLine) => string;
+export declare const formatProjectLine: (point: Vec3, u: number, v: number) => string;
+export type RectPose = {
+    position: Vec3;
+    w: number;
+    h: number;
+    rotation: {
+        x: number;
+        y: number;
+        z: number;
+        order: 'YXZ';
+    };
+    rmsPx?: number;
+    note: string;
+};
+/** A rectangle lying flat on the plane y = planeY, any size. Corners TL;TR;BR;BL in pixels. */
+export declare const rectOnPlane: (camera: Camera, size: ImageSize, corners: ReadonlyArray<[number, number]>, planeY: number) => {
+    pose: RectPose;
+    behind: boolean;
+};
+/**
+ * Free pose (position, lean, roll) of a rectangle whose width, height or both
+ * are known. With one side known the other is solved, less well determined.
+ */
+export declare const rectOfKnownSize: (camera: Camera, size: ImageSize, corners: ReadonlyArray<[number, number]>, known: {
+    width?: number;
+    height?: number;
+}) => RectPose;
+//# sourceMappingURL=placement.d.ts.map

package/dist/scene/placement.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"placement.d.ts","sourceRoot":"","sources":["../../src/scene/placement.ts"],"names":[],"mappings":"AAAA,OAAO,EAKL,KAAK,MAAM,EACX,KAAK,SAAS,EACd,KAAK,IAAI,EACV,MAAM,mBAAmB,CAAC;AAK3B;;;;;GAKG;AAEH,MAAM,MAAM,OAAO,GAAG;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,CAAC,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,IAAI,CAAC;IAAC,MAAM,EAAE,OAAO,CAAA;CAAE,CAAC;AAE7E,eAAO,MAAM,QAAQ,GACnB,QAAQ,MAAM,EACd,MAAM,SAAS,EACf,QAAQ,aAAa,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACvC,QAAQ,MAAM,KACb,OAAO,EAAiF,CAAC;AAE5F,eAAO,MAAM,aAAa,GAAI,MAAM,OAAO,KAAG,MAG7C,CAAC;AAEF,eAAO,MAAM,iBAAiB,GAAI,OAAO,IAAI,EAAE,GAAG,MAAM,EAAE,GAAG,MAAM,KAAG,MAErE,CAAC;AAEF,MAAM,MAAM,QAAQ,GAAG;IACrB,QAAQ,EAAE,IAAI,CAAC;IACf,CAAC,EAAE,MAAM,CAAC;IACV,CAAC,EAAE,MAAM,CAAC;IACV,QAAQ,EAAE;QAAE,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,KAAK,CAAA;KAAE,CAAC;IAC5D,KAAK,CAAC,EAAE,MAAM,CAAC;IACf,IAAI,EAAE,MAAM,CAAC;CACd,CAAC;AAYF,+FAA+F;AAC/F,eAAO,MAAM,WAAW,GACtB,QAAQ,MAAM,EACd,MAAM,SAAS,EACf,SAAS,aAAa,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACxC,QAAQ,MAAM,KACb;IAAE,IAAI,EAAE,QAAQ,CAAC;IAAC,MAAM,EAAE,OAAO,CAAA;CAoBnC,CAAC;AAEF;;;GAGG;AACH,eAAO,MAAM,eAAe,GAC1B,QAAQ,MAAM,EACd,MAAM,SAAS,EACf,SAAS,aAAa,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACxC,OAAO;IAAE,KAAK,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,MAAM,CAAA;CAAE,KACzC,QA+DF,CAAC"}

package/dist/scene/placement.js

@@ -0,0 +1,107 @@
+import { planeHit, pixelRay, projectPoints, rotationMatrix, } from './camera-model.js';
+import { multistart } from './least-squares.js';
+import { fixed, roundTo, SceneArgumentError } from './numbers.js';
+import { seededRandom } from './seeded-random.js';
+export const castRays = (camera, size, pixels, planeY) => pixels.map(([u, v]) => ({ u, v, ...planeHit(camera, size, u, v, planeY) }));
+export const formatRayLine = (line) => (`${fixed(line.u, 0)},${fixed(line.v, 0)} -> x=${fixed(line.point[0], 4)}`
+    + ` y=${fixed(line.point[1], 4)} z=${fixed(line.point[2], 4)}`);
+export const formatProjectLine = (point, u, v) => (`${fixed(point[0], 4)},${fixed(point[1], 4)},${fixed(point[2], 4)} -> px ${fixed(u, 1)},${fixed(v, 1)}`);
+const FLAT_NOTE = 'mesh face is local xy facing +z; rotation.x=-pi/2 lays it flat facing up,'
+    + ' local +y points away from the camera';
+const AMBIGUITY_NOTE = 'rectangles seen from one view are ambiguous in lean versus depth; if it lands inside'
+    + ' another object, constrain it with what you know (rests on a surface, leans on a wall)';
+const SOLVED_SIDE_NOTE = 'one dimension was solved, so lean, roll and that dimension trade against each other;'
+    + ' prefer --plane-y for anything lying flat. ';
+const subtract = (a, b) => [a[0] - b[0], a[1] - b[1], a[2] - b[2]];
+const norm = (a) => Math.hypot(...a);
+/** A rectangle lying flat on the plane y = planeY, any size. Corners TL;TR;BR;BL in pixels. */
+export const rectOnPlane = (camera, size, corners, planeY) => {
+    const hits = corners.map(([u, v]) => planeHit(camera, size, u, v, planeY));
+    const P = hits.map((hit) => hit.point);
+    const centre = [0, 0, 0];
+    for (const point of P)
+        for (let axis = 0; axis < 3; axis += 1)
+            centre[axis] += point[axis] / P.length;
+    const top = subtract(P[1], P[0]);
+    const bottom = subtract(P[2], P[3]);
+    const width = (norm(top) + norm(bottom)) / 2;
+    const height = (norm(subtract(P[3], P[0])) + norm(subtract(P[2], P[1]))) / 2;
+    const yaw = Math.atan2(-(top[2] + bottom[2]), top[0] + bottom[0]);
+    return {
+        pose: {
+            position: [roundTo(centre[0], 4), roundTo(centre[1], 4), roundTo(centre[2], 4)],
+            w: roundTo(width, 4),
+            h: roundTo(height, 4),
+            rotation: { x: roundTo(-Math.PI / 2, 6), y: roundTo(yaw, 6), z: 0, order: 'YXZ' },
+            note: FLAT_NOTE,
+        },
+        behind: hits.some((hit) => hit.behind),
+    };
+};
+/**
+ * Free pose (position, lean, roll) of a rectangle whose width, height or both
+ * are known. With one side known the other is solved, less well determined.
+ */
+export const rectOfKnownSize = (camera, size, corners, known) => {
+    const { width, height } = known;
+    if (width === undefined && height === undefined) {
+        throw new SceneArgumentError('rect needs --plane-y, or --width and/or --height');
+    }
+    const seed = Math.max(...[width, height].filter((value) => value !== undefined));
+    const solvedSide = width === undefined || height === undefined ? 1 : 0;
+    const centroid = corners.reduce((total, [u, v]) => [total[0] + u / corners.length, total[1] + v / corners.length], [0, 0]);
+    const ray = pixelRay(camera, size, centroid[0], centroid[1]);
+    const sides = (p) => [
+        width ?? p[6],
+        height ?? p[6],
+    ];
+    const residuals = (p) => {
+        const [w, h] = sides(p);
+        const local = [[-w / 2, h / 2, 0], [w / 2, h / 2, 0], [w / 2, -h / 2, 0], [-w / 2, -h / 2, 0]];
+        // Euler YXZ with rotation.set(x = p[4], y = p[3], z = p[5]); lean is rotation.x.
+        const r = rotationMatrix(p[3], p[4], p[5]);
+        const world = local.map((point) => [
+            r[0][0] * point[0] + r[0][1] * point[1] + r[0][2] * point[2] + p[0],
+            r[1][0] * point[0] + r[1][1] * point[1] + r[1][2] * point[2] + p[1],
+            r[2][0] * point[0] + r[2][1] * point[1] + r[2][2] * point[2] + p[2],
+        ]);
+        const projected = projectPoints(camera, world, size);
+        const out = new Float64Array(8);
+        projected.forEach((projection, index) => {
+            const penalty = Math.max(projection.depth + 0.01, 0) * 1e4;
+            out[2 * index] = projection.u - corners[index][0] + penalty;
+            out[2 * index + 1] = projection.v - corners[index][1] + penalty;
+        });
+        return out;
+    };
+    const random = seededRandom(2);
+    const starts = [];
+    for (let index = 0; index < 300; index += 1) {
+        const reach = random.uniform(0.2, 3) * seed * 3;
+        const start = [
+            camera.pos[0] + ray[0] * reach,
+            camera.pos[1] + ray[1] * reach,
+            camera.pos[2] + ray[2] * reach,
+            random.uniform(-1.2, 1.2),
+            random.uniform(-1.4, 1.4),
+            random.uniform(-0.6, 0.6),
+        ];
+        if (solvedSide)
+            start.push(seed * random.uniform(0.5, 2));
+        starts.push(Float64Array.from(start));
+    }
+    const low = [-1e4, -1e4, -1e4, -4, -1.6, -3.2, ...(solvedSide ? [seed * 0.05] : [])];
+    const high = [1e4, 1e4, 1e4, 4, 1.6, 3.2, ...(solvedSide ? [seed * 20] : [])];
+    const best = multistart(residuals, starts, low, high);
+    const p = best.parameters;
+    const [w, h] = sides(p);
+    return {
+        position: [roundTo(p[0], 4), roundTo(p[1], 4), roundTo(p[2], 4)],
+        w: roundTo(w, 4),
+        h: roundTo(h, 4),
+        rotation: { x: roundTo(p[4], 6), y: roundTo(p[3], 6), z: roundTo(p[5], 6), order: 'YXZ' },
+        rmsPx: roundTo(Math.sqrt(best.cost / 8), 2),
+        note: (solvedSide ? SOLVED_SIDE_NOTE : '') + AMBIGUITY_NOTE,
+    };
+};
+//# sourceMappingURL=placement.js.map

package/dist/scene/placement.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"placement.js","sourceRoot":"","sources":["../../src/scene/placement.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,QAAQ,EACR,QAAQ,EACR,aAAa,EACb,cAAc,GAIf,MAAM,mBAAmB,CAAC;AAC3B,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClE,OAAO,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AAWlD,MAAM,CAAC,MAAM,QAAQ,GAAG,CACtB,MAAc,EACd,IAAe,EACf,MAAuC,EACvC,MAAc,EACH,EAAE,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;AAE5F,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,IAAa,EAAU,EAAE,CAAC,CACtD,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,SAAS,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE;MACvE,MAAM,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAC/D,CAAC;AAEF,MAAM,CAAC,MAAM,iBAAiB,GAAG,CAAC,KAAW,EAAE,CAAS,EAAE,CAAS,EAAU,EAAE,CAAC,CAC9E,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CACxG,CAAC;AAWF,MAAM,SAAS,GAAG,2EAA2E;MACzF,uCAAuC,CAAC;AAC5C,MAAM,cAAc,GAAG,sFAAsF;MACzG,wFAAwF,CAAC;AAC7F,MAAM,gBAAgB,GAAG,sFAAsF;MAC3G,6CAA6C,CAAC;AAElD,MAAM,QAAQ,GAAG,CAAC,CAAO,EAAE,CAAO,EAAQ,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACrF,MAAM,IAAI,GAAG,CAAC,CAAO,EAAU,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;AAEnD,+FAA+F;AAC/F,MAAM,CAAC,MAAM,WAAW,GAAG,CACzB,MAAc,EACd,IAAe,EACf,OAAwC,EACxC,MAAc,EACuB,EAAE;IACvC,MAAM,IAAI,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC;IAC3E,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IACvC,MAAM,MAAM,GAAS,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/B,KAAK,MAAM,KAAK,IAAI,CAAC;QAAE,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,IAAI,CAAC;YAAE,MAAM,CAAC,IAAI,CAAE,IAAI,KAAK,CAAC,IAAI,CAAE,GAAG,CAAC,CAAC,MAAM,CAAC;IACxG,MAAM,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;IACtC,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;IAC7C,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACjF,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,OAAO;QACL,IAAI,EAAE;YACJ,QAAQ,EAAE,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC/E,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;YACpB,CAAC,EAAE,OAAO,CAAC,MAAM,EAAE,CAAC,CAAC;YACrB,QAAQ,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE;YACjF,IAAI,EAAE,SAAS;SAChB;QACD,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,CAAC;KACvC,CAAC;AACJ,CAAC,CAAC;AAEF;;;GAGG;AACH,MAAM,CAAC,MAAM,eAAe,GAAG,CAC7B,MAAc,EACd,IAAe,EACf,OAAwC,EACxC,KAA0C,EAChC,EAAE;IACZ,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,GAAG,KAAK,CAAC;IAChC,IAAI,KAAK,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;QAChD,MAAM,IAAI,kBAAkB,CAAC,kDAAkD,CAAC,CAAC;IACnF,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC,MAAM,CAAC,CAAC,KAAK,EAAmB,EAAE,CAAC,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC;IAClG,MAAM,UAAU,GAAG,KAAK,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvE,MAAM,QAAQ,GAAG,OAAO,CAAC,MAAM,CAC7B,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC,EACjF,CAAC,CAAC,EAAE,CAAC,CAAC,CACP,CAAC;IACF,MAAM,GAAG,GAAG,QAAQ,CAAC,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,MAAM,KAAK,GAAG,CAAC,CAAoB,EAAoB,EAAE,CAAC;QACxD,KAAK,IAAI,CAAC,CAAC,CAAC,CAAE;QACd,MAAM,IAAI,CAAC,CAAC,CAAC,CAAE;KAChB,CAAC;IACF,MAAM,SAAS,GAAG,CAAC,CAAe,EAAgB,EAAE;QAClD,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,KAAK,GAAW,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACvG,iFAAiF;QACjF,MAAM,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;QAC9C,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,KAAK,EAAQ,EAAE,CAAC;YACvC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE;YACpE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE;YACpE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE;SACrE,CAAC,CAAC;QACH,MAAM,SAAS,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;QACrD,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAChC,SAAS,CAAC,OAAO,CAAC,CAAC,UAAU,EAAE,KAAK,EAAE,EAAE;YACtC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,KAAK,GAAG,IAAI,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;YAC3D,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,UAAU,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAE,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;YAC7D,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAE,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QACnE,CAAC,CAAC,CAAC;QACH,OAAO,GAAG,CAAC;IACb,CAAC,CAAC;IACF,MAAM,MAAM,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,MAAM,GAAmB,EAAE,CAAC;IAClC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,GAAG,EAAE,KAAK,IAAI,CAAC,EAAE,CAAC;QAC5C,MAAM,KAAK,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;QAChD,MAAM,KAAK,GAAG;YACZ,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK;YAC9B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK;YAC9B,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK;YAC9B,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC;YACzB,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC;YACzB,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC;SAC1B,CAAC;QACF,IAAI,UAAU;YAAE,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1D,MAAM,CAAC,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;IACxC,CAAC;IACD,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACrF,MAAM,IAAI,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC9E,MAAM,IAAI,GAAG,UAAU,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;IACtD,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;IAC1B,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACxB,OAAO;QACL,QAAQ,EAAE,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC;QACnE,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,CAAC,EAAE,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,QAAQ,EAAE,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE;QAC5F,KAAK,EAAE,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAC3C,IAAI,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,cAAc;KAC5D,CAAC;AACJ,CAAC,CAAC"}

package/dist/scene/raster.d.ts

@@ -0,0 +1,18 @@
+/**
+ * 8-bit RGB pixels in memory: the one image shape every scene helper works on.
+ * Alpha is dropped without compositing, as the Python helpers' convert("RGB")
+ * did, and before any resize so a transparent edge never bleeds into colour.
+ */
+export type RgbImage = {
+    data: Buffer;
+    width: number;
+    height: number;
+};
+export declare const readRgb: (file: string) => Promise<RgbImage>;
+/** Lanczos resize to an exact size; the same size is returned untouched. */
+export declare const resizeRgb: (image: RgbImage, width: number, height: number) => Promise<RgbImage>;
+export declare const readRgbAt: (file: string, width: number, height: number) => Promise<RgbImage>;
+/** PNG bytes of an in-memory image. */
+export declare const encodePng: (image: RgbImage) => Promise<Buffer>;
+export declare const writePng: (image: RgbImage, file: string) => Promise<void>;
+//# sourceMappingURL=raster.d.ts.map

package/dist/scene/raster.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"raster.d.ts","sourceRoot":"","sources":["../../src/scene/raster.ts"],"names":[],"mappings":"AAEA;;;;GAIG;AACH,MAAM,MAAM,QAAQ,GAAG;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAAC;AAEvE,eAAO,MAAM,OAAO,GAAU,MAAM,MAAM,KAAG,OAAO,CAAC,QAAQ,CAQ5D,CAAC;AAMF,4EAA4E;AAC5E,eAAO,MAAM,SAAS,GAAU,OAAO,QAAQ,EAAE,OAAO,MAAM,EAAE,QAAQ,MAAM,KAAG,OAAO,CAAC,QAAQ,CAOhG,CAAC;AAEF,eAAO,MAAM,SAAS,GAAU,MAAM,MAAM,EAAE,OAAO,MAAM,EAAE,QAAQ,MAAM,KAAG,OAAO,CAAC,QAAQ,CAE7F,CAAC;AAEF,uCAAuC;AACvC,eAAO,MAAM,SAAS,GAAU,OAAO,QAAQ,KAAG,OAAO,CAAC,MAAM,CAE/D,CAAC;AAEF,eAAO,MAAM,QAAQ,GAAU,OAAO,QAAQ,EAAE,MAAM,MAAM,KAAG,OAAO,CAAC,IAAI,CAE1E,CAAC"}